title: Mapping Tractography Across Subjects

abstract: Diffusion magnetic resonance imaging (dMRI) and tractography provide means to study the anatomical structures within the white matter of the brain. When studying tractography data across subjects, it is usually necessary to align, i.e. to register, tractographies together. This registration step is most often performed by applying the transformation resulting from the registration of other volumetric images (T1, FA). In contrast with registration methods that "transform" tractographies, in this work, we try to find which streamline in one tractography correspond to which streamline in the other tractography, without any transformation. In other words, we try to find a "mapping" between the tractographies. We propose a graph-based solution for the tractography mapping problem and we explain similarities and differences with the related well-known graph matching problem. Specifically, we define a loss function based on the pairwise streamline distance and reformulate the mapping problem as combinatorial optimization of that loss function. We show preliminary promising results where we compare the proposed method, implemented with simulated annealing, against a standard registration techniques in a task of segmentation of the corticospinal tract.


title: Proximal-gradient algorithms for fractional programming

abstract: In this paper we propose two proximal gradient algorithms for fractional programming problems in real Hilbert spaces, where the numerator is a proper, convex and lower semicontinuous function and the denominator is a smooth function, either concave or convex. In the iterative schemes, we perform a proximal step with respect to the nonsmooth numerator and a gradient step with respect to the smooth denominator. The algorithm in case of a concave denominator has the particularity that it generates sequences which approach both the (global) optimal solutions set and the optimal objective value of the underlying fractional programming problem. In case of a convex denominator the numerical scheme approaches the set of critical points of the objective function, provided the latter satisfies the Kurdyka-\L{}ojasiewicz property.


title: Turbulent reconnection of magnetic bipoles in stratified turbulence

abstract: We consider strongly stratified forced turbulence in a plane-parallel layer with helicity and corresponding large-scale dynamo action in the lower part and non-helical turbulence in the upper. The magnetic field is found to develop strongly concentrated bipolar structures near the surface. They form elongated bands with a sharp interface between opposite polarities. Unlike earlier experiments with imposed magnetic field, the inclusion of rotation does not strongly suppress the formation of these structures. We perform a systematic numerical study of this phenomenon by varying magnetic Reynolds number, scale separation ratio, and Coriolis number. We focus on the formation of a current sheet between bipolar regions where reconnection of oppositely oriented field lines occurs. We determine the reconnection rate by measuring either the inflow velocity in the vicinity of the current sheet or by measuring the electric field in the reconnection region. We demonstrate that for large Lundquist numbers, S>10^3, the reconnection rate is nearly independent of S in agreement with results of recent numerical simulations performed by other groups in simpler settings.


title: Lower-Vietoris-type Topologies on Hyperspaces

abstract: We introduce a new lower-Vietoris-type hypertopology in a way similar to that with which a new upper-Vietoris-type hypertopology was introduced in G. Dimov and D. Vakarelov, "On Scott consequence systems", Fundamenta Informaticae, 33 (1998), 43-70. (it was called there {\em Tychonoff-type hypertopology}). We study this new hypertopology and, in particular, we generalize many results from E. Cuchillo-Ibanez, M. A. Moron and F. R. Ruiz del Portal, "Lower semifinite topology in hyperspaces", Topology Proceedings, 17 (1992), 29-39. As a corollary, we get that for every continuous map $f:X\longrightarrow X$, where $X$ is a continuum, there exist a subcontinuum $K$ of $X$ such that $f(K)=K.$


title: Kernels for sequentially ordered data

abstract: We present a novel framework for kernel learning with sequential data of any kind, such as time series, sequences of graphs, or strings. Our approach is based on signature features which can be seen as an ordered variant of sample (cross-)moments; it allows to obtain a "sequentialized" version of any static kernel. The sequential kernels are efficiently computable for discrete sequences and are shown to approximate a continuous moment form in a sampling sense. A number of known kernels for sequences arise as "sequentializations" of suitable static kernels: string kernels may be obtained as a special case, and alignment kernels are closely related up to a modification that resolves their open non-definiteness issue. Our experiments indicate that our signature-based sequential kernel framework may be a promising approach to learning with sequential data, such as time series, that allows to avoid extensive manual pre-processing.

categories: stat.ML cs.DM cs.LG math.ST stat.ME stat.TH


title: Hard core repulsion and supersolid cluster crystals

abstract: We study the effect of a short-ranged hard-core repulsion on the stability and superfluid properties of the cluster crystal phase of two-dimensional (2D) soft core bosons. Results of Quantum Monte Carlo simulations on a cogent test case suggest that the main physical properties of the phase remain unaltered if the range d of the inner repulsive core is sufficiently short, even if the strength of the repulsion is several orders of magnitude greater than the outer soft core barrier. Only if d is an appreciable fraction of the size of the clusters (> 5%) does a sufficiently strong hard core repulsion cause the crystal to break down into a homogeneous superfluid; a moderate inner core repulsion enhances the superfluid response of the crystalline phase.


title: Mass transfer between debris discs during close stellar encounters

abstract: We study mass transfers between debris discs during stellar encounters. We carried out numerical simulations of close flybys of two stars, one of which has a disc of planetesimals represented by test particles. We explored the parameter space of the encounters, varying the mass ratio of the two stars, their pericentre and eccentricity of the encounter, and its geometry. We find that particles are transferred to the other star from a restricted radial range in the disc and the limiting radii of this transfer region depend on the parameters of the encounter. We derive an approximate analytic description of the inner radius of the region. The efficiency of the mass transfer generally decreases with increasing encounter pericentre and increasing mass of the star initially possessing the disc. Depending on the parameters of the encounter, the transfer particles have a specific distributions in the space of orbital elements (semimajor axis, eccentricity, inclination, and argument of pericentre) around their new host star. The population of the transferred particles can be used to constrain the encounter through which it was delivered. We expect that many stars experienced transfer among their debris discs and planetary systems in their birth environment. This mechanism presents a formation channel for objects on wide orbits of arbitrary inclinations, typically having high eccentricity but possibly also close-to-circular (eccentricities of about 0.1). Depending on the geometry, such orbital elements can be distinct from those of the objects formed around the star.


title: Isometries of nilpotent metric groups

abstract: We consider Lie groups equipped with arbitrary distances. We only assume that the distance is left-invariant and induces the manifold topology. For brevity, we call such object metric Lie groups. Apart from Riemannian Lie groups, distinguished examples are sub-Riemannian Lie groups and, in particular, Carnot groups equipped with Carnot-Carath\'eodory distances. We study the regularity of isometries, i.e., distance-preserving homeomorphisms. Our first result is the analyticity of such maps between metric Lie groups. The second result is that if two metric Lie groups are connected and nilpotent then every isometry between the groups is the composition of a left translation and an isomorphism. There are counterexamples if one does not assume the groups to be either connected or nilpotent. The first result is based on a solution of the Hilbert 5th problem by Montgomery and Zippin. The second result is proved, via the first result, considering the Riemannian case, which for self-isometries was solved by Wolf.


title: On the Vinogradov mean value

abstract: A discussion of recent work of C.Demeter, L.Guth and the author of the proof of the Vinogradov Main Conjecture using the decoupling theory for curves


title: On Multi-step MLE-process for Markov Sequences

abstract: We consider the problem of the construction of the estimator-process of the unknown finite-dimensional parameter in the case of the observations of nonlinear autoregressive process. The estimation is done in two or three steps. First we estimate the unknown parameter by a learning relatively short part of observations and then we use the one-step MLE idea to construct an-estimator process which is asymptotically equivalent to the MLE. To have the learning interval shorter we introduce the two-step procedure which leads to the asymptotically efficient estimator-process too. The presented results are illustrated with the help of two numerical examples.


title: A simple nonsingular inflationary quintessential model

abstract: In the framework of a flat Friedmann-Lema{\^\i}tre-Robertson-Walker (FLRW) geometry, we present a nonsingular model (no big bang singularity at finite time) of our universe describing its evolution starting from its early inflationary era up to the present accelerating phase. We found that a hydrodynamical fluid with nonlinear equation of state could result in such a non singular scenario, which after the end of this inflationary stage, suffers a sudden phase transition and enters into the stiff matter dominated era, and the universe becomes reheated due to a huge amount of particle production. Finally, it asymptotically enters into the de Sitter phase concluding the present accelerated expansion. We show that the background providing an inflationary potential leads to a power spectrum of the cosmological perturbations which fit well with the latest Planck estimations. At the end, we compared our viable potential with some known inflationary quintessential potentials, which shows that, our model is an improved version of them because it contains an analytic solution that allows us to perform analytic calculations.


title: Weakly Mutually Uncorrelated Codes

abstract: We introduce the notion of weakly mutually uncorrelated (WMU) sequences, motivated by applications in DNA-based storage systems and synchronization protocols. WMU sequences are characterized by the property that no sufficiently long suffix of one sequence is the prefix of the same or another sequence. In addition, WMU sequences used in DNA-based storage systems are required to have balanced compositions of symbols and to be at large mutual Hamming distance from each other. We present a number of constructions for balanced, error-correcting WMU codes using Dyck paths, Knuth's balancing principle, prefix synchronized and cyclic codes.


title: On Asanov's Finsleroid-Finsler metrics as the solutions of a conformal rigidity problem

abstract: Finsleroid-Finsler metrics form an important class of singular (y-local) Finsler metrics. They were introduced by G. S. Asanov [2] in 2006. As the special case of the general construction Asanov produced singular (y - local) examples of Landsberg spaces of dimension at least three that are not of Berwald type. The existence of regular (y - global) Landsberg metrics that are not of Berwald type is an open problem up to this day; for a detailed exposition of the so-called unicorn problem in Finsler geometry see D. Bao [3]. In this paper we are going to characterize the Finsleroid-Finsler metrics as the solutions of a conformal rigidity problem. We are looking for (non-Riemannian) Finsler metrics admitting a (non-homothetic) conformal change such that the mixed curvature tensor of the Berwald connection contracted by the derivatives of the logarithmic scale function is invariant. We prove that the solutions of class at least $\mathcal{C}^2$ on the complement of the zero section are conformal to Finsleroid-Finsler metrics.


title: The orthogonality of Al-Salam-Carlitz polynomials for complex parameters

abstract: In this contribution, we study the orthogonality conditions satisfied by Al-Salam-Carlitz polynomials $U^{(a)}_n(x;q)$ when the parameters $a$ and $q$ are not necessarily real nor `classical', i.e., the linear functional $\bf u$ with respect to such polynomial sequence is quasi-definite and not positive definite. We establish orthogonality on a simple contour in the complex plane which depends on the parameters. In all cases we show that the orthogonality conditions characterize the Al-Salam-Carlitz polynomials $U_n^{(a)}(x;q)$ of degree $n$ up to a constant factor. We also obtain a generalization of the unique generating function for these polynomials.


title: Factorizing the factorization - a spectral-element solver for elliptic equations with linear operation count

abstract: High-order methods gain more and more attention in computational fluid dynamics. However, the potential advantage of these methods depends critically on the availability of efficient elliptic solvers. With spectral-element methods, static condensation is a common approach to reduce the number of degree of freedoms and to improve the condition of the algebraic equations. The resulting system is block-structured and the face-based operator well suited for matrix-matrix multiplications. However, a straight-forward implementation scales super-linearly with the number of unknowns and, therefore, prohibits the application to high polynomial degrees. This paper proposes a novel factorization technique, which yields a linear operation count of just 13N multiplications, where N is the total number of unknowns. In comparison to previous work it saves a factor larger than 3 and clearly outpaces unfactored variants for all polynomial degrees. Using the new technique as a building block for a preconditioned conjugate gradient method resulted in a runtime scaling linearly with N for polynomial degrees $2 \leq p \leq 32$ . Moreover the solver proved remarkably robust for aspect ratios up to 128.


title: The exponential map in non-commutative probability

abstract: The wrapping transformation $W$ is a homomorphism from the semigroup of probability measures on the real line, with the convolution operation, to the semigroup of probability measures on the circle, with the multiplicative convolution operation. We show that on a large class $\mathcal{L}$ of measures, $W$ also transforms the three non-commutative convolutions---free, Boolean, and monotone---to their multiplicative counterparts. Moreover, the restriction of $W$ to $\mathcal{L}$ preserves various qualitative properties of measures and triangular arrays. We use these facts to give short proofs of numerous known, and new, results about multiplicative convolutions.


title: Local mirror symmetry and the sunset Feynman integral

abstract: We study the sunset Feynman integral defined as the scalar two-point self-energy at two-loop order in a two dimensional space-time. We firstly compute the Feynman integral, for arbitrary internal masses, in terms of the regulator of a class in the motivic cohomology of a 1-parameter family of open elliptic curves. Using an Hodge theoretic (B-model) approach, we show that the integral is given by a sum of elliptic dilogarithms evaluated at the divisors determined by the punctures. Secondly we associate to the sunset elliptic curve a local non-compact Calabi-Yau 3-fold, obtained as a limit of elliptically fibered compact Calabi-Yau 3-folds. By considering the limiting mixed Hodge structure of the Batyrev dual A-model, we arrive at an expression for the sunset Feynman integral in terms of the local Gromov-Witten prepotential of the del Pezzo surface of degree 6. This expression is obtained by proving a strong form of local mirror symmetry which identifies this prepotential with the second regulator period of the motivic cohomology class.

categories: hep-th math-ph math.AG math.KT math.MP math.NT


title: Casimir entropy and internal energy of the objects in fluctuating scalar and electromagnetic fields

abstract: Casimir entropy is an important aspect of casimir effect.In this paper,we employ the path integral method to derive the total relation for casimir entropy and internal energy of arbitrary shaped objects in the presence of two,three and four dimensions scalar fields and electromagnetic field.We obtain the casimir entropy and internal energy of two nanoribbon immersed in scalar field and two nanospheres immersed in scalar field and electromagnetic field.The casmir entropy of two nanospheres immersed in the electromagnetic field in small interval of temperature variations,shown a different behavior.


title: Observation of Quantized Conductance in the Surface Bands of Bismuth Nanowires

abstract: We report the experimental observation of quantized conductance by Rashba bands at the surface of 50-nm Bi nanowires. With increasing magnetic fields along the wire axis, the wires exhibit a stepwise increase of conductance with as many as four distinct plateaus together with oscillatory thermopower. The observations can be accounted for by the increase of the number of propagating surface modes. The modes have very high mobility and are associated with Aharonov-Bohm interference around the perimeter consistently with theory of quasiballistic one-dimensional nanowires considering that Lorentz forces decouple the modes from scattering at the nanowire surface.


title: An eccentric binary millisecond pulsar with a helium white dwarf companion in the Galactic Field

abstract: Low-mass white dwarfs (LMWDs) are believed to be exclusive products of binary evolution, as the Universe is not yet old enough to produce them from single stars. Because of the strong tidal forces operating during the binary interaction phase, the remnant host systems observed today are expected to have negligible eccentricities. Here, we report on the first unambiguous identification of a LMWD in an eccentric (e=0.13) orbit with a millisecond pulsar, which directly contradicts this picture. We use our spectra and radio-timing solution (derived elsewhere) to infer the WD temperature T_eff = 8600 +/- 190 K) and 3D systemic velocity (179.5 km\s). We also place model-independent constraints on the WD radius (R_WD = 0.024+/- 0.004/0.002 R_sun) and surface gravity (log g = 7.11 +/- 0.08/0.16 dex). The WD and kinematic properties are consistent with the expectations for low-mass X-ray binary evolution and disfavour a three-body formation channel. In the case of the high eccentricity being the result of a spontaneous phase transition, we infer a mass of 1.6 M_sun for the progenitor of the pulsar, which is too low for the quark-nova mechanism proposed by Jiang et al. (2015). Similarly, the scenario of Freire & Tauris (2014), in which a WD collapses onto a neutron star via an rotationally-delayed accretion-induced collapse, requires both a high-mass differentially rotating progenitor and a significant momentum kick at birth under our constraints. Contrarily, we find that eccentricity pumping via interaction with a transient circumbinary disk is consistent with all inferred properties. Finally, we report tentative evidence for pulsations which, if confirmed, would transform the star into an unprecedented laboratory for WD physics and stellar convection.


title: Proof Lengths for Instances of the Paris-Harrington Principle

abstract: As Paris and Harrington have famously shown, Peano Arithmetic does not prove that for all numbers $k,m,n$ there is an $N$ which satisfies the statement $\operatorname{PH}(k,m,n,N)$: For any $k$-colouring of its $n$-element subsets the set $\{0,\dots,N-1\}$ has a large homogeneous subset of size $\geq m$. At the same time very weak theories can establish the $\Sigma_1$-statement $\exists_N\operatorname{PH}(\overline k,\overline m,\overline n,N)$ for any fixed parameters $k,m,n$. Which theory, then, does it take to formalize natural proofs of these instances? It is known that $\forall_m\exists_N\operatorname{PH}(\overline k,m,\overline n,N)$ has a natural and short proof (relative to $n$ and $k$) by $\Sigma_{n-1}$-induction. In contrast, we show that there is an elementary function $e$ such that any proof of $\exists_N\operatorname{PH}(\overline{e(n)},\overline{n+1},\overline n,N)$ by $\Sigma_{n-2}$-induction is ridiculously long. In order to establish this result on proof lengths we give a computational analysis of slow provability, a notion introduced by Sy-David Friedman, Rathjen and Weiermann. We will see that slow uniform $\Sigma_1$-reflection is related to a function that has a considerably lower growth rate than $F_{\varepsilon_0}$ but dominates all functions $F_\alpha$ with $\alpha<\varepsilon_0$ in the fast-growing hierarchy.


title: On the geometry of the LLSvS eightfold

abstract: In this note we make a few remarks about the geometry of the holomorphic symplectic manifold Z constructed by C.Lehn, M.Lehn, C.Sorger and D. van Straten as a two-step contraction of the variety of twisted cubic curves on a cubic fourfold Y in P^5. We show that Z is birational to a component of a moduli space of stable sheaves in the Calabi-Yau subcategory of the derived category of Y. Using this description we deduce that the twisted cubics contained in a hyperplane section Y_H of Y give rise to a Lagrangian subvariety Z_H in Z. For a generic choice of the hyperplane, Z_H is birational to the theta-divisor in the intermediate Jacobian of Y_H.


title: Carbon Phosphide Monolayer with Superior Carrier Mobility

abstract: Two dimensional (2D) materials with a finite band gap and high carrier mobility are sought after materials from both fundamental and technological perspectives. In this paper, we present the results based on the particle swarm optimization method and density functional theory which predict three geometrically different phases of carbon phosphide (CP) monolayer consisted of sp2 hybridized C atoms and sp3 hybridized P atoms in hexagonal networks. Two of the phases, referred to as {\alpha}-CP and \b{eta}-CP with puckered and buckled surfaces, respectively are semiconducting with highly anisotropic electronic and mechanical properties. More remarkably, they have lightest electrons and holes among the known 2D semiconductors, yielding superior carrier mobility. The {\gamma}-CP has a distorted hexagonal network and exhibits a semi-metallic behavior with Dirac cones. These theoretical findings suggest the binary CP monolayer to be yet unexplored 2D materials holding great promises for applications in high-performance electronics and optoelectronics.


title: Lipreading with Long Short-Term Memory

abstract: Lipreading, i.e. speech recognition from visual-only recordings of a speaker's face, can be achieved with a processing pipeline based solely on neural networks, yielding significantly better accuracy than conventional methods. Feed-forward and recurrent neural network layers (namely Long Short-Term Memory; LSTM) are stacked to form a single structure which is trained by back-propagating error gradients through all the layers. The performance of such a stacked network was experimentally evaluated and compared to a standard Support Vector Machine classifier using conventional computer vision features (Eigenlips and Histograms of Oriented Gradients). The evaluation was performed on data from 19 speakers of the publicly available GRID corpus. With 51 different words to classify, we report a best word accuracy on held-out evaluation speakers of 79.6% using the end-to-end neural network-based solution (11.6% improvement over the best feature-based solution evaluated).


title: The Right Way to Search Evolving Graphs

abstract: Evolving graphs arise in problems where interrelations between data change over time. We present a breadth first search (BFS) algorithm for evolving graphs that computes the most direct influences between nodes at two different times. Using simple examples, we show that naive unfoldings of adjacency matrices miscount the number of temporal paths. By mapping an evolving graph to an adjacency matrix of an equivalent static graph, we prove that our generalization of the BFS algorithm correctly accounts for paths that traverse both space and time. Finally, we demonstrate how the BFS over evolving graphs can be applied to mine citation networks.


title: On MBR codes with replication

abstract: An early paper by Rashmi et. al. presented the construction of an $(n,k,d=n-1)$ MBR regenerating code featuring the inherent double replication of all code symbols and repair-by-transfer (RBT), both of which are important in practice. We first show that no MBR code can contain even a single code symbol that is replicated more than twice. We then go on to present two new families of MBR codes which feature double replication of all systematic message symbols. The codes also possess a set of $d$ nodes whose contents include the message symbols and which can be repaired through help-by-transfer (HBT). As a corollary, we obtain systematic RBT codes for the case $d=(n-1)$ that possess inherent double replication of all code symbols and having a field size of $O(n)$ in comparison with the general, $O(n^2)$ field size requirement of the earlier construction by Rashmi et. al. For the cases $(k=d=n-2)$ or $(k+1=d=n-2)$, the field size can be reduced to $q=2$, and hence the codes can be binary. We also give a necessary and sufficient condition for the existence of MBR codes having double replication of all code symbols and also suggest techniques which will enable an arbitrary MBR code to be converted to one with double replication of all code symbols.


title: Local Sobolev Constant Estimate for Integral Ricci Curvature Bounds

abstract: We obtain a local Sobolev constant estimate for integral Ricci curvature, which enables us to extend several important tools such as the maximal principle, the gradient estimate, the heat kernel estimate and the $L^2$ Hessian estimate to manifolds with integral Ricci lower bounds, without the non-collapsing conditions.


title: Supergravity background of lambda-deformed model for AdS2 x S2 supercoset

abstract: Starting with the F/G supercoset model corresponding to the AdS_n x S^n superstring one can define the lambda-model of arXiv:1409.1538 either as a deformation of the F/F gauged WZW model or as an integrable one-parameter generalization of the non-abelian T-dual of the AdS_n x S^n superstring sigma model with respect to the whole supergroup F. Here we consider the case of n=2 and find the explicit form of the 4d target space background for the lambda-model for the PSU(1,1|2)/[SO(1,1) x SO(2)] supercoset. We show that this background represents a solution of type IIB 10d supergravity compactified on a 6-torus with only metric, dilaton Phi and the RR 5-form (represented by a 2-form F in 4d) being non-trivial. This implies that the lambda-model is Weyl invariant at the quantum level and thus defines a consistent superstring sigma model. The supergravity solution we find is different from the one in arXiv:1410.1886 which should correspond to a version of the lambda-model where only the bosonic subgroup of F is gauged. Still, the two solutions have equivalent scaling limit of arXiv:1504.07213 leading to the isometric background for the metric and e^Phi F which is related to the eta-deformed AdS_2 x S^2 sigma model of arXiv:1309.5850. Similar results are expected in the AdS_3 x S^3 and AdS_5 x S^5 cases.


title: Probing top quark neutral couplings in the Standard Model Effective Field Theory at NLO QCD

abstract: Top quark pair production in association with a $Z$-boson or a photon at the LHC directly probes neutral top-quark couplings. We present predictions for these two processes in the Standard Model (SM) Effective Field Theory (EFT) at next-to-leading (NLO) order in QCD. We include the full set of CP-even dimension-six operators that enter the top-quark interactions with the SM gauge bosons. For comparison, we also present predictions in the SMEFT for top loop-induced $HZ$ production at the LHC and for $t\bar{t}$ production at the ILC at NLO in QCD. Results for total cross sections and differential distributions are obtained and uncertainties coming from missing higher orders in the strong coupling and in the EFT expansions are discussed. NLO results matched to the parton shower are available, allowing for event generation to be directly employed in an experimental analyses. Our framework provides a solid basis for the interpretation of current and future measurements in the SMEFT, with improved accuracy and precision.


title: Ordered groupoid quotients and congruences on inverse semigroups

abstract: We introduce a preorder on an inverse semigroup $S$ associated to any normal inverse subsemigroup $N$, that lies between the natural partial order and Green's ${\mathscr J}$-relation. The corresponding equivalence relation $\simeq_N$ is not necessarily a congruence on $S$, but the quotient set does inherit a natural ordered groupoid structure. We show that this construction permits the factorisation of any inverse semigroup homomorphism into a composition of a quotient map and a star-injective functor, and that this decomposition implies a classification of congruences on $S$. We give an application to the congruence and certain normal inverse subsemigroups associate to an inverse monoid presentation.


title: Weak-coupling superconductivity in a strongly correlated iron pnictide

abstract: Iron-based superconductors have been found to exhibit an intimate interplay of orbital, spin, and lattice degrees of freedom, dramatically affecting their low-energy electronic properties, including superconductivity. Albeit the precise pairing mechanism remains unidentified, several candidate interactions have been suggested to mediate the superconducting pairing, both in the orbital and in the spin channel. Here, we employ optical spectroscopy (OS), angle-resolved photoemission spectroscopy (ARPES), ab initio band-structure, and Eliashberg calculations to show that nearly optimally doped NaFe$_{0.978}$Co$_{0.022}$As exhibits some of the strongest orbitally selective electronic correlations in the family of iron pnictides. Unexpectedly, we find that the mass enhancement of itinerant charge carriers in the strongly correlated band is dramatically reduced near the $\Gamma$ point and attribute this effect to orbital mixing induced by pronounced spin-orbit coupling. Embracing the true band structure allows us to describe all low-energy electronic properties obtained in our experiments with remarkable consistency and demonstrate that superconductivity in this material is rather weak and mediated by spin fluctuations.


title: Two-parameter quantum general linear supergroups

abstract: The universal R-matrix of two-parameter quantum general linear supergroups is computed explicitly based on the RTT realization of Faddeev--Reshetikhin--Takhtajan.


title: Sequential double cross-validation for assessment of added predictive ability in high-dimensional omic applications

abstract: Enriching existing predictive models with new biomolecular markers is an important task in the new multi-omic era. Clinical studies increasingly include new sets of omic measurements which may prove their added value in terms of predictive performance. We introduce a two-step approach for the assessment of the added predictive ability of omic predictors, based on sequential double cross-validation and regularized regression models. We propose several performance indices to summarize the two-stage prediction procedure and a permutation test to formally assess the added predictive value of a second omic set of predictors over a primary omic source. The performance of the test is investigated through simulations. We illustrate the new method through the systematic assessment and comparison of the performance of transcriptomics and metabolomics sources in the prediction of body mass index (BMI) using longitudinal data from the Dietary, Lifestyle, and Genetic determinants of Obesity and Metabolic syndrome (DILGOM) study, a population-based cohort from Finland.


title: Rotational directionality via symmetry-breaking in an electrostatic motor

abstract: We theoretically investigate how one can achieve a preferred rotational direction for the case of a simple electrostatic motor. The motor is composed by a rotor and two electronic reservoirs. Electronic islands on the rotor can exchange electrons with the reservoirs. An electrostatic field exerts a force on the occupied islands. The charge dynamics and the electrostatic field drive rotations of the rotor. Coupling to an environment lead to damping on the rotational degree of freedom. We use two different approaches to the charge dynamics in the electronic islands: hopping process and mean-field. The hopping process approach takes into account charge fluctuations, which can appear along Coulomb blockade effects in nanoscale systems. The mean-field approach neglects the charge fluctuations on the islands, which is typically suitable for larger systems. We show that for a system described by the mean-field equations one can in principle prepare initial conditions to obtain a desired rotational direction. In contrast, this is not possible in the stochastic description. However, for both cases one can achieve rotational directionality by changing the geometry of the rotor. By scanning the space formed by the relevant geometric parameters we find optimal geometries, while fixing the dissipation and driving parameters. Remarkably, in the hopping process approach perfect rotational directionality is possible for a large range of geometries.


title: On the toric ideals of matroids of a fixed rank

abstract: In $1980$ White conjectured that every element of the toric ideal of a matroid is generated by quadratic binomials corresponding to symmetric exchanges. We prove White's conjecture for high degrees with respect to the rank. This extends our result arXiv:1302.5236 confirming White's conjecture `up to saturation'. Furthermore, we study degrees of Gr\"{o}bner bases and Betti tables of the toric ideals of matroids of a fixed rank.


title: The radiative zone of the Sun and the tachocline: stability of baroclinic patterns of differential rotation

abstract: Barotropic rotation and radiative equilibrium are mutually incompatible in stars. The issue is often addressed by allowing for a meridional circulation, but this is not devoid of theoretical complications. Models of rotation in the Sun which maintain strict radiative equilibrium, making use of the observation that the Sun is not in a state of barotropic rotation, have recently been suggested. To investigate the dynamical behaviour of these solutions, we study the local stability of stratified, weakly magnetized, differentially rotating fluids to non-axisymmetric perturbations. Finite heat conductivity, kinematic viscosity, and resistivity are present. The evolution of local embedded perturbations is governed by a set of coupled, ordinary differential equations with time-dependent coefficients. Two baroclinic models of rotation for the upper radiative zone and tachocline are studied: (i) an interpolation based on helioseismology data, (ii) a theoretical solution directly compatible with radiative equilibrium. The growth of the local Goldreich-Schubert-Fricke instability appears to be suppressed, largely because of the viscosity. An extensive exploration of wavenumber space is carried out, with and without a magnetic field. Although we easily find classical local instabilities when they ought formally to be present, for the Sun the analysis reveals neither unstable solutions, nor even solutions featuring a large transient growth. We have not ruled out larger scale or nonlinear instabilities, nor have we rigorously proven local stability. But rotational configurations in close agreement with observations, generally thought to be vulnerable to the classic local Goldreich-Schubert-Fricke instability, do appear to be locally stable under rather general circumstances.


title: Spectrally Grouped Total Variation Reconstruction for Scatter Imaging Using ADMM

abstract: We consider X-ray coherent scatter imaging, where the goal is to reconstruct momentum transfer profiles (spectral distributions) at each spatial location from multiplexed measurements of scatter. Each material is characterized by a unique momentum transfer profile (MTP) which can be used to discriminate between different materials. We propose an iterative image reconstruction algorithm based on a Poisson noise model that can account for photon-limited measurements as well as various second order statistics of the data. To improve image quality, previous approaches use edge-preserving regularizers to promote piecewise constancy of the image in the spatial domain while treating each spectral bin separately. Instead, we propose spectrally grouped regularization that promotes piecewise constant images along the spatial directions but also ensures that the MTPs of neighboring spatial bins are similar, if they contain the same material. We demonstrate that this group regularization results in improvement of both spectral and spatial image quality. We pursue an optimization transfer approach where convex decompositions are used to lift the problem such that all hyper-voxels can be updated in parallel and in closed-form. The group penalty introduces a challenge since it is not directly amendable to these decompositions. We use the alternating directions method of multipliers (ADMM) to replace the original problem with an equivalent sequence of sub-problems that are amendable to convex decompositions, leading to a highly parallel algorithm. We demonstrate the performance on real data.


title: Identifying differences in the rules of interaction between individuals in moving animal groups

abstract: Collective movement can be achieved when individuals respond to the local movements and positions of their neighbours. Some individuals may disproportionately influence group movement if they occupy particular spatial positions in the group, for example, positions at the front of the group. We asked, therefore, what led individuals in moving pairs of fish (Gambusia holbrooki) to occupy a position in front of their partner. Individuals adjusted their speed and direction differently in response to their partner's position, resulting in individuals occupying different positions in the group. Individuals that were found most often at the front of the pair had greater mean changes in speed than their partner, and were less likely to turn towards their partner, compared to those individuals most often found at the back of the pair. The pair moved faster when led by the individual that was usually at the front. Our results highlight how differences in the social responsiveness between individuals can give rise to leadership in free moving groups. They also demonstrate how the movement characteristics of groups depend on the spatial configuration of individuals within them.


title: Emergent properties of nuclei from ab initio coupled-cluster calculations

abstract: Emergent properties such as nuclear saturation and deformation, and the effects on shell structure due to the proximity of the scattering continuum and particle decay channels are fascinating phenomena in atomic nuclei. In recent years, ab initio approaches to nuclei have taken the first steps towards tackling the computational challenge of describing these phenomena from Hamiltonians with microscopic degrees of freedom. This endeavor is now possible due to ideas from effective field theories, novel optimization strategies for nuclear interactions, ab initio methods exhibiting a soft scaling with mass number, and ever-increasing computational power. This paper reviews some of the recent accomplishments. We also present new results. The recently optimized chiral interaction NNLO$_{\rm sat}$ is shown to provide an accurate description of both charge radii and binding energies in selected light- and medium-mass nuclei up to $^{56}$Ni. We derive an efficient scheme for including continuum effects in coupled-cluster computations of nuclei based on chiral nucleon-nucleon and three-nucleon forces, and present new results for unbound states in the neutron-rich isotopes of oxygen and calcium. The coupling to the continuum impacts the energies of the $J^\pi = {1/2}^-,{3/2}^-,{7/2}^-,{3/2}^+$ states in $^{17,23,25}$O, and - contrary to naive shell-model expectations - the level ordering of the $J^\pi = {3/2}^+,{5/2}^+,{9/2}^+$ states in $^{53,55,61}$Ca.


title: Quantum Walks with Dynamical Control: Graph Engineering, Initial State Preparation and State Transfer

abstract: Quantum walks are a well-established model for the study of coherent transport phenomena and provide a universal platform in quantum information theory. Dynamically influencing the walker's evolution gives a high degree of flexibility for studying various applications. Here, we present time-multiplexed finite quantum walks of variable size, the preparation of non-localized input states and their dynamical evolution. As a further application, we implement a state transfer scheme for an arbitrary input state to two different output modes. The presented experiments rely on the full dynamical control of a time-multiplexed quantum walk, which includes adjustable coin operation as well as the possibility to flexibly configure the underlying graph structures.


title: Why state of quantum system is fully defined by density matrix

abstract: We show that probabilities of results of all possible measurements performing on a quantum system depend on the system's state only through its density matrix. Therefore all experimentally available information about the state contains in the density matrix. In this study, we do not postulate that measurements obey some given formalism (such as observables, positive-operator valued measures, etc.), and do not use Born rule. The process of measurement is considered in a fully operational manner---as an interaction of a measured system with some black-box apparatus. The key point of our approach is the proof that, for improper mixtures, the expected value of any measurement depends linearly on the reduced density function. Such a proof is achieved by considering appropriate thought experiments. We demonstrate that Born rule can be derived as a natural consequence of our results.


title: Expansion of polynomial Lie group integrals in terms of certain maps on surfaces, and factorizations of permutations

abstract: Using the diagrammatic approach to integrals over Gaussian random matrices, we find a representation for polynomial Lie group integrals as infinite sums over certain maps on surfaces. The maps involved satisfy a specific condition: they have some marked vertices, and no closed walks that avoid these vertices. We also formulate our results in terms of permutations, arriving at new kinds of factorization problems.


title: Active and Reactive Energy Balance Equations in Active and Reactive Time

abstract: Electrical networks, and physical systems in general, are known to satisfy a power balance equation which states that the rate of change of the energy in time equals the power at the port of the network minus the power dissipated. However, when complex power is considered, there does not seem to exist a similar statement for the imaginary power, either in the time-domain or the frequency-domain. Recently, in the context of electromagnetic fields, it has been shown by complexifying the time to t+js and interpreting s as reactive time, that it is possible to set up an imaginary power balance in terms of the rate of change of reactive energy in reactive time. Here these ideas are specialized to linear and time-invariant RLC networks. For non-sinusoidal waveforms it is shown that the rate of change of reactive energy in reactive time contains all the essential properties and features of the commonly accepted definition of reactive power under sinusoidal conditions. We believe that this provides an unambiguous and physically motivated resolution to the longstanding debate on how to generalize reactive power to non-sinusoidal waveforms.


title: Spectral triples for nested fractals

abstract: It is shown that, for nested fractals [T.Lindstrom, Mem. Amer. Math. Soc. 420, 1990], the main structural data, such as the Hausdorff dimension and measure, the geodesic distance (when it exists) induced by the immersion in $R^n$, and the self-similar energy can all be recovered by the description of the fractals in terms of the spectral triples considered in [D.Guido, T.Isola, in "Advances in Operator Algebras and Mathematical Physics", Theta Series in Advanced Mathematics, Bucharest 2005].


title: Evolution of segmented strings

abstract: I explain how to evolve segmented strings in de Sitter and anti-de Sitter spaces of any dimension in terms of forward-directed null displacements. The evolution is described entirely in terms of discrete hops which do not require a continuum spacetime. Moreover, the evolution rule is purely algebraic, so it can be defined not only on ordinary real de Sitter and anti-de Sitter, but also on the rational points of the quadratic equations that define these spaces. For three-dimensional anti-de Sitter space, a simpler evolution rule is possible that descends from the Wess-Zumino-Witten equations of motion. In this case, one may replace three-dimensional anti-de Sitter space by a non-compact discrete subgroup of SL(2,R) whose structure is related to the Pell equation. A discrete version of the BTZ black hole can be constructed as a quotient of this subgroup. This discrete black hole avoids the firewall paradox by a curious mechanism: even for large black holes, there are no points inside the horizon until one reaches the singularity.


title: Physics postgraduate teaching assistants' grading approaches: Conflicting goals and practices

abstract: Grading can shape students' learning and encourage use of effective problem solving practices. Teaching assistants (TAs) are often responsible for grading student solutions and providing feedback, thus, their perceptions of grading may impact grading practices in the physics classroom. Understanding TAs' perceptions of grading is instrumental for curriculum developers as well as professional development leaders interested in improving grading practices. In order to identify TAs' perceptions of grading, we used a data collection tool designed to elicit TAs' considerations when making grading decisions as well as elicit possible conflicts between their stated goals and actual grading practices. The tool was designed to explicate TAs' attitudes towards research-based grading practices that promote effective problem solving approaches. TAs were first asked to state their goals for grading in general. Then, TAs graded student solutions in a simulated setting while explicating and discussing their underlying considerations. The data collection tool was administered at the beginning of TAs' first postgraduate teaching appointment and again after one semester of teaching experience. We found that almost all of the TAs stated that the purpose of grading was formative, i.e., grading should encourage students to learn from their mistakes as well as inform the instructor of common student difficulties. However, when making grading decisions in a simulated setting, the majority of TAs' grading considerations focused on correctness and they did not assign grades in a way that encourages use of effective problem solving approaches. TAs' perceptions of grading did not change significantly during one semester of teaching experience.


title: Functional Renormalisation Group analysis of Tensorial Group Field Theories on $\mathbb{R}^d$

abstract: Rank-d Tensorial Group Field Theories are quantum field theories defined on a group manifold $G^{\times d}$, which represent a non-local generalization of standard QFT, and a candidate formalism for quantum gravity, since, when endowed with appropriate data, they can be interpreted as defining a field theoretic description of the fundamental building blocks of quantum spacetime. Their renormalisation analysis is crucial both for establishing their consistency as quantum field theories, and for studying the emergence of continuum spacetime and geometry from them. In this paper, we study the renormalisation group flow of two simple classes of TGFTs, defined for the group $G=\mathbb{R}$ for arbitrary rank, both without and with gauge invariance conditions, by means of functional renormalisation group techniques. The issue of IR divergences is tackled by the definition of a proper thermodynamic limit for TGFTs. We map the phase diagram of such models, in a simple truncation, and identify both UV and IR fixed points of the RG flow. Encouragingly, for all the models we study, we find evidence for the existence of a phase transition of condensation type.


title: Switching of chiral magnetic skyrmions by picosecond magnetic field pulses via transient topological states

abstract: Magnetic chiral skyrmions are vortex like spin structures that appear as stable or meta-stable states in magnetic materials due to the interplay between the symmetric and antisymmetric exchange interactions, applied magnetic field and/or uniaxial anisotropy. Their small size and internal stability make them prospective objects for data storage but for this, the controlled switching between skyrmion states of opposite polarity and topological charge is essential. Here we present a study of magnetic skyrmion switching by an applied magnetic field pulse based on a discrete model of classical spins and atomistic spin dynamics. We found a finite range of coupling parameters corresponding to the coexistence of two degenerate isolated skyrmions characterized by mutually inverted spin structures with opposite polarity and topological charge. We demonstrate how for a wide range of material parameters a short inclined magnetic field pulse can initiate the reliable switching between these states at GHz rates. Detailed analysis of the switching mechanism revealed the complex path of the system accompanied with the excitation of a chiral-achiral meron pair and the formation of an achiral skyrmion.


title: Renormalizable Tensor Field Theories

abstract: Extending tensor models at the field theoretical level, tensor field theories are nonlocal quantum field theories with Feynman graphs identified with simplicial complexes. They become relevant for addressing quantum topology and geometry in any dimension and therefore form an interesting class of models for studying quantum gravity. We review the class of perturbatively renormalizable tensor field theories and some of their features.


title: Slow Reflection

abstract: We describe a "slow" version of the hierarchy of uniform reflection principles over Peano Arithmetic ($\mathbf{PA}$). These principles are unprovable in Peano Arithmetic (even when extended by usual reflection principles of lower complexity) and introduce a new provably total function. At the same time the consistency of $\mathbf{PA}$ plus slow reflection is provable in $\mathbf{PA}+\operatorname{Con}(\mathbf{PA})$. We deduce a conjecture of S.-D. Friedman, Rathjen and Weiermann: Transfinite iterations of slow consistency generate a hierarchy of precisely $\varepsilon_0$ stages between $\mathbf{PA}$ and $\mathbf{PA}+\operatorname{Con}(\mathbf{PA})$ (where $\operatorname{Con}(\mathbf{PA})$ refers to the usual consistency statement).


title: Supersymmetric plasma systems and their nonsupersymmetric counterparts

abstract: In this thesis a systematic comparison of supersymmetric plasma systems and their nonsupersymmetric counterparts is presented. The work is motivated by the AdS/CFT correspondence and the main aim is to check how much the plasma governed by the N=4 super Yang-Mills theory resembles the quark-gluon plasma studied experimentally in relativistic heavy-ion collisions. The analysis is done in a weak coupling regime where perturbative methods are applicable. Since the Keldysh-Schwinger approach is used, not only equilibrium but also nonequilibrium plasmas, which are assumed to be ultrarelativistic, are under consideration. First, using the functional techniques we introduce Faddeev-Popov ghosts into the Keldysh-Schwinger formalism of nonAbelian gauge theories. Next the collective excitations of the N=1 SUSY QED plasma are considered and compared to those of the usual QED system. The analysis is repeated to confront with each other the plasmas governed by the N=4 super Yang-Mills and QCD theories. Finally, transport properties of the systems are also discussed. Since both collective excitations and transport characteristics of supersymmetric and nonsupersymmetric plasmas are very similar to each other we conclude that supersymmetric plasma systems are qualitatively the same as their nonsupersymmetric partners in weak coupling regime. The quantitative differences mostly reflect different numbers of degrees of freedom.


title: Very high energy Fermi/LAT detection of HESS J0632+057

abstract: We report on the results of ~7 yrs of the very-high energy (10-600 GeV) observations of HESS J0632+057 with FERMI/LAT. In the highest energy band, 200-600 GeV, the source is clearly detected with the statistical significance >3.6sigma at orbital phases 0.2-0.4 and 0.6-0.8 at which HESS J0632+057 is known to demonstrate enhanced emission in TeV energy band. The analysis did not reveal the emission from HESS J0632+057 at lower energies and different orbital phases. Using the upper limits on source's flux we locate the break of the spectrum to >80 GeV and low-energy slope <2.0 (2sigma statistical significance).


title: New calibrations for abundance determinations in HII regions

abstract: Simple relations for deriving the oxygen abundance in HII regions with intensities of the three strong emission lines R_2, R_3, and N_2 (R calibration) or S_2, R_3, and N_2 (S calibration) in their spectra are suggested. A sample of 313 reference HII regions of the counterpart method is used as calibrating data points. Relations for the determination of nitrogen abundances, the R calibration, are also constructed. We find that the oxygen and nitrogen abundances in high-metallicity HII regions can be estimated using the intensities of the two strong lines R_2 and N_2 (or S_2 and N_2 for oxygen) only. The corresponding two-dimensional relations are provided. There are considerable advantages of the suggested calibration relations as compared to the existing ones. First, the oxygen and nitrogen abundances estimated through the suggested calibrations agree with the Te-based abundances within ~0.1 dex over the whole metallicity range, i.e., the relative accuracy of the calibration-based abundances is 0.1 dex. Although we constructed distinct relations for high- and low-metallicity objects, the separation between these two can be simply obtained from the intensity of the N_2 line. Moreover, the applicability ranges of the high- and low-metallicity relations overlap for adjacent metallicities, i.e., the transition zone disappears. Second, the oxygen abundances produced by the two suggested calibrations are in remarkable agreement with each other. In fact, the R-based and S-based oxygen abundances agree within ~0.05 dex in the majority of cases for more than three thousand HII region spectra.


title: Gauge Theory by canonical Transformations

abstract: Electromagnetism, the strong and the weak interaction are commonly formulated as gauge theories in a Lagrangian description. In this paper we present an alternative formal derivation of U(1)-gauge theory in a manifestly covariant Hamilton formalism. We make use of canonical transformations as our guiding tool to formalize the gauging procedure. The introduction of the gauge field, its transformation behaviour and a dynamical gauge field Lagrangian/Hamiltonian are unavoidable consequences of this formalism, whereas the form of the free gauge Lagrangian/Hamiltonian depends on the selection of the gauge dependence of the canonically conjugate gauge fields.


title: Random walks in Dirichlet environment: an overview

abstract: Random Walks in Dirichlet Environment (RWDE) correspond to Random Walks in Random Environment (RWRE) on $\Bbb{Z}^d$ where the transition probabilities are i.i.d. at each site with a Dirichlet distribution. Hence, the model is parametrized by a family of positive weights $(\alpha_i)_{i=1, \ldots, 2d}$, one for each direction of $\Bbb{Z}^d$. In this case, the annealed law is that of a reinforced random walk, with linear reinforcement on directed edges. RWDE have a remarkable property of statistical invariance by time reversal from which can be inferred several properties that are still inaccessible for general environments, such as the equivalence of static and dynamic points of view and a description of the directionally transient and ballistic regimes. In this paper we give a state of the art on this model and several sketches of proofs presenting the core of the arguments. We also present new computation of the large deviation rate function for one dimensional RWDE.


title: Uncertainty quantification of effective nuclear interactions

abstract: We give a brief review on the development of phenomenological NN interactions and the corresponding quantification of statistical uncertainties. We look into the uncertainty of effective interactions broadly used in mean field calculations through the Skyrme parameters and effective field theory counter-terms by estimating both statistical and systematic uncertainties stemming from the NN interaction. We also comment on the role played by different fitting strategies on the light of recent developments.


title: WiSeDB: A Learning-based Workload Management Advisor for Cloud Databases

abstract: Workload management for cloud databases must deal with the tasks of resource provisioning, query placement and query scheduling in a manner that meets the application's performance goals while minimizing the cost of using cloud resources. Existing solutions have approached these three challenges in isolation, and with only a particular type of performance goal in mind. In this paper, we introduce WiSeDB, a learning-based framework for generating holistic workload management solutions customized to application-defined performance metrics and workload characteristics. Our approach relies on supervised learning to train cost-effective decision tree models for guiding query placement, scheduling, and resource provisioning decisions. Applications can use these models for both batch and online scheduling of incoming workloads. A unique feature of our system is that it can adapt its offline model to stricter/looser performance goals with minimal re-training. This allows us to present alternative workload management strategies that address the typical performance vs. cost trade-off of cloud services. Experimental results show that our approach has very low training overhead while offering low cost strategies for a variety of performance goals and workload characteristics.


title: Numerical Simulation of Tidal Evolution of a Viscoelastic Body Modelled with a Mass-Spring Network

abstract: We use a damped mass-spring model within an N-body code to simulate the tidal evolution of the spin and orbit of a self-gravitating viscoelastic spherical body moving around a point-mass perturber. The damped mass-spring model represents a Kelvin-Voigt viscoelastic solid. We measure the tidal quality function (the dynamical Love number $\,k_2\,$ divided by the tidal quality factor $\,Q\,$) from the numerically computed tidal drift of the semimajor axis of the binary. The shape of $\,k_2/Q\,$, as a function of the principal tidal frequency, reproduces the kink shape predicted by Efroimsky (2012a; CeMDA 112$\,:\,$283) for the tidal response of near-spherical homogeneous viscoelastic rotators. We demonstrate that we can directly simulate the tidal evolution of spinning viscoelastic objects. In future, the mass-spring N-body model can be generalised to inhomogeneous and/or non-spherical bodies.


title: Consistent quantization of massless fields of any spin and the generalized Maxwell's equations

abstract: A simplified formalism of first quantized massless fields of any spin is presented. The angular momentum basis for particles of zero mass and finite spin s of the D^(s-1/2,1/2) representation of the Lorentz group is used to describe the wavefunctions. The advantage of the formalism is that by equating to zero the s-1 components of the wave functions, the 2s-1 subsidiary conditions (needed to eliminate the non-forward and non-backward helicities) are automatically satisfied. Probability currents and Lagrangians are derived allowing a first quantized formalism. A simple procedure is derived for connecting the wave functions with potentials and gauge conditions. The spin 1 case is of particular interest and is described with the D^(1/2,1/2) vector representation of the well known self-dual representation of the Maxwell's equations. This representation allows us to generalize Maxwell's equations by adding the E_0 and B_0 components to the electric and magnetic four-vectors. Restrictions on their existence are discussed.


title: New classes of degree sequences with fast mixing swap Markov chain sampling

abstract: In network modeling of complex systems one is often required to sample random realizations of networks that obey a given set of constraints, usually in form of graph measures. A much studied class of problems targets uniform sampling of simple graphs with given degree sequence or also with given degree correlations expressed in the form of a joint degree matrix. One approach is to use Markov chains based on edge switches (swaps) that preserve the constraints, are irreducible (ergodic) and fast mixing. In 1999, Kannan, Tetali and Vempala (KTV) proposed a simple swap Markov chain for sampling graphs with given degree sequence and conjectured that it mixes rapidly (in poly-time) for arbitrary degree sequences. While the conjecture is still open, it was proven for special degree sequences, in particular, for those of undirected and directed regular simple graphs, of half-regular bipartite graphs, and of graphs with certain bounded maximum degrees. Here we prove the fast mixing KTV conjecture for novel, exponentially large classes of irregular degree sequences. Our method is based on a canonical decomposition of degree sequences into split graph degree sequences, a structural theorem for the space of graph realizations and on a factorization theorem for Markov chains. After introducing bipartite splitted degree sequences, we also generalize the canonical split graph decomposition for bipartite and directed graphs.


title: Twisted Interferometry: the topological perspective

abstract: Three manifold topology is used to analyze the effect of anyonic interferometers in which the probe anyons' path along an arm crosses itself, leading to a "twisted" or braided space-time trajectory for the probe anyons. In the case of Ising non-Abelian anyons, twisted interferometry is shown to be able to generate a topologically protected $\pi/8$-phase gate, which cannot be generated from quasiparticle braiding.


title: Differences in chemical doping matter - Superconductivity in Ti1-xTaxSe2 but not in Ti1-xNbxSe2

abstract: We report that 1T-TiSe2, an archetypical layered transition metal dichalcogenide, becomes superconducting when Ta is substituted for Ti but not when Nb is substituted for Ti. This is unexpected because Nb and Ta should be chemically equivalent electron donors. Superconductivity emerges near x = 0.02 for Ti1-xTaxSe2, while for Ti1-xNbxSe2, no superconducting transitions are observed above 0.4 K. The equivalent chemical nature of the dopants is confirmed by X-ray photoelectron spectroscopy. ARPES and Raman scattering studies show similarities and differences between the two systems, but the fundamental reasons why the Nb and Ta dopants yield such different behavior are unknown. We present a comparison of the electronic phase diagrams of many electron-doped 1T-TiSe2 systems, showing that they behave quite differently, which may have broad implications in the search for new superconductors. We propose that superconducting Ti0.8Ta0.2Se2 will be suitable for devices and other studies based on exfoliated crystal flakes.


title: Using Reed-Solomon codes in the $\left( U\mid U+V\right)$ construction and an application to cryptography

abstract: In this paper we present a modification of Reed-Solomon codes that beats the Guruwami-Sudan $1-\sqrt{R}$ decoding radius of Reed-Solomon codes at low rates $R$. The idea is to choose Reed-Solomon codes $U$ and $V$ with appropriate rates in a $\left( U\mid U+V\right)$ construction and to decode them with the Koetter-Vardy soft information decoder. We suggest to use a slightly more general version of these codes (but which has the same decoding performances as the $\left( U\mid U+V\right)$-construction) for code-based cryptography, namely to build a McEliece scheme. The point is here that these codes not only perform nearly as well (or even better in the low rate regime) as Reed-Solomon codes, their structure seems to avoid the Sidelnikov-Shestakov attack which broke a previous McEliece proposal based on generalized Reed-Solomon codes.


title: The Void Galaxy Survey: Star Formation Properties

abstract: We study the star formation properties of 59 void galaxies as part of the Void Galaxy Survey (VGS). Current star formation rates are derived from $\rm{H\alpha}$ and recent star formation rates from near-UV imaging. In addition, infrared 3.4 $\rm{\mu m}$, 4.6 $\rm{\mu m}$, 12 $\rm{\mu m}$ and 22 $\rm{\mu m}$ WISE emission is used as star formation and mass indicator. Infrared and optical colours show that the VGS sample displays a wide range of dust and metallicity properties. We combine these measurements with stellar and HI masses to measure the specific SFRs ($\rm{SFR/M_{*}}$) and star formation efficiencies ($\rm{SFR/M_{HI}}$). We compare the star formation properties of our sample with galaxies in the more moderate density regions of the cosmic web, 'the field'. We find that specific SFRs of the VGS galaxies as a function of stellar and HI mass are similar to those of the galaxies in these field regions. Their $\rm{SFR\alpha}$ is slightly elevated than the galaxies in the field for a given total HI mass. In the global star formation picture presented by Kennicutt-Schmidt, VGS galaxies fall into the regime of low average star formation and correspondingly low HI surface density. Their mean $\rm{SFR\alpha/M_{HI}}$ and $\rm{SFR\alpha/M_{*}}$ are of the order of $\rm{10^{-9.9}}$ $\rm{yr^{-1}}$. We conclude that while the large scale underdense environment must play some role in galaxy formation and growth through accretion, we find that even with respect to other galaxies in the more mildly underdense regions, the increase in star formation rate is only marginal.


title: On the geometry of border rank algorithms for matrix multiplication and other tensors with symmetry

abstract: We establish basic information about border rank algorithms for the matrix multiplication tensor and other tensors with symmetry. We prove that border rank algorithms for tensors with symmetry (such as matrix multiplication and the determinant polynomial) come in families that include representatives with normal forms. These normal forms will be useful both to develop new efficient algorithms and to prove lower complexity bounds. We derive a border rank version of the substitution method used in proving lower bounds for tensor rank. We use this border-substitution method and a normal form to improve the lower bound on the border rank of matrix multiplication by one, to 2n^2- n+1. We also point out difficulties that will be formidable obstacles to future progress on lower complexity bounds for tensors because of the "wild" structure of the Hilbert scheme of points.


title: Hairy black holes in the XX-th and XXI-st centuries

abstract: This is a brief summary of the most important hairy black hole solutions in 3+1 spacetime dimensions discovered over the last 25 years. These were first of all the Einstein-Yang-Mills black holes and their various generalizations including the Higgs field, the dilaton and the curvature corrections, and also the Skyrme black holes. More recently, these were black holes supporting a scalar field violating the energy conditions or non-minimally coupled to gravity, and also spinning black holes with massive complex scalar hair. Finally, these were black holes with massive graviton hair.


title: Generalized nil-Coxeter algebras over discrete complex reflection groups

abstract: We define and study generalized nil-Coxeter algebras associated to Coxeter groups. Motivated by a question of Coxeter (1957), we construct the first examples of such finite-dimensional algebras that are not the 'usual' nil-Coxeter algebras: a novel 2-parameter type $A$ family that we call $NC_A(n,d)$. We explore several combinatorial properties of $NC_A(n,d)$, including its Coxeter word basis, length function, and Hilbert-Poincare series, and show that the corresponding generalized Coxeter group is not a flat deformation of $NC_A(n,d)$. These algebras yield symmetric semigroup module categories that are necessarily not monoidal; we write down their Tannaka-Krein duality. Further motivated by the Broue-Malle-Rouquier (BMR) freeness conjecture [J. reine angew. math. 1998], we define generalized nil-Coxeter algebras over all discrete real or complex reflection groups $W$, finite or infinite. We provide a complete classification of all such algebras that are finite-dimensional. Remarkably, these turn out to be either the usual nil-Coxeter algebras, or the algebras $NC_A(n,d)$. This proves as a special case - and strengthens - the lack of equidimensional nil-Coxeter analogues for finite complex reflection groups. In particular, generic Hecke algebras are not flat deformations of $NC_W$ for $W$ complex.


title: Interaction of chiral rafts in self-assembled colloidal membranes

abstract: Colloidal membranes are monolayer assemblies of rodlike particles that capture the long-wavelength properties of lipid bilayer membranes on the colloidal scale. Recent experiments on colloidal membranes formed by chiral rodlike viruses showed that introducing a second species of virus with different length and opposite chirality leads to the formation of rafts --- micron-sized domains of one virus species floating in a background of the other viruses [Sharma et al., Nature 513, 77 (2014)]. In this article we study the interaction of such rafts using liquid crystal elasticity theory. By numerically minimizing the director elastic free energy, we predict the tilt angle profile for both a single raft and two rafts in a background membrane, and the interaction between two rafts as a function of their separation. We find that the chiral penetration depth in the background membrane sets the scale for the range of the interaction. We compare our results with the experimental data and find good agreement for the strength and range of the interaction. Unlike the experiments, however, we do not observe a complete collapse of the data when rescaled by the tilt angle at the raft edge.


title: Predicting Bulk Metallic Glass Forming Ability with the Thermodynamic Density of Competing Crystalline States

abstract: Despite two decades of studies, the formation of metallic glasses, very promising systems for industrial applications, still remains mostly unexplained. This lack of knowledge hinders the search for new systems, still performed with combinatorial trial and error. In the past, it was speculated that some sort of "confusion" during crystallization could play a key role during their formation. In this article, we propose a heuristic descriptor quantifying such confusion. It is based on the "thermodynamic density of competing crystalline states", parameterized from high-throughput ab-initio calculations. The existence of highly enthalpy-degenerate but geometrically-different phases frustrates the crystallization process and promotes glass formation. Two test beds are considered. A good and a bad glass-former, CuZr and NiZr, are experimentally characterized with high-throughput synthesis. The experimental results corroborate the capability of the heuristic descriptor in predicting glass forming ability through the compositional space. Our analysis is expected to deepen the understanding of the underlying mechanisms and to accelerate the discovery of novel metallic glasses.


title: Low-energy neutrinos

abstract: There exist several kinds of sources emitting neutrinos in the MeV energy range. These low-energy neutrinos from different sources can be often detected by the same multipurpose detectors. The status-of-art of the feld of solar neutrinos, geoneutrinos, and the search for sterile neutrino with artifcial neutrino sources is provided here; other neutrino sources, as for example reactor or high-energy neutrinos, are described elsewhere. For each of these three felds, the present-day motivation and open questions, as well as the latest experimental results and future perspectives are discussed.


title: Stable Gapless Bose Liquid Phases without any Symmetry

abstract: It is well-known that a stable algebraic spin liquid state (or equivalently an algebraic Bose liquid (ABL) state) with emergent gapless photon excitations can exist in quantum spin ice systems, or in a quantum dimer model on a bipartite $3d$ lattice. This photon phase is stable against any weak perturbation without assuming any symmetry. Further works concluded that certain lattice models give rise to more exotic stable algebraic Bose liquid phases with graviton-like excitations. In this paper we will show how these algebraic Bose liquid states can be generalized to stable phases with even more exotic types of gapless excitations and then argue that these new phases are stable against weak perturbations. We also explicitly show that these theories have an (algebraic) topological ground state degeneracy on a torus, and construct the corresponding topological invariants.


title: THz conductivity of Sr$_{1-x}$Ca$_x$RuO$_3$

abstract: We investigate the optical conductivity of Sr$_{1-x}$Ca$_x$RuO$_3$ across the ferromagnetic to paramagnetic transition that occurs at $x=0.8$. The thin films were grown by metalorganic aerosol deposition with $0 \leq x \leq 1$ onto NdGaO$_3$ substrates. We performed THz frequency domain spectroscopy in a frequency range from 3~cm$^{-1}$ to 40~cm$^{-1}$ (100~GHz to 1.4~THz) and at temperatures ranging from 5~K to 300~K, measuring transmittivity and phase shift through the films. From this we obtained real and imaginary parts of the optical conductivity. The end-members, ferromagnetic SrRuO$_3$ and paramagnetic CaRuO$_3$, show a strongly frequency-dependent metallic response at temperatures below 20~K. Due to the high quality of these samples we can access pronounced intrinsic electronic contributions to the optical scattering rate, which at 1.4~THz exceeds the residual scattering rate by more than a factor of three. Deviations from a Drude response start at about 0.7~THz for both end-members in a remarkably similar way. For the intermediate members a higher residual scattering originating in the compositional disorder leads to a featureless optical response, instead. The relevance of low-lying interband transitions is addressed by a calculation of the optical conductivity within density functional theory in the local density approximation (LDA).


title: The omega-inequality problem for concatenation hierarchies of star-free languages

abstract: The problem considered in this paper is whether an inequality of omega-terms is valid in a given level of a concatenation hierarchy of star-free languages. The main result shows that this problem is decidable for all (integer and half) levels of the Straubing-Th\'erien hierarchy.


title: On Elliptic Algebras and Large-n Supersymmetric Gauge Theories

abstract: In this note we further develop the duality between supersymmetric gauge theories in various dimensions and elliptic integrable systems such as Ruijsenaars-Schneider model and periodic intermediate long wave hydrodynamics. These models arise in instanton counting problems and are described by certain elliptic algebras. We discuss the correspondence between the two types of models by employing the large-n limit of the dual gauge theory. In particular we provide non-Abelian generalization of our previous result on the intermediate long wave model.


title: A Statistical Framework for the Utilization of Simultaneous Pupil Plane and Focal Plane Telemetry for Exoplanet Imaging, Part I: Accounting for Aberrations in Multiple Planes

abstract: A new generation of telescopes with mirror diameters of 20 m or more, called extremely large telescopes (ELTs) has the potential to provide unprecedented imaging and spectroscopy of exo-planetary systems, if the difficulties in achieving the extremely high dynamic range required to differentiate the planetary signal from the star can be overcome to a sufficient degree. Fully utilizing the potential of ELTs for exoplanet imaging will likely require simultaneous and self-consistent determination of both the planetary image and the unknown aberrations in multiple planes of the optical system, using statistical inference based on the wavefront sensor and science camera data streams. This approach promises to overcome the most important systematic errors inherent in the various schemes based on differential imaging, such as ADI and SDI. This paper is the first in a series on this subject, in which a formalism is established for the exoplanet imaging problem, setting the stage for the statistical inference methods to follow in the future. Every effort has been made to be rigorous and complete, so that validity of approximations to be made later can be assessed. Here, the polarimetric image is expressed in terms of aberrations in the various planes of a polarizing telescope with an adaptive optics system. Further, it is shown that current methods that utilize focal plane sensing to correct the speckle field, e.g., electric field conjugation, rely on the tacit assumption that aberrations on multiple optical surfaces can be represented as aberration on a single optical surface, ultimately limiting their potential effectiveness for ground-based astronomy.


title: Doubling Constructions for Covering Groups and Tensor Product L-Functions

abstract: This is a research announcement concerning a series of constructions obtained by applying the "doubling method" from the theory of automorphic forms to covering groups. Using these constructions, we obtain partial tensor product L-functions attached to generalized Shimura lifts, which may be defined in a natural way since at almost all places the representations are unramified principal series.


title: Homotopical Complexity of a 3D Billiard Flow

abstract: In this paper we study the homotopical rotation vectors and the homotopical rotation sets for the billiard flow on the unit flat torus with three, mutually intersecting and mutually orthogonal cylindrical scatterers removed from it. The natural habitat for these objects is the infinite cone erected upon the Cantor set $\text{Ends}(\textbf{F}_3)$ of all "ends" of the hyperbolic group $\textbf{F}_3=\pi_1(\mathbf{Q})$. An element of $\text{Ends}(\textbf{F}_3)$ describes the direction in (the Cayley graph of) the group $\textbf{F}_3$ in which the considered trajectory escapes to infinity, whereas the height function $s$ ($s\ge 0$) of the cone gives us the average speed at which this escape takes place. The main results of this paper claim that the orbits can only escape to infinity at a speed not exceeding $\sqrt{3}$, and in any direction $e\in\text{Ends}(\mathbf{F}_3)$ the escape is feasible with any prescribed speed $s$, $0\leq s\leq 1/3$. This means that the radial upper and lower bounds for the rotation set $R$ are actually pretty close to each other. Furthermore, we prove the convexity of the set $AR$ of constructible rotation vectors, and that the set of rotation vectors of periodic orbits is dense in $AR$. We also provide effective lower and upper bounds for the topological entropy of the studied billiard flow.


title: Implicit Filtered PN for High-Energy Density Thermal Radiation Transport using Discontinuous Galerkin Finite Elements

abstract: In this work, we provide a fully-implicit implementation of the time-dependent, filtered spherical harmonics (FPN) equations for non-linear, thermal radiative transfer. We investigate local filtering strategies and analyze the effect of the filter on the conditioning of the system, showing in particular that the filter improves the convergence properties of the iterative solver. We also investigate numerically the rigorous error estimates derived in the linear setting, to determine whether they hold also for the non-linear case. Finally, we simulate a standard test problem on an unstructured mesh and make comparisons with implicit Monte-Carlo (IMC) calculations.


title: Thermal Transport in Graphene, Few-Layer Graphene and Graphene Nanoribbons

abstract: The discovery of unusual heat conduction properties of graphene has led to a surge of theoretical and experimental studies of phonon transport in two-dimensional material systems. The rapidly developing graphene thermal field spans from theoretical physics to practical engineering applications. In this invited review we outline different theoretical approaches developed for describing phonon transport in graphene and provide comparison with available experimental thermal conductivity data. A special attention is given to analysis of the recent theoretical results for the phonon thermal conductivity of graphene and few-layer graphene, the effects of the strain, defects, isotopes and edge scattering on the acoustic phonon transport in these material systems.


title: Categorical Judgment with a Variable Decision Rule

abstract: A new Thurstonian rating scale model uses a variable decision rule (VDR) that incorporates three previously formulated, distinct decision rules. The model includes probabilities for choosing each rule, along with Gaussian representation and criterion densities. Numerical optimisation techniques were validated through demonstrating that the model fits simulated data tightly. For simulations with 400 trials per stimulus (tps), useful information emerged about the generating parameters. However, larger experiments (e.g. 4000 tps) proved desirable for better recovery of generating parameters and to support trustworthy choices between competing models by the Akaike Information Criterion. In reanalyses of experiments by others, the VDR model explained most of the data better than did classical signal detection theory models.


title: Kalman-Filter-Based Particle Tracking on Parallel Architectures at Hadron Colliders

abstract: Power density constraints are limiting the performance improvements of modern CPUs. To address this we have seen the introduction of lower-power, multi-core processors such as GPGPU, ARM and Intel MIC. To stay within the power density limits but still obtain Moore's Law performance/price gains, it will be necessary to parallelize algorithms to exploit larger numbers of lightweight cores and specialized functions like large vector units. Track finding and fitting is one of the most computationally challenging problems for event reconstruction in particle physics. At the High-Luminosity Large Hadron Collider (HL-LHC), for example, this will be by far the dominant problem. The need for greater parallelism has driven investigations of very different track finding techniques such as Cellular Automata or Hough Transforms. The most common track finding techniques in use today, however, are those based on the Kalman Filter. Significant experience has been accumulated with these techniques on real tracking detector systems, both in the trigger and offline. They are known to provide high physics performance, are robust, and are in use today at the LHC. We report on porting these algorithms to new parallel architectures. Our previous investigations showed that, using optimized data structures, track fitting with Kalman Filter can achieve large speedups both with Intel Xeon and Xeon Phi. We report here our progress towards an end-to-end track reconstruction algorithm fully exploiting vectorization and parallelization techniques in a realistic experimental environment.


title: Millimeter-Wave Broadband Anti-Reflection Coatings Using Laser Ablation of Sub-Wavelength Structures

abstract: We report on the first use of laser ablation to make sub-millimeter, broad-band, anti-reflection coatings (ARC) based on sub-wavelength structures (SWS) on alumina and sapphire. We used a 515 nm laser to produce pyramid-shaped structures with pitch of about 320 $\mu$m and total height of near 800 $\mu$m. Transmission measurements between 70 and 140 GHz are in agreement with simulations using electromagnetic propagation software. The simulations indicate that SWS ARC with the fabricated shape should have a fractional bandwidth response of $\Delta \nu / \nu_{center}\ = 0.55$ centered on 235 GHz for which reflections are below 3%. Extension of the bandwidth to both lower and higher frequencies, between few tens of GHz and few THz, should be straightforward with appropriate adjustment of laser ablation parameters.


title: A nonlinear model for rotationally constrained convection with Ekman pumping

abstract: It is a well established result of linear theory that the influence of differing mechanical boundary conditions, i.e., stress-free or no-slip, on the primary instability in rotating convection becomes asymptotically small in the limit of rapid rotation. This is accounted for by the diminishing impact of the viscous stresses exerted within Ekman boundary layers and the associated vertical momentum transport by Ekman pumping. By contrast, in the nonlinear regime recent experiments and supporting simulations are now providing evidence that the efficiency of heat transport remains strongly influenced by Ekman pumping in the rapidly rotating limit. In this paper, a reduced model is developed for the case of low Rossby number convection in a plane layer geometry with no-slip upper and lower boundaries held at fixed temperatures. A complete description of the dynamics requires the existence of three distinct regions within the fluid layer: a geostrophically balanced interior where fluid motions are predominately aligned with the axis of rotation, Ekman boundary layers immediately adjacent to the bounding plates, and thermal wind layers driven by Ekman pumping in between. The reduced model uses a classical Ekman pumping parameterization to alleviate the need for spatially resolving the Ekman boundary layers. Results are presented for both linear stability theory and a special class of nonlinear solutions described by a single horizontal spatial wavenumber. It is shown that Ekman pumping allows for significant enhancement in the heat transport relative to that observed in simulations with stress-free boundaries. Without the intermediate thermal wind layer the nonlinear feedback from Ekman pumping would be able to generate a heat transport that diverges to infinity. This layer arrests this blowup resulting in finite heat transport at a significantly enhanced value.


title: A fully discrete BEM-FEM scheme for transient acoustic waves

abstract: We study a symmetric BEM-FEM coupling scheme for the scattering of transient acoustic waves by bounded inhomogeneous anisotropic obstacles in a homogeneous field. An incident wave in free space interacts with the obstacles and produces a combination of transmission and scattering. The transmitted part of the wave is discretized in space by finite elements while the scattered wave is reduced to two fields defined on the boundary of the obstacles and is discretized in space with boundary elements. We choose a coupling formulation that leads to a symmetric system of integro-differential equations. The retarded boundary integral equations are discretized in time by Convolution Quadrature, and the interior field is discretized in time with the trapezoidal rule. We show that the scattering problem generates a C_0 group of isometries in a Hilbert space, and use associated estimates to derive stability and convergence results. We provide numerical experiments and simulations to validate our results and demonstrate the flexibility of the method.


title: Plasmon excitations in layered high-Tc cuprates

abstract: Motivated by the recent resonant inelastic x-ray scattering (RIXS) experiment for the electron-doped cuprates Nd_{2-x}Ce_{x}CuO_{4} with x around 0.15, we compute the density-density correlation function in the t-J model on a square lattice by including interlayer hopping and the long-range Coulomb interaction. We find that collective charge excitations are realized not inside the particle-hole continuum, but above the continuum as plasmons. The plasmon mode has a rather flat dispersion near the in-plane momentum q=(0,0) with a typical excitation energy of the order of the intralayer hopping t when the out-of-plane momentum qz is zero. However, when qz becomes finite, the plasmon dispersion changes drastically near q=(0,0), leading to a strong dispersive feature with an excitation gap scaled by the interlayer hopping tz. We discuss that the mode recently observed by RIXS near q=(0,0) can be associated with the plasmon mode with a finite qz.


title: Anomalous increase in nematic-isotropic transition temperature in dimer molecules induced by magnetic field

abstract: We have determined the nematic-isotropic transition temperature as a function of applied magnetic field in three different thermotropic liquid crystalline dimers. These molecules are comprised of two rigid calamitic moieties joined end to end by flexible spacers with odd numbers of methylene groups. They show an unprecedented magnetic field enhancement of nematic order in that the transition temperature is increased by up to 15K when subjected to 22T magnetic field. The increase is conjectured to be caused by a magnetic field-induced decrease of the average bend angle in the aliphatic spacers connecting the rigid mesogenic units of the dimers.


title: Rectifiability, interior approximation and Harmonic Measure

abstract: We prove a structure theorem for any $n$-rectifiable set $E\subset \mathbb{R}^{n+1}$, $n\ge 1$, satisfying a weak version of the lower ADR condition, and having locally finite $H^n$ ($n$-dimensional Hausdorff) measure. Namely, that $H^n$-almost all of $E$ can be covered by a countable union of boundaries of bounded Lipschitz domains contained in $\mathbb{R}^{n+1}\setminus E$. As a consequence, for harmonic measure in the complement of such a set $E$, we establish a non-degeneracy condition which amounts to saying that $H^n|_E$ is "absolutely continuous" with respect to harmonic measure in the sense that any Borel subset of $E$ with strictly positive $H^n$ measure has strictly positive harmonic measure in some connected component of $\mathbb{R}^{n+1}\setminus E$. We also provide some counterexamples showing that our result for harmonic measure is optimal. Moreover, we show that if, in addition, a set $E$ as above is the boundary of a connected domain $\Omega \subset \mathbb{R}^{n+1}$ which satisfies an infinitesimal interior thickness condition, then $H^n|_{\partial\Omega}$ is absolutely continuous (in the usual sense) with respect to harmonic measure for $\Omega$. Local versions of these results are also proved: if just some piece of the boundary is $n$-rectifiable then we get the corresponding absolute continuity on that piece. As a consequence of this and recent results by Azzam-Hofmann-Martell-Mayboroda-Mourgoglou-Tolsa-Volberg, we can decompose the boundary of any open connected set satisfying the previous conditions in two disjoint pieces: one that is $n$-rectifiable where Hausdorff measure is absolutely continuous with respect to harmonic measure and another purely $n$-unrectifiable piece having vanishing harmonic measure.


title: A New Numerical Method for Fast Solution of Partial Integro-Differential Equations

abstract: A new method of numerical solution for partial differential equations is proposed. The method is based on a fast matrix multiplication algorithm. Two-dimensional Poison equation is used for comparison of the proposed method with conventional numerical methods. It was shown that the new method allows for linear growth in the number of elementary addition and multiplication operations with the growth of grid size, as contrasted with quadratic growth necessitated by the standard numerical methods. The proposed method can be easily generalized for any differential equations.


title: The abundance and spatial distribution of ultra-diffuse galaxies in nearby galaxy clusters

abstract: Recent observations have highlighted a significant population of faint but large (r_eff>1.5 kpc) galaxies in the Coma cluster. The origin of these Ultra Diffuse Galaxies (UDGs) remains puzzling, as the interpretation of the observational results has been hindered by the subjective selection of UDGs, and the limited study of only the Coma (and some examples in the Virgo-) cluster. We extend the study of UDGs using 8 clusters in the redshift range 0.044<z<0.063 with deep g- and r-band imaging data taken with MegaCam at the CFHT. We describe an automatic selection pipeline for quantitative identification, tested for completeness using image simulations of these galaxies. We find that the abundance of the UDGs we can detect increases with cluster mass, reaching ~200 in typical haloes of M200~10^15 Msun. The cluster UDGs have colours consistent with the cluster red sequence, and have a steep size distribution that declines as n~r_eff^-3.4. Their radial distribution is significantly steeper than NFW in the outskirts, and is significantly shallower in the inner parts. They follow the same radial distribution as the more massive quiescent galaxies in the clusters, except within the core region of r<0.15XR200 (or <300 kpc). Within this region the number density of UDGs drops and is consistent with zero. These diffuse galaxies can only resist tidal forces down to this cluster-centric distance if they are highly centrally dark-matter dominated. The observation that the radial distribution of more compact dwarf galaxies (r_eff<1.0 kpc) with similar luminosities follows the same distribution as the UDGs, but exist down to a smaller distance of 100kpc from the cluster centres, indicates that they may have similarly massive sub-haloes as the UDGs. Although several scenarios can give rise to the UDG population, our results point to differences in the formation history as the most plausible explanation.


title: Resolution-independent modeling of environmental effects in semi-analytic models of galaxy formation that include ram-pressure stripping of both hot and cold gas

abstract: The quenching of star formation in satellite galaxies is observed over a wide range of dark matter halo masses and galaxy environments. In the recent Guo et al (2011) and Fu et al (2013) semi-analytic + N-body models, the gaseous environment of the satellite galaxy is governed by the properties of the dark matter subhalo in which it resides. This quantity depends of the resolution of the N-body simulation, leading to a divergent fraction of quenched satellites in high- and low-resolution simulations. Here, we incorporate an analytic model to trace the subhaloes below the resolution limit. We demonstrate that we then obtain better converged results between the Millennium I and II simulations, especially for the satellites in the massive haloes ($\rm log M_{halo}=[14,15]$). We also include a new physical model for the ram-pressure stripping of cold gas in satellite galaxies. However, we find very clear discrepancies with observed trends in quenched satellite galaxy fractions as a function of stellar mass at fixed halo mass. At fixed halo mass, the quenched fraction of satellites does not depend on stellar mass in the models, but increases strongly with mass in the data. In addition to the over-prediction of low-mass passive satellites, the models also predict too few quenched central galaxies with low stellar masses, so the problems in reproducing quenched fractions are not purely of environmental origin. Further improvements to the treatment of the gas-physical processes regulating the star formation histories of galaxies are clearly necessary to resolve these problems.


title: Di-photon excess at LHC and the gamma ray excess at the Galactic Centre

abstract: Motivated by the recent indications for a 750 GeV resonance in the di-photon final state at the LHC, in this work we analyse the compatibility of the excess with the broad photon excess detected at the Galactic Centre. Intriguingly, by analysing the parameter space of an effective models where a 750 GeV pseudoscalar particles mediates the interaction between the Standard Model and a scalar dark sector, we prove the compatibility of the two signals. We show, however, that the LHC mono-jet searches and the Fermi LAT measurements strongly limit the viable parameter space. We comment on the possible impact of cosmic antiproton flux measurement by the AMS-02 experiment.