title: Sharing Means Renting?: An Entire-marketplace Analysis of Airbnb

abstract: Airbnb, an online marketplace for accommodations, has experienced a staggering growth accompanied by intense debates and scattered regulations around the world. Current discourses, however, are largely focused on opinions rather than empirical evidences. Here, we aim to bridge this gap by presenting the first large-scale measurement study on Airbnb, using a crawled data set containing 2.3 million listings, 1.3 million hosts, and 19.3 million reviews. We measure several key characteristics at the heart of the ongoing debate and the sharing economy. Among others, we find that Airbnb has reached a global yet heterogeneous coverage. The majority of its listings across many countries are entire homes, suggesting that Airbnb is actually more like a rental marketplace rather than a spare-room sharing platform. Analysis on star-ratings reveals that there is a bias toward positive ratings, amplified by a bias toward using positive words in reviews. The extent of such bias is greater than Yelp reviews, which were already shown to exhibit a positive bias. We investigate a key issue---commercial hosts who own multiple listings on Airbnb---repeatedly discussed in the current debate. We find that their existence is prevalent, they are early-movers towards joining Airbnb, and their listings are disproportionately entire homes and located in the US. Our work advances the current understanding of how Airbnb is being used and may serve as an independent and empirical reference to inform the debate.


title: Optical nanoscopy via quantum control

abstract: We present a scheme for nanoscopic imaging of a quantum mechanical two-level system using an optical probe in the far-field. Existing super-resolution schemes require more than two-levels and depend on an incoherent response to the lasers. Here, quantum control of the two states proceeds via rapid adiabatic passage. We implement this scheme on an array of semiconductor self-assembled quantum dots. Each quantum dot results in a bright spot in the image with extents down to 30 nm ({\lambda}/31). Rapid adiabatic passage is established as a versatile tool in the super-resolution toolbox.


title: Rayleigh-B\`enard convection in the generalized Oberbeck-Boussinesq system

abstract: The original Oberbeck-Boussinesq (OB) equations which are the coupled two dimensional Navier-Stokes(NS) and heat conduction equations have been investigated by E.N. Lorenz half a century ago with Fourier series and opened the way to the paradigm of chaos. In our former study-Chaos, Solitons and Fractals78, 249 (2015)-we presented fully analytic solutions for the velocity, pressure and temperature fields with the aim of the self-similar Ansatz and gave a possible explanation of the Rayleigh--B\`enard convection cells. Now we generalize the Oberbeck-Boussinesq hydrodynamical system, going beyond the first order Boussinesq approximation and consider a non-linear temperature coupling. We investigate more general, power law dependent fluid viscosity or heat conduction material equations as well. Our analytic results obtained via the self-similar Ansatz may attract the interest of various fields like meteorology, oceanography or climate studies.


title: Locality Sim: Cloud Simulator with Data Locality

abstract: Cloud Computing (CC) is a model for enabling on-demand access to a shared pool of configurable computing resources. Testing and evaluating the performance of the cloud environment for allocating, provisioning, scheduling, and data allocation policy have great attention to be achieved. Therefore, using cloud simulator would save time and money, and provide a flexible environment to evaluate new research work. Unfortunately, the current simulators (e.g., CloudSim, NetworkCloudSim, GreenCloud, etc..) deal with the data as for size only without any consideration about the data allocation policy and locality. On the other hand, the NetworkCloudSim simulator is considered one of the most common used simulators because it includes different modules which support needed functions to a simulated cloud environment, and it could be extended to include new extra modules. According to work in this paper, the NetworkCloudSim simulator has been extended and modified to support data locality. The modified simulator is called LocalitySim. The accuracy of the proposed LocalitySim simulator has been proved by building a mathematical model. Also, the proposed simulator has been used to test the performance of the three-tire data center as a case study with considering the data locality feature.


title: Existence of Some Signed Magic Arrays

abstract: We consider the notion of a signed magic array, which is an $m \times n$ rectangular array with the same number of filled cells $s$ in each row and the same number of filled cells $t$ in each column, filled with a certain set of numbers that is symmetric about the number zero, such that every row and column has a zero sum. We attempt to make progress toward a characterization of for which $(m, n, s, t)$ there exists such an array. This characterization is complete in the case where $n = s$ and in the case where $n = m$; we also characterize three-fourths of the cases where $n = 2m$.


title: Sparse Identification for Nonlinear Optical Communication Systems: SINO Method

abstract: We introduce low complexity machine learning based approach for mitigating nonlinear impairments in optical fiber communications systems. The immense intricacy of the problem calls for the development of "smart" methodology, simplifying the analysis without losing the key features that are important for recovery of transmitted data. The proposed sparse identification method for optical systems (SINO) allows to determine the minimal (optimal) number of degrees of freedom required for adaptive mitigation of detrimental nonlinear effects. We demonstrate successful application of the SINO method both for standard fiber communication links and for few-mode spatial-division-multiplexing systems.


title: Some new gradient estimates for two nonlinear parabolic equations under Ricci flow

abstract: In this paper, by maximum principle and cutoff function, we investigate gradient estimates for positive solutions to two nonlinear parabolic equations under Ricci flow. The related Harnack inequalities are deduced. An result about positive solutions on closed manifolds under Ricci flow is abtained. As applications, gradient estimates and Harnack inequalities for positive solutions to the heat equation under Ricci flow are derived. These results in the paper can be regard as generalizing the gradient estimates of Li-Yau, J. Y. Li, Hamilton and Li-Xu to the Ricci flow. Our results also improve the estimates of S. P. Liu and J. Sun to the nonlinear parabolic equation under Ricci flow.


title: Eikonal instability of Gauss-Bonnet-(anti-)de Sitter black holes

abstract: Here we have shown that asymptotically anti-de Sitter (AdS) black holes in the Einstein-Gauss-Bonnet (GB) theory are unstable under linear perturbations of spacetime in some region of parameters. This (eikonal) instability develops at high multipole numbers. We found the exact parametric regions of the eikonal instability and extended this consideration to asymptotically flat and de Sitter cases. The approach to the threshold of instability is driven by purely imaginary quasinormal modes, which are similar to those found recently in [5] for the higher curvature corrected black hole with the planar horizon. The found instability may indicate limits of holographic applicability of the GB-AdS backgrounds. Recently, through the analysis of critical behavior in AdS space-time in the presence of Gauss-Bonnet term, it was shown [19] that that, if the total energy content of the AdS space-time is small, then no black holes can be formed with mass less than some critical value. A similar mass gap was also found when considering collapse of mass shells in asymptotically flat Gauss-Bonnet theories [20]. The found instability of all sufficiently small Einstein-Gauss-Bonnet-AdS, -dS and asymptotically flat black holes may explain the existing mass gaps in their formation.


title: Bimeromorphic geometry of K\"ahler threefolds

abstract: We describe the recently established minimal model program for (non-algebraic) K\"ahler threefolds as well as the abundance theorem for these spaces.


title: Application of Fuzzy Logic in Design of Smart Washing Machine

abstract: Washing machine is of great domestic necessity as it frees us from the burden of washing our clothes and saves ample of our time. This paper will cover the aspect of designing and developing of Fuzzy Logic based, Smart Washing Machine. The regular washing machine (timer based) makes use of multi-turned timer based start-stop mechanism which is mechanical as is prone to breakage. In addition to its starting and stopping issues, the mechanical timers are not efficient with respect of maintenance and electricity usage. Recent developments have shown that merger of digital electronics in optimal functionality of this machine is possible and nowadays in practice. A number of international renowned companies have developed the machine with the introduction of smart artificial intelligence. Such a machine makes use of sensors and smartly calculates the amount of run-time (washing time) for the main machine motor. Realtime calculations and processes are also catered in optimizing the run-time of the machine. The obvious result is smart time management, better economy of electricity and efficiency of work. This paper deals with the indigenization of FLC (Fuzzy Logic Controller) based Washing Machine, which is capable of automating the inputs and getting the desired output (wash-time).


title: In vitro Demonstration of Cancer Inhibiting Properties from Stratified Self-Organized Micro-Discharge Plasma-Liquid Interface

abstract: Experiments have revealed a nontrivial cancer-inhibiting capability of liquid media treated by the plasma jet capable of forming thinly stratified self-organized patterns at a plasma-liquid interface. A pronounced cancer depressing activity towards at least two kinds of human cancer cells, namely breast cancer MDA-MB-231 and human glioblastoma U87 cancer lines, was demonstrated. After a short treatment at the thinly stratified self-organized plasma-liquid interface pattern, the cancer inhibiting media demonstrate well pronounced depression and apoptosis activities towards tumor cells, not achievable without interfacial stratification of plasma jet to thin (of several um) current filaments, which therefore play a pivotal (yet still not completely clear) role in building up the cancer inhibition properties. Moreover, thinly stratified, self-organized interfacial discharge is capable to efficiently control the ROS and RNS concentrations in the cancer-inhibiting media, and in particular, abnormal ROS/RNS ratios not achievable in discharges which do not form stratified thin-filament patterns could be obtained. These results were explained in terms of interaction of thin plasma filaments of the self-organized pattern with gas and liquid, where the unusual interaction conditions (i.e., high surface-to-volume ratios etc.) cause accumulation of ROS, RNS and other species in unusual ratios and concentrations, thus forming potentially efficient anti-cancer cocktail. Our funding could be extremely important for handling the cancer proliferation problem, and hence, it should be brought to light to attract due attention of the researchers and explore the possible potential of this approach in tackling the challenging problem of high cancer-induced mortality and rising morbidity trends.


title: A non-increasing tree growth process for recursive trees and applications

abstract: We introduce a non-increasing tree growth process $((T_n,\sigma_n),\, n\ge 1)$, where $T_n$ is a rooted labeled tree on $n$ vertices and ${\sigma}_n$ is a permutation of the vertex labels. The construction of $(T_{n},{\sigma}_n)$ from $(T_{n-1},{\sigma}_{n-1})$ involves rewiring a random (possibly empty) subset of edges in $T_{n-1}$ towards the newly added vertex; as a consequence $T_{n-1} \not\subset T_n$ with positive probability. The key feature of the process is that the shape of $T_n$ has the same law as that of a random recursive tree, while the degree distribution of any given vertex is not monotonous in the process. We present two applications. First, while couplings between Kingman's coalescent and random recursive trees where known for any fixed $n$, this new process provides a non-standard coupling of all finite Kingman's coalescents. Second, we use the new process and the Chen-Stein method to extend the well-understood properties of degree distribution of random recursive trees to extremal-range cases. Namely, we obtain convergence rates on the number of vertices with degree at least $c\ln n$, $c\in (1,2)$, in trees with $n$ vertices. Further avenues of research are discussed.


title: Autonomous Multirobot Excavation for Lunar Applications

abstract: In this paper, a control approach called Artificial Neural Tissue (ANT) is applied to multirobot excavation for lunar base preparation tasks including clearing landing pads and burying of habitat modules. We show for the first time, a team of autonomous robots excavating a terrain to match a given 3D blueprint. Constructing mounds around landing pads will provide physical shielding from debris during launch/landing. Burying a human habitat modules under 0.5 m of lunar regolith is expected to provide both radiation shielding and maintain temperatures of -25 $^{o}$C. This minimizes base life-support complexity and reduces launch mass. ANT is compelling for a lunar mission because it doesn't require a team of astronauts for excavation and it requires minimal supervision. The robot teams are shown to autonomously interpret blueprints, excavate and prepare sites for a lunar base. Because little pre-programmed knowledge is provided, the controllers discover creative techniques. ANT evolves techniques such as slot-dozing that would otherwise require excavation experts. This is critical in making an excavation mission feasible when it is prohibitively expensive to send astronauts. The controllers evolve elaborate negotiation behaviors to work in close quarters. These and other techniques such as concurrent evolution of the controller and team size are shown to tackle problem of antagonism, when too many robots interfere reducing the overall efficiency or worse, resulting in gridlock. While many challenges remain with this technology our work shows a compelling pathway for field testing this approach.


title: The next-to-minimal weights of binary projective Reed-Muller codes

abstract: Projective Reed-Muller codes were introduced by Lachaud, in 1988 and their dimension and minimum distance were determined by Serre and S{\o}rensen in 1991. In coding theory one is also interested in the higher Hamming weights, to study the code performance. Yet, not many values of the higher Hamming weights are known for these codes, not even the second lowest weight (also known as next-to-minimal weight) is completely determined. In this paper we determine all the values of the next-to-minimal weight for the binary projective Reed-Muller codes, which we show to be equal to the next-to-minimal weight of Reed-Muller codes in most, but not all, cases.


title: Towards an optofluidic pump?

abstract: Most of the vibrating mechanisms of optofluidicsystems are based on local heating of membranes that in-duces liquid flow. We report here a new type of diaphragmpump in a liquid film based on the optical radiation pres-sure force. We modulate a low power laser that gener-ates, at resonance, a symmetric vibration of a free standingsoap film. The film lifetime strongly varies from 56 s at lowpower (2 mW) to 2 s at higher power (70 mW). Since thelaser beam only acts mechanically on the interfaces, sucha pump could be easily implemented on delicate micro-equipment on chips or in biological systems.


title: Generic coexistence of Fermi arcs and Dirac cones on the surface of time-reversal invariant Weyl semimetals

abstract: The hallmark of Weyl semimetals is the existence of open constant-energy contours on their surface -- the so-called Fermi arcs -- connecting Weyl points. Here, we show that for time-reversal symmetric realizations of Weyl semimetals these Fermi arcs in many cases coexist with closed Fermi pockets originating from surface Dirac cones pinned to time-reversal invariant momenta. The existence of Fermi pockets is required for certain Fermi-arc connectivities due to additional restrictions imposed by the six $\mathbb{Z}_2$ topological invariants characterizing a generic time-reversal invariant Weyl semimetal. We show that a change of the Fermi-arc connectivity generally leads to a different topology of the surface Fermi surface, and identify the half-Heusler compound LaPtBi under in-plane compressive strain as a material that realizes this surface Lifshitz transition. We also discuss universal features of this coexistence in quasi-particle interference spectra.


title: Formation and evolution of exoplanets in different environments

abstract: The ultimate goal of exoplanetologists is to discover life outside our Earth and to fully understand our place in the Universe. Even though we have never been closer to attaining this goal, we still need to understand how and where the planets (efficiently) form. In this manuscript I briefly discuss the important role of stellar metallicity and chemistry on the formation and evolution of exoplanets.


title: Up- and down-quark contributions to the nucleon electromagnetic form factors at low $Q^2$

abstract: The spatial distribution of charge and magnetization within the nucleon (proton and neutron) is encoded in the elastic electromagnetic form factors $G_E^{(p,n)}$ and $G_M^{(p,n)}$. These form factors have been precisely measured utilizing elastic electron scattering, and the combination of proton and neutron form factors allows for the separation of the up- and down-quark contributions to the nucleon form factors. We expand on our original analyses and extract the up- and down-quark contributions to the nucleon electromagnetic form factors from worldwide data with an emphasis on precise new data covering the low-momentum region, which is sensitive to the large-scale structure of the nucleon. From these, we construct the flavor-separated Dirac and Pauli form factors and their ratios, and compare the results to recent extractions and theoretical calculations and models.


title: On the next-to-minimal weight of projective Reed-Muller codes

abstract: In this paper we present several values for the next-to-minimal weights of projective Reed-Muller codes. We work over $\mathbb{F}_q$ with $q \geq 3$ since in IEEE-IT 62(11) p. 6300-6303 (2016) we have determined the complete values for the next-to-minimal weights of binary projective Reed-Muller codes. As in loc. cit. here we also find examples of codewords with next-to-minimal weight whose set of zeros is not in a hyperplane arrangement.


title: An Integrated Conceptual Model for Information System Security Risk Management and Enterprise Architecture Management based on TOGAF, ArchiMate, IAF and DoDAF

abstract: Risk management is today a major steering tool for any organization wanting to deal with Information System (IS) security. However, IS Security Risk Management (ISSRM) remains difficult to establish and maintain, mainly in a context of multi-regulations with complex and inter-connected IS. We claim that a connection with Enterprise Architecture Management (EAM) contributes to deal with these issues. A first step towards a better integration of both domains is to define an integrated EAM-ISSRM conceptual model. Among the steps of the research method followed to define such an integrated EAM-ISSRM conceptual, this technical report presents the whole outputs (through alignment tables) of the conceptual alignment between concepts used to model EA (based on ArchiMate, TOGAF, IAF and DoDAF) and concepts of the ISSRM domain model.


title: Markov State Models from short non-Equilibrium Simulations - Analysis and Correction of Estimation Bias

abstract: Many state of the art methods for the thermodynamic and kinetic characterization of large and complex biomolecular systems by simulation rely on ensemble approaches, where data from large numbers of relatively short trajectories are integrated. In this context, Markov state models (MSMs) are extremely popular because they can be used to compute stationary quantities and long-time kinetics from ensembles of short simulations, provided that these short simulations are in "local equilibrium" within the MSM states. However, in the last over 15 years since the inception of MSMs, it has been controversially discussed and not yet been answered how deviations from local equilibrium can be detected, whether these deviations induce a practical bias in MSM estimation, and how to correct for them. In this paper, we address these issues: We systematically analyze the estimation of Markov state models (MSMs) from short non-equilibrium simulations, and we provide an expression for the error between unbiased transition probabilities and the expected estimate from many short simulations. We show that the unbiased MSM estimate can be obtained even from relatively short non-equilibrium simulations in the limit of long lag times and good discretization. Further, we exploit observable operator model (OOM) theory to derive an unbiased estimator for the MSM transition matrix that corrects for the effect of starting out of equilibrium, even when short lag times are used. Finally, we show how the OOM framework can be used to estimate the exact eigenvalues or relaxation timescales of the system without estimating an MSM transition matrix, which allows us to practically assess the discretization quality of the MSM. Applications to model systems and molecular dynamics simulation data of alanine dipeptide are included for illustration. The improved MSM estimator is implemented in PyEMMA as of version 2.3.

categories: physics.data-an cond-mat.stat-mech math.PR physics.bio-ph physics.chem-ph


title: From Foucault's Pendulum to the Gauss--Bonnet Theorem

abstract: We present a self-contained proof of the Gauss-Bonnet theorem for two-dimensional surfaces embedded in $R^3$ using just classical vector calculus. The exposition should be accessible to advanced undergraduate and non-expert graduate students. It may be viewed as an illustration and exercise in multivariate calculus and a motivation to go deeper into the fields of geometry and topology.


title: The geometry of a critical percolation cluster on the UIPT

abstract: We consider a critical Bernoulli site percolation on the uniform infinite planar triangulation. We study the tail distributions of the peeling time, perimeter, and volume of the hull of a critical cluster. The exponents obtained here differs by a factor 2 from those computed previously by Angel and Curien (2015) in the case of critical site percolation on the uniform infinite half-plane triangulation.


title: Disease Progression Modeling and Prediction through Random Effect Gaussian Processes and Time Transformation

abstract: The development of statistical approaches for the joint modelling of the temporal changes of imaging, biochemical, and clinical biomarkers is of paramount importance for improving the understanding of neurodegenerative disorders, and for providing a reference for the prediction and quantification of the pathology in unseen individuals. Nonetheless, the use of disease progression models for probabilistic predictions still requires investigation, for example for accounting for missing observations in clinical data, and for accurate uncertainty quantification. We tackle this problem by proposing a novel Gaussian process-based method for the joint modeling of imaging and clinical biomarker progressions from time series of individual observations. The model is formulated to account for individual random effects and time reparameterization, allowing non-parametric estimates of the biomarker evolution, as well as high flexibility in specifying correlation structure, and time transformation models. Thanks to the Bayesian formulation, the model naturally accounts for missing data, and allows for uncertainty quantification in the estimate of evolutions, as well as for probabilistic prediction of disease staging in unseen patients. The experimental results show that the proposed model provides a biologically plausible description of the evolution of Alzheimer's pathology across the whole disease time-span as well as remarkable predictive performance when tested on a large clinical cohort with missing observations.


title: Diffusion of new products with recovering consumers

abstract: We consider the diffusion of new products in the discrete Bass-SIR model, in which consumers who adopt the product can later "recover" and stop influencing their peers to adopt the product. To gain insight into the effect of the social network structure on the diffusion, we focus on two extreme cases. In the "most-connected" configuration where all consumers are inter-connected (complete network), averaging over all consumers leads to an aggregate model, which combines the Bass model for diffusion of new products with the SIR model for epidemics. In the "least-connected" configuration where consumers are arranged on a circle and each consumer can only be influenced by his left neighbor (one-sided 1D network), averaging over all consumers leads to a different aggregate model which is linear, and can be solved explicitly. We conjecture that for any other network, the diffusion is bounded from below and from above by that on a one-sided 1D network and on a complete network, respectively. When consumers are arranged on a circle and each consumer can be influenced by his left and right neighbors (two-sided 1D network), the diffusion is strictly faster than on a one-sided 1D network. This is different from the case of non-recovering adopters, where the diffusion on one-sided and on two-sided 1D networks is identical. We also propose a nonlinear model for recoveries, and show that consumers' heterogeneity has a negligible effect on the aggregate diffusion.


title: sTools - a data reduction pipeline for the GREGOR Fabry-P\'erot Interferometer and the High-resolution Fast Imager at the GREGOR solar telescope

abstract: A huge amount of data has been acquired with the GREGOR Fabry-P\'erot Interferometer (GFPI), large-format facility cameras, and since 2016 with the High-resolution Fast Imager (HiFI). These data are processed in standardized procedures with the aim of providing science-ready data for the solar physics community. For this purpose, we have developed a user-friendly data reduction pipeline called "sTools" based on the Interactive Data Language (IDL) and licensed under creative commons license. The pipeline delivers reduced and image-reconstructed data with a minimum of user interaction. Furthermore, quick-look data are generated as well as a webpage with an overview of the observations and their statistics. All the processed data are stored online at the GREGOR GFPI and HiFI data archive of the Leibniz Institute for Astrophysics Potsdam (AIP). The principles of the pipeline are presented together with selected high-resolution spectral scans and images processed with sTools.


title: Multi-level Compressed Sensing Petrov-Galerkin discretization of high-dimensional parametric PDEs

abstract: We analyze a novel multi-level version of a recently introduced compressed sensing (CS) Petrov-Galerkin (PG) method from [H. Rauhut and Ch. Schwab: Compressive Sensing Petrov-Galerkin approximation of high-dimensional parametric operator equations, Math. Comp. 304(2017) 661-700] for the solution of many-parametric partial differential equations. We propose to use multi-level PG discretizations, based on a hierarchy of nested finite dimensional subspaces, and to reconstruct parametric solutions at each level from level-dependent random samples of the high-dimensional parameter space via CS methods such as weighted l1-minimization. For affine parametric, linear operator equations, we prove that our approach allows to approximate the parametric solution with (almost) optimal convergence order as specified by certain summability properties of the coefficient sequence in a general polynomial chaos expansion of the parametric solution and by the convergence order of the PG discretization in the physical variables. The computations of the parameter samples of the PDE solution is "embarrassingly parallel", as in Monte-Carlo Methods. Contrary to other recent approaches, and as already noted in [A. Doostan and H. Owhadi: A non-adapted sparse approximation of PDEs with stochastic inputs. JCP 230(2011) 3015-3034] the optimality of the computed approximations does not require a-priori assumptions on ordering and structure of the index sets of the largest gpc coefficients (such as the "downward closed" property). We prove that under certain assumptions work versus accuracy of the new algorithms is asymptotically equal to that of one PG solve for the corresponding nominal problem on the finest discretization level up to a constant.


title: Detecting changes in slope with an $L_0$ penalty

abstract: Whilst there are many approaches to detecting changes in mean for a univariate time-series, the problem of detecting multiple changes in slope has comparatively been ignored. Part of the reason for this is that detecting changes in slope is much more challenging. For example, simple binary segmentation procedures do not work for this problem, whilst efficient dynamic programming methods that work well for the change in mean problem cannot be directly used for detecting changes in slope. We present a novel dynamic programming approach, CPOP, for finding the "best" continuous piecewise-linear fit to data. We define best based on a criterion that measures fit to data using the residual sum of squares, but penalises complexity based on an $L_0$ penalty on changes in slope. We show that using such a criterion is more reliable at estimating changepoint locations than approaches that penalise complexity using an $L_1$ penalty. Empirically CPOP has good computational properties, and can analyse a time-series with over 10,000 observations and over 100 changes in a few minutes. Our method is used to analyse data on the motion of bacteria, and provides fits to the data that both have substantially smaller residual sum of squares and are more parsimonious than two competing approaches.


title: Formal Proof of the Weak Goodstein Theorem

abstract: For many years, I have been interested in introducing students to the development of complex systems by means of modelling and refinement. To this end, I did not find anything better than presenting many examples of system developments. However, I figured out that my examples were not explicit enough on how (mechanical) proofs are performed. So, besides courses presenting these examples and also some courses in various forms of proofs (propositional calculus, first order predicate calculus, set theory), I decided to study the work of professional mathematicians, thinking that it could be good examples for students. Among the works I already studied and reconstructed are the theorem of Zermelo, the theorem of Cantor-Bernstein, the planar graph theorem of Kuratowski, the topological proof of the infinity of primes of Furstenberg, the intermediate value theorem of Bolzano, the Archimedean property of the set of Real numbers, and others. More recently, I found that the Goodstein theorem was also very interesting. The purpose of this short note is to give some information about this theorem and the way I introduce a weak form of it to students.


title: Existence and non-existence of minimal graphs

abstract: We study the Dirichlet problem for minimal surface systems in arbitrary dimension and codimension via mean curvature flow, and obtain the existence of minimal graphs over arbitrary mean convex bounded $C^2$ domains for a large class of prescribed boundary data. This result can be seen as a natural generalization of the classical sharp criterion for solvability of the minimal surface equation by Jenkins-Serrin. In contrast, we also construct a class of prescribed boundary data on just mean convex domains for which the Dirichlet problem in codimension 2 is not solvable. Moreover, we study existence and the uniqueness of minimal graphs by perturbation.


title: Pareto Efficient Multi Objective Optimization for Local Tuning of Analogy Based Estimation

abstract: Analogy Based Effort Estimation (ABE) is one of the prominent methods for software effort estimation. The fundamental concept of ABE is closer to the mentality of expert estimation but with an automated procedure in which the final estimate is generated by reusing similar historical projects. The main key issue when using ABE is how to adapt the effort of the retrieved nearest neighbors. The adaptation process is an essential part of ABE to generate more successful accurate estimation based on tuning the selected raw solutions, using some adaptation strategy. In this study we show that there are three interrelated decision variables that have great impact on the success of adaptation method: (1) number of nearest analogies (k), (2) optimum feature set needed for adaptation, and (3) adaptation weights. To find the right decision regarding these variables, one need to study all possible combinations and evaluate them individually to select the one that can improve all prediction evaluation measures. The existing evaluation measures usually behave differently, presenting sometimes opposite trends in evaluating prediction methods. This means that changing one decision variable could improve one evaluation measure while it is decreasing the others. Therefore, the main theme of this research is how to come up with best decision variables that improve adaptation strategy and thus, the overall evaluation measures without degrading the others. The impact of these decisions together has not been investigated before, therefore we propose to view the building of adaptation procedure as a multi-objective optimization problem. The Particle Swarm Optimization Algorithm (PSO) is utilized to find the optimum solutions for such decision variables based on optimizing multiple evaluation measures


title: SD-CPS: Taming the Challenges of Cyber-Physical Systems with a Software-Defined Approach

abstract: Cyber-Physical Systems (CPS) revolutionize various application domains with integration and interoperability of networking, computing systems, and mechanical devices. Due to its scale and variety, CPS faces a number of challenges and opens up a few research questions in terms of management, fault-tolerance, and scalability. We propose a software-defined approach inspired by Software-Defined Networking (SDN), to address the challenges for a wider CPS adoption. We thus design a middleware architecture for the correct and resilient operation of CPS, to manage and coordinate the interacting devices centrally in the cyberspace whilst not sacrificing the functionality and performance benefits inherent to a distributed execution.


title: The Shapley Value of Digraph Games

abstract: In this paper the Shapley value of digraph (directed graph) games are considered. Digraph games are transferable utility (TU) games with limited cooperation among players, where players are represented by nodes. A restrictive relation between two adjacent players is established by a directed line segment. Directed path, connecting the initial player with the terminal player, form the coalition among players. A dominance relation is established between players and this relation determines whether or not a player wants to cooperate. To cooperate, we assume that a player joins a coalition where he/she is not dominated by any other players.The Shapley value is defined as the average of marginal contribution vectors corresponding to all permutations that do not violate the subordination of players. The Shapley value for cyclic digraph games is calculated and analyzed. For a given family of characteristic functions, a quick way to calculate Shapley values is formulated.


title: An Optics-Based Approach to Thermal Management of Photovoltaics: Selective-Spectral and Radiative Cooling

abstract: For commercial one-sun solar modules, up to 80% of the incoming sunlight may be dissipated as heat, potentially raising the temperature 20 C - 30 C higher than the ambient. In the long term, extreme self-heating erodes efficiency and shortens lifetime, thereby dramatically reducing the total energy output. Therefore, it is critically important to develop effective and practical (and preferably passive) cooling methods to reduce operating temperature of PV modules. In this paper, we explore two fundamental (but often overlooked) origins of PV self-heating, namely, sub-bandgap absorption and imperfect thermal radiation. The analysis suggests that we redesign the optical properties of the solar module to eliminate parasitic absorption (selective-spectral cooling) and enhance thermal emission (radiative cooling). Our Comprehensive opto-electro-thermal simulation shows that the proposed techniques would cool the one-sun and low-concentrated terrestrial solar modules up to 10 C and 20 C, respectively. This self-cooling would substantially extend the lifetime for solar modules, with The corresponding increase in energy yields and reduced LCOE.


title: Multigraded linear series and recollement

abstract: Given a scheme $Y$ equipped with a collection of globally generated vector bundles $E_1, \dots, E_n$, we study the universal morphism from $Y$ to a fine moduli space $\mathcal{M}(E)$ of cyclic modules over the endomorphism algebra of $E:=\mathcal{O}_Y\oplus E_1\oplus\cdots \oplus E_n$. This generalises the classical morphism to the linear series of a basepoint-free line bundle on a scheme. We describe the image of the morphism and present necessary and sufficient conditions for surjectivity in terms of a recollement of a module category. When the morphism is surjective, this gives a fine moduli space interpretation of the image, and as an application we show that for a small, finite subgroup $G\subset \text{GL}(2,k)$, every sub-minimal partial resolution of $\mathbb{A}^2_k/G$ is isomorphic to a fine moduli space $\mathcal{M}(E_C)$ where $E_C$ is a summand of the bundle $E$ defining the reconstruction algebra. We also consider applications to Gorenstein affine threefolds, where Reid's recipe sheds some light on the classes of algebra from which one can reconstruct a given crepant resolution.


title: $\frak{g}$-quasi-Frobenius Lie algebras

abstract: A Lie version of Turaev's $\overline{G}$-Frobenius algebras from 2-dimensional homotopy quantum field theory is proposed. The foundation for this Lie version is a structure we call a \textit{$\frak{g}$-quasi-Frobenius Lie algebra} for $\frak{g}$ a finite dimensional Lie algebra. The latter consists of a quasi-Frobenius Lie algebra $(\frak{q},\beta)$ together with a left $\frak{g}$-module structure which acts on $\frak{q}$ via derivations and for which $\beta$ is $\frak{g}$-invariant. Geometrically, $\frak{g}$-quasi-Frobenius Lie algebras are the Lie algebra structures associated to symplectic Lie groups with an action by a Lie group $G$ which acts via symplectic Lie group automorphisms. In addition to geometry, $\frak{g}$-quasi-Frobenius Lie algebras can also be motivated from the point of view of category theory. Specifically, $\frak{g}$-quasi Frobenius Lie algebras correspond to \textit{quasi Frobenius Lie objects} in $\mathbf{Rep}(\frak{g})$. If $\frak{g}$ is now equipped with a Lie bialgebra structure, then the categorical formulation of $\overline{G}$-Frobenius algebras given in \cite{KP} suggests that the Lie version of a $\overline{G}$-Frobenius algebra is a quasi-Frobenius Lie object in $\mathbf{Rep}(D(\frak{g}))$, where $D(\frak{g})$ is the associated (semiclassical) Drinfeld double. We show that if $\frak{g}$ is a quasitriangular Lie bialgebra, then every $\frak{g}$-quasi-Frobenius Lie algebra has an induced $D(\frak{g})$-action which gives it the structure of a $D(\frak{g})$-quasi-Frobenius Lie algebra.


title: Coupling between Ion-Acoustic Waves and Neutrino Oscillations

abstract: The work investigates the coupling between ion-acoustic waves and neutrino flavor oscillations in a non-relativistic electron-ion plasma under the influence of a mixed neutrino beam. Neutrino oscillations are mediated by the flavor polarization vector dynamics in a material medium. The linear dispersion relation around homogeneous static equilibria is developed. When resonant with the ion-acoustic mode, the neutrino flavor oscillations can transfer energy to the plasma exciting a new fast unstable mode in extreme astrophysical scenarios. The growth rate and the unstable wavelengths are determined in typical type II supernovae parameters. The predictions can be useful for a new indirect probe on neutrino oscillations in nature.


title: A STEREO Survey of Magnetic Cloud Coronal Mass Ejections Observed at Earth in 2008--2012

abstract: We identify coronal mass ejections (CMEs) associated with magnetic clouds (MCs) observed near Earth by the Wind spacecraft from 2008 to mid-2012, a time period when the two STEREO spacecraft were well positioned to study Earth-directed CMEs. We find 31 out of 48 Wind MCs during this period can be clearly connected with a CME that is trackable in STEREO imagery all the way from the Sun to near 1 AU. For these events, we perform full 3-D reconstructions of the CME structure and kinematics, assuming a flux rope morphology for the CME shape, considering the full complement of STEREO and SOHO imaging constraints. We find that the flux rope orientations and sizes inferred from imaging are not well correlated with MC orientations and sizes inferred from the Wind data. However, velocities within the MC region are reproduced reasonably well by the image-based reconstruction. Our kinematic measurements are used to provide simple prescriptions for predicting CME arrival times at Earth, provided for a range of distances from the Sun where CME velocity measurements might be made. Finally, we discuss the differences in the morphology and kinematics of CME flux ropes associated with different surface phenomena (flares, filament eruptions, or no surface activity).


title: High Speed Elephant Flow Detection Under Partial Information

abstract: In this paper we introduce a new framework to detect elephant flows at very high speed rates and under uncertainty. The framework provides exact mathematical formulas to compute the detection likelihood and introduces a new flow reconstruction lemma under partial information. These theoretical results lead to the design of BubbleCache, a new elephant flow detection algorithm designed to operate near the optimal tradeoff between computational scalability and accuracy by dynamically tracking the traffic's natural cutoff sampling rate. We demonstrate on a real world 100 Gbps network that the BubbleCache algorithm helps reduce the computational cost by a factor of 1000 and the memory requirements by a factor of 100 while detecting the top flows on the network with very high probability.


title: Full Momentum and Energy Resolved Spectral Function of a 2D Electronic System

abstract: The single-particle spectral function measures the density of electronic states (DOS) in a material as a function of both momentum and energy, providing central insights into phenomena such as superconductivity and Mott insulators. While scanning tunneling microscopy (STM) and other tunneling methods have provided partial spectral information, until now only angle-resolved photoemission spectroscopy (ARPES) has permitted a comprehensive determination of the spectral function of materials in both momentum and energy. However, ARPES operates only on electronic systems at the material surface and cannot work in the presence of applied magnetic fields. Here, we demonstrate a new method for determining the full momentum and energy resolved electronic spectral function of a two-dimensional (2D) electronic system embedded in a semiconductor. In contrast with ARPES, the technique remains operational in the presence of large externally applied magnetic fields and functions for electronic systems with zero electrical conductivity or with zero electron density. It provides a direct high-resolution and high-fidelity probe of the dispersion and dynamics of the interacting 2D electron system. By ensuring the system of interest remains under equilibrium conditions, we uncover delicate signatures of many-body effects involving electron-phonon interactions, plasmons, polarons, and a novel phonon analog of the vacuum Rabi splitting in atomic systems.


title: Imaging magnetic vortex configurations in ferromagnetic nanotubes

abstract: We image the remnant magnetization configurations of CoFeB and permalloy nanotubes (NTs) using x-ray magnetic circular dichroism photo-emission electron microscopy. The images provide direct evidence for flux-closure configurations, including a global vortex state, in which magnetization points circumferentially around the NT axis. Furthermore, micromagnetic simulations predict and measurements confirm that vortex states can be programmed as the equilibrium remnant magnetization configurations by reducing the NT aspect ratio.


title: Trion and Biexciton in Monolayer Transition Metal Dichalcogenides

abstract: We study the trion and biexciton in transition metal dichalcogenides monolayers within the framework of a nonrelativistic potential model using the method of hyperspherical harmonics (HH). We solve the three- and four-body Schr\"odinger equations with the Keldysh potential by expanding the wave functions of a trion and biexciton in terms of the antisymmetrized HH. Results of the calculations for the ground state energies are in good agreement with similar calculations for the Keldysh potential and in reasonable agreement with experimental measurements of trion and biexciton binding energies.


title: Deep Convolutional Denoising of Low-Light Images

abstract: Poisson distribution is used for modeling noise in photon-limited imaging. While canonical examples include relatively exotic types of sensing like spectral imaging or astronomy, the problem is relevant to regular photography now more than ever due to the booming market for mobile cameras. Restricted form factor limits the amount of absorbed light, thus computational post-processing is called for. In this paper, we make use of the powerful framework of deep convolutional neural networks for Poisson denoising. We demonstrate how by training the same network with images having a specific peak value, our denoiser outperforms previous state-of-the-art by a large margin both visually and quantitatively. Being flexible and data-driven, our solution resolves the heavy ad hoc engineering used in previous methods and is an order of magnitude faster. We further show that by adding a reasonable prior on the class of the image being processed, another significant boost in performance is achieved.


title: A Multiscale-Analysis of Stochastic Bistable Reaction-Diffusion Equations

abstract: A multiscale analysis of 1D stochastic bistable reaction-diffusion equations with additive noise is carried out w.r.t. travelling waves within the variational approach to stochastic partial differential equations. It is shown with explicit error estimates on appropriate function spaces that up to lower order w.r.t. the noise amplitude, the solution can be decomposed into the orthogonal sum of a travelling wave moving with random speed and into Gaussian fluctuations. A stochastic differential equation describing the speed of the travelling wave and a linear stochastic partial differential equation describing the fluctuations are derived in terms of the coefficients. Our results extend corresponding results obtained for stochastic neural field equations to the present class of stochastic dynamics.


title: Collective optomechanical effects in cavity quantum electrodynamics

abstract: We investigate a cavity quantum electrodynamic effect, where the alignment of two-dimensional freely rotating optical dipoles is driven by their collective coupling to the cavity field. By exploiting the formal equivalence of a set of rotating dipoles with a polymer we calculate the partition function of the coupled light-matter system and demonstrate it exhibits a second order phase transition between a bunched state of isotropic orientations and a stretched one with all the dipoles aligned. Such a transition manifests itself as an intensity-dependent shift of the polariton mode resonance. Our work, lying at the crossroad between cavity quantum electrodynamics and quantum optomechanics, is a step forward in the on-going quest to understand how strong coupling can be exploited to influence matter internal degrees of freedom.


title: Parametric decay of plasma waves near the upper-hybrid resonance

abstract: An intense X wave propagating perpendicularly to dc magnetic field is unstable with respect to a parametric decay into an electron Bernstein wave and a lower-hybrid wave. A modified theory of this effect is proposed that extends to the high-intensity regime, where the instability rate $\gamma$ ceases to be a linear function of the incident-wave amplitude. An explicit formula for $\gamma$ is derived and expressed in terms of cold-plasma parameters. Theory predictions are in reasonable agreement with the results of the particle-in-cell simulations presented in a separate publication.


title: Magnetization in narrow ribbons: curvature effects

abstract: A ribbon is a surface swept out by a line segment turning as it moves along a central curve. For narrow magnetic ribbons, for which the length of the line segment is much less than the length of the curve, the anisotropy induced by the magnetostatic interaction is biaxial, with hard axis normal to the ribbon and easy axis along the central curve. The micromagnetic energy of a narrow ribbon reduces to that of a one-dimensional ferromagnetic wire, but with curvature, torsion and local anisotropy modified by the rate of turning. These general results are applied to two examples, namely a helicoid ribbon, for which the central curve is a straight line, and a M\"obius ribbon, for which the central curve is a circle about which the line segment executes a $180^\circ$ twist. In both examples, for large positive tangential anisotropy, the ground state magnetization lies tangent to the central curve. As the tangential anisotropy is decreased, the ground state magnetization undergoes a transition, acquiring an in-surface component perpendicular to the central curve. For the helicoid ribbon, the transition occurs at vanishing anisotropy, below which the ground state is uniformly perpendicular to the central curve. The transition for the M\"obius ribbon is more subtle; it occurs at a positive critical value of the anisotropy, below which the ground state is nonuniform. For the helicoid ribbon, the dispersion law for spin wave excitations about the tangential state is found to exhibit an asymmetry determined by the geometric and magnetic chiralities.


title: To Boost or Not to Boost? On the Limits of Boosted Trees for Object Detection

abstract: We aim to study the modeling limitations of the commonly employed boosted decision trees classifier. Inspired by the success of large, data-hungry visual recognition models (e.g. deep convolutional neural networks), this paper focuses on the relationship between modeling capacity of the weak learners, dataset size, and dataset properties. A set of novel experiments on the Caltech Pedestrian Detection benchmark results in the best known performance among non-CNN techniques while operating at fast run-time speed. Furthermore, the performance is on par with deep architectures (9.71% log-average miss rate), while using only HOG+LUV channels as features. The conclusions from this study are shown to generalize over different object detection domains as demonstrated on the FDDB face detection benchmark (93.37% accuracy). Despite the impressive performance, this study reveals the limited modeling capacity of the common boosted trees model, motivating a need for architectural changes in order to compete with multi-level and very deep architectures.


title: Emergent stochastic oscillations and signal detection in tree networks of excitable elements

abstract: We study the stochastic dynamics of strongly-coupled excitable elements on a tree network. The peripheral nodes receive independent random inputs which may induce large spiking events propagating through the branches of the tree and leading to global coherent oscillations in the network. This scenario may be relevant to action potential generation in certain sensory neurons, which possess myelinated distal dendritic tree-like arbors with excitable nodes of Ranvier at peripheral and branching nodes and exhibit noisy periodic sequences of action potentials. We focus on the spiking statistics of the central node, which fires in response to a noisy input at peripheral nodes. We show that, in the strong coupling regime, relevant to myelinated dendritic trees, the spike train statistics can be predicted from an isolated excitable element with rescaled parameters according to the network topology. Furthermore, we show that by varying the network topology the spike train statistics of the central node can be tuned to have a certain firing rate and variability, or to allow for an optimal discrimination of inputs applied at the peripheral nodes.

categories: physics.bio-ph cond-mat.dis-nn nlin.AO q-bio.NC


title: Speedmeter scheme for gravitational-wave detectors based on EPR quantum entanglement

abstract: We propose a new implementation of a quantum speed meter QND measurement scheme. It employs two independent optical readouts of the interferometer test masses, featuring strongly different values of the bandwidths $\gamma_{1,2}$ and of the optical circulating power $I_{1,2}$, with the special relationship of $I_1/I_2=\gamma_1/\gamma_2$. The outputs of these two position meters have to be combined by an additional beamsplitter. In this scheme, signals at the common and the differential outputs of the interferometer setup are proportional to the position and the velocity of the test masses, respectively. The influence of the position meter-like back action force associated with the position signal can be cancelled using the EPR approach by measuring the amplitude quadrature of the beamsplitter common output correlated with this force. In the standard signal-recycled Michelson interferometer topology of the modern gravitational-wave detectors, two independent optical position meters can be implemented by two orthogonal polarisations of the probe light. Our analysis shows that the EPR speedmeter provides significantly improved sensitivity for all frequencies below $\sim 30\,{\rm Hz}$ compared to an equivalent signal recycled Michelson interferometer. We believe the EPR speedmeter scheme to be very attractive for future upgrades of gravitational wave detectors, because it requires only minor changes to be implemented in the interferometer hardware and allows to switch between the position meter and the speed meter modes within short time-scales and without any changes to the hardware.


title: Star formation driven galactic winds in UGC 10043

abstract: We study the galactic wind in the edge-on spiral galaxy UGC 10043 with the combination of the CALIFA integral field spectroscopy data, scanning Fabry-Perot interferometry (FPI), and multiband photometry. We detect ionized gas in the extraplanar regions reaching a relatively high distance, up to ~ 4 kpc above the galactic disk. The ionized gas line ratios ([N ii]/Ha, [S ii]/Ha and [O i]/Ha) present an enhancement along the semi minor axis, in contrast with the values found at the disk, where they are compatible with ionization due to H ii-regions. These differences, together with the biconic symmetry of the extra-planar ionized structure, makes UGC 10043 a clear candidate for a galaxy with gas outflows ionizated by shocks. From the comparison of shock models with the observed line ratios, and the kinematics observed from the FPI data, we constrain the physical properties of the observed outflow. The data are compatible with a velocity increase of the gas along the extraplanar distances up to < 400 km/s and the preshock density decreasing in the same direction. We also observe a discrepancy in the SFR estimated based on Ha (0.36 Msun/yr ) and the estimated with the CIGALE code, being the latter 5 times larger. Nevertheless, this SFR is still not enough to drive the observed galactic wind if we do not take into account the filling factor. We stress that the combination of the three techniques of observation with models is a powerful tool to explore galactic winds in the Local Universe.


title: Flux tubes at Finite Temperature

abstract: In this work, we show the flux tubes of the quark-antiquark and quark-quark at finite temperature for SU(3) Lattice QCD. The chromomagnetic and chromoelectric fields are calculated above and below the phase transition.


title: Quantum-continuum simulation of the electrochemical response of pseudocapacitor electrodes from realistic conditions

abstract: Pseudocapacitors are energy-storage devices characterized by fast and reversible redox reactions that enable them to store large amounts of electrical energy at high rates. We simulate the response of pseudocapacitive electrodes under realistic conditions to identify the microscopic factors that determine their performance, focusing on ruthenia (RuO2) as a prototypical electrode material. Electronic-structure methods are used together with a self-consistent continuum solvation (SCCS) model to build a complete dataset of free energies as the surface of the charged electrode is gradually covered with protons under applied voltage. The resulting dataset is exploited to compute hydrogen-adsorption isotherms and charge-voltage responses by means of grand-canonical sampling, finding close agreement with experimental voltammetry. These simulations reveal that small changes on the order of 5 {\mu}F/cm2 in the intrinsic double-layer capacitance of the electrode-electrolyte interface can induce variations of up to 40 {\mu}F/cm2 in the overall pseudocapacitance.


title: Deep Class Aware Denoising

abstract: The increasing demand for high image quality in mobile devices brings forth the need for better computational enhancement techniques, and image denoising in particular. At the same time, the images captured by these devices can be categorized into a small set of semantic classes. However simple, this observation has not been exploited in image denoising until now. In this paper, we demonstrate how the reconstruction quality improves when a denoiser is aware of the type of content in the image. To this end, we first propose a new fully convolutional deep neural network architecture which is simple yet powerful as it achieves state-of-the-art performance even without being class-aware. We further show that a significant boost in performance of up to $0.4$ dB PSNR can be achieved by making our network class-aware, namely, by fine-tuning it for images belonging to a specific semantic class. Relying on the hugely successful existing image classifiers, this research advocates for using a class-aware approach in all image enhancement tasks.


title: Finite time singularity formation for moving interface Euler equations

abstract: This paper proposes a new general methodology for finite-time singularity formation for moving interface problems involving the incompressible Euler equations in the plane. The first problem considered is the two-phase Euler vortex sheets problem with surface tension for which is proved the finite time singularity of the natural norm of the problem for suitable initial data. This is in striking contrast with the case of finite time splash and splat singularity formation for the one phase Euler equations introduced in [4] and studied in a more general context in [8], for which the natural norm (in the one phase fluid) stays finite all the way until contact. The second problem considered involves the presence of a heavier rigid body moving in the perfect fluid. We first establish that the rigid body will hit the bottom of the fluid domain in finite time (in a more general context and very different methodology than first done for this problem in [19]). We next establish that a surface energy blows up, and characterize the acceleration of the rigid body at contact, depending on the nature of the contact zone: Acceleration opposes the motion at contact, and is either positive finite if the contact zone contains a curve, or infinite otherwise.


title: Exploring the squark flavour structure of the MSSM

abstract: We present an extensive study of the MSSM parameter space allowing for general generation mixing in the squark sector. Employing an MCMC algorithm, we establish the parameter ranges which are allowed with respect to various experimental and theoretical constraints. Based on this analysis, we propose benchmark scenarios for future studies. Moreover, we discuss aspects of signatures at the LHC.


title: Consistency relations for the Lagrangian halo bias and their implications

abstract: The protohalo patches from which halos form are defined by a number of constraints imposed on the Lagrangian dark matter density field. Each of these constraints contributes to biasing the spatial distribution of the protohalos relative to the matter. We show how measurements of this spatial distribution -- linear combinations of protohalo bias factors -- can be used to make inferences about the physics of halo formation. Our analysis exploits the fact that halo bias factors satisfy consistency relations which encode this physics, and that these relations are the same even for sub-populations in which assembly bias has played a role. We illustrate our methods using a model in which three parameters matter: a density threshold, the local slope and the curvature of the smoothed density field. The latter two are nearly degenerate; our approach naturally allows one to build an accurate effective two-parameter model for which the consistency relations still apply. This, with an accurate description of the smoothing window, allows one to describe the protohalo-matter cross-correlation very well, both in Fourier and configuration space. We then use our determination of the large scale bias parameters together with the consistency relations, to estimate the enclosed density and mean slope on the Lagrangian radius scale of the protohalos. Direct measurements of these quantities, made on smaller scales than those on which the bias parameters are typically measured, are in good agreement.


title: Virtual Network Migration on the GENI Wide-Area SDN-Enabled Infrastructure

abstract: A virtual network (VN) contains a collection of virtual nodes and links assigned to underlying physical resources in a network substrate. VN migration is the process of remapping a VN's logical topology to a new set of physical resources to provide failure recovery, energy savings, or defense against attack. Providing VN migration that is transparent to running applications is a significant challenge. Efficient migration mechanisms are highly dependent on the technology deployed in the physical substrate. Prior work has considered migration in data centers and in the PlanetLab infrastructure. However, there has been little effort targeting an SDN-enabled wide-area networking environment - an important building block of future networking infrastructure. In this work, we are interested in the design, implementation and evaluation of VN migration in GENI as a working example of such a future network. We identify and propose techniques to address key challenges: the dynamic allocation of resources during migration, managing hosts connected to the VN, and flow table migration sequences to minimize packet loss. We find that GENI's virtualization architecture makes transparent and efficient migration challenging. We suggest alternatives that might be adopted in GENI and are worthy of adoption by virtual network providers to facilitate migration.


title: Energetic instability of passive states in thermodynamics

abstract: Passivity is a fundamental concept in thermodynamics that demands a quantum system's energy cannot be lowered by any reversible, unitary process acting on the system. In the limit of many such systems, passivity leads in turn to the concept of complete passivity, thermal states, and the emergence of a thermodynamic temperature. In contrast, here we need only consider a single system and show that every passive state except the thermal state is unstable under a weaker form of reversibility. More precisely, we show that given a single copy of any athermal quantum state we may extract a maximal amount of energy from the state when we can use a machine that operates in a reversible cycle and whose state is left unchanged. This means that for individual systems the only form of passivity that is stable under general reversible processes is complete passivity, and thus provides a single-shot and more physically motivated identification of thermal states and the emergence of temperature. The machine which extracts work from passive states exploits the fact that one can find a subspace which acts as a virtual hot reservoir, and a subspace which acts as a virtual cold reservoir. We show that an optimal amount of work can be extracted, and that the machine operates at the Carnot efficiency between pairs of virtual reservoirs.


title: Near Horizon Geometries and Black Hole Holograph

abstract: Two quasi local approaches to black holes are combined: Near Horizon Geometries (NHG) and stationary Black Hole Holographs (BHH). Necessary and sufficient conditions on BHH data for the emergence of NHGs as resulting vacuum solutions to Einstein's equations are found.


title: Moving planes, Jacobi curves and the dynamical approach to Finsler geometry

abstract: We express invariants of Finsler manifolds in a geometrical way by means of using moving planes and their associated Jacobi curves, which are curves in a fixed homogeneous Grassmann manifold. Some applications are given.


title: Effects on the CMB from Compactification Before Inflation

abstract: Many theories beyond the Standard Model include extra dimensions, though these have yet to be directly observed. In this work we consider the possibility of a compactification mechanism which both allows extra dimensions and is compatible with current observations. This compactification is predicted to leave a signature on the CMB by altering the amplitude of the low l multipoles, dependent on the amount of inflation. Recently discovered CMB anomalies at low multipoles may be evidence for this. In our model we assume the spacetime is the product of a four-dimensional spacetime and flat extra dimensions. Before the compactification, both the four-dimensional space- time and the extra dimensions can either be expanding or contracting independently. Taking into account physical constraints, we explore the observational consequences and the plausibility of these different models.


title: On Non-Linear Markov Operators: surjectivity vs orthogonal preserving property

abstract: In the present paper, we consider nonlinear Markov operators, namely polynomial stochastic operators. We introduce a notion of orthogonal preserving polynomial stochastic operators. The purpose of this study is to show that surjectivity of nonlinear Markov operators is equivalent to their orthogonal preserving property.


title: Higher order multipoles in metamaterial homogenization

abstract: The higher order multipoles above the electric quadrupole are commonly neglected in metamaterial homogenization. We show that they nevertheless can be significant when second order spatial dispersive effects, such as the magnetic response, are considered. In this respect, they can be equally important as the magnetization and quadrupole terms, and should not automatically be neglected.


title: A new approximation method for geodesics on the space of K\"ahler metrics using complexified symplectomorphisms and Gr\"obner Lie series

abstract: It has been shown that the Cauchy problem for geodesics in the space of K\"ahler metrics with a fixed cohomology class on a compact complex manifold $M$ can be effectively reduced to the problem of finding the flow of a related hamiltonian vector field $X_H$, followed by analytic continuation of the time to complex time. This opens the possibility of expressing the geodesic $\omega_t$ in terms of Gr\"obner Lie series of the form $\exp(\sqrt{-1} \, tX_H)(f)$, for local holomorphic functions $f$. The main goal of this paper is to use truncated Lie series as a new way of constructing approximate solutions to the geodesic equation. For the case of an elliptic curve and $H$ a certain Morse function squared, we approximate the relevant Lie series by their first twelve terms, calculated with the help of Mathematica. This leads to approximate geodesics which hit the boundary of the space of K\"ahler metrics in finite geodesic time. For quantum mechanical applications, one is interested also on the non-K\"ahler polarizations that one obtains by crossing the boundary of the space of K\"ahler structures. Properties of the approximate geodesics and its extensions are also studied using Mathematica.


title: An inexact strategy for the projected gradient algorithm in vector optimization problems on variable ordered spaces

abstract: Variable order structures model situations in which the comparison between two points depends on a point-to-cone map. In this paper, an inexact projected gradient method for solving smooth constrained vector optimization problems on variable ordered spaces is presented. It is shown that every accumulation point of the generated sequence satisfies the first order necessary optimality condition. The convergence of all accumulation points to a weakly efficient point is established under suitable convexity assumptions for the objective function. The convergence results are also derived in the particular case in which the problem is unconstrained and if exact directions are taken as descent directions. Furthermore, we investigate the application of the proposed method to optimization models where the domain of the variable order map and the objective function are the same. In this case, similar concepts and convergence results are presented. Finally, some computational experiments designed to illustrate the behavior of the proposed inexact methods versus the exact ones (in terms of CPU time) are performed.


title: Functions on surfaces and constructions of manifolds in dimensions three, four and five

abstract: We offer a new proof that two closed oriented 4-manifolds are cobordant if their signatures agree, in the spirit of Lickorish's proof that all closed oriented 3-manifolds bound 4-manifolds. Where Lickorish uses Heegaard splittings we use trisections. In fact we begin with a subtle recasting of Lickorish's argument: Instead of factoring the gluing map for a Heegaard splitting as a product of Dehn twists, we encode each handlebody in a Heegaard splitting in terms of a Morse function on the surface and build the 4-manifold from a generic homotopy between the two functions. This extends up a dimension by encoding a trisection of a closed 4-manifold as a triangle (circle) of functions and constructing an associated 5-manifold from an extension to a 2-simplex (disk) of functions. This borrows ideas from Hatcher and Thurston's proof that the mapping class group of a surface is finitely presented.


title: Internal delensing of Planck CMB temperature and polarization

abstract: We present a first internal delensing of CMB maps, both in temperature and polarization, using the public foreground-cleaned (SMICA) Planck 2015 maps. After forming quadratic estimates of the lensing potential, we use the corresponding displacement field to undo the lensing on the same data. We build differences of the delensed spectra to the original data spectra specifically to look for delensing signatures. After taking into account reconstruction noise biases in the delensed spectra, we find an expected sharpening of the power spectrum acoustic peaks with a delensing efficiency of $29\,\%$ ($TT$) $25\,\%$ ($TE$) and $22\,\%$ ($EE$). The detection significance of the delensing effects is very high in all spectra: $12\,\sigma$ in $EE$ polarization; $18\,\sigma$ in $TE$; and $20\,\sigma$ in $TT$. The null hypothesis of no lensing in the maps is rejected at $26\,\sigma$. While direct detection of the power in lensing $B$-modes themselves is not possible at high significance at Planck noise levels, we do detect (at $4.5\,\sigma$ under the null hypothesis) delensing effects in the $B$-mode map, with $7\,\%$ reduction in lensing power. Our results provide a first demonstration of polarization delensing, and generally of internal CMB delensing, and stand in agreement with the baseline $\Lambda$CDM Planck 2015 cosmology expectations.


title: Notes on "Einstein metrics on compact simple Lie groups attached to standard triples"

abstract: In the paper "Einstein metrics on compact simple Lie groups attached to standard triples", the authors introduced the definition of standard triples and proved that every compact simple Lie group $G$ attached to a standard triple $(G,K,H)$ admits a left-invariant Einstein metric which is not naturally reductive except the standard triple $(\Sp(4),2\Sp(2),4\Sp(1))$. For the triple $(\Sp(4),2\Sp(2),4\Sp(1))$, we find there exists an involution pair of $\sp(4)$ such that $4\sp(1)$ is the fixed point of the pair, and then give the decomposition of $\sp(4)$ as a direct sum of irreducible $\ad(4\sp(1))$-modules. But $\Sp(4)/4\Sp(1)$ is not a generalized Wallach space. Furthermore we give left-invariant Einstein metrics on $\Sp(4)$ which are non-naturally reductive and $\Ad(4\Sp(1))$-invariant. For the general case $(\Sp(2n_1n_2),2\Sp(n_1n_2),2n_2\Sp(n_1))$, there exist $2n_2-1$ involutions of $\sp(2n_1n_2)$ such that $2n_2\sp(n_1))$ is the fixed point of these $2n_2-1$ involutions, and it follows the decomposition of $\sp(2n_1n_2)$ as a direct sum of irreducible $\ad(2n_2\sp(n_1))$-modules. In order to give new non-naturally reductive and $\Ad(2n_2\Sp(n_1)))$-invariant Einstein metrics on $\Sp(2n_1n_2)$, we prove a general result, i.e. $\Sp(2k+l)$ admits at least two non-naturally reductive Einstein metrics which are $\Ad(\Sp(k)\times\Sp(k)\times\Sp(l))$-invariant if $k<l$. It implies that every compact simple Lie group $\Sp(n)$ for $n\geq 4$ admits at least $2[\frac{n-1}{3}]$ non-naturally reductive left-invariant Einstein metrics.


title: Comment on "Minimal Basis Iterative Stockholder: Atoms in Molecules for Force-Field Development"

abstract: Verstraelen et al. (J. Chem. Theory Comput. 12 (2016) 3894-3912) recently introduced a new method for partitioning the electron density of a material into constituent atoms. Their approach falls within the class of atomic population analysis methods called stockholder charge partitioning methods in which a material electron distribution is divided into overlapping atoms. The Minimal Basis Iterative Stockholder (MBIS) method proposed by Verstraelen et al. composes the pro-atom density as a sum of exponential functions, where the number of exponential functions equals that element's row in the Periodic Table. Specifically, one exponential function is used for H and He, two for Li through Ne, three for Na through Ar, etc. In the MBIS method, the exponential functions parameters defining the pro-atom density are optimized in a self-consistent iterative procedure. Close examination reveals some important anomalies in the article by Verstaelen et al. The purpose of this comment article is to bring these important issues to readers' attention and to start a discussion of them.


title: Atomic vapor as a source of tunable, non-Gaussian self-reconstructing optical modes

abstract: In this manuscript, we demonstrate the ability of nonlinear light-atom interactions to produce tunably non-Gaussian, partially self-healing optical modes. Gaussian spatial-mode light tuned near to the atomic resonances in hot rubidium vapor is shown to result in non-Gaussian output mode structures that may be controlled by varying either the input beam power or the temperature of the atomic vapor. We show that the output modes exhibit a degree of self-reconstruction after encountering an obstruction in the beam path. The resultant modes are similar to truncated Bessel-Gauss modes that exhibit the ability to self-reconstruct earlier upon propagation than Gaussian modes. The ability to generate tunable, self-reconstructing beams has potential applications to a variety of imaging and communication scenarios.


title: Perturbative Expansion of Irreversible Work in Fokker-Planck Equation a la Quantum Mechanics

abstract: We discuss the systematic expansion of the solution of the Fokker-Planck equation with the help of the eigenfunctions of the time-dependent Fokker-Planck operator. The expansion parameter is the time derivative of the external parameter which controls the form of an external potential. Our expansion corresponds to the perturbative calculation of the adiabatic motion in quantum mechanics. With this method, we derive a new formula to calculate the irreversible work order by order, which is expressed as the expectation value with a pseudo density matrix. Applying this method to the case of the harmonic potential, we show that the first order term of the expansion gives the exact result. Because we do not need to solve the coupled differential equations of moments, our method simplifies the calculations of various functions such as the fluctuation of the irreversible work per unit time. We further investigate the exact optimized protocol to minimize the irreversible work by calculating its variation with respect to the control parameter itself.

categories: cond-mat.stat-mech math-ph math.MP nucl-th quant-ph


title: Towards an algebraic natural proofs barrier via polynomial identity testing

abstract: We observe that a certain kind of algebraic proof - which covers essentially all known algebraic circuit lower bounds to date - cannot be used to prove lower bounds against VP if and only if what we call succinct hitting sets exist for VP. This is analogous to the Razborov-Rudich natural proofs barrier in Boolean circuit complexity, in that we rule out a large class of lower bound techniques under a derandomization assumption. We also discuss connections between this algebraic natural proofs barrier, geometric complexity theory, and (algebraic) proof complexity.


title: Inelastic Neutrino-Nucleus Interactions within the Spectral Function Formalism

abstract: We report the results of a study of neutrino-carbon interactions at beam energies ranging between few hundreds MeV and few tens of GeV, carried out within the framework of the impulse approximation using a realistic spectral function. The contributions of quasi elastic scattering, resonance production and deep inelastic scattering---consistently obtained, for first time, from a model based on a realistic description of the nuclear ground state---are compared and analyzed.


title: Jupiter's formation and its primordial internal structure

abstract: The composition of Jupiter and the primordial distribution of the heavy elements are determined by its formation history. As a result, in order to constrain the primordial internal structure of Jupiter the growth of the core and the deposition and settling of accreted planetesimals must be followed in detail. In this paper we determine the distribution of the heavy elements in proto-Jupiter and determine the mass and composition of the core. We find that while the outer envelope of proto-Jupiter is typically convective and has an homogeneous composition, the innermost regions have compositional gradients. In addition, the existence of heavy elements in the envelope leads to much higher internal temperatures (several times 10$^4$ K) than in the case of a hydrogen-helium envelope. The derived core mass depends on the actual definition of the core: if the core is defined as the region in which the heavy-element mass fraction is above some limit (say 0.5), then it can be much more massive ($\sim$ 15 M$_{\oplus}$) and more extended (10\% of the planet's radius) than in the case where the core is just the region with 100\% heavy elements. In the former case Jupiter's core also consists of hydrogen and helium. Our results should be taken into account when constructing internal structure models of Jupiter and when interpreting the upcoming data from the {\it Juno} (NASA) mission.


title: Amoebas of curves and the Lyashko-Looijenga map

abstract: For any curve $\mathcal{V}$ in a toric surface $X$, we study the critical locus $S(\mathcal{V})$ of the moment map $\mu$ from $\mathcal{V}$ to its compactified amoeba $\mu(\mathcal{V})$. We show that for curves $\mathcal{V}$ in a fixed complete linear system, the critical locus $S(\mathcal{V})$ is smooth apart from some real codimension $1$ walls. We then investigate the topological classification of pairs $(\mathcal{V},S(\mathcal{V}))$ when $\mathcal{V}$ and $S(\mathcal{V})$ are smooth. As a main tool, we use the Lyashko-Looijenga mapping ($\mathcal{LL}$) relative to the logarithmic Gauss map $\gamma : \mathcal{V} \rightarrow \mathbb{C}P^1$. We prove two statements concerning $\mathcal{LL}$ that are crucial for our study: the map $\mathcal{LL}$ is algebraic; the map $\mathcal{LL}$ extends to nodal curves. It allows us to construct many examples of pairs $(\mathcal{V},S(\mathcal{V}))$ by perturbing nodal curves.


title: Fluctuation theorems for discrete kinetic models of molecular motors

abstract: Motivated by discrete kinetic models for non-cooperative molecular motors on periodic tracks, we consider random walks (also not Markov) on quasi one dimensional (1d) lattices, obtained by gluing several copies of a fundamental graph in a linear fashion. We show that, for a suitable class of quasi 1d lattices, the large deviation rate function associated to the position of the walker satisfies a Gallavotti-Cohen symmetry for any choice of the dynamical parameters defining the stochastic walk. This class includes the linear model considered in \cite{LLM1}. We also derive fluctuation theorems for the time-integrated cycle currents and discuss how the matrix approach of \cite{LLM1} can be extended to derive the above Gallavotti-Cohen symmetry for any Markov random walk on $\mathbb{Z}$ with periodic jump rates. Finally, we review in the present context some large deviation results of \cite{FS1} and give some specific examples with explicit computations.

categories: cond-mat.stat-mech math-ph math.MP math.PR


title: Follow the Compressed Leader: Faster Online Learning of Eigenvectors and Faster MMWU

abstract: The online problem of computing the top eigenvector is fundamental to machine learning. In both adversarial and stochastic settings, previous results (such as matrix multiplicative weight update, follow the regularized leader, follow the compressed leader, block power method) either achieve optimal regret but run slow, or run fast at the expense of loosing a $\sqrt{d}$ factor in total regret where $d$ is the matrix dimension. We propose a $\textit{follow-the-compressed-leader (FTCL)}$ framework which achieves optimal regret without sacrificing the running time. Our idea is to "compress" the matrix strategy to dimension 3 in the adversarial setting, or dimension 1 in the stochastic setting. These respectively resolve two open questions regarding the design of optimal and efficient algorithms for the online eigenvector problem.


title: In situ upgrade of quantum simulators to universal computers

abstract: Quantum simulators, machines that can replicate the dynamics of quantum systems, are being built as useful devices and are seen as a stepping stone to universal quantum computers. A key difference between the two is that computers have the ability to perform the logic gates that make up algorithms. We propose a method for learning how to construct these gates efficiently by using the simulator to perform optimal control on itself. This bypasses two major problems of purely classical approaches to the control problem: the need to have an accurate model of the system, and a classical computer more powerful than the quantum one to carry out the required simulations. Strong evidence that the scheme scales polynomially in the number of qubits, for systems of up to 9 qubits with Ising interactions, is presented from numerical simulations carried out in different topologies. This suggests that this in situ approach is a practical way of upgrading quantum simulators to computers.


title: DeepStack: Expert-Level Artificial Intelligence in No-Limit Poker

abstract: Artificial intelligence has seen several breakthroughs in recent years, with games often serving as milestones. A common feature of these games is that players have perfect information. Poker is the quintessential game of imperfect information, and a longstanding challenge problem in artificial intelligence. We introduce DeepStack, an algorithm for imperfect information settings. It combines recursive reasoning to handle information asymmetry, decomposition to focus computation on the relevant decision, and a form of intuition that is automatically learned from self-play using deep learning. In a study involving 44,000 hands of poker, DeepStack defeated with statistical significance professional poker players in heads-up no-limit Texas hold'em. The approach is theoretically sound and is shown to produce more difficult to exploit strategies than prior approaches.


title: Effective Potential for Cellular Size Control

abstract: For various species of biological cells, experimental observations indicate the existence of universal distributions of the cellular size, scaling relations between the cell-size moments and simple rules for the cell-size control. We address a class of models for the control of cell division, and present the steady state distributions. By introducing concepts such as effective force and potential, we are able to address the appearance of scaling collapse of different distributions and the connection between various moments of the cell-size. Our approach allows us to derive strict bounds which a potential cell-size control scenario must meet in order to yield a steady state distribution. The so-called "adder" model for cell-size control exhibits the weakest control that still enables the existence of stable size distribution, a fact that might explain the relative "popularity" of this scenario for different cells.


title: Online characterization of planetary surfaces: PlanetServer, an open-source analysis and visualization tool

abstract: The lack of open-source tools for hyperspectral data visualization and analysiscreates a demand for new tools. In this paper we present the new PlanetServer,a set of tools comprising a web Geographic Information System (GIS) and arecently developed Python Application Programming Interface (API) capableof visualizing and analyzing a wide variety of hyperspectral data from differentplanetary bodies. Current WebGIS open-source tools are evaluated in orderto give an overview and contextualize how PlanetServer can help in this mat-ters. The web client is thoroughly described as well as the datasets availablein PlanetServer. Also, the Python API is described and exposed the reason ofits development. Two different examples of mineral characterization of differenthydrosilicates such as chlorites, prehnites and kaolinites in the Nili Fossae areaon Mars are presented. As the obtained results show positive outcome in hyper-spectral analysis and visualization compared to previous literature, we suggestusing the PlanetServer approach for such investigations.


title: Open quantum generalisation of Hopfield neural networks

abstract: We propose a new framework to understand how quantum effects may impact on the dynamics of neural networks. We implement the dynamics of neural networks in terms of Markovian open quantum systems, which allows us to treat thermal and quantum coherent effects on the same footing. In particular, we propose an open quantum generalisation of the celebrated Hopfield neural network, the simplest toy model of associative memory. We determine its phase diagram and show that quantum fluctuations give rise to a qualitatively new non-equilibrium phase. This novel phase is characterised by limit cycles corresponding to high-dimensional stationary manifolds that may be regarded as a generalisation of storage patterns to the quantum domain.


title: Evidence for azimuthal variations of the oxygen abundance gradient tracing the spiral structure of the galaxy HCG91c

abstract: Context. The distribution of elements in galaxies forms an important diagnostic tool to characterize the system's formation and evolution. This tool is however complex to use in practice, as galaxies are subject to a range of simultaneous physical processes active from pc to kpc scales. This renders observations of the full optical extent of galaxies down to sub-kpc scales essential. Aims. Using the WiFeS integral field spectrograph, we previously detected abrupt and localized variations in the gas-phase oxygen abundance of the spiral galaxy HCG91c. Here, we follow-up on these observations to map HCG91c's disk out to ~2Re at a resolution of 600pc, and characterize the non-radial variations of the gas-phase oxygen abundance in the system. Methods. We obtained deep MUSE observations of the target under ~0.6 arcsec seeing conditions. We perform both a spaxel-based and aperture-based analysis of the data to map the spatial variations of 12+log(O/H) across the disk of the galaxy. Results. We confirm the presence of rapid variations of the oxygen abundance across the entire extent of the galaxy previously detected with WiFeS, for all azimuths and radii. The variations can be separated in two categories: a) localized and associated with individual HII regions, and b) extended over kpc scales, and occurring at the boundaries of the spiral structures in the galaxy. Conclusions. Our MUSE observations suggest that the enrichment of the interstellar medium in HGC91c has proceeded preferentially along spiral structures, and less efficiently across them. Our dataset highlights the importance of distinguishing individual star-forming regions down to scales of a few 100pc when using integral field spectrographs to spatially resolve the distribution of oxygen abundances in a given system, and accurately characterize azimuthal variations and intrinsic scatter.


title: SXP7.92: A Recently Rediscovered Be/X-ray Binary in the Small Magellanic Cloud, Viewed Edge On

abstract: We present a detailed optical and X-ray study of the 2013 outburst of the Small Magellanic Cloud Be/X-ray binary SXP7.92, as well as an overview of the last 18 years of observations from OGLE, RXTE, Chandra and XMM-Newton. We revise the position of this source to RA(J2000)=00:57:58.4, Dec(J2000)=-72:22:29.5 with a $1\sigma$ uncertainty of 1.5 arcsec, correcting the previously reported position by Coe et al. (2009) by more than 20 arcminutes. We identify and spectrally classify the correct counterpart as a B1Ve star. The optical spectrum is distinguished by an uncharacteristically deep narrow Balmer series, with the H$\alpha$ line in particular having a distinctive shell profile, i.e. a deep absorption core embedded in an emission line. We interpret this as evidence that we are viewing the system edge on and are seeing self obscuration of the circumstellar disc. We derive an optical period for the system of 40.0$\pm0.3$ days, which we interpret as the orbital period, and present several mechanisms to describe the X-ray/Optical behaviour in the recent outburst, in particular the "flares" and "dips" seen in the optical light curve, including a transient accretion disc and an elongated precessing disc.


title: Houdini for Astrophysical Visualization

abstract: The rapid growth in scale and complexity of both computational and observational astrophysics over the past decade necessitates efficient and intuitive methods for examining and visualizing large datasets. Here we discuss some newly developed tools to import and manipulate astrophysical data into the three dimensional visual effects software, {\it Houdini}. This software is widely used by visual effects artists, but a recently implemented Python API now allows astronomers to more easily use Houdini as a visualization tool. This paper includes a description of features, work flow, and various example visualizations. The project website, www.ytini.com, contains Houdini tutorials and links to the Python script Bitbucket repository aimed at a scientific audience to simplify the process of importing and rendering astrophysical data.


title: Direct evidence of hierarchical assembly at low masses from isolated dwarf galaxy groups

abstract: The demographics of dwarf galaxy populations have long been in tension with predictions from the Cold Dark Matter (CDM) paradigm. If primordial density fluctuations were scale-free as predicted, dwarf galaxies should themselves host dark matter subhaloes, the most massive of which may have undergone star formation resulting in dwarf galaxy groups. Ensembles of dwarf galaxies are observed as satellites of more massive galaxies, and there is observational and theoretical evidence to suggest that these satellites at z=0 were captured by the massive host halo as a group. However, the evolution of dwarf galaxies is highly susceptible to environment making these satellite groups imperfect probes of CDM in the low mass regime. We have identified one of the clearest examples to date of hierarchical structure formation at low masses: seven isolated, spectroscopically confirmed groups with only dwarf galaxies as members. Each group hosts 3-5 known members, has a baryonic mass of ~4.4 x 10^9 to 2 x 10^10 Msun, and requires a mass-to-light ratio of <100 to be gravitationally bound. Such groups are predicted to be rare theoretically and found to be rare observationally at the current epoch and thus provide a unique window into the possible formation mechanism of more massive, isolated galaxies.


title: Selection bias in dynamically-measured super-massive black hole samples: dynamical masses and dependence on S\'ersic index

abstract: We extend the comparison between the set of local galaxies having dynamically measured black holes with galaxies in the Sloan Digital Sky Survey (SDSS). We first show that the most up-to-date local black hole samples of early-type galaxies with measurements of effective radii, luminosities, and S\'ersic indices of the bulges of their host galaxies, have dynamical mass and S\'ersic index distributions consistent with those of SDSS early-type galaxies of similar bulge stellar mass. The host galaxies of local black hole samples thus do not appear structurally different from SDSS galaxies, sharing similar dynamical masses, light profiles and light distributions. Analysis of the residuals reveals that velocity dispersion is more fundamental than S\'ersic index n in the scaling relations between black holes and galaxies. Indeed, residuals with S\'ersic index could be ascribed to the (weak) correlation with bulge mass or even velocity dispersion. Finally, targetted Monte Carlo simulations that include the effects of the sphere of influence of the black hole, and tuned to reproduce the observed residuals and scaling relations in terms of velocity dispersion and stellar mass, show that, at least for galaxies with Mbulge > 1e10 Msun and n>5, the observed mean black hole mass at fixed S\'ersic index is biased significantly higher than the intrinsic value.


title: On the formation of hot and warm Jupiters via secular high-eccentricity migration in stellar triples

abstract: Hot Jupiters (HJs) are Jupiter-like planets orbiting their host star in tight orbits of a few days. They are commonly believed not to have formed in situ, requiring inwards migration towards the host star. One of the proposed migration scenarios is secular high-eccentricity or high-$e$ migration, in which the orbit of the planet is perturbed to high eccentricity by secular processes, triggering strong tidal evolution and orbital migration. Previous theoretical studies have considered secular excitation in stellar binaries. Recently, a number of HJs have been observed in stellar triple systems. In the latter, the secular dynamics are much richer compared to stellar binaries, and HJs could potentially be formed more efficiently. Here, we investigate this possibility by modeling the secular dynamical and tidal evolution of planets in two hierarchical configurations in stellar triple systems. We find that the HJ formation efficiency is higher compared to stellar binaries, but only by at most a few tens of per cent. The orbital properties of the HJs formed in the simulations are very similar to HJs formed in stellar binaries, and similarly to studies of the latter we find no significant number of warm Jupiters. HJs are only formed in our simulations for triples with specific orbital configurations, and our constraints are approximately consistent with current observations. In future, this allows to rule out high-$e$ migration in stellar triples if a HJ is detected in a triple grossly violating these constraints.


title: $^{13}$CO/C$^{18}$O Gradients Across the Disks of Nearby Spiral Galaxies

abstract: We use the IRAM Large Program EMPIRE and new high-resolution ALMA data to measure 13CO(1-0)/C18O(1-0) intensity ratios across nine nearby spiral galaxies. These isotopologues of CO are typically optically thin across most of the area in galaxy disks, and this ratio allows us to gauge their relative abundance due to chemistry or stellar nucleosynthesis effects. Resolved 13CO/C18O gradients across normal galaxies have been rare due to the faintness of these lines. We find a mean 13CO/C18O ratio of 6.0$\pm$0.9 for the central regions of our galaxies. This agrees well with results in the Milky Way, but differs from results for starburst galaxies (3.4$\pm$0.9) and ultraluminous infrared galaxies (1.1$\pm$0.4). In our sample, the 13CO/C18O ratio consistently increases with increasing galactocentric radius and decreases with increasing star formation rate surface density. These trends qualitatively agree with expectations for carbon and oxygen isotopic abundance variations due to stellar nucleosynthesis, with a possible effect of fractionation.