title: Inferring Passenger Type from Commuter Eigentravel Matrices

abstract: A sufficient knowledge of the demographics of a commuting public is essential in formulating and implementing more targeted transportation policies, as commuters exhibit different ways of traveling. With the advent of the Automated Fare Collection system (AFC), probing the travel patterns of commuters has become less invasive and more accessible. Consequently, numerous transport studies related to human mobility have shown that these observed patterns allow one to pair individuals with locations and/or activities at certain times of the day. However, classifying commuters using their travel signatures is yet to be thoroughly examined. Here, we contribute to the literature by demonstrating a procedure to characterize passenger types (Adult, Child/Student, and Senior Citizen) based on their three-month travel patterns taken from a smart fare card system. We first establish a method to construct distinct commuter matrices, which we refer to as eigentravel matrices, that capture the characteristic travel routines of individuals. From the eigentravel matrices, we build classification models that predict the type of passengers traveling. Among the models explored, the gradient boosting method (GBM) gives the best prediction accuracy at 76%, which is 84% better than the minimum model accuracy (41%) required vis-\`a-vis the proportional chance criterion. In addition, we find that travel features generated during weekdays have greater predictive power than those on weekends. This work should not only be useful for transport planners, but for market researchers as well. With the awareness of which commuter types are traveling, ads, service announcements, and surveys, among others, can be made more targeted spatiotemporally. Finally, our framework should be effective in creating synthetic populations for use in real-world simulations that involve a metropolitan's public transport system.

categories: physics.soc-ph cs.CY physics.data-an stat.AP stat.ML


title: A modular modelling framework for hypotheses testing in the simulation of urbanisation

abstract: In this paper, we present a modelling experiment developed to study systems of cities and processes of urbanisation in large territories over long time spans. Building on geographical theories of urban evolution, we rely on agent-based models to 1/ formalise complementary and alternative hypotheses of urbanisation and 2/ explore their ability to simulate observed patterns in a virtual laboratory. The paper is therefore divided into two sections : an overview of the mechanisms implemented to represent competing hypotheses used to simulate urban evolution; and an evaluation of the resulting model structures in their ability to simulate - efficiently and parsimoniously - a system of cities (the Former Soviet Union) over several periods of time (before and after the crash of the USSR). We do so using a modular framework of model-building and evolutionary algorithms for the calibration of several model structures. This project aims at tackling equifinality in systems dynamics by confronting different mechanisms with similar evaluation criteria. It enables the identification of the best-performing models with respect to the chosen criteria by scanning automatically the parameter along with the space of model structures (as combinations of modelled dynamics).


title: Search for Pulsed Emission in Archival VERITAS Data

abstract: Since the 2011 VERITAS discovery of very high energy (VHE; E>100 GeV) gamma rays from the Crab pulsar, there has been concerted effort by the gamma-ray astrophysics community to detect other pulsars in the VHE band in order to place better constraints on emission models. Pulsar modelling demonstrates that much of the magnetosphere is opaque to VHE photons, limiting emission regions to the outer magnetosphere or beyond the light cylinder. The locations of 19 known pulsars have been observed by VERITAS since full observations began in 2007 with 11 locations having more than 20 hours of observations. Observations of VHE emission from more sources could provide key data to help constrain current models of emission location and mechanisms. We present the status of the ongoing VERITAS program searching for pulsed emission in archival data.


title: Photometry of the Stingray Nebula (V839 Ara) from 1889-2015 Across the Ionization of Its Planetary Nebula

abstract: Up until around 1980, the Stingray was an ordinary B1 post-AGB star, but then it suddenly sprouted bright emission lines like in a planetary nebula (PN), and soon after this the Hubble Space Telescope (HST) discovered a small PN around the star, so apparently we have caught a star in the act of ionizing a PN. We report here on a well-sampled light curve from 1889 to 2015, with unique coverage of the prior century plus the entire duration of the PN formation plus three decades of its aftermath. Surprisingly, the star anticipated the 1980's ionization event by declining from B=10.30 in 1889 to B=10.76 in 1980. Starting in 1980, the central star faded fast, at a rate of 0.20 mag/year, reaching B=14.64 in 1996. This fast fading is apparently caused by the central star shrinking in size. From 1994-2015, the V-band light curve is almost entirely from the flux of two bright [OIII] emission lines from the unresolved nebula, and it shows a consistent decline at a rate of 0.090 mag/year. This steady fading (also seen in the radio and infrared) has a time scale equal to that expected for ordinary recombination within the nebula, immediately after a short-duration ionizing event in the 1980s. We are providing the first direct measure of the rapidly changing luminosity of the central star on both sides of a presumed thermal pulse in 1980, with this providing a strong and critical set of constraints, and these are found to sharply disagree with theoretical models of PN evolution.


title: Non-Hermitian extension of gauge theories and implications for neutrino physics

abstract: An extension of QED is considered in which the Dirac fermion has both Hermitian and anti-Hermitian mass terms, as well as both vector and axial-vector couplings to the gauge field. Gauge invariance is restored when the Hermitian and anti-Hermitian masses are of equal magnitude, and the theory reduces to that of a single massless Weyl fermion. An analogous non-Hermitian Yukawa theory is considered, and it is shown that this model can explain the smallness of the light-neutrino masses and provide an additional source of leptonic CP violation.


title: On the feasibility of ab initio electronic structure calculations for Cu using a single s orbital basis

abstract: The accuracy of a single s-orbital representation of Cu towards enabling multi-thousand atom ab initio calculations of electronic structure is evaluated in this work. If an electrostatic compensation charge of approximately 0.3 electrons per atom is used in this basis representation of copper, the electronic transmission in bulk and nanocrystalline Cu compares accurately to that obtained with a Double Zeta Polarized basis set. The use of this representation is analogous to the use of single band effective mass representation for semiconductor electronic structure. With a basis of just one s-orbital per Cu atom, the representation is extremely computationally efficient and can be used to provide much needed ab initio insight into electronic transport in nanocrystalline Cu interconnects at realistic dimensions.


title: Performance Comparisons of Geographic Routing Protocols in Mobile Ad Hoc Networks

abstract: Geographic routing protocols greatly reduce the requirements of topology storage and provide flexibility in the accommodation of the dynamic behavior of mobile ad hoc networks. This paper presents performance evaluations and comparisons of two geographic routing protocols and the popular AODV protocol. The tradeoffs among the average path reliabilities, average conditional delays, average conditional numbers of hops, and area spectral efficiencies and the effects of various parameters are illustrated for finite ad hoc networks with randomly placed mobiles. This paper uses a dual method of closed-form analysis and simple simulation that is applicable to most routing protocols and provides a much more realistic performance evaluation than has previously been possible. Some features included in the new analysis are shadowing, exclusion and guard zones, distance-dependent fading, and interference correlation.


title: Inadequate experimental methods and erroneous epilepsy diagnostic criteria result in confounding acquired focal epilepsy with genetic absence epilepsy

abstract: Here we provide a thorough discussion of the study conducted by Rodgers et al. (J Neurosci. 2015; 35(24):9194-204. doi: 10.1523/JNEUROSCI.0919-15.2015) to investigate focal seizures and acquired epileptogenesis induced by head injury in the rat. This manuscript serves as supplementary document for our letter to the Editor to appear in the Journal of Neuroscience. We find that the subject article suffers from poor experimental design, very selective consideration of antecedent literature, and application of inappropriate epilepsy diagnostic criteria which, together, lead to unwarranted conclusions.


title: General and exact approach to percolation on random graphs

abstract: We present a comprehensive and versatile theoretical framework to study site and bond percolation on clustered and correlated random graphs. Our contribution can be summarized in three main points. (i) We introduce a set of iterative equations that solve the exact distribution of the size and composition of components in finite size quenched or random multitype graphs. (ii) We define a very general random graph ensemble that encompasses most of the models published to this day, and also that permits to model structural properties not yet included in a theoretical framework. Site and bond percolation on this ensemble is solved exactly in the infinite size limit using probability generating functions [i.e., the percolation threshold, the size and the composition of the giant (extensive) and small components]. Several examples and applications are also provided. (iii) Our approach can be adapted to model interdependent graphs---whose most striking feature is the emergence of an extensive component via a discontinuous phase transition---in an equally general fashion. We show how a graph can successively undergo a continuous then a discontinuous phase transition, and preliminary results suggest that clustering increases the amplitude of the discontinuity at the transition.


title: Fast Clustering and Topic Modeling Based on Rank-2 Nonnegative Matrix Factorization

abstract: The importance of unsupervised clustering and topic modeling is well recognized with ever-increasing volumes of text data. In this paper, we propose a fast method for hierarchical clustering and topic modeling called HierNMF2. Our method is based on fast Rank-2 nonnegative matrix factorization (NMF) that performs binary clustering and an efficient node splitting rule. Further utilizing the final leaf nodes generated in HierNMF2 and the idea of nonnegative least squares fitting, we propose a new clustering/topic modeling method called FlatNMF2 that recovers a flat clustering/topic modeling result in a very simple yet significantly more effective way than any other existing methods. We implement highly optimized open source software in C++ for both HierNMF2 and FlatNMF2 for hierarchical and partitional clustering/topic modeling of document data sets. Substantial experimental tests are presented that illustrate significant improvements both in computational time as well as quality of solutions. We compare our methods to other clustering methods including K-means, standard NMF, and CLUTO, and also topic modeling methods including latent Dirichlet allocation (LDA) and recently proposed algorithms for NMF with separability constraints. Overall, we present efficient tools for analyzing large-scale data sets, and techniques that can be generalized to many other data analytics problem domains.


title: Numerical study of the optical nonlinearity of doped and gapped graphene: From weak to strong field excitation

abstract: Numerically solving the semiconductor Bloch equations within a phenomenological relaxation time approximation, we extract both the linear and nonlinear optical conductivities of doped graphene and gapped graphene under excitation by a laser pulse. We discuss in detail the dependence of second harmonic generation, third harmonic generation, and the Kerr effects on the doping level, the gap, and the electric field amplitude. The numerical results for weak electric fields agree with those calculated from available analytic perturbation formulas. For strong electric fields when saturation effects are important, all the effective third order nonlinear response coefficients show a strong field dependence.


title: Pitt inequalities and restriction theorems for the Fourier transform

abstract: We prove new Pitt inequalities for the Fourier transforms with radial and non-radial weights using weighted restriction inequalities for the Fourier transform on the sphere. We also prove new Riemann-Lebesgue estimates and versions of the uncertainty principle for the Fourier transform.


title: Control principles of metabolic networks

abstract: Deciphering the control principles of metabolism and its interaction with other cellular functions is central to biomedicine and biotechnology. Yet, understanding the efficient control of metabolic fluxes remains elusive for large-scale metabolic networks. Existing methods either require specifying a cellular objective or are limited to small networks due to computational complexity. Here we develop an efficient computational framework for flux control by introducing a complete set of flux coupling relations. We analyze 23 metabolic networks from all kingdoms of life, and identify the driver reactions facilitating their control on a large scale. We find that most unicellular organisms require less extensive control than multicellular organisms. The identified driver reactions are under strong transcriptional regulation in Escherichia coli. In human cancer cells driver reactions play pivotal roles in tumor development, representing potential therapeutic targets. The proposed framework helps us unravel the regulatory principles of complex diseases and design novel engineering strategies at the interface of gene regulation, signaling, and metabolism.


title: Stochastic Frontier I & D of fractal dimensions for technological innovation

abstract: This paper presents an analysis of the study variables such as gdp, employment levels, the level of R & D and technology that will serve as the basis for stochastic modeling of production possibilities frontier in the goodness of fractal dimensions Ex Ante and Ex Post a priori to determine the levels of causality immediately and check its accuracy and power of indexing, using high frequency data and thus address the response this assumption of stochastic frontiers with level N of partitions in time.


title: Impact of Artificial Intelligence on Economic Theory

abstract: Artificial intelligence has impacted many aspects of human life. This paper studies the impact of artificial intelligence on economic theory. In particular we study the impact of artificial intelligence on the theory of bounded rationality, efficient market hypothesis and prospect theory.


title: The Poker-Litigation Game

abstract: Is litigation a serious search for truth or simply a game of skill or luck? Although the process of litigation has been modeled as a Prisoner's Dilemma, as a War of Attrition, as a Game of Chicken and even as a simple coin toss, no one has formally modeled litigation as a game of poker. This paper is the first to do so. We present a simple "poker-litigation game" and find the optimal strategy for playing this game.


title: Assessing Consistency of Consumer Confidence Data using Dynamic Latent Class Analysis

abstract: In many countries information on expectations collected through consumer confidence surveys are used in macroeconomic policy formulation. Unfortunately, before doing so, the consistency of responses is often not taken into account, leading to biases creeping in and affecting the reliability of the indices hence created. This paper describes how latent class analysis may be used to check the consistency of responses and ensure a parsimonious questionnaire. In particular, we examine how temporal changes may be incorporated into the model. Our methodology is illustrated using three rounds of Consumer Confidence Survey (CCS) conducted by Reserve Bank of India (RBI).


title: Dynamic Model of the Price Dispersion of Homogeneous Goods

abstract: Presented is an analytic microeconomic model of the temporal price dispersion of homogeneous goods in polypoly markets. This new approach is based on the idea that the price dispersion has its origin in the dynamics of the purchase process. The price dispersion is determined by the chance that demanded and supplied product units meet in a given price interval. It can be characterized by a fat-tailed Laplace distribution for short and by a lognormal distribution for long time horizons. Taking random temporal variations of demanded and supplied units into account both the mean price and also the standard deviation of the price dispersion are governed by a lognormal distribution. A comparison with empirical investigations confirms the model statements.


title: Wealth distribution across communities of adaptive financial agents

abstract: This paper studies the trading volumes and wealth distribution of a novel agent-based model of an artificial financial market. In this model, heterogeneous agents, behaving according to the Von Neumann and Morgenstern utility theory, may mutually interact. A Tobin-like tax (TT) on successful investments and a flat tax are compared to assess the effects on the agents' wealth distribution. We carry out extensive numerical simulations in two alternative scenarios: i) a reference scenario, where the agents keep their utility function fixed, and ii) a focal scenario, where the agents are adaptive and self-organize in communities, emulating their neighbours by updating their own utility function. Specifically, the interactions among the agents are modelled through a directed scale-free network to account for the presence of community leaders, and the herding-like effect is tested against the reference scenario. We observe that our model is capable of replicating the benefits and drawbacks of the two taxation systems and that the interactions among the agents strongly affect the wealth distribution across the communities. Remarkably, the communities benefit from the presence of leaders with successful trading strategies, and are more likely to increase their average wealth. Moreover, this emulation mechanism mitigates the decrease in trading volumes, which is a typical drawback of TTs.


title: Tax Bond Creation Using a Structural Model and its Extensions

abstract: This article describes and explores taxes and debt in finance. Here a situation is thought about, where tax payments would qualify to be considered as debt. Using this principle we can infer that it is possible to create and price a type of bond (Tax Normalization Guarantee) for companies, which would allow them to enter in temporary tax breaks to allow them to free capital. Finally it is explored a way to structure these bonds in financial products and valuate them.


title: Imaginary Time Correlations for a High-Density two-dimensional Electron Gas

abstract: We evaluate imaginary time density-density correlation functions for a two-dimensional homogeneous electron gas using the phaseless auxiliary field quantum Monte Carlo method. We show that such methodology, once equipped with suitable numerical stabilization techniques necessary to deal with exponentials, products and inversions of large matrices, gives access to the calculation of imaginary time correlation functions for medium-sized systems; we present simulations of a number up to 42 correlated fermions in the continuum, using up to 300 plane waves as basis set elements. We discuss the numerical stabilization techniques and the computational complexity of the methodology. We perform the inverse Laplace transform of the obtained density-density correlation functions, assessing the ability of the phaseless auxiliary field quantum Monte Carlo method to evaluate dynamical properties of many-fermion systems.


title: Light Efficient Flutter Shutter

abstract: Flutter shutter is a technique in which the exposure is chopped into segments and light is only integrated part of the time. By carefully selecting the chopping sequence it is possible to better condition the data for reconstruction problems such as motion deblurring, focal sweeping, and compressed sensing. The partial exposure trades better conditioning for less energy. In problems such as motion deblurring the available energy is what caused the problem in the first place (as strong illumination allows short exposure thus eliminates motion blur). It is still beneficial because the benefit from the better conditioning outweighs the cost in energy. This documents is focused on light efficient flutter shutter that provides better conditioning and better energy utilization than conventional flutter shutter.


title: Optimal searching of gapped repeats in a word

abstract: Following (Kolpakov et al., 2013; Gawrychowski and Manea, 2015), we continue the study of {\em $\alpha$-gapped repeats} in strings, defined as factors $uvu$ with $|uv|\leq \alpha |u|$. Our main result is the $O(\alpha n)$ bound on the number of {\em maximal} $\alpha$-gapped repeats in a string of length $n$, previously proved to be $O(\alpha^2 n)$ in (Kolpakov et al., 2013). For a closely related notion of maximal $\delta$-subrepetition (maximal factors of exponent between $1+\delta$ and $2$), our result implies the $O(n/\delta)$ bound on their number, which improves the bound of (Kolpakov et al., 2010) by a $\log n$ factor. We also prove an algorithmic time bound $O(\alpha n+S)$ ($S$ size of the output) for computing all maximal $\alpha$-gapped repeats. Our solution, inspired by (Gawrychowski and Manea, 2015), is different from the recently published proof by (Tanimura et al., 2015) of the same bound. Together with our bound on $S$, this implies an $O(\alpha n)$-time algorithm for computing all maximal $\alpha$-gapped repeats.


title: Exploring the thermodynamics of non-commutative scalar fields

abstract: We study the thermodynamic properties of the Bose-Einstein condensate (BEC) in the context of the quantum field theory with non-commutative target space. Our main goal is to investigate in which temperature and/or energy regimes the non-commutativity can characterize some influence in the BEC properties described by a relativistic massive non-commutative boson gas. The non-commutativity parameters play a key role in the modified dispersion relations of the non-commutative fields, leading to a new phenomenology. We have obtained the condensate fraction, internal energy, pressure and specific heat of the system and taken ultra-relativistic (UR) and non-relativistic limits (NR). The non-commutative effects in the thermodynamic properties of the system are discussed. We found that there appear interesting signatures around the critical temperature.


title: A test of electric charge conservation with Borexino

abstract: Borexino is a liquid scintillation detector located deep underground at the Laboratori Nazionali del Gran Sasso (LNGS, Italy). Thanks to the unmatched radio-purity of the scintillator, and to the well understood detector response at low energy, a new limit on the stability of the electron for decay into a neutrino and a single mono-energetic photon was obtained. This new bound, tau > 6.6 10**28 yr at 90 % C.L., is two orders of magnitude better than the previous limit.


title: The role of multistability and transient trajectories in networked dynamical systems: Connectomic dynamics of C. elegans and behavioral assays

abstract: The neural dynamics of the nematode C. elegans are experimentally low-dimensional and correspond to discrete behavioral states, where previous modeling work has found neural proxies for some of these states. Experimental results further suggest that dynamics may be understood as long-timescale transitions between multiple low-dimensional attractors. To identify multistable regimes of our model, we develop a method for systematic generation of bifurcation diagrams and their analysis in an interpretable low-dimensional subspace, showing the existence and nature of multistable input responses at a glance. Stimulation of the PLM neuron pair, experimentally associated with forward movement and shown in simulation to drive a limit cycle, defines our low-dimensional projection space. We then obtain bifurcation diagrams for single-neuron excitation over a range of amplitudes and which classify whether the dynamics in this projection space are associated with a limit cycle, fixed point, or multiple states. In the specific case of compound input into both the PLM pair and ASK pair we discover bistability of a limit cycle and a fixed point, with transitional timescales between different states being much longer than other timescales in the system. This suggests consistency of our model with the characterization of dynamics in neural systems as long-timescale transitions between discrete, low-dimensional attractors corresponding to behavioral states. Our methodology thus prescribes a method for identifying these states and transitions in response to arbitrary input.


title: Evolution of cooperation in a particular case of the infinitely repeated Prisoner's Dilemma with three strategies

abstract: We will study a population of individuals playing the infinitely repeated Prisoner's Dilemma under replicator dynamics. The population consists of three kinds of individuals using the following reactive strategies: ALLD (individuals which always defect), ATFT (almost tit-for-tat: individuals which almost always repeat the opponent's last move) and G (generous individuals, which always cooperate when the opponent cooperated in the last move and have a positive probability $q$ of cooperating when they are defected). Our aim is studying in a mathematically rigorous fashion the dynamics of a simplified version for the computer experiment in [Nowak, Sigmund, Nature, 355, pp. 250--53, 1992] involving 100 reactive strategies. We will see that as the generosity degree of the G individuals varies, equilibria (rest points) of the dynamics appear or disappear, and the dynamics changes accordingly. Not only we will prove that the results of the experiment are true in our simplified version, but we will have complete control on the existence or non-existence of the equilbria for the dynamics for all possible values of the parameters, given that ATFT individuals are close enough to TFT. For most values of the parameters the dynamics will be completely determined.


title: Analog Computing Using Graphene-based Metalines

abstract: We introduce the new concept of "metalines" for manipulating the amplitude and phase profile of an incident wave locally and independently. Thanks to the highly confined graphene plasmons, a transmit-array of graphene-based metalines is used to realize analog computing on an ultra-compact, integrable and planar platform. By employing the general concepts of spatial Fourier transformation, a well-designed structure of such meta-transmit-array combined with graded index lenses can perform two mathematical operations; i.e. differentiation and integration, with high efficiency. The presented configuration is about 60 times shorter than the recent structure proposed by Silva et al.(Science, 2014, 343, 160-163); moreover, our simulated output responses are in more agreement with the desired analytic results. These findings may lead to remarkable achievements in light-based plasmonic signal processors at nanoscale instead of their bulky conventional dielectric lens-based counterparts.


title: Quantum Phase Imaging using Spatial Entanglement

abstract: Entangled photons have the remarkable ability to be more sensitive to signal and less sensitive to noise than classical light. Joint photons can sample an object collectively, resulting in faster phase accumulation and higher spatial resolution, while common components of noise can be subtracted. Even more, they can accomplish this while physically separate, due to the nonlocal properties of quantum mechanics. Indeed, nearly all quantum optics experiments rely on this separation, using individual point detectors that are scanned to measure coincidence counts and correlations. Scanning, however, is tedious, time consuming, and ill-suited for imaging. Moreover, the separation of beam paths adds complexity to the system while reducing the number of photons available for sampling, and the multiplicity of detectors does not scale well for greater numbers of photons and higher orders of entanglement. We bypass all of these problems here by directly imaging collinear photon pairs with an electron-multiplying CCD camera. We show explicitly the benefits of quantum nonlocality by engineering the spatial entanglement of the illuminating photons and introduce a new method of correlation measurement by converting time-domain coincidence counting into spatial-domain detection of selected pixels. We show that classical transport-of-intensity methods are applicable in the quantum domain and experimentally demonstrate nearly optimal (Heisenberg-limited) phase measurement for the given quantum illumination. The methods show the power of direct imaging and hold much potential for more general types of quantum information processing and control.


title: Machine Learning Model of the Swift/BAT Trigger Algorithm for Long GRB Population Studies

abstract: To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift/BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien 2014 is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of $\gtrsim97\%$ ($\lesssim 3\%$ error), which is a significant improvement on a cut in GRB flux which has an accuracy of $89.6\%$ ($10.4\%$ error). These models are then used to measure the detection efficiency of Swift as a function of redshift $z$, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of $n_0 \sim 0.48^{+0.41}_{-0.23} \ {\rm Gpc}^{-3} {\rm yr}^{-1}$ with power-law indices of $n_1 \sim 1.7^{+0.6}_{-0.5}$ and $n_2 \sim -5.9^{+5.7}_{-0.1}$ for GRBs above and below a break point of $z_1 \sim 6.8^{+2.8}_{-3.2}$. This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting. The code used in this is analysis is publicly available online (https://github.com/PBGraff/SwiftGRB_PEanalysis).


title: A Theory of Solving TAP Equations for Ising Models with General Invariant Random Matrices

abstract: We consider the problem of solving TAP mean field equations by iteration for Ising model with coupling matrices that are drawn at random from general invariant ensembles. We develop an analysis of iterative algorithms using a dynamical functional approach that in the thermodynamic limit yields an effective dynamics of a single variable trajectory. Our main novel contribution is the expression for the implicit memory term of the dynamics for general invariant ensembles. By subtracting these terms, that depend on magnetizations at previous time steps, the implicit memory terms cancel making the iteration dependent on a Gaussian distributed field only. The TAP magnetizations are stable fixed points if an AT stability criterion is fulfilled. We illustrate our method explicitly for coupling matrices drawn from the random orthogonal ensemble.


title: Kelvin Helmholtz instability in an ultrathin air film causes drop splashing on smooth surfaces

abstract: When a fast-moving drop impacts onto a smooth substrate, splashing will be produced at the edge of the expanding liquid sheet. This ubiquitous phenomenon lacks a fundamental understanding. Combining experiment with model, we illustrate that the ultrathin air film trapped under the expanding liquid front triggers splashing. Because this film is thinner than the mean free path of air molecules, the interior airflow transfers momentum with an unusually high velocity comparable to the speed of sound and generates a stress 10 times stronger than the airflow in common situations. Such a large stress initiates Kelvin Helmholtz instabilities at small length scales and effectively produces splashing. Our model agrees quantitatively with experimental verifications and brings a fundamental understanding to the ubiquitous phenomenon of drop splashing on smooth surfaces.


title: Constructing Dualities from Quantum State Manifolds

abstract: The thesis develops a systematic procedure to construct semi-classical gravitational duals from quantum state manifolds. Though the systems investigated are simple quantum mechanical systems without gauge symmetry many familiar concepts from the conventional gauge/gravity duality come about in a very natural way. The investigation of the low-dimensional manifolds link existing results in the $AdS_2/CFT_1$ literature. We are able to extend these in various ways and provide an explicit dictionary. The higher dimensional investigation is also concluded with a simple dictionary, but this dictionary requires the inclusion of many bulk coordinates. Consequently further work is needed to relate these results to existing literature. Possible ways to achieve this are discussed.


title: A major electronics upgrade for the H.E.S.S. Cherenkov telescopes 1-4

abstract: The High Energy Stereoscopic System (H.E.S.S.) is an array of imaging atmospheric Cherenkov telescopes (IACTs) located in the Khomas Highland in Namibia. It consists of four 12-m telescopes (CT1-4), which started operations in 2003, and a 28-m diameter one (CT5), which was brought online in 2012. It is the only IACT system featuring telescopes of different sizes, which provides sensitivity for gamma rays across a very wide energy range, from ~30 GeV up to ~100 TeV. Since the camera electronics of CT1-4 are much older than the one of CT5, an upgrade is being carried out; first deployment was in 2015, full operation is planned for 2016. The goals of this upgrade are threefold: reducing the dead time of the cameras, improving the overall performance of the array and reducing the system failure rate related to aging. Upon completion, the upgrade will assure the continuous operation of H.E.S.S. at its full sensitivity until and possibly beyond the advent of CTA. In the design of the new components, several CTA concepts and technologies were used and are thus being evaluated in the field: The upgraded read-out electronics is based on the NECTAR readout chips; the new camera front- and back-end control subsystems are based on an FPGA and an embedded ARM computer; the communication between subsystems is based on standard Ethernet technologies. These hardware solutions offer good performance, robustness and flexibility. The design of the new cameras is reported here.


title: Current fluctuations and large deviations for periodic TASEP on the relaxation scale

abstract: The one-dimensional totally asymmetric simple exclusion process (TASEP) with $N$ particles on a periodic lattice of $L$ sites is an interacting particle system with hopping rates breaking detailed balance. The total time-integrated current of particles $Q$ between time $0$ and time $T$ is studied for this model in the thermodynamic limit $L,N\to\infty$ with finite density of particles $\overline{\rho}=N/L$. The current $Q$ takes at leading order a deterministic value which follows from the hydrodynamic evolution of the macroscopic density profile by the inviscid Burgers' equation. Using asymptotics of Bethe ansatz formulas for eigenvalues and eigenvectors, an exact expression for the probability distribution of the fluctuations of $Q$ is derived on the relaxation time scale $T\sim L^{3/2}$ for an evolution conditioned on simple initial and final states. For flat initial and final states, a large deviation function expressed simply in terms of the Airy function is obtained at small rescaled time $T/L^{3/2}$.

categories: cond-mat.stat-mech math-ph math.MP nlin.SI


title: Topology On BCK-Modules

abstract: In this paper, we introduce the notion of a BCK-topological module in a natural way and establish that every decreasing sequence of submodules on a BCK-module M over bounded commutative BCK-algebra X is indeed a BCK- topological module. We have defined the notion of compatible and strict BCK- module homomorphisms, and establish that a strict BCK-module homomorphism is an open as well as a continuous mapping. Also, we establish the necessary and sufficient condition for a compatible mapping to be strict.


title: Macroscopic Einstein Equations for a Cosmological Model with $\lambda$-term

abstract: Through averaging the Einstein equations over transverse gravitational perturbations it is obtained a closed system of two ordinary differential equations describing macroscopic cosmological evolution of the isotropic space-flat Universe filled with gravitational radiation. It is found an asymptotic solution of evolution equation for gravitational perturbation amplitude. Making the substitution of this solution into Einstein equation averaged over gravitational perturbations, the single evolution non-linear ordinary differential second-order equation relative to macroscopic scale factor is obtained. It is also found a solution of evolution equation for scale factor in WKB-approximation which analytically describes the process of transformation from ultrarelativistic regime of cosmological extension to inflationary one.


title: New mapping properties of the Time Domain Electric Field Integral Equation

abstract: We show some improved mapping properties of the Time Domain Electric Field Integral Equation and of its Galerkin semidiscretization in space. We relate the weak distributional framework with a stronger class of solutions using a group of strongly continuous operators. The stability and error estimates we derive are sharper than those in the literature.


title: Existence of traversable wormholes in the spherical stellar systems

abstract: Potentiality of the presence of traversable wormholes in the outer/inner regions of the halos of galaxies, situated on the Navarro-Frenk-White (NFW) density profile and Universal Rotation Curve (URC) dark matter models have been investigated recently. Since this covers our own galaxy also as a possible home for traversable wormholes it prompts us to further the subject by considering alternative density distributions. From this token herein we make use of the Einasto model to describe the density profiles for the same purpose. Our choice for the latter is based on the fact that theoretical dark matter halos produced in computer simulations are best described by such a profile. For technical reasons we trim the number of parameters in the Einasto profile to a possible minimum. Based on such a model it is shown that traversable wormholes in the outer regions of spiral galaxies are possible while the inner part regions prohibit such formations.


title: Characterization of transiting exoplanets by way of differential photometry

abstract: This paper describes a simple activity for plotting and characterizing the light curve from an exoplanet transit event by way of differential photometry analysis. Using free digital imaging software, participants analyse a series of telescope images with the goal of calculating various exoplanet parameters, including its size, orbital radius and habitability. The activity has been designed for a high school or undergraduate university level and introduces fundamental concepts in astrophysics and an understanding of the basis for exoplanetary science, the transit method and digital photometry.


title: Universal set of Dynamically Protected Gates for Bipartite Qubit Networks II: Soft Pulse Implementation of the [[5,1,3]] Quantum Error Correcting Code

abstract: We model repetitive quantum error correction (QEC) with the single-error-correcting five-qubit code on a network of individually-controlled qubits with always-on Ising couplings, using our previously designed universal set of quantum gates based on sequences of shaped decoupling pulses. In addition to serving as accurate quantum gates, the sequences also provide dynamical decoupling (DD) of low-frequency phase noise. The simulation involves integrating unitary dynamics of six qubits over the duration of tens of thousands of control pulses, using classical stochastic phase noise as a source of decoherence. The combined DD/QEC protocol dramatically improves the coherence, with the QEC alone responsible for more than an order of magnitude infidelity reduction.


title: Train faster, generalize better: Stability of stochastic gradient descent

abstract: We show that parametric models trained by a stochastic gradient method (SGM) with few iterations have vanishing generalization error. We prove our results by arguing that SGM is algorithmically stable in the sense of Bousquet and Elisseeff. Our analysis only employs elementary tools from convex and continuous optimization. We derive stability bounds for both convex and non-convex optimization under standard Lipschitz and smoothness assumptions. Applying our results to the convex case, we provide new insights for why multiple epochs of stochastic gradient methods generalize well in practice. In the non-convex case, we give a new interpretation of common practices in neural networks, and formally show that popular techniques for training large deep models are indeed stability-promoting. Our findings conceptually underscore the importance of reducing training time beyond its obvious benefit.


title: Factorization of Temperley--Lieb diagrams

abstract: The Temperley--Lieb algebra is a finite dimensional associative algebra that arose in the context of statistical mechanics and occurs naturally as a quotient of the Hecke algebra arising from a Coxeter group of type $A$. It is often realized in terms of a certain diagram algebra, where every diagram can be written as a product of "simple diagrams." These factorizations correspond precisely to factorizations of the so-called fully commutative elements of the Coxeter group that index a particular basis. Given a reduced factorization of a fully commutative element, it is straightforward to construct the corresponding diagram. On the other hand, it is generally difficult to reconstruct the factorization given an arbitrary diagram. We present an efficient algorithm for obtaining a reduced factorization for a given diagram.


title: Transport coefficients of the Gribov-Zwanziger plasma

abstract: We study dynamic features of a plasma consisting of gluons whose infrared dynamics is improved by the Gribov-Zwanziger quantization. This approach embodies essential features of color confinement which set the plasma apart from conventional quasiparticle systems in several aspects. Our study focusses on a boost-invariant expansion for in- and out-of-equilibrium settings, which at late times can be characterized by the sound velocity, $c_s$, and the shear, $\eta$, and bulk, $\zeta$, viscosities. We obtain explicit expressions for the transport coefficients $\eta$ and $\zeta$ and check that they are consistent with the numerical solutions of the kinetic equation. At high temperature, keeping both the Gribov parameter and the relaxation time constant, we find a scaling $\zeta/\eta \propto 1/3 - c_s^2$ which manifests strong breaking of conformal symmetry in contrast to the case of weakly coupled plasmas.


title: General Relativity

abstract: This is an English translation of the Italian version of an encyclopedia chapter that appeared in the Italian Encyclopedia of the Physical Sciences, edited by Bruno Bertotti (1994). Following requests from colleagues we have decided to make it available to a more general readership. We present the motivation for constructing General Relativity, provide a short discussion of tensor algebra, and follow the set up of Einstein equations. We discuss briefly the initial value problem, the linear approximation and how should non gravitational physics be described in curved spacetime.


title: Finding matrix product state representations of highly-excited eigenstates of many-body localized Hamiltonians

abstract: A key property of many-body localization, the localization of quantum particles in systems with both quenched disorder and interactions, is the area law entanglement of even highly excited eigenstates of many-body localized Hamiltonians. Matrix Product States (MPS) can be used to efficiently represent low entanglement (area law) wave functions in one dimension. An important application of MPS is the widely used Density Matrix Renormalization Group (DMRG) algorithm for finding ground states of one dimensional Hamiltonians. Here, we describe two algorithms, the Shift and Invert MPS (SIMPS) and excited state DMRG which finds highly-excited eigenstates of many-body localized Hamiltonians. Excited state DMRG uses a modified sweeping procedure to identify eigenstates whereas SIMPS is a shift-inverse procedure that applies the inverse of the shifted Hamiltonian to a MPS multiple times to project out the targeted eigenstate. To demonstrate the power of these methods we verify the breakdown of the Eigenstate Thermalization Hypothesis (ETH) in the many-body localized phase of the random field Heisenberg model, show the saturation of entanglement in the MBL phase and generate local excitations.


title: Chiral Heat Wave and mixing of Magnetic, Vortical and Heat waves in chiral media

abstract: We show that a hot rotating fluid of relativistic chiral fermions possesses a new gapless collective mode associated with coherent propagation of energy density and chiral density waves along the axis of rotation. This mode, which we call the Chiral Heat Wave, emerges due to a mixed gauge-gravitational anomaly. At finite density the Chiral Heat Wave couples to the Chiral Vortical Wave while in the presence of an external magnetic field it mixes with the Chiral Magnetic Wave. The coupling of the Chiral Magnetic and Chiral Vortical Waves is also demonstrated. We find that the coupled waves - which are coherent fluctuations of the vector, axial and energy currents - have generally different velocities compared to the velocities of the individual waves.


title: Landau levels of Majorana fermions in a spin liquid

abstract: Majorana fermions were originally proposed as elementary particles acting as their own antiparticles. In recent years, it has become clear that Majorana fermions can instead be realized in condensed-matter systems as emergent quasiparticles, a situation often accompanied by topological order. Here we propose a physical system which realizes Landau levels - highly degenerate single-particle states usually resulting from an orbital magnetic field acting on charged particles - for Majorana fermions. This is achieved in a variant of a quantum spin system due to Kitaev which is distorted by triaxial strain. This strained Kitaev model displays a spin-liquid phase with charge-neutral Majorana-fermion excitations whose spectrum corresponds to that of Landau levels, here arising from a tailored pseudo-magnetic field. We show that measuring the dynamic spin susceptibility reveals the Landau-level structure by a remarkable mechanism of probe-induced bound-state formation.


title: Non-local gravity in D-dimensions: Propagator, entropy and bouncing Cosmology

abstract: We present the graviton propagator for an infinite derivative, $D$-dimensional, non-local action, up to quadratic order in curvature around a Minkowski background, and discuss the conditions required for this class of gravity theory to be ghost-free. We then study the gravitational entropy for de-Sitter and Anti-de Sitter backgrounds, before comparing with a recently derived result for a Schwarzschild blackhole, generalised to arbitrary $D$-dimensions, whereby the entropy is given simply by the area law. A novel approach of decomposing the entropy into its $(r, t)$ and spherical components is adopted in order to illustrate the differences more clearly. We conclude with a discussion of de-Sitter entropy in the framework of a non-singular bouncing cosmology.


title: A Soft X-Ray Spectral Episode for the Clocked Burster, GS 1826-24 as Measured by Swift and NuSTAR

abstract: We report on NuSTAR and Swift observations of a soft state of the neutron star low-mass X-ray binary GS 1826-24, commonly known as the "clocked" burster. The transition to the soft state was recorded in 2014 June through an increase of the 2-20 keV source intensity measured by MAXI, simultaneous with a decrease of the 15-50 keV intensity measured by Swift/BAT. The episode lasted approximately two months, after which the source returned to its usual hard state. We analyze the broad-band spectrum measured by Swift/XRT and NuSTAR, and estimate the accretion rate during the soft episode to be about 13% of Eddington, within the range of previous observations. However, the best fit spectral model, adopting the double Comptonization used previously, exhibits significantly softer components. We detect seven type-I X-ray bursts, all significantly weaker (and with shorter rise and decay times) than observed previously. The burst profiles and recurrence times vary significantly, ruling out the regular bursts that are typical for this source. One burst exhibited photospheric radius expansion, and we estimate the source distance at about (5.7 / xi_b^1/2) kpc, where xi_b parameterizes the possible anisotropy of the burst emission. Interpreting the soft state as a transition from an optically thin inner flow to an optically thick flow passing through a boundary layer, as is commonly observed in similar systems, is contradicted by the lower optical depth measured for the double Comptonization model we find for this soft state. The effect of a change in disk geometry on the burst behavior remains unclear.


title: Lepton Flavor Non-Universality in B-meson Decays from a U(2) Flavor Model

abstract: We address the recent anomalies in semi-leptonic $B$-meson decays using a model of fermion masses based on the $U(2)$ flavor symmetry. The new contributions to $b \to s \ell \ell$ transitions arise due to a tree-level exchange of a $Z^\prime$ vector boson gauging a $U(1)$ subgroup of the flavor symmetry. They are controlled by a single parameter and are approximately aligned to the Standard Model prediction, with constructive interference in the $e$-channel and destructive interference in the $\mu$-channel. The current experimental data on semi-leptonic $B$-meson decays can be very well reproduced without violating existing constraints from flavor violation in the quark and lepton sectors. Our model will be tested by new measurements of $b \to s \ell \ell$ transitions and also by future electroweak precision tests, direct $Z^\prime$ searches, and $\mu$-$e$ conversion in nuclei.


title: Push it to the limit: Local Group constraints on high-redshift stellar mass functions for Mstar > 10^5 Msun

abstract: We constrain the evolution of the galaxy stellar mass function from 2 < z < 5 for galaxies with stellar masses as low as 10^5 Msun by combining star formation histories of Milky Way satellite galaxies derived from deep Hubble Space Telescope observations with merger trees from the ELVIS suite of N-body simulations. This approach extends our understanding more than two orders of magnitude lower in stellar mass than is currently possible by direct imaging. We find the faint end slopes of the mass functions to be alpha= -1.42(+0.07/-0.05) at z = 2 and alpha = -1.57^(+0.06/-0.06) at z = 5, and show the slope only weakly evolves from z = 5 to z = 0. Our findings are in stark contrast to a number of direct detection studies that suggest slopes as steep as alpha = -1.9 at these epochs. Such a steep slope would result in an order of magnitude too many luminous Milky Way satellites in a mass regime that is observationally complete (Mstar > 2*10^5 Msun at z = 0). The most recent studies from ZFOURGE and CANDELS also suggest flatter faint end slopes that are consistent with our results, but with a lower degree of precision. This work illustrates the strong connections between low and high-z observations when viewed through the lens of LCDM numerical simulations.


title: Spin-mediated dissipation and frequency shifts of a cantilever at milliKelvin temperatures

abstract: We measure the dissipation and frequency shift of a magnetically coupled cantilever in the vicinity of a silicon chip, down to $25$ mK. The dissipation and frequency shift originates from the interaction with the unpaired electrons, associated with the dangling bonds in the native oxide layer of the silicon, which form a two dimensional system of electron spins. We approach the sample with a $3.43$ $\mu$m-diameter magnetic particle attached to an ultrasoft cantilever, and measure the frequency shift and quality factor as a function of temperature and the distance. Using a recent theoretical analysis [J. M. de Voogd et al., arXiv:1508.07972 (2015)] of the dynamics of a system consisting of a spin and a magnetic resonator, we are able to fit the data and extract the relaxation time $T_1=0.39\pm0.08$ ms and spin density $\sigma=0.14\pm0.01$ spins per nm$^2$. Our analysis shows that at temperatures $\leq500$ mK magnetic dissipation is an important source of non-contact friction.


title: Chemistry of the Most Metal-poor Stars in the Bulge and the z > 10 Universe

abstract: Metal-poor stars in the Milky Way are local relics of the epoch of the first stars and the first galaxies. However, a low metallicity does not prove that a star formed in this ancient era, as metal-poor stars form over a range of redshift in different environments. Theoretical models of Milky Way formation have shown that at constant metallicity, the oldest stars are those closest to the center of the Galaxy on the most tightly-bound orbits. For that reason, the most metal-poor stars in the bulge of the Milky Way provide excellent tracers of the chemistry of the high-redshift universe. We report the dynamics and detailed chemical abundances of three stars in the bulge with [Fe/H] $\lesssim-2.7$, two of which are the most metal-poor stars in the bulge in the literature. We find that with the exception of scandium, all three stars follow the abundance trends identified previously for metal-poor halo stars. These three stars have the lowest [Sc II/Fe] abundances yet seen in $\alpha$-enhanced giant stars in the Galaxy. Moreover, all three stars are outliers in the otherwise tight [Sc II/Fe]-[Ti II/Fe] relation observed among metal-poor halo stars. Theoretical models predict that there is a 30% chance that at least one of these stars formed at $z\gtrsim15$, while there is a 70% chance that at least one formed at $10 \lesssim z \lesssim 15$. These observations imply that by $z\sim10$, the progenitor galaxies of the Milky Way had both reached [Fe/H] $\sim-3.0$ and established the abundance pattern observed in extremely metal-poor stars.


title: The IMF-sensitive 1.14-micron Na I doublet in early-type galaxies

abstract: We present J-band spectroscopy of passive galaxies focusing on the Na I doublet at 1.14 {\mu}m. Like the Na I 0.82 {\mu}m doublet, this feature is strong in low-mass stars and hence may provide a useful probe of the initial mass function (IMF). From high signal-to-noise composite spectra, we find that Na I 1.14 {\mu}m increases steeply with increasing velocity dispersion, {\sigma}, and for the most massive galaxies (\sigma > 300 km/s) is much stronger than predicted from synthetic spectra with Milky-Way-like IMFs and solar abundances. Reproducing Na I 1.14 {\mu}m at high {\sigma} likely requires either a very high [Na/H], or a bottom-heavy IMF, or a combination of both. Using the Na D line to break the degeneracy between IMF and abundance, we infer [Na/H] $\approx$ +0.5 and a steep IMF (single-slope-equivalent x $\approx$ 3.2, where x = 2.35 for Salpeter), for the high-\sigma galaxies. At lower mass ({\sigma} = 50-100 km/s), the line strengths are compatible with MW-like IMFs and near-solar [Na/H]. We highlight two galaxies in our sample where strong gravitational lensing masses favour MW-like IMFs. Like the high-{\sigma} sample on average, these galaxies have strong Na I 1.14 \mu m; taken in isolation their sodium indices imply bottom-heavy IMFs which are hard to reconcile with the lensing masses. An alternative full-spectrum-fitting approach, applied to the high-\sigma sample, recovers an IMF less heavy than Salpeter, but under-predicts the Na I 1.14 \mu m line at the 5{\sigma} level. We conclude that current models struggle to reproduce this feature in the most massive galaxies without breaking other constraints, and caution against over-reliance on the sodium lines in spectroscopic IMF studies.


title: The environments of high redshift radio galaxies and quasars: probes of protoclusters

abstract: We use the GALFORM semi-analytical model to study high density regions traced by radio galaxies and quasars at high redshifts. We explore the impact that baryonic physics has upon the properties of galaxies in these environments. Star-forming emission-line galaxies (Ly{\alpha} and H{\alpha} emitters) are used to probe the environments at high redshifts. Radio galaxies are predicted to be hosted by more massive haloes than quasars, and this is imprinted on the amplitude of galaxy overdensities and cross-correlation functions. We find that Ly{\alpha} radiative transfer and AGN feedback indirectly affect the clustering on small scales and also the stellar masses, star- formation rates and gas metallicities of galaxies in dense environments. We also investigate the relation between protoclusters associated with radio galaxies and quasars, and their present- day cluster descendants. The progenitors of massive clusters associated with radio galaxies and quasars allow us to determine an average protocluster size in a simple way. Overdensities within the protoclusters are found to correlate with the halo descendant masses. We present scaling relations that can be applied to observational data. By computing projection effects due to the wavelength resolution of modern spectrographs and narrow-band filters we show that the former have enough spectral resolution to map the structure of protoclusters, whereas the latter can be used to measure the clustering around radio galaxies and quasars over larger scales to determine the mass of dark matter haloes hosting them.


title: The Caterpillar Project: A Large Suite of Milky Way Sized Halos

abstract: We present the largest number of Milky Way sized dark matter halos simulated at very high mass ($\sim$$10^4$ M$_\odot$/particle) and temporal resolution ($\sim$5 Myrs/snapshot) done to date, quadrupling what is currently available in the literature. This initial suite consists of the first 24 halos of the $Caterpillar$ $Project$ (www.caterpillarproject.org) whose project goal of 60 - 70 halos will be made public when complete. We resolve $\sim$20,000 gravitationally bound subhalos within the virial radius of each host halo. Over the ranges set by our spatial resolution our convergence is excellent and improvements were made upon current state-of-the-art halo finders to better identify substructure at such high resolutions (e.g., on average we recover $\sim$4 subhalos in each host halo above 10$^8$ M$_\odot$ which would have otherwise not been found using conventional methods). For our relaxed halos, the inner profiles are reasonably fit by Einasto profiles ($\alpha$ = 0.169 $\pm$ 0.023) though this depends on the relaxed nature and assembly history of a given halo. Averaging over all halos, the substructure mass fraction is $f_{m,subs} = 0.121 \pm 0.041$, and mass function slope is d$N$/d$M\propto M^{-1.88 \pm 0.10}$ though we find scatter in the normalizations for fixed halo mass due to more concentrated hosts having less subhalos at fixed subhalo mass. There are no biases stemming from Lagrangian volume selection as all Lagrangian volume types are included in our sample. Our detailed contamination study of 264 low resolution halos has resulted in obtaining very large and unprecedented, high-resolution regions around our host halos for our target resolution (sphere of radius $\sim$$1.4 \pm 0.4$ Mpc) allowing for accurate studies of low mass dwarf galaxies at large galactocentric radii and the very first stellar systems at high redshift ($z \geq$ 10).


title: Permutation Orbifolds in the large N Limit

abstract: The space of permutation orbifolds is a simple landscape of two dimensional CFTs, generalizing the well-known symmetric orbifolds. We consider constraints which a permutation orbifold with large central charge must obey in order to be holographically dual to a weakly coupled (but possibly stringy) theory of gravity in AdS. We then construct explicit examples of permutation orbifolds which obey these constraints. In our constructions the spectrum remains finite at large N, but differs qualitatively from that of symmetric orbifolds. We also discuss under what conditions the correlation functions factorize at large N and thus reduce to those of a generalized free field in AdS. We show that this happens not just for symmetric orbifolds, but also for permutation groups which act "democratically" in a sense which we define.


title: The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

abstract: We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.


title: Special Quasirandom Structures: a selection approach for stochastic homogenization

abstract: We adapt and study a variance reduction approach for the homogenization of elliptic equations in divergence form. The approach, borrowed from atomistic simulations and solid-state science [von Pezold et al, Physical Review B 2010; Wei et al, Physical Review B 1990; Zunger et al, Physical Review Letters 1990], consists in selecting random realizations that best satisfy some statistical properties (such as the volume fraction of each phase in a composite material) usually only obtained asymptotically. We study the approach theoretically in some simplified settings (one-dimensional setting, perturbative setting in higher dimensions), and numerically demonstrate its efficiency in more general cases.


title: Does the obscured AGN fraction really depend on luminosity?

abstract: We use a sample of 151 local non-blazar AGN selected from the INTEGRAL all-sky hard X-ray survey to investigate if the observed declining trend of the fraction of obscured (i.e. showing X-ray absorption) AGN with increasing luminosity is mostly an intrinsic or selection effect. Using a torus-obscuration model, we demonstrate that in addition to negative bias, due to absorption in the torus, in finding obscured AGN in hard X-ray flux limited surveys, there is also positive bias in finding unobscured AGN, due to Compton reflection in the torus. These biases can be even stronger taking into account plausible intrinsic collimation of hard X-ray emission along the axis of the obscuring torus. Given the AGN luminosity function, which steepens at high luminosities, these observational biases lead to a decreasing observed fraction of obscured AGN with increasing luminosity even if this fraction has no intrinsic luminosity dependence. We find that if the central hard X-ray source in AGN is isotropic, the intrinsic (i.e. corrected for biases) obscured AGN fraction still shows a declining trend with luminosity, although the intrinsic obscured fraction is significantly larger than the observed one: the actual fraction is larger than $\sim 85$% at $L\lesssim 10^{42.5}$ erg/s (17--60 keV), and decreases to $\lesssim 60$% at $L\gtrsim 10^{44}$ erg/s. In terms of the half-opening angle, $\theta$, of an obscuring torus, this implies that $\theta\lesssim 30$ deg in lower-luminosity AGN, and $\theta\gtrsim 45$ deg in higher-luminosity ones. If, however, the emission from the central SMBH is collimated as $dL/d\Omega\propto\cos\alpha$, the intrinsic dependence of the obscured AGN fraction is consistent with a luminosity-independent torus half-opening angle $\theta\sim 30$ deg.


title: Variational wave function for an anisotropic single-hole-doped $t$-$J$ ladder

abstract: Based on three general guiding principles, i.e., no double occupancy constraint, accurate description of antiferromagnetism at half-filling, and the precise sign structure of the $t$-$J$ model, a new ground state wave function has been constructed recently [Weng, New J. Phys. 13, 103039 (2011)]. In this paper, we specifically study such kind of variational ground state for the one-hole-doped anisotropic two-leg $t$-$J$ ladder using variational Monte Carlo (VMC) method. The results are then systematically compared with those recently obtained by density matrix renormalization group (DMRG) simulation. An excellent agreement is found between the VMC and DMRG results, including a "quantum critical point" at the anisotropy parameter $\alpha=\alpha_c\approx0.7$ (with the parameters $t/J=3$), and the emergence of charge modulation and momentum (Fermi point) reconstruction at $\alpha>\alpha_c$ due to the quantum interference of the sign structure. In particular, the wave function indicates that a Landau's quasiparticle description remains valid at $\alpha<\alpha_c$ but fails at $\alpha>\alpha_c$ due to the breakdown of the one-to-one correspondence of momentum and translational symmetry of the hole. The explicit form of the wave function provides a direct understanding on how the many-body strong correlation effect takes place non-perturbatively in a doped Mott insulator, which sheds interesting light on the two-dimensional case where the same type of wave function was proposed to describe the cuprate superconductor.


title: The spatial clustering of ROSAT All-Sky Survey Active Galactic Nuclei IV. More massive black holes reside in more massive dark matter halos

abstract: This is the fourth paper in a series that reports on our investigation of the clustering properties of active galactic nuclei (AGN) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this paper we investigate the cause of the X-ray luminosity dependence of the clustering of broad-line, luminous AGN at 0.16<z<0.36. We fit the H-alpha line profile in the SDSS spectra for all X-ray and optically-selected broad-line AGN, determine the mass of the super-massive black hole (SMBH), M_BH, and infer the accretion rate relative to Eddington (L/L_EDD). Since M_BH and L/L_EDD are correlated, we create AGN subsamples in one parameter while maintaining the same distribution in the other parameter. In both the X-ray and optically-selected AGN samples we detect a weak clustering dependence with M_BH and no statistically significant dependence on L/L_EDD. We find a difference of up to 2.7sigma when comparing the objects that belong to the 30% least and 30% most massive M_BH subsamples, in that luminous broad-line AGN with more massive black holes reside in more massive parent dark matter halos at these redshifts. These results provide evidence that higher accretion rates in AGN do not necessarily require dense galaxy environments in which more galaxy mergers and interactions are expected to channel large amounts of gas onto the SMBH. We also present semi-analytic models which predict a positive M_DMH dependence on M_BH, which is most prominent at M_BH ~ 10^{8-9} M_SUN.


title: Spectral correlation in down-converted photon pairs at telecom wavelength

abstract: Sources of photon pairs based on the spontaneous parametric down conversion process are commonly used for long distance quantum communication. The key feature for improving the range of transmission is engineering their spectral properties. Following two experimental papers [Opt. Lett., 38, 697 (2013)] and [Opt. Lett., 39, 1481 (2014)] we analytically and numerically analyze the characteristics of a source. It is based on a $\beta$-barium borate (BBO) crystal cut for type II phase matching at the degenerated frequencies 755 nm $\rightarrow$ 1550 nm + 1550 nm. Our analysis shows a way for full control of spectral correlation within a fiber-coupled photon pair simultaneously with optimal brightness.


title: Gravitational Encounters and the Evolution of Galactic Nuclei. III. Anomalous Relaxation

abstract: This paper is the third in a series presenting the results of direct numerical integrations of the Fokker-Planck equation for stars orbiting a supermassive black hole (SBH) at the center of a galaxy. The algorithm of Paper II included diffusion coefficients that described the effects of random ("classical") and correlated ("resonant") relaxation. In this paper, the diffusion coefficients of Paper II have been generalized to account for the effects of "anomalous relaxation," the qualitatively different way in which eccentric orbits evolve in the regime of rapid relativistic precession. Two functional forms for the anomalous diffusion coefficients are investigated, based on power-law or exponential modifications of the resonant diffusion coefficients. The parameters defining the modified coefficients are first constrained by comparing the results of Fokker-Planck integrations with previously-published N-body integrations. Steady-state solutions are then obtained via the Fokker-Planck equation for models with properties similar to those of the Milky Way nucleus. Inclusion of anomalous relaxation leads to the formation of less prominent cores than in the case of resonant relaxation alone, due to the lengthening of diffusion timescales for eccentric orbits. Steady-state capture rates of stars by the SBH are found to always be less, by as much as an order of magnitude, than capture rates in the presence of resonant relaxation alone.


title: Long-lived states with well-defined spins in spin-$1/2$ homogeneous Bose gases

abstract: Many-body eigenfunctions of the total spin operator can be constructed from the spin and spatial wavefunctions with non-trivial permutation symmetries. Spin-dependent interactions can lead to relaxation of the spin eigenstates to the thermal equilibrium. A mechanism that stabilizes the many-body entangled states is proposed here. Surprisingly, in spite coupling with the chaotic motion of the spatial degrees of freedom, the spin relaxation can be suppressed by destructive quantum interference due to spherical vector and tensor terms of the spin-dependent interactions. Tuning the scattering lengths by the method of Feshbach resonances, readily available in cold atomic labs, can enhance the relaxation timescales by several orders of magnitude.


title: On entropy production in the Madelung fluid and the role of Bohm's potential in classical diffusion

abstract: The Madelung equations map the non-relativistic time-dependent Schrodinger equation into hydrodynamic equations of a virtual fluid. Here we show that an increase of the Boltzmann entropy of this Madelung fluid is proportional to the expectation value of its velocity divergence. Hence, entropy growth is accompanied by expansion resulting from the ability of the Madelung fluid to be compressible. The compressibility itself reflects superposition of solutions of the Schrodinger equation. Thus, in unitary processes where the Madelung fluid expands and then shrinks, the Boltzmann entropy may, correspondingly, grow and then decrease. The notion of entropy growth due to expansion is common in diffusive processes, however in the latter the process is irreversible. Much unlike the Boltzmann entropy, the von Neumann entropy, does not vary with time. To elucidate the physical underpinning of the Boltzmann entropy, we examine several specific examples. We demonstrate that, for classical diffusive processes, the "force" accelerating diffusion has the form of the positive gradient of the quantum Bohm potential. In the Madelung fluid, the advective and the diffusive velocities correspond respectively to the the real and imaginary parts of the complex momentum. We find that the diffusion coefficient provides a lower bound of Heisenberg uncertainty type product between the gas mean free path and the Brownian momentum.


title: All-Optical Vector Measurement of Spin-Orbit-Induced Torques Using Both Polar and Quadratic Magneto-Optic Kerr Effects

abstract: We demonstrate that the magneto-optic-Kerr effect with normal light incidence can be used to obtain quantitative optical measurements of both components of spin-orbit-induced torque (both the antidamping and effective-field components) in heavy-metal/ferromagnet bilayers. This is achieved by analyzing the quadratic Kerr effect as well as the polar Kerr effect. The two effects can be distinguished by properly selecting the polarization of the incident light. We use this all-optical technique to determine the spin-orbit torques generated by a series of Pt/Permalloy samples, finding values in excellent agreement with spin-torque ferromagnetic resonance measurements.


title: On positive viscosity solutions of fractional Lane-Emden systems

abstract: We discuss on existence, nonexistence and uniqueness of positive viscosity solutions for Lane-Emden systems involving the fractional Laplacian on bounded domains. As a byproduct, we obtain the critical hyperbole associated to the these systems.


title: On the Concept of Cryptographic Quantum Hashing

abstract: In the paper we define a notion of quantum resistant ($(\epsilon,\delta)$-resistant) hash function which combine together a notion of pre-image (one-way) resistance ($\epsilon$-resistance) property we define in the paper and the notion of collision resistance ($\delta$-resistance) properties. We show that in the quantum setting a one-way resistance property and collision resistance property are correlated: the "more" a quantum function is one-way resistant the "less" it collision resistant and vice versa. We present an explicit quantum hash function which is "balanced" one-way resistant and collision resistant and demonstrate how to build a large family quantum hash functions. Balanced quantum hash functions need a high degree of entanglement between the qubits. We use a "phase constructions" technique to express quantum hashing constructions, which is good to map hash states to coherent states in a superposition of time-bin modes. The later is ready to be implemented with current optical technology.


title: Nanomechanical detection of the spin Hall effect

abstract: The spin Hall effect creates a spin current in response to a charge current in a material that has strong spin-orbit coupling. The size of the spin Hall effect in many materials is disputed, requiring independent measurements of the effect. We develop a novel mechanical method to measure the size of the spin Hall effect, relying on the equivalence between spin and angular momentum. The spin current carries angular momentum, so the flow of angular momentum will result in a mechanical torque on the material. We determine the size and geometry of this torque and demonstrate that it can be measured using a nanomechanical device. Our results show that measurement of the spin Hall effect in this manner is possible and also opens possibilities for actuating nanomechanical systems with spin currents.


title: Machine Learning Methods to Analyze Arabidopsis Thaliana Plant Root Growth

abstract: One of the challenging problems in biology is to classify plants based on their reaction on genetic mutation. Arabidopsis Thaliana is a plant that is so interesting, because its genetic structure has some similarities with that of human beings. Biologists classify the type of this plant to mutated and not mutated (wild) types. Phenotypic analysis of these types is a time-consuming and costly effort by individuals. In this paper, we propose a modified feature extraction step by using velocity and acceleration of root growth. In the second step, for plant classification, we employed different Support Vector Machine (SVM) kernels and two hybrid systems of neural networks. Gated Negative Correlation Learning (GNCL) and Mixture of Negatively Correlated Experts (MNCE) are two ensemble methods based on complementary feature of classical classifiers; Mixture of Expert (ME) and Negative Correlation Learning (NCL). The hybrid systems conserve of advantages and decrease the effects of disadvantages of NCL and ME. Our Experimental shows that MNCE and GNCL improve the efficiency of classical classifiers, however, some SVM kernels function has better performance than classifiers based on neural network ensemble method. Moreover, kernels consume less time to obtain a classification rate.


title: Probabilistic Neural Network Training for Semi-Supervised Classifiers

abstract: In this paper, we propose another version of help-training approach by employing a Probabilistic Neural Network (PNN) that improves the performance of the main discriminative classifier in the semi-supervised strategy. We introduce the PNN-training algorithm and use it for training the support vector machine (SVM) with a few numbers of labeled data and a large number of unlabeled data. We try to find the best labels for unlabeled data and then use SVM to enhance the classification rate. We test our method on two famous benchmarks and show the efficiency of our method in comparison with pervious methods.


title: Noise-assisted energy transport in electrical oscillator networks with off-diagonal dynamical disorder

abstract: Noise is generally thought as detrimental for energy transport in coupled oscillator networks. However, it has been shown that for certain coherently evolving systems, the presence of noise can enhance, somehow unexpectedly, their transport efficiency; a phenomenon called environment-assisted quantum transport (ENAQT) or dephasing-assisted transport. Here, we report on the experimental observation of such effect in a network of coupled electrical oscillators. We demonstrate that by introducing stochastic fluctuations in one of the couplings of the network, a relative enhancement in the energy transport efficiency of $22.5 \pm 3.6\,\%$ can be observed.


title: Subshifts with Slowly Growing Numbers of Follower Sets

abstract: For any subshift, define $F_X(n)$ to be the collection of distinct follower sets of words of length $n$ in $X$. Based on a similar result of the second and third authors, we conjecture that if there exists an $n$ for which $|F_X(n)| \leq n$, then $X$ is sofic. In this paper, we prove several results related to this conjecture, including verifying it for $n \leq 3$, proving that the conjecture is true for a large class of coded subshifts, and showing that if there exists $n$ for which $|F_X(n)| \leq \log_2(n+1)$, then $X$ is sofic.


title: Plasma Viscosity with Mass Transport in Spherical ICF Implosion Simulations

abstract: The effects of viscosity and small-scale atomic-level mixing on plasmas in inertial confinement fusion (ICF) currently represent challenges in ICF research. Many current ICF hydrodynamic codes ignore the effects of viscosity though recent research indicates viscosity and mixing by classical transport processes may have a substantial impact on implosion dynamics. We have implemented a Lagrange hydrodynamic code in one-dimensional spherical geometry with plasma viscosity and mass transport and including a three temperature model for ions, electrons, and radiation treated in a gray radiation diffusion approximation. The code is used to study ICF implosion differences with and without plasma viscosity and to determine the impacts of viscosity on temperature histories and neutron yield. It was found that plasma viscosity has substantial impacts on ICF shock dynamics characterized by shock burn timing, maximum burn temperatures, convergence ratio, and time history of neutron production rates. Plasma viscosity reduces the need for artificial viscosity to maintain numerical stability in the Lagrange formulation and also modifies the flux-limiting needed for electron thermal conduction.


title: Modeling Long-term Outcomes and Treatment Effects After Androgen Deprivation Therapy for Prostate Cancer

abstract: Analyzing outcomes in long-term cancer survivor studies can be complex. The effects of predictors on the failure process may be difficult to assess over longer periods of time, as the commonly used assumption of proportionality of hazards holding over an extended period is often questionable. In this manuscript, we compare seven different survival models that estimate the hazard rate and the effects of proportional and non-proportional covariates. In particular, we focus on an extension of the the multi-resolution hazard (MRH) estimator, combining a non-proportional hierarchical MRH approach with a data-driven pruning algorithm that allows for computational efficiency and produces robust estimates even in times of few observed failures. Using data from a large-scale randomized prostate cancer clinical trial, we examine patterns of biochemical failure and estimate the time-varying effects of androgen deprivation therapy treatment and other covariates. We compare the impact of different modeling strategies and smoothness assumptions on the estimated treatment effect. Our results show that the benefits of treatment diminish over time, possibly with implications for future treatment protocols.


title: The Stripe 82 Massive Galaxy Project I: Catalog Construction

abstract: The Stripe 82 Massive Galaxy Catalog (S82-MGC) is the largest-volume stellar mass-limited sample of galaxies beyond z~1 constructed to date. Spanning 139.4 deg2, the S82-MGC includes a mass-limited sample of 41,770 galaxies with log Mstar > 11.2 to z~0.7, sampling a volume of 0.3 Gpc3, roughly equivalent to the volume of the Sloan Digital Sky Survey-I/II (SDSS-I/II) z < 0.15 MAIN sample. The catalog is built on three pillars of survey data: the SDSS Stripe 82 Coadd photometry which reaches r-band magnitudes of 23.5 AB, YJHK photometry at depths of 20th magnitude (AB) from the UK Infrared Deep Sky Survey (UKIDSS) Large Area Survey, and over 70,000 spectroscopic galaxy redshifts from SDSS-I/II and the Baryon Oscillation Spectroscopic Survey (BOSS). We describe the catalog construction and verification, the production of 9-band matched aperture photometry, tests of existing and newly estimated photometric redshifts required to supplement spectroscopic redshifts for 55% of the log Mstar > 11.2 sample, and geometric masking. We provide near-IR based stellar mass estimates and compare these to previous estimates. All catalog products are made publicly available. The S82-MGC not only addresses previous statistical limitations in high-mass galaxy evolution studies but begins tackling inherent data challenges in the coming era of wide-field imaging surveys.


title: A Dataset for Improved RGBD-based Object Detection and Pose Estimation for Warehouse Pick-and-Place

abstract: An important logistics application of robotics involves manipulators that pick-and-place objects placed in warehouse shelves. A critical aspect of this task corre- sponds to detecting the pose of a known object in the shelf using visual data. Solving this problem can be assisted by the use of an RGB-D sensor, which also provides depth information beyond visual data. Nevertheless, it remains a challenging problem since multiple issues need to be addressed, such as low illumination inside shelves, clutter, texture-less and reflective objects as well as the limitations of depth sensors. This paper provides a new rich data set for advancing the state-of-the-art in RGBD- based 3D object pose estimation, which is focused on the challenges that arise when solving warehouse pick- and-place tasks. The publicly available data set includes thousands of images and corresponding ground truth data for the objects used during the first Amazon Picking Challenge at different poses and clutter conditions. Each image is accompanied with ground truth information to assist in the evaluation of algorithms for object detection. To show the utility of the data set, a recent algorithm for RGBD-based pose estimation is evaluated in this paper. Based on the measured performance of the algorithm on the data set, various modifications and improvements are applied to increase the accuracy of detection. These steps can be easily applied to a variety of different methodologies for object pose detection and improve performance in the domain of warehouse pick-and-place.


title: Adiabatic and Hamiltonian computing on a 2D lattice with simple 2-qubit interactions

abstract: We show how to perform universal Hamiltonian and adiabatic computing using a time-independent Hamiltonian on a 2D grid describing a system of hopping particles which string together and interact to perform the computation. In this construction, the movement of one particle is controlled by the presence or absence of other particles, an effective quantum field effect transistor that allows the construction of controlled-NOT and controlled-rotation gates. The construction translates into a model for universal quantum computation with time-independent 2-qubit ZZ and XX+YY interactions on an (almost) planar grid. The effective Hamiltonian is arrived at by a single use of first-order perturbation theory avoiding the use of perturbation gadgets. The dynamics and spectral properties of the effective Hamiltonian can be fully determined as it corresponds to a particular realization of a mapping between a quantum circuit and a Hamiltonian called the space-time circuit-to-Hamiltonian construction. Because of the simple interactions required, and because no higher-order perturbation gadgets are employed, our construction is potentially realizable using superconducting or other solid-state qubits.


title: A Keck Survey of Gravitationally-Lensed Star-Forming Galaxies: High Spatial Resolution Studies of Kinematics and Metallicity Gradients

abstract: We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 gravitationally lensed star-forming galaxies at a mean redshift of $z\simeq2$ based on Keck laser-assisted adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (IFU) spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY survey, we sample these sub-L$^{\ast}$ galaxies with source-plane resolutions of a few hundred parsecs ensuring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high spatial resolution data offers a critical check on the structural properties of larger samples derived with coarser sampling using multiple-IFU instruments. We demonstrate how kinematic complexities essential to understanding the maturity of an early star-forming galaxy can often only be revealed with better sampled data. Although we include four sources from our earlier work, the present study provides a more representative sample unbiased with respect to emission line strength. Contrary to earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior inconsistent with a simple picture of well-ordered rotation developing concurrently with established steep metal gradients in all but merging systems. Comparing our observations with the predictions of hydrodynamical simulations suggests that gas and metals have been mixed by outflows or other strong feedback processes, flattening the metal gradients in early star-forming galaxies.


title: The Kozai-Lidov Mechanism in Hydrodynamical Disks. III. Effects of Disk Mass and Self-Gravity

abstract: Previously we showed that a substantially misaligned viscous accretion disk with pressure that orbits around one component of a binary system can undergo global damped Kozai-Lidov (KL) oscillations. These oscillations produce periodic exchanges of the disk eccentricity with inclination. The disk KL mechanism is quite robust and operates over a wide range of binary and disk parameters. However, the effects of self-gravity, which are expected to suppress the KL oscillations for sufficiently massive disks, were ignored. Here, we analyze the effects of disk self-gravity by means of hydrodynamic simulations and compare the results with the expectations of analytic theory. The disk mass required for suppression in the simulations is a few percent of the mass of the central star and this roughly agrees with an analytical estimate. The conditions for suppression of the KL oscillations in the simulations are close to requiring that the disk be gravitationally unstable. We discuss some implications of our results for the dynamics of protoplanetary disks and the related planet formation.


title: Beam engineering for selective and enhanced coupling to multipolar resonances

abstract: Multipolar electromagnetic phenomena in sub-wavelength resonators are at the heart of metamaterial science and technology. In this letter, we demonstrate selective and enhanced coupling to specific multipole resonances via beam engineering. We first derive an analytical method for determining the scattering and absorption of spherical nanoparticles (NPs) that depends only on the local electromagnetic field quantities within an inhomogeneous beam. Using this analytical technique, we demonstrate the ability to drastically manipulate the scattering properties of a spherical NP by varying illumination properties and demonstrate the excitation of a longitudinal quadrupole mode that cannot be accessed with conventional illumination. This work enhances the understanding of fundamental light-matter interactions in metamaterials, and lays the foundation for researchers to identify, quantify, and manipulate multipolar light-matter interactions through optical beam engineering.


title: Catenaries in viscous fluid

abstract: This work explores a simple model of a slender, flexible structure in a uniform flow, providing analytical solutions for the translating, axially flowing equilibria of strings subjected to a uniform body force and drag forces linear in the velocities. The classical catenaries are extended to a five-parameter family of curves. A sixth parameter affects the tension in the curves. Generic configurations are planar, represented by a single first order equation for the tangential angle. The effects of varying parameters on representative shapes, orbits in angle-curvature space, and stress distributions are shown. As limiting cases, the solutions include configurations corresponding to "lariat chains" and the towing, reeling, and sedimentation of flexible cables in a highly viscous fluid. Regions of parameter space corresponding to infinitely long, semi-infinite, and finite length curves are delineated. Almost all curves subtend an angle less than $\pi$ radians, but curious special cases with doubled or infinite range occur on the borders between regions. Separate transitions in the tension behavior, and counterintuitive results regarding finite towing tensions for infinitely long cables, are presented. Several physically inspired boundary value problems are solved and discussed.


title: Optimal control of nonlinear systems governed by Dirichlet fractional Laplacian in the minimax framework

abstract: We consider an optimal control problem governed by a class of boundary value problem with the spectral Dirichlet fractional Laplacian. Some sufficient condition for the existence of optimal processes is stated. The proof of the main result relies on variational structure of the problem. To show that boundary value problem with the Dirichlet fractional Laplacian has a weak solution we employ the renowned Ky Fan Theorem.


title: Inca Foams

abstract: We study a certain class of embedded two-foams that arise from gluing discs into ribbon torus knots along nonintersecting torus meridians. We exhibit several equivalent diagrammatic formalisms for these objects and identify several of their invariants, including a unique prime decomposition.


title: Counting non-standard binary representations

abstract: Let $\mathcal{A}$ be a finite subset of $\mathbb{N}$ including $0$ and $f_\mathcal{A}(n)$ be the number of ways to write $n=\sum_{i=0}^{\infty}\epsilon_i2^i$, where $\epsilon_i\in\mathcal{A}$. We consider asymptotics of the summatory function $s_\mathcal{A}(r,m)$ of $f_\mathcal{A}(n)$ from $m2^r$ to $m2^{r+1}-1$ and show that $s_{\mathcal{A}}(r,m)\approx c(\mathcal{A},m)\left|\mathcal{A}\right|^r$ for some $c(\mathcal{A},m)\in\mathbb{Q}$.


title: The fate of Hamilton's Hodograph in Special and General Relativity

abstract: The hodograph of a non-relativistic particle motion in Euclidean space is the curve described by its momentum vector. For a general central orbit problem the hodograph is the inverse of the pedal curve of the orbit, (i.e. its polar reciprocal), rotated through a right angle. Hamilton showed that for the Kepler/Coulomb problem, the hodograph is a circle whose centre is in the direction of a conserved eccentricity vector. The addition of an inverse cube law force induces the eccentricity vector to precess and with it the hodograph. The same effect is produced by a cosmic string. If one takes the relativistic momentum to define the hodograph, then for the Sommerfeld (i.e. the special relativistic Kepler/Coulomb problem) there is an effective inverse cube force which causes the hodograph to precess. If one uses Schwarzschild coordinates one may also define a a hodograph for timelike or null geodesics moving around a black hole. Iheir pedal equations are given. In special cases the hodograph may be found explicitly. For example the orbit of a photon which starts from the past singularity, grazes the horizon and returns to future singularity is a cardioid, its pedal equation is Cayley's sextic the inverse of which is Tschirhausen's cubic. It is also shown that that provided one uses Beltrami coordinates, the hodograph for the non-relativistic Kepler problem on hyperbolic space is also a circle. An analogous result holds for the the round 3-sphere. In an appendix the hodograph of a particle freely moving on a group manifold equipped with a left-invariant metric is defined.


title: Image Classification with Rejection using Contextual Information

abstract: We introduce a new supervised algorithm for image classification with rejection using multiscale contextual information. Rejection is desired in image-classification applications that require a robust classifier but not the classification of the entire image. The proposed algorithm combines local and multiscale contextual information with rejection, improving the classification performance. As a probabilistic model for classification, we adopt a multinomial logistic regression. The concept of rejection with contextual information is implemented by modeling the classification problem as an energy minimization problem over a graph representing local and multiscale similarities of the image. The rejection is introduced through an energy data term associated with the classification risk and the contextual information through an energy smoothness term associated with the local and multiscale similarities within the image. We illustrate the proposed method on the classification of images of H&E-stained teratoma tissues.