title: On multi-scale percolation behaviour of the effective conductivity for the lattice model

abstract: Macroscopic properties of heterogeneous media are frequently modelled by regular lattice models, which are based on a relatively small basic cluster of lattice sites. Here, we extend one of such models to any cluster's size kxk. We also explore its modified form. The focus is on the percolation behaviour of the effective conductivity of random two- and three-phase systems. We consider only the influence of geometrical features of local configurations at different length scales k. At scales accessible numerically, we find that an increase in the size of the basic cluster leads to characteristic displacements of the percolation threshold. We argue that the behaviour is typical of materials, whose conductivity is dominated by a few linear, percolation-like, conducting paths. Such a system can be effectively treated as one-dimensional medium. We also develop a simplified model that permits of an analysis at any scale. It is worth mentioning that the latter approach keeps the same thresholds predicted by the former one. We also briefly discuss a three-phase system, where the double-thresholds paths appear on model surfaces.


title: Atomic decompositions for Hardy spaces related to Schr\"odinger operators

abstract: Let L_U = -Delta+U be a Schr\"odinger operator on R^d, where U\in L^1_{loc}(R^d) is a non-negative potential and d\geq 3. The Hardy space H^1(L_U) is defined in terms of the maximal function for the semigroup K_{t,U} = exp(-t L_U), namely H^1(L_U) = {f\in L^1(R^d): \|f\|_{H^1(L_U)}:= \|sup_{t>0} |K_{t,U} f| \|_{L^1(R^d)} < \infty. Assume that U=V+W, where V\geq 0 satisfies the global Kato condition sup_{x\in R^d} \int_{R^d} V(y)|x-y|^{2-d} < \infty. We prove that, under certain assumptions on W\geq 0, the space H^1(L_U) admits an atomic decomposition of local type. An atom a for H^1(L_U) is either of the form a(x)=|Q|^{-1}\chi_Q(x), where Q are special cubes determined by W, or a satisfies the cancellation condition \int a(x)w(x) dx = 0, where w is an (-Delta+V)-harmonic function given by w(x) = lim_{t\to \infty} K_{t,V} 1(x). Furthermore, we show that, in some cases, the cancellation condition \int_{R^d} a(x)w(x) dx = 0 can be replaced by the classical one \int_{R^d} a(x) dx = 0. However, we construct another example, such that the atomic spaces with these two cancellation conditions are not equivalent as Banach spaces.


title: Similarity of closed polygonal curves in Frechet metric

abstract: The article analyzes similarity of closed polygonal curves in Frechet metric, which is stronger than the well-known Hausdorff metric and therefore is more appropriate in some applications. An algorithm that determines whether the Frechet distance between two closed polygonal curves with m and n vertices is less than a given number is described. The described algorithm takes O(mn) time whereas the previously known algorithms take O(mn log(mn)) time.


title: Lexical Normalisation of Twitter Data

abstract: Twitter with over 500 million users globally, generates over 100,000 tweets per minute . The 140 character limit per tweet, perhaps unintentionally, encourages users to use shorthand notations and to strip spellings to their bare minimum "syllables" or elisions e.g. "srsly". The analysis of twitter messages which typically contain misspellings, elisions, and grammatical errors, poses a challenge to established Natural Language Processing (NLP) tools which are generally designed with the assumption that the data conforms to the basic grammatical structure commonly used in English language. In order to make sense of Twitter messages it is necessary to first transform them into a canonical form, consistent with the dictionary or grammar. This process, performed at the level of individual tokens ("words"), is called lexical normalisation. This paper investigates various techniques for lexical normalisation of Twitter data and presents the findings as the techniques are applied to process raw data from Twitter.


title: Non-Archimedean Whittaker functions as characters: a probabilistic approach to the Shintani-Casselman-Shalika formula

abstract: For a reductive group $G$ over a non-Archimedean local field (e.g $GL_n( \mathbb{Q}_p )$ ), Jacquet's Whittaker function is essentially proportional to a character of an irreducible representation of the Langlands dual group $G^\vee( \mathbb{C} )$ ( a Schur function if $G = GL_n( \mathbb{Q}_p )$). We propose a probabilistic approach to this claim, known as the Shintani-Casselman-Shalika formula, when the group $G$ has at least one minuscule cocharacter in the coweight lattice. Our presentation goes along the following lines. Thanks to a minuscule random walk $W^{(z)}$ on the coweight lattice and a related random walk on the Borel subgroup, we establish a Poisson kernel formula for the non-Archimedean Whittaker function. The expression and its ingredients are similar to the one previously obtained by the author in the Archimedean case. A simple manipulation reduces the problem to evaluating the probability of $W^{(z)}$ never exiting the Weyl chamber. Then, an implementation of the reflection principle forces the appearance of the Weyl character formula and therefore retrieves characters of $G^\vee\left( \mathbb{C} \right)$. The construction of the random walk on the Borel subgroup requires some care. It is extracted from a spherical random walk whose increments have a distribution that can be understood as elements from the spherical Hecke algebra.


title: Hodge integrals and tau-symmetric integrable hierarchies of Hamiltonian evolutionary PDEs

abstract: For an arbitrary semisimple Frobenius manifold we construct Hodge integrable hierarchy of Hamiltonian partial differential equations. In the particular case of quantum cohomology the tau-function of a solution to the hierarchy generates the intersection numbers of the Gromov--Witten classes and their descendents along with the characteristic classes of Hodge bundles on the moduli spaces of stable maps. For the one-dimensional Frobenius manifold the Hodge hierarchy is a deformation of the Korteweg--de Vries hierarchy depending on an infinite number of parameters. Conjecturally this hierarchy is a universal object in the class of scalar Hamiltonian integrable hierarchies possessing tau-functions.


title: The Role of Emotions in Propagating Brands in Social Networks

abstract: A key aspect of word of mouth marketing are emotions. Emotions in texts help propagating messages in conventional advertising. In word of mouth scenarios, emotions help to engage consumers and incite to propagate the message further. While the function of emotions in offline marketing in general and word of mouth marketing in particular is rather well understood, online marketing can only offer a limited view on the function of emotions. In this contribution we seek to close this gap. We therefore investigate how emotions function in social media. To do so, we collected more than 30,000 brand marketing messages from the Google+ social networking site. Using state of the art computational linguistics classifiers, we compute the sentiment of these messages. Starting out with Poisson regression-based baseline models, we seek to replicate earlier findings using this large data set. We extend upon earlier research by computing multi-level mixed effects models that compare the function of emotions across different industries. We find that while the well known notion of activating emotions propagating messages holds in general for our data as well. But there are significant differences between the observed industries.


title: Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements

abstract: We propose an effective and flexible way to assemble finite element stiffness and mass matrices in MATLAB. We apply this for problems discretized by edge finite elements. Typical edge finite elements are Raviart-Thomas elements used in discretizations of H(div) spaces and Nedelec elements in discretizations of H(curl) spaces. We explain vectorization ideas and comment on a freely available MATLAB code which is fast and scalable with respect to time.


title: Measurement of the CP-violating phase $\phi_s$ in $\bar{B}^{0}_{s}\to D_{s}^{+}D_{s}^{-}$ decays

abstract: We present a measurement of the $CP$-violating weak mixing phase $\phi_s$ using the decay $\bar{B}^{0}_{s}\to D_{s}^{+}D_{s}^{-}$ in a data sample corresponding to $3.0$ fb$^{-1}$ of integrated luminosity collected with the LHCb detector in $pp$ collisions at centre-of-mass energies of 7 and 8 TeV. An analysis of the time evolution of the system, which does not constrain $|\lambda|=1$ to allow for the presence of $CP$ violation in decay, yields $\phi_s = 0.02 \pm 0.17$ (stat) $\pm 0.02$ (syst) rad, $|\lambda| = 0.91^{+0.18}_{-0.15}$ (stat) $\pm0.02$ (syst). This result is consistent with the Standard Model expectation.


title: Linear spectro-polarimetry: a new diagnostic tool for the classification and characterisation of asteroids

abstract: We explore the use of spectro-polarimetry as a remote sensing tool for asteroids in addition to traditional reflectance measurements. In particular we are interested in possible relationships between the wavelength-dependent variation of linear polarization and the properties of the surfaces, including albedo and composition. We have obtained optical spectro-polarimetric measurements of a dozen asteroids of different albedo and taxonomic classes and of two small regions at the limb of the Moon. We found that objects with marginally different relative reflectance spectra (in the optical) may have totally different polarization spectra. This suggests that spectro-polarimetry may be used to refine the classification of asteroids. We also found that in some cases the Umov law may be violated, that is, in contrast to what is expected from basic physical considerations, the fraction of linear polarization and the reflectance may be positively correlated. In agreement with a few previous studies based on multi-colour broadband polarimetry, we found that the variation of linear polarization with wavelength and with phase-angle is correlated with the albedo and taxonomic class of the objects. Finally, we have serendipitously discovered that spinel-rich asteroid (599) Luisa, located very close to the Watsonia family, is a member of the rare class of Barbarian asteroids. We suggest that future modelling attempts of the surface structure of asteroids should be aimed at explaining both reflectance and polarization spectra.


title: Approximability of Guarding Weak Visibility Polygons

abstract: The art gallery problem enquires about the least number of guards that are sufficient to ensure that an art gallery, represented by a polygon $P$, is fully guarded. In 1998, the problems of finding the minimum number of point guards, vertex guards, and edge guards required to guard $P$ were shown to be APX-hard by Eidenbenz, Widmayer and Stamm. In 1987, Ghosh presented approximation algorithms for vertex guards and edge guards that achieved a ratio of $\mathcal{O}(\log n)$, which was improved upto $\mathcal{O}(\log\log OPT)$ by King and Kirkpatrick in 2011. It has been conjectured that constant-factor approximation algorithms exist for these problems. We settle the conjecture for the special class of polygons that are weakly visible from an edge and contain no holes by presenting a 6-approximation algorithm for finding the minimum number of vertex guards that runs in $\mathcal{O}(n^2)$ time. On the other hand, for weak visibility polygons with holes, we present a reduction from the Set Cover problem to show that there cannot exist a polynomial time algorithm for the vertex guard problem with an approximation ratio better than $((1 - \epsilon)/12)\ln n$ for any $\epsilon>0$, unless NP=P. We also show that, for the special class of polygons without holes that are orthogonal as well as weakly visible from an edge, the approximation ratio can be improved to 3. Finally, we consider the Point Guard problem and show that it is NP-hard in the case of polygons weakly visible from an edge.


title: Optimal two-qubit tomography based on local and global measurements: Maximal robustness against errors as described by condition numbers

abstract: We present an error analysis of various tomographic protocols based on the linear inversion for the reconstruction of an unknown two-qubit state. We solve the problem of finding a tomographic protocol which is the most robust against errors in terms of the lowest value (i.e., equal to 1) of a condition number, as required by the Gastinel-Kahan theorem. In contrast, standard tomographic protocols, including those based on mutually unbiased bases, are nonoptimal for determining all 16 elements of an unknown two-qubit density matrix. Our method is based on the measurements of the 16 generalized Pauli operators, where twelve of them can be locally measured, and the other four require nonlocal Bell measurements. Our method corresponds to selectively measuring, one by one, all of the real and imaginary elements of an unknown two-qubit density matrix. We describe two experimentally feasible setups of this protocol for the optimal reconstruction of two photons in an unknown polarization state using conventional detectors and linear-optical elements. Moreover, we define the operators for the optimal reconstruction of the states of multiqubit or multilevel (qudit) systems.


title: Can the Abraham Light Momentum and Energy in a Medium Constitute a Lorentz Four-Vector?

abstract: By analyzing the Einstein-box thought experiment with the principle of relativity, it is shown that Abraham's light momentum and energy in a medium cannot constitute a Lorentz four-vector, and they consequentially break global momentum and energy conservation laws. In contrast, Minkowski's momentum and energy always constitute a Lorentz four-vector no matter whether in a medium or in vacuum, and the Minkowski's momentum is the unique correct light momentum. A momentum-associated photon mass in a medium is exposed, which explains why only the Abraham's momentum is derived in the traditional "center-of-mass-energy" approach. The EM boundary-condition matching approach, combined with Einstein light-quantum hypothesis, is proposed to analyze this thought experiment, and it is found for the first time that only from Maxwell equations without resort to the relativity, the correctness of light momentum definitions cannot be identified. Optical pulling effect is studied as well.


title: Decomposition of Differential Games

abstract: This paper provides a decomposition technique for the purpose of simplifying the solution of certain zero-sum differential games. The games considered terminate when the state reaches a target, which can be expressed as the union of a collection of target subsets; the decomposition consists of replacing the original target by each of the target subsets. The value of the original game is then obtained as the lower envelope of the values of the collection of games resulting from the decomposition, which can be much easier to solve than the original game. Criteria are given for the validity of the decomposition. The paper includes examples, illustrating the application of the technique to pursuit/evasion games, where the decomposition arises from considering the interaction of individual pursuer/evader pairs.


title: The effect of anchoring on nematic flow in channels

abstract: Understanding the flow of liquid crystals in microfluidic environments plays an important role in many fields, including device design and microbiology. We perform hybrid lattice-Boltzmann simulations of a nematic liquid crystal flowing under an applied pressure gradient in two-dimensional channels with various anchoring boundary conditions at the substrate walls. We investigate the relation between flow rate and pressure gradient and the corresponding profile of the nematic director, and find significant departures from the linear Poiseuille relation. We also identify a morphological transition in the director profile and explain this in terms of an instability in the dynamical equations. We examine the qualitative and quantitative effects of changing the type and strength of the anchoring. Understanding such effects may provide a useful means of quantifying the anchoring of a substrate by measuring its flow properties.


title: Laboratory Test Bench for Research Network and Cloud Computing

abstract: At present moment, there is a great interest in development of information systems operating in cloud infrastructures. Generally, many of tasks remain unresolved such as tasks of optimization of large databases in a hybrid cloud infrastructure, quality of service (QoS) at different levels of cloud services, dynamic control of distribution of cloud resources in application systems and many others. Research and development of new solutions can be limited in case of using emulators or international commercial cloud services, due to the closed architecture and limited opportunities for experimentation. Article provides answers to questions on the establishment of a pilot cloud practically "at home" with the ability to adjust the width of the emulation channel and delays in data transmission. It also describes architecture and configuration of the experimental setup. The proposed modular structure can be expanded by available computing power.


title: DISA at ImageCLEF 2014 Revised: Search-based Image Annotation with DeCAF Features

abstract: This paper constitutes an extension to the report on DISA-MU team participation in the ImageCLEF 2014 Scalable Concept Image Annotation Task as published in [3]. Specifically, we introduce a new similarity search component that was implemented into the system, report on the results achieved by utilizing this component, and analyze the influence of different similarity search parameters on the annotation quality.


title: Managing HIV/AIDS in Malawi

abstract: The epidemic of HIV in Malawi started early and at its peak 15% of all adults were infected with HIV. Malawi is a low-income country and the cost of putting all HIV-positive people in Malawi onto ART, expressed as a percentage of the gross domestic product, is the highest in the world. Nevertheless, Malawi has made great progress and the greatly reduced cost of potent anti-retroviral therapy (ART) makes it possible to contemplate ending the epidemic of HIV/AIDS. Here we consider what would have happened without ART, the No ART counterfactual, the impact if the current level of roll-out of ART is maintained, the Current Programme, and the likely impact if treatment is made available to everyone who is eligible under the 2013 guidelines of the World Health Organization reaching full coverage by 2020, the Expanded Programme. The Current Programme has substantially reduced the epidemic of HIV and the number of people dying of AIDS. The Expanded Programme has the potential to avert more infections, save more lives and end the epidemic. The annual cost of managing HIV will increase from about US$132 million in 2014 to about US$155 million in 2020 but will fall after that. If the Expanded Programme is implemented several key areas must be addressed. Testing services will need to be expanded and supported by mass testing campaigns, so as to diagnose people with HIV and enrol them in treatment and care as early as possible. A regular and uninterrupted supply of drugs will have to be assured. The quantity and quality of existing health staff will need to be strengthened. Community health workers will need to be mobilized and trained to encourage people to be tested and accept treatment, to monitor progress and to support people on treatment; this in turn will help to reduce stigma and discrimination, loss to follow up of people diagnosed with HIV, and improve adherence for those on treatment.


title: Resolute: An Assurance Case Language for Architecture Models

abstract: Arguments about the safety, security, and correctness of a complex system are often made in the form of an assurance case. An assurance case is a structured argument, often represented with a graphical interface, that presents and supports claims about a system's behavior. The argument may combine different kinds of evidence to justify its top level claim. While assurance cases deliver some level of guarantee of a system's correctness, they lack the rigor that proofs from formal methods typically provide. Furthermore, changes in the structure of a model during development may result in inconsistencies between a design and its assurance case. Our solution is a framework for automatically generating assurance cases based on 1) a system model specified in an architectural design language, 2) a set of logical rules expressed in a domain specific language that we have developed, and 3) the results of other formal analyses that have been run on the model. We argue that the rigor of these automatically generated assurance cases exceeds those of traditional assurance case arguments because of their more formal logical foundation and direct connection to the architectural model.


title: Contact potential instability in the path-integral description of itinerant ferromagnetism

abstract: It has long been predicted that a two-component non-localized Fermi gas will exhibit spontaneous polarization for sufficiently strong repulsive interactions, a phenomenon which is called itinerant ferromagnetism. Recent experiments with ultracold atomic gases have reached the interaction strength for which theoretical models have predicted the occurrence of the normal-to-itinerant-ferromagnetic phase transition, but so far this transition has not been observed. The instability of the repulsive branch of the Feshbach resonance prevents the formation of the itinerant ferromagnetic state, but it is not clear whether this is the only instability impeding its experimental realization. In this article, we use the path-integral formalism with density fields in the Hubbard-Stratonovich transformation to study the stability of a homogeneous two-component Fermi gas with contact interactions. Within the saddle-point approximation we show that none of the extrema of the action are minima, meaning all extrema are unstable to small density fluctuations. This implies a more general mechanical instability of the polarized (itinerant ferromagnetic) and normal states of the system in the path-integral formalism. We find that it is important to consider the stability of the system when studying itinerant ferromagnetism. Since (mechanical) stability may be influenced by the details of the interaction potential, we suggest the use of a more realistic potential than the contact potential in future theoretical descriptions.


title: The geometry of dimer models

abstract: This is an expanded version of a three-hour minicourse given at the winterschool Winterbraids IV held in Dijon in February 2014. The aim of these lectures was to present some aspects of the dimer model to a geometrically minded audience. We spoke neither of braids nor of knots, but tried to show how several geometrical tools that we know and love (e.g. (co)homology, spin structures, real algebraic curves) can be applied to very natural problems in combinatorics and statistical physics. These lecture notes do not contain any new results, but give a (relatively original) account of the works of Kasteleyn, Cimasoni-Reshetikhin and Kenyon-Okounkov-Sheffield.


title: A characterization of $K_{2,4}$-minor-free graphs

abstract: We provide a complete structural characterization of $K_{2,4}$-minor-free graphs. The $3$-connected $K_{2,4}$-minor-free graphs consist of nine small graphs on at most eight vertices, together with a family of planar graphs that contains $K_4$ and, for each $n \ge 5$, $2n-8$ nonisomorphic graphs of order $n$. To describe the $2$-connected $K_{2,4}$-minor-free graphs we use $xy$-outerplanar graphs, graphs embeddable in the plane with a Hamilton $xy$-path so that all other edges lie on one side of this path. We show that, subject to an appropriate connectivity condition, $xy$-outerplanar graphs are precisely the graphs that have no rooted $K_{2,2}$-minor where $x$ and $y$ correspond to the two vertices on one side of the bipartition of $K_{2,2}$. Each $2$-connected $K_{2,4}$-minor-free graph is then (i) outerplanar, (ii) the union of three $xy$-outerplanar graphs and possibly the edge $xy$, or (iii) obtained from a $3$-connected $K_{2,4}$-minor-free graph by replacing each edge $x_iy_i$ in a set $\{x_1 y_1, x_2 y_2, \ldots, x_k y_k\}$ satisfying a certain condition by an $x_i y_i$-outerplanar graph.


title: Global Existence and Regularity Results for Strongly Coupled Nonregular Parabolic Systems via Iterative Methods

abstract: The global existence of classical solutions to strongly coupled parabolic systems is shown to be equivalent to the availability of an iterative scheme producing a sequence of solutions with uniform continuity in the BMO norms. Amann's results on global existence of classical solutions still hold under much weaker condition that their BMO norms do not blow up in finite time. The proof makes use of some new global and local weighted Gagliardo-Nirenberg inequalities involving BMO norms.


title: Scattering of a Single Plasmon by Three Non-equally Spaced Quantum Dots System Coupled to One-Dimensional Waveguide

abstract: Scattering properties of a single plasm on interacting with three non-equally spaced quantum dots coupled to one-dimensional surface plasmonic waveguide is investigated theoretically via the real-space approach. It is demonstrated that the transmission and reflection of a single plasmon can be switched on or off by controlling the detuning and changing the interparticle distances between the quantum dots. By controlling the transition frequencies and interparticle distances of QDs, one can construct a half-transmitting mirror with three QDs system. We also showed that controlling the transition frequencies and interparticle distances of QDs results in the complete transmission peak near the zero detuning.


title: Exterior powers of F-zips

abstract: An F-zip over a field of positive characteristic is a vector space together with two filtrations whose subquotients are related in a certain way. We will define the category of F-zips and some basic constructions in it, especially exterior powers. If the ground field is algebraically closed, one can give a classification of F-Zips in terms of combinatorics. However, the way constructions and concepts in the category of F-zips manifest themselves in terms of the classifying invariant, is yet to be fully understood. The theory of F-crystals suggests that another invariant might be useful in trying to improve the understanding of F-zips. Given an F-zip, we calculate for every $1$-dimensional F-zip (of which there is essentially one for every integer $d$) $\mathbf{1}(d)$ and every $r\in\mathbb{Z}$ the dimension of the space of F-zip morphisms from $\mathbf{1}(d)$ into the $r$-th exterior power of the given F-zip. To make sense of this however, we will have to canonically decompose these spaces of morphisms each into two subspaces that are finite-dimensional over the prime field and its prime field respectively. One result will then be a way to calculate these numbers for a given isomorphism type. Our main result however, is a negative one: The invariant does not classify F-zips.


title: Switching of a Single Photon by Two {\Lambda}-type Three-Level Quantum Dots Embedded in Cavities Coupling to One-Dimensional Waveguide

abstract: Switching of a single photon interacting with two {\Lambda}-type three-level quantum dots embedded in cavities coupled to one-dimensional waveguide is investigated theoretically via the real-space approach. We demonstrated that switching of a single photon can be achieved by tuning the classic driving field on or off, and by controlling the QD-cavity coupling strength, Rabi frequency and the cavity-waveguide coupling rate. The transmission properties of a single photon by such a nanosystem discussed here could find the applications in the design of next-generation quantum devices and quantum information.


title: Automatic Error Localization for Software using Deductive Verification

abstract: Even competent programmers make mistakes. Automatic verification can detect errors, but leaves the frustrating task of finding the erroneous line of code to the user. This paper presents an automatic approach for identifying potential error locations in software. It is based on a deductive verification engine, which detects errors in functions annotated with pre- and post-conditions. Using an automatic theorem prover, our approach finds expressions in the code that can be modified such that the program satisfies its specification. Scalability is achieved by analyzing each function in isolation. We have implemented our approach in the widely used Frama-C framework and present first experimental results. This is an extended version of [8], featuring an additional appendix.


title: Suppressed epidemics in multi-relational networks

abstract: A two-state epidemic model in networks with links mimicking two kinds of relationships between connected nodes is introduced. Links of weights w1 and w0 occur with probabilities p and 1-p, respectively. The fraction of infected nodes rho(p) shows a non-monotonic behavior, with rho drops with p for small p and increases for large p. For small to moderate w1/w0 ratios, rho(p) exhibits a minimum that signifies an optimal suppression. For large w1/w0 ratios, the suppression leads to an absorbing phase consisting only of healthy nodes within a range p_L =< p =< p_R, and an active phase with mixed infected and healthy nodes for p < p_L and p>p_R. A mean field theory that ignores spatial correlation is shown to give qualitative agreement and capture all the key features. A physical picture that emphasizes the intricate interplay between infections via w0 links and within clusters formed by nodes carrying the w1 links is presented. The absorbing state at large w1/w0 ratios results when the clusters are big enough to disrupt the spread via w0 links and yet small enough to avoid an epidemic within the clusters. A theory that uses the possible local environments of a node as variables is formulated. The theory gives results in good agreement with simulation results, thereby showing the necessity of including longer spatial correlations.


title: Ultimate boundedness of droop controlled Microgrids with secondary loops

abstract: In this paper we study theoretical properties of inverter-based microgrids controlled via primary and secondary loops. Stability of these microgrids has been the subject of a number of recent studies. Conventional approaches based on standard hierarchical control rely on time-scale separation between primary and secondary control loops to show local stability of equilibria. In this paper we show that (i) frequency regulation can be ensured without assuming time-scale separation and, (ii) ultimate boundedness of the trajectories starting inside a region of the state space can be guaranteed under a condition on the inverters power injection errors. The trajectory ultimate bound can be computed by simple iterations of a nonlinear mapping and provides a certificate of the overall performance of the controlled microgrid.


title: Distribution of zeros of polynomials with positive coefficients

abstract: We describe the limit zero distributions of sequences of polynomials with positive coefficients.


title: Exact WKB analysis and cluster algebras II: Simple poles, orbifold points, and generalized cluster algebras

abstract: This is a continuation of developing mutation theory in exact WKB analysis using the framework of cluster algebras. Here we study the Schrodinger equation on a compact Riemann surface with turning points of simple-pole type. We show that the orbifold triangulations by Felikson, Shapiro, and Tumarkin provide a natural framework of describing the mutation of Stokes graphs, where simple poles correspond to orbifold points. We then show that under the mutation of Stokes graphs around simple poles the Voros symbols mutate as the variables of generalized cluster algebras introduced by Chekhov and Shapiro.


title: The strong representation for the nonparametric estimation of length-biased and right-censored data

abstract: In this article, we consider a useful product-limit estimator of distribution function proposed by Huang & Qin(2011) when the observations are subject to length-biased and right-censored data. The estimator retains the simple closed-form expression of the truncation product-limit estimator with some good properties. An almost sure representation for the estimator is obtained which can be used to derive many properties of functional statistics based on this product-limit estimator. The rate for the remainder in the representation is of order $O(n^{-3/4}(\log(n)^{3/4}) $ a.s.


title: Ytterbium-driven strong enhancement of electron-phonon coupling in graphene

abstract: We present high-resolution angle-resolved photoemission spectroscopy study in conjunction with first principles calculations to investigate how the interaction of electrons with phonons in graphene is modified by the presence of Yb. We find that the transferred charges from Yb to the graphene layer hybridize with the graphene $\pi$ bands, leading to a strong enhancement of the electron-phonon interaction. Specifically, the electron-phonon coupling constant is increased by as much as a factor of 10 upon the introduction of Yb with respect to as grown graphene ($\leq$0.05). The observed coupling constant constitutes the highest value ever measured for graphene and suggests that the hybridization between graphene and the adatoms might be a critical parameter in realizing superconducting graphene.


title: Gorenstein formats, canonical and Calabi-Yau threefolds

abstract: We extend the known classification of threefolds of general type that are complete intersections to various classes of non-complete intersections, and find other classes of polarised varieties, including Calabi-Yau threefolds with canonical singularities, that are not complete intersections. Our methods apply more generally to construct orbifolds described by equations in given Gorenstein formats.


title: Shape co-existence and parity doublet in Zr isotopes

abstract: We studied the ground and excited states properties for Zr isotopes starting^M from proton to neutron drip-lines using the relativistic and non-relativistic mean field formalisms with BCS and Bogolyubov pairing. The celebrity NL3 and SLy4 parameter sets are used in the calculations. We find spherical ground and low-lying^M superdeformed excited states in most of the isotopes. Several couples of^M $\Omega^{\pi}=1/2^{\pm}$ parity doublets configurations are found,^M while analyzing the single-particle energy levels of the superdeformed ^M configurations.


title: Electron Spin Decoherence in Silicon Carbide Nuclear Spin Bath

abstract: In this paper, we study the electron spin decoherence of single defects in silicon carbide (SiC) nuclear spin bath. We find that, although the natural abundance of $^{29}\rm{Si}$ ($p_{\rm{Si}}=4.7\%$) is about 4 times larger than that of $^{13}{\rm C}$ ($p_{\rm{C}}=1.1\%$), the electron spin coherence time of defect centers in SiC nuclear spin bath in strong magnetic field ($B>300~\rm{Gauss}$) is longer than that of nitrogen-vacancy (NV) centers in $^{13}{\rm C}$ nuclear spin bath in diamond. The reason for this counter-intuitive result is the suppression of heteronuclear-spin flip-flop process in finite magnetic field. Our results show that electron spin of defect centers in SiC are excellent candidates for solid state spin qubit in quantum information processing.


title: The thermodynamics of creating correlations: Limitations and optimal protocols

abstract: We establish a rigorous connection between fundamental resource theories at the quantum scale. Correlations and entanglement constitute indispensable resources for numerous quantum information tasks. However, their establishment comes at the cost of energy, the resource of thermodynamics, and is limited by the initial entropy. Here, the optimal conversion of energy into correlations is investigated. Assuming the presence of a thermal bath, we establish general bounds for arbitrary systems and construct a protocol saturating them. The amount of correlations, quantified by the mutual information, can increase at most linearly with the available energy, and we determine where the linear regime breaks down. We further consider the generation of genuine quantum correlations, focusing on the fundamental constituents of our universe: fermions and bosons. For fermionic modes, we find the optimal entangling protocol. For bosonic modes, we show that while Gaussian operations can be outperformed in creating entanglement, their performance is optimal for high energies.


title: Fermi-normal, optical, and wave-synchronous coordinates for spacetime with a plane gravitational wave

abstract: Fermi normal coordinates provide a standardized way to describe the effects of gravitation from the point of view of an inertial observer. These coordinates have always been introduced via perturbation expansions and were usually limited to distances much less than the characteristic length scale set by the curvature of spacetime. For a plane gravitational wave this scale is given by its wavelength which defines the domain of validity for these coordinates known as the long-wavelength regime. The symmetry of this spacetime, however, allows us to extend Fermi normal coordinates far beyond the long-wavelength regime. Here we present an explicit construction for this long-range Fermi normal coordinate system based on the unique solution of the boundary-value problem for spacelike geodesics. The resulting formulae amount to summation of the infinite series for Fermi normal coordinates previously obtained with perturbation expansions. We also consider two closely related normal coordinate systems: optical coordinates which are built from null geodesics and wave-synchronous coordinates which are built from spacelike geodesics locked in phase with the propagating gravitational wave. The wave-synchronous coordinates yield the exact solution of Peres and Ehlers-Kundt which is globally defined. In this case, the limitation of the long-wavelength regime is completely overcome, and the system of wave-synchronous coordinates becomes valid for arbitrarily large distances. Comparison of the different coordinate systems is done by considering the motion of an inertial test mass in the field of a plane gravitational wave.


title: Functoriality and duality in Morse-Conley-Floer homology

abstract: In~\cite{rotvandervorst} a homology theory --Morse-Conley-Floer homology-- for isolated invariant sets of arbitrary flows on finite dimensional manifolds is developed. In this paper we investigate functoriality and duality of this homology theory. As a preliminary we investigate functoriality in Morse homology. Functoriality for Morse homology of closed manifolds is known~\cite{abbondandoloschwarz, aizenbudzapolski,audindamian, kronheimermrowka, schwarz}, but the proofs use isomorphisms to other homology theories. We give direct proofs by analyzing appropriate moduli spaces. The notions of isolating map and flow map allows the results to generalize to local Morse homology and Morse-Conley-Floer homology. We prove Poincar\'e type duality statements for local Morse homology and Morse-Conley-Floer homology.


title: Constraining dark sector perturbations I: cosmic shear and CMB lensing

abstract: We present current and future constraints on equations of state for dark sector perturbations. The equations of state considered are those corresponding to a generalized scalar field model and time-diffeomorphism invariant $L(g)$ theories that are equivalent to models of a relativistic elastic medium and also Lorentz violating massive gravity. We develop a theoretical understanding of the observable impact of these models. In order to constrain these models we use CMB temperature data from Planck, BAO measurements, CMB lensing data from Planck and the South Pole Telescope, and weak galaxy lensing data from CFHTLenS. We find non-trivial exclusions on the range of parameters, although the data remains compatible with $w=-1$. We gauge how future experiments will help to constrain the parameters. This is done via a likelihood analysis for CMB experiments such as CoRE and PRISM, and tomographic galaxy weak lensing surveys, focussing in on the potential discriminatory power of Euclid on mildly non-linear scales.


title: Accretion of Dark Energy onto Higher Dimensional Charged BTZ Black Hole

abstract: In this work, we have studied the accretion of $(n+2)$-dimensional charged BTZ black hole (BH). The critical point and square speed of sound have been obtained. The mass of the BTZ BH has been calculated and we have observed that the mass of the BTZ BH is related with square root of the energy density of dark energy which accretes onto BH in our accelerating FRW universe. We have assumed modified Chaplygin gas (MCG) as a candidate of dark energy which accretes onto BH and we have found the expression of BTZ BH mass. Since in our solution of MCG, this model generates only quintessence dark energy (not phantom) and so BTZ BH mass increases during the whole evolution of the accelerating universe. Next we have assumed 5 kinds of parametrizations of well known dark energy models. These models generate both quintessence and phantom scenarios i.e., phantom crossing models. So if these dark energies accrete onto the BTZ BH, then in quintessence stage, BH mass increases upto a certain value (finite value) and then decreases to a certain finite value for phantom stage during whole evolution of the universe. We have shown these results graphically.


title: Diffusion of hidden charm mesons in hadronic medium

abstract: The drag and diffusion coefficients of a hot hadronic medium have been evaluated by using hidden charm mesons as probes. The matrix elements for the evaluation of these coefficients are calculated using an effective theory as well as from scattering lengths. Although the transport coefficients show a significant rise with temperature its effects on the suppression of $J/\psi$ in hadronic matter is not significant.


title: Offline Trace Checking of Quantitative Properties of Service-Based Applications

abstract: Service-based applications are often developed as compositions of partner services. A service integrator needs precise methods to specify the quality attributes expected by each partner service, as well as effective techniques to verify these attributes. In previous work, we identified the most common specification patterns related to provisioning service-based applications and developed an expressive specification language (SOLOIST) that supports them. SOLOIST is an extension of metric temporal logic with aggregate temporal modalities that can be used to write quantitative temporal properties. In this paper we address the problem of performing offline checking of service execution traces against quantitative requirements specifications written in SOLOIST. We present a translation of SOLOIST into CLTLB(D), a variant of linear temporal logic, and reduce the trace checking of SOLOIST to bounded satisfiability checking of CLTLB(D), which is supported by ZOT, an SMT-based verification toolkit. We detail the results of applying the proposed offline trace checking procedure to different types of traces, and compare its performance with previous work.


title: Inverse problem for extragalactic transport of ultra-high energy cosmic rays

abstract: The energy spectra and composition of ultra-high energy cosmic rays are changing in a course of propagation in the expanding Universe filled with background radiation. We developed a numerical code for solution of inverse problem for cosmic-ray transport equations that allows the determination of average source spectra of different nuclei from the cosmic ray spectra observed at the Earth. Employing this approach, the injection spectra of protons and Iron nuclei in extragalactic sources are found assuming that only these species are accelerated at the source. The data from the Auger experiment and the combined data from the Telescope Array + HiRes experiments are used to illustrate the method.


title: Jet Production in p-Pb Collisions

abstract: One of the major results from the study of high energy heavy ion collisions is the observation of jet quenching. The suppression of the number of jets observed in heavy ion collisions relative to pp collisions at the same energy scaled by the number of binary collisions, is attributed to partonic energy loss in the quark gluon plasma (QGP). However, cold nuclear matter effects due to the presence of a nucleus in the initial state could also influence this measurement. To disentangle these effects p-Pb collisions are studied, where QGP formation is not expected to occur and only cold nuclear matter effects are present. In addition to being an important baseline for understanding jet quenching, jets in p-Pb collisions may also be used to provide constraints on the nuclear parton distribution functions. Fully reconstructed jets measured using the ALICE tracking system and electro-magnetic calorimeter in p-Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV are reported. In addition to the spectra, studies of the jet fragmentation behavior in p-Pb collisions are also presented.


title: Embedded Markov chain approximations in Skorokhod topologies

abstract: In order to approximate a continuous time stochastic process by discrete time Markov chains one has several options to embed the Markov chains into continuous time processes. On the one hand there is the Markov embedding, which uses exponential waiting times. On the other hand each Skorokhod topology naturally suggests a certain embedding. These are the step function embedding for $J_1$, the linear interpolation embedding for $M_1$, the multi step embedding for $J_2$ and a more general embedding for $M_2$. We show that the convergence of the step function embedding in $J_1$ implies the convergence of the other embeddings in the corresponding topologies, respectively. For the converse statement a $J_1$-tightness condition for embedded Markov chains is given. The result relies on various representations of the Skorokhod topologies. Additionally it is shown that $J_1$ convergence is equivalent to the joint convergence in $M_1$ and $J_2$.


title: Four-Quark Effective Operators at Hadron Colliders

abstract: The robustness of translating effective operator constraints to BSM theories crucially depends on the mass and coupling of BSM particles. This is especially relevant for hadron colliders where the partonic centre of mass energy is around the typical energy scales of natural BSM theories. The caveats in applying the limits are discussed using Z' and G' models, illustrating the effects for a large class of models. This analysis shows that the applicability of effective operators mainly depends on the ratio of the transfer energy in the events and the mass scale of the full theory. Moreover, based on these results a method is developed to recast existing experimental limits on four-quark effective operators to the full theory parameter space.


title: Shockingly low water abundances in Herschel / PACS observations of low-mass protostars in Perseus

abstract: Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in the immediate protostellar environment. Herschel/PACS spectral maps of 22 objects in the Perseus molecular cloud were obtained as part of the `William Herschel Line Legacy' survey. Line emission from H$_\mathrm{2}$O, CO, and OH is tested against shock models from the literature. Observed line ratios are remarkably similar and do not show variations with source physical parameters. Observations show good agreement with the shock models when line ratios of the same species are compared. Ratios of various H$_\mathrm{2}$O lines provide a particularly good diagnostic of pre-shock gas densities, $n_\mathrm{H}\sim10^{5}$ cm$^{-3}$, in agreement with typical densities obtained from observations of the post-shock gas. The corresponding shock velocities, obtained from comparison with CO line ratios, are above 20 km\,s$^{-1}$. However, the observations consistently show one-to-two orders of magnitude lower H$_\mathrm{2}$O-to-CO and H$_\mathrm{2}$O-to-OH line ratios than predicted by the existing shock models. The overestimated model H$_\mathrm{2}$O fluxes are most likely caused by an overabundance of H$_\mathrm{2}$O in the models since the excitation is well-reproduced. Illumination of the shocked material by ultraviolet photons produced either in the star-disk system or, more locally, in the shock, would decrease the H$_\mathrm{2}$O abundances and reconcile the models with observations. Detections of hot H$_\mathrm{2}$O and strong OH lines support this scenario.


title: Equi-homogeneity, Assouad Dimension and Non-autonomous Dynamics

abstract: We show that self-similar sets arising from iterated function systems that satisfy the Moran open-set condition, a canonical class of fractal sets, are `equi-homogeneous'. This is a regularity property that, roughly speaking, means that at each fixed length-scale any two neighbourhoods of the set have covers of approximately equal cardinality. Self-similar sets are notable in that they are Ahlfors-David regular, which implies that their Assouad and box-counting dimensions coincide. More generally, attractors of non-autonomous iterated functions systems (where maps are allowed to vary between iterations) can have distinct Assouad and box-counting dimensions. Consequently the familiar notion of Ahlfors-David regularity is too strong to be useful in the analysis of this important class of sets, which include generalised Cantor sets and possess different dimensional behaviour at different length-scales. We further develop the theory of equi-homogeneity showing that it is a weaker property than Ahlfors-David regularity and distinct from any previously defined notion of dimensional equivalence. However, we show that if the upper and lower box-counting dimensions of an equi-homogeneous set are equal and `attained' in a sense we make precise then the lower Assouad, Hausdorff, packing, lower box-counting, upper box-counting and Assouad dimensions coincide. Our main results provide conditions under which the attractor of a non-autonomous iterated function system is equi-homogeneous and we use this to compute the Assouad dimension of a certain class of these highly non-trivial sets.


title: Probing Individual Star Forming Regions Within Strongly Lensed Galaxies at z > 1

abstract: Star formation occurs on physical scales corresponding to individual star forming regions, typically of order ~100 parsecs in size, but current observational facilities cannot resolve these scales within field galaxies beyond the local universe. However, the magnification from strong gravitational lensing allows us to measure the properties of these discrete star forming regions within galaxies in the distant universe. New results from multi-wavelength spectroscopic studies of a sample of extremely bright, highly magnified lensed galaxies are revealing the complexity of star formation on sub-galaxy scales during the era of peak star formation in the universe. We find a wide range of properties in the rest-frame UV spectra of individual galaxies, as well as in spectra that originate from different star forming regions within the same galaxy. Large variations in the strengths and velocity structure of Lyman-alpha and strong P Cygni lines such as C IV, and MgII provide new insights into the astrophysical relationships between extremely massive stars, the elemental abundances and physical properties of the nebular gas those stars ionize, and the galactic-scale outflows they power.


title: Hydrogen Diffusion and Stabilization in Single-crystal VO2 Micro/nanobeams by Direct Atomic Hydrogenation

abstract: We report measurements of the diffusion of atomic hydrogen in single crystalline VO2 micro/nanobeams by direct exposure to atomic hydrogen, without catalyst. The atomic hydrogen is generated by a hot filament, and the doping process takes place at moderate temperature (373 K). Undoped VO2 has a metal-to-insulator phase transition at ~340 K between a high-temperature, rutile, metallic phase and a low-temperature, monoclinic, insulating phase with a resistance exhibiting a semiconductor-like temperature dependence. Atomic hydrogenation results in stabilization of the metallic phase of VO2 micro/nanobeams down to 2 K, the lowest point we could reach in our measurement setup. Based on observing the movement of the hydrogen diffusion front in single crystalline VO2 beams, we estimate the diffusion constant for hydrogen along the c-axis of the rutile phase to be 6.7 x 10^{-10} cm^2/s at approximately 373 K, exceeding the value in isostructural TiO2 by ~ 38x. Moreover, we find that the diffusion constant along the c-axis of the rutile phase exceeds that along the equivalent a-axis of the monoclinic phase by at least three orders of magnitude. This remarkable change in kinetics must originate from the distortion of the "channels" when the unit cell doubles along this direction upon cooling into the monoclinic structure. Ab initio calculation results are in good agreement with the experimental trends in the relative kinetics of the two phases. This raises the possibility of a switchable membrane for hydrogen transport.


title: An explicit iterative method to solve generalized mixed equilibrium problem, variational inequality problem and hierarchical fixed point problem for a nearly nonexpansive mapping

abstract: In this paper, we introduce a new iterative method to find a common solution of a generalized mixed equilibrium problem, a variational inequality problem and a hierarchical fixed point problem for a demicontinuous nearly nonexpansive mapping. We prove that the proposed method converges strongly to a common solution of above problems under the suitable conditions. It is also noted that the main theorem is proved without usual demiclosedness condition. Also, under the appropriate assumptions on the control sequences and operators, our iterative method can be reduced to recent methods. So, the results here improve and extend some recent corresponding results given by many other authors.


title: The Panchromatic High-Resolution Spectroscopic Survey of Local Group Star Clusters - I. General Data Reduction Procedures for the VLT/X-shooter UVB and VIS arm

abstract: Our dataset contains spectroscopic observations of 29 globular clusters in the Magellanic Clouds and the Milky Way performed with VLT/X-shooter. Here we present detailed data reduction procedures for the VLT/X-shooter UVB and VIS arm. These are not restricted to our particular dataset, but are generally applicable to different kinds of X-shooter data without major limitation on the astronomical object of interest. ESO's X-shooter pipeline (v1.5.0) performs well and reliably for the wavelength calibration and the associated rectification procedure, yet we find several weaknesses in the reduction cascade that are addressed with additional calibration steps, such as bad pixel interpolation, flat fielding, and slit illumination corrections. Furthermore, the instrumental PSF is analytically modeled and used to reconstruct flux losses at slit transit and for optimally extracting point sources. Regular observations of spectrophotometric standard stars allow us to detect instrumental variability, which needs to be understood if a reliable absolute flux calibration is desired. A cascade of additional custom calibration steps is presented that allows for an absolute flux calibration uncertainty of less than ten percent under virtually every observational setup provided that the signal-to-noise ratio is sufficiently high. The optimal extraction increases the signal-to-noise ratio typically by a factor of 1.5, while simultaneously correcting for resulting flux losses. The wavelength calibration is found to be accurate to an uncertainty level of approximately 0.02 Angstrom. We find that most of the X-shooter systematics can be reliably modeled and corrected for. This offers the possibility of comparing observations on different nights and with different telescope pointings and instrumental setups, thereby facilitating a robust statistical analysis of large datasets.


title: Optical transmittance of multilayer graphene

abstract: We study the optical transmittance of multilayer graphene films up to 65 layers thick. By combing large-scale tight-binding simulation and optical measurement on CVD multilayer graphene, the optical transmission through graphene films in the visible region is found to be solely determined by the number of graphene layers. We argue that the optical transmittance measurement is more reliable in the determination of the number of layers than the commonly used Raman Spectroscopy. Moreover, optical transmittance measurement can be applied also to other 2D materials with weak van der Waals interlayer interaction.


title: On the p-regularized trust region subproblem

abstract: The $p$-regularized subproblem (p-RS) is a regularisation technique in computing a Newton-like step for unconstrained optimization, which globally minimizes a local quadratic approximation of the objective function while incorporating with a weighted regularisation term $\frac{\sigma}{p} \|x\|^p$. The global solution of the $p$-regularized subproblem for $p=3$, also known as the cubic regularization, has been characterized in literature. In this paper, we resolve both the global and the local non-global minimizers of (p-RS) for $p>2$ with necessary and sufficient optimality conditions. Moreover, we prove a parallel result of Mart\'{\i}nez \cite{Mar} that the (p-RS) for $p>2$, analogous to the trust region subproblem, can have at most one local non-global minimizer. When the (p-RS) is subject to a fixed number $m$ additional linear inequality constraints, we show that the uniqueness of the local solution of the (p-RS) (if exists at all), especially for $p=4$, can be applied to solve such an extension in polynomial time.


title: Solutions to the Navier-Stokes Equations with Mixed Boundary Conditions in Two-Dimensional Bounded Domains

abstract: In this paper we consider the system of the non-steady Navier-Stokes equations with mixed boundary conditions. We study the existence and uniqueness of a solution of this system. We define Banach spaces $X$ and $Y$, respectively, to be the space of "possible" solutions of this problem and the space of its data. We define the operator $\mathcal{N}:X\rightarrow Y$ and formulate our problem in terms of operator equations. Let $\mathbf{u}\in X$ and ${{\mathcal G}_{\mathcal P}}_{\mathbf{u}}: X\rightarrow Y$ be the Frechet derivative of $\mathcal{N}$ at $\mathbf{u}$. We prove that ${{\mathcal G}_{\mathcal P}}_{\mathbf{u}}$ is one-to-one and onto $Y$. Consequently, suppose that the system is solvable with some given data (the initial velocity and the right hand side). Then there exists a unique solution of this system for data which are small perturbations of the previous ones. Next result proved in the Appendix of this paper is $W^{2,2}$- regularity of solutions of steady Stokes system with mixed boundary condition for sufficiently smooth data.


title: Computable Stochastic Processes

abstract: The aim of this paper is to present an elementary computable theory of probability, random variables and stochastic processes. The probability theory is baed on existing approaches using valuations and lower integrals. Various approaches to random variables are discussed, including the approach based on completions in a Polish space. We apply the theory to the study of stochastic dynamical systems in discrete-time, and give a brief exposition of the Wiener process as a foundation for stochastic differential equations. The theory is based within the framework of type-two effectivity, so has an explicit direct link with Turing computation, and is expressed in a system of computable types and operations, so has a clean mathematical description.


title: High-resolution, H band Spectroscopy of Be Stars with SDSS-III/APOGEE: I. New Be Stars, Line Identifications, and Line Profiles

abstract: APOGEE has amassed the largest ever collection of multi-epoch, high-resolution (R~22,500), H-band spectra for B-type emission line (Be) stars. The 128/238 APOGEE Be stars for which emission had never previously been reported serve to increase the total number of known Be stars by ~6%. We focus on identification of the H-band lines and analysis of the emission peak velocity separations (v_p) and emission peak intensity ratios (V/R) of the usually double-peaked H I and non-hydrogen emission lines. H I Br11 emission is found to preferentially form in the circumstellar disks at an average distance of ~2.2 stellar radii. Increasing v_p toward the weaker Br12--Br20 lines suggests these lines are formed interior to Br11. By contrast, the observed IR Fe II emission lines present evidence of having significantly larger formation radii; distinctive phase lags between IR Fe II and H I Brackett emission lines further supports that these species arise from different radii in Be disks. Several emission lines have been identified for the first time including ~16895, a prominent feature in the spectra for almost a fifth of the sample and, as inferred from relatively large v_p compared to the Br11-Br20, a tracer of the inner regions of Be disks. Unlike the typical metallic lines observed for Be stars in the optical, the H-band metallic lines, such as Fe II 16878, never exhibit any evidence of shell absorption, even when the H I lines are clearly shell-dominated. The first known example of a quasi-triple-peaked Br11 line profile is reported for HD 253659, one of several stars exhibiting intra- and/or extra-species V/R and radial velocity variation within individual spectra. Br11 profiles are presented for all discussed stars, as are full APOGEE spectra for a portion of the sample.


title: Covariant asymmetric wave packet for a field-theoretical description of neutrino oscillations

abstract: We consider a class of models for the relativistic covariant wave packets which can be used as asymptotically free in and out states in the quantum field theoretical formalisms for description of the neutrino flavor oscillation phenomenon. We demonstrate that the new "asymmetric" wave packet (AWP) is an appropriate alternative for the more convenient "symmetric" wave packets, like the so-called relativistic Gaussian packet (RGP) widely used in the QFT-based approaches to neutrino oscillations. We show that RGP is not a particular case of AWP, although many properties of these models are almost identical in the quasistable regime. We discuss some features of AWP distinguishing it from RGP.


title: Class polynomials for some affine Hecke algebras

abstract: Class polynomials attached to affine Hecke algebras were first introduced by X.~He in \cite{He1}. They play an important role in the study of affine Deligne-Lusztig varieties. Motivated by \cite{He2}, we compute the class polynomials attached to an affine Hecke algebra of type (twisted) $\widetilde{A}_2$. Using these class polynomials we prove a conjecture of G\"{o}rtz-Haines-Kottwitz-Reuman for the general linear group, unitary group and division algebra of semisimple rank 2. Furthermore, we discuss some interesting patterns on affine Deligne-Lusztig varieties.


title: Efficient Coordinated Recovery of Sparse Channels in Massive MIMO

abstract: This paper addresses the problem of estimating sparse channels in massive MIMO-OFDM systems. Most wireless channels are sparse in nature with large delay spread. In addition, these channels as observed by multiple antennas in a neighborhood have approximately common support. The sparsity and common support properties are attractive when it comes to the efficient estimation of large number of channels in massive MIMO systems. Moreover, to avoid pilot contamination and to achieve better spectral efficiency, it is important to use a small number of pilots. We present a novel channel estimation approach which utilizes the sparsity and common support properties to estimate sparse channels and require a small number of pilots. Two algorithms based on this approach have been developed which perform Bayesian estimates of sparse channels even when the prior is non-Gaussian or unknown. Neighboring antennas share among each other their beliefs about the locations of active channel taps to perform estimation. The coordinated approach improves channel estimates and also reduces the required number of pilots. Further improvement is achieved by the data-aided version of the algorithm. Extensive simulation results are provided to demonstrate the performance of the proposed algorithms.


title: Quark number susceptibilities at finite chemical potential from fugacity expansion

abstract: Generalized quark number susceptibilities are expected to be good probes for the phase transitions in QCD and the search of a possible critical point. However, their computation in lattice QCD is plagued by the complex action problem which appears at finite chemical potential mu. In this work we explore the possibilities of an expansion in the fugacity parameter exp(mu beta) which has features that make, in particular quark number related bulk observables easily accessible. We present results at finite chemical potential for generalized susceptibilities up to the 4th order as well as their ratios and compare them to model calculations.


title: High Temperature Superconductivity in the Cuprates

abstract: The discovery of high temperature superconductivity in the cuprates in 1986 triggered a spectacular outpouring of creative and innovative scientific inquiry. Much has been learned over the ensuing 28 years about the novel forms of quantum matter that are exhibited in this strongly correlated electron system. This progress has been made possible by improvements in sample quality, coupled with the development and refinement of advanced experimental techniques. In part, avenues of inquiry have been motivated by theoretical developments, and in part new theoretical frameworks have been conceived to account for unanticipated experimental observations. An overall qualitative understanding of the nature of the superconducting state itself has been achieved, while profound unresolved issues have come into increasingly sharp focus concerning the astonishing complexity of the phase diagram, the unprecedented prominence of various forms of collective fluctuations, and the simplicity and insensitivity to material details of the "normal" state at elevated temperatures. New conceptual approaches, drawing from string theory, quantum information theory, and various numerically implemented approximate approaches to problems of strong correlations are being explored as ways to come to grips with this rich tableaux of interrelated phenomena.


title: The not-so-nonlinear nonlinearity of Einstein's equation

abstract: Many of the technical complications associated with the general theory of relativity ultimately stem from the nonlinearity of Einstein's equation. It is shown here that an appropriate choice of dynamical variables may be used to eliminate all such nonlinearities beyond a particular order: Both Landau-Lifshitz and tetrad formulations of Einstein's equation are obtained which involve only finite products of the unknowns and their derivatives. Considerable additional simplifications arise in physically-interesting cases where metrics becomes approximately Kerr or, e.g., plane waves, suggesting that the variables described here can be used to efficiently reformulate perturbation theory in a variety of contexts. In all cases, these variables are shown to have simple geometrical interpretations which directly relate the local causal structure associated with the metric of interest to the causal structure associated with a prescribed background. A new method to search for exact solutions is outlined as well.


title: Heavy-flavour production and multiplicity dependence in pp and p--Pb collisions with ALICE

abstract: The production of heavy quarks in pp collisions provides a precision test of perturbative QCD calculations at the LHC energies. More complex collision systems like p--Pb collisions allow studies of cold nuclear matter effects, such as modifications of the parton distribution functions at small x and of the $\kt$ broadening effect. We present the ALICE results of prompt D-meson production as a function of the charged-particle multiplicity, in pp and p--Pb collisions at $\sqrt{s}=7$ TeV and $\sqrt{s_{NN}}=5.02$ TeV respectively. The per-event yield of D mesons in different multiplicity and $\pt$ intervals are compared for pp and p--Pb collisions to study the contribution of multi-parton interactions to open-charm production. Angular correlations of prompt D mesons and heavy-flavour decay electrons with charged hadrons in pp and p-Pb collisions are also shown in different kinematic ranges and compared to pQCD models. These measurements provide information on the charm fragmentation processes, on cold nuclear matter effects on charm production, and on the participation of charm in the collective motion arising in small collision systems like p-Pb.


title: Quantized Method Solution for Various Fluid-Solid Reaction Models

abstract: Fluid-solid reactions exist in many chemical and metallurgical process industries. Several models describe these reactions such as volume reaction model, grain model, random pore model and nucleation model. These models give two nonlinear coupled partial differential equations (CPDE) that must be solved numerically. A new approximate solution technique (quantized method) has been introduced for some of these models in recent years. In this work, the various fluid-solid reaction models with their quantized and numerical solutions have been discussed.


title: Quark spin-orbit correlations

abstract: The proton spin puzzle issue focused the attention on the parton spin and orbital angular momentum contributions to the proton spin. However, a complete characterization of the proton spin structure requires also the knowledge of the parton spin-orbit correlation. We showed that this quantity can be expressed in terms of moments of measurable parton distributions. Using the available phenomenological information about the valence quarks, we concluded that this correlation is negative, meaning that the valence quark spin and kinetic orbital angular momentum are, in average, opposite. The quark spin-orbit correlation can also be expressed more intuitively in terms of relativistic phase-space distributions, which can be seen as the mother distributions of the standard generalized and transverse-momentum dependent parton distributions. We present here for the first time some examples of the general multipole decomposition of these phase-space distributions.


title: Characterization of curves in $C^{(2)}$

abstract: In this paper we characterize the irreducible curves lying in $C^{(2)}$. We prove that a curve $B$ has a degree one morphism to $C^{(2)}$ with image a curve of degree $d$ with irreducible preimage in $C\times C$ if and only if there exists an irreducible smooth curve $D$ and morphisms from $D$ to $C$ and $B$ of degrees $d$ and $2$ respectively forming a diagram which does not reduce.


title: On a model of a population with variable motility

abstract: We study a reaction-diffusion equation with a nonlocal reaction term that models a population with variable motility. We establish a global supremum bound for solutions of the equation. We investigate the asymptotic (long-time and long-range) behavior of the population. We perform a certain rescaling and prove that solutions of the rescaled problem converge locally uniformly to zero in a certain region and stay positive (in some sense) in another region. These regions are determined by two viscosity solutions of a related Hamilton-Jacobi equation.


title: Constraining the gravitational wave energy density of the Universe in the Range 0.1 Hz to 1 Hz using the Apollo Seismic Array

abstract: In this paper, we describe an analysis of Apollo era lunar seismic data that places an upper limit on an isotropic stochastic gravitational-wave background integrated over a year in the frequency range 0.1\,Hz -- 1\,Hz. We find that because the Moon's ambient noise background is much quieter than that of the Earth, significant improvements over an Earth based analysis were made. We find an upper limit of $\Omega_{\rm GW}<1.2\times 10^{5}$, which is three orders of magnitude smaller than a similar analysis of a global network of broadband seismometers on Earth and the best limits in this band to date. We also discuss the benefits of a potential Earth-Moon correlation search and compute the time-dependent overlap reduction function required for such an analysis. For this search, we find an upper limit an order of magnitude larger than the Moon-Moon search.


title: Galaxy And Mass Assembly (GAMA): The dependence of the galaxy luminosity function on environment, redshift and colour

abstract: We use 80922 galaxies in the Galaxy And Mass Assembly (GAMA) survey to measure the galaxy luminosity function (LF) in different environments over the redshift range 0.04<z<0.26. The depth and size of GAMA allows us to define samples split by colour and redshift to measure the dependence of the LF on environment, redshift and colour. We find that the LF varies smoothly with overdensity, consistent with previous results, with little environmental dependent evolution over the last 3 Gyrs. The modified GALFORM model predictions agree remarkably well with our LFs split by environment, particularly in the most overdense environments. The LFs predicted by the model for both blue and red galaxies are consistent with GAMA for the environments and luminosities at which such galaxies dominate. Discrepancies between the model and the data seen in the faint end of the LF suggest too many faint red galaxies are predicted, which is likely to be due to the over-quenching of satellite galaxies. The excess of bright blue galaxies predicted in underdense regions could be due to the implementation of AGN feedback not being sufficiently effective in the lower mass halos.


title: Continuous representation for shell models of turbulence

abstract: In this work we construct and analyze continuous hydrodynamic models in one space dimension, which are induced by shell models of turbulence. After Fourier transformation, such continuous models split into an infinite number of uncoupled subsystems, which are all identical to the same shell model. The two shell models, which allow such a construction, are considered: the dyadic (Desnyansky--Novikov) model with the intershell ratio $\lambda = 2^{3/2}$ and the Sabra model of turbulence with $\lambda = \sqrt{2+\sqrt{5}} \approx 2.058$. The continuous models allow understanding various properties of shell model solutions and provide their interpretation in physical space. We show that the asymptotic solutions of the dyadic model with Kolmogorov scaling correspond to the shocks (discontinuities) for the induced continuous solutions in physical space, and the finite-time blowup together with its viscous regularization follow the scenario similar to the Burgers equation. For the Sabra model, we provide the physical space representation for blowup solutions and intermittent turbulent dynamics.


title: A short proof of the multilinear Kakeya inequality

abstract: We give a short proof of a slightly weaker version of the multilinear Kakeya inequality proven by Bennett, Carbery, and Tao.


title: Dynamic message-passing approach for kinetic spin models with reversible dynamics

abstract: A method to approximately close the dynamic cavity equations for synchronous reversible dynamics on a locally tree-like topology is presented. The method builds on $(a)$ a graph expansion to eliminate loops from the normalizations of each step in the dynamics, and $(b)$ an assumption that a set of auxilary probability distributions on histories of pairs of spins mainly have dependencies that are local in time. The closure is then effectuated by projecting these probability distributions on $n$-step Markov processes. The method is shown in detail on the level of ordinary Markov processes ($n=1$), and outlined for higher-order approximations ($n>1$). Numerical validations of the technique are provided for the reconstruction of the transient and equilibrium dynamics of the kinetic Ising model on a random graph with arbitrary connectivity symmetry.


title: Perturbation-based inference for diffusion processes: Obtaining effective models from multiscale data

abstract: We consider the inference problem for parameters in stochastic differential equation models from discrete time observations (e.g. experimental or simulation data). Specifically, we study the case where one does not have access to observations of the model itself, but only to a perturbed version which converges weakly to the solution of the model. Motivated by this perturbation argument, we study the convergence of estimation procedures from a numerical analysis point of view. More precisely, we introduce appropriate consistency, stability, and convergence concepts and study their connection. It turns out that standard statistical techniques, such as the maximum likelihood estimator, are not convergent methodologies in this setting, since they fail to be stable. Due to this shortcoming, we introduce and analyse a novel inference procedure for parameters in stochastic differential equation models which turns out to be convergent. As such, the method is particularly suited for the estimation of parameters in effective (i.e. coarse-grained) models from observations of the corresponding multiscale process. We illustrate these theoretical findings via several numerical examples.


title: Repr\'esentations modulo $p$ de $\mathrm{GL}(m,D)$, $D$ alg\`ebre \`a division sur un corps local

abstract: This Ph.D. thesis belongs to the realm of mod $p$ representation theory of $p$-adic groups. The main object of study is the inner form of the general linear group $\mathrm{GL}(m,D)$ where $D$ is a division algebra over a non-Archimedean local field. The first part focuses on the rank $1$ case ($m=2$): it develops in detail the irreducibility criterion for parabolically induced representations and a classification of irreducible admissible smooth representations "\`a la Barthel-Livn\'e". Then one turns to the generalized Steinberg representations: building on ideas of Grosse-Kl\"onne and Herzig, this chapter manages to get a result that is valid for any reductive group over a non-Archimedean local field. Finally the last part exploits the two previous ones, and by overcoming some combinatorial difficulties that arise in the rank $>1$ case, one gets a classification "\`a la Herzig" for $\mathrm{GL}(3,D)$ and $\mathrm{GL}(m,D)$, $m>3$, with some technical assumptions on $D$.


title: Position Auctions with Externalities and Brand Effects

abstract: This paper presents models for predicted click-through rates in position auctions that take into account two possibilities that are not normally considered---that the identities of ads shown in other positions may affect the probability that an ad in a particular position receives a click (externalities) and that some ads may be less adversely affected by being shown in a lower position than others (brand effects). We present a general axiomatic methodology for how click probabilities are affected by the qualities of the ads in the other positions, and illustrate that using these axioms will increase revenue as long as higher quality ads tend to be ranked ahead of lower quality ads. We also present appropriate algorithms for selecting the optimal allocation of ads when predicted click-through rates are governed by either the models of externalities or brand effects that we consider. Finally, we analyze the performance of a greedy algorithm of ranking the ads by their expected cost-per-1000-impressions bids when the true click-through rates are governed by our model of predicted click-through rates with brand effects and illustrate that such an algorithm will potentially cost as much as half of the total possible social welfare.


title: Cofree coalgebras over operads and representative functions

abstract: We give a recursive formula to compute the cofree coalgebra P^\vee(C) over any colored operad P in Set, CGHaus or (dg)Vect. The construction is closed to that of Smith but different. We use a more conceptual approach to simplify the proofs that P^\vee is the cofree P-coalgebra functor and also the comonad generating P-coalgebras. In a second part, when P is a linear or dg-operad over a field, we generalize the notion of representative functions of Block & Leroux and prove that P^\vee(C) is simply the subobject of representative elements in the "completed P-algebra" P^\wedge(C). This says that our recursion (as well as that of Smith) stops at the first step.


title: Compute Less to Get More: Using ORC to Improve Sparse Filtering

abstract: Sparse Filtering is a popular feature learning algorithm for image classification pipelines. In this paper, we connect the performance of Sparse Filtering with spectral properties of the corresponding feature matrices. This connection provides new insights into Sparse Filtering; in particular, it suggests early stopping of Sparse Filtering. We therefore introduce the Optimal Roundness Criterion (ORC), a novel stopping criterion for Sparse Filtering. We show that this stopping criterion is related with pre-processing procedures such as Statistical Whitening and demonstrate that it can make image classification with Sparse Filtering considerably faster and more accurate.


title: Imaging an Event Horizon: Mitigation of Scattering Toward Sagittarius A*

abstract: The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.


title: Evolution of Dark Matter Halos and their Radio Emissions

abstract: Radio synchrotron emission is expected as a natural by-product of the self-annihilation of super-symmetric dark matter particles. In this work we discuss the general properties of the radio emission expected in a wide range of dark matter halos, from local dwarf spheroidal galaxies to large and distant galaxy clusters with the aim to determine the neutralino dark matter detection prospects of the Square Kilometre Array (SKA). The analysis of the SKA detection of dark matter(DM)-induced radio emission is presented for structures spanning a wide range of masses and redshifts, and we also analyze the limits that the SKA can set on the thermally averaged neutralino annihilation cross-section in the event of non-detection. To this aim, we construct a model of the redshift evolution of the radio emissions of dark matter halos and apply it to generate predicted fluxes from a range of neutralino masses and annihilation channels for the dark matter halos surrounding dwarf galaxies, galaxies and galaxy clusters. Using the available SKA performance predictions and its ability to determine an independent measure of the magnetic field in cosmic structures, we explore both the detailed detection prospects and the upper-bounds that might be placed on the neutralino annihilation cross-section in the event of non-detection. We find that the SKA can access a neutralino parameter space far larger than that of any preceding indirect-detection experiment, also improving on the realistic CTA detection prospects, with the possibility of setting cross-section upper-bounds up to four orders of magnitude below the thermal relic density bound. Additionally, we find that neutralino radio emissions carry redshift-independent signatures of the dominant annihilation channel and of neutralino mass, offering therefore a means to identify such non-thermal emissions within the observing frequency range of the SKA.


title: Integrable turbulence and formation of rogue waves

abstract: In the framework of the focusing Nonlinear Schrodinger (NLS) equation we study numerically the nonlinear stage of the modulation instability (MI) of the condensate. As expected, the development of the MI leads to formation of "integrable turbulence" [V.E. Zakharov, Turbulence in integrable systems, Stud. in Appl. Math. 122, no. 3, 219-234, (2009)]. We study the time evolution of it's major characteristics averaged across realizations of initial data - the condensate solution seeded by small random noise with fixed statistical properties. The measured quantities are: (1) wave-action spectrum and spatial correlation function, (2) the probability density function (PDF) of wave amplitudes and their momenta, and (3) kinetic and potential energies.


title: Conformal symmetries of gravity from asymptotic methods: further developments

abstract: In this thesis, the symmetry structure of gravitational theories at null infinity is studied further, in the case of pure gravity in four dimensions and also in the case of Einstein-Yang-Mills theory in $d$ dimensions with and without a cosmological constant. The first part of this thesis is devoted to the presentation of asymptotic methods (symmetries, solution space and surface charges) applied to gravity in the case of the BMS gauge in three and four spacetime dimensions. The second part of this thesis contains the original contributions. Firstly, it is shown that the enhancement from Lorentz to Virasoro algebra also occurs for asymptotically flat spacetimes defined in the sense of Newman-Unti. As a first application, the transformation laws of the Newman-Penrose coefficients characterizing solution space of the Newman-Unti approach are worked out, focusing on the inhomogeneous terms that contain the information about central extensions of the theory. These transformations laws make the conformal structure particularly transparent, and constitute the main original result of the thesis. Secondly, asymptotic symmetries of the Einstein-Yang-Mills system with or without cosmological constant are explicitly worked out in a unified manner in $d$ dimensions. In agreement with a recent conjecture, a Virasoro-Kac-Moody type algebra is found not only in three dimensions but also in the four dimensional asymptotically flat case. These two parts of the thesis are supplemented by appendices.


title: The chemical composition of red giants in 47 Tucanae I: Fundamental parameters and chemical abundance patterns

abstract: Context: The study of chemical abundance patterns in globular clusters is of key importance to constrain the different candidates for intra-cluster pollution of light elements. Aims: We aim at deriving accurate abundances for a large range of elements in the globular cluster 47 Tucanae (NGC 104) to add new constraints to the pollution scenarios for this particular cluster, expanding the range of previously derived element abundances. Methods: Using tailored 1D LTE atmospheric models together with a combination of equivalent width measurements, LTE, and NLTE synthesis we derive stellar parameters and element abundances from high-resolution, high signal-to-noise spectra of 13 red giant stars near the tip of the RGB. Results: We derive abundances of a total 27 elements (O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ru, Ba, La, Ce, Pr, Nd, Eu, Dy). Departures from LTE were taken into account for Na, Al and Ba. We find a mean [Fe/H] = $-0.78\pm0.07$ and $[\alpha/{\rm Fe}]=0.34\pm0.03$ in good agreement with previous studies. The remaining elements show good agreement with the literature, but the inclusion of NLTE for Al has a significant impact on the behaviour of this key element. Conclusions: We confirm the presence of an Na-O anti-correlation in 47 Tucanae found by several other works. Our NLTE analysis of Al shifts the [Al/Fe] to lower values, indicating that this may be overestimated in earlier works. No evidence for an intrinsic variation is found in any of the remaining elements.


title: Lurking in Social Networks: Topology-based Analysis and Ranking Methods

abstract: The massive presence of silent members in online communities, the so-called lurkers, has long attracted the attention of researchers in social science, cognitive psychology, and computer-human interaction. However, the study of lurking phenomena represents an unexplored opportunity of research in data mining, information retrieval and related fields. In this paper, we take a first step towards the formal specification and analysis of lurking in social networks. We address the new problem of lurker ranking and propose the first centrality methods specifically conceived for ranking lurkers in social networks. Our approach utilizes only the network topology without probing into text contents or user relationships related to media. Using Twitter, Flickr, FriendFeed and GooglePlus as cases in point, our methods' performance was evaluated against data-driven rankings as well as existing centrality methods, including the classic PageRank and alpha-centrality. Empirical evidence has shown the significance of our lurker ranking approach, and its uniqueness in effectively identifying and ranking lurkers in an online social network.


title: Differentially Private Exponential Random Graphs

abstract: We propose methods to release and analyze synthetic graphs in order to protect privacy of individual relationships captured by the social network. Proposed techniques aim at fitting and estimating a wide class of exponential random graph models (ERGMs) in a differentially private manner, and thus offer rigorous privacy guarantees. More specifically, we use the randomized response mechanism to release networks under $\epsilon$-edge differential privacy. To maintain utility for statistical inference, treating the original graph as missing, we propose a way to use likelihood based inference and Markov chain Monte Carlo (MCMC) techniques to fit ERGMs to the produced synthetic networks. We demonstrate the usefulness of the proposed techniques on a real data example.


title: Higher order recurrence relation for exceptional Charlier, Meixner, Hermite and Laguerre orthogonal polynomials

abstract: In this paper we prove in a constructing way that exceptional Charlier, Meixner, Hermite and Laguerre polynomials satisfy higher order recurrence relations. Our conjecture is that the recurrence relations provided in this paper have minimal order.


title: A Mixtures-of-Experts Framework for Multi-Label Classification

abstract: We develop a novel probabilistic approach for multi-label classification that is based on the mixtures-of-experts architecture combined with recently introduced conditional tree-structured Bayesian networks. Our approach captures different input-output relations from multi-label data using the efficient tree-structured classifiers, while the mixtures-of-experts architecture aims to compensate for the tree-structured restrictions and build a more accurate model. We develop and present algorithms for learning the model from data and for performing multi-label predictions on future data instances. Experiments on multiple benchmark datasets demonstrate that our approach achieves highly competitive results and outperforms the existing state-of-the-art multi-label classification methods.


title: Fluid flow enhances the effectiveness of toxin export by aquatic microorganisms: a first-passage perspective on microvilli and the concentration boundary layer

abstract: A central challenge for organisms during development is determining a means to efficiently export toxic molecules from inside the developing embryo. For aquatic microorganisms, the strategies employed should be robust with respect to the variable ocean environment and limit the chances that exported toxins are reabsorbed. As a result, the problem of toxin export is closely related to the physics of mass transport in a fluid. In this paper we consider a model first-passage problem for the uptake of exported toxins by a spherical embryo. By considering how macroscale fluid turbulence manifests itself on the microscale of the embryo, we determine that fluid flow enhances the effectiveness of toxin export as compared to the case of diffusion-limited transport. In the regime of large P\'eclet number, a perturbative solution of the advection-diffusion equation reveals that a concentration boundary layer forms at the surface of the embryo. The model results suggest a functional role for cell surface roughness in the export process, with the thickness of the concentration boundary layer setting the length scale for cell membrane protrusions known as microvilli. We highlight connections between the model results and experiments on the development of sea urchin embryos.


title: Can Sol's Explanation for the Evolution of Animal Innovation Account for Human Innovation?

abstract: Sol argues that innovation propensity is not a specialized adaptation resulting from targeted selection but an instance of exaptation because selection cannot act on situations that are only encountered once. In exaptation, a trait that originally evolved to solve one problem is co-opted to solve a new problem; thus the trait or traits in question must be necessary and sufficient to solve the new problem. Sol claims that traits such as persistence and neophilia, are necessary and sufficient for animal innovation, which is a matter of trial and error. We suggest that this explanation does not extend to human innovation, which involves strategy, logic, intuition, and insight, and requires traits that evolved, not as a byproduct of some other function, but for the purpose of coming up with adaptive responses to environmental variability itself. We point to an agent based model that indicates the feasibility of two such proposed traits: (1) chaining, the ability to construct complex thoughts from simple ones, and (2) contextual focus, the ability to shift between convergent and divergent modes of thought. We agree that there is a sense in which innovation is exaptation--it occurs when an existing object or behaviour is adapted to new needs or tastes--and refer to a mathematical model of biological and cultural exaltation. We conclude that much is gained by comparing and contrasting animal and human innovation.


title: Coexistence of superconductivity and density wave in $\mathrm{Ba_2Ti_2Fe_2As_4O}$: an optical spectroscopy study

abstract: We performed an optical spectroscopy measurement on single crystals of $\mathrm{Ba_2Ti_2Fe_2As_4O}$, which is a newly discovered superconductor showing a coexistence of superconductivity and density wave order. The study reveals a spectral change related to the formation density wave energy gap below $T_{DW}\approx$125 K, leading to the removal of about half of the multiple Fermi surface sheets. The ratio of 2$\Delta_{DW}$/$k_B T_{DW}\approx$ 11.9 is considerably larger than the mean-field value based on the weak-coupling BCS theory. At the lowest temperature in the superconducting state, we observed opening of superconducting energy gaps $\Delta_1(0) =3.4$ meV and $\Delta_2(0)=7.9$ meV. The properties of superconducting state in $\mathrm{Ba_2Ti_2Fe_2As_4O}$ are similar to that in $\mathrm{BaFe_{1.85}Co_{0.15}As_2}$.


title: Algebraic distance for anisotropic diffusion problems: multilevel results

abstract: In this paper we motivate, discuss the implementation and present the resulting numerics for a new definition of strength of connection which is based on the notion of algebraic distance. This algebraic distance measure, combined with compatible relaxation, is used to choose suitable coarse grids and accurate interpolation operators for algebraic multigrid algorithms. The main tool of the proposed measure is the least squares functional defined using a set of relaxed test vectors. The motivating application is the anisotropic diffusion problem, in particular problems with non-grid aligned anisotropy. We demonstrate numerically that the measure yields a robust technique for determining strength of connectivity among variables, for both two-grid and multigrid solvers. %We illustrate the use of the measure to construct, in addition, an adaptive aggregation form of interpolation for the targeted anisotropic problems. %Our approach is not a two-level approach -- we provide preliminary results that show its extendability to multigrid. The proposed algebraic distance measure can also be used in any other coarsening procedure, assuming a rich enough set of the near-kernel components of the matrix for the targeted system is known or computed.


title: Confinement with Perturbation Theory, after All?

abstract: I call attention to the possibility that QCD bound states (hadrons) could be derived using rigorous Hamiltonian, perturbative methods. Solving Gauss' law for $A^0$ with a non-vanishing boundary condition at spatial infinity gives an \order{\alpha_s^0} linear potential for color singlet $q\bar q$ and $qqq$ states. These states are Poincar\'e and gauge covariant and thus can serve as initial states of a perturbative expansion, replacing the conventional free $in$ and $out$ states. The coupling freezes at $\alpha_s(0)\simeq 0.5$, allowing reasonable convergence. The \order{\alpha_s^0} bound states have a sea of $q\bar q$ pairs, while transverse gluons contribute only at \order{\alpha_s}. Pair creation in the linear $A^0$ potential leads to string breaking and hadron loop corrections. These corrections give finite widths to excited states, as required by unitarity. Several of these features have been verified analytically in $D=1+1$ dimensions, and some in $D=3+1$.