title: Electronic friction near metal surfaces: a case where molecule-metal couplings depend on nuclear coordinates

abstract: We derive an explicit form for the electronic friction as felt by a molecule near a metal surface for the general case that molecule-metal couplings depend on nuclear coordinates. Our work generalizes a previous study by von Oppen et al [Beilstein Journal of Nanotechnology, 3, 144, 2012], where we now go beyond the Condon approximation (i.e. molecule-metal couplings are not held constant). Using a non-equilibrium Green's function formalism in the adiabatic limit, we show that fluctuating metal-molecule couplings lead to new frictional damping terms and random forces, plus a correction to the potential of mean force. Numerical tests are performed and compared with a modified classical master equation; our results indicate that violating the Condon approximation can have a large effect on dynamics.


title: Deviations from universality in the fluctuation behavior of a heterogeneous complex system reveal intrinsic properties of components: The case of the international currency market

abstract: Identifying behavior that is relatively invariant under different conditions is a challenging task in far-from-equilibrium complex systems. As an example of how the existence of a semi-invariant signature can be masked by the heterogeneity in the properties of the components comprising such systems, we consider the exchange rate dynamics in the international currency market. We show that the exponents characterizing the heavy tails of fluctuation distributions for different currencies systematically diverge from a putative universal form associated with the median value (~2) of the exponents. We relate the degree of deviation of a particular currency from such an "inverse square law" to fundamental macroscopic properties of the corresponding economy, viz., measures of per capita production output and diversity of export products. We also show that in contrast to uncorrelated random walks exhibited by the exchange rate dynamics for currencies belonging to developed economies, those of the less developed economies show characteristics of sub-diffusive processes which we relate to the anti-correlated nature of the corresponding fluctuations. Approaches similar to that presented here may help in identifying invariant features obscured by the heterogeneous nature of components in other complex systems.


title: The big-bang theory: construction, evolution and status

abstract: Over the past century, rooted in the theory of general relativity, cosmology has developed a very successful physical model of the universe: the {\em big-bang model}. Its construction followed different stages to incorporate nuclear processes, the understanding of the matter present in the universe, a description of the early universe and of the large scale structure. This model has been confronted to a variety of observations that allow one to reconstruct its expansion history, its thermal history and the structuration of matter. Hence, what we refer to as the big-bang model today is radically different from what one may have had in mind a century ago. This construction changed our vision of the universe, both on observable scales and for the universe as a whole. It offers in particular physical models for the origins of the atomic nuclei, of matter and of the large scale structure. This text summarizes the main steps of the construction of the model, linking its main predictions to the observations that back them up. It also discusses its weaknesses, the open questions and problems, among which the need for a dark sector including dark matter and dark energy.


title: Generation and precise control of dynamic biochemical gradients for cellular assays

abstract: Spatial gradients of diffusible signalling molecules play crucial roles in controlling diverse cellular behaviour such as cell differentiation, tissue patterning and chemotaxis. In this paper, we report the design and testing of a microfluidic device for diffusion-based gradient generation for cellular assays. A unique channel design of the device eliminates cross-flow between the source and sink channels, thereby stabilising gradients by passive diffusion. The platform also enables quick and flexible control of chemical concentration that makes highly dynamic gradients in diffusion chambers. A model with the first approximation of diffusion and surface adsorption of molecules recapitulates the experimentally observed gradients. Budding yeast cells cultured in a gradient of a chemical inducer expressed a reporter fluorescence protein in a concentration-dependent manner. This microfluidic platform serves as a versatile prototype applicable to a broad range of biomedical investigations.


title: Planar digraphs of digirth four are 2-colourable

abstract: Neumann-Lara conjectured in 1985 that every planar digraph with digirth at least three is 2-colourable, meaning that the vertices can be 2-coloured without creating any monochromatic directed cycles. We prove a relaxed version of this conjecture: every planar digraph of digirth at least four is 2-colourable.


title: Diffusion paths between product life-cycles in European phonographic markets

abstract: We have investigated the product life-cycles of almost 17 000 hit singles performed on the 12 biggest national phonographic markets in Europe including: Austria, Belgium, France, Germany, Ireland, Italy, Netherlands, Norway, Spain, Sweden, Switzerland and the United Kingdom. We have considered weekly singles charts from the last 50 years (1966-2015) in each country. We analyze the spread of hit singles popularity (chart topping) as an epidemiological process performed on various European countries. Popular hit singles are contagious from one country to another. Thus, we consider time delays between the initial hit single release and reaching the highest position on consecutive national singles charts. We create directed network of countries representing transitions of hit singles popularity between countries. It is obtained by simulating the most likely paths and picking up the most frequent links. A country of initial hit single release is considered as a source of infection. Our algorithm builds up spanning trees by attaching new nodes. The probability of attachment depends on: 1) new nodes immunity 2) infectivity of previous nodes from the tree. Thus we obtain network of popularity spread with a hub the UK, a bridge the Netherlands and outliers Italy and Spain. We have found a characteristic topology of hit singles popularity spread. The positive correlation between this network and geographic or cultural grid-map of Europe is also observed. However, the network of popularity spread has some typical properties of complex networks.


title: Exploring Arsenic danger awareness in the Polish Copper Basin via network simulation - preliminary results

abstract: Information spread plays crucial role in risk management in case of environmental danger. The structure of local society may be well described by social network analysis. We have investigated the role of hubs within that concept. In case of danger, there are two different strategies of information management: 1) information spread that leads to awareness of the whole society 2) keeping the whole information in secret under control that leads to a partial social awareness, available to a small number of people only. In our model, the probability of information spread between two nodes is inverse proportional to connectivity of the next node because people who have a lot of connections are more immune. We have applied agent-based modelling on Barabasi Albert networks to explore various scenarios of information spread. We have considered recent arsenic environmental danger in G{\l}og\'ow and Legnica (Copper Basin) according to the official available data (2015 - 2016) from social network analysis point of view. We have considered various levels of environmental danger. Despite blocking of information by hubs, the successful information spread is possible when levels of danger are high enough. The perception and impression of information spread by society is also.


title: Multi-stage high order semi-Lagrangian schemes for incompressible flows in Cartesian geometries

abstract: Efficient transport algorithms are essential to the numerical resolution of incompressible fluid flow problems. Semi-Lagrangian methods are widely used in grid based methods to achieve this aim. The accuracy of the interpolation strategy then determines the properties of the scheme. We introduce a simple multi-stage procedure which can easily be used to increase the order of accuracy of a code based on multi-linear interpolations. This approach is an extension of a corrective algorithm introduced by Dupont \& Liu (2003, 2007). This multi-stage procedure can be easily implemented in existing parallel codes using a domain decomposition strategy, as the communications pattern is identical to that of the multi-linear scheme. We show how a combination of a forward and backward error correction can provide a third-order accurate scheme, thus significantly reducing diffusive effects while retaining a non-dispersive leading error term.


title: Strong spin Seebeck effect in Kondo T-shaped double quantum dots

abstract: We theoretically investigate the thermoelectric and spin thermoelectric properties of a T-shaped double quantum dot strongly coupled to two ferromagnetic leads, focusing on transport regime where the system exhibits the two-stage Kondo effect. We study the dependence of the (spin) Seebeck coefficient, the corresponding power factor and the figure of merit on temperature, leads' spin polarization and dot level position. We show that the thermal conductance fulfills a modified Wiedemann-Franz law. We also demonstrate that the spin thermopower is enhanced at temperatures corresponding to the second stage of Kondo screening. Very interestingly, the spin-thermoelectric response of the system is found to be highly sensitive to the spin polarization of the leads. In some cases spin polarization of the order of 1% is sufficient for a strong spin Seebeck effect to occur. This is explained as a consequence of the interplay between the two-stage Kondo effect and the exchange field induced in the double quantum dot. All calculations are performed with the aid of numerical renormalization group technique.


title: On the existence of non-hyperbolic ergodic measures as the limit of periodic measures

abstract: [GIKN] and [BBD1] propose two very different ways for building non hyperbolic measures, [GIKN] building such a measure as the limit of periodic measures and [BBD1] as the $\omega$-limit set of a single orbit, with a uniformly vanishing Lyapunov exponent. The technique in [GIKN] was essentially used in a generic setting, as the periodic orbits were built by small perturbations. It is not known if the measures obtained by the technique in [BBD1] are accumulated by periodic measures. In this paper we use a shadowing lemma from [G]: $\bullet$for getting the periodic orbits in [GIKN] without perturbing the dynamics, $\bullet$for recovering the compact set in [BBD1] with a uniformly vanishing Lyapunov exponent by considering the limit of periodic orbits. As a consequence, we prove that there exists an open and dense subset $\mathcal{U}$ of the set of robustly transitive non-hyperbolic diffeomorphisms far from homoclinic tangencies, such that for any $f\in\mathcal{U}$, there exists a non-hyperbolic ergodic measure with full support and approximated by hyperbolic periodic measures. We also prove that there exists an open and dense subset $\mathcal{V}$ of the set of diffeomorphisms exhibiting a robust cycle, such that for any $f\in\mathcal{V}$, there exists a non-hyperbolic ergodic measure approximated by hyperbolic periodic measures.


title: Connected Lie Groupoids are Internally Connected and Integral Complete in Synthetic Differential Geometry

abstract: We extend some fundamental definitions and constructions in the established generalisation of Lie theory involving Lie groupoids by reformulating them in terms of groupoids internal to a well-adapted model of synthetic differential geometry. In particular we define internal counterparts of the definitions of source path and source simply connected groupoid and the integration of $A$-paths. The main results of this paper show that if a classical Hausdorff Lie groupoid satisfies one of the classical connectedness conditions it also satisfies its internal counterpart.


title: Quantifying and Reducing Stereotypes in Word Embeddings

abstract: Machine learning algorithms are optimized to model statistical properties of the training data. If the input data reflects stereotypes and biases of the broader society, then the output of the learning algorithm also captures these stereotypes. In this paper, we initiate the study of gender stereotypes in {\em word embedding}, a popular framework to represent text data. As their use becomes increasingly common, applications can inadvertently amplify unwanted stereotypes. We show across multiple datasets that the embeddings contain significant gender stereotypes, especially with regard to professions. We created a novel gender analogy task and combined it with crowdsourcing to systematically quantify the gender bias in a given embedding. We developed an efficient algorithm that reduces gender stereotype using just a handful of training examples while preserving the useful geometric properties of the embedding. We evaluated our algorithm on several metrics. While we focus on male/female stereotypes, our framework may be applicable to other types of embedding biases.


title: Infinite Product Exponents for Modular Forms

abstract: Recently, D. Choi obtained a description of the coefficients of the infinite product expansions of meromorphic modular forms over $\Gamma_0(N)$. Using this result, we provide some bounds on these infinite product coefficients for holomorphic modular forms. We give an exponential upper bound for the growth of these coefficients. We show that this bound is also a lower bound in the case that the genus of the associated modular curve $X_0(N)$ is $0$ or $1$.


title: A Boltzmann Machine Implementation for the D-Wave

abstract: The D-Wave is an adiabatic quantum computer. It is an understatement to say that it is not a traditional computer. It can be viewed as a computational accelerator or more precisely a computational oracle, where one asks it a relevant question and it returns a useful answer. The question is how do you ask a relevant question and how do you use the answer it returns. This paper addresses these issues in a way that is pertinent to machine learning. A Boltzmann machine is implemented with the D-Wave since the D-Wave is merely a hardware instantiation of a partially connected Boltzmann machine. This paper presents a prototype implementation of a 3-layered neural network where the D-Wave is used as the middle (hidden) layer of the neural network. This paper also explains how the D-Wave can be utilized in a multi-layer neural network (more than 3 layers) and one in which each layer may be multiple times the size of the D-Wave being used.


title: Prospects for joint observations of gravitational waves and gamma rays from merging neutron star binaries

abstract: The detection of the events GW150914 and GW151226, both consistent with the merger of a binary black hole system (BBH), opened the era of gravitational wave (GW) astronomy. Besides BBHs, the most promising GW sources are the coalescences of binary systems formed by two neutron stars or a neutron star and a black hole. These mergers are thought to be connected with short Gamma Ray Bursts (GRBs), therefore combined observations of GW and electromagnetic (EM) signals could definitively probe this association. We present a detailed study on the expectations for joint GW and high-energy EM observations of coalescences of binary systems of neutron stars with Advanced Virgo and LIGO and with the \emph{Fermi} gamma-ray telescope. To this scope, we designed a dedicated Montecarlo simulation pipeline for the multimessenger emission and detection by GW and gamma-ray instruments, considering the evolution of the GW detector sensitivities. We show that the expected rate of joint detection is low during the Advanced Virgo and Advanced LIGO 2016-2017 run; however, as the interferometers approach their final design sensitivities, the rate will increase by $\sim$ a factor of ten. Future joint observations will help to constrain the association between short GRBs and binary systems and to solve the puzzle of the progenitors of GWs. Comparison of the joint detection rate with the ones predicted in this paper will help to constrain the geometry of the GRB jet.


title: Introducing a Calculus of Effects and Handlers for Natural Language Semantics

abstract: In compositional model-theoretic semantics, researchers assemble truth-conditions or other kinds of denotations using the lambda calculus. It was previously observed that the lambda terms and/or the denotations studied tend to follow the same pattern: they are instances of a monad. In this paper, we present an extension of the simply-typed lambda calculus that exploits this uniformity using the recently discovered technique of effect handlers. We prove that our calculus exhibits some of the key formal properties of the lambda calculus and we use it to construct a modular semantics for a small fragment that involves multiple distinct semantic phenomena.


title: Bootstrapping with Models: Confidence Intervals for Off-Policy Evaluation

abstract: For an autonomous agent, executing a poor policy may be costly or even dangerous. For such agents, it is desirable to determine confidence interval lower bounds on the performance of any given policy without executing said policy. Current methods for exact high confidence off-policy evaluation that use importance sampling require a substantial amount of data to achieve a tight lower bound. Existing model-based methods only address the problem in discrete state spaces. Since exact bounds are intractable for many domains we trade off strict guarantees of safety for more data-efficient approximate bounds. In this context, we propose two bootstrapping off-policy evaluation methods which use learned MDP transition models in order to estimate lower confidence bounds on policy performance with limited data in both continuous and discrete state spaces. Since direct use of a model may introduce bias, we derive a theoretical upper bound on model bias for when the model transition function is estimated with i.i.d. trajectories. This bound broadens our understanding of the conditions under which model-based methods have high bias. Finally, we empirically evaluate our proposed methods and analyze the settings in which different bootstrapping off-policy confidence interval methods succeed and fail.


title: Cutting out the middleman: measuring nuclear area in histopathology slides without segmentation

abstract: The size of nuclei in histological preparations from excised breast tumors is predictive of patient outcome (large nuclei indicate poor outcome). Pathologists take into account nuclear size when performing breast cancer grading. In addition, the mean nuclear area (MNA) has been shown to have independent prognostic value. The straightforward approach to measuring nuclear size is by performing nuclei segmentation. We hypothesize that given an image of a tumor region with known nuclei locations, the area of the individual nuclei and region statistics such as the MNA can be reliably computed directly from the image data by employing a machine learning model, without the intermediate step of nuclei segmentation. Towards this goal, we train a deep convolutional neural network model that is applied locally at each nucleus location, and can reliably measure the area of the individual nuclei and the MNA. Furthermore, we show how such an approach can be extended to perform combined nuclei detection and measurement, which is reminiscent of granulometry.


title: Slice-quaternionic Hopf surfaces

abstract: We investigate slice-quaternionic Hopf surfaces. In particular, we construct new structures of slice-quaternionic manifold on $\mathbb{S}^1\times\mathbb{S}^7$, we study their group of automorphisms and their deformations.


title: Desain dan Implementasi Face Recognition dan Live Streaming pada Sistem Digital Assistant untuk Staf Medik Fungsional menggunakan Google Glass

abstract: DIANE is a digital assistant system that aims to fasten the doctor access to various informations at the hospital such as health care facility, medical records, and also human resource data. The fasten access could be achieved by implementing face recognition and live streaming as part of the digital assistant system. By using face recognition, DIANE can identify patient's identity. Moreover, DIANE could know patient's biodata and also patient's medical records. Live streaming in DIANE enables doctor to communicate with other doctors in different place. Through the increasing speed of doctor to access the various informations at the hospital, it is expected that doctor can make decisions more quickly. As a result, more and more patients can be treated by the doctor properly. ---- DIANE merupakan sebuah sistem digital assistant yang bertujuan untuk mempercepat akses staf medik fungsional sebagai seorang profesional dengan waktu terbatas terhadap informasi seperti layanan kesehatan, rekam medis, dan personalia pada SIRS tingkat rumah sakit. Percepatan akses terhadap informasi tersebut diwujudkan melalui penggunaan fitur face recognition dan live streaming. Dengan menggunakan face recognition, DIANE dapat mengidentifikasi pasien sehingga dapat memberikan tampilan data berupa biodata dan rekam medis pasien. Sedangkan menggunakan live streaming pada DIANE, seorang staf medik fungsional dapat berkomunikasi dengan staf medik fungsional di tempat lain. Dengan peningkatan kecepatan akses terhadap beragam informasi tersebut diharapkan seorang staf medik fungsional dapat mengambil keputusan dengan lebih cepat. Akibatnya, semakin banyak pasien yang dapat ditangani oleh staf medik fungsional dengan baik.


title: A new characterization of the jump rate for piecewise-deterministic Markov processes with discrete transitions

abstract: Piecewise-deterministic Markov processes form a general class of non-diffusion stochastic models that involve both deterministic trajectories and random jumps at random times. In this paper, we state a new characterization of the jump rate of such a process with discrete transitions. We deduce from this result a nonparametric technique for estimating this feature of interest. We state the uniform convergence in probability of the estimator. The methodology is illustrated on a numerical example.


title: Quantum processing by remote quantum control

abstract: Client-server models enable computations to be hosted remotely on quantum servers. We present a novel protocol for realizing this task, with practical advantages when using technology feasible in the near term. Client tasks are realized as linear combinations of operations implemented by the server, where the linear coefficients are hidden from the server. We report on an experimental demonstration of our protocol using linear optics, which realizes linear combination of two single-qubit operations by a remote single-qubit control. In addition, we explain when our protocol can remain efficient for larger computations, as well as some ways in which privacy can be maintained using our protocol.


title: Spectroscopy at the solar limb: II. Are spicules heated to coronal temperatures ?

abstract: Spicules of the so-called type II were suggested to be relevant for coronal heating because of their ubiquity on the solar surface and their eventual extension into the corona. We investigate whether solar spicules are heated to transition-region or coronal temperatures and reach coronal heights (>6 Mm) using multi-wavelength observations of limb spicules in different chromospheric spectral lines (Ca II H, Hepsilon, Halpha, Ca II IR at 854.2 nm, He I at 1083 nm). We determine the line width of individual spicules and throughout the field of view and estimate the maximal height that different types of off-limb features reach. We derive estimates of the kinetic temperature and the non-thermal velocity from the line width of spectral lines from different chemical elements. We find that most regular spicules reach a maximal height of about 6 Mm above the solar limb. The majority of features found at larger heights are irregularly shaped with a significantly larger lateral extension than spicules. Both individual and average line profiles in all spectral lines show a decrease in their line width with height above the limb with very few exceptions. Both the kinetic temperature and the non-thermal velocity decrease with height above the limb. We find no indications that the spicules in our data reach coronal heights or transition-region or coronal temperatures.


title: A Study of Energy and Locality Effects using Space-filling Curves

abstract: The cost of energy is becoming an increasingly important driver for the operating cost of HPC systems, adding yet another facet to the challenge of producing efficient code. In this paper, we investigate the energy implications of trading computation for locality using Hilbert and Morton space-filling curves with dense matrix-matrix multiplication. The advantage of these curves is that they exhibit an inherent tiling effect without requiring specific architecture tuning. By accessing the matrices in the order determined by the space-filling curves, we can trade computation for locality. The index computation overhead of the Morton curve is found to be balanced against its locality and energy efficiency, while the overhead of the Hilbert curve outweighs its improvements on our test system.


title: Single-diffractive Drell-Yan pair production at the LHC

abstract: We present predictions for single-diffractive low-mass Drell-Yan pair production in $pp$ collisions at the LHC at $\sqrt{s}=13$ TeV. Predictions are obtained adopting a factorised form for the relevant cross sections and are based on a new set of diffractive parton distributions resulting from the QCD analysis of combined HERA leading proton data. We discuss a number of observables useful to characterise the expected factorisation breaking effects.


title: The Minimum Cost Connected Subgraph Problem in Medical Image Analysis

abstract: Several important tasks in medical image analysis can be stated in the form of an optimization problem whose feasible solutions are connected subgraphs. Examples include the reconstruction of neural or vascular structures under connectedness constraints. We discuss the minimum cost connected subgraph (MCCS) problem and its approximations from the perspective of medical applications. We propose a) objective-dependent constraints and b) novel constraint generation schemes to solve this optimization problem exactly by means of a branch-and-cut algorithm. These are shown to improve scalability and allow us to solve instances of two medical benchmark datasets to optimality for the first time. This enables us to perform a quantitative comparison between exact and approximative algorithms, where we identify the geodesic tree algorithm as an excellent alternative to exact inference on the examined datasets.


title: Desain dan Implementasi Sistem Digital Assistant Berbasis Google Glass pada Rumah Sakit

abstract: In order to improve the performance of the hospital, one solution that is needed is a system that can facilitate functional medical staff in the discharge of their duties. This paper discusses the design and implementation of digital assistant system designed for functional medical staff at the hospital. This system is implemented in Google Glass, PC, tablet, and smart phones. Google Glass device is equipped with two modules, which are patient data module with face recognition feature and communication module with livestreaming feature. Tablet devices and mobile phones are equipped with the personnel module and display medical records module. Tablet devices are also equipped with prescription writing. The whole process in all four devices require webservice to be executed. The implementation results indicate that the four software has been built based on the design and integrated with HIS (Hospital Information System) database via web service. ---- Demi meningkatkan kinerja rumah sakit, salah satu solusi yang diperlukan adalah sebuah sistem yang dapat mempermudah staf medik fungsional dalam menunaikan tugasnya. Makalah ini membahas mengenai desain dan implementasi sistem digital assistant yang dirancang untuk staf medik fungsional di rumah sakit. Sistem ini diimplementasikan pada perangkat keras Google Glass, PC, tablet, dan ponsel pintar. Perangkat Google Glass dilengkapi dengan dua fitur yaitu face recognition dan live streaming. Perangkat tablet dan ponsel dilengkapi dengan modul personalia dan modul tampilan rekam medis. Perangkat tablet juga dilengkapi dengan fitur penulisan resep. Seluruh proses pada keempat perangkat tersebut membutuhkan web service untuk dapat dieksekusi. Hasil implementasi mengindikasikan bahwa keempat perangkat lunak telah sesuai dengan rancangan dan terintegrasikan dengan database SIRS tingkat rumah sakit melalui web service.


title: LSTM-Based Predictions for Proactive Information Retrieval

abstract: We describe a method for proactive information retrieval targeted at retrieving relevant information during a writing task. In our method, the current task and the needs of the user are estimated, and the potential next steps are unobtrusively predicted based on the user's past actions. We focus on the task of writing, in which the user is coalescing previously collected information into a text. Our proactive system automatically recommends the user relevant background information. The proposed system incorporates text input prediction using a long short-term memory (LSTM) network. We present simulations, which show that the system is able to reach higher precision values in an exploratory search setting compared to both a baseline and a comparison system.


title: A geometric inequality for convex free boundary hypersurfaces in the unit ball

abstract: We use the inverse mean curvature flow with a free boundary perpendicular to the sphere to prove a geometric inequality involving the Willmore energy for convex hypersurfaces of dimension $n\geq 3$ with boundary on the sphere.


title: Indicable groups and $p_c<1$

abstract: A conjecture of Benjamini & Schramm from 1996 states that any finitely generated group that is not a finite extension of Z has a non-trivial percolation phase. Our main results prove this conjecture for certain groups, and in particular prove that any group with a non-trivial homomorphism into the additive group of real numbers satisfies the conjecture. We use this to reduce the conjecture to the case of hereditary just-infinite groups. The novelty here is mainly in the methods used, combining the methods of EIT and evolving sets, and using the algebraic properties of the group to apply these methods.


title: Rating Effects on Social News Posts and Comments

abstract: At a time when information seekers first turn to digital sources for news and opinion, it is critical that we understand the role that social media plays in human behavior. This is especially true when information consumers also act as information producers and editors through their online activity. In order to better understand the effects that editorial ratings have on online human behavior, we report the results of a two large-scale in-vivo experiments in social media. We find that small, random rating manipulations on social media posts and comments created significant changes in downstream ratings resulting in significantly different final outcomes. We found positive herding effects for positive treatments on posts, increasing the final rating by 11.02% on average, but not for positive treatments on comments. Contrary to the results of related work, we found negative herding effects for negative treatments on posts and comments, decreasing the final ratings on average, of posts by 5.15% and of comments by 37.4%. Compared to the control group, the probability of reaching a high rating (>=2000) for posts is increased by 24.6% when posts receive the positive treatment and for comments is decreased by 46.6% when comments receive the negative treatment.


title: Doubly coupled matter fields in massive bigravity

abstract: In the context of massive (bi-)gravity non-minimal matter couplings have been proposed. These couplings are special in the sense that they are free of the Boulware-Deser ghost below the strong coupling scale and can be used consistently as an effective field theory. Furthermore, they enrich the phenomenology of massive gravity. We consider these couplings in the framework of bimetric gravity and study the cosmological implications for background and linear tensor, vector, and scalar perturbations. Previous works have investigated special branch of solutions. Here we perform a complete perturbation analysis for the general background equations of motion completing previous analysis.


title: Comparing the hierarchy of keywords in on-line news portals

abstract: The tagging of on-line content with informative keywords is a widespread phenomenon from scientific article repositories through blogs to on-line news portals. In most of the cases, the tags on a given item are free words chosen by the authors independently. Therefore, relations among keywords in a collection of news items is unknown. However, in most cases the topics and concepts described by these keywords are forming a latent hierarchy, with the more general topics and categories at the top, and more specialised ones at the bottom. Here we apply a recent, cooccurrence-based tag hierarchy extraction method to sets of keywords obtained from four different on-line news portals. The resulting hierarchies show substantial differences not just in the topics rendered as important (being at the top of the hierarchy) or of less interest (categorised low in the hierarchy), but also in the underlying network structure. This reveals discrepancies between the plausible keyword association frameworks in the studied news portals.


title: Vibrato and automatic differentiation for high order derivatives and sensitivities of financial options

abstract: This paper deals with the computation of second or higher order greeks of financial securities. It combines two methods, Vibrato and automatic differentiation and compares with other methods. We show that this combined technique is faster than standard finite difference, more stable than automatic differentiation of second order derivatives and more general than Malliavin Calculus. We present a generic framework to compute any greeks and present several applications on different types of financial contracts: European and American options, multidimensional Basket Call and stochastic volatility models such as Heston's model. We give also an algorithm to compute derivatives for the Longstaff-Schwartz Monte Carlo method for American options. We also extend automatic differentiation for second order derivatives of options with non-twice differentiable payoff. 1. Introduction. Due to BASEL III regulations, banks are requested to evaluate the sensitivities of their portfolios every day (risk assessment). Some of these portfolios are huge and sensitivities are time consuming to compute accurately. Faced with the problem of building a software for this task and distrusting automatic differentiation for non-differentiable functions, we turned to an idea developed by Mike Giles called Vibrato. Vibrato at core is a differentiation of a combination of likelihood ratio method and pathwise evaluation. In Giles [12], [13], it is shown that the computing time, stability and precision are enhanced compared with numerical differentiation of the full Monte Carlo path. In many cases, double sensitivities, i.e. second derivatives with respect to parameters, are needed (e.g. gamma hedging). Finite difference approximation of sensitivities is a very simple method but its precision is hard to control because it relies on the appropriate choice of the increment. Automatic differentiation of computer programs bypass the difficulty and its computing cost is similar to finite difference, if not cheaper. But in finance the payoff is never twice differentiable and so generalized derivatives have to be used requiring approximations of Dirac functions of which the precision is also doubtful. The purpose of this paper is to investigate the feasibility of Vibrato for second and higher derivatives. We will first compare Vibrato applied twice with the analytic differentiation of Vibrato and show that it is equivalent, as the second is easier we propose the best compromise for second derivatives: Automatic Differentiation of Vibrato. In [8], Capriotti has recently investigated the coupling of different mathematical methods -- namely pathwise and likelihood ratio methods -- with an Automatic differ


title: On the elliptic $\mathfrak{gl}_2$ solid-on-solid model: functional relations and determinants

abstract: In this work we study an elliptic solid-on-solid model with domain-wall boundaries having the elliptic quantum group $\mathcal{E}_{p, \gamma}[\widehat{\mathfrak{gl}_2}]$ as its underlying symmetry algebra. We elaborate on results previously presented by the author and extend our analysis to include continuous families of single determinantal representations for the model's partition function. Interestingly, our families of representations are parameterized by two continuous complex variables which can be arbitrarily chosen without affecting the partition function.

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


title: Vibrationally resolved UV/Vis spectroscopy with time-dependent density functional based tight binding

abstract: We report a time-dependent density functional based tight-binding (TD-DFTB) scheme for the calculation of UV/Vis spectra, explicitly taking into account the excitation of nuclear vibrations via the adiabatic Hessian Franck-Condon (AH|FC) method with a harmonic approximation for the nuclear wavefunction. The theory of vibrationally resolved UV/Vis spectroscopy is first summarized from the viewpoint of TD-DFTB. The method is benchmarked against time-dependent density functional theory (TD-DFT) calculations for strongly dipole allowed excitations in various aromatic and polar molecules. Using the recent 3ob:freq parameter set of Elstner's group, very good agreement with TD-DFT calculations using local functionals was achieved.


title: Modeling the distribution of distance data in Euclidean space

abstract: Phylogenetic inference-the derivation of a hypothesis for the common evolutionary history of a group of species- is an active area of research at the intersection of biology, computer science, mathematics, and statistics. One assumes the data contains a phylogenetic signal that will be recovered with varying accuracy due to the quality of the method used, and the quality of the data. The input for distance-based inference methods is an element of a Euclidean space with coordinates indexed by the pairs of organisms. For several algorithms there exists a subdivision of this space into polyhedral cones such that inputs in the same cone return the same tree topology. The geometry of these cones has been used to analyze the inference algorithms. In this chapter, we model how input data points drawn from DNA sequences are distributed throughout Euclidean space in relation to the space of tree metrics, which in turn can also be described as a collection of polyhedral cones.


title: Lorentz factor distribution of blazars from the optical Fundamental plane of black hole activity

abstract: Blazar radiation is dominated by a relativistic jet which can be modeled at first approximation using just two intrinsic parameters - the Lorentz factor $\Gamma$ and the viewing angle $\theta$. Blazar jet observations are often beamed due to relativistic effects, complicating the understanding of these intrinsic properties. The most common way to estimate blazar Lorentz factors needs the estimation of apparent jet speeds and Doppler beaming factors. We present a new and independent method of constructing the blazar Lorentz factor distribution, using the optical fundamental plane of black hole activity. The optical fundamental plane is a plane stretched out by both the supermassive black holes and the X-ray binaries, in the 3D space provided by their [OIII] line luminosity, radio luminosity and black hole mass. We use the intrinsic radio luminosity obtained from the optical fundamental plane to constrain the boosting parameters of the VLBA Imaging and Polarimetry Survey (VIPS) blazar sample. We find a blazar bulk Lorentz factor distribution in the form of a power law as $N(\Gamma) \propto \Gamma^{-2.1 \pm 0.4}$ for the $\Gamma$ range of 1 to 40. We also discuss the viewing angle distribution of the blazars and the dependence of our results on the input parameters.


title: LHCb pentaquarks in constituent quark models

abstract: The recently discovered $P_c(4380)^+$ and $P_c(4450)^+$ states at LHCb have masses close to the $\bar D\Sigma_c^*$ and $\bar D^*\Sigma_c$ thresholds, respectively, which suggest that they may have significant meson-baryon molecular components. We analyze these states in the framework of a constituent quark model which has been applied to a wide range of hadronic observables, being the model parameters, therefore, completely constrained. The $P_c(4380)^+$ and $P_c(4450)^+$ are studied as molecular states composed by charmed baryons and open charm mesons. Several bound states with the proper binding energy are found in the $\bar D\Sigma_c^*$ and $\bar D^*\Sigma_c$ channels. We discuss the possible assignments of these states from their decay widths. Moreover, two more states are predicted, associated with the $\bar D\Sigma_c$ and $\bar D^* \Sigma_c^*$ thresholds.


title: On an inequality of different metrics for algebraic polynomials

abstract: We establish an inequality of different metrics for algebraic polynomials.


title: Minimal quartet-metric gravity: beyond LCDM per scalar graviton as the unified dark matter and dark energy

abstract: In the framework of the minimal quartet-metric gravity/systogravity, a scalar graviton/systolon is stated as a universal dark component, with supplementary manifestations in the different contexts either as dark matter or dark energy. An ensuing extension to the standard {\Lambda}CDM model is developed. A modification of the late expansion of the Universe, with an attractor of a scalar master equation defining an effective cosmological constant, which supersedes the true one, is proposed. A new partial solution to the cosmological constant problem is discussed.


title: Classification of 8-dimensional rank two commutative semifields

abstract: We classify the rank two commutative semifields which are 8-dimensional over their center $\mathbb{F}_{q}$. This is done using computational methods utilizing the connection to linear sets in $\mathrm{PG}(2,q^{4})$. We then apply our methods to complete the classification of rank two commutative semifields which are 10-dimensional over $\mathbb{F}_{3}$. The implications of these results are detailed for other geometric structures such as semifield flocks, ovoids of parabolic quadrics, and eggs.


title: A Note on Efficiency of Downsampling and Color Transformation in Image Quality Assessment

abstract: Several existing and successful full reference image quality assessment (IQA) models use linear color transformation and downsampling before measuring similarity or quality of images. This paper indicates to the right order of these two procedures and that the existing models have not chosen the more efficient approach. In addition, efficiency of these metrics is not compared in a fair basis in the literature.


title: Classical analogue of the continuous transition between the Weisskopf-Wigner exponential decay and the Rabi oscillation

abstract: When a discrete state is coupled to a continuum, the dynamics can be described either by the Weisskopf-Wigner exponential decay or by the Rabi oscillation, depending on the relative magnitudes of the continuum width and of the Rabi frequency. A continuous transition between these two regimes exists, as demonstrated in 1977 by C. Cohen-Tannoudji and P. Avan. Here, we describe a fully analogous transition in classical mechanics, by studying the dynamics of two coupled mechanical oscillators in the presence of damping. By varying the relative magnitudes of the damping and coupling terms, we observe a continuous transition between a regime analogous to the Rabi oscillation and a regime analogous to the Weisskopf-Wigner exponential decay.


title: Closed Form Fractional Integration and Differentiation via Real Exponentially Spaced Pole-Zero Pairs

abstract: We derive closed-form expressions for the poles and zeros of approximate fractional integrator/differentiator filters, which correspond to spectral roll-off filters having any desired log-log slope to a controllable degree of accuracy over any bandwidth. The filters can be described as a uniform exponential distribution of poles along the negative-real axis of the s plane, with zeros interleaving them. Arbitrary spectral slopes are obtained by sliding the array of zeros relative to the array of poles, where each array maintains periodic spacing on a log scale. The nature of the slope approximation is close to Chebyshev optimal in the interior of the pole-zero array, approaching conjectured Chebyshev optimality over all frequencies in the limit as the order approaches infinity. Practical designs can arbitrarily approach the equal-ripple approximation by enlarging the pole-zero array band beyond the desired frequency band. The spectral roll-off slope can be robustly modulated in real time by varying only the zeros controlled by one slope parameter. Software implementations are provided in matlab and Faust.


title: Certain identities on derivatives of radial homogeneous and logarithmic functions

abstract: Let $k$ be a natural number and $s$ be real. In the 1-dimensional case, the $k$-th order derivatives of the functions $\lvert x\rvert^s$ and $\log \lvert x\rvert$ are multiples of $\lvert x\rvert^{s-k}$ and $\lvert x\rvert^{-k}$, respectively. In the present paper, we generalize this fact to higher dimensions by introducing a suitable norm of the derivatives, and give the exact values of the multiples.


title: Invariance in Quantum Walks

abstract: In this Chapter, we present some interesting properties of quantum walks on the line. We concentrate our attention in the emergence of invariance and provide some insights into the ultimate origin of the observed behavior. In the first part of the Chapter, we review the building blocks of the quantum-mechanical version of the standard random walk in one dimension. The most distinctive difference between random and quantum walks is the replacement of the random coin in the former by the action of a unitary operator upon some internal property of the later. We provide explicit expressions for the solution to the problem when the most general form for the homogeneous unitary operator is considered, and we analyze several key features of the system as the presence of symmetries or stationary limits. After that, we analyze the consequences of letting the properties of the coin operator change from site to site, and from time step to time step. In spite of this lack of homogeneity, the probabilistic properties of the motion of the walker can remain unaltered if the coin variability is chosen adequately. Finally, we show how this invariance can be connected to the gauge freedom of electromagnetism.


title: Linear and nonlinear fractional Voigt models

abstract: We consider fractional generalizations of the ordinary differential equation that governs the creep phenomenon. Precisely, two Caputo fractional Voigt models are considered: a rheological linear model and a nonlinear one. In the linear case, an explicit Volterra representation of the solution is found, involving the generalized Mittag-Leffler function in the kernel. For the nonlinear fractional Voigt model, an existence result is obtained through a fixed point theorem. A nonlinear example, illustrating the obtained existence result, is given.


title: New formulas for the spectral radius via Aluthge transform

abstract: In this paper we give several expressions of spectral radius of a bounded operator on a Hilbert space, in terms of iterates of Aluthge transformation, numerical radius and the asymptotic behavior of the powers of this operator. Also we obtain several characterizations of normaloid operators.


title: BiFold visualization of bipartite datasets

abstract: The emerging domain of data-enabled science necessitates development of algorithms and tools for knowledge discovery. Human interaction with data through well-constructed graphical representation can take special advantage of our visual ability to identify patterns. We develop a data visualization framework, called BiFold, for exploratory analysis of bipartite datasets that describe binary relationships between groups of objects. Typical data examples would include voting records, organizational memberships, and pairwise associations, or other binary datasets. BiFold provides a low dimensional embedding of data that represents similarity by visual nearness, analogous to Multidimensional Scaling (MDS). The unique and new feature of BiFold is its ability to simultaneously capture both within-group and between-group relationships among objects, enhancing knowledge discovery. We benchmark BiFold using the {\it Southern Women Dataset}, where social groups are now visually evident. We construct BiFold plots for two US voting datasets: For the presidential election outcomes since 1976, BiFold illustrates the evolving geopolitical structures that underlie these election results. For Senate congressional voting, BiFold identifies a partisan coordinate, separating senators into two parties while simultaneously visualizing a bipartisan-coalition coordinate which captures the ultimate fate of the bills (pass/fail). Finally, we consider a global cuisine dataset of the association between recipes and food ingredients. BiFold allows us to visually compare and contrast cuisines while also allowing identification of signature ingredients of individual cuisines.


title: DoReFa-Net: Training Low Bitwidth Convolutional Neural Networks with Low Bitwidth Gradients

abstract: We propose DoReFa-Net, a method to train convolutional neural networks that have low bitwidth weights and activations using low bitwidth parameter gradients. In particular, during backward pass, parameter gradients are stochastically quantized to low bitwidth numbers before being propagated to convolutional layers. As convolutions during forward/backward passes can now operate on low bitwidth weights and activations/gradients respectively, DoReFa-Net can use bit convolution kernels to accelerate both training and inference. Moreover, as bit convolutions can be efficiently implemented on CPU, FPGA, ASIC and GPU, DoReFa-Net opens the way to accelerate training of low bitwidth neural network on these hardware. Our experiments on SVHN and ImageNet datasets prove that DoReFa-Net can achieve comparable prediction accuracy as 32-bit counterparts. For example, a DoReFa-Net derived from AlexNet that has 1-bit weights, 2-bit activations, can be trained from scratch using 6-bit gradients to get 46.1\% top-1 accuracy on ImageNet validation set. The DoReFa-Net AlexNet model is released publicly.


title: Jordan product commuting maps with $\lambda$-Aluthge transform

abstract: Let H and K be two complex Hilbert spaces and B(H) be the algebra of bounded linear operators from H into itself. The main purpose in this paper is to obtain a characterization of bijective maps $\Phi$ : B(H) $\rightarrow$ B(K) satisfying the following condition $\Delta$ $\lambda$ ($\Phi$(A) $\bullet$ $\Phi$(B)) = $\Phi$($\Delta$ $\lambda$ (A $\bullet$ B)) for all A, B $\in$ B(H), where $\Delta$ $\lambda$ (T) stands the $\lambda$-Aluthge transform of the operator T $\in$ B(H) and A $\bullet$ B = 1 2 (AB + BA) is the Jordan product of A and B. We prove that a bijective map $\Phi$ satisfies the above condition, if and only if there exists an unitary operator U : H $\rightarrow$ K, such that $\Phi$ has the form $\Phi$(A) = UAU * for all A $\in$ B(H).


title: MetaZipf. (Re)producing knowledge about city size distributions

abstract: Zipf's law for cities is probably the most famous regularity in social sciences. So much that, a hundred years of publication later, its status is not clear: is it a law of social organisation? Is it an instrument of description of city size distributions? Is it an element of validation of geographical objects (cities and systems of cities in particular)? Empirical estimations of the rank-size parameters are very numerous and contradict each other. In this study, we present the results of a reproducible meta-analysis of the largest pool of papers regarding this issue, obtained from the collection of data made open and the construction of an online interactive application which allows the reader to explore this literature. We find that a large part of the variations observed in the measure of Zipf's coefficient is unnecessary as it comes from the choice of different technical specifications in the way cities are defined, whereas some of the current theories to explain the remaining share of variations are challenged by our results.


title: Decoherence and disentanglement of qubits detecting scalar fields in an expanded universe

abstract: We consider Unruh-Wald qubit detector model adopted for the far future region of an exactly solvable 1+1 dimensional scalar field theory in a Robertson-Walker expanding spacetime. It is shown that the expansion of the universe in its history enhances the decoherence of the qubit coupled with a scalar field. Moreover, we consider two entangled qubits, each locally coupled a scalar field. The expansion of the universe in its history degrades the entanglement between the qubits, and can lead to entanglement sudden death if the initial entanglement is small enough. The details depend on the parameters characterizing the expansion of the universe. This work, albeit on a toy model, suggests that the history of the universe might be probed through the coherent and entanglement behavior of future detectors of quantum fields.


title: Pragmatic factors in image description: the case of negations

abstract: We provide a qualitative analysis of the descriptions containing negations (no, not, n't, nobody, etc) in the Flickr30K corpus, and a categorization of negation uses. Based on this analysis, we provide a set of requirements that an image description system should have in order to generate negation sentences. As a pilot experiment, we used our categorization to manually annotate sentences containing negations in the Flickr30K corpus, with an agreement score of K=0.67. With this paper, we hope to open up a broader discussion of subjective language in image descriptions.


title: Product commuting maps with the $\lambda$-Aluthge transform

abstract: Let H and K be two Hilbert spaces and B(H) be the algebra of all bounded linear operators from H into itself. The main purpose of this paper is to obtain a characterization of bijective maps $\Phi$ : B(H) $\rightarrow$ B(K) satisfying the following condition $\Delta$ $\lambda$ ($\Phi$(A)$\Phi$(B)) = $\Phi$($\Delta$ $\lambda$ (AB)) f orall A, B $\in$ B(H), where $\Delta$ $\lambda$ (T) stands the $\lambda$-Aluthge transform of the operator T $\in$ B(H). More precisely, we prove that a bijective map $\Phi$ satisfies the above condition, if and only , if $\Phi$(A) = U AU * for all A $\in$ B(H), for some unitary operator U : H $\rightarrow$ K.


title: Scalar correlator, Higgs decay into quarks, and scheme variations of the QCD coupling

abstract: In this work, the perturbative QCD series of the scalar correlation function $\Psi(s)$ is investigated. Besides ${\rm Im}\Psi(s)$, which is relevant for Higgs decay into quarks, two other physical correlators, $\Psi^{"}(s)$ and $D^L(s)$, have been employed in QCD applications like quark mass determinations or hadronic $\tau$ decays. $D^L(s)$ suffers from large higher-order corrections and, by resorting to the large-$\beta_0$ approximation, it is shown that this is related to a spurious renormalon ambiguity at $u=1$. Hence, this correlator should be avoided in phenomenological analyses. Moreover, it turns out advantageous to express the quark mass factor, introduced to make the scalar current renormalisation group invariant, in terms of the renormalisation invariant quark mass $\hat m_q$. To further study the behaviour of the perturbative expansion, we introduce a QCD coupling $\hat\alpha_s$, whose running is explicitly renormalisation scheme independent. The scheme dependence of $\hat\alpha_s$ is parametrised by a single parameter $C$, being related to transformations of the QCD scale parameter $\Lambda$. It is demonstrated that appropriate choices of $C$ lead to a substantial improvement in the behaviour of the perturbative series for $\Psi^{"}(s)$ and ${\rm Im}\Psi(s)$.


title: Simulating nonlinear cosmological structure formation with massive neutrinos

abstract: We present a new method for simulating cosmologies that contain massive particles with thermal free streaming motion, such as massive neutrinos or warm/hot dark matter. This method combines particle and fluid descriptions of the thermal species to eliminate the shot noise known to plague conventional N-body simulations. We describe this method in detail, along with results for a number of test cases to validate our method, and check its range of applicability. Using this method, we demonstrate that massive neutrinos can produce a significant scale-dependence in the large-scale biasing of deep voids in the matter field. We show that this scale-dependence may be quantitatively understood using an extremely simple spherical expansion model which reproduces the behavior of the void bias for different neutrino parameters.


title: Cobordism categories and moduli spaces of odd dimensional manifolds

abstract: We prove that the stable moduli space of $(n-1)$-connected, $n$-parallelizable, $(2n+1)$-dimensional manifolds is homology equivalent to an infinite loopspace for $n \geq 4, n \neq 7$. The main novel ingredient is a version of the cobordism category incorporating surgery data in the form of Lagrangian subspaces.


title: Variable time flow as an alternative to dark energy

abstract: Time is a parameter playing a central role in our most fundamental modelling of natural laws. Relativity theory shows that the comparison of times measured by different clocks depends on their relative motions and on the strength of the gravitational field in which they are embedded. In standard cosmology, the time parameter is the one measured by fundamental clocks, i.e. clocks at rest with respect to the expanding space. This proper time is assumed to flow at a constant rate throughout the whole history of the Universe. We make the alternative hypothesis that the rate at which cosmological time flows depends on the global geometric curvature the Universe. Using a simple one-parameter model for the relation between proper time and curvature, we build a cosmological model that fits the Type Ia Supernovae data (the best cosmological standard candles) without the need for dark energy nor probably exotic dark matter.


title: Entanglement of superconducting qubits via acceleration radiation

abstract: We show that simulated relativistic motion can generate entanglement between artificial atoms and protect them from spontaneous emission. We consider a pair of superconducting qubits coupled to a resonator mode, where the modulation of the coupling strength can mimic the harmonic motion of the qubits at relativistic speeds, generating acceleration radiation. We find the optimal feasible conditions for generating a stationary entangled state between the qubits when they are initially prepared in their ground state. Furthermore, we analyze the effects of motion on the probability of spontaneous emission in the standard scenarios of single-atom and two-atom superradiance, where one or two excitations are initially present. Finally, we show that relativistic motion induces sub-radiance and can generate a Zeno-like effect, preserving the excitations from radiative decay.


title: Optical properties of superconducting EuFe2(As1-xPx)2

abstract: We present a broadband optical-conductivity study of superconducting single-crystalline EuFe2(As1-xPx)2 with three different substitutional levels. We analyze the normal-state electrodynamics by decomposing the conductivity spectra using a Drude-Lorentz model with two Drude terms representing two groups of carriers with different scattering rates. The analysis reveals that the scattering rate of at least one of the Drude components develops linearly with temperature for each doping level. This points towards strong electron-electron correlations and a non-Fermi-liquid behavior in the P-substituted superconducting Eu-122 pnictides. We also detect a transfer of the spectral weight from mid-infrared to higher frequencies and assign it to the Hund's-rule coupling between itinerant and localized carriers. The conductivity spectra below the superconducting transition show no sharp features to be associated with the dirty-limit superconducting BCS gaps. We interpret these results in terms of clean-limit superconductivity in EuFe2(As1-xPx)2. The best parametrization fit can be achieved using a two-gap model. We find that the larger gap at the hole pockets of the Fermi surface is likely to be isotropic, while the smaller gap at the electron pockets is anisotropic or even nodal.


title: A constant-time algorithm for middle levels Gray codes

abstract: For any integer $n\geq 1$ a middle levels Gray code is a cyclic listing of all $n$-element and $(n+1)$-element subsets of $\{1,2,\ldots,2n+1\}$ such that any two consecutive subsets differ in adding or removing a single element. The question whether such a Gray code exists for any $n\geq 1$ has been the subject of intensive research during the last 30 years, and has been answered affirmatively only recently [T. M\"utze. Proof of the middle levels conjecture. Proc. London Math. Soc., 112(4):677--713, 2016]. In a follow-up paper [T. M\"utze and J. Nummenpalo. An efficient algorithm for computing a middle levels Gray code. To appear in ACM Transactions on Algorithms, 2018] this existence proof was turned into an algorithm that computes each new set in the Gray code in time $\mathcal{O}(n)$ on average. In this work we present an algorithm for computing a middle levels Gray code in optimal time and space: each new set is generated in time $\mathcal{O}(1)$ on average, and the required space is $\mathcal{O}(n)$.


title: On pair correlation for generic diagonal forms

abstract: We establish new pair correlation results for certain generic homogenous diagonal forms evaluated on the integers. Methods are analytic leading to explicit quantitative statements.


title: Taming the higher power corrections in semileptonic B decays

abstract: We study the effect of dimension 7 and 8 operators on inclusive semileptonic B decays and the extraction of |Vcb|. Using moments of semileptonic B decay spectra and information based on the Lowest-Lying State saturation Approximation (LLSA) we perform a global fit of the nonperturbative parameters of the heavy quark expansion including for the first time the O(1/mb^{4,5}) contributions. Higher power corrections appear to have a very small effect on the extraction of |Vcb|, independently of the weight we attribute to the LLSA.


title: Scheme variations of the QCD coupling and hadronic $\tau$ decays

abstract: The Quantum Chromodynamics (QCD) coupling, $\alpha_s$, is not a physical observable of the theory since it depends on conventions related to the renormalization procedure. We introduce a definition of the QCD coupling, denoted by $\hat\alpha_s$, whose running is explicitly renormalization scheme invariant. The scheme dependence of the new coupling $\hat\alpha_s$ is parameterized by a single parameter $C$, related to transformations of the QCD scale $\Lambda$. It is demonstrated that appropriate choices of $C$ can lead to substantial improvements in the perturbative prediction of physical observables. As phenomenological applications, we study $e^+e^-$ scattering and decays of the $\tau$ lepton into hadrons, both being governed by the QCD Adler function.


title: Vortex reconnection in the three dimensional Navier-Stokes equations

abstract: We prove that the vortex structures of solutions to the 3D Navier-Stokes equations can change their topology without any loss of regularity. More precisely, we construct smooth high-frequency solutions to the Navier-Stokes equations where vortex lines and vortex tubes of arbitrarily complicated topologies are created and destroyed in arbitrarily small times. This instance of vortex reconnection is structurally stable and in perfect agreement with the existing computer simulations and experiments. We also provide a (non-structurally stable) scenario where the destruction of vortex structures is instantaneous.


title: Understanding Innovation to Drive Sustainable Development

abstract: Innovation is among the key factors driving a country's economic and social growth. But what are the factors that make a country innovative? How do they differ across different parts of the world and different stages of development? In this work done in collaboration with the World Economic Forum (WEF), we analyze the scores obtained through executive opinion surveys that constitute the WEF's Global Competitiveness Index in conjunction with other country-level metrics and indicators to identify actionable levers of innovation. The findings can help country leaders and organizations shape the policies to drive developmental activities and increase the capacity of innovation.


title: Resolving the Schwarzschild singularity in both classic and quantum gravity

abstract: The Schwarzschild singularity's resolution has key values in cracking the key mysteries related with black holes, the origin of their horizon entropy and the information missing puzzle involved in their evaporations. We provide in this work the general dynamic inner metric of collapsing stars with horizons and with non-trivial radial mass distributions. We find that static central singularities are not the final state of the system. Instead, the final state of the system is a periodically zero-cross breathing ball. Through 3+1 decomposed general relativity and its quantum formulation, we establish a functional Schr\"odinger equation controlling the micro-state of this breathing ball and show that, the system configuration with all the matter concentrating on the central point is not the unique eigen-energy-density solution. Using a Bohr-Sommerfield like "orbital" quantisation assumption, we show that for each black hole of horizon radius $r_h$, there are about $e^{r_h^2/\ell^2_\mathrm{pl}}$ allowable eigen-energy-density profile. This naturally leads to physic interpretations for the micro-origin of horizon entropy, as well as solutions to the information missing puzzle involved in Hawking radiations.


title: On the prediction loss of the lasso in the partially labeled setting

abstract: In this paper we revisit the risk bounds of the lasso estimator in the context of transductive and semi-supervised learning. In other terms, the setting under consideration is that of regression with random design under partial labeling. The main goal is to obtain user-friendly bounds on the off-sample prediction risk. To this end, the simple setting of bounded response variable and bounded (high-dimensional) covariates is considered. We propose some new adaptations of the lasso to these settings and establish oracle inequalities both in expectation and in deviation. These results provide non-asymptotic upper bounds on the risk that highlight the interplay between the bias due to the mis-specification of the linear model, the bias due to the approximate sparsity and the variance. They also demonstrate that the presence of a large number of unlabeled features may have significant positive impact in the situations where the restricted eigenvalue of the design matrix vanishes or is very small.


title: Semi-classical resonances associated with a periodic orbit of hyperbolic type

abstract: We consider in this Note resonances for a $h$-Pseudo-Differential Operator $H(x,hD_x;h)$ on $L^2(M)$ induced by a periodic orbit of hyperbolic type, as arises for Schr\"odinger operator with AC Stark effect when $M={\bf R}^n$, or the geodesic flow on an axially symmetric manifold $M$, extending Poincar\'e example of Lagrangian systems with 2 degrees of freedom. We generalize the framework of [G\'eSj], in the sense that we allow for hyperbolic and elliptic eigenvalues of Poincar\'e map, and look for so-called semi-excited resonances with imaginary part of magnitude $-h\log h$, or $h^s$, with $0<s<1$.


title: Chirality modifies the interaction between knots

abstract: In this study we consider an idealization of a typical optical tweezers experiment involving a semiflexible double-knotted polymer, with steric hindrance and persistence length matching those of dsDNA in high salt concentration, under strong stretching. Using exhaustive Molecular Dynamics simulations we show that not only does a double-knotted dsDNA filament under tension possess a free energy minimum when the two knots are intertwined, but also that the depth of this minimum depends on the relative chirality of the two knots. We rationalize this dependence of the effective interaction on the chirality in terms of a competition between chain entropy and bending energy.


title: The Panchromatic Hubble Andromeda Treasury XV. The BEAST: Bayesian Extinction and Stellar Tool

abstract: We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys.


title: Asymptotically Optimal Approximation Algorithms for Coflow Scheduling

abstract: Many modern datacenter applications involve large-scale computations composed of multiple data flows that need to be completed over a shared set of distributed resources. Such a computation completes when all of its flows complete. A useful abstraction for modeling such scenarios is a {\em coflow}, which is a collection of flows (e.g., tasks, packets, data transmissions) that all share the same performance goal. In this paper, we present the first approximation algorithms for scheduling coflows over general network topologies with the objective of minimizing total weighted completion time. We consider two different models for coflows based on the nature of individual flows: circuits, and packets. We design constant-factor polynomial-time approximation algorithms for scheduling packet-based coflows with or without given flow paths, and circuit-based coflows with given flow paths. Furthermore, we give an $O(\log n/\log \log n)$-approximation polynomial time algorithm for scheduling circuit-based coflows where flow paths are not given (here $n$ is the number of network edges). We obtain our results by developing a general framework for coflow schedules, based on interval-indexed linear programs, which may extend to other coflow models and objective functions and may also yield improved approximation bounds for specific network scenarios. We also present an experimental evaluation of our approach for circuit-based coflows that show a performance improvement of at least 22% on average over competing heuristics.


title: Geometric approach to entanglement quantification with polynomial measures

abstract: We show that the quantification of entanglement of any rank-2 state with any polynomial entanglement measure can be recast as a geometric problem on the corresponding Bloch sphere. This approach provides novel insight into the properties of entanglement and allows us to relate different polynomial measures to each other, simplifying their quantification. In particular, unveiling and exploiting the geometric structure of the concurrence for two qubits, we show that the convex roof of any polynomial measure of entanglement can be quantified exactly for all rank-2 states of an arbitrary number of qubits which have only one or two unentangled states in their range. We give explicit examples by quantifying the three-tangle exactly for several representative classes of three-qubit states. We further show how our methods can be used to obtain analytical results for entanglement of more complex states if one can exploit symmetries in their geometric representation.


title: An high definition view of the COSMOS Wall at z~0.73

abstract: We present a study of a large filamentary structure at z~0.73 in the field of the COSMOS survey, the so-called COSMOS Wall. This structure encompasses a comprehensive range of environments from a dense cluster and a number of galaxy groups to filaments, less dense regions, and adjacent voids. It thus provides a valuable laboratory for the accurate mapping of environmental effects on galaxy evolution at a look-back time of ~6.5 Gyr, when the Universe was roughly half its present age. We performed deep spectroscopic observations with VIMOS at VLT of a K-band selected sample of galaxies in this complex structure, building a sample of galaxies complete in galaxy stellar mass down to a lower limit of log(M/M_sun)~ 9.8, which is significantly deeper than previously available data. Thanks to its location within the COSMOS survey, each galaxy benefits from a wealth of ancillary information. In this paper we detail the survey strategy and weighting scheme adopted to account for the biases introduced by the photometric pre-selection of our targets. We present our galaxy stellar mass and rest-frame magnitudes estimates together with a group catalog obtained with our new data and their member galaxies color/mass distribution. Owing to to our new sample we can perform a detailed, high definition mapping of the complex COSMOS Wall structure. The sharp environmental information, coupled with high quality spectroscopic information and rich ancillary data available in the COSMOS field, enables a detailed study of galaxy properties as a function of local environment in a redshift slice where environmental effects are important, and in a stellar mass range where mass and environment driven effects are both at work.


title: Fundamental Implications of Intergalactic Magnetic Field Observations

abstract: Helical intergalactic magnetic fields at the $\sim 10^{-14}~{\rm G}$ level on $\sim 10~{\rm Mpc}$ length scales are indicated by current gamma ray observations. The existence of magnetic fields in cosmic voids and their non-trivial helicity suggest that they must have originated in the early universe and thus have implications for the fundamental interactions. I combine present knowledge of the observational constraints and the dynamics of cosmological magnetic fields to derive characteristics that would need to be explained by the magnetic field generation mechanism. The importance of CP violation and a possible crucial role for chiral effects in the early universe are pointed out.


title: Extended homologous series of Sn-O layered systems: a first-principles study

abstract: Apart from the most studied tin-oxide compounds, SnO and SnO2, intermediate states have been claimed to exist for more than a hundred years. In addition to the known homologous series (Seko et al., Phys. Rev. Lett. 100, 045702 (2008)), we here predict the existence of several new compounds with an O concentration between 50 % (SnO) and 67 % (SnO2). All these intermediate compounds are constructed from removing one or more (101) oxygen layers of SnO2. Since the van der Waals (vdW) interaction is known to be important for the Sn-Sn interlayer distances, we use a vdW-corrected functional, and compare these results with results obtained with PBE and hybrid functionals. We present the electronic properties of the intermediate structures and we observe a decrease of the band gap when (i) the O concentration increases and (ii) more SnO-like units are present for a given concentration. The contribution of the different atoms to the valence and conduction band is also investigated.


title: Analysis of uncertainties in alpha-particle optical-potential assessment below the Coulomb barrier

abstract: Background: Recent high-precision measurements of alpha-induced reaction data below the Coulomb barrier have pointed out questions of the alpha-particle optical-model potential (OMP) which are yet open within various mass ranges. Purpose: The applicability of a previous optical potential and eventual uncertainties and/or systematic errors of the OMP assessment at low energies can be further considered on this basis. Method: Nuclear model parameters based on the analysis of recent independent data, particularly gamma-ray strength functions, have been involved within statistical model calculation of the (alpha,x) reaction cross sections. Results: The above-mentioned potential provides a consistent description of the recent alpha-induced reaction data with no empirical rescaling factors of the and/or nucleon widths. Conclusions: A suitable assessment of alpha-particle optical potential below the Coulomb barrier should involve the statistical-model parameters beyond this potential on the basis of a former analysis of independent data.


title: Role of vertex corrections in the matrix formulation of the random phase approximation for the multiorbital Hubbard model

abstract: In the framework of a multiorbital Hubbard model description of superconductivity, a matrix formulation of the superconducting pairing interaction that has been widely used is designed to treat spin, charge and orbital fluctuations within a random phase approximation (RPA). In terms of Feynman diagrams, this takes into account particle-hole ladder and bubble contributions as expected. It turns out, however, that this matrix formulation also generates additional terms which have the diagrammatic structure of vertex corrections. Here we examine these terms and discuss the relationship between the matrix-RPA superconducting pairing interaction and the Feynman diagrams that it sums.


title: Essential self-adjointness of powers of first-order differential operators on non-compact manifolds with low-regularity metrics

abstract: We consider first-order differential operators with locally bounded measurable coefficients on vector bundles with measurable coefficient metrics. Under a mild set of assumptions, we demonstrate the equivalence between the essential self-adjointness of such operators to a negligible boundary property. When the operator possesses higher regularity coefficients, we show that higher powers are essentially self-adjoint if and only if this condition is satisfied. In the case that the low-regularity Riemannian metric induces a complete length space, we demonstrate essential self-adjointness of the operator and its higher powers up to the regularity of its coefficients. We also present applications to Dirac operators on Dirac bundles when the metric is non-smooth.


title: Revisiting on-shell renormalization conditions in theories with flavour mixing

abstract: In this review, we present a derivation of the on-shell renormalization conditions for scalar and fermionic fields in theories with and without parity conservation. We also discuss the specifics of Majorana fermions. Our approach only assumes a canonical form for the renormalized propagators and exploits the fact that the inverse propagators are non-singular in $\varepsilon = p^2 - m_n^2$, where $p$ is the external four-momentum and $m_n$ is a pole mass. In this way, we obtain full agreement with commonly used on-shell conditions. We also discuss how they are implemented in renormalization.


title: A Novel Quasi-One-Dimensional Topological Insulator in Bismuth Iodide $\beta$-Bi$_4$I$_4$

abstract: Recent progress in the field of topological states of matter(1,2) has largely been initiated by the discovery of bismuth and antimony chalcogenide bulk topological insulators (TIs)(3-6), followed by closely related ternary compounds(7-16) and predictions of several weak TIs(17-19). However, both the conceptual richness of Z$_2$ classification of TIs as well as their structural and compositional diversity are far from being fully exploited. Here, a new Z$_2$ topological insulator is theoretically predicted and experimentally confirmed in the $\beta$-phase of quasi-one-dimensional bismuth iodide Bi$_4$I$_4$. The electronic structure of $\beta$-Bi$_4$I$_4$, characterized by Z$_2$ invariants (1;110), is in proximity of both the weak TI phase (0;001) and the trivial insulator phase (0;000). Our angle-resolved photoemission spectroscopy measurements on the (001) surface reveal a highly anisotropic band-crossing feature located at the point of the surface Brillouin zone and showing no dispersion with the photon energy, thus being fully consistent with the theoretical prediction.


title: Single quark entropy and the Polyakov loop

abstract: We study Quantum Chromodynamics (QCD) with 2+1 flavors with almost physical quark masses using the highly improved staggered quark action (HISQ). We calculate the Polyakov loop in a wide temperature range, obtain the free energy and the entropy of a single static quark and discuss the QCD crossover region in detail. We show that the entropy has a peak close to the chiral crossover and consider the consequences for the deconfinement aspects of the crossover phenomena. We study the renormalized Polyakov loop susceptibilities and place them into the context of the crossover. We also obtain a quantitative result for the onset of weak coupling behavior at high temperatures.


title: Gauge-invariant description of Higgs phenomenon and quark confinement

abstract: We propose a novel description for the Higgs mechanism by which a gauge boson acquires the mass. We do not assume spontaneous breakdown of gauge symmetry signaled by a non-vanishing vacuum expectation value of the scalar field. In fact, we give a manifestly gauge-invariant description of the Higgs mechanism in the operator level, which does not rely on spontaneous symmetry breaking. This enables us to discuss the confinement-Higgs complementarity from a new perspective. The "Abelian" dominance in quark confinement of the Yang-Mills theory is understood as a consequence of the gauge-invariant Higgs phenomenon for the relevant Yang-Mills-Higgs model.


title: Dirac fermions at high-index surfaces of bismuth chalcogenide topological insulator nanostructures

abstract: Binary bismuth chalcogenides Bi$_2$Se$_3$, Bi$_2$Te$_3$, and related materials are currently being extensively investigated as the reference topological insulators (TIs) due to their simple surface-state band dispersion (single Dirac cone) and relatively large bulk band gaps. Nanostructures of TIs are of particular interest as an increased surface-to-volume ratio enhances the contribution of surfaces states, meaning they are promising candidates for potential device applications. So far, the vast majority of research efforts have focused on the low-energy (0001) surfaces, which correspond to natural cleavage planes in these layered materials. However, the surfaces of low-dimensional nanostructures (nanoplatelets, nanowires, nanoribbons) inevitably involve higher-index facets. We perform a systematic ab initio investigation of the surfaces of bismuth chalcogenide TI nanostructures characterized by different crystallographic orientations, atomic structures and stoichiometric compositions. We find several stable terminations of high-index surfaces, which can be realized at different values of the chemical potential of one of constituent elements. For the uniquely defined stoichiometric termination, the topological Dirac fermion states are shown to be strongly anisotropic with a clear dependence of Fermi velocities and spin polarization on the surface orientation. Self-doping effects and the presence of topologically trivial mid-gap states are found to characterize the non-stoichiometric surfaces. The results of our study pave the way towards experimental control of topologically protected surface states in bismuth chalcogenide nanostructures.


title: Gender Representation on Journal Editorial Boards in the Mathematical Sciences

abstract: We study gender representation on the editorial boards of 435 journals in the mathematical sciences. Women are known to comprise approximately 15% of tenure-stream faculty positions in doctoral-granting mathematical sciences departments in the United States. Compared to this pool, the likely source of journal editorships, we find that 8.9% of the 13067 editorships in our study are held by women. We describe group variations within the editorships by identifying specific journals, subfields, publishers, and countries that significantly exceed or fall short of this average. To enable our study, we develop a semi-automated method for inferring gender that has an estimated accuracy of 97.5%. Our findings provide the first measure of gender distribution on editorial boards in the mathematical sciences, offer insights that suggest future studies in the mathematical sciences, and introduce new methods that enable large-scale studies of gender distribution in other fields.


title: Polymetric Rhythmic Feel for a Cognitive Drum Computer

abstract: This paper addresses a question about music cognition: how do we derive polymetric structures. A preference rule system is presented which is implemented into a drum computer. The preference rule system allows inferring local polymetric structures, like two-over-three and three-over-two. By analyzing the micro-timing of West African percussion music a timing pattern consisting of six pulses was discovered. It integrates binary and ternary rhythmic feels. The presented drum computer integrates the discovered superimposed polymetric swing (timing and velocity) appropriate to the rhythmic sequence the user inputs. For binary sequences, the amount of binary swing is increased and for ternary sequences, the ternary swing is increased.


title: Evidence for collectivity in pp collisions at the LHC

abstract: Measurements of two- and multi-particle angular correlations in pp collisions at sqrt(s) = 5, 7, and 13 TeV are presented as a function of charged-particle multiplicity. The data, corresponding to integrated luminosities of 1.0 inverse picobarn (5 TeV), 6.2 inverse picobarns (7 TeV), and 0.7 inverse picobarns (13 TeV), were collected using the CMS detector at the LHC. The second-order (v[2]) and third-order (v[3]) azimuthal anisotropy harmonics of unidentified charged particles, as well as v[2] of K0 short and Lambda/anti-Lambda particles, are extracted from long-range two-particle correlations as functions of particle multiplicity and transverse momentum. For high-multiplicity pp events, a mass ordering is observed for the v[2] values of charged hadrons (mostly pions), K0 short, and Lambda/anti-Lambda, with lighter particle species exhibiting a stronger azimuthal anisotropy signal below pt of about 2 GeV/c. For 13 TeV data, the v[2] signals are also extracted from four- and six-particle correlations for the first time in pp collisions, with comparable magnitude to those from two-particle correlations. These observations are similar to those seen in pPb and PbPb collisions, and support the interpretation of a collective origin for the observed long-range correlations in high-multiplicity pp collisions.


title: A variational H(div) finite element discretisation approach for perfect incompressible fluids

abstract: We propose a finite element discretisation approach for the incompressible Euler equations which mimics their geometric structure and their variational derivation. In particular, we derive a finite element method that arises from a nonholonomic variational principle and an appropriately defined Lagrangian, where finite element H(div) vector fields are identified with advection operators; this is the first successful extension of the structure-preserving discretisation of Pavlov et al. (2009) to the finite element setting. The resulting algorithm coincides with the energy-conserving scheme presented in Guzm\'an et al. (2016). Through the variational derivation, we discover that it also satisfies a discrete analogous of Kelvin's circulation theorem. Further, we propose an upwind-stabilised version of the scheme which dissipates enstrophy whilst preserving energy conservation and the discrete Kelvin's theorem. We prove error estimates for this version of the scheme, and we study its behaviour through numerical tests.


title: Electron Two-stream Instability and Its Application in Solar and Heliophysics

abstract: It is well known that electron beams accelerated in solar flares can drive two-stream instability and produce radio bursts in the solar corona as well as in the interplanetary medium. Recent observations show that the solar wind likely originates from nanoflare-like events near the surface of the Sun where locally heated plasma escapes along open field lines into space. Recent numerical simulations and theoretical studies show that electron two-stream instability (ETSI) driven by nanoflare-accelerated electron beams can produce the observed nanoflare-type radio bursts, the non-Maxwellian electron velocity distribution function of the solar wind, and the kinetic scale turbulence in solar wind. This brief review focus on the basic theoretical framework and recent progress in the nonlinear evolution of ETSI, including the formation of electron holes, Langmuir wave generation in warm plasma, and the nonlinear modulation instability and Langmuir collapse. Potential applications in heliophysics and astrophysics are discussed.


title: Primal-dual interior-point multigrid method for topology optimization

abstract: An interior point method for the structural topology optimization is proposed. The linear systems arising in the method are solved by the conjugate gradient method preconditioned by geometric multigrid. The resulting method is then compared with the so-called optimality condition method, an established technique in topology optimization. This method is also equipped with the multigrid preconditioned conjugate gradient algorithm. We conclude that, for large scale problems, the interior point method with an inexact iterative linear solver is superior to any other variant studied in the paper.


title: Regular global attractors for wave equations with degenerate memory

abstract: We consider the wave equation with degenerate viscoelastic dissipation recently examined in Cavalcanti, Fatori, and Ma, Attractors for wave equations with degenerate memory, J. Differential Equations (2016). Under some additional assumptions, we show the existence of a compact absorbing set. This result provides further regularity for the global attractor and shows that it consists of regular solutions.


title: Application of the generalized Kirchhoff's law to calculation of photoluminescence spectra of one-dimensional photonic crystals

abstract: The approach based on the generalized Kirchhoff's law for calculating photoluminescence (PL) spectra of one-dimensional (1D) multi-layered structures, in particular, 1D photonic crystals has been developed. It is valid in the local thermodynamic equilibrium approximation and leads to simple and explicit expression for the photoluminescence intensity. In the framework of the present theory the analytical expression for the spontaneous emission intensity enhancement factor (IEF) for a 1D photonic crystal has been derived. It takes a particularly simple form in the case of a sufficiently large number of the layers and is well suitable for analysis; in particular, it explains the difference in the emission intensity at frequencies near different edges of photonic band-gaps (PBGs), where the intensity is relatively high, and specificity of suppression of the emission in a given frequency range. Also, the developed approach is discussed in connection with the standard method using the Fermi's golden rule and the concept of the local density of (optical) states (LDOS).


title: Parallel Priority-Flood Depression Filling For Trillion Cell Digital Elevation Models On Desktops Or Clusters

abstract: Algorithms for extracting hydrologic features and properties from digital elevation models (DEMs) are challenged by large datasets, which often cannot fit within a computer's RAM. Depression filling is an important preconditioning step to many of these algorithms. Here, I present a new, linearly-scaling algorithm which parallelizes the Priority-Flood depression-filling algorithm by subdividing a DEM into tiles. Using a single-producer, multi-consumer design, the new algorithm works equally well on one core, multiple cores, or multiple machines and can take advantage of large memories or cope with small ones. Unlike previous algorithms, the new algorithm guarantees a fixed number of memory access and communication events per subdivision of the DEM. In comparison testing, this results in the new algorithm running generally faster while using fewer resources than previous algorithms. For moderately sized tiles, the algorithm exhibits ~60% strong and weak scaling efficiencies up to 48 cores, and linear time scaling across datasets ranging over three orders of magnitude. The largest dataset on which I run the algorithm has 2 trillion (2*10^12) cells. With 48 cores, processing required 4.8 hours wall-time (9.3 compute-days). This test is three orders of magnitude larger than any previously performed in the literature. Complete, well-commented source code and correctness tests are available for download from a repository.


title: Exact special twist method for quantum Monte Carlo simulations

abstract: We present a systematic investigation of the special twist method introduced by Rajagopal $\textit{et al.}$ [ Phys. Rev. B 51, 10591 (1995) ] for reducing finite-size effects in correlated calculations of periodic extended systems with Coulomb interactions and Fermi statistics. We propose a procedure for finding special twist values which, at variance with previous applications of this method, reproduce the energy of the mean-field infinite-size limit solution within an adjustable (arbitrarily small) numerical error. This choice of the special twist is shown to be the most accurate single-twist solution for curing one-body finite-size effects in correlated calculations. For these reasons we dubbed our procedure "exact special twist" (EST). EST only needs a fully converged independent-particles or mean-field calculation within the primitive cell and a simple fit to find the special twist along a specific direction in the Brillouin zone. We first assess the performances of EST in a simple correlated model such as the 3D electron gas. Afterwards, we test its efficiency within $\textit{ab initio}$ quantum Monte Carlo simulations of metallic elements of increasing complexity. We show that EST displays an overall good performance in reducing finite-size errors comparable to the widely used twist average technique but at a much lower computational cost, since it involves the evaluation of just one wavefunction. We also demonstrate that the EST method shows similar performances in the calculation of correlation functions, such as the ionic forces for structural relaxation and the pair radial distribution function in liquid hydrogen. Our conclusions point to the usefulness of EST for correlated supercell calculations, our method will be particularly relevant when the physical problem under consideration requires large periodic cells.


title: ArbiLoMod: Local Solution Spaces by Random Training in Electrodynamics

abstract: The simulation method ArbiLoMod has the goal to provide users of Finite Element based simulation software with quick re-simulation after localized changes to the model under consideration. It generates a Reduced Order Model (ROM) for the full model without ever solving the full model. To this end, a localized variant of the Reduced Basis method is employed, solving only small localized problems in the generation of the reduced basis. The key to quick re-simulation lies in recycling most of the localized basis vectors after a localized model change. In this publication, ArbiLoMod's local training algorithm is analyzed numerically for the non-coercive problem of time harmonic Maxwell's equations in 2D, formulated in H(curl).