title: Green function for linearized Navier-Stokes around a boundary layer profile: near critical layers

abstract: This is a continuation and completion of the program (initiated in \cite{GrN1,GrN2}) to derive pointwise estimates on the Green function and sharp bounds on the semigroup of linearized Navier-Stokes around a generic stationary boundary layer profile. This is done via a spectral analysis approach and a careful study of the Orr-Sommerfeld equations, or equivalently the Navier-Stokes resolvent operator $(\lambda - L)^{-1}$. The earlier work (\cite{GrN1,GrN2}) treats the Orr-Sommerfeld equations away from critical layers: this is the case when the phase velocity is away from the range of the background profile or when $\lambda$ is away from the Euler continuous spectrum. In this paper, we study the critical case: the Orr-Sommerfeld equations near critical layers, providing pointwise estimates on the Green function as well as carefully studying the Dunford's contour integral near the critical layers. As an application, we obtain pointwise estimates on the Green function and sharp bounds on the semigroup of the linearized Navier-Stokes problem near monotonic boundary layers that are spectrally stable to the Euler equations, complementing \cite{GrN1,GrN2} where unstable profiles are considered.


title: Parameterised Post-Newtonian Expansion in Screened Regions

abstract: The parameterised post-Newtonian (PPN) formalism has enabled stringent tests of static weak-field gravity in a theory-independent manner. Here we incorporate screening mechanisms of modified gravity theories into the framework by introducing an effective gravitational coupling and defining the PPN parameters as functions of position. To determine these functions we develop a general method for efficiently performing the post-Newtonian expansion in screened regimes. For illustration, we derive all the PPN functions for a cubic galileon and a chameleon model. We also analyse the Shapiro time delay effect for these two models and find no deviations from General Relativity insofar as the signal path and the perturbing mass reside in a screened region of space.


title: Spectral Scaling Technique to Determine Extreme Carrington-level Geomagnetically Induced Currents Effects

abstract: Space weather events produce variations in the electric current in the Earth's magnetosphere and ionosphere. From these high altitude atmospheric regions, resulting geomagnetically induced currents (GICs) can lead to fluctuations in ground currents that affect the electric power grid and potentially overload transformers during extreme storms. The most extreme geomagnetic storm on record, known as the 1859 Carrington event, was so intense that ground-based magnetometers were saturated at high magnetic latitudes. The most reliable, un-saturated observation is the hour-resolution data from the Colaba Magnetic Observatory in India. However, higher frequency components -- fluctuations at second through minute time cadence -- to the magnetic field can play a significant role in GIC-related effects. We present a new method for scaling higher frequency observations to create a realistic Carrington-like event magnetic field model, using modern magnetometer observations. Using the magnetic field model and ground conductivity models, we produce an electric field model. This method can be applied to create similar magnetic and electric field models for studies of GIC effects on power-grids.


title: Constrained Bayesian Networks: Theory, Optimization, and Applications

abstract: We develop the theory and practice of an approach to modelling and probabilistic inference in causal networks that is suitable when application-specific or analysis-specific constraints should inform such inference or when little or no data for the learning of causal network structure or probability values at nodes are available. Constrained Bayesian Networks generalize a Bayesian Network such that probabilities can be symbolic, arithmetic expressions and where the meaning of the network is constrained by finitely many formulas from the theory of the reals. A formal semantics for constrained Bayesian Networks over first-order logic of the reals is given, which enables non-linear and non-convex optimisation algorithms that rely on decision procedures for this logic, and supports the composition of several constrained Bayesian Networks. A non-trivial case study in arms control, where few or no data are available to assess the effectiveness of an arms inspection process, evaluates our approach. An open-access prototype implementation of these foundations and their algorithms uses the SMT solver Z3 as decision procedure, leverages an open-source package for Bayesian inference to symbolic computation, and is evaluated experimentally.


title: Plankton: the paradox and the power law

abstract: Two basic features of assemblages of unicellular plankton: (1) their high biodiversity and (2) the power-law structure of their abundance, can be explained by an allometric scaling of cell growth and mortality with respect to cell size. To show this, we describe a numerical study of a size-structured, multispecies, population-dynamic model; the model has a single resource, supporting an arbitrary number of phytoplankton and zooplankton species. If the number of plankton species is large enough, the death rate of prey and cell growth rate of predators have approximate allometric scalings with cell size. Together, these scalings give rise to an equilibrium distribution of abundance near the power law, on which many species can coexist. Scalings of this kind cannot be achieved if the number of species is small. This suggests that the conjunction of species-richness and power-law structures in plankton communities is more than a coincidence. Although the exact allometric scalings used here should not be expected in practice, exclusion of species should be relatively slow if they lie close to the power law. Thus the forces needed to achieve coexistence could be effective, even if they are relatively weak.


title: Multipole analysis of redshift-space distortions around cosmic voids

abstract: We perform a comprehensive redshift-space distortion analysis based on cosmic voids in the large-scale distribution of galaxies observed with the Sloan Digital Sky Survey. To this end, we measure multipoles of the void-galaxy cross-correlation function and compare them with standard model predictions in cosmology. Merely considering linear-order theory allows us to accurately describe the data on the entire available range of scales and to probe void-centric distances down to about $2h^{-1}{\rm Mpc}$. Common systematics, such as the Fingers-of-God effect, scale-dependent galaxy bias, and nonlinear clustering do not seem to play a significant role in our analysis. We constrain the growth rate of structure via the redshift-space distortion parameter $\beta$ at two median redshifts, $\beta(\bar{z}=0.32)=0.599^{+0.134}_{-0.124}$ and $\beta(\bar{z}=0.54)=0.457^{+0.056}_{-0.054}$, with a precision that is competitive with state-of-the-art galaxy-clustering results. While the high-redshift constraint perfectly agrees with model expectations, we observe a mild $2\sigma$ deviation at $\bar{z}=0.32$, which increases to $3\sigma$ when the data is restricted to the lowest available redshift range of $0.15<z<0.33$.


title: Galactic Orbits of Globular Clusters in the Region of the Galactic Bulge

abstract: Galactic orbits have been constructed over long time intervals for ten globular clusters located near the Galactic center. A model with an axially symmetric gravitational potential for the Galaxy was initially applied, after which a non-axially symmetric potential corresponding to the central bar was added. Variations in the trajectories of all these globular clusters in the XY plane due to the influence of the bar were detected. These were greatest for the cluster Terzan 4 in the meridional (RZ) plane. The globular clusters Terzan 1, Terzan 2, Terzan 4, Terzan 9, NGC 6522, and NGC 6558 always remained within the Galactic bulge, no farther than 4 kpc from the Galactic center.


title: A new length scale, and modified Einstein-Cartan-Dirac equations for a point mass

abstract: We have recently proposed a new action principle for combining Einstein equations and the Dirac equation for a point mass. We used a length scale $L_{CS}$, dubbed the Compton-Schwarzschild length, to which the Compton wavelength and Schwarzschild radius are small mass and large mass approximations, respectively. Here we write down the field equations which follow from this action. We argue that the large mass limit yields Einstein equations, provided we assume wave function collapse and localisation for large masses. The small mass limit yields the Dirac equation. We explain why the Kerr-Newman black hole has the same gyromagnetic ratio as the Dirac electron, both being twice the classical value. The small mass limit also provides compelling reasons for introducing torsion, which is sourced by the spin density of the Dirac field. There is thus a symmetry between torsion and gravity: torsion couples to quantum objects through Planck's constant $\hbar$ (but not $G$) and is important in the microscopic limit. Whereas gravity couples to classical matter, as usual, through Newton's gravitational constant $G$ (but not $\hbar$), and is important in the macroscopic limit. We construct the Einstein-Cartan-Dirac equations which include the length $L_{CS}$. We find a potentially significant change in the coupling constant of the torsion driven cubic non-linear self-interaction term in the Dirac-Hehl-Datta equation. We speculate on the possibility that gravity is not a fundamental interaction, but emerges as a consequence of wave function collapse, and that the gravitational constant maybe expressible in terms of Planck's constant and the parameters of dynamical collapse models.


title: The denominators of power sums of arithmetic progressions

abstract: In a recent paper the authors studied the denominators of polynomials that represent power sums by Bernoulli's formula. Here we extend our results to power sums of arithmetic progressions. In particular, we obtain a simple explicit criterion for integrality of the coefficients of these polynomials. As applications, we obtain new results on the sequence of denominators of the Bernoulli polynomials. A consequence is that certain quotients of successive denominators are infinitely often integers, which we characterize.


title: The Milky Way's Central Molecular Zone

abstract: This review compiles the results of recent studies of molecular gas conditions in the central six hundred parsecs of our Galaxy. The review begins by placing our Galactic center into context with the rest of our galaxy. It next discusses the wealth of previous research on the Galactic center, before focusing on what is known about the molecular interstellar medium in this region. It focuses especially on a surge in interest in this region and new studies conducted in the last five years. It concludes by highlighting open questions that remain, and the potential for new facilities such as ALMA to make progress in resolving these uncertainties.


title: Nanoplatelets as material system between strong confinement and weak confinement

abstract: Recently, the fabrication of CdSe nanoplatelets became an important research topic. Nanoplatelets are often described as having a similar electronic structure as 2D dimensional quantum wells and are promoted as colloidal quantum wells with monolayer precision width. In this paper, we show, that nanoplatelets are not ideal quantum wells, but cover depending on the size: the strong confinement regime, an intermediate regime and a Coulomb dominated regime. Thus, nanoplatelets are an ideal platform to study the physics in these regimes. Therefore, the exciton states of the nanoplatelets are numerically calculated by solving the full four dimensional Schr\"odinger equation. We compare the results with approximate solutions from semiconductor quantum well and quantum dot theory. The paper can also act as review of these concepts for the colloidal nanoparticle community.


title: Evolutionary dynamics of the cryptocurrency market

abstract: The cryptocurrency market surpassed the barrier of \$100 billion market capitalization in June 2017, after months of steady growth. Despite its increasing relevance in the financial world, however, a comprehensive analysis of the whole system is still lacking, as most studies have focused exclusively on the behaviour of one (Bitcoin) or few cryptocurrencies. Here, we consider the history of the entire market and analyse the behaviour of 1,469 cryptocurrencies introduced between April 2013 and June 2017. We reveal that, while new cryptocurrencies appear and disappear continuously and their market capitalization is increasing (super-)exponentially, several statistical properties of the market have been stable for years. These include the number of active cryptocurrencies, the market share distribution and the turnover of cryptocurrencies. Adopting an ecological perspective, we show that the so-called neutral model of evolution is able to reproduce a number of key empirical observations, despite its simplicity and the assumption of no selective advantage of one cryptocurrency over another. Our results shed light on the properties of the cryptocurrency market and establish a first formal link between ecological modelling and the study of this growing system. We anticipate they will spark further research in this direction.

categories: physics.soc-ph cs.SI nlin.AO q-bio.PE q-fin.GN


title: Horn Maximum Satisfiability: Reductions, Algorithms & Applications

abstract: Recent years have witness remarkable performance improvements in maximum satisfiability (MaxSAT) solvers. In practice, MaxSAT algorithms often target the most generic MaxSAT formulation, whereas dedicated solvers, which address specific subclasses of MaxSAT, have not been investigated. This paper shows that a wide range of optimization and decision problems are either naturally formulated as MaxSAT over Horn formulas, or permit simple encodings using Horn MaxSAT. Furthermore, the paper also shows how linear time decision procedures for Horn formulas can be used for developing novel algorithms for the Horn MaxSAT problem.


title: Spectral estimates for Schr\"odinger operators on periodic discrete graphs

abstract: We consider normalized Laplacians and their perturbations by periodic potentials (Schr\"odinger operators) on periodic discrete graphs. The spectrum of the operators consists of an absolutely continuous part (a union of a finite number of non-degenerate bands) and a finite number of flat bands, i.e., eigenvalues of infinite multiplicity. We obtain estimates of the Lebesgue measure of the spectrum in terms of geometric parameters of the graphs and show that they become identities for some class of graphs. We determine two-sided estimates on the length of the first spectral band and on the effective mass at the bottom of the spectrum of the Laplace and Schr\"odinger operators. In particular, these estimates yield that the first spectral band of Schr\"odinger operators is non-degenerate.


title: QCD Corrections to the Color-singlet $J/\psi$ Production in Deeply Inelastic Scattering at HERA

abstract: We present the first study of the QCD corrections to the color-singlet (CS) $J/\psi$ production in deeply inelastic $ep$ scattering at HERA. The $K$-factor ranges from $0.85$ to $2.38$ in the kinematic regions we study. In low transverse momentum regions, the $K$-factors is even smaller, and close to 1, which indicates good convergence of the perturbative expansion. With the QCD corrections, the CS cross section is still below the data. At least at QCD next-to-leading order, the color-octet mechanism is necessary to describe the data.


title: Complex Organic Molecules Towards Embedded Low-Mass Protostars

abstract: Complex organic molecules (COMs) have been observed towards several low-mass young stellar objects (LYSOs). Small and heterogeneous samples have so far precluded conclusions on typical COM abundances, as well as the origin(s) of abundance variations between sources. We present observations towards 16 deeply embedded (Class 0/I) low-mass protostars using the IRAM 30m telescope. We detect CH$_2$CO, CH$_3$CHO, CH$_3$OCH$_3$, CH$_3$OCHO, CH$_3$CN, HNCO, and HC$_3$N towards 67%, 37%, 13%, 13%, 44%, 81%, and 75% of sources respectively. Median column densities derived using survival analysis range between 6.0x10$^{10}$ cm$^{-2}$ (CH$_3$CN) and 2.4x10$^{12}$ cm$^{-2}$ (CH$_3$OCH$_3$) and median abundances range between 0.48% (CH$_3$CN) and 16% (HNCO) with respect to CH$_3$OH. Column densities for each molecule vary by about one order of magnitude across the sample. Abundances with respect to CH$_3$OH are more narrowly distributed, especially for oxygen-bearing species. We compare observed median abundances with a chemical model for low-mass protostars and find fair agreement, although some modeling work remains to bring abundances higher with respect to CH$_3$OH. Median abundances with respect to CH$_3$OH in LYSOs are also found to be generally comparable to observed abundances in hot cores, hot corinos, and massive young stellar objects. Compared with comets, our sample is comparable for all molecules except HC$_3$N and CH$_2$CO, which likely become depleted at later evolutionary stages.


title: A Modular Regularized Variational Multiscale Proper Orthogonal Decomposition for Incompressible Flows

abstract: In this paper, we propose, analyze and test a post-processing implementation of a projection-based variational multiscale (VMS) method with proper orthogonal decomposition (POD) for the incompressible Navier-Stokes equations. The projection-based VMS stabilization is added as a separate post-processing step to the standard POD approximation, and since the stabilization step is completely decoupled, the method can easily be incorporated into existing codes, and stabilization parameters can be tuned independent from the time evolution step. We present a theoretical analysis of the method, and give results for several numerical tests on benchmark problems which both illustrate the theory and show the proposed method's effectiveness.


title: The Real No-Boundary Wave Function in Lorentzian Quantum Cosmology

abstract: It is shown that the standard no-boundary wave function has a natural expression in terms of a Lorentzian path integral with its contour defined by Picard-Lefschetz theory. The wave function is real, satisfies the Wheeler-DeWitt equation and predicts an ensemble of asymptotically classical, inflationary universes with nearly-Gaussian fluctuations and with a smooth semiclassical origin.


title: Multiple Dirichlet series associated to additive and Dirichlet characters

abstract: In this article we study analytic properties of the multiple Dirichlet series associated to additive and Dirichlet characters. For the multiple Dirichlet series associated to additive characters, the meromorphic continuation is established via obtaining translation formulas satisfied by these multiple Dirichlet series. While it seems difficult to obtain such a translation formula for the multiple Dirichlet series associated to Dirichlet characters, we rely on their intrinsic connection with the multiple Dirichlet series associated to additive characters in order to investigate their analytic characteristics. We are also able to determine the exact set of singularities of the multiple Dirichlet series associated to additive characters.


title: Local and global strong solutions for SQG in bounded domains

abstract: We prove local well-posedness for the inviscid surface quasigeostrophic (SQG) equation in bounded domains of $\mathbb{R}^2$. When fractional Dirichlet Laplacian dissipation is added, global existence of strong solutions is obtained for small data for critical and supercritical cases. Global existence of strong solutions with arbitrary data is obtained in the subcritical cases.


title: Scalable Mobile Crowdsensing via Peer-to-Peer Data Sharing

abstract: Mobile crowdsensing (MCS) is a new paradigm of sensing by taking advantage of the rich embedded sensors of mobile user devices. However, the traditional server-client MCS architecture often suffers from the high operational cost on the centralized server (e.g., for storing and processing massive data), hence the poor scalability. Peer-to-peer (P2P) data sharing can effectively reduce the server's cost by leveraging the user devices' computation and storage resources. In this work, we propose a novel P2P-based MCS architecture, where the sensing data is saved and processed in user devices locally and shared among users in a P2P manner. To provide necessary incentives for users in such a system, we propose a quality-aware data sharing market, where the users who sense data can sell data to others who request data but not want to sense the data by themselves. We analyze the user behavior dynamics from the game-theoretic perspective, and characterize the existence and uniqueness of the game equilibrium. We further propose best response iterative algorithms to reach the equilibrium with provable convergence. Our simulations show that the P2P data sharing can greatly improve the social welfare, especially in the model with a high transmission cost and a low trading price.


title: GP-ILQG: Data-driven Robust Optimal Control for Uncertain Nonlinear Dynamical Systems

abstract: As we aim to control complex systems, use of a simulator in model-based reinforcement learning is becoming more common. However, it has been challenging to overcome the Reality Gap, which comes from nonlinear model bias and susceptibility to disturbance. To address these problems, we propose a novel algorithm that combines data-driven system identification approach (Gaussian Process) with a Differential-Dynamic-Programming-based robust optimal control method (Iterative Linear Quadratic Control). Our algorithm uses the simulator's model as the mean function for a Gaussian Process and learns only the difference between the simulator's prediction and actual observations, making it a natural hybrid of simulation and real-world observation. We show that our approach quickly corrects incorrect models, comes up with robust optimal controllers, and transfers its acquired model knowledge to new tasks efficiently.


title: Frequency-Multiplexed bias and readout of a 16-pixel Superconducting Nanowire Single-Photon Detector Array

abstract: We demonstrate a 16-pixel array of radio-frequency superconducting nanowire single-photon detectors with an integrated and scalable frequency-division multiplexing architecture, reducing the required bias and readout lines to a single microwave feed line. The electrical behavior of the photon-sensitive nanowires, embedded in a resonant circuit, as well as the optical performance and timing jitter of the single detectors is discussed. Besides the single pixel measurements we also demonstrate the operation of the 16-pixel array with a temporal, spatial and photon-number resolution.


title: Liouville theorems for a family of very degenerate elliptic non linear operators

abstract: We prove nonexistence results of Liouville type for nonnegative viscosity solutions of some equations involving the fully nonlinear degenerate elliptic operators ${\cal P}^\pm_k$, defined respectively as the sum of the largest and the smallest $k$ eigenvalues of the Hessian matrix. For the operator ${\cal P}^+_k$ we obtain results analogous to those which hold for the Laplace operator in space dimension $k$. Whereas, owing to the stronger degeneracy of the operator ${\cal P}^-_k$, we get totally different results.


title: Software Vulnerability Analysis Using CPE and CVE

abstract: In this paper, we analyze the Common Platform Enumeration (CPE) dictionary and the Common Vulnerabilities and Exposures (CVE) feeds. These repositories are widely used in Vulnerability Management Systems (VMSs) to check for known vulnerabilities in software products. The analysis shows, among other issues, a lack of synchronization between both datasets that can lead to incorrect results output by VMSs relying on those datasets. To deal with these problems, we developed a method that recommends to a user a prioritized list of CPE identifiers for a given software product. The user can then assign (and, if necessary, adapt) the most suitable CPE identifier to the software so that regular (e.g., daily) checks can find known vulnerabilities for this software in the CVE feeds. Our evaluation of this method shows that this interaction is indeed necessary because a fully automated CPE assignment is prone to errors due to the CPE and CVE shortcomings. We implemented an open-source VMS that employs the proposed method and published it on GitHub.


title: Exploiting the Structure via Sketched Gradient Algorithms

abstract: Sketched gradient algorithms have been recently introduced for efficiently solving the large-scale constrained Least-squares regressions. In this paper we provide novel convergence analysis for the basic method {\it Gradient Projection Classical Sketch} (GPCS) to reveal the fast linear convergence rate of GPCS towards a vicinity of the solution thanks to the intrinsic low-dimensional geometric structure of the solution prompted by constraint set. Similar to our analysis we observe computational and sketch size trade-offs in numerical experiments. Hence we justify that the combination of gradient methods and the sketching technique is a way of designing efficient algorithms which can actively exploit the low-dimensional structure to accelerate computation in large scale data regression and signal processing applications.


title: About a non-standard interpolation problem

abstract: Using algebraic methods, and motivated by the one variable case, we study a multipoint interpolation problem in the setting of several complex variables. The duality realized by the residue generator associated with an underlying Gorenstein algebra, using the Lagrange interpolation polynomial, plays a key role in the arguments.


title: Resonant single-parameter pumping in graphene

abstract: We present results for non-adiabatic single-parameter pumping in a ballistic graphene field-effect transistor. We investigate how scattering from an ac-driven top gate results in dc charge current from source to drain in an asymmetric setup caused either by geometry of the device or different doping of leads. Charge current is computed using Floquet scattering matrix approach in Landauer-B\"uttiker operator formalism. We single out two mechanisms contributing to the pumped current: Fabry-P\'erot interference in open channels and quasibound state resonant scattering through closed channels. We identify two distinct parameter regimes based on the quasibound state scattering mechanism: high and low doping of contacts compared to the frequency of the ac drive. We show that the latter regime results in a stronger peak pump current. We discuss how back gate potential and temperature dependence can be used to change the direction of the pumped current, operating the device as a switch.


title: Interplay of fast and slow dynamics in rare transition pathways: the disk-to-slab transition in the 2d Ising model

abstract: Rare transitions between long-lived stable states are often analyzed in terms of free energy landscapes computed as functions of a few collective variables. Here, using transitions between geometric phases as example, we demonstrate that the effective dynamics of a system along these variables are an essential ingredient in the description of rare events and that the static perspective provided by the free energy alone may be misleading. In particular, we investigate the disk-to-slab transition in the two-dimensional Ising model starting with a calculation of a two-dimensional free energy landscape and the distribution of committor probabilities. While at first sight it appears that the committor is incompatible with the free energy, they can be reconciled with each other using a two-dimensional Smoluchowski equation that combines the free energy landscape with state dependent diffusion coefficients. These results illustrate that dynamical information is not only required to calculate rate constants but that neglecting dynamics may also lead to an inaccurate understanding of the mechanism of a given process.


title: Electrostatic and induction effects in the solubility of water in alkanes

abstract: Experiments show that at 298~K and 1 atm pressure the transfer free energy, $\mu^{\rm ex}$, of water from its vapor to liquid normal alkanes $C_nH_{2n+2}$ ($n=5\ldots12$) is negative. Earlier it was found that with the united-atom TraPPe model for alkanes and the SPC/E model for water, one had to artificially enhance the attractive alkane-water cross interaction to capture this behavior. Here we revisit the calculation of $\mu^{\rm ex}$ using the polarizable AMOEBA and the non-polarizable Charmm General (CGenFF) forcefields. We test both the AMOEBA03 and AMOEBA14 water models; the former has been validated with the AMOEBA alkane model while the latter is a revision of AMOEBA03 to better describe liquid water. We calculate $\mu^{\rm ex}$ using the test particle method. With CGenFF, $\mu^{\rm ex}$ is positive and the error relative to experiments is about 1.5 $k_{\rm B}T$. With AMOEBA, $\mu^{\rm ex}$ is negative and deviations relative to experiments are between 0.25 $k_{\rm B}T$ (AMOEBA14) and 0.5 $k_{\rm B}T$ (AMOEBA03). Quantum chemical calculations in a continuum solvent suggest that zero point effects may account for some of the deviation. Forcefield limitations notwithstanding, electrostatic and induction effects, commonly ignored in considerations of water-alkane interactions, appear to be decisive in the solubility of water in alkanes.


title: An improved tree-graph bound

abstract: This report discusses the improved bound of the cluster expansion, recently proposed by Procacci and Yuhjtman (Lett. Math. Phys. 107, 31, 2017). Brydges and Helmuth noticed the relevance of Kruskal's algorithm, which allows to streamline and simplify the argument. These clarifications have proved useful in arXiv:1705.04385.


title: A return to eddy viscosity model for epistemic UQ in RANS closures

abstract: For the purpose of Uncertainty Quantification (UQ) of Reynolds-Averaged Navier-Stokes closures, we introduce a framework in which perturbations in the eigenvalues of the anisotropy tensor are made in order to bound a Quantity-of-Interest based on limiting states of turbulence. To make the perturbations representative of local flow features, we introduce two additional transport equations for linear combinations of these aforementioned eigenvalues. The location, magnitude and direction of the eigenvalue perturbations are now governed by the model transport equations. The general behavior of our discrepancy model is determined by two coefficients, resulting in a low-dimensional UQ problem. We will furthermore show that the behavior of the model is intuitive and rooted in the physical interpretation of misalignment between the mean strain and Reynolds stresses.


title: Probabilistic Matrix Factorization for Automated Machine Learning

abstract: In order to achieve state-of-the-art performance, modern machine learning techniques require careful data pre-processing and hyperparameter tuning. Moreover, given the ever increasing number of machine learning models being developed, model selection is becoming increasingly important. Automating the selection and tuning of machine learning pipelines consisting of data pre-processing methods and machine learning models, has long been one of the goals of the machine learning community. In this paper, we tackle this meta-learning task by combining ideas from collaborative filtering and Bayesian optimization. Using probabilistic matrix factorization techniques and acquisition functions from Bayesian optimization, we exploit experiments performed in hundreds of different datasets to guide the exploration of the space of possible pipelines. In our experiments, we show that our approach quickly identifies high-performing pipelines across a wide range of datasets, significantly outperforming the current state-of-the-art.


title: Multistage Electronic Nematic Transitions in Cuprate Superconductors: Functional-Renormalization-Group Analysis

abstract: Recently, complex phase transitions accompanied by the rotational symmetry breaking have been discovered experimentally in cuprate superconductors. To find the realized order parameters, we study various charge susceptibilities in an unbiased way, by applying the functional-renormalization-group method to the realistic $d$-$p$ Hubbard model. Without assuming the wavevector of the order parameter, we reveal that the most dominant instability is the uniform ($q = 0$) charge modulation on the $p_x$ and $p_y$ orbitals, which possesses the d-symmetry. This uniform nematic order triggers another nematic p-orbital density wave along the axial (Cu-Cu) direction at $Q_a = (\pi/2,0)$. It is predicted that uniform nematic order is driven by the spin fluctuations in the pseudogap region, and another nematic density-wave order at $q = Q_a$ is triggered by the uniform order. The predicted multistage nematic transitions are caused by the Aslamazov-Larkin-type fluctuation-exchange processes.


title: On some families of smooth affine spherical varieties of full rank

abstract: Let G be a complex connected reductive group. I. Losev has shown that a smooth affine spherical G-variety X is uniquely determined by its weight monoid, which is the set of irreducible representations of G that occur in the coordinate ring of X. In this paper we use a combinatorial characterization of the weight monoids of smooth affine spherical varieties to classify: (a) all such varieties when G is $\mathrm{SL}(2) \times \mathbb{C}^{\times}$ and (b) all such varieties for G simple which have a G-saturated weight monoid of full rank. We also use the characterization and F. Knop's classification theorem for multiplicity free Hamiltonian manifolds to give a new proof of C. Woodward's result that every reflective Delzant polytope is the moment polytope of such a manifold.


title: Construction of lepton mass matrices and TeV-scale phenomenology in the minimal left-right symmetric model

abstract: We develop a systematic procedure of constructing lepton mass matrices that satisfy all the experimental constraints in the light lepton sector of the minimal left-right symmetric model with type-I seesaw dominance. This method is unique since it is applicable to the most general cases of type-I seesaw with complex electroweak vacuum expectation values in the model. With this method, we investigate the TeV-scale phenomenology in the normal hierarchy without fine-tuning of model parameters, focusing on the charged lepton flavour violation, neutrinoless double beta decay, and electric dipole moments of charged leptons. We examine the predictions for typical ranges of associated observables such as branching ratios of rare lepton decays, and study how those experimental constraints affect the model parameter space. The most notable result is that the regions of parameter space that allow small light neutrino masses have been constrained by the present experimental bounds from charged lepton flavour violation. Furthermore, we also find that the mass of the lightest heavy neutrino should be relatively small in order to satisfy those experimental constraints.


title: Remarks on regular black holes

abstract: Recently it has been claimed by Chinaglia and Zerbini that the curvature singularity is present even in the so-called regular black hole solutions of the Einstein equations. In this brief note we show that this criticism is devoid of any physical content.


title: The intergalactic electromagnetic cascade solution for the anomalies from $\gamma$-ray blazar observations

abstract: Recent progress in very high energy (VHE, E >100 GeV) $\gamma$-ray observations, together with advances in the extragalactic background light (EBL) modelling, allows to search for new phenomena such as $\gamma$-axion-like particle ($\gamma \rightarrow$ ALP) oscillation and to explore the extragalactic magnetic field (EGMF) strength and structure. These studies are usually performed by searching for some deviation from the so-called absorption-only model, that accounts for only primary photon absorption on the EBL and adiabatic losses. In fact, there exist several indications that the absorption-only model is incomplete. We present and discuss the intergalactic electromagnetic cascade model (IECM) --- the simplest model that allows to coherently explain all known anomalies. This model has a number of robust signatures that could be searched for with present and future instruments. The IECM model may serve as a new background template, allowing to make future searches for $\gamma \rightarrow$ ALP oscillation more robust. A detailed account of our calculations is available in astro-ph/1609.01013v2 (A\&A,in print).


title: The Monster, the Baby Monster and Traces of Singular Moduli

abstract: In this note, we provide evidence for new (super) moonshines relating the Monster and the Baby monster to some weakly holomorphic weight 1/2 modular forms defined by Zagier in his work on traces of singular moduli. They are similar in spirit to the recently discovered Thompson moonshine.


title: The Mellin Formalism for Boundary CFT$_d$

abstract: We extend the Mellin representation of conformal field theory (CFT) to allow for conformal boundaries and interfaces. We consider the simplest holographic setup dual to an interface CFT - a brane filling an $AdS_{d}$ subspace of $AdS_{d+1}$ - and perform a systematic study of Witten diagrams in this setup. As a byproduct of our analysis, we show that geodesic Witten diagrams in this geometry reproduce interface CFT$_d$ conformal blocks, generalizing the analogous statement for CFTs with no defects.


title: Curiosity-driven Exploration by Self-supervised Prediction

abstract: In many real-world scenarios, rewards extrinsic to the agent are extremely sparse, or absent altogether. In such cases, curiosity can serve as an intrinsic reward signal to enable the agent to explore its environment and learn skills that might be useful later in its life. We formulate curiosity as the error in an agent's ability to predict the consequence of its own actions in a visual feature space learned by a self-supervised inverse dynamics model. Our formulation scales to high-dimensional continuous state spaces like images, bypasses the difficulties of directly predicting pixels, and, critically, ignores the aspects of the environment that cannot affect the agent. The proposed approach is evaluated in two environments: VizDoom and Super Mario Bros. Three broad settings are investigated: 1) sparse extrinsic reward, where curiosity allows for far fewer interactions with the environment to reach the goal; 2) exploration with no extrinsic reward, where curiosity pushes the agent to explore more efficiently; and 3) generalization to unseen scenarios (e.g. new levels of the same game) where the knowledge gained from earlier experience helps the agent explore new places much faster than starting from scratch. Demo video and code available at https://pathak22.github.io/noreward-rl/


title: Boundary regularity of stochastic PDEs

abstract: The boundary behaviour of solutions of stochastic PDEs with Dirichlet boundary conditions can be surprisingly - and in a sense, arbitrarily - bad: as shown by Krylov, for any $\alpha>0$ one can find a simple $1$-dimensional constant coefficient linear equation whose solution at the boundary is not $\alpha$-H\"older continuous. We obtain a positive counterpart of this: under some mild regularity assumptions on the coefficients, solutions of semilinear SPDEs on $C^1$ domains are proved to be $\alpha$-H\"older continuous up to the boundary with some $\alpha>0$.


title: Fault-tolerant quantum computation with few qubits

abstract: Reliable qubits are difficult to engineer, but standard fault-tolerance schemes use seven or more physical qubits to encode each logical qubit, with still more qubits required for error correction. The large overhead makes it hard to experiment with fault-tolerance schemes with multiple encoded qubits. The 15-qubit Hamming code protects seven encoded qubits to distance three. We give fault-tolerant procedures for applying arbitrary Clifford operations on these encoded qubits, using only two extra qubits, 17 total. In particular, individual encoded qubits within the code block can be targeted. Fault-tolerant universal computation is possible with four extra qubits, 19 total. The procedures could enable testing more sophisticated protected circuits in small-scale quantum devices. Our main technique is to use gadgets to protect gates against correlated faults. We also take advantage of special code symmetries, and use pieceable fault tolerance.


title: Maximum Selection and Ranking under Noisy Comparisons

abstract: We consider $(\epsilon,\delta)$-PAC maximum-selection and ranking for general probabilistic models whose comparisons probabilities satisfy strong stochastic transitivity and stochastic triangle inequality. Modifying the popular knockout tournament, we propose a maximum-selection algorithm that uses $\mathcal{O}\left(\frac{n}{\epsilon^2}\log \frac{1}{\delta}\right)$ comparisons, a number tight up to a constant factor. We then derive a general framework that improves the performance of many ranking algorithms, and combine it with merge sort and binary search to obtain a ranking algorithm that uses $\mathcal{O}\left(\frac{n\log n (\log \log n)^3}{\epsilon^2}\right)$ comparisons for any $\delta\ge\frac1n$, a number optimal up to a $(\log \log n)^3$ factor.


title: Evaluating XMPP Communication in IEC 61499-based Distributed Energy Applications

abstract: The IEC 61499 reference model provides an international standard developed specifically for supporting the creation of distributed event-based automation systems. Functionality is abstracted into function blocks which can be coded graphically as well as via a text-based method. As one of the design goals was the ability to support distributed control applications, communication plays a central role in the IEC 61499 specification. In order to enable the deployment of functionality to distributed platforms, these platforms need to exchange data in a variety of protocols. IEC 61499 realizes the support of these protocols via "Service Interface Function Blocks" (SIFBs). In the context of smart grids and energy applications, IEC 61499 could play an important role, as these applications require coordinating several distributed control logics. Yet, the support of grid-related protocols is a pre-condition for a wide-spread utilization of IEC 61499. The eXtensible Messaging and Presence Protocol (XMPP) on the other hand is a well-established protocol for messaging, which has recently been adopted for smart grid communication. Thus, SIFBs for XMPP facilitate distributed control applications, which use XMPP for exchanging all control relevant data, being realized with the help of IEC 61499. This paper introduces the idea of integrating XMPP into SIFBs, demonstrates the prototypical implementation in an open source IEC 61499 platform and provides an evaluation of the feasibility of the result.


title: DeepRT: deep learning for peptide retention time prediction in proteomics

abstract: Accurate predictions of peptide retention times (RT) in liquid chromatography have many applications in mass spectrometry-based proteomics. Herein, we present DeepRT, a deep learning based software for peptide retention time prediction. DeepRT automatically learns features directly from the peptide sequences using the deep convolutional Neural Network (CNN) and Recurrent Neural Network (RNN) model, which eliminates the need to use hand-crafted features or rules. After the feature learning, principal component analysis (PCA) was used for dimensionality reduction, then three conventional machine learning methods were utilized to perform modeling. Two published datasets were used to evaluate the performance of DeepRT and we demonstrate that DeepRT greatly outperforms previous state-of-the-art approaches ELUDE and GPTime.


title: On Channel State Feedback Model and Overhead in Theoretical and Practical Views

abstract: Channel state feedback plays an important role to the improvement of link performance in current wireless communication systems, and even more in the next generation. The feedback information, however, consumes the uplink bandwidth and thus generates overhead. In this paper, we investigate the impact of channel state feedback and propose an improved scheme to reduce the overhead in practical communication systems. Compared with existing schemes, we introduce a more accurate channel model to describe practical wireless channels and obtain the theoretical lower bounds of overhead for the periodical and aperiodical feedback schemes. The obtained theoretical results provide us the guidance to optimise the design of feedback systems, such as the number of bits used for quantizing channel states. We thus propose a practical feedback scheme that achieves low overhead and improved performance over currently widely used schemes such as zero holding. Simulation experiments confirm its advantages and suggest its potentially wide applications in the next generation of wireless systems.


title: Window for preferred axion models

abstract: We discuss phenomenological criteria for defining "axion windows", namely regions in the parameter space of the axion-photon coupling where realistic models live. Currently, the boundaries of this region depend on somewhat arbitrary criteria, and it would be highly desirable to specify them in terms of precise phenomenological requirements. We first focus on hadronic axion models within post-inflationary scenarios, in which the initial abundance of the new vectorlike quarks $Q$ is thermal. We classify their representations $R_Q$ by requiring that $i)$ the $Q$ are sufficiently short lived to avoid issues with long-lived strongly interacting relics, $ii)$ the theory remains weakly coupled up to the Planck scale. The more general case of multiple $R_Q$ is also studied, and the absolute upper and lower bounds on the axion-photon coupling as a function of the axion mass is identified. Pre-inflationary scenarios in which the axion decay constant remains bounded as $f_a\leq 5\cdot 10^{11}\,$GeV allow for axion-photon couplings only about 20% larger. Realistic Dine-Fischler-Srednicki-Zhitnitsky type of axion models also remain encompassed within the hadronic axion window. Some mechanisms that can allow to enhance the axion-photon coupling to values sizeably above the preferred window are discussed.


title: A Catalog of GALEX Ultraviolet Emission from Asymptotic Giant Branch Stars

abstract: We have performed a comprehensive study of the UV emission detected from AGB stars by the Galaxy Evolution Explorer (GALEX). Of the 468 AGB stars in our sample, 316 were observed by GALEX. In the NUV bandpass ($\lambda_{\rm eff} \sim 2310$z A), 179 AGB stars were detected and 137 were not detected. Only 38 AGB stars were detected in the FUV bandpass ($\lambda_{\rm eff} \sim1528$ A). We find that NUV emission is correlated with optical to near infrared emission leading to higher detection fractions among the brightest, hence closest, AGB stars. Comparing the AGB time-variable visible phased light curves to corresponding GALEX NUV phased light curves we find evidence that for some AGB stars the NUV emission varies in phase with the visible light curves. We also find evidence that the NUV emission, and possibly, the FUV emission are anti-correlated with the circumstellar envelope density. These results suggest that the origin of the GALEX-detected UV emission is an inherent characteristic of the AGB stars that can most likely be traced to a combination of photospheric and chromospheric emission. In most cases, UV detections of AGB stars are not likely to be indicative of the presence of binary companions.


title: Construction of Hamiltonians by supervised learning of energy and entanglement spectra

abstract: Correlated many-body problems ubiquitously appear in various fields of physics such as condensed matter physics, nuclear physics, and statistical physics. However, due to the interplay of the large number of degrees of freedom, it is generically impossible to treat these problems from first principles. Thus the construction of a proper model, namely effective Hamiltonian, is essential. Here, we propose a simple scheme of constructing Hamiltonians from given energy or entanglement spectra with machine learning. We apply the proposed scheme to the Hubbard model at the half-filling, and compare the obtained effective low-energy spin-1/2 model with several analytic results based on the high order perturbation theory which have been inconsistent with each other. We also show that our approach can be used to construct the entanglement Hamiltonian of a quantum many-body state from its entanglement spectrum as well. We exemplify this using the ground states of the $S=1/2$ two-leg Heisenberg ladders. We observe a qualitative difference between the entanglement Hamiltonians of the two phases (the Haldane phase and the Rung Singlet phase) of the model due to the different origin of the entanglement. In the Haldane phase, we find that the entanglement Hamiltonian is non-local by nature, and the locality can be restored by introducing the anisotropy and turning the system into the large-$D$ phase. Possible applications to the study of strongly-correlated systems and the model construction from experimental data are discussed.


title: EMERGE - An empirical model for the formation of galaxies since $z\sim10$

abstract: We present EMERGE, an Empirical ModEl for the foRmation of GalaxiEs, describing the evolution of individual galaxies in large volumes from $z\sim10$ to the present day. We assign a star formation rate to each dark matter halo based on its growth rate, which specifies how much baryonic material becomes available, and the instantaneous baryon conversion efficiency, which determines how efficiently this material is converted to stars, thereby capturing the baryonic physics. Satellites are quenched following the delayed-then-rapid model, and they are tidally disrupted once their subhalo has lost a significant fraction of its mass. The model is constrained with observed data extending out to high redshift. The empirical relations are very flexible, and the model complexity is increased only if required by the data, assessed by several model selection statistics. We find that for the same final halo mass galaxies can have very different star formation histories. Nevertheless, the average star formation and accretion rates are in good agreement with models following an abundance matching strategy. Galaxies that are quenched at $z=0$ typically have a higher peak star formation rate compared to their star-forming counterparts. The accretion of stars can dominate the total mass of massive galaxies, but is insignificant for low-mass systems, independent of star-formation activity. EMERGE predicts stellar-to-halo mass ratios for individual galaxies and introduces scatter self-consistently. We find that at fixed halo mass, passive galaxies have a higher stellar mass on average. The intra-cluster-mass in massive haloes can be up to 8 times larger than the mass of the central galaxy. Clustering for star-forming and quenched galaxies is in good agreement with observational constraints, indicating a realistic assignment of galaxies to haloes.


title: Efficient adiabatic hydrodynamical simulations of the high-redshift intergalactic medium

abstract: We present a post-processing tool for GADGET-2 adiabatic simulations to model various observed properties of the Ly$\alpha$ forest at $2.5 \leq z \leq 4$ that enables an efficient parameter estimation. In particular, we model the thermal and ionization histories that are not computed self-consistently by default in GADGET-2. We capture the effect of pressure smoothing by running GADGET-2 at an elevated temperature floor and using an appropriate smoothing kernel. We validate our procedure by comparing different statistics derived from our method with those derived using self-consistent simulations with GADGET-3. These statistics are: line of sight density field power spectrum, flux probability distribution function, flux power spectrum, wavelet statistics, curvature statistics, HI column density (${\rm N_{HI}}$) distribution function, linewidth ($b$) distribution and $b$ versus $\log {\rm N_{HI}}$ scatter. For the temperature floor of $10^4$ K and typical signal-to-noise of 25, the results agree well within 20 percent of the self-consistent GADGET-3 simulation. However, this difference is smaller than the expected $1\sigma$ sample variance for an absorption path length of $\sim 5.35$ at $z=3$. Moreover for a given cosmology, we gain a factor of $\sim N$ in computing time for modelling the intergalactic medium under $N \gg 1$ different thermal histories. In addition, our method allows us to simulate the non-equilibrium evolution of thermal and ionization state of the gas and include heating due to non-standard sources like cosmic rays and high energy $\gamma$-rays from Blazars.


title: The Energy Distribution of Subjets and the Jet Shape

abstract: We present a framework that describes the energy distribution of subjets of radius $r$ within a jet of radius $R$. We consider both an inclusive sample of subjets as well as subjets centered around a predetermined axis, from which the jet shape can be obtained. For $r \ll R$ we factorize the physics at angular scales $r$ and $R$ to resum the logarithms of $r/R$. For central subjets, we consider both the standard jet axis and the winner-take-all axis, which involve double and single logarithms of $r/R$, respectively. All relevant one-loop matching coefficients are given, and an inconsistency in some previous results for cone jets is resolved. Our results for the standard jet shape differ from previous calculations at next-to-leading logarithmic order, because we account for the recoil of the standard jet axis due to soft radiation. Numerical results are presented for an inclusive subjet sample for $pp \to {\rm jet}+X$ at next-to-leading order plus leading logarithmic order.


title: Detection of virial shocks in stacked Fermi-LAT clusters

abstract: Galaxy clusters are thought to grow by accreting mass through large-scale, strong, yet elusive, virial shocks. Such a shock is expected to accelerate relativistic electrons, thus generating a spectrally-flat leptonic virial-ring. However, until now, only the nearby Coma cluster has shown evidence for a $\gamma$-ray virial ring. We stack Fermi-LAT data for the 112 most massive, high latitude, extended clusters, enhancing the ring sensitivity by rescaling clusters to their virial radii and utilizing the expected flat energy spectrum. In addition to a central unresolved, hard signal (detected at the $\sim 5.8\sigma$ confidence level), probably dominated by AGN, we identify (at the $5.8\sigma$ confidence level) a bright, spectrally-flat $\gamma$-ray ring at the expected virial shock position. The ring signal implies that the shock deposits $\sim 0.6\%$ (with an interpretation uncertainty factor $\sim2$) of the thermal energy in relativistic electrons over a Hubble time. This result, consistent with the Coma signal, validates and calibrates the virial shock model, and indicates that the cumulative emission from such shocks significantly contributes to the diffuse extragalactic $\gamma$-ray and low-frequency radio backgrounds.


title: Is a top-heavy initial mass function needed to reproduce the submillimeter galaxy number counts?

abstract: Matching the number counts and redshift distribution of submillimeter galaxies (SMGs) without invoking modifications to the initial mass function (IMF) has proved challenging for semi-analytic models (SAMs) of galaxy formation. We adopt a previously developed SAM that is constrained to match the z=0 galaxy stellar mass function and makes various predictions that agree well with observational constraints; we do not recalibrate the SAM for this work. We implement three prescriptions to predict the submillimeter flux densities of the model galaxies; two depend solely on star formation rate, whereas the other also depends on the dust mass. By comparing the predictions of the models, we find that taking into account the dust mass, which affects the dust temperature and thus influences the far-infrared spectral energy distribution, is crucial for matching the number counts and redshift distribution of SMGs. Moreover, despite using a standard IMF, our model can match the observed SMG number counts and redshift distribution reasonably well, which contradicts the conclusions of some previous studies that a top-heavy IMF, in addition to taking into account the effect of dust mass, is needed to match these observations. Although we have not identified the key ingredient that is responsible for our model matching the observed SMG number counts and redshift distribution without IMF variation- which is challenging given the different prescriptions for physical processes employed in the SAMs of interest-our results demonstrate that in SAMs, IMF variation is degenerate with other physical processes, such as stellar feedback.


title: Quantized charge transport in chiral Majorana edge modes

abstract: Majorana fermions can be realized as quasiparticles in topological superconductors, with potential applications in topological quantum computing. Recently, lattices of magnetic adatoms deposited on the surface of s-wave superconductors - Shiba lattices - have been proposed as a new platform for topological superconductivity. These systems possess the great advantage that they are accessible via scanning-probe techniques, and thus enable the local manipulation and detection of Majorana modes. Realizing quantum bits on this basis will rely on the creation of nanoscopic Shiba lattices, so-called Shiba islands. Here, we demonstrate that the topological Majorana edge modes of such islands display universal electronic and transport properties. Most remarkably, these Majorana modes possess a quantized charge conductance that is proportional to the topological Chern number, C, and carry a supercurrent whose chirality reflects the sign of C. These results establish nanoscopic Shiba islands as promising components in future topology-based devices.


title: Dissipative Landau-Zener problem and thermally assisted Quantum Annealing

abstract: We revisit the issue of thermally assisted quantum annealing by a detailed study of the dissipative Landau-Zener problem in presence of a Caldeira-Leggett bath of harmonic oscillators, using both a weak-coupling quantum master equation and a quasi-adiabatic path-integral approach. Building on the known zero-temperature exact results (Wubs et al., PRL 97, 200404 (2006)), we show that a finite temperature bath can have a beneficial effect on the ground-state probability only if it couples also to a spin-direction that is transverse with respect to the driving field, while no improvement is obtained for the more commonly studied purely longitudinal coupling. In particular, we also highlight that, for a transverse coupling, raising the bath temperature further improves the ground-state probability in the fast-driving regime. We discuss the relevance of these findings for the current quantum-annealing flux qubit chips.


title: Sub-Riemannian interpolation inequalities

abstract: We prove that ideal sub-Riemannian manifolds (i.e., admitting no non-trivial abnormal minimizers) support interpolation inequalities for optimal transport. A key role is played by sub-Riemannian Jacobi fields and distortion coefficients, whose properties are remarkably different with respect to the Riemannian case. As a byproduct, we characterize the cut locus as the set of points where the squared sub-Riemannian distance fails to be semiconvex, answering to a question raised by Figalli and Rifford in [Geom. Funct. Anal. (2010) 20: 124]. As an application, we deduce sharp and intrinsic Borell-Brascamp-Lieb and geodesic Brunn-Minkowski inequalities in the aforementioned setting. For the case of the Heisenberg group, we recover in an intrinsic way the results recently obtained by Balogh, Krist\'aly and Sipos in [Calc. Var. PDE (2018) 57: 61], and we extend them to the class of generalized H-type Carnot groups. Our results do not require the distribution to have constant rank, yielding for the particular case of the Grushin plane a sharp measure contraction property and a sharp Brunn-Minkowski inequality.


title: Emergence of Chiral Spin Liquids via Quantum Melting of Non-Coplanar Magnetic Orders

abstract: Quantum spin liquids (QSLs) are long-range entangled states of quantum magnets which lie beyond the Landau paradigm of classifying phases of matter via broken symmetries. A physical route to arriving at QSLs is via frustration-induced quantum melting of ordered states such as valence bond crystals or magnetic orders. Here, we show, using extensive exact diagonalization (ED) and density-matrix renormalization group (DMRG) studies of concrete $SU(2)$ invariant spin models on honeycomb, triangular and square lattices, that chiral spin liquids (CSLs) emerge as descendants of triple-$Q$ spin crystals with tetrahedral magnetic order and a large scalar spin chirality. Such ordered-to-CSL melting transitions may yield lattice realizations of effective Chern-Simons-Higgs field theories. Our work provides a distinct unifying perspective on the emergence of CSLs, and suggests that materials with magnetic skyrmion crystal order might provide a good starting point to search for CSLs.


title: Rapid growth of black holes accompanied with hot or warm outflows exposed to anisotropic super-Eddington radiation

abstract: We perform two-dimensional radiation hydrodynamical simulations of accretion flows onto a black hole (BH) with a mass of $10^3\leq M_{\rm BH}/M_{\odot} \lesssim 10^6$ in order to study rapid growth of BHs in the early Universe. For spherically symmetric flows, hyper-Eddington accretion onto the BH from outside the Bondi radius can occur unimpeded by radiation feedback only when the BH mass is higher than $\simeq 10^4~M_{\odot}(n_\infty/10^5~{\rm cm}^{-3})^{-1}(T_\infty/10^4~{\rm K})^{3/2}$, where $n_\infty$ and $T_\infty$ are the density and temperature of ambient gas. Here, we study the properties of accretion flows exposed to anisotropic radiation from a nuclear accretion disk with a luminosity higher than the Eddington value ($L_{\rm Edd}$) due to collimation toward the bipolar directions. We find that, unlike the spherically symmetric case, even less massive BHs with $M_{\rm BH} < 10^4~M_{\odot}$ can be fed by surrounding gas at high accretion rates of $\gtrsim L_{\rm Edd}/c^2$ through the equatorial plane, while ionized regions expand to the polar directions producing hot outflows with $T\sim 10^5$K. For more massive BHs with $M_{\rm BH}\gtrsim 5\times 10^5~M_{\odot}$, neutral gas through the equatorial plane totally covers the central radiating region due to the non-radial gas motions, and thus the emergent radiation in all directions is blocked. Because of efficient recombination by hydrogen, the entire flow results in neutral and warm gas with $T \simeq 8000~{\rm K}$ . The central BH is fed through the equator at the averaged rate of $\sim 5\times 10^4~L_{\rm Edd}/c^2$, which corresponds to $\sim 50~\%$ of the inflow rate from the Bondi radius. Moreover, radiation momentum absorbed by neutral hydrogen produces warm outflows toward the bipolar directions at $\sim 30~\%$ of the BH feeding rate and with a typical velocity of $\simeq 50~{\rm km~s}^{-1}$.


title: An extremization principle for the entropy of rotating BPS black holes in AdS$_5$

abstract: We show that the Bekenstein-Hawking entropy of a class of BPS electrically charged rotating black holes in AdS$_5\times S^5$ can be obtained by a simple extremization principle. We expect that this extremization corresponds to the attractor mechanism for BPS rotating black holes in five-dimensional gauged supergravity, which is still unknown. The expression to be extremized has a suggestive resemblance to anomaly polynomials and the supersymmetric Casimir energy recently studied for $\mathcal{N}=4$ super Yang-Mills.


title: Orbital superconductivity, defects and pinned nematic fluctuations in the doped iron chalcogenide FeSe$_{0.45}$Te$_{0.55}$

abstract: We demonstrate that the differential conductance, $dI/dV$, measured via spectroscopic imaging scanning tunneling microscopy in the doped iron chalcogenide FeSe$_{0.45}$Te$_{0.55}$, possesses a series of characteristic features that allow one to extract the orbital structure of the superconducting gaps. This yields nearly isotropic superconducting gaps on the two hole-like Fermi surfaces, and a strongly anisotropic gap on the electron-like Fermi surface. Moreover, we show that the pinning of nematic fluctuations by defects can give rise to a dumbbell-like spatial structure of the induced impurity bound states, and explains the related $C_2$-symmetry in the Fourier transformed differential conductance.


title: Decay of Ultralight Axion Condensates

abstract: Axion particles can form macroscopic condensates, whose size can be galactic in scale for models with very small axion masses $m\sim10^{-22}$ eV, and which are sometimes referred to under the name of Fuzzy Dark Matter. Many analyses of these condensates are done in the non-interacting limit, due to the weakness of the self-interaction coupling of axions. We investigate here how certain results change upon inclusion of these interactions, finding a decreased maximum mass and a modified mass-radius relationship. Further, these condensates are, in general, unstable to decay through number-changing interactions. We analyze the stability of galaxy-sized condensates of axion-like particles, and sketch the parameter space of stable configurations as a function of a binding energy parameter. We find a strong lower bound on the size of Fuzzy Dark Matter condensates which are stable to decay, with lifetimes longer than the age of the universe.


title: The redshift distribution of cosmological samples: a forward modeling approach

abstract: Determining the redshift distribution $n(z)$ of galaxy samples is essential for several cosmological probes including weak lensing. For imaging surveys, this is usually done using photometric redshifts estimated on an object-by-object basis. We present a new approach for directly measuring the global $n(z)$ of cosmological galaxy samples, including uncertainties, using forward modeling. Our method relies on image simulations produced using UFig (Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation) within the $MCCL$ (Monte-Carlo Control Loops) framework. The galaxy population is modeled using parametric forms for the luminosity functions, spectral energy distributions, sizes and radial profiles of both blue and red galaxies. We apply exactly the same analysis to the real data and to the simulated images, which also include instrumental and observational effects. By adjusting the parameters of the simulations, we derive a set of acceptable models that are statistically consistent with the data. We then apply the same cuts to the simulations that were used to construct the target galaxy sample in the real data. The redshifts of the galaxies in the resulting simulated samples yield a set of $n(z)$ distributions for the acceptable models. We demonstrate the method by determining $n(z)$ for a cosmic shear like galaxy sample from the 4-band Subaru Suprime-Cam data in the COSMOS field. We also complement this imaging data with a spectroscopic calibration sample from the VVDS survey. We compare our resulting posterior $n(z)$ distributions to the one derived from photometric redshifts estimated using 36 photometric bands in COSMOS and find good agreement. This offers good prospects for applying our approach to current and future large imaging surveys.


title: Generalised Proca Theories

abstract: In this work we will summarise the recent progress made in constructing consistent theories for a massive vector field with derivative self-interactions. The construction is such that only the three desired polarisations of the Proca field propagate. We apply a systematic construction of the interactions by using the anti-symmetric Levi-Civita tensors. This finite family of allowed derivative self-interactions can be also obtained from the decoupling limit by imposing that the St\"uckelberg field only contains second derivatives both with itself and with the transverse modes. These interactions can be generalised to curved backgrounds, which relies on the presence of non-minimal couplings and constitutes a general family of vector-tensor interactions. We discuss also some extensions of these interactions by alleviating the restriction of second order nature of the equations or by imposing global symmetries. We will also comment on their interesting cosmological applications.


title: Neutrino Masses and Absence of Flavor Changing Interactions in the 2HDM from Gauge Principles

abstract: We propose several Two Higgs Doublet Models with the addition of an Abelian gauge group which free the usual framework from flavor changing neutral interactions and explain neutrino masses through the seesaw mechanism. We discuss the kinetic and mass-mixing gripping phenomenology which encompass several constraints coming from atomic parity violation, the muon anomalous magnetic moment, rare meson decays, Higgs physics, LEP precision data, neutrino-electron scattering, low energy accelerators and LHC probes.


title: Entanglement enabled telescopic arrays in the presence of decoherence

abstract: We consider quantum enhancement of direct-detection interferometric measurements to increase telescope resolution. We propose a protocol of measuring interferometric visibility function using imperfectly entangled states shared between remote telescopes. We show how errors in visibility measurement, and in turn, errors in intensity distribution of a distant object depend on the degree of entanglement of the shared quantum resource. We determine that these errors are sufficiently small over a wide range of resource states which makes our technique feasible in practical environments.


title: Intertwined Orders in Holography: Pair and Charge Density Waves

abstract: Building on [1], we examine a holographic model in which a U(1) symmetry and translational invariance are broken spontaneously at the same time. The symmetry breaking is realized through the St\"{u}ckelberg mechanism, and leads to a scalar condensate and a charge density which are spatially modulated and exhibit unidirectional stripe order. Depending on the choice of parameters, the oscillations of the scalar condensate can average out to zero, with a frequency which is half of that of the charge density. In this case the system realizes some of the key features of pair density wave order. The model also admits a phase with co-existing superconducting and charge density wave orders, in which the scalar condensate has a uniform component. In our construction the various orders are intertwined with each other and have a common origin. The fully backreacted geometry is computed numerically, including for the case in which the theory contains axions. The latter can be added to explicitly break translational symmetry and mimic lattice-type effects.


title: Optimal Rates and Tradeoffs in Multiple Testing

abstract: Multiple hypothesis testing is a central topic in statistics, but despite abundant work on the false discovery rate (FDR) and the corresponding Type-II error concept known as the false non-discovery rate (FNR), a fine-grained understanding of the fundamental limits of multiple testing has not been developed. Our main contribution is to derive a precise non-asymptotic tradeoff between FNR and FDR for a variant of the generalized Gaussian sequence model. Our analysis is flexible enough to permit analyses of settings where the problem parameters vary with the number of hypotheses $n$, including various sparse and dense regimes (with $o(n)$ and $\mathcal{O}(n)$ signals). Moreover, we prove that the Benjamini-Hochberg algorithm as well as the Barber-Cand\`{e}s algorithm are both rate-optimal up to constants across these regimes.


title: Carving Out the End of the World or (Superconformal Bootstrap in Six Dimensions)

abstract: We bootstrap ${\cal N}=(1,0)$ superconformal field theories in six dimensions, by analyzing the four-point function of flavor current multiplets. Assuming $E_8$ flavor group, we present universal bounds on the central charge $C_T$ and the flavor central charge $C_J$. Based on the numerical data, we conjecture that the rank-one E-string theory saturates the universal lower bound on $C_J$, and numerically determine the spectrum of long multiplets in the rank-one E-string theory. We comment on the possibility of solving the higher-rank E-string theories by bootstrap and thereby probing M-theory on AdS${}_7\times{\rm S}^4$/$\mathbb{Z}_2$.


title: A class of exponential neighbourhoods for the quadratic travelling salesman problem

abstract: The Quadratic Travelling Salesman Problem (QTSP) is to find a least-cost Hamiltonian cycle in an edge-weighted graph, where costs are defined on all pairs of edges such that each edge in the pair is contained in the Hamiltonian cycle. This is a more general version than the one that appears in the literature as the QTSP, denoted here as the \emph{adjacent quadratic} TSP, which only considers costs for pairs of adjacent edges. Major directions of research work on the linear TSP include exact algorithms, heuristics, approximation algorithms, polynomially solvable special cases and exponential neighbourhoods among others. In this paper we explore the complexity of searching exponential neighbourhoods for QTSP, the fixed-rank QTSP, and the adjacent quadratic TSP. The fixed-rank QTSP is introduced as a restricted version of the QTSP where the cost matrix has fixed rank $p$. It is shown that fixed-rank QTSP is solvable in pseudopolynomial time and admits an FPTAS for each of the special cases studied, except for the case of matching edge ejection tours. The adjacent quadratic TSP is shown to be polynomially-solvable in many of the cases for which the linear TSP is polynomially-solvable. Interestingly, optimizing over the matching edge ejection tour neighbourhood is shown to be pseudopolynomial for the rank 1 case without a linear term in the objective function, but NP-hard for the adjacent quadratic TSP case. We also show that the quadratic shortest path problem on an acyclic digraph can be solved in pseudopolynomial time and by an FPTAS when the rank of the associated cost matrix is fixed.


title: Probabilistically Safe Policy Transfer

abstract: Although learning-based methods have great potential for robotics, one concern is that a robot that updates its parameters might cause large amounts of damage before it learns the optimal policy. We formalize the idea of safe learning in a probabilistic sense by defining an optimization problem: we desire to maximize the expected return while keeping the expected damage below a given safety limit. We study this optimization for the case of a robot manipulator with safety-based torque limits. We would like to ensure that the damage constraint is maintained at every step of the optimization and not just at convergence. To achieve this aim, we introduce a novel method which predicts how modifying the torque limit, as well as how updating the policy parameters, might affect the robot's safety. We show through a number of experiments that our approach allows the robot to improve its performance while ensuring that the expected damage constraint is not violated during the learning process.


title: Spin and tunneling dynamics in an asymmetrical double quantum dot with spin - orbit coupling

abstract: In this article we study the spin and tunneling dynamics as a function of magnetic field in a one-dimensional GaAs double quantum dot with both the Dresselhaus and Rashba spin-orbit coupling. In particular we consider different spatial widths for the spin-up and spin-down electronic states. We find that the spin dynamics is a superposition of slow as well as fast Rabi oscillations. It is found that the Rashba interaction strength as well as the external magnetic field strongly modifies the slow Rabi oscillations which is particularly useful for single qubit manipulation for possible quantum computer applications.


title: Learning Probabilistic Programs Using Backpropagation

abstract: Probabilistic modeling enables combining domain knowledge with learning from data, thereby supporting learning from fewer training instances than purely data-driven methods. However, learning probabilistic models is difficult and has not achieved the level of performance of methods such as deep neural networks on many tasks. In this paper, we attempt to address this issue by presenting a method for learning the parameters of a probabilistic program using backpropagation. Our approach opens the possibility to building deep probabilistic programming models that are trained in a similar way to neural networks.


title: Quantum fluctuation theorems, contextuality and work quasi-probabilities

abstract: We discuss the role of contextuality within quantum fluctuation theorems, in the light of a recent no-go result by Perarnau \emph{et al}. We show that any fluctuation theorem reproducing the two-point-measurement scheme for classical states either admits a notion of work quasi-probability or fails to describe protocols exhibiting contextuality. Conversely, we describe a protocol that smoothly interpolates between the two-point measurement work distribution for projective measurements and Allahverdyan's work quasi-probability for weak measurements, and show that the negativity of the latter is a direct signature of contextuality.


title: Constraining the contribution of active galactic nuclei to reionisation

abstract: Recent results have suggested that active galactic nuclei (AGN) could provide enough photons to reionise the Universe. We assess the viability of this scenario using a semi-numerical framework for modeling reionisation, to which we add a quasar contribution by constructing a Quasar Halo Occupation Distribution (QHOD) based on Giallongo et al. observations. Assuming a constant QHOD, we find that an AGN-only model cannot simultaneously match observations of the optical depth $\tau_e$, neutral fraction, and ionising emissivity. Such a model predicts $\tau_e$ too low by $\sim 2\sigma$ relative to Planck constraints, and reionises the Universe at $z\lesssim 5$. Arbitrarily increasing the AGN emissivity to match these results yields a strong mismatch with the observed ionising emissivity at $z\sim 5$. If we instead assume a redshift-independent AGN luminosity function yielding an emissivity evolution like that assumed in Madau & Haardt model, then we can match $\tau_e$ albeit with late reionisation, however such evolution is inconsistent with observations at $z\sim 4-6$ and poorly motivated physically. These results arise because AGN are more biased towards massive halos than typical reionising galaxies, resulting in stronger clustering and later formation times. AGN-dominated models produce larger ionising bubbles that are reflected in $\sim\times 2$ more 21cm power on all scales. A model with equal parts galaxies and AGN contribution is still (barely) consistent with observations, but could be distinguished using next-generation 21cm experiments HERA and SKA-low. We conclude that, even with recent claims of more faint AGN than previously thought, AGN are highly unlikely to dominate the ionising photon budget for reionisation.


title: Gauge invariant one-loop corrections to Higgs boson couplings in non-minimal Higgs models

abstract: We comprehensively evaluate renormalized Higgs boson couplings at one-loop level in non-minimal Higgs models such as the Higgs Singlet Model (HSM) and the four types of Two Higgs Doublet Models (THDMs) with a softly-broken $Z_2$ symmetry. The renormalization calculation is performed in the on-shell scheme improved by using the pinch technique to eliminate the gauge dependence in the renormalized couplings. We first review the pinch technique for scalar boson two-point functions in the Standard Model (SM), the HSM and the THDMs. We then discuss the difference in the results of the renormalized Higgs boson couplings between the improved on-shell scheme and the ordinal one with a gauge dependence appearing in mixing parameters of scalar bosons. Finally, we widely investigate how we can identify the HSM and the THDMs focusing on the pattern of deviations in the renormalized Higgs boson couplings from predictions in the SM.


title: Aggregation-fragmentation-diffusion model for trail dynamics

abstract: We investigate statistical properties of trails formed by a random process incorporating aggregation, fragmentation, and diffusion. In this stochastic process, which takes place in one spatial dimension, two neighboring trails may combine to form a larger one and also, one trail may split into two. In addition, trails move diffusively. The model is defined by two parameters which quantify the fragmentation rate and the fragment size. In the long-time limit, the system reaches a steady state, and our focus is the limiting distribution of trail weights. We find that the density of trail weight has power-law tail $P(w) \sim w^{-\gamma}$ for small weight $w$. We obtain the exponent $\gamma$ analytically, and find that it varies continuously with the two model parameters. The exponent $\gamma$ can be positive or negative, so that in one range of parameters small-weight tails are abundant, and in the complementary range, they are rare.


title: General Bootstrap Equations in 4D CFTs

abstract: We provide a framework for generic 4D conformal bootstrap computations. It is based on the unification of two independent approaches, the covariant (embedding) formalism and the non-covariant (conformal frame) formalism. We construct their main ingredients (tensor structures and differential operators) and establish a precise connection between them. We supplement the discussion by additional details like classification of tensor structures of n-point functions, normalization of 2-point functions and seed conformal blocks, Casimir differential operators and treatment of conserved operators and permutation symmetries. Finally, we implement our framework in a Mathematica package and make it freely available.


title: Radiation Pressure Driven Ion Weibel Instability and Collisionless Shocks

abstract: The Weibel instability from counterstreaming plasma flows is a basic process highly relevant for collisionless shock formation in astrophysics. In this Letter we investigate, via two- and three- dimensional simulations, suitable configurations for laboratory investigations of the ion Weibel instability (IWI) driven by a fast quasi-neutral plasma flow launched into the target via the radiation pressure of an ultra-high-intensity (UHI) laser pulse ('Hole-Boring' process). The use of S-polarized light at oblique incidence is found to be an optimal configuration for driving IWI, as it prevents the development of surface rippling observed at normal incidence, that would lead to strong electron heating and favors competing instabilities. Conditions for the evolution of IWI into a collisionless shock are also investigated.


title: A statistical physics approach to learning curves for the Inverse Ising problem

abstract: Using methods of statistical physics, we analyse the error of learning couplings in large Ising models from independent data (the inverse Ising problem). We concentrate on learning based on local cost functions, such as the pseudo-likelihood method for which the couplings are inferred independently for each spin. Assuming that the data are generated from a true Ising model, we compute the reconstruction error of the couplings using a combination of the replica method with the cavity approach for densely connected systems. We show that an explicit estimator based on a quadratic cost function achieves minimal reconstruction error, but requires the length of the true coupling vector as prior knowledge. A simple mean field estimator of the couplings which does not need such knowledge is asymptotically optimal, i.e. when the number of observations is much large than the number of spins. Comparison of the theory with numerical simulations shows excellent agreement for data generated from two models with random couplings in the high temperature region: a model with independent couplings (Sherrington-Kirkpatrick model), and a model where the matrix of couplings has a Wishart distribution.


title: UVI colour gradients of 0.4<z<1.4 star-forming main sequence galaxies in CANDELS: dust extinction and star formation profiles

abstract: This paper uses radial colour profiles to infer the distributions of dust, gas and star formation in z=0.4-1.4 star-forming main sequence galaxies. We start with the standard UVJ-based method to estimate dust extinction and specific star formation rate (sSFR). By replacing J with I band, a new calibration method suitable for use with ACS+WFC3 data is created (i.e. UVI diagram). Using a multi-wavelength multi-aperture photometry catalogue based on CANDELS, UVI colour profiles of 1328 galaxies are stacked in stellar mass and redshift bins. The resulting colour gradients, covering a radial range of 0.2--2.0 effective radii, increase strongly with galaxy mass and with global $A_V$. Colour gradient directions are nearly parallel to the Calzetti extinction vector, indicating that dust plays a more important role than stellar population variations. With our calibration, the resulting $A_V$ profiles fall much more slowly than stellar mass profiles over the measured radial range. sSFR gradients are nearly flat without central quenching signatures, except for $M_*>10^{10.5} M_{\odot}$, where central declines of 20--25 per cent are observed. Both sets of profiles agree well with previous radial sSFR and (continuum) $A_V$ measurements. They are also consistent with the sSFR profiles and, if assuming a radially constant gas-to-dust ratio, gas profiles in recent hydrodynamic models. We finally discuss the striking findings that SFR scales with stellar mass density in the inner parts of galaxies, and that dust content is high in the outer parts despite low stellar-mass surface densities there.


title: The Three-Dimensional Expansion of the Ejecta from Tycho's Supernova Remnant

abstract: We present the first three-dimensional measurements of the velocity of various ejecta knots in Tycho's supernova remnant, known to result from a Type Ia explosion. Chandra X-ray observations over a 12-year baseline from 2003 to 2015 allow us to measure the proper motion of nearly 60 "tufts" of Si-rich ejecta, giving us the velocity in the plane of the sky. For the line of sight velocity, we use two different methods: a non-equilibrium ionization model fit to the strong Si and S lines in the 1.2-2.8 keV regime, and a fit consisting of a series of Gaussian lines. These methods give consistent results, allowing us to determine the red or blue shift of each of the knots. Assuming a distance of 3.5 kpc, we find total velocities that range from 2400 to 6600 km s$^{-1}$, with a mean of 4430 km s$^{-1}$. We find several regions where the ejecta knots have overtaken the forward shock. These regions have proper motions in excess of 6000 km s$^{-1}$. Some Type Ia supernova explosion models predict a velocity asymmetry in the ejecta. We find no such velocity asymmetries in Tycho, and discuss our findings in light of various explosion models, favoring those delayed detonation models with relatively vigorous and symmetrical deflagrations. Finally, we compare measurements with models of the remnant's evolution that include both smooth and clumpy ejecta profiles, finding that both ejecta profiles can be accommodated by the observations.


title: Community structure of copper supply networks in the prehistoric Balkans: An independent evaluation of the archaeological record from the 7th to the 4th millennium BC

abstract: Complex networks analyses of many physical, biological and social phenomena show remarkable structural regularities, yet, their application in studying human past interaction remains underdeveloped. Here, we present an innovative method for identifying community structures in the archaeological record that allow for independent evaluation of the copper using societies in the Balkans, from c. 6200 to c. 3200 BC. We achieve this by exploring modularity of networked systems of these societies across an estimated 3000 years. We employ chemical data of copper-based objects from 79 archaeological sites as the independent variable for detecting most densely interconnected sets of nodes with a modularity maximization method. Our results reveal three dominant modular structures across the entire period, which exhibit strong spatial and temporal significance. We interpret patterns of copper supply among prehistoric societies as reflective of social relations, which emerge as equally important as physical proximity. Although designed on a variable isolated from any archaeological and spatiotemporal information, our method provides archaeologically and spatiotemporally meaningful results. It produces models of human interaction and cooperation that can be evaluated independently of established archaeological systematics, and can find wide application on any quantitative data from archaeological and historical record.


title: Discontinuous Galerkin algorithms for fully kinetic plasmas

abstract: We present a new algorithm for the discretization of the Vlasov-Maxwell system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin finite element method for the spatial discretization, we obtain a high order accurate solution for the plasma's distribution function. Time stepping for the distribution function is done explicitly with a third order strong-stability preserving Runge-Kutta method. Since the Vlasov equation in the Vlasov-Maxwell system is a high dimensional transport equation, up to six dimensions plus time, we take special care to note various features we have implemented to reduce the cost while maintaining the integrity of the solution, including the use of a reduced high-order basis set. A series of benchmarks, from simple wave and shock calculations, to a five dimensional turbulence simulation, are presented to verify the efficacy of our set of numerical methods, as well as demonstrate the power of the implemented features.

categories: physics.plasm-ph astro-ph.SR math.NA physics.comp-ph physics.space-ph


title: A criticality result for polycycles in a family of quadratic reversible centers

abstract: We consider the family of dehomogenized Loud's centers $X_{\mu}=y(x-1)\partial_x+(x+Dx^2+Fy^2)\partial_y,$ where $\mu=(D,F)\in\mathbb{R}^2,$ and we study the number of critical periodic orbits that emerge or dissapear from the polycycle at the boundary of the period annulus. This number is defined exactly the same way as the well-known notion of cyclicity of a limit periodic set and we call it criticality. The previous results on the issue for the family $\{X_{\mu},\mu\in\mathbb{R}^2\}$ distinguish between parameters with criticality equal to zero (regular parameters) and those with criticality greater than zero (bifurcation parameters). A challenging problem not tackled so far is the computation of the criticality of the bifurcation parameters, which form a set $\Gamma_{B}$ of codimension 1 in $\mathbb{R}^2$. In the present paper we succeed in proving that a subset of $\Gamma_{B}$ has criticality equal to one.


title: The nuclear force imprints revealed on the elastic scattering of protons with $^{10}$C

abstract: How does nature hold together protons and neutrons to form the wide variety of complex nuclei in the universe? Describing many-nucleon systems from the fundamental theory of quantum chromodynamics has been the greatest challenge in answering this question. The chiral effective field theory description of the nuclear force now makes this possible but requires certain parameters that are not uniquely determined. Defining the nuclear force needs identification of observables sensitive to the different parametrizations. From a measurement of proton elastic scattering on $^{10}$C at TRIUMF and ab initio nuclear reaction calculations we show that the shape and magnitude of the measured differential cross section is strongly sensitive to the nuclear force prescription.


title: Isolating the chiral magnetic effect from backgrounds by pair invariant mass

abstract: Topological gluon configurations in quantum chromodynamics induce quark chirality imbalance in local domains, which can result in the chiral magnetic effect (CME)--an electric charge separation along a strong magnetic field. Experimental searches for the CME in relativistic heavy ion collisions via the charge-dependent azimuthal correlator ($\Delta\gamma$) suffer from large backgrounds arising from particle correlations (e.g. due to resonance decays) coupled with the elliptic anisotropy. We propose differential measurements of the $\Delta\gamma$ as a function of the pair invariant mass ($m_{\rm inv}$), by restricting to high $m_{\rm inv}$ thus relatively background free, and by studying the $m_{\rm inv}$ dependence to separate the possible CME signal from backgrounds. We demonstrate by model studies the feasibility and effectiveness of such measurements for the CME search.


title: Simulation of electrochemical processes during oxygen evolution on $\mathrm{Pb-MnO_2}$ composite electrodes

abstract: The geometric properties of $\mathrm{Pb-MnO_2}$ composite electrodes are studied, and a general formula is presented for the length of the triple phase boundary (TPB) on two dimensional (2D) composite electrodes using sphere packing and cutting simulations. The difference in the geometrical properties of 2D (or compact) and 3D (or porous) electrodes is discussed. It is found that the length of the TPB is the only reasonable property of a 2D electrode that follows a 1/r particle radius relationship. Subsequently, sphere packing cuts are used to derive a statistical electrode surface that is the basis for the earlier proposed simulations of different electrochemical mechanisms. It is shown that two of the proposed mechanisms (conductivity and a two-step-two-material kinetic mechanism) can explain the current increase at $\mathrm{Pb-MnO_2}$ anodes compared to standard lead anodes. The results show that although $\mathrm{MnO_2}$ has low conductivity, when combined with Pb as the metal matrix, the behaviour of the composite is not purely ohmic but is also affected by activation overpotentials, increasing the current density close to the TPB. Current density is inversely proportional to the radius of the catalyst particles, matching with earlier experimental results. A hypothetical two-step-two-material mechanism with intermediate $\mathrm{H_2O_2}$ that reacts on both the Pb matrix and $\mathrm{MnO_2}$ catalyst is studied. It was found that assuming quasi-reversible generation of $\mathrm{H_2O_2}$ followed by its chemical decomposition on $\mathrm{MnO_2}$, results are obtained that agree with the experiments. It is further emphasised that both the Pb matrix and $\mathrm{MnO_2}$ catalyst are necessary and their optimum ratio depends on the used current density. Yet, additional experimental evidence is needed to support the postulated mechanism.


title: Achieving Large, Tunable Strain in Monolayer Transition-Metal Dichalcogenides

abstract: We describe a facile technique based on polymer encapsulation to apply several percent controllable strains to monolayer and few-layer Transition Metal Dichalcogenides (TMDs). We use this technique to study the lattice response to strain via polarized Raman spectroscopy in monolayer WSe2 and WS2. The application of strain causes mode-dependent redshifts, with larger shift rates observed for in-plane modes. We observe a splitting of the degeneracy of the in-plane E' modes in both materials and measure the Gruneisen parameters. At large strain, we observe that the reduction of crystal symmetry can lead to a change in the polarization response of the A' mode in WS2. While both WSe2 and WS2 exhibit similar qualitative changes in the phonon structure with strain, we observe much larger changes in mode positions and intensities with strain in WS2. These differences can be explained simply by the degree of iconicity of the metal-chalcogen bond.


title: Modelling high resolution ALMA observations of strongly lensed highly star forming galaxies detected by Herschel

abstract: We have modelled high resolution ALMA imaging of six strong gravitationally lensed galaxies detected by the Herschel Space Observatory. Our modelling recovers mass properties of the lensing galaxies and, by determining magnification factors, intrinsic properties of the lensed sub-millimetre sources. We find that the lensed galaxies all have high ratios of star formation rate to dust mass, consistent with or higher than the mean ratio for high redshift sub-millimetre galaxies and low redshift ultra-luminous infra-red galaxies. Source reconstruction reveals that most galaxies exhibit disturbed morphologies. Both the cleaned image plane data and the directly observed interferometric visibilities have been modelled, enabling comparison of both approaches. In the majority of cases, the recovered lens models are consistent between methods, all six having mass density profiles that are close to isothermal. However, one system with poor signal to noise shows mildly significant differences.