title: Global instability in the Ghil-Sellers model

abstract: The Ghil-Sellers energy balance model of Earth's climate, features -- for a considerable range of the solar intensity -- two stable climate states (a warm and a cold snowball Earth), where the bistability results from the celebrated ice-albedo feedback. The unstable solution is obtained and characterized in this paper. We find such unstable states by applying for the first time in a geophysical context the so-called edge tracking method that has been used for studying multiple coexisting states in shear flows. We examine robustness, efficiency, and accuracy properties of the edge tracking algorithm. We find that the procedure is the most efficient when taking a single bisection per cycle. Due to the strong diffusivity of the system trajectories of transient dynamics, initialized between the stable states with respect to the mean temperature, are confined to the heteroclininc trajectory, one which connects the fixed unstable and stable states, after relatively short transient times. This constraint dictates a functional relationship between observables. We characterize such a relationship between the global average temperature and a descriptor of nonequilibrium thermodynamics, the large scale temperature gradient between low and high latitudes. We find that a maximum of the temperature gradient is realized at the same value of the average temperature, about 270 K, largely independent of the strength of incoming solar radiation. We also examine the structural properties of the system defined by the equilibria depending on the solar strength. We construct new bifurcation diagrams in terms of quantities relevant for describing the thermodynamic disequilibrium, such as the temperature gradient and the material entropy production due to heat transport. We compare our results for the EBM to results for the intermediate complexity GCM PlaSim and find an interesting qualitative agreement.


title: On monodromy representations in Denham-Suciu fibrations

abstract: We study the monodromy representation corresponding to a fibration introduced by G. Denham and A. Suciu, which involves polyhedral products given in Definition 2.2. Algebraic and geometric descriptions for these monodromy representations are given. In particular, we study the case of a product of two finite cyclic groups and obtain representations into $Out(F_n)$ and $SL_n(\mathbb Z)$. We give algebraic descriptions of monodromy for the case of a product of any two finite groups . Finally we give a geometric description for monodromy representations of a product of 2 or more finite groups to $Out(F_n)$, as well as some algebraic properties. The geometric description does not rely on choosing a basis for the fundamental group of the fibre in terms of commutators, hence avoids this delicate question.


title: Excluding the Light Dark Matter Window of a 331 Model Using LHC and Direct Dark Matter Detection Data

abstract: We sift the impact of the recent Higgs precise measurements, and recent dark matter direct detection results, on the dark sector of an electroweak extension of the Standard Model that has a complex scalar as dark matter. We find that in this model the Higgs decays with a large branching ratio into dark matter particles, and charged scalars when these are kinematically available, for any coupling strength differently from the so called Higgs portal. Moreover, we compute the abundance and spin-independent WIMP-nucleon scattering cross section, which are driven by the Higgs and $Z^{\prime}$ boson processes. We decisively exclude the $1-500$~GeV dark matter window and find the most stringent lower bound in the literature on the scale of symmetry breaking of the model namely $10$~TeV, after applying the LUX-2013 limit. Interestingly, the projected XENON1T constraint will be able to rule out the entire $1$~GeV-$1000$~GeV dark matter mass range. Lastly, for completeness, we compute the charged scalar production cross section at the LHC and comment on the possibility of detection at current and future LHC runnings.


title: Effective-range approximations for resonant scattering of cold atoms

abstract: Studies of cold atom collisions and few-body interactions often require the energy dependence of the scattering phase shift, which is usually expressed in terms of an effective-range expansion. We use accurate coupled-channel calculations on $^{6}$Li, $^{39}$K and $^{133}$Cs to explore the behavior of the effective range in the vicinity of both broad and narrow Feshbach resonances. We show that commonly used expressions for the effective range break down dramatically for narrow resonances and near the zero-crossings of broad resonances. We present an alternative parametrization of the effective range that is accurate through both the pole and the zero-crossing for both broad and narrow resonances. However, the effective range expansion can still fail at quite low collision energies, particularly around narrow resonances. We demonstrate that an analytical form of an energy and magnetic field-dependent phase shift, based on multichannel quantum defect theory, gives accurate results for the energy-dependent scattering length.


title: On the Construction of Doubly Even Order Magic Squares

abstract: Magic squares have been an enthralling topic in mathematics for centuries. They are formed by filling in all the cells of a square matrix with the numbers starting from one so that the sum of all rows, columns, and diagonals is the same. Magic squares have applications in entertainment, music and even cryptography. This paper focuses on a simple and easy method to construct doubly even magic squares.


title: On the Spatial Ambiguity Function

abstract: The ambiguity function of a spatial signal in the "signal spatial frequency displacement vs antenna linear displacement" coordinates is considered in this paper; in case of a monochromatic signal with known carrier frequency, spatial signal vector angular displacement is considered as the second coordinate. The analysis is made taking into account an external (reradiated or thermal) background which contributes considerable peculiarities in the properties of the spatial ambiguity function. Accuracy of measuring the signal linear or angular displacement and resolution are determined. Examples for the Fresnel and Fraunhofer zones are given.


title: Automorphism groups of simplicial complexes of infinite type surfaces

abstract: Let S be any orientable surface of infinite genus with a finite number of boundary components. In this work we consider the curve complex C(S), the nonseparating curve complex N(S) and the Schmutz graph G(S) of S. When all the topological ends of S carry genus, we show that all elements in the automorphism groups Aut(C(S)), Aut(N(S)) and Aut(G(S)) are geometric, i.e. these groups are naturally isomorphic to the extended mapping class group MCG*(S) of the infinite surface S. Finally, we study rigidity phenomena within Aut(C(S)) and Aut(N(S))


title: Thermodynamics with continuous information flow

abstract: We provide a unified thermodynamic formalism describing information transfers in autonomous as well as nonautonomous systems described by stochastic thermodynamics. We demonstrate how information is continuously generated in an auxiliary system and then transferred to a relevant system that can utilize it to fuel otherwise impossible processes. Indeed, while the joint system satisfies the second law, the entropy balance for the relevant system is modified by an information term related to the mutual information rate between the two systems. We show that many important results previously derived for nonautonomous Maxwell demons can be recovered from our formalism and use a cycle decomposition to analyze the continuous information flow in autonomous systems operating at steady-state. A model system is used to illustrate our findings.


title: Separating Regular Languages with First-Order Logic

abstract: Given two languages, a separator is a third language that contains the first one and is disjoint from the second one. We investigate the following decision problem: given two regular input languages of finite words, decide whether there exists a first-order definable separator. We prove that in order to answer this question, sufficient information can be extracted from semigroups recognizing the input languages, using a fixpoint computation. This yields an EXPTIME algorithm for checking first-order separability. Moreover, the correctness proof of this algorithm yields a stronger result, namely a description of a possible separator. Finally, we generalize this technique to answer the same question for regular languages of infinite words.


title: An enlargement of filtration formula with application to progressive enlargement with multiple random times

abstract: Given a reference filtration $\mathbb{F}$, we develop in this work a generic method for computing the semimartingale decomposition of $\mathbb{F}$-martingales in some specific enlargements of $\mathbb{F}$. This method is then applied to the study of progressive enlargement with multiple non-ordered random times, for which explicit decompositions can be obtained under the absolute continuity condition of Jacod.


title: Variations on Classical and Quantum Extractors

abstract: Many constructions of randomness extractors are known to work in the presence of quantum side information, but there also exist extractors which do not [Gavinsky {\it et al.}, STOC'07]. Here we find that spectral extractors $\psi$ with a bound on the second largest eigenvalue $\lambda_{2}(\psi^{\dagger}\circ\psi)$ are quantum-proof. We then discuss fully quantum extractors and call constructions that also work in the presence of quantum correlations decoupling. As in the classical case we show that spectral extractors are decoupling. The drawback of classical and quantum spectral extractors is that they always have a long seed, whereas there exist classical extractors with exponentially smaller seed size. For the quantum case, we show that there exists an extractor with extremely short seed size $d=O(\log(1/\epsilon))$, where $\epsilon>0$ denotes the quality of the randomness. In contrast to the classical case this is independent of the input size and min-entropy and matches the simple lower bound $d\geq\log(1/\epsilon)$.


title: $K$-theory and homotopies of 2-cocycles on transformation groups

abstract: This paper constitutes a first step in the author's program to investigate the question of when a homotopy of 2-cocycles $\omega = \{\omega_t\}_{t \in [0,1]}$ on a locally compact Hausdorff groupoid $\mathcal{G}$ induces an isomorphism of the $K$-theory groups of the reduced twisted groupoid $C^*$-algebras: $K_*(C^*_r(\mathcal{G}, \omega_0)) \cong K_*(C^*_r(\mathcal{G}, \omega_1)).$ Generalizing work of Echterhoff, L\"uck, Phillips, and Walters from 2010, we show that if $\mathcal{G} = G \ltimes X$ is a second countable locally compact transformation group, then whenever $G$ satisfies the Baum-Connes conjecture with coefficients, a homotopy $\omega = \{\omega_t\}_{t \in [0,1]}$ of 2-cocycles on $G \ltimes X$ gives rise to an isomorphism $K_*(C^*_r(G \ltimes X, \omega_0)) \cong K_*(C^*_r(G \ltimes X, \omega_1)).$


title: Partitioning Complex Networks via Size-constrained Clustering

abstract: The most commonly used method to tackle the graph partitioning problem in practice is the multilevel approach. During a coarsening phase, a multilevel graph partitioning algorithm reduces the graph size by iteratively contracting nodes and edges until the graph is small enough to be partitioned by some other algorithm. A partition of the input graph is then constructed by successively transferring the solution to the next finer graph and applying a local search algorithm to improve the current solution. In this paper, we describe a novel approach to partition graphs effectively especially if the networks have a highly irregular structure. More precisely, our algorithm provides graph coarsening by iteratively contracting size-constrained clusterings that are computed using a label propagation algorithm. The same algorithm that provides the size-constrained clusterings can also be used during uncoarsening as a fast and simple local search algorithm. Depending on the algorithm's configuration, we are able to compute partitions of very high quality outperforming all competitors, or partitions that are comparable to the best competitor in terms of quality, hMetis, while being nearly an order of magnitude faster on average. The fastest configuration partitions the largest graph available to us with 3.3 billion edges using a single machine in about ten minutes while cutting less than half of the edges than the fastest competitor, kMetis.


title: High-Energy Continuum of Magnetic Excitations in the Two-Dimensional Quantum Antiferromagnet Sr2CuO2Cl2

abstract: We have measured the magnetic excitation spectrum of the model square-lattice spin-1/2 antiferromagnet $\rm Sr_2CuO_2Cl_2$ over a broad range of energy and momentum using high-resolution inelastic neutron scattering (INS). The magnon dispersion along the zone boundary was accurately measured to be a 43 meV between (1/2,0) and (3/4,1/4) indicating the importance of coupling beyond nearest neighbors in the spin Hamiltonian. We observe a strong momentum dependent damping of the zone-boundary magnons at (1/2,0) revealing a high energy continuum of magnetic excitations. A direct comparison between our INS measurements and resonant inelastic X-ray scattering (RIXS) measurements shows that the RIXS spectrum contains significant contributions from higher energy excitations not previously considered. Our observations demonstrate that this high-energy continuum of magnetic fluctuations is a ubiquitous feature of insulating monolayer cuprates, apparent in both inelastic neutron and light scattering measurements.


title: Threshold state and a conjecture of Poghosyan, Poghosyan, Priezzhev and Ruelle

abstract: We prove a precise relationship between the threshold state of the fixed-energy sandpile and the stationary state of Dhar's abelian sandpile: In the limit as the initial condition tends to negative infinity, the former is obtained by size-biasing the latter according to burst size, an avalanche statistic. The question of whether and how these two states are related has been a subject of some controversy since 2000. The size-biasing in our result arises as an instance of a Markov renewal theorem, and implies that the threshold and stationary distributions are not equal even in the negative infinity limit. We prove that nevertheless in this limit the total amount of sand in the threshold state converges in distribution to the total amount of sand in the stationary state, confirming a conjecture of Poghosyan, Poghosyan, Priezzhev and Ruelle.


title: Cavatappi 2.0: More of the same but better

abstract: Innocent musing on geodesics on the surface of helical pasta shapes leads to a single continuous 4-parameter family of surfaces invariant under at least a 1-parameter symmetry group and which contains as various limits spheres, tori, helical tubes, and cylinders, all useful for illustrating various aspects of geometry in a visualizable setting that are important in special and general relativity. In this family the most aesthetically pleasing surfaces come from screw-rotating a plane cross-sectional curve perpendicular to itself, i.e., orthogonal to the tangent vector to a helix. If we impose instead this orthogonality in the Lorentzian geometry of 3-dimensional Minkowski spacetime with a timelike helical "central" world line representing a circular orbit, we can model the Fermi Born rigid model of the classical electron in such an orbit around the nucleus, and visualize the Fermi coordinate grid and its intersection with the world tube of the equatorial circle of the spherical surface of the electron (suppressing one spatial dimension). This leads what we might playfully term "relativistic pasta". This is a useful 2-dimensional stationary spacetime with closed spatial slices on which to illustrate the slicing and threading splittings of general relativity relative to a Killing congruence, like stationary axisymmetric spacetimes including the rotating black hole family.


title: Can the clicks of the detectors provide a complete description of Nature?

abstract: No matter how counterintuitive they are, quantum phenomena are all simple consequences of the laws of Quantum Mechanics. It is not needed to extend the theory with hidden mechanisms or additional principles to explain what Quantum Mechanics already predicts. This indubitable fact is often taken as supporting the view that all we can know about the universe comes from the outcomes of the quantum observations. According to this view, we can even learn the physical laws, in particular the properties of the space, particles, fields, and interactions, solely from the outcomes of the quantum observations. In this article it is shown that the unitary symmetry of the laws of Quantum Mechanics imposes severe restrictions in learning the physical laws of the universe, if we know only the observables and their outcomes.


title: More on Universal Superconformal Attractors

abstract: We define a general class of superconformal inflationary attractor models studied in our previous papers [1,2], which include, among others, inflationary models of 'induced gravity' which were argued to retain perturbative unitarity up to the Planck scale [3]. The choice of the function $f(\phi)$ in superconformal attractors for the models introduced in [3] is $f(\phi)= \phi^n - \xi^{-1}$. We present both superconformal and supergravity versions of these models, which were derived in [1,2] for arbitrary $f(\phi)$, together with the universal conditions required for stabilization of the extra three moduli present in the superconformal attractors, in addition to the inflaton.


title: Time flat surfaces and the monotonicity of the spacetime Hawking mass II

abstract: In this sequel paper we give a shorter, second proof of the monotonicity of the Hawking mass for time flat surfaces under spacelike uniformly area expanding flows in spacetimes that satisfy the dominant energy condition. We also include a third proof which builds on a known formula and describe a class of sufficient conditions of divergence type for the monotonicity of the Hawking mass. These flows of surfaces may have connections to the problem in general relativity of bounding the total mass of a spacetime from below by the quasi-local mass of spacelike 2-surfaces in the spacetime.


title: Two Steps to Obfuscation

abstract: This note addresses the historical antecedents of the 1998 PageRank measure of centrality. An identity relation links it to 1990-1991 models of Friedkin and Johnsen.


title: Quantum interference between independent reservoirs in open quantum systems

abstract: When a quantum system interacts with multiple reservoirs, the environmental effects are usually treated in an additive manner. We show that this assumption breaks down for non-Markovian environments that have finite memory times. Specifically, we demonstrate that quantum interferences between independent environments can qualitatively modify the dynamics of the physical system. We illustrate this effect with a two level system coupled to two structured photonic reservoirs, discuss its origin using a non-equilibrium diagrammatic technique, and show an example when the application of this interference can result in an improved dark state preparation in a $\lambda$ system.


title: On a digital quantum description of physical space versus the continuum description

abstract: The continuum of real numbers has served well as a model for physical space in mechanics and field theories. However it is a well-motivated and popular idea that at the fundamental Planck scale the combination of gravitational and quantum effects forces us to re-evaluate the concept of space (and time), and some sort of discreteness or granularity is expected. Here we study a cubic grid of points representing fundamental volumes of space and derive the spectra of momenta and wave number, which are not generally proportional to each other. The momentum is bounded and discrete for a finite grid. The uncertainty principle must be modified to reflect small scale spatial non-locality but it is consistent with the standard uncertainty principle in the continuum limit. We only consider kinematics in this work, and do not discuss time or energy or dynamical evolution.


title: Statistical quasi Cauchyness in two normed spaces

abstract: A function $f$ defined on a subset $E$ of a two normed space $X$ is statistically ward continuous if it preserves statistically quasi-Cauchy sequences of points in $E$ where a sequence $(x_n)$ is statistically quasi-Cauchy if $(\Delta x_{n})$ is a statistically null sequence. A subset $E$ of $X$ is statistically ward compact if any sequence of points in $E$ has a statistically quasi-Cauchy subsequence. In this paper, new kinds of continuities are investigated in two normed spaces. It turns out that uniform limit of statistically ward continuous functions is again statistically ward continuous.


title: Relativistic transfer ionization and the Breit interaction

abstract: We consider correlated transfer ionization in relativistic collisions between a highly charged ion and a light atom. In this process two quasi-free electrons of the atom interact with each other during the short collision time that results in capture of one of them by the ion and emission of the other. We show that this process is strongly influenced by the generalized Breit interaction already at modest relativistic impact energies.


title: Testing the equal-time angular-averaged consistency relation of the gravitational dynamics in N-body simulations

abstract: We explicitly test the equal-time consistency relation between the angular-averaged bispectrum and the power spectrum of the matter density field, employing a large suite of cosmological $N$-body simulations. This is the lowest-order version of the relations between $(\ell+n)$-point and $n$-point polyspectra, where one averages over the angles of $\ell$ soft modes. This relation depends on two wave numbers, $k'$ in the soft domain and $k$ in the hard domain. We show that it holds up to a good accuracy, when $k'/k\ll 1$ and $k'$ is in the linear regime, while the hard mode $k$ goes from linear ($0.1\,h\mathrm{Mpc}^{-1}$) to nonlinear ($1.0\,h\mathrm{Mpc}^{-1}$) scales. On scales $k\lesssim 0.4\,h\mathrm{Mpc}^{-1}$, we confirm the relation within the statistical error of the simulations (typically a few percent depending on the wave number), even though the bispectrum can already deviate from leading-order perturbation theory by more than $30\%$. We further examine the relation on smaller scales with higher resolution simulations. We find that the relation holds within the statistical error of the simulations at $z=1$, whereas we find deviations as large as $\sim 7\%$ at $k \sim 1.0\,h\mathrm{Mpc}^{-1}$ at $z=0.35$. We show that this can be explained partly by the breakdown of the approximation $\Omega_\mathrm{m}/f^2\simeq1$ with supplemental simulations done in the Einstein-de Sitter background cosmology. We also estimate the impact of this approximation on the power spectrum and bispectrum.


title: New Type IIB Backgrounds and Aspects of Their Field Theory Duals

abstract: In this paper we study aspects of geometries in Type IIA and Type IIB String theory and elaborate on their field theory dual pairs. The backgrounds are associated with reductions to Type IIA of solutions with $G_2$ holonomy in eleven dimensions. We classify these backgrounds according to their G-structure, perform a non-Abelian T-duality on them and find new Type IIB configurations presenting {\it dynamical} $SU(2)$-structure. We study some aspects of the associated field theories defined by these new backgrounds. Various technical details are clearly spelled out.


title: Dynamic Optical Superlattices with Topological Bands

abstract: We introduce an all-optical approach to producing high-flux synthetic magnetic fields for neutral atoms or molecules by designing intrinsically time-periodic optical superlattices. A single laser source, modulated to generate two frequencies, suffices to create dynamically modulated interference patterns which have topological Floquet energy bands. In particular, we propose a simple laser setup that realizes a tight-binding model with uniform flux per plaquette and well-separated Chern bands. Our method relies only on the particles' scalar polarizability and far detuned light.


title: Chemo-Archaeological Downsizing in a Hierarchical Universe: Impact of a Top Heavy IGIMF

abstract: We make use of a semi-analytical model of galaxy formation to investigate the origin of the observed correlation between [a/Fe] abundance ratios and stellar mass in elliptical galaxies. We implement a new galaxy-wide stellar initial mass function (Top Heavy Integrated Galaxy Initial Mass Function, TH-IGIMF) in the semi-analytic model SAG and evaluate its impact on the chemical evolution of galaxies. The SFR-dependence of the slope of the TH-IGIMF is found to be key to reproducing the correct [a/Fe]-stellar mass relation. Massive galaxies reach higher [a/Fe] abundance ratios because they are characterized by more top-heavy IMFs as a result of their higher SFR. As a consequence of our analysis, the value of the minimum embedded star cluster mass and of the slope of the embedded cluster mass function, which are free parameters involved in the TH-IGIMF theory, are found to be as low as 5 solar masses and 2, respectively. A mild downsizing trend is present for galaxies generated assuming either a universal IMF or a variable TH-IGIMF. We find that, regardless of galaxy mass, older galaxies (with formation redshifts > 2) are formed in shorter time-scales (< 2 Gyr), thus achieving larger [a/Fe] values. Hence, the time-scale of galaxy formation alone cannot explain the slope of the [a/Fe]-galaxy mass relation, but is responsible for the big dispersion of [a/Fe] abundance ratios at fixed stellar mass.We further test the hyphothesis of a TH-IGIMF in elliptical galaxies by looking into mass-to-light ratios, and luminosity functions. Models with a TH-IGIMF are also favoured by these constraints. In particular, mass-to-light ratios agree with observed values for massive galaxies while being overpredicted for less massive ones; this overprediction is present regardless of the IMF considered.


title: Superconducting proximity effect in topological metals

abstract: Much interest in the superconducting proximity effect in three-dimensional (3D) topological insulators (TIs) has been driven by the potential to induce Majorana bound states at the interface. Most candidate materials for 3D TI, however, are bulk metals, with bulk states at the Fermi level coexisting with well-defined surface states exhibiting spin-momentum locking. In such topological metals, the proximity effect can differ qualitatively from that in TIs. By studying a model topological metal-superconductor (TM-SC) heterostructure within the Bogoliubov-de Gennes formalism, we show that the pair amplitude reaches the naked surface, unlike in a topological insulator-superconductor (TI-SC) heterostructure where it is confined to the interface. Furthermore, we predict vortex-bound-state spectra to contain a Majorana zero-mode localized at the naked surface, separated from the bulk vortex-bound-state spectra by a finite energy gap in such a TM-SC heterostructure. These naked-surface-bound modes are amenable to experimental observation and manipulation, presenting advantages of TM-SC over TI-SC.


title: A simple technique for combining simplified models and its application to direct stop production

abstract: The results of many LHC searches for supersymmetric particles are interpreted using simplified models, in which one fixes the masses and couplings of most sparticles then scans over a few remaining masses of interest. We present a new technique for combining multiple simplified models (that requires no additional simulation) thereby highlighting the utility and limitations of simplified models in general, and demonstrating a simple way of improving LHC search strategies. The technique is used to derive limits on the stop mass that are model independent, modulo some reasonably generic assumptions which are quantified precisely. We find that current ATLAS and CMS results exclude stop masses up to 340 GeV for neutralino masses up to 120 GeV, provided that the total branching ratio into channels other than top-neutralino and bottom-chargino is small, and that there is no mass difference smaller than 10 GeV in the mass spectrum. In deriving these limits we place upper bounds on the branching ratios for complete stop pair decay processes for many values of the stop, neutralino and chargino masses. These are available with this paper.


title: General Features of the Relaxation Dynamics of Interacting Quantum Systems

abstract: We study numerically and analytically isolated interacting quantum systems that are taken out of equilibrium instantaneously (quenched). The probability of finding the initial state in time, the so-called fidelity, decays fastest for systems described by full random matrices, where simultaneous many-body interactions are implied. In the realm of realistic systems with two-body interactions, the dynamics is slower and depends on the interplay between the initial state and the Hamiltonian characterizing the system. The fastest fidelity decay in this case is Gaussian and can persist until saturation. A simple general picture, in which the fidelity plays a central role, is also achieved for the short-time dynamics of few-body observables. It holds for initial states that are eigenstates of the observables. We also discuss the need to reassess analytical expressions that were previously proposed to describe the evolution of the Shannon entropy. Our analyses are mainly developed for initial states that can be prepared in experiments with cold atoms in optical lattices.


title: The four-loop remainder function and multi-Regge behavior at NNLLA in planar N=4 super-Yang-Mills theory

abstract: We present the four-loop remainder function for six-gluon scattering with maximal helicity violation in planar N=4 super-Yang-Mills theory, as an analytic function of three dual-conformal cross ratios. The function is constructed entirely from its analytic properties, without ever inspecting any multi-loop integrand. We employ the same approach used at three loops, writing an ansatz in terms of hexagon functions, and fixing coefficients in the ansatz using the multi-Regge limit and the operator product expansion in the near-collinear limit. We express the result in terms of multiple polylogarithms, and in terms of the coproduct for the associated Hopf algebra. From the remainder function, we extract the BFKL eigenvalue at next-to-next-to-leading logarithmic accuracy (NNLLA), and the impact factor at NNNLLA. We plot the remainder function along various lines and on one surface, studying ratios of successive loop orders. As seen previously through three loops, these ratios are surprisingly constant over large regions in the space of cross ratios, and they are not far from the value expected at asymptotically large orders of perturbation theory.


title: Privacy and National Security Issues in Social Networks: The Challenges

abstract: Online social networks are becoming a major growth point of the internet, as individuals, companies and governments constantly desire to interact with one another, the ability of the internet to deliver this networking capabilities grows stronger. In this paper, we looked at the structure and components of the member profile and the challenges of privacy issues faced by individuals and governments that participate in social networking. We also looked at how it can be used to distort national security, how it became the new weapons of mass mobilization and also how social networks have became the rallying forces for revolutions and social justice.


title: Detection of Stacked Filament Lensing Between SDSS Luminous Red Galaxies

abstract: We search for the lensing signal of massive filaments between 135,000 pairs of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey. We develop a new estimator that cleanly removes the much larger shear signal of the neighboring LRG halos, relying only on the assumption of spherical symmetry. We consider two models: a thick filament model constructed from ray-tracing simulations for $\Lambda$CDM model, and a thin filament model which models the filament by a string of halos along the line connecting the two LRGs. We show that the filament lensing signal is in nice agreement with the thick simulation filament, while strongly disfavoring the thin model. The magnitude of the lensing shear due to the filament is below $10^{-4}$. Employing the likelihood ratio test, we find a 4.5$\sigma$ significance for the detection of the filament lensing signal, corresponding to a null hypothesis fluctuation probability of $3 \times 10^{-6}$. We also carried out several null tests to verify that the residual shear signal from neighboring LRGs and other shear systematics are minimized.


title: Novel Computational Methods for High-Dimensional Stochastic Sensitivity Analysis

abstract: This paper presents three new computational methods for calculating design sensitivities of statistical moments and reliability of high-dimensional complex systems subject to random input. The first method represents a novel integration of polynomial dimensional decomposition (PDD) of a multivariate stochastic response function and score functions. Applied to the statistical moments, the method provides mean-square convergent analytical expressions of design sensitivities of the first two moments of a stochastic response. The second and third methods, relevant to probability distribution or reliability analysis, exploit two distinct combinations built on PDD: the PDD-SPA method, entailing the saddlepoint approximation (SPA) and score functions; and the PDD-MCS method, utilizing the embedded Monte Carlo simulation (MCS) of the PDD approximation and score functions. For all three methods developed, the statistical moments or failure probabilities and their design sensitivities are both determined concurrently from a single stochastic analysis or simulation. Numerical examples, including a 100-dimensional mathematical problem, indicate that the new methods developed provide not only theoretically convergent or accurate design sensitivities, but also computationally efficient solutions. A practical example involving robust design optimization of a three-hole bracket illustrates the usefulness of the proposed methods.


title: The Fates of Mutant Lineages and the Distribution of Fitness Effects of Beneficial Mutations in Laboratory Budding Yeast Populations

abstract: The outcomes of evolution are determined by which mutations occur and fix. In rapidly adapting microbial populations, this process is particularly hard to predict because lineages with different beneficial mutations often spread simultaneously and interfere with one another's fixation. Hence to predict the fate of any individual variant, we must know the rate at which new mutations create competing lineages of higher fitness. Here, we directly measured the effect of this interference on the fates of specific adaptive variants in laboratory Saccharomyces cerevisiae populations and used these measurements to infer the distribution of fitness effects of new beneficial mutations. To do so, we seeded marked lineages with different fitness advantages into replicate populations and tracked their subsequent frequencies for hundreds of generations. Our results illustrate the transition between strongly advantageous lineages which decisively sweep to fixation and more moderately advantageous lineages that are often outcompeted by new mutations arising during the course of the experiment. We developed an approximate likelihood framework to compare our data to simulations and found that the effects of these competing beneficial mutations were best approximated by an exponential distribution, rather than one with a single effect size. We then used this inferred distribution of fitness effects to predict the rate of adaptation in a set of independent control populations. Finally, we discuss how our experimental design can serve as a screen for rare, large-effect beneficial mutations.


title: Real-Time Notification for Resource Synchronization

abstract: Web applications frequently leverage resources made available by remote web servers. As resources are created, updated, deleted, or moved, these applications face challenges to remain in lockstep with the server's change dynamics. Several approaches exist to help meet this challenge for use cases where "good enough" synchronization is acceptable. But when strict resource coverage or low synchronization latency is required, commonly accepted Web-based solutions remain elusive. This paper details characteristics of an approach that aims at decreasing synchronization latency while maintaining desired levels of accuracy. The approach builds on pushing change notifications and pulling changed resources and it is explored with an experiment based on a DBpedia Live instance.


title: Homotopy Representations over the Orbit Category

abstract: Let G be a finite group. The unit sphere in a finite-dimensional orthogonal G-representation motivates the definition of homotopy representations, due to tom Dieck. We introduce an algebraic analogue, and establish its basic properties including the Borel-Smith conditions and realization by finite G-CW-complexes.


title: Space-time S-matrix and Flux-tube S-matrix III. The two-particle contributions

abstract: We consider light-like Wilson loops with hexagonal geometry in the planar limit of N=4 Super-Yang-Mills theory. Within the Operator-Product-Expansion framework these loops receive contributions from all states that can propagate on top of the colour flux tube sourced by any two opposite edges of the loops. Of particular interest are the two-particle contributions. They comprise virtual effects like the propagation of a pair of scalars, fermions, and gluons, on top of the flux tube. Each one of them is thoroughly discussed in this paper. Our main result is the prediction of all the twist-2 corrections to the expansion of the dual 6-gluons MHV amplitude in the near-collinear limit at finite coupling. At weak coupling, our result was recently used by Dixon, Drummond, Duhr and Pennington to predict the full amplitude at four loops. At strong coupling, it allows us to make contact with the classical string description and to recover the (previously elusive) AdS(3) mode from the continuum of two-fermion states. More generally, the two-particle contributions serve as an exemplar for all the multi-particle corrections.


title: On the evolution of the CO snow line in protoplanetary disks

abstract: CO is thought to be a vital building block for prebiotic molecules that are necessary for life. Thus, understanding where CO existed in a solid phase within the solar nebula is important for understanding the origin of life. We model the evolution of the CO snow line in a protoplanetary disk. We find that the current observed location of the CO snow line in our solar system, and in the solar system analogue TW Hydra, cannot be explained by a fully turbulent disk model. With time-dependent disk models we find that the inclusion of a dead zone (a region of low turbulence) can resolve this problem. Furthermore, we obtain a fully analytic solution for the CO snow line radius for late disk evolutionary times. This will be useful for future observational attempts to characterize the demographics and predict the composition and habitability of exoplanets.


title: Ballistic performance comparison of monolayer transition metal dichalcogenide MX2 (M = Mo, W; X = S, Se, Te) MOSFETs

abstract: We study the transport properties of monolayer MX2 (M = Mo, W; X = S, Se, Te) n- and p- channel metal-oxide-semiconductor field effect transistors (MOSFETs) using full-band ballistic non-equilibrium Green's function simulations with an atomistic tight-binding Hamiltonian with hopping potentials obtained from density functional theory. We discuss the subthreshold slope, drain-induced barrier lowering (DIBL), as well as gate-induced drain leakage (GIDL) for different monolayer MX2 MOSFETs. We also report the possibility of negative differential resistance behavior in the output characteristics of nanoscale monolayer MX2 MOSFETs.


title: An Energy-Minimization Finite-Element Approach for the Frank-Oseen Model of Nematic Liquid Crystals: Continuum and Discrete Analysis

abstract: This paper outlines an energy-minimization finite-element approach to the computational modeling of equilibrium configurations for nematic liquid crystals under free elastic effects. The method targets minimization of the system free energy based on the Frank-Oseen free-energy model. Solutions to the intermediate discretized free elastic linearizations are shown to exist generally and are unique under certain assumptions. This requires proving continuity, coercivity, and weak coercivity for the accompanying appropriate bilinear forms within a mixed finite-element framework. Error analysis demonstrates that the method constitutes a convergent scheme. Numerical experiments are performed for problems with a range of physical parameters as well as simple and patterned boundary conditions. The resulting algorithm accurately handles heterogeneous constant coefficients and effectively resolves configurations resulting from complicated boundary conditions relevant in ongoing research.


title: Why do we change whatever amount we found in the first envelope: the Wikipedia "two envelopes problem" commented

abstract: We analyze the main arguments that attempt to explain why there is no point in changing the envelope. Most people confuse estimation and calculation, conditional and unconditional probabilities, random and non-random variables, modelling and reality, chance and uncertainty (or credence), whatever X known and whatever X unknown, personal and interpersonal estimates. Such confusions suggest that there would be no point in swapping the envelope. Our analysis is that if the first envelope is opened, whatever the amount, it is rational and consistent to think that it is worth changing it. If the first envelope is not opened, it is rational and consistent to believe there is no point in swapping the envelope.


title: Veronese Geometry and the Electroweak Vacuum Moduli Space

abstract: We explain the origin of the Veronese surface in the vacuum moduli space geometry of the MSSM electroweak sector. While this result appeared many years ago using techniques of computational algebraic geometry, it has never been demonstrated analytically. Here, we present an analytical derivation of the vacuum geometry of the electroweak theory by understanding how the F- and D-term relations lead to the Veronese surface. We moreover give a detailed description of this geometry, realising an extra branch as a zero-dimensional point when quadratic Higgs lifting deformations are incorporated into the superpotential.


title: Ground state and excitation properties of soft-core bosons

abstract: We study the physics of soft-core bosons at zero temperature in two dimensions for a class of potentials that could be realised in experiments with Rydberg dressed Bose-Einstein condensates. We analyze the ground state properties of the system in detail and provide a complete description of the excitation spectra in both superfluid, supersolid and crystalline phase for a wide range of interaction strengths and densities. In addition we describe a method to extract the transverse gapless excitation modes in the phases with broken translational symmetry within the framework of path integral Monte Carlo methods.


title: Distributed synthesis for acyclic architectures

abstract: The distributed synthesis problem is about constructing cor- rect distributed systems, i.e., systems that satisfy a given specification. We consider a slightly more general problem of distributed control, where the goal is to restrict the behavior of a given distributed system in order to satisfy the specification. Our systems are finite state machines that communicate via rendez-vous (Zielonka automata). We show decidability of the synthesis problem for all omega-regular local specifications, under the restriction that the communication graph of the system is acyclic. This result extends a previous decidability result for a restricted form of local reachability specifications.


title: On the spheroidized carbide dissolution and elemental partitioning in a high carbon bearing steel 100Cr6

abstract: We report on the characterization of high carbon bearing steel 100Cr6 using electron microscopy and atom probe tomography in combination with multi-component diffusion simulations (DICTRA). Scanning electron micrographs show that around 14 vol.% spheroidized carbides are formed during soft annealing and only 3 vol.% remain after dissolution into the austenitic matrix by austenitization at 1123 K (850 {\deg}C) for 300 s. The spheroidized particles are identified as (Fe, Cr)3C by transmission electron microscopy. Atom probe analyses reveal the redistribution and partitioning behaviors of elements, i.e. C, Si, Mn, Cr, Fe in both, the spheroidized carbides and the bainitic matrix in the sample isothermally heat-treated at 773 K (500 {\deg}C) after austenitization. A homogeneous distribution of C and gradual gradient of Cr was detected within the spheroidized carbides. Due to its limited diffusivity in (Fe, Cr)3C, Cr exhibits a maximum concentration at the surface of spheroidized carbides (16 at.%) and decreases gradually from surface towards the core down to a level of about 2 at.%. The atom probe results also indicate that the partially dissolved spheroidized carbides during austenitization may serve as nucleation sites for intermediate temperature cementite within bainite, which results in a relatively softer surface and harder core in spheroidized particles. This microstructure may contribute to the good wear resistance and fatigue properties


title: More on loops in reheating: Non-gaussianities and tensor power spectrum

abstract: We consider the single field chaotic m^2\phi^2 inflationary model with a period of preheating, where the inflaton decays to another scalar field \chi\ in the parametric resonance regime. In a recent work, one of us has shown that the \chi\ modes circulating in the loops during preheating notably modify the <\zeta\zeta> correlation function. We first rederive this result using a different gauge condition hence reconfirm that superhorizon \zeta\ modes are affected by the loops in preheating. Further, we examine how \chi\ loops give rise to non-gaussianity and affect the tensor perturbations. For that, all cubic and some higher order interactions involving two \chi\ fields are determined and their contribution to the non-gaussianity parameter f_{NL} and the tensor power spectrum are calculated at one loop. Our estimates for these corrections show that while a large amount of non-gaussianity can be produced during reheating, the tensor power spectrum receive moderate corrections. We observe that the loop quantum effects increase with more \chi\ fields circulating in the loops indicating that the perturbation theory might be broken down. These findings demonstrate that the loop corrections during reheating are significant and they must be taken into account for precision inflationary cosmology.


title: Multiple Window Moving Horizon Estimation

abstract: Long horizon lengths in Moving Horizon Estimation are desirable to reach the performance limits of the full information estimator. However, the conventional MHE technique suffers from a number of deficiencies in this respect. First, the problem complexity scales at least linearly with the horizon length selected, which restrains from selecting long horizons if computational limitations are present. Second, there is no monitoring of constraint activity/inactivity which results in conducting redundant constrained minimizations even when no constraints are active. In this study we develop a Multiple-Window Moving Horizon Estimation strategy (MW-MHE) that exploits constraint inactivity to reduce the problem size in long horizon estimation problems. The arrival cost is approximated using the unconstrained full information estimator arrival cost to guarantee stability of the technique. A new horizon length selection criteria is developed based on maximum sensitivity between remote states in time. The development will be in terms of general causal descriptor systems, which includes the standard state space representation as a special case. The potential of the new estimation algorithm will be demonstrated with an example showing a significant reduction in both computation time and numerical errors compared to conventional MHE.


title: Element-Resolved Corrosion Analysis of Stainless-Type Glass-Forming Steels

abstract: Ultrathin passive films effectively prevent the chemical attack of stainless steel grades in corrosive environments; their stability critically depends on the interplay between structure and chemistry of the constituents Fe-Cr-Mo. In particular, nanoscale inhomogeneities along the surface can have a tremendous impact on material failure, but are yet barely understood. Addressing a stainless-type glass-forming Fe50Cr15Mo14C15B6 alloy and utilizing a combination of complementary high-resolution analytical techniques, we relate near-atomistic insight into different gradual nanostructures with time- and element-resolved dissolution behavior. The progressive elemental segregation on the nanoscale is followed in its influence on the concomitant degree of passivity. A detrimental transition from Cr-controlled passivity to Mo-controlled breakdown is dissected atom-by-atom demonstrating the importance of nanoscale knowledge for understanding corrosion.


title: Flow-based reputation with uncertainty: Evidence-Based Subjective Logic

abstract: The concept of reputation is widely used as a measure of trustworthiness based on ratings from members in a community. The adoption of reputation systems, however, relies on their ability to capture the actual trustworthiness of a target. Several reputation models for aggregating trust information have been proposed in the literature. The choice of model has an impact on the reliability of the aggregated trust information as well as on the procedure used to compute reputations. Two prominent models are flow-based reputation (e.g., EigenTrust, PageRank) and Subjective Logic based reputation. Flow-based models provide an automated method to aggregate trust information, but they are not able to express the level of uncertainty in the information. In contrast, Subjective Logic extends probabilistic models with an explicit notion of uncertainty, but the calculation of reputation depends on the structure of the trust network and often requires information to be discarded. These are severe drawbacks. In this work, we observe that the `opinion discounting' operation in Subjective Logic has a number of basic problems. We resolve these problems by providing a new discounting operator that describes the flow of evidence from one party to another. The adoption of our discounting rule results in a consistent Subjective Logic algebra that is entirely based on the handling of evidence. We show that the new algebra enables the construction of an automated reputation assessment procedure for arbitrary trust networks, where the calculation no longer depends on the structure of the network, and does not need to throw away any information. Thus, we obtain the best of both worlds: flow-based reputation and consistent handling of uncertainties.


title: Rate of Convergence for Large Coupling Limits in Sobolev Spaces

abstract: We estimate the rate of convergence, in the so-called large coupling limit, for Schr\"odinger type operators on bounded domains. The Schr\"odinger we deal with have "interaction potentials" supported in a compact inclusion. We show that if the boundary of the inclusion is sufficiently smooth, one essentially recovers the "free Hamiltonian" in the exterior domain with Dirichlet boundary conditions. In addition, we obtain a convergence rate, in $L^2$, that is $\mathcal{O}(\lambda^{-\frac{1}{4}})$ where $\lambda$ is the coupling parameter. Our methods include energy estimates, trace estimates, interpolation and duality.


title: On the entanglement of formation of two-mode Gaussian states: a compact form

abstract: We define an EPR-like uncertainty by using the Duan et al.'s inequality which gives a necessary and sufficient condition for the separability of any two-mode Gaussian state. We show that for a given amount of entanglement, the uncertainty is minimized by pure two-mode squeezed states. Using this fact, we write the optimal pure-state decomposition and derive a compact form for the entanglement of formation of two-mode Gaussian states. For illustration purposes, we consider symmetric and squeezed thermal states as special cases and evaluate their entanglement of formation explicitly. For the symmetric states, our result is in agreement with the Giedke et al.'s one. To our knowledge, our work is the first one which gives the exact entanglement of formation of two-mode squeezed thermal states explicitly.


title: $p$-adic quasi Gibbs measure for Vannimenus model on a Cayley tree

abstract: In this paper we study the boundedness of the $p$-adic quasi Gibbs measures for the Vannimenus model on a Cayley tree of order two.


title: Pre-mainsequence stellar evolution in N-body models

abstract: We provide a set of analytic fits to the radii of pre-mainsequence stars in the mass range 0.1 < M/Msun < 8.0. We incorporate the formulae in N-body cluster models for evolution from the beginning of pre-main sequence. In models with 1,000 stars and high initial cluster densities, pre-mainsequence evolution causes roughly twice the number of collisions between stars than in similar models with evolution begun only from the zero-age main sequence. The collisions are often all part of a runaway sequence that creates one relatively massive star.


title: Covariance among independent variables determines the overfitting and underfitting problems in variation partitioning methods: with a special focus on the mixed co-variation

abstract: The effectiveness and validity of applying variation partitioning methods in community ecology has been questioned. Here, using mathematical deduction and numerical simulation, we made an attempt to uncover the underlying mechanisms determining the effectiveness of variation partitioning techniques. The covariance among independent variables determines the under-fitting and over-fitting problem with the variation partitioning process. Ideally, it is assumed that the covariance among independent variables will be zero (no correlation at all), however, typically there will be some colinearities. Therefore, we analyzed the role of slight covariance on influencing species variation partitioning. We concluded that when the covariance between spatial and environmental predictors is positive, all the three components-pure environmental, spatial variations and mixed covariation were over-fitted, with the sign of the true covariation being negative. In contrast, when the covariance is negative, all the three components were under-fitted with the sign of true covariation being positive. Other factors, including extra noise levels, the strengths of variable coefficients and the patterns of landscape gradients, could reduce the fitting problems caused by the covariance of variables. The conventional calculation of mixed covariation is incorrect and misleading, as the true and estimated covariations are always sign-opposite. In conclusion, I challenge the conventional three-step procedure of variation partitioning, suggesting that a full regression model with all variables together is robust enough to correctly partition variations.


title: A Keplerian-based Hamiltonian Splitting for Gravitational $N$-body Simulations

abstract: We developed a Keplerian-based Hamiltonian splitting for solving the gravitational $N$-body problem. This splitting allows us to approximate the solution of a general $N$-body problem by a composition of multiple, independently evolved $2$-body problems. While the Hamiltonian splitting is exact, we show that the composition of independent $2$-body problems results in a non-symplectic non-time-symmetric first-order map. A time-symmetric second-order map is then constructed by composing this basic first-order map with its self-adjoint. The resulting method is precise for each individual $2$-body solution and produces quick and accurate results for near-Keplerian $N$-body systems, like planetary systems or a cluster of stars that orbit a supermassive black hole. The method is also suitable for integration of $N$-body systems with intrinsic hierarchies, like a star cluster with primordial binaries. The superposition of Kepler solutions for each pair of particles makes the method excellently suited for parallel computing; we achieve $\gtrsim 64\%$ efficiency for only $8$ particles per core, but close to perfect scaling for $16384$ particles on a $128$ core distributed-memory computer. We present several implementations in \texttt{Sakura}, one of which is publicly available via the AMUSE framework.


title: On presimplifiable group rings

abstract: A ring A is called presimplifiable if whenever a; b belongs to A and a = ab, then either a = 0 or b is a unit in A. Let A be a commutative ring and G be an abelian torsion group. For the group ring A[G], we prove that A[G] is presimplifiable if and only if A is presimplifiable and G is a p-group with p belongs to the Jacobson radical of A, and it is shown that A[G] is domainlike (i.e all zero divisors are nilpotents) if and only if A is domainlike and G is a p-group and p is a nilpotent in A. Furthermore, whenever the group ring A[G] is presimplifiable we prove that A[H] is presimplifiable for any subgroup H of G. Also, for a torsion free group G we prove that A[G] is domainlike if and only if A[G] is integral domain.


title: Lie Group Spectral Variational Integrators

abstract: We present a new class of high-order variational integrators on Lie groups. We show that these integrators are symplectic, momentum preserving, and can be constructed to be of arbitrarily high-order, or can be made to converge geometrically. Furthermore, these methods are stable and accurate for very large time steps. We demonstrate the construction of one such variational integrator for the rigid body, and discuss how this construction could be generalized to other related Lie group problems. We close with several numerical examples which demonstrate our claims, and discuss further extensions of our work.


title: The starbursts in the Milky Way

abstract: High-mass stars are major players in the chemical and dynamical evolution of galaxies, and young massive clusters are the natural laboratories to study their evolution and their impact on star formation processes. Only in recent years have we become aware of the existence of numerous massive (M_cl > 10000M_sun) clusters in our Galaxy. Here I give a review, rather biased towards my own research interests, of the observational and theoretical efforts that have led to a description of their properties, and present an overview of the two (perhaps three) starburst regions known outside the Galactic Centre neighbourhood: the Scutum Complex, its putative counterpart on the far side of the Long Bar, and the starburst cluster Westerlund 1.


title: Differential Privacy: An Economic Method for Choosing Epsilon

abstract: Differential privacy is becoming a gold standard for privacy research; it offers a guaranteed bound on loss of privacy due to release of query results, even under worst-case assumptions. The theory of differential privacy is an active research area, and there are now differentially private algorithms for a wide range of interesting problems. However, the question of when differential privacy works in practice has received relatively little attention. In particular, there is still no rigorous method for choosing the key parameter $\epsilon$, which controls the crucial tradeoff between the strength of the privacy guarantee and the accuracy of the published results. In this paper, we examine the role that these parameters play in concrete applications, identifying the key questions that must be addressed when choosing specific values. This choice requires balancing the interests of two different parties: the data analyst and the prospective participant, who must decide whether to allow their data to be included in the analysis. We propose a simple model that expresses this balance as formulas over a handful of parameters, and we use our model to choose $\epsilon$ on a series of simple statistical studies. We also explore a surprising insight: in some circumstances, a differentially private study can be more accurate than a non-private study for the same cost, under our model. Finally, we discuss the simplifying assumptions in our model and outline a research agenda for possible refinements.


title: Adaptive-Sparse Polynomial Dimensional Decomposition Methods for High-Dimensional Stochastic Computing

abstract: This article presents two novel adaptive-sparse polynomial dimensional decomposition (PDD) methods for solving high-dimensional uncertainty quantification problems in computational science and engineering. The methods entail global sensitivity analysis for retaining important PDD component functions, and a full- or sparse-grid dimension-reduction integration or quasi Monte Carlo simulation for estimating the PDD expansion coefficients. A unified algorithm, endowed with two distinct ranking schemes for grading component functions, was created for their numerical implementation. The fully adaptive-sparse PDD method is comprehensive and rigorous, leading to the second-moment statistics of a stochastic response that converges to the exact solution when the tolerances vanish. A partially adaptive-sparse PDD method, obtained through regulated adaptivity and sparsity, is economical and is, therefore, expected to solve practical problems with numerous variables. Compared with past developments, the adaptive-sparse PDD methods do not require its truncation parameter(s) to be assigned \emph{a priori} or arbitrarily. The numerical results reveal that an adaptive-sparse PDD method achieves a desired level of accuracy with considerably fewer coefficients compared with existing PDD approximations. For a required accuracy in calculating the probabilistic response characteristics, the new bivariate adaptive-sparse PDD method is more efficient than the existing bivariately truncated PDD method by almost an order of magnitude. Finally, stochastic dynamic analysis of a disk brake system was performed, demonstrating the ability of the new methods to tackle practical engineering problems.


title: L-infinity Norm Design of Linear-phase Robust Broadband Beamformers using Constrained Optimization

abstract: A new method for the design of linear-phase robust far-field broadband beamformers using constrained optimization is proposed. In the method, the maximum passband ripple and minimum stopband attenuation are ensured to be within prescribed levels, while at the same time maintaining a good linear-phase characteristic at a prescribed group delay in the passband. Since the beamformer is intended primarily for small-sized microphone arrays where the microphone spacing is small relative to the wavelength at low frequencies, the beamformer can become highly sensitive to spatial white noise and array imperfections if a direct minimization of the error is performed. Therefore, to limit the sensitivity of the beamformer the optimization is carried out by constraining a sensitivity parameter, namely, the white noise gain (WNG) to be above prescribed levels across the frequency band. Two novel design variants have been developed. The first variant is formulated as a convex optimization problem where the maximum error in the passband is minimized, while the second variant is formulated as an iterative optimization problem and has the advantage of significantly improving the linear-phase characteristics of the beamformer under any prescribed group delay or linear-array configuration. In the second variant, the passband group-delay deviation is minimized while ensuring that the maximum passband ripple and stopband attenuation are within prescribed levels. To reduce the computational effort in carrying out the optimization, a nonuniform variable sampling approach over the frequency and angular dimensions is used to compute the required parameters. Experiment results show that beamformers designed using the proposed methods have much smaller passband group-delay deviation for similar passband ripple and stopband attenuation than a modified version of an existing method.


title: Elements of Trust in Named-Data Networking

abstract: In contrast to today's IP-based host-oriented Internet architecture, Information-Centric Networking (ICN) emphasizes content by making it directly addressable and routable. Named Data Networking (NDN) architecture is an instance of ICN that is being developed as a candidate next-generation Internet architecture. By opportunistically caching content within the network (in routers), NDN appears to be well-suited for large-scale content distribution and for meeting the needs of increasingly mobile and bandwidth-hungry applications that dominate today's Internet. One key feature of NDN is the requirement for each content object to be digitally signed by its producer. Thus, NDN should be, in principle, immune to distributing fake (aka "poisoned") content. However, in practice, this poses two challenges for detecting fake content in NDN routers: (1) overhead due to signature verification and certificate chain traversal, and (2) lack of trust context, i.e., determining which public keys are trusted to verify which content. Because of these issues, NDN does not force routers to verify content signatures, which makes the architecture susceptible to content poisoning attacks. This paper explores root causes of, and some cures for, content poisoning attacks in NDN. In the process, it becomes apparent that meaningful mitigation of content poisoning is contingent upon a network-layer trust management architecture, elements of which we construct while carefully justifying specific design choices. This work represents the initial effort towards comprehensive trust management for NDN.


title: Finding Eigenvalues the Rupert Way

abstract: We develop a stability index for the travelling waves of non-linear reaction diffusion equations using the geometric phase induced on the Hopf bundle $S^{2n-1} \subset \mathbb{C}^n$. This can be viewed as an alternative formulation of the winding number calculation of the Evans function, whose zeroes correspond to the eigenvalues of the linearization of reaction diffusion operators about the wave. The stability of a travelling wave can be determined by the existence of eigenvalues of positive real part for the linear operator. Our method of geometric phase for locating and counting eigenvalues is inspired by the numerical results in Way's Dynamics in the Hopf bundle, the geometric phase and implications for dynamical systems. We provide a detailed proof of the relationship between the phase and eigenvalues for dynamical systems defined on $\mathbb{C}^2$ and sketch the proof of the method of geometric phase for $\mathbb{C}^n$ and its generalization to boundary-value problems. Implementing the numerical method, modified from Way's work, we conclude with open questions inspired from the results.


title: A minimal-variable symplectic integrator on spheres

abstract: We construct a symplectic, globally defined, minimal-coordinate, equivariant integrator on products of 2-spheres. Examples of corresponding Hamiltonian systems, called spin systems, include the reduced free rigid body, the motion of point vortices on a sphere, and the classical Heisenberg spin chain, a spatial discretisation of the Landau-Lifschitz equation. The existence of such an integrator is remarkable, as the sphere is neither a vector space, nor a cotangent bundle, has no global coordinate chart, and its symplectic form is not even exact. Moreover, the formulation of the integrator is very simple, and resembles the geodesic midpoint method, although the latter is not symplectic.


title: Resistive flow in a weakly interacting Bose-Einstein condensate

abstract: We report the direct observation of resistive flow through a weak link in a weakly interacting atomic Bose-Einstein condensate. Two weak links separate our ring-shaped superfluid atomtronic circuit into two distinct regions, a source and a drain. Motion of these weak links allows for creation of controlled flow between the source and the drain. At a critical value of the weak link velocity, we observe a transition from superfluid flow to superfluid plus resistive flow. Working in the hydrodynamic limit, we observe a conductivity that is 4 orders of magnitude larger than previously reported conductivities for a Bose-Einstein condensate with a tunnel junction. Good agreement with zero-temperature Gross-Pitaevskii simulations and a phenomenological model based on phase slips indicate that the creation of excitations plays an important role in the resulting conductivity. Our measurements of resistive flow elucidate the microscopic origin of the dissipation and pave the way for more complex atomtronic devices.


title: Permutation Patterns in Latin Squares

abstract: In this paper we study pattern avoidance in Latin Squares, which gives us a two dimensional analogue of the well studied notion of pattern avoidance in permutations. Our main results include enumerating and characterizing the Latin Squares which avoid patterns of length three and a generalization of the Erd\H{o}s-Szekeres theorem. We also discuss equivalence classes among longer patterns, and conclude by describing open questions of interest both in light of pattern avoidance and their potential to reveal information about the structure of Latin Squares. Along the way, we show that classical results need not trivially generalize, and demonstrate techniques that may help answer future questions.


title: Zero-bias autoencoders and the benefits of co-adapting features

abstract: Regularized training of an autoencoder typically results in hidden unit biases that take on large negative values. We show that negative biases are a natural result of using a hidden layer whose responsibility is to both represent the input data and act as a selection mechanism that ensures sparsity of the representation. We then show that negative biases impede the learning of data distributions whose intrinsic dimensionality is high. We also propose a new activation function that decouples the two roles of the hidden layer and that allows us to learn representations on data with very high intrinsic dimensionality, where standard autoencoders typically fail. Since the decoupled activation function acts like an implicit regularizer, the model can be trained by minimizing the reconstruction error of training data, without requiring any additional regularization.


title: On non-local representations of the ageing algebra in $d\geq 1$ dimensions

abstract: Non-local representations of the ageing algebra for generic dynamical exponents $z$ and for any space dimension $d\geq 1$ are constructed. The mechanism for the closure of the Lie algebra is explained. The Lie algebra generators contain higher-order differential operators or the Riesz fractional derivative. Co-variant two-time response functions are derived. An application to phase-separation in the conserved spherical model is described.


title: Young Star Clusters In The Circumnuclear Region Of NGC 2110

abstract: High-resolution observations in the near infrared show star clusters around the active galactic nucleus (AGN) of the Seyfert 1 NGC2110, along with a 90 x 35 pc bar of shocked gas material around its nucleus. These are seen for the first time in our imaging and gas kinematics of the central 100pc with the Keck OSIRIS instrument with adaptive optics. Each of these clusters is 2-3 times brighter than the Arches cluster close to the centre of the Milky Way. The core star formation rate (SFR) is 0.3 M$_\odot$/yr. The photoionized gas (He I) dynamics imply an enclosed mass of 3-4 x 10$^8$ M$_\odot$. These observations demonstrate the physical linkage between AGN feedback, which triggers star formation in massive clusters, and the resulting stellar (and SNe) winds, which cause the observed [Fe II] emission and feed the black hole.


title: Inferring microscopic kinetics of a Markov process using maximum caliber

abstract: We present a principled approach for estimating the matrix of microscopic rates among states of a Markov process, given only its stationary state population distribution and a single average global kinetic observable. We adapt Maximum Caliber, a variational principle in which a path entropy is maximized over the distribution of all the possible trajectories, subject to basic kinetic constraints and some average dynamical observables. We show that this approach leads, under appropriate conditions, to the continuous-time master equation and a Smoluchowski-like equation that is valid for both equilibrium and non-equilibrium stationary states. We illustrate the method by computing the solvation dynamics of water molecules from molecular dynamics trajectories.


title: On the Spectrum of Wenger Graphs

abstract: Let $q=p^e$, where $p$ is a prime and $e\geq 1$ is an integer. For $m\geq 1$, let $P$ and $L$ be two copies of the $(m+1)$-dimensional vector spaces over the finite field $\mathbb{F}_q$. Consider the bipartite graph $W_m(q)$ with partite sets $P$ and $L$ defined as follows: a point $(p)=(p_1,p_2,\ldots,p_{m+1})\in P$ is adjacent to a line $[l]=[l_1,l_2,\ldots,l_{m+1}]\in L$ if and only if the following $m$ equalities hold: $l_{i+1} + p_{i+1}=l_{i}p_1$ for $i=1,\ldots, m$. We call the graphs $W_m(q)$ Wenger graphs. In this paper, we determine all distinct eigenvalues of the adjacency matrix of $W_m(q)$ and their multiplicities. We also survey results on Wenger graphs.


title: Self-Completeness in Alternative Theories of Gravity

abstract: It has recently been shown via an equivalence of gravitational radius and Compton wavelength in four dimensions that the trans-Planckian regime of gravity may by semi-classical, and that this point is defined by a minimum horizon radius commensurate with the Planck mass. We extend the formalism to gravity in the context of Randall-Sundrum and the generalized uncertainty principle.


title: Regeneration tests of a room temperature magnetic refrigerator and heat pump

abstract: A magnetic heat pump apparatus consisting of a solid magnetic refrigerant, gadolinium, and a liquid regenerator column of ethanol and water has been tested. Utilizing a 7T field, it produced a maximum temperature span of 80 K, and in separate tests, a lowest temperature of 241 K and a highest temperature of 328 K. Thermocouples, placed at intervals along the regenerator tube, permitted measurement of the temperature distribution in the regenerator fluid. No attempt was made to extract refrigeration from the device, but analysis of the temperature distributions shows that 34 watts of refrigeration was produced.


title: Intrinsically Motivated Learning of Visual Motion Perception and Smooth Pursuit

abstract: We extend the framework of efficient coding, which has been used to model the development of sensory processing in isolation, to model the development of the perception/action cycle. Our extension combines sparse coding and reinforcement learning so that sensory processing and behavior co-develop to optimize a shared intrinsic motivational signal: the fidelity of the neural encoding of the sensory input under resource constraints. Applying this framework to a model system consisting of an active eye behaving in a time varying environment, we find that this generic principle leads to the simultaneous development of both smooth pursuit behavior and model neurons whose properties are similar to those of primary visual cortical neurons selective for different directions of visual motion. We suggest that this general principle may form the basis for a unified and integrated explanation of many perception/action loops.


title: Nonlinearity management of photonic composites and observation of spatial-modulation instability due to quintic nonlinearity

abstract: We present a procedure for nonlinearity management of metal-dielectric composites. Varying the volume fraction occupied by silver nanoparticles suspended in acetone we could cancel the refractive index related to the third-order susceptibility, $\chi_{eff}^{(3)}$, and the nonlinear refraction behavior was due to the fifth-order susceptibility, $\chi_{eff}^{(5)}$. Hence, in a cross-phase modulation experiment, we demonstrated for the first time the effect of spatial-modulation- instability due to $\chi_{eff}^{(5)}$. The results are corroborated with numerical calculations based on a generalized Maxwell-Garnet model.


title: Geometry and Expressive Power of Conditional Restricted Boltzmann Machines

abstract: Conditional restricted Boltzmann machines are undirected stochastic neural networks with a layer of input and output units connected bipartitely to a layer of hidden units. These networks define models of conditional probability distributions on the states of the output units given the states of the input units, parametrized by interaction weights and biases. We address the representational power of these models, proving results their ability to represent conditional Markov random fields and conditional distributions with restricted supports, the minimal size of universal approximators, the maximal model approximation errors, and on the dimension of the set of representable conditional distributions. We contribute new tools for investigating conditional probability models, which allow us to improve the results that can be derived from existing work on restricted Boltzmann machine probability models.


title: Electronic relaxation rates in metallic ferromagnets

abstract: We show that the single-particle and transport relaxation rates in ferromagnetic metals, which determine the thermal and electrical conductivity, respectively, at asymptotically low temperature do not obey a power law as previously thought, but rather show an exponential temperature dependence. The reason is the splitting of the conduction band that inevitably results from a nonzero magnetization. At higher temperatures there is a sizable temperature window where the transport rate shows a T^2 temperature dependence, in accord with prior results. This window is separated from the asymptotic regime by a temperature scale that is estimated to range from tens of mK to tens of K for typical ferromagnets. We motivate and derive a very general effective theory for metallic magnets that we then use to derive these results. Comparisons with existing experiments are discussed, and predictions for future experiments at low temperatures are made.


title: Measuring the transmission matrix for microwave radiation propagating through random waveguides: fundamentals and applications

abstract: This thesis describes the measurement and analysis of the transmission matrix (TM) for microwave radiation propagating through multichannel random waveguides in the crossover to Anderson localization. Eigenvalues of the transmission matrix and the associated eigenchannels are obtained via a singular value decomposition of the TM. The sum of the transmission eigenvalues yields the transmittance $\it T$, which is the classical analog of the dimensionless conductance $\textsl g$. The dimensionless conductance $\textsl g$ is the electronic conductance in units of the quantum conductance, $G/(e^2/h)$.


title: Quantum Walks of Two Interacting Particles in One Dimension

abstract: We investigate continuous-time quantum walks of two indistinguishable particles (bosons, fermions or hard-core bosons) in one-dimensional lattices with nearest-neighbour interactions. The two interacting particles can undergo independent- and/or co-walking dependent on both quantum statistics and interaction strength. We find that two strongly interacting particles may form a bound state and then co-walk like a single composite particle with statistics-dependent propagation speed. Such an effective single-particle picture of co-walking is analytically derived in the context of degenerate perturbation and the analytical results are well consistent with direct numerical simulation. In addition to implementing universal quantum computation and observing bound states, two-particle quantum walks offer a novel route to detecting quantum statistics. Our theoretical results can be examined in experiments of light propagations in two-dimensional waveguide arrays or spin-impurity dynamics of ultracold atoms in one-dimensional optical lattices.


title: Constant Rank Bimatrix Games are PPAD-hard

abstract: The rank of a bimatrix game (A,B) is defined as rank(A+B). Computing a Nash equilibrium (NE) of a rank-$0$, i.e., zero-sum game is equivalent to linear programming (von Neumann'28, Dantzig'51). In 2005, Kannan and Theobald gave an FPTAS for constant rank games, and asked if there exists a polynomial time algorithm to compute an exact NE. Adsul et al. (2011) answered this question affirmatively for rank-$1$ games, leaving rank-2 and beyond unresolved. In this paper we show that NE computation in games with rank $\ge 3$, is PPAD-hard, settling a decade long open problem. Interestingly, this is the first instance that a problem with an FPTAS turns out to be PPAD-hard. Our reduction bypasses graphical games and game gadgets, and provides a simpler proof of PPAD-hardness for NE computation in bimatrix games. In addition, we get: * An equivalence between 2D-Linear-FIXP and PPAD, improving a result by Etessami and Yannakakis (2007) on equivalence between Linear-FIXP and PPAD. * NE computation in a bimatrix game with convex set of Nash equilibria is as hard as solving a simple stochastic game. * Computing a symmetric NE of a symmetric bimatrix game with rank $\ge 6$ is PPAD-hard. * Computing a (1/poly(n))-approximate fixed-point of a (Linear-FIXP) piecewise-linear function is PPAD-hard. The status of rank-$2$ games remains unresolved.


title: Discrete branching laws for minimal holomorphic representations

abstract: We find the explicit branching laws for the restriction of minimal holomorphic representations to symmetric subgroups in the case where the restriction is discretely decomposable. For holomorphic pairs the minimal holomorphic representation decomposes into a direct sum of lowest weight representations which is made explicit. For non-holomorphic pairs the restriction is shown to be irreducible and identified with a known representation. We further study a conjecture by Kobayashi on the behaviour of associated varieties under restriction and confirm this conjecture in the setting of this paper.


title: Improved Design Method for Nearly Linear-Phase IIR Filters Using Constrained Optimization

abstract: A new optimization method for the design of nearly linear-phase IIR digital filters that satisfy prescribed specifications is proposed. The group-delay deviation is minimized under the constraint that the passband ripple and stopband attenuation are within the prescribed specifications and either a prescribed or an optimized group delay can be achieved. By representing the filter in terms of a cascade of second-order sections, a non-restrictive stability constraint characterized by a set of linear inequality constraints can be incorporated in the optimization algorithm. An additional feature of the method, which is very useful in certain applications, is that it provides the capability of constraining the maximum gain in transition bands to be below a prescribed level. Experimental results show that filters designed using the proposed method have much lower group-delay deviation for the same passband ripple and stopband attenuation when compared with corresponding filters designed with several state-of-the-art competing methods.


title: Numerical solutions of a boundary value problem on the sphere using radial basis functions

abstract: Boundary value problems on the unit sphere arise naturally in geophysics and oceanography when scientists model a physical quantity on large scales. Robust numerical methods play an important role in solving these problems. In this article, we construct numerical solutions to a boundary value problem defined on a spherical sub-domain (with a sufficiently smooth boundary) using radial basis functions (RBF). The error analysis between the exact solution and the approximation is provided. Numerical experiments are presented to confirm theoretical estimates.


title: Adaptive Search Algorithms for Discrete Stochastic Optimization: A Smooth Best-Response Approach

abstract: This paper considers simulation-based optimization of the performance of a regime-switching stochastic system over a finite set of feasible configurations. Inspired by the stochastic fictitious play learning rules in game theory, we propose an adaptive simulation-based search algorithm that uses a smooth best-response sampling strategy and tracks the set of global optima, yet distributes the search so that most of the effort is spent on simulating the system performance at the global optima. The algorithm converges weakly to the set of global optima even when the observation data is correlated (as long as a weak law of large numbers holds). Numerical examples show that the proposed scheme yields a faster convergence for finite sample lengths compared with several existing random search and pure exploration methods in the literature.


title: Robust acoustic cloaking with density-near-zero materials

abstract: Using the phenomenon of extraordinary sound transmission (EST) through artificially constructed density-near-zero (DNZ) materials, we propose an alternative approach to realize a very robust acoustic cloaking, which includes one absorber and one projector connected by an energy channel. The elementary unit cell is made of one single-piece homogeneous copper, but it well maintains both the original planar wavefront and the nearly perfect one-dimensional transmission, in the presence of any inserted objects. Furthermore, the overall cloaked space can be increased, designed, and distributed by arbitrarily assembling additional DNZ cells without any upper limit in the total cloaking volume. We demonstrate acoustic cloaking by different combinations of DNZ cells in free space as well as in acoustic waveguides, which enables any objects inside the cells imperceptible along undistorted sound paths.


title: A Unified Framework of Topological Phases with Symmetry

abstract: In topological phases in $2+1$ dimensions, anyons fall into representations of quantum group symmetries. As proposed in our work (arXiv:1308.4673), physics of a symmetry enriched phase can be extracted by the Mathematics of (hidden) quantum group symmetry breaking of a "parent phase". This offers a unified framework and classification of the symmetry enriched (topological) phases, including symmetry protected trivial phases as well. In this paper, we extend our investigation to the case where the "parent" phases are non-Abelian topological phases. We show explicitly how one can obtain the topological data and symmetry transformations of the symmetry enriched phases from that of the "parent" non-Abelian phase. Two examples are computed: (1) the $\text{Ising}\times\overline{\text{Ising}}$ phase breaks into the $\mathbb{Z}_2$ toric code with $\mathbb{Z}_2$ global symmetry; (2) the $SU(2)_8$ phase breaks into the chiral Fibonacci $\times$ Fibonacci phase with a $\mathbb{Z}_2$ symmetry, a first non-Abelian example of symmetry enriched topological phase beyond the gauge theory construction.


title: Parameter convexity and concavity of generalized trigonometric functions

abstract: We study the convexity properties of the generalized trigonometric functions considered as functions of parameter. We show that $p\to\sin_p(y)$ and $p\to\cos_p(y)$ are log-concave on the appropriate intervals while $p\to\tan_p(y)$ is log-convex. We also prove similar facts about the generalized hyperbolic functions. In particular, our results settle the major part of a conjecture put forward in a recent paper by Baricz, Bhayo and Vuorinen.


title: Block patterns in Stirling permutations

abstract: We introduce and study a new notion of patterns in Stirling and $k$-Stirling permutations, which we call block patterns. We prove a general result which allows us to compute generating functions for the occurrences of various block patterns in terms of generating functions for the occurrences of patterns in permutations. This result yields a number of applications involving, among other things, Wilf equivalence of block patterns and a new interpretation of Bessel polynomials. We also show how to interpret our results for a certain class of labeled trees, which are in bijection with Stirling permutations.


title: Defect-Mediated Relaxation in the Random Tiling Phase of a Binary Mixture: Birth, Death and Mobility of an Atomic Zipper

abstract: This paper describes the mechanism of defect-mediated relaxation in a dodecagonal square-triangle random tiling phase exhibited by a simulated binary mixture of soft discs in 2D. We examine the internal transitions within the elementary mobile defect (christened the 'zipper') that allow it to move, as well as the mechanisms by which the zipper is created and annihilated. The structural relaxation of the random tiling phase is quantified and we show that this relaxation is well described by a model based on the distribution of waiting times for each atom to be visited by the diffusing zipper. This system, representing one of the few instances where a well defined mobile defect is capable of structural relaxation, can provide a valuable test case for general theories of relaxation in complex and disordered materials.


title: The nonequivariant coherent-constructible correspondence and tilting

abstract: The coherent-constructible correspondence is a relationship between coherent sheaves on a toric variety X, and constructible sheaves on a real torus T. This was discovered by Bondal, and explored in the equivariant setting by Fang, Liu, Treumann and Zaslow. In this paper we collect partial results towards a proof of the non equivariant coherent-constructible correspondence. Also, we give applications to the construction of tilting complexes in the derived category of toric DM stacks.


title: On realizations of polynomial algebras with three generators via deformed oscillator algebras

abstract: We present the most general polynomial Lie algebra generated by a second order integral of motion and one of order M, construct the Casimir operator, and show how the Jacobi identity provides the existence of a realization in terms of deformed oscillator algebra. We also present the classical analog of this construction for the most general Polynomial Poisson algebra. Two specific classes of such polynomial algebras are discussed that include the symmetry algebras observed for various 2D superintegrable systems.


title: Magnetic field control of the intraband optical absorption in two-dimensional quantum rings

abstract: Linear and nonlinear optical absorption coefficients of the two-dimensional semiconductor ring in the perpendicular magnetic field $\bf B$ are calculated within independent electron approximation. Characteristic feature of the energy spectrum are crossings of the levels with adjacent nonpositive magnetic quantum numbers as the intensity $B$ changes. It is shown that the absorption coefficient of the associated optical transition is drastically decreased at the fields corresponding to the crossing. Proposed model of the Volcano disc allows to get simple mathematical analytical results, which provide clear physical interpretation. An interplay between positive linear and intensity-dependent negative cubic absorption coefficients is discussed; in particular, critical light intensity at which additional resonances appear in the total absorption dependence on the light frequency, is calculated as a function of the magnetic field and levels' broadening.


title: Orientational Ordering of Passivating Ligands on CdS Nanorods in Solution Generates Strong Rod-Rod Interactions

abstract: We present the first nearly atomistic molecular dynamics study of nanorod-nanorod association in explicit solvent, showing that inter-rod forces can be dominated by microscopic factors absent in common continuum descriptions. Specifically, we find that alkane ligands on faceted CdS nanorods in n-hexane undergo a temperature-dependent order-disorder transition akin to that of self-assembled monolayers on macroscopic substrates. This collective ligand alignment organizes nearby solvent molecules, strongly influencing the statistics of rod-rod separation. The strong temperature-dependence of this mechanism could be exploited in the laboratory to manipulate and optimize the assembly of ordered structures.


title: Metric tree-like structures in real-life networks: an empirical study

abstract: Based on solid theoretical foundations, we present strong evidences that a number of real-life networks, taken from different domains like Internet measurements, biological data, web graphs, social and collaboration networks, exhibit tree-like structures from a metric point of view. We investigate few graph parameters, namely, the tree-distortion and the tree-stretch, the tree-length and the tree-breadth, the Gromov's hyperbolicity, the cluster-diameter and the cluster-radius in a layering partition of a graph, which capture and quantify this phenomenon of being metrically close to a tree. By bringing all those parameters together, we not only provide efficient means for detecting such metric tree-like structures in large-scale networks but also show how such structures can be used, for example, to efficiently and compactly encode approximate distance and almost shortest path information and to fast and accurately estimate diameters and radii of those networks. Estimating the diameter and the radius of a graph or distances between its arbitrary vertices are fundamental primitives in many data and graph mining algorithms.


title: Regularization Parameter Estimation for Underdetermined problems by the $\chi^2$ principle with application to $2D$ focusing gravity inversion

abstract: The $\chi^2$-principle generalizes the Morozov discrepancy principle (MDP) to the augmented residual of the Tikhonov regularized least squares problem. Weighting of the data fidelity by a known Gaussian noise distribution on the measured data, when the regularization term is weighted by unknown inverse covariance information on the model parameters, the minimum of the Tikhonov functional is a random variable following a $\chi^2$-distribution with $m+p-n$ degrees of freedom, model matrix $G:$ $m \times n$ and regularizer $L:$ $p\times n$. It is proved that the result holds also for $m<n$ when $m+p\ge n$. A Newton root-finding algorithm is used to find the regularization parameter $\alpha$ which yields the optimal inverse covariance weighting in the case of a white noise assumption on the mapped model data. It is implemented for small-scale problems using the generalized singular value decomposition. Numerical results verify the algorithm for the case of regularizers approximating zero to second order derivative approximations, contrasted with the methods of generalized cross validation and unbiased predictive risk estimation. The inversion of underdetermined $2D$ focusing gravity data produces models with non-smooth properties, for which typical implementations in this field use the iterative minimum support (MS) stabilizer and both regularizer and regularizing parameter are updated each iteration. For a simulated data set with noise, the regularization parameter estimation methods for underdetermined data sets are used in this iterative framework, also contrasted with the L-curve and MDP. Experiments demonstrate efficiency and robustness of the $\chi^2$-principle, moreover the L-curve and MDP are generally outperformed. Furthermore, the MS is of general use for the $\chi^2$-principle when implemented without the knowledge of a mean value of the model.


title: Quantum Phases and Collective Excitations in Bose-Hubbard Models with Staggered Magnetic Flux

abstract: We study the quantum phases of a Bose-Hubbard model with staggered magnetic flux in two dimensions, as has been realized recently [Aidelsburger {\it et al.}, PRL, {\bf 107}, 255301 (2011)]. Within mean field theory, we show how the structure of the condensates evolves from weak to strong coupling limit, exhibiting a tricritical point at the Mott-superfluid transition. Non-trivial topological structures (Dirac points) in the quasi-particle (hole) excitations in the Mott state are found within random phase approximation and we discuss how interaction modifies their structures. Excitation gap in the Mott state closes at different ${\bf k}$ points when approaching the superfluid states, which is consistent with the findings of mean field theory.


title: An Electrostatics Problem on the Sphere Arising from a Nearby Point Charge

abstract: For a positively charged insulated d-dimensional sphere we investigate how the distribution of this charge is affected by proximity to a nearby positive or negative point charge when the system is governed by a Riesz s-potential 1/r^s, s>0, where r denotes Euclidean distance between point charges. Of particular interest are those distances from the point charge to the sphere for which the equilibrium charge distribution is no longer supported on the whole of the sphere (i.e. spherical caps of negative charge appear). Arising from this problem attributed to A. A. Gonchar are sequences of polynomials of a complex variable that have some fascinating properties regarding their zeros.