title: MOND Prediction for the Velocity Dispersion of the `Feeble Giant' Crater II

abstract: Crater II is an unusual object among the dwarf satellite galaxies of the Local Group in that it has a very large size for its small luminosity. This provides a strong test of MOND, as Crater II should be in the deep MOND regime ($g_{in} \approx 34\;\mathrm{km}^2\,\mathrm{s}^{-2}\,\mathrm{kpc}^{-1} \ll a_0 = 3700\;\mathrm{km}^2\,\mathrm{s}^{-2}\,\mathrm{kpc}^{-1}$). Despite its great distance ($\approx 120$ kpc) from the Milky Way, the external field of the host ($g_{ex} \approx 282\; \mathrm{km}^2\,\mathrm{s}^{-2}\,\mathrm{kpc}^{-1}$) comfortably exceeds the internal field. Consequently, Crater II should be subject to the external field effect, a feature unique to MOND. This leads to the prediction of a very low velocity dispersion: $\sigma_{efe} = 2.1^{+0.9}_{-0.6}\;\mathrm{km}\,\mathrm{s}^{-1}$.


title: Simulations of the Fe K-alpha Energy Spectra from Gravitationally Microlensed Quasars

abstract: The analysis of the Chandra X-ray observations of the gravitationally lensed quasar RX J1131-1231 revealed the detection of multiple and energy-variable spectral peaks. The spectral variability is thought to result from the microlensing of the Fe K-alpha emission, selectively amplifying the emission from certain regions of the accretion disk with certain effective frequency shifts of the Fe K-alpha line emission. In this paper, we combine detailed simulations of the emission of Fe K-alpha photons from the accretion disk of a Kerr black hole with calculations of the effect of gravitational microlensing on the observed energy spectra. The simulations show that microlensing can indeed produce multiply peaked energy spectra. We explore the dependence of the spectral characteristics on black hole spin, accretion disk inclination, corona height, and microlensing amplification factor, and show that the measurements can be used to constrain these parameters. We find that the range of observed spectral peak energies of QSO RX J1131-1231 can only be reproduced for black hole inclinations exceeding 70 degree and for lamppost corona heights of less than 30 gravitational radii above the black hole. We conclude by emphasizing the scientific potential of studies of the microlensed Fe K$\alpha$ quasar emission and the need for more detailed modeling that explores how the results change for more realistic accretion disk and corona geometries and microlensing magnification patterns. A full analysis should furthermore model the signal-to-noise ratio of the observations and the resulting detection biases.


title: Interplay of non-symmorphic symmetry and spin-orbit coupling in hyperkagome spin liquids: Applications to Na$_4$Ir$_3$O$_8$

abstract: Na$_4$Ir$_3$O$_8$ provides a material platform to study three-dimensional quantum spin liquids in the geometrically frustrated hyperkagome lattice of Ir$^{4+}$ ions. In this work, we consider quantum spin liquids on hyperkagome lattice for generic spin models, focusing on the effects of anisotropic spin interactions. In particular, we classify possible $\mathbb{Z}_2$ and $U(1)$ spin liquid states, following the projective symmetry group analysis in the slave-fermion representation. There are only three distinct $\mathbb{Z}_2$ spin liquids, together with 2 different $U(1)$ spin liquids. The non-symmorphic space group symmetry of hyperkagome lattice plays a vital role in simplifying the classification, forbidding "$\pi$-flux" or "staggered-flux" phases in contrast to symmorphic space groups. We further prove that both $U(1)$ states and one $Z_2$ state among all 3 are symmetry-protected gapless spin liquids, robust against any symmetry-preserving perturbations. Motivated by the "spin-freezing" behavior recently observed in Na$_4$Ir$_3$O$_8$ at low temperatures, we further investigate the nearest-neighbor spin model with dominant Heisenberg interaction subject to all possible anisotropic perturbations from spin-orbit couplings. We found a $U(1)$ spin liquid ground state with spinon fermi surfaces is energetically favored over $Z_2$ states. Among all spin-orbit coupling terms, we show that only Dzyaloshinskii-Moriya (DM) interaction can induce spin anisotropy in the ground state when perturbing from the isotropic Heisenberg limit. Our work paves the way for a systematic study of quantum spin liquids in various materials with a hyperkagome crystal structure.


title: Initial conditions for inflation

abstract: Within the $\alpha$-attractors framework we investigate scalar potentials with the same pole as the one featured in the kinetic term. We show that, in field space, this leads to directions without a plateau. Using this, we present a proposal, which manages to overcome the initial conditions problem of inflation with a plateau. An earlier period of proto-inflation, beginning at Planck scale, accounts for the Universe expansion and arranges the required initial conditions for inflation on the plateau to commence. We show that, if proto-inflation is power-law, it does not suffer from a sub-Planckian eternal inflationary stage, which would otherwise be a problem. A simple model realisation is constructed in the context of $\alpha$-attractors, which can both generate the inflationary plateau and the exponential slopes around it, necessary for the two inflation stages. Our mechanism allows to assume chaotic initial conditions at the Planck scale for proto-inflation, it is generic and it is shown to work without fine-tuning.


title: Breit-Wigner approximation for propagators of mixed unstable states

abstract: For systems of unstable particles that mix with each other, an approximation of the fully momentum-dependent propagator matrix is presented in terms of a sum of simple Breit-Wigner propagators that are multiplied with finite on-shell wave function normalisation factors. The latter are evaluated at the complex poles of the propagators. The pole structure of general propagator matrices is carefully analysed, and it is demonstrated that in the proposed approximation imaginary parts arising from absorptive parts of loop integrals are properly taken into account. Applying the formalism to the neutral MSSM Higgs sector with complex parameters, very good numerical agreement is found between cross sections based on the full propagators and the corresponding cross sections based on the described approximation. The proposed approach does not only technically simplify the treatment of propagators with non-vanishing off-diagonal contributions, it is shown that it can also facilitate an improved theoretical prediction of the considered observables via a more precise implementation of the total widths of the involved particles. It is also well-suited for the incorporation of interference effects arising from overlapping resonances.


title: Robust and Parallel Bayesian Model Selection

abstract: Effective and accurate model selection is an important problem in modern data analysis. One of the major challenges is the computational burden required to handle large data sets that cannot be stored or processed on one machine. Another challenge one may encounter is the presence of outliers and contaminations that damage the inference quality. The parallel "divide and conquer" model selection strategy divides the observations of the full data set into roughly equal subsets and perform inference and model selection independently on each subset. After local subset inference, this method aggregates the posterior model probabilities or other model/variable selection criteria to obtain a final model by using the notion of geometric median. This approach leads to improved concentration in finding the "correct" model and model parameters and also is provably robust to outliers and data contamination.


title: Nonstationary POT modelling of air pollution concentrations: Statistical analysis of the traffic and meteorological impact

abstract: Predicting the occurrence, level and duration of high air pollution concentrations exceeding a given critical level enables researchers to study the health impact of road traffic on local air quality and to inform public policy action. Precise estimates of the probabilities of occurrence and level of extreme concentrations are formidable due to the combination of complex physical and chemical processes involved. This underpins the need for developing sophisticated extreme value models, in particular allowing for non-stationarity of environmental time series. In this paper, extremes of nitrogen oxide (NO), nitrogen dioxide (NO$_2$) and ozone (O$_3$) concentrations are investigated using two models. Model I is based on an extended peaks-over-threshold (POT) approach developed by A. C. Davison and R. L. Smith, whereby the parameters of the underlying generalized Pareto distribution (GPD) are treated as functions of covariates (i.e., traffic and meteorological factors). The new Model II resolves the lack of threshold stability in the Davison--Smith model by constructing a special functional form for the GPD parameters. For each of the models, the effects of trafic and meteorological factors on the frequency and size of extreme values are estimated using Markov chain Monte Carlo methods. Finally, appropriate goodness-of-fit tests and model selection criteria confirm that Model II significantly outperforms Model I in estimation and forecasting of extremes.


title: Structure of radiation dominated gravitoturbulent quasar discs

abstract: Self-gravitating accretion discs in a gravitoturbulent state, including radiation and gas pressures, are studied using a set of new analytical solutions. While the Toomre parameter of the disc remains close to its critical value for the onset of gravitational instability, the dimensionless stress parameter is uniquely determined from the thermal energy reservoir of the disc and its cooling rate. Our solutions are applicable to the accretion discs with dynamically important radiation pressure like in the quasars discs. We show that physical quantities of a gravitoturbulent disc in the presence of radiation are significantly modified compared to solutions with only gas pressure. We show that the dimensionless stress parameter is an increasing function of the radial distance so that its steepness strongly depends on the accretion rate. In a disc without radiation its slope is 4.5, however, we show that in the presence of radiation, it varies between 2 and 4.5 depending on the accretion rate and the central mass. As for the surface density, we find a shallower profile with an exponent -2 in a disc with sub-Eddington accretion rate compared to a similar disc, but without radiation, where its surface density slope is -3 independent of the accretion rate. We then investigate gravitational stability of the disc when the stress parameter reaches to its critical value. In order to self-consistently determine the fragmentation boundary, however, it is shown that the critical value of the stress parameter is a power-law function of the ratio of gas pressure and the total pressure and its exponent is around 1.7. We also estimate the maximum mass of the central black hole using our analytical solutions.


title: Helicity Evolution at Small $x$: Flavor Singlet and Non-Singlet Observables

abstract: We extend our earlier results for the quark helicity evolution at small $x$ to derive the small-$x$ asymptotics of the flavor singlet and flavor non-singlet quark helicity TMDs and PDFs and of the $g_1$ structure function. In the flavor singlet case we re-derive the evolution equations obtained in our previous paper on the subject, performing additional cross-checks of our results. In the flavor non-singlet case we construct new small-$x$ evolution equations by employing the large-$N_c$ limit. All evolution equations resum double-logarithmic powers of $\alpha_s \, \ln^2 (1/x)$ in the polarization-dependent evolution along with the single-logarithmic powers of $\alpha_s \, \ln (1/x)$ in the unpolarized evolution which includes saturation effects. We solve the linearized flavor non-singlet equation analytically, obtaining an intercept which agrees with the one calculated earlier by Bartels, Ermolaev and Ryskin using the infra-red evolution equations. Our numerical solution of the linearized large-$N_c$ evolution equations for the flavor singlet case is presented in the accompanying Letter and is further discussed here.


title: Localization of interacting Fermi gases in quasiperiodic potentials

abstract: We investigate the zero-temperature metal-insulator transition in a one-dimensional two-component Fermi gas in the presence of a quasi-periodic potential resulting from the superposition of two optical lattices of equal intensity but incommensurate periods. A mobility edge separating (low energy) Anderson localized and (high energy) extended single-particle states appears in this continuous-space model beyond a critical intensity of the quasi-periodic potential. In order to discern the metallic phase from the insulating phase in the interacting many-fermion system, we employ unbiased quantum Monte Carlo (QMC) simulations combined with the many-particle localization length familiar from the modern theory of the insulating state. In the noninteracting limit, the critical optical-lattice intensity for the metal-insulator transition predicted by the QMC simulations coincides with the Anderson localization transition of the single-particle eigenstates. We show that weak repulsive interactions induce a shift of this critical point towards larger intensities, meaning that repulsion favors metallic behavior. This shift appears to be linear in the interaction parameter, suggesting that even infinitesimal interactions can affect the position of the critical point.


title: Better Streaming Algorithms for the Maximum Coverage Problem

abstract: We study the classic NP-Hard problem of finding the maximum $k$-set coverage in the data stream model: given a set system of $m$ sets that are subsets of a universe $\{1,\ldots,n \}$, find the $k$ sets that cover the most number of distinct elements. The problem can be approximated up to a factor $1-1/e$ in polynomial time. In the streaming-set model, the sets and their elements are revealed online. The main goal of our work is to design algorithms, with approximation guarantees as close as possible to $1-1/e$, that use sublinear space $o(mn)$. Our main results are: Two $(1-1/e-\epsilon)$ approximation algorithms: One uses $O(\epsilon^{-1})$ passes and $\tilde{O}(\epsilon^{-2} k)$ space whereas the other uses only a single pass but $\tilde{O}(\epsilon^{-2} m)$ space. We show that any approximation factor better than $(1-(1-1/k)^k)$ in constant passes requires $\Omega(m)$ space for constant $k$ even if the algorithm is allowed unbounded processing time. We also demonstrate a single-pass, $(1-\epsilon)$ approximation algorithm using $\tilde{O}(\epsilon^{-2} m \cdot \min(k,\epsilon^{-1}))$ space. We also study the maximum $k$-vertex coverage problem in the dynamic graph stream model. In this model, the stream consists of edge insertions and deletions of a graph on $N$ vertices. The goal is to find $k$ vertices that cover the most number of distinct edges. We show that any constant approximation in constant passes requires $\Omega(N)$ space for constant $k$ whereas $\tilde{O}(\epsilon^{-2}N)$ space is sufficient for a $(1-\epsilon)$ approximation and arbitrary $k$ in a single pass. For regular graphs, we show that $\tilde{O}(\epsilon^{-3}k)$ space is sufficient for a $(1-\epsilon)$ approximation in a single pass. We generalize this to a $(\kappa-\epsilon)$ approximation when the ratio between the minimum and maximum degree is bounded below by $\kappa$.


title: The Tungsten Project: Dielectronic Recombination Data For Collisional-Radiative Modelling In ITER

abstract: Tungsten is an important metal in nuclear fusion reactors. It will be used in the divertor component of ITER (Latin for 'the way'). The Tungsten Project aims to calculate partial and total DR rate coefficients for the isonuclear sequence of Tungsten. The calculated data will be made available as and when they are produced via the open access database OPEN- ADAS in the standard adf09 and adf48 file formats. We present our progress thus far, detailing calculational methods, and showing comparisons with other available data. We conclude with plans for the future.


title: Two constructions of virtually contact structures

abstract: Motivated by recent developments in proving the Weinstein conjecture we introduce the notion of covering contact connected sum for virtually contact manifolds and construct virtually contact structures on boundaries of subcritical handle bodies.


title: Plancherel decomposition of Howe duality and Euler factorization of automorphic functionals

abstract: There are several global functionals on irreducible automorphic representations which are Eulerian, that is: pure tensors of local functionals, when the representation is written as an Euler product $\pi = \otimes'_v \pi_v$ of local representations. The precise factorization of such functionals is of interest to number theorists and is -- naturally -- very often related to special values of $L$-functions. The purpose of this paper is to develop in full generality the Plancherel formula for the Weil or oscillator representation, considered as a unitary representation of a reductive dual pair, and to use it in order to demonstrate a very general principle of Euler factorization: local factors are determined via the Langlands correspondence by a local Plancherel formula. This pattern has already been observed and conjectured in the author's prior work with Venkatesh in the case of period integrals. Here, it is shown that the Rallis inner product formula amounts to the same principle in the setting of global Howe duality.


title: A high-order staggered meshless method for elliptic problems

abstract: We present a new meshless method for scalar diffusion equations which is motivated by their compatible discretizations on primal-dual grids. Unlike the latter though, our approach is truly meshless because it only requires the graph of nearby neighbor connectivity of the discretization points $\bm{x}_i$. This graph defines a local primal-dual grid complex with a \emph{virtual} dual grid, in the sense that specification of the dual metric attributes is implicit in the method's construction. Our method combines a topological gradient operator on the local primal grid with a Generalized Moving Least Squares approximation of the divergence on the local dual grid. We show that the resulting approximation of the div-grad operator maintains polynomial reproduction to arbitrary orders and yields a meshless method, which attains $O(h^{m})$ convergence in both $L^2$ and $H^1$ norms, similar to mixed finite element methods. We demonstrate this convergence on curvilinear domains using manufactured solutions. Application of the new method to problems with discontinuous coefficients reveals solutions that are qualitatively similar to those of compatible mesh-based discretizations.


title: A Reinforcement Learning Approach to the View Planning Problem

abstract: We present a Reinforcement Learning (RL) solution to the view planning problem (VPP), which generates a sequence of view points that are capable of sensing all accessible area of a given object represented as a 3D model. In doing so, the goal is to minimize the number of view points, making the VPP a class of set covering optimization problem (SCOP). The SCOP is NP-hard, and the inapproximability results tell us that the greedy algorithm provides the best approximation that runs in polynomial time. In order to find a solution that is better than the greedy algorithm, (i) we introduce a novel score function by exploiting the geometry of the 3D model, (ii) we model an intuitive human approach to VPP using this score function, and (iii) we cast VPP as a Markovian Decision Process (MDP), and solve the MDP in RL framework using well-known RL algorithms. In particular, we use SARSA, Watkins-Q and TD with function approximation to solve the MDP. We compare the results of our method with the baseline greedy algorithm in an extensive set of test objects, and show that we can out-perform the baseline in almost all cases.


title: Natural generalized mirage mediation

abstract: In the supersymmetric scenario known as mirage mediation (MM), the soft SUSY breaking terms receive comparable anomaly-mediation and moduli-mediation contributions leading to the phenomenon of mirage unification. The simplest MM SUSY breaking models which are consistent with the measured Higgs mass and sparticle mass constraints are strongly disfavoured by fine-tuning considerations. However, while MM makes robust predictions for gaugino masses, the scalar sector is quite sensitive to specific mechanisms for moduli stabilization and potential uplifting. We suggest here a broader setup of generalized mirage mediation (GMM), where heretofore discrete parameters are allowed as continuous to better parametrize these other schemes. We find that natural SUSY spectra consistent with both the measured value of m(h). as well as LHC lower bounds on superpartner masses are then possible. We explicitly show that models generated from natural GMM may be beyond the reach of even high-luminosity LHC searches. In such a case, the proposed International Linear e^+e^- Collider (ILC) will be required for natural SUSY discovery via higgsino pair production reactions. We also outline prospects for detection of higgsino-like WIMPs from natural GMM.


title: HERUS: A CO Atlas from SPIRE Spectroscopy of local ULIRGs

abstract: We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671\mu$m with photometry in the SPIRE bands at 250$\mu$m, 350$\mu$m and 500$\mu$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $\mu$m, [CI] 609 $\mu$m, and [NII] 205 $\mu$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}


title: Elliptic Polylogarithms and Basic Hypergeometric Functions

abstract: Multiple elliptic polylogarithms can be written as (multiple) integrals of products of basic hypergeometric functions. The latter are computable, to arbitrary precision, using a q-difference equation and q-contiguous relations.


title: A Loomis-Sikorski theorem and functional calculus for a generalized Hermitian algebra

abstract: A generalized Hermitian (GH-) algebra is a generalization of the partially ordered Jordan algebra of all Hermitian operators on a Hilbert space. We introduce the notion of a gh-tribe, which is a commutative GH-algebra of functions on a nonempty set $X$ with pointwise partial order and operations, and we prove that every commutative GH-algebra is the image of a gh-tribe under a surjective GH-morphism. Using this result, we prove each element $a$ of a GH-algebra $A$ corresponds to a real observable $\xi_a$ on the $\sigma$-orthomodular lattice of projections in $A$ and that $\xi_a$ determines the spectral resolution of $a$. Also, if $f$ is a continuous function defined on the spectrum of $a$, we formulate a definition of $f(a)$, thus obtaining a continuous functional calculus for $A$.


title: Structured adaptive and random spinners for fast machine learning computations

abstract: We consider an efficient computational framework for speeding up several machine learning algorithms with almost no loss of accuracy. The proposed framework relies on projections via structured matrices that we call Structured Spinners, which are formed as products of three structured matrix-blocks that incorporate rotations. The approach is highly generic, i.e. i) structured matrices under consideration can either be fully-randomized or learned, ii) our structured family contains as special cases all previously considered structured schemes, iii) the setting extends to the non-linear case where the projections are followed by non-linear functions, and iv) the method finds numerous applications including kernel approximations via random feature maps, dimensionality reduction algorithms, new fast cross-polytope LSH techniques, deep learning, convex optimization algorithms via Newton sketches, quantization with random projection trees, and more. The proposed framework comes with theoretical guarantees characterizing the capacity of the structured model in reference to its unstructured counterpart and is based on a general theoretical principle that we describe in the paper. As a consequence of our theoretical analysis, we provide the first theoretical guarantees for one of the most efficient existing LSH algorithms based on the HD3HD2HD1 structured matrix [Andoni et al., 2015]. The exhaustive experimental evaluation confirms the accuracy and efficiency of structured spinners for a variety of different applications.


title: A Theme-Rewriting Approach for Generating Algebra Word Problems

abstract: Texts present coherent stories that have a particular theme or overall setting, for example science fiction or western. In this paper, we present a text generation method called {\it rewriting} that edits existing human-authored narratives to change their theme without changing the underlying story. We apply the approach to math word problems, where it might help students stay more engaged by quickly transforming all of their homework assignments to the theme of their favorite movie without changing the math concepts that are being taught. Our rewriting method uses a two-stage decoding process, which proposes new words from the target theme and scores the resulting stories according to a number of factors defining aspects of syntactic, semantic, and thematic coherence. Experiments demonstrate that the final stories typically represent the new theme well while still testing the original math concepts, outperforming a number of baselines. We also release a new dataset of human-authored rewrites of math word problems in several themes.


title: Diffractive processes at the LHC within $k_t$ -factorization approach

abstract: We discuss the single diffractive production of $c \bar c$ pairs and charmed mesons at the LHC. In addition to standard collinear approach, for a first time we propose a $k_t$-factorization approach to the diffractive processes. The transverse momentum dependent (the unintegrated diffractive parton distributions) in proton are obtained with the help of the Kimber-Martin-Ryskin prescription where collinear diffractive PDFs are used as input. In this calculation the transverse momentum of the pomeron is neglected with respect to transverse momentum of partons entering the hard process. The results of the new approach are compared with those of the standard collinear one. Significantly larger cross sections are obtained in the $k_t$-factorization approach where some part of higher-order effects is effectively included. Some correlation observables, like azimuthal angle correlation between $c$ and $\bar c$, and $c \bar c$ pair transverse momentum distribution were obtained for the first time.


title: RadioHound: A Pervasive Sensing Network for Sub-6 GHz Dynamic Spectrum Monitoring

abstract: We design a custom spectrum sensing network, called RadioHound, capable of tuning from 25 MHz to 6 GHz, which covers nearly all widely-deployed wireless activity. We describe the system hardware and network infrastructure in detail with a view towards driving the cost, size, and power usage of the sensors as low as possible. The system estimates the spatial variation of radio-frequency power from an unknown random number of sources. System performance is measured by computing the mean square error against a simulated radio-frequency environment. We find that the system performance depends heavily on the deployment density of the sensors. Consequently, we derive an expression for the sensor density as a function of environmental characteristics and confidence in measurement quality.


title: Direct Photon Elliptic Flow at BNL Relativistic Heavy Ion Collider (RHIC) and CERN Large Hadron Collider (LHC)

abstract: We use an event-by-event hydrodynamical description of the heavy-ion collision process with Glauber initial conditions to calculate the thermal emission of photons. The photon rates in the hadronic phase follow from a spectral function approach and a density expansion, while in the partonic phase they follow from the AMY perturbative rates. The calculated photon elliptic flows are lower than those reported recently by both the ALICE and PHENIX collaborations.


title: A model of infant speech perception and learning

abstract: Infant speech perception and learning is modeled using Echo State Network classification and Reinforcement Learning. Ambient speech for the modeled infant learner is created using the speech synthesizer Vocaltractlab. An auditory system is trained to recognize vowel sounds from a series of speakers of different anatomies in Vocaltractlab. Having formed perceptual targets, the infant uses Reinforcement Learning to imitate his ambient speech. A possible way of bridging the problem of speaker normalisation is proposed, using direct imitation but also including a caregiver who listens to the infants sounds and imitates those that sound vowel-like.


title: Density-dependent clustering: I. Pulling back the curtains on motions of the BAO peak

abstract: The most common statistic used to analyze large-scale structure surveys is the correlation function, or power spectrum. Here, we show how `slicing' the correlation function on local density brings sensitivity to interesting non-Gaussian features in the large-scale structure, such as the expansion or contraction of baryon acoustic oscillations (BAO) according to the local density. The sliced correlation function measures the large-scale flows that smear out the BAO, instead of just correcting them as reconstruction algorithms do. Thus, we expect the sliced correlation function to be useful in constraining the growth factor, and modified gravity theories that involve the local density. Out of the studied cases, we find that the run of the BAO peak location with density is best revealed when slicing on a $\sim 40$ Mpc/$h$ filtered density. But slicing on a $\sim100$ Mpc/$h$ filtered density may be most useful in distinguishing between underdense and overdense regions, whose BAO peaks are separated by a substantial $\sim 5$ Mpc/$h$ at $z=0$. We also introduce `curtain plots' showing how local densities drive particle motions toward or away from each other over the course of an $N$-body simulation.


title: A Versatile Technique to Enable sub-milli-Kelvin Instrument Stability for Precise Radial Velocity Measurements: Tests with the Habitable-zone Planet Finder

abstract: Insufficient instrument thermo-mechanical stability is one of the many roadblocks for achieving 10cm/s Doppler radial velocity (RV) precision, the precision needed to detect Earth-twins orbiting Solar-type stars. Highly temperature and pressure stabilized spectrographs allow us to better calibrate out instrumental drifts, thereby helping in distinguishing instrumental noise from astrophysical stellar signals. We present the design and performance of the Environmental Control System (ECS) for the Habitable-zone Planet Finder (HPF), a high-resolution (R=50,000) fiber-fed near infrared (NIR) spectrograph for the 10m Hobby Eberly Telescope at McDonald Observatory. HPF will operate at 180K, driven by the choice of an H2RG NIR detector array with a 1.7micron cutoff. This ECS has demonstrated 0.6mK RMS stability over 15 days at both 180K and 300K, and maintained high quality vacuum (<$10^{-7}$Torr) over months, during long-term stability tests conducted without a planned passive thermal enclosure surrounding the vacuum chamber. This control scheme is versatile and can be applied as a blueprint to stabilize future NIR and optical high precision Doppler instruments over a wide temperature range from ~77K to elevated room temperatures. A similar ECS is being implemented to stabilize NEID, the NASA/NSF NN-EXPLORE spectrograph for the 3.5m WIYN telescope at Kitt Peak, operating at 300K. A full SolidWorks 3D-CAD model and a comprehensive parts list of the HPF ECS are included with this manuscript to facilitate the adaptation of this versatile environmental control scheme in the broader astronomical community.


title: Generalized K-Shift Forbidden Substrings in Permutations

abstract: In this note we continue the analysis of permutations that avoid substrings j(j+k), 1 <= j <= n-k, k < n, as well as substrings j(j+k) (mod n), 1 <= j <= n. In the first case the number of such permutations can be obtained from recursions starting from derangement numbers, while in the (mod n) case the number of permutations depends on whether n is prime, and more generally, on whether n and k are relatively prime.


title: Dynamical Behavior Investigation and Analysis of Novel Mechanism for Simulated Spherical Robot named "RollRoller"

abstract: This paper introduces a simulation study of fluid actuated multi-driven closed system as spherical mobile robot called "RollRoller". Robot's mechanism design consists of two essential parts: tubes to lead a core and mechanical controlling parts to correspond movements. Our robot gets its motivation force by displacing the spherical movable mass known as core in curvy manners inside certain pipes. This simulation investigates by explaining the mechanical and structural features of the robot for creating hydraulic-base actuation via force and momentum analysis. Next, we categorize difficult and integrated 2D motions to omit unstable equilibrium points through derived nonlinear dynamics. We propose an algorithmic position control in forward direction that creates hybrid model as solution for motion planning problem in spherical robot. By deriving nonlinear dynamics of the spherical robot and implementing designed motion planning, we show how RollRoller can be efficient in high speed movements in comparison to the other pendulum-driven models. Then, we validate the results of this position control obtained by nonlinear dynamics via Adams/view simulation which uses the imported solid model of RollRoller. Lastly, We have a look to the circular maneuver of this robot by the same simulator.


title: Wave-particle Interactions in a Resonant System of Photons and Ion-solvated Water

abstract: We investigate a laser model for a resonant system of photons and ion cluster-solvated rotating water molecules in which ions in the cluster are identical and have very low, non-relativistic velocities and direction of motion parallel to a static electric field induced in a single direction. This model combines Dicke superradiation with wave-particle interaction. As the result, we find that the equations of motion of the system are expressed in terms of a conventional free electron laser system. This result leads to a mechanism for dynamical coherence, induced by collective instability in the wave-particle interaction.


title: On the recent seismic activity at Ioannina (Greece): Pre-earthquake electromagnetic emissions with critical and tricritical behavior

abstract: In this short paper we report that the dynamics of the earthquake (EQ) preparation processes, as embedded in the observed fracture-induced electromagnetic emissions (EME), were characterized by critical and then tricritical behavior a few days before the occurrence of the recent Ioannina (Greece) earthquake [($39.79^o N, 20.69^o E$), 15 October 2016, $M_L=5.3$]. Specifically, an excerpt of the MHz EME recorded by our remote telemetric station near the city of Ioannina five days before the EQ was found to possess critical characteristics, followed by an excerpt possessing tricritical characteristics. The analysis was performed by means of the method of critical fluctuations (MCF) revealing critical.


title: Confinement from Center Vortices: A review of old and new results

abstract: I briefly review the numerical evidence, some old and some quite recent, in favor of the center vortex theory of confinement.


title: Arc-transitive digraphs with quasiprimitive local actions

abstract: Let $\Gamma$ be a finite $G$-vertex-transitive digraph. The in-local action of $(\Gamma,G)$ is the permutation group $L_-$ induced by the vertex-stabiliser on the set of in-neighbours of $v$. The out-local action $L_+$ is defined analogously. Note that $L_-$ and $L_+$ may not be isomorphic. We thus consider the problem of determining which pairs $(L_-,L_+)$ are possible. We prove some general results, but pay special attention to the case when $L_-$ and $L_+$ are both quasiprimitive. (Recall that a permutation group is quasiprimitive if each of its nontrivial normal subgroups is transitive.) Along the way, we prove a structural result about pairs of finite quasiprimitive groups of the same degree, one being (abstractly) isomorphic to a proper quotient of the other.


title: Constraints on the Early Terrestrial Surface UV Environment Relevant to Prebiotic Chemistry

abstract: The UV environment is a key boundary condition for the origin of life. However, considerable uncertainty exists as to planetary conditions and hence surface UV at abiogenesis. Here, we present two-stream multi-layer clear-sky calculations of the UV surface radiance on Earth at 3.9 Ga to constrain the UV surface fluence as a function of albedo, solar zenith angle (SZA), and atmospheric composition. Variation in albedo and latitude (through SZA) can affect maximum photoreaction rates by a factor of >10.4; for the same atmosphere, photoreactions can proceed an order of magnitude faster at the equator of a snowball Earth than at the poles of a warmer world. Surface conditions are important considerations when computing prebiotic UV fluences. For climatically reasonable levels of CO2, fluence shortward of 189 nm is screened out, meaning that prebiotic chemistry is robustly shielded from variations in UV fluence due to solar flares or variability. Strong shielding from CO2 also means that the UV surface fluence is insensitive to plausible levels of CH4, O2, and O3. At scattering wavelengths, UV fluence drops off comparatively slowly with increasing CO2 levels. However, if SO2 and/or H2S can build up to the 1-100 ppm level as hypothesized by some workers, then they can dramatically suppress surface fluence and hence prebiotic photoprocesses. H2O is a robust UV shield for <198 nm. This means that regardless of the levels of other atmospheric gases, fluence <198 nm is only available for cold, dry atmospheres, meaning sources with emission <198 nm (e.g. ArF eximer lasers) can only be used in simulations of cold environments with low abundance of volcanogenic gases. On the other hand, fluence at 254 nm is unshielded by H2O and is available across a broad range of CO2 columns, meaning that mercury lamps are suitable for initial studies regardless of the uncertainty in primordial H2O and CO2 levels.


title: Quantum-mechanical calculations of cross sections for electron collisions with atoms and molecules

abstract: An overview of quantum-mechanical methods to generate cross-section data for electron collisions with atoms and molecules is presented. Particular emphasis is placed on the time-independent close-coupling approach, since it is particularly suitable for low-energy collisions and also allows for systematic improvements as well as uncertainty estimates. The basic ideas are illustrated with examples for electron collisions with argon atoms and methane. For many atomic systems, such as e-Ar collisions, highly reliable cross sections can now be computed with quantified uncertainties. On the other hand, while electron collision calculations with molecules do provide key input data for plasma models, the methods and computer codes presently used require further development to make these inputs robust.


title: Relaxation of a Simulated Lipid Bilayer Vesicle Compressed by an AFM

abstract: Using Coarse-Grained Molecular Dynamics simulations, we study the relaxation of bilayer vesicles, uniaxially compressed by an Atomic Force Microscope (AFM) cantilever. The relaxation time exhibits a strong force-dependence. Force-compression curves are very similar to recent experiments wherein giant unilamellar vesicles were compressed in a nearly identical manner.


title: At the Frontier of Precision QCD in the LHC Era

abstract: This thesis discusses recent advances in precision calculations of quantum chromodynamics and their application to the Large Hadron Collider (LHC) physics program and beyond.


title: Cross-Lingual Syntactic Transfer with Limited Resources

abstract: We describe a simple but effective method for cross-lingual syntactic transfer of dependency parsers, in the scenario where a large amount of translation data is not available. The method makes use of three steps: 1) a method for deriving cross-lingual word clusters, which can then be used in a multilingual parser; 2) a method for transferring lexical information from a target language to source language treebanks; 3) a method for integrating these steps with the density-driven annotation projection method of Rasooli and Collins (2015). Experiments show improvements over the state-of-the-art in several languages used in previous work, in a setting where the only source of translation data is the Bible, a considerably smaller corpus than the Europarl corpus used in previous work. Results using the Europarl corpus as a source of translation data show additional improvements over the results of Rasooli and Collins (2015). We conclude with results on 38 datasets from the Universal Dependencies corpora.


title: Observation of a non-Abelian Yang Monopole: From New Chern Numbers to a Topological Transition

abstract: Because global topological properties are robust against local perturbations, understanding and manipulating the topological properties of physical systems is essential in advancing quantum science and technology. For quantum computation, topologically protected qubit operations can increase computational robustness, and for metrology the quantized Hall effect directly defines the von Klitzing constant. Fundamentally, topological order is generated by singularities called topological defects in extended spaces, and is quantified in terms of Chern numbers, each of which measures different sorts of fields traversing surfaces enclosing these topological singularities. Here, inspired by high energy theories, we describe our synthesis and characterization of a singularity present in non-Abelian gauge theories - a Yang monopole - using atomic Bose-Einstein condensates in a five-dimensional space, and quantify the monopole in terms of Chern numbers measured on enclosing manifolds. While the well-known 1st Chern number vanished, the 2nd Chern number, measured for the first time in any physical settings, did not. By displacing the manifold, we then observed a phase transition from "topological" to "trivial" as the monopole left the manifold.


title: The DRAT format and DRAT-trim checker

abstract: This document describes the DRAT format for clausal proofs and the DRAT-trim proof checker.


title: Tetrahedral Order in Liquid Crystals

abstract: We review the impact of tetrahedral order on the macroscopic dynamics of bent-core liquid crystals. We discuss tetrahedral order comparing with other types of orientational order, like nematic, polar nematic, polar smectic, and active polar order. In particular, we present hydrodynamic equations for phases, where only tetrahedral order exists or tetrahedral order is combined with nematic order. Among the latter we discriminate between three cases, where the nematic director (a) orients along a 4-fold, or (b) along a 3-fold symmetry axis of the tetrahedral structure. For the optically isotropic T_d phase, which only has tetrahedral order, we focus on the coupling of flow with e.g. temperature gradients and on the specific orientation behavior in external electric fields. For the transition to the nematic phase, electric fields lead to a temperature shift that is linear in the field strength. Electric fields induce nematic order, again linear in the field strength. If strong enough, electric fields can change the tetrahedral structure and symmetry leading to a polar phase. We briefly deal with the T phase that arises when tetrahedral order occurs in a system of chiral molecules. To case (a) belong (i) the non-polar, achiral, optically uniaxial D2d phase with ambidextrous helicity and with orientational frustration in external fields, (ii) the non-polar tetragonal S4 phase, (iii) the non-polar, orthorhombic D2 phase that is structurally chiral featuring ambidextrous chirality, (iv) the polar orthorhombic C2v phase, and (v) the polar monoclinic C2 phase. Case (b) results in a trigonal C3v phase that behaves like a biaxial polar nematic phase. Finally we discuss some experiments that show typical effects related to the existence of tetrahedral order.


title: Measurements of charged-particle distributions with the ATLAS detector

abstract: Inclusive charged-particle measurements probe the low-energy region of the non-perturbative quantum chromodynamics. The ATLAS collaboration has recently measured the charged-particle multiplicity and its dependence on transverse momentum and pseudorapidity in special data sets with low LHC beam currents, recorded at centre-of-mass energies of 8 TeV and 13 TeV. The measurements at 8 TeV cover a wide spectrum using charged-particle selections with minimum transverse momentum of both 100 MeV and 500 MeV and in various phase space regions of low and high charged-particle multiplicities, some of which are studied for the first time by ATLAS. The measurements at 13 TeV also present detailed studies with a minimum transverse momentum of both 100 MeV and 500 MeV. The measurements are compared with predictions of various tuned Monte Carlo generators and are found to provide strong constraints on these. None of the Monte Carlo generators with their respective tunes are able to reproduce all the features of the data.


title: Rotors in triangles and tethrahedra

abstract: A polytope $P$ is circumscribed about a convex body $\Phi\subset \mathbb{R}^n$ if $\Phi\subset P$ and each facet of $P$ is contained in a support hyperplane of $\Phi$. We say that a convex body $\Phi\subset \mathbb{R}^n$ is a rotor of a polytope $P$ if for each rotation $\rho$ of $\mathbb{R}^n$ there exist a translation $\tau$ so that $P$ is circumscribed about $\tau\rho\Phi$. In this paper we shall prove that if $P$ is a triangle, then there is a baricentric formula that describes the curvature of bd$\Phi$ at the contact points, $\{A_1, A_2,A_3\}$. We prove also that if $\Phi\subset \mathbb{R}^3$ is a convex body which is a rotor in a tetrahedron $T$ and if $\Phi$ intersects the faces of $T$ at the points $\{x_1, \dots, x_4\}$, then the normal lines of $\Phi$ at the contact points with $T$, $\{x_1, \dots, x_4\}$ generically belong to one ruling of a quadric surface.


title: Novel Electron-Phonon Relaxation Pathway in Graphite Revealed by Time-Resolved Raman Scattering and Angle-Resolved Photoemission Spectroscopy

abstract: Time dynamics of photoexcited electron-hole pairs is important for a number of technologies, in particular solar cells. We combined ultrafast pump-probe Raman scattering and photoemission to directly follow electron-hole excitations as well as the G-phonon in graphite after an excitation by an intense laser pulse. This phonon is known to couple relatively strongly to electrons. Cross-correlating effective electronic and phonon temperatures places new constraints on model-based fits. The accepted two-temperature model predicts that G-phonon population should start to increase as soon as excited electron-hole pairs are created and that the rate of increase should not depend strongly on the pump fluence. Instead we found that the increase of the G-phonon population occurs with a delay of $\sim$65 fs. This time-delay is also evidenced by the absence of the so-called self-pumping for G phonons. It decreases with increased pump fluence. We show that these observations imply a new relaxation pathway: Instead of hot carriers transferring energy to G-phonons directly, the energy is first transferred to optical phonons near the zone boundary K-points, which then decay into G-phonons via phonon-phonon scattering. Our work demonstrates that phonon-phonon interactions must be included in any calculations of hot carrier relaxation in optical absorbers even when only short timescales are considered.


title: Gamma-Ray Emission from Supernova Remnants and Surrounding Molecular Clouds

abstract: Galactic cosmic rays are believed to be accelerated at supernova remnant shocks. Gamma-ray observations of both supernova remnants and associated molecular clouds have been used in several occasions to test (so far quite successfully) this popular hypothesis. Despite that, a conclusive solution to the problem of cosmic ray origin is still missing, and further observational and theoretical efforts are needed. In this paper, the current status of these investigations is briefly reviewed.


title: Enhancing ICA Performance by Exploiting Sparsity: Application to FMRI Analysis

abstract: Independent component analysis (ICA) is a powerful method for blind source separation based on the assumption that sources are statistically independent. Though ICA has proven useful and has been employed in many applications, complete statistical independence can be too restrictive an assumption in practice. Additionally, important prior information about the data, such as sparsity, is usually available. Sparsity is a natural property of the data, a form of diversity, which, if incorporated into the ICA model, can relax the independence assumption, resulting in an improvement in the overall separation performance. In this work, we propose a new variant of ICA by entropy bound minimization (ICA-EBM)-a flexible, yet parameter-free algorithm-through the direct exploitation of sparsity. Using this new SparseICA-EBM algorithm, we study the synergy of independence and sparsity through simulations on synthetic as well as functional magnetic resonance imaging (fMRI)-like data.


title: Holographic Theory of Gravity and Cosmology

abstract: According to the holographic principle, the maximum amount of information stored in a region of space scales as the area of its two-dimensional surface, like a hologram. We show that the holographic principle can be understood heuristically as originated from quantum fluctuations of spacetime. Applied to cosmology, this consideration leads to a dynamical cosmological constant $\Lambda$ of the observed magnitude, in agreement with the result obtained for the present and recent cosmic eras, by using unimodular gravity and causal-set theory. By generalizing the concept of entropic gravity, we find a critical acceleration parameter related to $\Lambda$ in galactic dynamics, and we construct a phenomenological model of dark matter which we call "modified dark matter" (MDM). We provide successful observational tests of MDM at both the galactic and cluster scales. We also discuss the possibility that the quanta of both dark energy and dark matter obey the quantum Boltzmann statistics or infinite statistics as described by a curious average of the bosonic and fermionic algebras.


title: The Spread of Cooperative Strategies on Grids with Random Asynchronous Updating

abstract: The Prisoner's Dilemma Process on a graph $G$ is an iterative process where each vertex, with a fixed strategy (cooperate or defect), plays the game with each of its neighbours. At the end of a round each vertex may change its strategy to that of its neighbour with the highest pay-off. Here we study the spread of cooperative and selfish behaviours on a toroidal grid, where each vertex is initially a cooperator with probability $p$. When vertices are permitted to change their strategies via a randomized asynchronous update scheme, we find that for some values of $p$ the limiting density of cooperators may be modelled as a polynomial in $p$. Theoretical bounds for this density are confirmed via simulation.


title: Many-body localization and delocalization from the perspective of Integrals of Motion

abstract: We study many-body localization (MBL) and delocalization from the perspective of integrals of motion (IOMs). MBL can be understood phenomenologically through the existence of macroscopically many localized IOMs. However, IOMs exist for all many-body systems, and non-localized IOMs determine properties on the ergodic side of the MBL transition too. Here we explore their properties using our method of displacement transformations. We show how different quantities can be calculated using the IOMs as an expansion in the number of operators. For all values of disorder the typical IOMs are localized, suggesting the importance of rare fluctuations in understanding the delocalization transition.


title: Relative weights approach to dynamical fermions at finite densities

abstract: The method of relative weights, coupled with mean field theory, is applied to the problem of simulating gauge theories with dynamical staggered fermions at finite densities. We present initial results and discuss issues so far encountered.


title: On the effects of alternative optima in context-specific metabolic model predictions

abstract: Recent methodological developments have facilitated the integration of high-throughput data into genome-scale models to obtain context-specific metabolic reconstructions. A unique solution to this data integration problem often may not be guaranteed, leading to a multitude of context-specific predictions equally concordant with the integrated data. Yet, little attention has been paid to the alternative optima resulting from the integration of context-specific data. Here we present computational approaches to analyze alternative optima for different context-specific data integration instances. By using these approaches on metabolic reconstructions for the leaf of Arabidopsis thaliana and the human liver, we show that the analysis of alternative optima is key to adequately evaluating the specificity of the predictions in particular cellular contexts. While we provide several ways to reduce the ambiguity in the context-specific predictions, our findings indicate that the existence of alternative optimal solutions warrant caution in detailed context-specific analyses of metabolism.


title: Integral operator approach over octonions to solution of nonlinear PDE

abstract: Integration of nonlinear partial differential equations with the help of the non-commutative integration over octonions is studied. An apparatus permitting to take into account symmetry properties of PDOs is developed. For this purpose formulas for calculations of commutators of integral and partial differential operators are deduced. Transformations of partial differential operators and solutions of partial differential equations are investigated. Theorems providing solutions of nonlinear PDEs are proved. Examples are given. Applications to PDEs of hydrodynamics and other types PDEs are described.


title: Finding Online Extremists in Social Networks

abstract: Online extremists in social networks pose a new form of threat to the general public. These extremists range from cyberbullies who harass innocent users to terrorist organizations such as the Islamic State of Iraq and Syria (ISIS) that use social networks to recruit and incite violence. Currently social networks suspend the accounts of such extremists in response to user complaints. The challenge is that these extremist users simply create new accounts and continue their activities. In this work we present a new set of operational capabilities to deal with the threat posed by online extremists in social networks. Using data from several hundred thousand extremist accounts on Twitter, we develop a behavioral model for these users, in particular what their accounts look like and who they connect with. This model is used to identify new extremist accounts by predicting if they will be suspended for extremist activity. We also use this model to track existing extremist users as they create new accounts by identifying if two accounts belong to the same user. Finally, we present a model for searching the social network to efficiently find suspended users' new accounts based on a variant of the classic Polya's urn setup. We find a simple characterization of the optimal search policy for this model under fairly general conditions. Our urn model and main theoretical results generalize easily to search problems in other fields.


title: Detection of faint broad emission lines in type 2 AGN: I. Near infrared observations and spectral fitting

abstract: We present medium resolution near infrared spectroscopic observations of 41 obscured and intermediate class AGN (type 2, 1.9 and 1.8; AGN2) with redshift $z \lesssim$0.1, selected from the Swift/BAT 70-month catalogue. The observations have been carried out in the framework of a systematic study of the AGN2 near infrared spectral properties and have been executed using ISAAC/VLT, X-shooter/VLT and LUCI/LBT, reaching an average S/N ratio of $\sim$30 per resolution element. For those objects observed with X-shooter we also obtained simultaneous optical and UV spectroscopy. We have identified a component from the broad line region in 13 out of 41 AGN2, with FWHM ${\rm > 800 }$ km/s. We have verified that the detection of the broad line region components does not significantly depend on selection effects due to the quality of the spectra, the X-ray or near infrared fluxes, the orientation angle of the host galaxy or the hydrogen column density measured in the X-ray band. The average broad line region components found in AGN2 has a significantly (a factor 2) smaller FWHM if compared with a control sample of type 1 AGN.


title: Measuring neutrino-induced exclusive charge-current final states on hydrogen at T2K

abstract: By taking advantage of symmetries with respect to the plane containing the directions of the neutrino and outgoing lepton, it is possible to isolate neutrino interactions on hydrogen in composite nuclear targets. This technique enables us to study the `primary' neutrino-nucleon interaction and therefore gain access to fundamental model parameters free from nuclear effects. Using T2K Monte Carlo equivalent to $\sim7\times10^{21}$ POT, we present an update on the measurement of the exclusive charged-current $\mu^-$, p, $\pi^+$ final state on hydrogen.


title: Commensurability for certain right-angled Coxeter groups and geometric amalgams of free groups

abstract: We give explicit necessary and sufficient conditions for the abstract commensurability of certain families of 1-ended, hyperbolic groups, namely right-angled Coxeter groups defined by generalized theta-graphs and cycles of generalized theta-graphs, and geometric amalgams of free groups whose JSJ graphs are trees of diameter at most 4. We also show that if a geometric amalgam of free groups has JSJ graph a tree, then it is commensurable to a right-angled Coxeter group, and give an example of a geometric amalgam of free groups which is not quasi-isometric (hence not commensurable) to any group which is finitely generated by torsion elements. Our proofs involve a new geometric realization of the right-angled Coxeter groups we consider, such that covers corresponding to torsion-free, finite-index subgroups are surface amalgams.


title: The Hodge theory of the Hecke category

abstract: Ideas from Hodge theory have found important applications in representation theory. We give a survey of joint work with Ben Elias which uncovers Hodge theoretic structure in the Hecke category ("Soergel bimodules"). We also outline similarities and differences to other combinatorial Hodge theories.


title: Charge transfer at hybrid inorganic-organic interfaces

abstract: Organic dopants are frequently used to surface-dope inorganic semiconductors. The resulted hybrid inorganic-organic materials have a crucial role in advanced functional materials and semiconductor devices. In this article, we study charge transfer at hybrid silicon-molecule interfaces theoretically. The idea is to filter out the best molecular acceptors to dope silicon with hole densities as high as 1013 cm-2. Here, we use a combinatorial algorithm merging chemical hardness method and first-principles DFT and GW calculations. We start by using the chemical hardness method which is simple and fast to narrow down our search for molecular dopants. Then, for the most optimistic candidates, we perform first-principles DFT calculations and discuss the necessity of GW corrections. This screening approach is quite general and applicable to other hybrid interfaces.


title: Pairing between zeros and critical points of random polynomials with independent roots

abstract: Let $p_n$ be a random, degree $n$ polynomial whose roots are chosen independently according to the probability measure $\mu$ on the complex plane. For a deterministic point $\xi$ lying outside the support of $\mu$, we show that almost surely the polynomial $q_n(z):=p_n(z)(z - \xi)$ has a critical point at distance $O(1/n)$ from $\xi$. In other words, conditioning the random polynomials $p_n$ to have a root at $\xi$, almost surely forces a critical point near $\xi$. More generally, we prove an analogous result for the critical points of $q_n(z):=p_n(z)(z - \xi_1)\cdots (z - \xi_k)$, where $\xi_1, \ldots, \xi_k$ are deterministic. In addition, when $k=o(n)$, we show that the empirical distribution constructed from the critical points of $q_n$ converges to $\mu$ in probability as the degree tends to infinity, extending a recent result of Kabluchko.


title: Multilevel Anomaly Detection for Mixed Data

abstract: Anomalies are those deviating from the norm. Unsupervised anomaly detection often translates to identifying low density regions. Major problems arise when data is high-dimensional and mixed of discrete and continuous attributes. We propose MIXMAD, which stands for MIXed data Multilevel Anomaly Detection, an ensemble method that estimates the sparse regions across multiple levels of abstraction of mixed data. The hypothesis is for domains where multiple data abstractions exist, a data point may be anomalous with respect to the raw representation or more abstract representations. To this end, our method sequentially constructs an ensemble of Deep Belief Nets (DBNs) with varying depths. Each DBN is an energy-based detector at a predefined abstraction level. At the bottom level of each DBN, there is a Mixed-variate Restricted Boltzmann Machine that models the density of mixed data. Predictions across the ensemble are finally combined via rank aggregation. The proposed MIXMAD is evaluated on high-dimensional realworld datasets of different characteristics. The results demonstrate that for anomaly detection, (a) multilevel abstraction of high-dimensional and mixed data is a sensible strategy, and (b) empirically, MIXMAD is superior to popular unsupervised detection methods for both homogeneous and mixed data.


title: Quark Helicity and Transversity Distributions

abstract: The quark transversity distribution inside nucleon is less understood than the quark unpolarized and helicity distributions inside nucleon. In particular, it is important to know clearly why the quark helicity and transversity distributions are different. We investigate the origin of their discrepancy.


title: DeepGraph: Graph Structure Predicts Network Growth

abstract: The topological (or graph) structures of real-world networks are known to be predictive of multiple dynamic properties of the networks. Conventionally, a graph structure is represented using an adjacency matrix or a set of hand-crafted structural features. These representations either fail to highlight local and global properties of the graph or suffer from a severe loss of structural information. There lacks an effective graph representation, which hinges the realization of the predictive power of network structures. In this study, we propose to learn the represention of a graph, or the topological structure of a network, through a deep learning model. This end-to-end prediction model, named DeepGraph, takes the input of the raw adjacency matrix of a real-world network and outputs a prediction of the growth of the network. The adjacency matrix is first represented using a graph descriptor based on the heat kernel signature, which is then passed through a multi-column, multi-resolution convolutional neural network. Extensive experiments on five large collections of real-world networks demonstrate that the proposed prediction model significantly improves the effectiveness of existing methods, including linear or nonlinear regressors that use hand-crafted features, graph kernels, and competing deep learning methods.


title: Gd$^{3+}$ - Gd$^{3+}$ distances exceeding 3 nm determined by very high frequency continuous wave electron paramagnetic resonance

abstract: Electron paramagnetic resonance spectroscopy in combination with site-directed spin-labeling is a very powerful tool for elucidating the structure and organization of biomolecules. Gd$^{3+}$ complexes have recently emerged as a new class of spin labels for distance determination by pulsed EPR spectroscopy at Q- and W-band. We present CW EPR measurements at 240 GHz (8.6 Tesla) on a series of Gd-rulers of the type Gd-PyMTA---spacer---Gd-PyMTA, with Gd-Gd distances ranging from 1.2 nm to 4.3 nm. CW EPR measurements of these Gd-rulers show that significant dipolar broadening of the central $|-1/2\rangle\rightarrow|1/2\rangle$ transition occurs at 30 K for Gd-Gd distances up to $\sim$ 3.4 nm with Gd-PyMTA as the spin label. This represents a significant extension for distances accessible by CW EPR, as nitroxide-based spin labels at X-band frequencies can typically only access distances up to $\sim$ 2 nm. We show that this broadening persists at biologically relevant temperatures above 200 K, and that this method is further extendable up to room temperature by immobilizing the sample in glassy trehalose. We show that the peak-to-peak broadening of the central transition follows the expected $~1/r^3$ dependence for the electron-electron dipolar interaction, from cryogenic temperatures up to room temperature. A simple procedure for simulating the dependence of the lineshape on interspin distance is presented, in which the broadening of the central transition is modeled as an $S=1/2$ spin whose CW EPR lineshape is broadened through electron-electron dipolar interactions with a neighboring $S=7/2$ spin.


title: First-principles analysis of electron transport in BaSnO$_3$

abstract: BaSnO$_3$ (BSO) is a promising transparent conducting oxide (TCO) with reported room-temperature (RT) Hall mobility exceeding 320 cm$^{2}$V$^{-1}$s$^{-1}$. Among perovskite oxides, it has the highest RT mobility, about 30 times higher than that of the prototypical SrTiO$_3$. Using first-principles calculations based on hybrid density functional theory, we elucidate the physical mechanisms that govern the mobility by studying the details of LO-phonon and ionized impurity scattering. A careful numerical analysis to obtain converged results within the relaxation-time approximation of Boltzmann transport theory is presented. The ${\bf k}$ dependence of the relaxation time is fully taken into account. We find that the high RT mobility in BSO originates not only from a small effective mass, but also from a significant reduction in the phonon scattering rate compared to other perovskite oxides; the origins of this reduction are identified. Ionized impurity scattering influences the total mobility even at RT for dopant densities larger than $5\times10^{18}$ cm$^{-3}$, and becomes comparable to LO-phonon scattering for $1\times10^{20}$ cm$^{-3}$ doping, reducing the drift mobility from its intrinsic LO-phonon-limited value of $\sim$594 cm$^{2}$V$^{-1}$s$^{-1}$ to less than 310 cm$^{2}$V$^{-1}$s$^{-1}$. We suggest pathways to avoid impurity scattering via modulation doping or polar discontinuity doping. We also explicitly calculate the Hall factor and Hall mobility, allowing a direct comparison to experimental reports for bulk and thin films and providing insights into the nature of the dominant mechanisms that limit mobility in state-of-the art samples.


title: QCD and $\gamma\,\gamma$ studies at FCC-ee

abstract: The Future Circular Collider (FCC) is a post-LHC project aiming at searches for physics beyond the SM in a new 80--100~km tunnel at CERN. Running in its first phase as a very-high-luminosity electron-positron collider (FCC-ee), it will provide unique possibilities for indirect searches of new phenomena through high-precision tests of the SM. In addition, by collecting tens of ab$^{-1}$ integrated luminosity in the range of center-of-mass energies $\sqrt{s}$~=90--350~GeV, the FCC-ee also offers unique physics opportunities for precise measurements of QCD phenomena and of photon-photon collisions through, literally, billions of hadronic final states as well as unprecedented large fluxes of quasireal $\gamma$'s radiated from the $\rm e^+e^-$ beams. We succinctly summarize the FCC-ee perspectives for high-precision extractions of the QCD coupling, for detailed analyses of parton radiation and fragmentation, and for SM and BSM studies through $\gamma\gamma$ collisions.


title: Witten indices of abelian M5 brane on $\mathbb{R}\times S^5$

abstract: Witten indices and partition functions are computed for abelian 6d tensor and hypermultiplets on $\mathbb{R}\times S^5$ in Lorentzian signature in an R gauge field background which preserves some supersymmetry. We consider a generic supersymmetric squashing that also admits squashing of the Hopf fiber. Wick rotation to Euclidean M5 brane amounts to Wick rotation of squashing parameters and the hypermultiplet mass parameter. We compute Casimir energies for tensor and hypermultiplets separately for general squashing, and match these with the corresponding gravitational anomaly polynomials. We extract Witten indices on $\mathbb{R}\times \mathbb{CP}^2$ and find that this is zero, again matching with the vanishing anomaly polynomial on an odd dimensional space.


title: Distinguishing the nonjet azimuth quadrupole from QCD jets and hydrodynamic flows via 2D angular correlations and quadrupole spectrum analysis

abstract: According to the flow narrative commonly applied to high-energy nuclear collisions a 1D cylindrical-quadrupole component of 2D angular correlations conventionally denoted by quantity $v_2$ is interpreted to represent elliptic flow: azimuth modulation of transverse or radial flow in noncentral nucleus-nucleus (A-A) collisions. The nonjet (NJ) quadrupole component exhibits various properties inconsistent with a flow or hydro interpretation, including the observation that NJ-quadrupole centrality variation in $A$-$A$ collisions has no relation to strongly-varying jet modification ("jet quenching") in those collisions commonly attributed to jet interaction with a dense flowing medium. In the present study I report isolation of quadrupole spectra from $p_t$-differential $v_2(p_t)$ data obtained at the relativistic heavy ion collider (RHIC) and large hadron collider (LHCr). I demonstrate that NJ quadrupole spectra have characteristics very different from the single-particle spectra for most hadrons, that quadrupole spectra indicate a common boosted hadron source for a small minority of hadrons that "carry" the quadrupole structure, that the narrow source-boost distribution is characteristic of an expanding thin cylindrical shell (also strongly contradicting a hydro interpretation), and that in the boost frame a single universal quadrupole spectrum (L\'evy distribution) on transverse mass $m_t$ accurately describes data for several hadron species scaled according to their statistical-model abundances. The quadrupole spectrum shape changes very little from RHIC to LHC energies. Taken in combination those characteristics strongly suggest a unique {\em nonflow} (and nonjet) QCD mechanism for the NJ quadrupole conventionally represented by $v_2$.


title: Robust quantum state transfer between two superconducting qubits via partial measurement

abstract: We develop a potentially practical proposal for robust quantum state transfer (QST) between two superconducting qubits coupled by a coplanar waveguide (CPW) resonator. We show that the partial measurement could drastically enhance the fidelity even when the dissipation of qubits and CPW is considered. Unlike many other schemes for QST, our proposal does not require the couplings between the qubits and the CPW resonator to be strong. In fact, our method works much more efficiently in the weak coupling regime. The underlying mechanism is attributed to the probabilistic nature of partial measurement.


title: Using Fast Weights to Attend to the Recent Past

abstract: Until recently, research on artificial neural networks was largely restricted to systems with only two types of variable: Neural activities that represent the current or recent input and weights that learn to capture regularities among inputs, outputs and payoffs. There is no good reason for this restriction. Synapses have dynamics at many different time-scales and this suggests that artificial neural networks might benefit from variables that change slower than activities but much faster than the standard weights. These "fast weights" can be used to store temporary memories of the recent past and they provide a neurally plausible way of implementing the type of attention to the past that has recently proved very helpful in sequence-to-sequence models. By using fast weights we can avoid the need to store copies of neural activity patterns.


title: Understanding experimentally-observed fluctuations

abstract: We discuss two topics on the experimental measurements of fluctuation observables in relativistic heavy-ion collisions. First, we discuss the effects of the thermal blurring, i.e. the blurring effect arising from the experimental measurement of fluctuations in momentum space as a proxy of the thermal fluctuations defined in coordinate space, on higher order cumulants. Second, we discuss the effect of imperfect efficiency of detectors on the measurement of higher order cumulants. We derive effective formulas which can carry out the correction of this effect for higher order cumulants based on the binomial model.


title: Detecting macroscopic quantum coherence with a cavity optomechanical system

abstract: The rigorous resource framework of quantum coherence has been set up recently and excited a wide variety of interests. Here we show that a quantum cavity optomechanical system, as an emerging platform, can behave with a certain value of quantum coherence at a macroscopic scale. We also find that the difference between the total optomechanical coherence and the sum of the optical and the mechanical coherence just equals their mutual information. Motivated by the detection of the optomechanical entanglement, an experimentally feasible scheme to probe the optomechanical coherence is proposed.


title: Algebraic representations and constructible sheaves

abstract: These are notes for my Takagi lecture at the University of Tokyo in November, 2016. I survey what is known about simple modules for reductive algebraic groups. The emphasis is on characteristic p>0 and Lusztig's character formula. I explain ideas connecting representations and constructible sheaves (Finkelberg-Mirkovic conjecture) in the spirit of the Kazhdan-Lusztig conjecture. I also discuss a conjecture with S. Riche (a theorem for GL_n) which should eventually make computations more feasible.


title: There are asymptotically the same number of Latin squares of each parity

abstract: A Latin square is reduced if its first row and column are in natural order. For Latin squares of a particular order $n$ there are four possible different parities. We confirm a conjecture of Stones and Wanless by showing asymptotic equality between the numbers of reduced Latin squares of each possible parity as the order $n\rightarrow\infty$.


title: Cousin groups and Hodge structures

abstract: We consider the geometric properties of Hodge Cousin groups, introduced in an unpublished paper \cite{OVV}, emphasizing the case of Hodge Cousin groups corresponding to polarized $\mathbb{Q}$-Hodge structures. Basing on this consideration, we introduce the class of abelian Cousin groups and prove an analogue of Poincar\'e complete reducibility theorem for them.


title: Foundations of Modern Query Languages for Graph Databases

abstract: We survey foundational features underlying modern graph query languages. We first discuss two popular graph data models: edge-labelled graphs, where nodes are connected by directed, labelled edges; and property graphs, where nodes and edges can further have attributes. Next we discuss the two most fundamental graph querying functionalities: graph patterns and navigational expressions. We start with graph patterns, in which a graph-structured query is matched against the data. Thereafter we discuss navigational expressions, in which patterns can be matched recursively against the graph to navigate paths of arbitrary length; we give an overview of what kinds of expressions have been proposed, and how they can be combined with graph patterns. We also discuss several semantics under which queries using the previous features can be evaluated, what effects the selection of features and semantics has on complexity, and offer examples of such features in three modern languages that are used to query graphs: SPARQL, Cypher and Gremlin. We conclude by discussing the importance of formalisation for graph query languages; a summary of what is known about SPARQL, Cypher and Gremlin in terms of expressivity and complexity; and an outline of possible future directions for the area.


title: Molecular Beam Epitaxy of 2D-layered Gallium Selenide on GaN substrates

abstract: Large area epitaxy of two-dimensional (2D) layered materials with high material quality is a crucial step in realizing novel device applications based on 2D materials. In this work, we report high-quality, crystalline, large-area gallium selenide (GaSe) films grown on bulk substrates such as c-plane sapphire and gallium nitride (GaN) using a valved cracker source for Se. (002)-oriented GaSe with random in-plane orientation of domains was grown on sapphire and GaN substrates at a substrate temperature of 350-450 C with complete surface coverage and smooth surface morphology. Higher growth temperature (575 C) resulted in the formation of single-crystalline {\epsilon}-GaSe triangular domains with six-fold symmetry confirmed by in-situ reflection high electron energy diffraction (RHEED) and off-axis x-ray diffraction (XRD). A two-step growth method involving high temperature nucleation of single crystalline domains and low temperature growth to enhance coalescence was adopted to obtain continuous (002)-oriented GaSe with an epitaxial relationship with the substrate. While six-fold symmetry was maintained in the two step growth, {\beta}-GaSe phase was observed in addition to the dominant {\epsilon}-GaSe in cross-sectional scanning transmission electron microscopy images. This work demonstrates the potential of growing high quality 2D-layered materials using molecular beam epitaxy and can be extended to the growth of other transition metal chalcogenides.


title: Adaptive Substring Extraction and Modified Local NBNN Scoring for Binary Feature-based Local Mobile Visual Search without False Positives

abstract: In this paper, we propose a stand-alone mobile visual search system based on binary features and the bag-of-visual words framework. The contribution of this study is three-fold: (1) We propose an adaptive substring extraction method that adaptively extracts informative bits from the original binary vector and stores them in the inverted index. These substrings are used to refine visual word-based matching. (2) A modified local NBNN scoring method is proposed in the context of image retrieval, which considers the density of binary features in scoring each feature matching. (3) In order to suppress false positives, we introduce a convexity check step that imposes a convexity constraint on the configuration of a transformed reference image. The proposed system improves retrieval accuracy by 11% compared with a conventional method without increasing the database size. Furthermore, our system with the convexity check does not lead to false positive results.


title: Auxiliary-field quantum Monte Carlo calculations with multiple-projector pseudopotentials

abstract: We have implemented recently developed multiple-projector pseudopotentials into the planewave based auxiliary-field quantum Monte Carlo (pw-AFQMC) method. Multiple-projector pseudopotentials can yield smaller planewave cut-offs while maintaining or improving transferability. This reduces the computational cost of pw-AFQMC, increasing its reach to larger and more complicated systems. We discuss the use of non-local pseudopotentials in the separable Kleinman-Bylander form, and the implementation in pw-AFQMC of the multiple-projector optimized norm-conserving pseudopotential ONCVPSP of Hamann. The accuracy of the method is first demonstrated by equation-of-state calculations of the ionic insulator NaCl and more strongly correlated metal Cu. The method is then applied to calibrate the accuracy of density functional theory (DFT) predictions of the phase stability of recently discovered high temperature and pressure superconducting sulfur hydride systems. We find that DFT results are in good agreement with pw-AFQMC, due to near cancellation of electron-electron correlation effects between different structures.


title: Online Training of an Opto-Electronic Reservoir Computer Applied to Real-Time Channel Equalisation

abstract: Reservoir Computing is a bio-inspired computing paradigm for processing time dependent signals. The performance of its analogue implementation are comparable to other state of the art algorithms for tasks such as speech recognition or chaotic time series prediction, but these are often constrained by the offline training methods commonly employed. Here we investigated the online learning approach by training an opto-electronic reservoir computer using a simple gradient descent algorithm, programmed on an FPGA chip. Our system was applied to wireless communications, a quickly growing domain with an increasing demand for fast analogue devices to equalise the nonlinear distorted channels. We report error rates up to two orders of magnitude lower than previous implementations on this task. We show that our system is particularly well-suited for realistic channel equalisation by testing it on a drifting and a switching channels and obtaining good performances


title: Embodiment of Learning in Electro-Optical Signal Processors

abstract: Delay-coupled electro-optical systems have received much attention for their dynamical properties and their potential use in signal processing. In particular it has recently been demonstrated, using the artificial intelligence algorithm known as reservoir computing, that photonic implementations of such systems solve complex tasks such as speech recognition. Here we show how the backpropagation algorithm can be physically implemented on the same electro-optical delay-coupled architecture used for computation with only minor changes to the original design. We find that, compared when the backpropagation algorithm is not used, the error rate of the resulting computing device, evaluated on three benchmark tasks, decreases considerably. This demonstrates that electro-optical analog computers can embody a large part of their own training process, allowing them to be applied to new, more difficult tasks.


title: Safeguarding Decentralized Wireless Networks Using Full-Duplex Jamming Receivers

abstract: In this paper, we study the benefits of full-duplex (FD) receiver jamming in enhancing the physical-layer security of a two-tier decentralized wireless network with each tier deployed with a large number of pairs of a single-antenna transmitter and a multi-antenna receiver. In the underlying tier, the transmitter sends unclassified information, and the receiver works in the halfduplex (HD) mode receiving the desired signal. In the overlaid tier, the transmitter deliveries confidential information in the presence of randomly located eavesdroppers, and the receiver works in the FD mode radiating jamming signals to confuse eavesdroppers and receiving the desired signal simultaneously. We provide a comprehensive performance analysis and network design under a stochastic geometry framework. Specifically, we consider the scenarios where each FD receiver uses single- and multi-antenna jamming, and analyze the connection probability and the secrecy outage probability of a typical FD receiver by providing accurate expressions and more tractable approximations for the two metrics. We further determine the optimal deployment of the FD-mode tier in order to maximize networkwide secrecy throughput subject to constraints including the given dual probabilities and the network-wide throughput of the HD-mode tier. Numerical results are demonstrated to verify our theoretical findings, and show that network-wide secrecy throughput is significantly improved by properly deploying the FD-mode tier.


title: Dual atomic interferometer with a tunable point of minimum magnetic sensitivity

abstract: Atomic interferometers are often affected by magnetic field fluctuations. Using the clock transition at zero magnetic field minimizes the effect of these fluctuations. There is another transition in rubidium that minimizes the magnetic sensitivity at 3.2 G.We combine the previous two transitions to obtain minimum magnetic sensitivity at a tunable magnetic field between 2.2 and 3.2 G. The two interferometers evolve independently from each other and we control the magnetic sensitivity by changing the population in both transitions with a microwave pulse.


title: Lexicon Integrated CNN Models with Attention for Sentiment Analysis

abstract: With the advent of word embeddings, lexicons are no longer fully utilized for sentiment analysis although they still provide important features in the traditional setting. This paper introduces a novel approach to sentiment analysis that integrates lexicon embeddings and an attention mechanism into Convolutional Neural Networks. Our approach performs separate convolutions for word and lexicon embeddings and provides a global view of the document using attention. Our models are experimented on both the SemEval'16 Task 4 dataset and the Stanford Sentiment Treebank, and show comparative or better results against the existing state-of-the-art systems. Our analysis shows that lexicon embeddings allow to build high-performing models with much smaller word embeddings, and the attention mechanism effectively dims out noisy words for sentiment analysis.


title: Prototype Filter Design for FBMC in Massive MIMO Channels

abstract: We perform an asymptotic study on the performance of filter bank multicarrier (FBMC) in the context of massive multi-input multi-output (MIMO). We show that the signal-to-interference-plus-noise ratio (SINR) cannot grow unboundedly by increasing the number of base station (BS) antennas, and is upper bounded by a certain deterministic value. This is a result of the correlation between the multi-antenna combining tap values and the channel impulse responses between the terminals and the BS antennas. To solve this problem, we introduce a simple FBMC prototype filter design method that removes this correlation, enabling us to achieve arbitrarily large SINR values by increasing the number of BS antennas.


title: Design of Provably Physical-Constraint-Preserving Methods for General Relativistic Hydrodynamics

abstract: The paper develops high-order physical-constraint-preserving (PCP) methods for general relativistic hydrodynamic (GRHD) equations, equipped with a general equation of state. Here the physical constraints, describing the admissible states of GRHD, are referred to the subluminal constraint on the fluid velocity and the positivity of the density, pressure and specific internal energy. Preserving these constraints is very important for robust computations, otherwise violating one of them will lead to the ill-posed problem and numerical instability. To overcome the difficulties arising from the inherent strong nonlinearity contained in the constraints, we derive an equivalent definition of the admissible states. Using this definition, we prove the convexity, scaling invariance and Lax-Friedrichs (LxF) splitting property of the admissible state set $\mathcal G$, and discover the dependence of $\mathcal G$ on the spacetime metric. Unfortunately, such dependence yields the non-equivalence of $\mathcal G$ at different points in curved spacetime, and invalidates the convexity of $\mathcal G$ in analyzing PCP schemes. This obstacle is effectively overcame by introducing a new formulation of the GRHD equations. Based on this formulation and the above theories, a first-order LxF scheme is designed on general unstructured mesh and rigorously proved to be PCP under a CFL condition. With two types of PCP limiting procedures, we design high-order, {\em provably} (not probably) PCP methods under discretization on the proposed new formulation. These high-order methods include the PCP finite difference, finite volume and discontinuous Galerkin methods.


title: Comment on "Anomalous Edge State in a Non-Hermitian Lattice"

abstract: In this comment, we criticize three main conclusions of the letter\cite{Lee2016}. We show that the concept of fractional winding number(FWN) is factitious, Lee's conclusions on Fig. 3 are finite-size effect and the breakdown of bulk-boundary correspondence (BBBC) cannot be explained by "defective".


title: Modeling Scalability of Distributed Machine Learning

abstract: Present day machine learning is computationally intensive and processes large amounts of data. It is implemented in a distributed fashion in order to address these scalability issues. The work is parallelized across a number of computing nodes. It is usually hard to estimate in advance how many nodes to use for a particular workload. We propose a simple framework for estimating the scalability of distributed machine learning algorithms. We measure the scalability by means of the speedup an algorithm achieves with more nodes. We propose time complexity models for gradient descent and graphical model inference. We validate our models with experiments on deep learning training and belief propagation. This framework was used to study the scalability of machine learning algorithms in Apache Spark.


title: The near-surface methane humidity on Titan

abstract: We retrieve vertical and meridional variations of methane mole fraction in Titan's lower troposphere by re-analyzing near-infrared ground-based observations from 17 July 2014 UT (Adamkovics et al., 2016). We generate synthetic spectra using atmospheric methane profiles that do not contain supersaturation or discontinuities to fit the observations, and thereby retrieve minimum saturation altitudes and corresponding specific humidities in the boundary layer. We relate these in turn to surface-level relative humidities using independent surface temperature measurements. We also compare our results with general circulation model simulations to interpret and constrain the relationship between humidities and surface liquids. The results show that Titan's lower troposphere is undersaturated at latitudes south of 60N, consistent with a dry surface there, but increases in humidity toward the north pole indicate appreciable surface liquid coverage. While our observations are consistent with considerably more liquid methane existing at the north pole than is present in observed lakes, a degeneracy between low-level methane and haze leads to substantial uncertainty in determining the extent of the source region.


title: Arc-like continua, Julia sets of entire functions, and Eremenko's Conjecture

abstract: A hyperbolic transcendental entire function with connected Fatou set is said to be "of disjoint type". It is known that a disjoint-type function provides a model for the dynamics near infinity of all maps in the same parameter space; hence a good understanding of these functions has implications in wider generality. Our goal is to study the topological properties of the Julia sets of entire functions of disjoint type. In particular, we give a detailed description of the topology of their connected components. More precisely, consider a "Julia continuum" C of such a function, i.e. the closure in the Riemann sphere of a component of the Julia set. We show that infinity is a terminal point of C, and that C has span zero in the sense of Lelek; under a mild geometric assumption on the function C is arc-like. (Whether every span zero continuum is also arc-like was a famous question in continuum theory, only recently resolved in the negative.) Conversely, we construct a single disjoint-type entire function with the remarkable property that each arc-like continuum with at least one terminal point is realised as a Julia continuum. The class of arc-like continua with terminal points is uncountable. It includes, in particular, the sin(1/x)-curve, the Knaster buckethandle and the pseudo-arc, so these can all occur as Julia continua of a disjoint-type entire function. We give similar descriptions of the possible topology of Julia continua that contain periodic points or points with bounded orbits, and answer a question of Bara\'nski and Karpi\'nska by showing that Julia continua need not contain points that are accessible from the Fatou set. Furthermore, we construct an entire function whose Julia set has connected components on which the iterates tend to infinity pointwise, but not uniformly. This is related to a famous conjecture of Eremenko concerning escaping sets of entire functions.


title: Linear Process Bootstrap Unit Root Test

abstract: One of the most widely applied unit root test, Phillips-Perron test, enjoys in general highpowers, but suffers from size distortions when moving average noise exists. As a remedy, thispaper proposes a nonparametric bootstrap unit root test that specifically targets moving aver-age noise. Via a bootstrap functional central limit theorem, the consistency of this bootstrapapproach is established under general assumptions which allows a large family of non-linear timeseries. In simulation, this bootstrap test alleviates the size distortions of the Phillips-Perrontest while preserving its high powers.


title: Semiclassical asymptotics of quantum weighted Hurwitz numbers

abstract: This work concerns the semiclassical asymptotics of quantum weighted double Hurwitz numbers. We compute the leading term of the partition function for three versions of the quantum weighted Hurwitz numbers, as well as lower order semiclassical corrections. The classical limit $\hbar \to 0$ is shown to reproduce the simple Hurwitz numbers studied by Pandharipande and Okounkov. The KP-Toda $\tau$-function serving as generating function for the quantum Hurwitz numbers is shown to converge in the classical limit to the generating function of Pandharipande and Okounkov and, with suitable scaling, so do the partition function, the weights and expectations of Hurwitz numbers.


title: Effects of Substrate Heating and Wettability on Evaporation Dynamics and Deposition Patterns for a Sessile Water Droplet Containing Colloidal Particles

abstract: Effects of substrate temperature, substrate wettability and particles concentration are experimentally investigated for evaporation of a sessile water droplet containing colloidal particles. Time-varying droplet shapes and temperature of the liquid-gas interface are measured using high-speed visualization and infrared thermography, respectively. The motion of the particles inside the evaporating droplet is qualitatively visualized by an optical microscope and profile of final particle deposit is measured by an optical profilometer. On a non-heated hydrophilic substrate, a ring-like deposit forms after the evaporation, as reported extensively in the literature; while on a heated hydrophilic substrate, a thinner ring with an inner deposit is reported in the present work. The latter is attributed to Marangoni convection and recorded motion of the particles as well as measured temperature gradient across the liquid-gas interface confirms this hypothesis. The thinning of the ring scales with the substrate temperature and is reasoned to stronger Marangoni convection at larger substrate temperature. In case of a non-heated hydrophobic substrate, an inner deposit forms due to very early depinning of the contact line. On the other hand, in case of a heated hydrophobic substrate, the substrate heating as well as larger particle concentration helps in the pinning of the contact line, which results in a thin ring with an inner deposit. We propose a regime map for predicting three types of deposits namely, ring, thin ring with inner deposit and inner deposit - for varying substrate temperature, substrate wettability and particles concentration. A first-order model corroborates the liquid-gas interface temperature measurements and variation in the measured ring profile with the substrate temperature.