title: Robust topology optimisation of microstructural details without length scale separation - using a spectral coarse basis preconditioner

abstract: This paper applies topology optimisation to the design of structures with periodic microstructural details without length scale separation, i.e. considering the complete macroscopic structure and its response, while resolving all microstructural details, as compared to the often used homogenisation approach. The approach takes boundary conditions into account and ensures connected and macroscopically optimised microstructures regardless of the difference in micro- and macroscopic length scales. This results in microstructures tailored for specific applications rather than specific properties. Dealing with the complete macroscopic structure and its response is computationally challenging as very fine discretisations are needed in order to resolve all microstructural details. Therefore, this article shows the benefits of applying a contrast-independent spectral preconditioner based on the multiscale finite element method (MsFEM) to large structures with fully-resolved microstructural details. The density-based topology optimisation approach combined with a Heaviside projection filter and a stochastic robust formulation is used on various problems, with both periodic and layered microstructures. The presented approach is shown to allow for the topology optimisation of very large problems in \textsc{Matlab}, specifically a problem with 26 million displacement degrees of freedom in 26 hours using a single computational thread.


title: Sign of Green's function of Paneitz operators and the Q curvature

abstract: In a conformal class of metrics with positive Yamabe invariant, we derive a necessary and sufficient condition for the existence of metrics with positive Q curvature. The condition is conformally invariant. We also prove some inequalities between the Green's functions of the conformal Laplacian operator and the Paneitz operator.


title: Fatigue Estimation Methods Comparison for Wind Turbine Control

abstract: Fatigue is a critical factor in structures as wind turbines exposed to harsh operating conditions, both in the design stage and control during their operation. In the present paper the most recognized approaches to estimate the damage caused by fatigue are discussed and compared, with special focus on their applicability for wind turbine control. The aim of this paper is to serve as a guide among the vast literature on fatigue and shed some light on the underlying relationships between these methods.


title: Q curvature on a class of manifolds with dimension at least 5

abstract: For a smooth compact Riemannian manifold with positive Yamabe invariant, positive Q curvature and dimension at least 5, we prove the existence of a conformal metric with constant Q curvature. Our approach is based on the study of extremal problem for a new functional involving the Paneitz operator.


title: Confining Strings in Supersymmetric Theories with Higgs Branches

abstract: We study flux tubes (strings) on the Higgs branches in supersymmetric gauge theories. In generic vacua on the Higgs branches strings were shown to develop long-range "tails" associated with massless fields, a characteristic feature of the Higgs branch (the only exception is the vacuum at the base of the Higgs branch). A natural infrared regularization for the above tails is provided by a finite string length L. We perform a numerical study of these strings in generic vacua. We focus on the simplest example of strings in N=1 supersymmetric QED with the Fayet-Iliopoulos term. In particular, we examine the accuracy of a logarithmic approximation (proposed earlier by Evlampiev and Yung) for the tension of such string solutions. In the Evlampiev-Yung formula the dependence of tension on the string length is logarithmic and the dependence on the geodesic length from the base of the Higgs branch is quadratic. We observe a remarkable agreement of our numerical results for the string tension with the Evlampiev-Yung analytic expression.


title: Dark Bogolon-Excitons in a Linear Atomic Super-Lattice

abstract: Dark and bright excitons are shown to appear naturally in a linear atomic super-lattice with two atoms per unit cell. In bringing the super-lattice into a strong coupling regime with a one-dimensional nanophotonic waveguide, bright excitons and photons are coherently mixed to form polaritons. Treating excitons as bosons implies a mechanism that forbids two excitations from being at the same atomic state, which is included here through a bosonization procedure with kinematic interactions. Interestingly these interactions couple dark and bright excitons, and which we exploit as a new tool for exciting dark states in a controllable way. We suggest a pump-probe experiment where two polaritons scatter into two dark excitons that found to be correlated and are represented as dark bogolon-excitons. The results can be adapted for any super-lattice of active materials, e.g., of organic molecules.


title: Analog Signal Processing Solution for Image Alignment

abstract: Imaging and Image sensors is a field that is continuously evolving. There are new products coming into the market every day. Some of these have very severe Size, Weight and Power constraints whereas other devices have to handle very high computational loads. Some require both these conditions to be met simultaneously. Current imaging architectures and digital image processing solutions will not be able to meet these ever increasing demands. There is a need to develop novel imaging architectures and image processing solutions to address these requirements. In this work we propose analog signal processing as a solution to this problem. The analog processor is not suggested as a replacement to a digital processor but it will be used as an augmentation device which works in parallel with the digital processor, making the system faster and more efficient. In order to show the merits of analog processing the highly computational Normalized Cross Correlation algorithm is implemented. We propose two novel modifications to the algorithm and a new imaging architecture which, significantly reduces the computation time.


title: 1974: the discovery of the first binary pulsar

abstract: The 1974 discovery, by Russell A. Hulse and Joseph H. Taylor, of the first binary pulsar PSR B1913+16, opened up new possibilities for the study of relativistic gravity. PSR B1913+16, as well as several other binary pulsars, provided {\it direct} observational proofs that gravity propagates at the velocity of light and has a quadrupolar structure. Binary pulsars also provided accurate tests of the strong-field regime of relativistic gravity. General Relativity has passed all the binary pulsar tests with flying colors. The discovery of binary pulsars had also very important consequences for astrophysics: accurate measurement of neutron star masses, improved understanding of the possible evolution scenarios for the co-evolution of binary stars, proof of the existence of binary neutron stars emitting gravitational waves for hundreds of millions of years, before coalescing in catastrophic events radiating intense gravitational-wave signals, and probably leading also to important emissions of electromagnetic radiation and neutrinos. This article reviews the history of the discovery of the first binary pulsar, and describes both its immediate impact, and its longer-term effect on theoretical and experimental studies of relativistic gravity.


title: The muonic hydrogen Lamb shift and the proton radius

abstract: We obtain a model independent expression for the muonic hydrogen Lamb shift up to $\mathcal O (m_\mu \alpha^6, m_\mu \alpha^5 \frac{m_\mu^2}{m_\rho^2})$. The hadronic effects are controlled by the chiral theory, which allows for their model independent determination. We give their complete expression including the pion and Delta particles. Out of this analysis and the experimental measurement of the muonic hydrogen Lamb shift we determine the electromagnetic proton radius: $r_p$=0.8412(15) fm. This number is at 6.8$\sigma$ variance with respect to the CODATA value. The parametric control of the uncertainties allows us to obtain a model independent determination of the error, which is dominated by hadronic effects.


title: Electric Field Induced Topological Phase Transition in Two-Dimensional Few-layer Black Phosphorus

abstract: Phosphorene is a novel two-dimensional material that can be isolated through mechanical exfoliation from layered black phosphorus, but unlike graphene and silicene, monolayer phosphorene has a large band gap. It was thus unsuspected to exhibit band inversion and the ensuing topological insulator behavior. It has recently attracted interest because of its proposed application as field effect transistors. Using first-principles calculations with applied perpendicular electric field F we predict a continuous transition from the normal insulator to a topological insulator and eventually to a metal as a function of F. The continuous tuning of topological behavior with electric field would lead to spin-separated, gapless edge states, i.e., quantum spins Hall effect. This finding opens the possibility of converting normal insulating materials into topological ones via electric field, and making a multi-functional field effect topological transistor that could manipulate simultaneously both spins and charge carrier.


title: Cut and conjugate points of the exponential map, with applications

abstract: The goal of this thesis is to study the singularities of the exponential map of Riemannian and Finsler manifolds (a concept related to caustics and catastrophes), and the object known as the cut locus (aka ridge, medial axis or skeleton), to improve existing results about its structure, to look at it in new ways, and to derive applications to the Ambrose conjecture and the Hamilton-Jacobi equations.


title: SOPHIE velocimetry of Kepler transit candidates XIV. A joint photometric, spectroscopic, and dynamical analysis of the Kepler-117 system

abstract: As part of our follow-up campaign of Kepler planets, we observed Kepler-117 with the SOPHIE spectrograph at the Observatoire de Haute-Provence. This F8-type star hosts two transiting planets in non-resonant orbits. The planets, Kepler-117 b and c, have orbital periods $\simeq 18.8$ and $\simeq 50.8$ days, and show transit-timing variations (TTVs) of several minutes. We performed a combined Markov chain Monte Carlo (MCMC) fit on transits, radial velocities, and stellar parameters to constrain the characteristics of the system. We included the fit of the TTVs in the MCMC by modeling them with dynamical simulations. In this way, consistent posterior distributions were drawn for the system parameters. According to our analysis, planets b and c have notably different masses ($0.094 \pm 0.033$ and $1.84 \pm 0.18$ M$_{\rm J}$) and low orbital eccentricities ($0.0493 \pm 0.0062$ and $0.0323 \pm 0.0033$). The uncertainties on the derived parameters are strongly reduced if the fit of the TTVs is included in the combined MCMC. The TTVs allow measuring the mass of planet b, although its radial velocity amplitude is poorly constrained. Finally, we checked that the best solution is dynamically stable.


title: Overlapping Community Discovery Methods: A Survey

abstract: The detection of overlapping communities is a challenging problem which is gaining increasing interest in recent years because of the natural attitude of individuals, observed in real-world networks, to participate in multiple groups at the same time. This review gives a description of the main proposals in the field. Besides the methods designed for static networks, some new approaches that deal with the detection of overlapping communities in networks that change over time, are described. Methods are classified with respect to the underlying principles guiding them to obtain a network division in groups sharing part of their nodes. For each of them we also report, when available, computational complexity and web site address from which it is possible to download the software implementing the method.


title: Modeling and simulations of quantum phase slips in ultrathin superconducting wires

abstract: We study quantum phase slips (QPS) in ultrathin superconducting wires. Starting from an effective one-dimensional microscopic model, which includes electromagnetic fluctuations, we map the problem to a (1+1)-dimensional gas of interacting instantons. We introduce a method to calculate the tunneling amplitude of quantum phase slips directly from Monte Carlo simulations. This allows us to go beyond the dilute instanton gas approximation and study the problem without any limitations of the density of QPS. We find that the tunneling amplitude shows a characteristic scaling behavior near the superconductor-insulator transition. We also calculate the voltage-charge relation of the insulating state, which is the dual of the Josephson current-phase relation in ordinary superconducting weak links. This evolves from a sinusoidal form in the regime of dilute QPS to more exotic shapes for higher QPS densities, where interactions are important.


title: Entanglement Generation and Dynamics for a Bose-Hubbard model in a Double-Well Potential

abstract: The study of entanglement between bosonic systems is of primary importance for establishing feasible resources needed for implementing quantum information protocols, both in their interacting atomic or photonic realizations. Atomic systems are particularly efficient in the production of large amounts of entanglement, providing higher information density than conventional qubit entangled states. Such increased quantum resources pave the way to novel fundamental tests of nature and efficient applications in quantum information, metrology and sensing. We consider a basic setup made up of two parties A and B, each one populated by a single level bosonic variable. The bosons are interacting and can hop between A and B, thus describing a two-site Bose-Hubbard Hamiltonian. We consider the generation of quantum states in several situations that cover the majority of physical realizations: ground state, finite temperature, unitary dynamics, dissipation through dephasing and loss of particles. The system is analyzed through truncated exact diagonalization, as a function of the microscopic parameters. The non separability of the obtained quantum states is estimated by means of the negativity, which has recently been proven to be a suitable measure of entanglement. Finally, we calculate lower bounds of the entanglement of formation, an indicator that quantifies the minimal amount of entanglement resources needed to build up such states.


title: Reliable approximation of separatrix manifolds in competition models with safety niches

abstract: In dynamical systems saddle points partition the domain into basins of attractions of the remaining locally stable equilibria. This situation is rather common especially in population dynamics models, like prey-predator or competition systems. Focusing on squirrels population models with niche, in this paper we design algorithms for the detection and the refinement of points lying on the separatrix manifold partitioning the phase space. We consider both the two populations and the three populations cases. To reconstruct the separatrix curve and surface, we apply the Partition of Unity method, which makes use of Wendland's functions as local approximants.


title: Long-time behavior, invariant measures and regularizing effects for stochastic scalar conservation laws

abstract: We study the long-time behavior and the regularity of pathwise entropy solutions to stochastic scalar conservation laws with random in time spatially homogeneous fluxes and periodic initial data. We prove that the solutions converge to their spatial average, which is the unique invariant measure, and provide a rate of convergence, the latter being new even in the deterministic case for dimensions higher than two. The main tool is a new regularization result in the spirit of averaging lemmata for scalar conservation laws, which, in particular, implies a regularization by noise-type result for pathwise quasi-solutions.


title: Entanglement entropy in dynamic quantum-coherent conductors

abstract: We investigate the entanglement and the R\'enyi entropies of two electronic leads connected by a quantum point contact. For non-interacting electrons, the entropies can be related to the cumulants of the full counting statistics of transferred charge which in principle are measurable. We consider the entanglement entropy generated by operating the quantum point contact as a quantum switch which is opened and closed in a periodic manner. Using a numerically exact approach we analyze the conditions under which a logarithmic growth of the entanglement entropy predicted by conformal field theory should be observable in an electronic conductor. In addition, we consider clean single-particle excitations on top of the Fermi sea (levitons) generated by applying designed pulses to the leads. We identify a Hong-Ou-Mandel-like suppression of the entanglement entropy by interfering two levitons on a quantum point contact tuned to half transmission.


title: The replica method in liquid theory: from the basics to explicit computations

abstract: In these notes we introduce briefly the fundamentals of the replica method in the context of liquid theory and the structural glass problem. In particular, we explain and show its usefulness as a computation framework in the context of the Random First Order Transition (RFOT) theory of the glass transition, whose defining points the reader is assumed to know. We shall give the intuitive idea of how and why the replica method is suitable for the description of the glass transition (the dynamical glass transition in particular) in real liquids, and then show how it can be used to make explicit computations and predictions that can be compared to experiments and numerical simulations.


title: Signatures of anomalous Higgs couplings in angular asymmetries of H --> Zl+l- and e+e- --> HZ

abstract: Parametrizing beyond Standard Model physics by the SU(3) x SU(2)_L x U(1)_Y dimension-six effective Lagrangian, we study the impact of anomalous Higgs couplings in angular asymmetries of the crossing-symmetric processes H --> Zl+l- and e+e- --> HZ. In the light of present bounds on d=6 couplings, we show that some asymmetries can reveal BSM effects that would otherwise be hidden in other observables. The d=6 HZgamma couplings as well as (to a lesser extent) HZll contact interactions can generate asymmetries at the several percent level, albeit having less significant effects on the di-lepton invariant mass spectrum of the decay H --> Zl+l-. The higher di-lepton invariant mass probed in e+e- --> HZ can lead to complementary anomalous coupling searches at e+e- colliders.


title: Discrete solitons and vortices on two-dimensional lattices of $\mathcal{PT}$-symmetric couplers

abstract: We introduce a 2D network built of $\mathcal{PT}$-symmetric dimers with on-site cubic nonlinearity, the gain and loss elements of the dimers being linked by parallel square-shaped lattices. The system may be realized as a set of $\mathcal{PT}$-symmetric dual-core waveguides embedded into a photonic crystal. The system supports $\mathcal{PT}$-symmetric and antisymmetric fundamental solitons (FSs) and on-site-centered solitary vortices (OnVs). Stability of these discrete solitons is the central topic of the consideration. Their stability regions in the underlying parameter space are identified through the computation of stability eigenvalues, and verified by direct simulations. Symmetric FSs represent the system's ground state, being stable at lowest values of the power, while anti-symmetric FSs and OnVs are stable at higher powers. Symmetric OnVs, which are also stable at lower powers, are remarkably robust modes: on the contrary to other PT-symmetric states, unstable OnVs do not blow up, but spontaneously rebuild themselves into stable FSs.


title: Discrete solitons and scattering of lattice waves in guiding arrays with a nonlinear $\mathcal{PT}$-symmetric defect

abstract: Discrete fundamental and dipole solitons are constructed, in an exact analytical form, in an array of linear waveguides with an embedded $\mathcal{PT}$-symmetric dimer, which is composed of two nonlinear waveguides carrying equal gain and loss. Fundamental solitons in tightly knit lattices, as well as all dipole modes, exist above a finite threshold value of the total power. However, the threshold vanishes for fundamental solitons in loosely knit lattices. The stability of the discrete solitons is investigated analytically by means of the Vakhitov-Kolokolov (VK) criterion, and, in the full form, via the computation of eigenvalues for perturbation modes. Fundamental and dipole solitons tend to be stable at smaller and larger values of the total power (norm), respectively. The increase of the strength of the coupling between the two defect-forming sites leads to strong expansion of the stability areas. The scattering problem for linear lattice waves impinging upon the defect is considered too.


title: Detecting nonlocal Cooper pair entanglement by optical Bell inequality violation

abstract: Based on the Bardeen Cooper Schrieffer (BCS) theory of superconductivity, the coherent splitting of Cooper pairs from a superconductor to two spatially separated quantum dots has been predicted to generate nonlocal pairs of entangled electrons. In order to test this hypothesis, we propose a scheme to transfer the spin state of a split Cooper pair onto the polarization state of a pair of optical photons. We show that the produced photon pairs can be used to violate a Bell inequality, unambiguously demonstrating the entanglement of the split Cooper pairs.


title: Level density of a Bose gas: beyond the saddle point approximation

abstract: The present article is concerned with the use of approximations in the calculation of the many-body density of states (MBDS) of a system with total energy E, composed by N bosons. In the mean-field framework, an integral expression for MBDS, which is proper to be performed by asymptotic expansions, can be derived. However, the standard second order steepest descent method cannot be applied to this integral when the ground-state is sufficiently populated. Alternatively, we derive a uniform formula for MBDS, which is potentially able to deal with this regime. In the case of the one-dimensional harmonic oscillator, using results found in the number theory literature, we show that the uniform formula improves the standard expression achieved by means of the second order method.


title: Incorporating Views on Market Dynamics in Options Hedging

abstract: We examine the possibility of incorporating information or views of market movements during the holding period of a portfolio, in the hedging of European options with respect to the underlying. Given a fixed holding period interval, we explore whether it is possible to adjust the number of shares needed to effectively hedge our position to account for views on market dynamics from present until the end of our interval, to account for the time-dependence of the options' sensitivity to the underlying. We derive an analytical expression for the number of shares needed by adjusting the standard Black-Scholes-Merton $\Delta$ quantity, in the case of an arbitrary process for implied volatility, and we present numerical results.


title: Designing for the Don't Cares: A story about a sociotechnical system

abstract: This article discusses the difficulties that arose when attempting to specify and design a large scale digital learning environment for Scottish schools. This had a potential user base of about 1 million users and was intended to replace an existing, under-used system. We found that the potential system users were not interested in engaging with the project and that there were immense problems with system governance. The only technique that we found to be useful were user stories, presenting scenarios of how the system might be used by students and their teachers. The designed architecture was based around a layered set of replaceable services.


title: Routing Diverse Crowds in Emergency with Dynamic Grouping

abstract: Evacuee routing algorithms in emergency typically adopt one single criterion to compute desired paths and ignore the specific requirements of users caused by different physical strength, mobility and level of resistance to hazard. In this paper, we present a quality of service (QoS) driven multi-path routing algorithm to provide diverse paths for different categories of evacuees. This algorithm borrows the concept of Cognitive Packet Network (CPN), which is a flexible protocol that can rapidly solve optimal solution for any user-defined goal function. Spatial information regarding the location and spread of hazards is taken into consideration to avoid that evacuees be directed towards hazardous zones. Furthermore, since previous emergency navigation algorithms are normally insensitive to sudden changes in the hazard environment such as abrupt congestion or injury of civilians, evacuees are dynamically assigned to several groups to adapt their course of action with regard to their on-going physical condition and environments. Simulation results indicate that the proposed algorithm which is sensitive to the needs of evacuees produces better results than the use of a single metric. Simulations also show that the use of dynamic grouping to adjust the evacuees' category and routing algorithms with regard for their on-going health conditions and mobility, can achieve higher survival rates.


title: Spin-orbit coupling and chaotic rotation for eccentric coorbital bodies

abstract: The presence of a co-orbital companion induces the splitting of the well known Keplerian spin-orbit resonances. It leads to chaotic rotation when those resonances overlap.


title: The $L^1$ gradient flow of a generalized scale invariant Willmore energy for radially non increasing functions

abstract: We use the minimizing movement theory to study the gradient flow associated with a non-regular relaxation of a geometric functional derived from the Willmore energy. Thanks to the coarea formula, one can define a Willmore energy on regular functions of L 1 (R d). This functional is extended to every L 1 function by taking its lower semi-continuous envelope. We study the flow generated by this relaxed energy for radially non-increasing functions, i.e. functions with balls as level sets. In the first part of the paper, we prove a coarea formula for the relaxed energy of such functions. Then we show that the flow consists on an erosion of the initial data. The erosion speed is given by a first order ordinary equation.


title: Dynamical screening in La2CuO4

abstract: We show that the dynamical screening of the Coulomb interaction among Cu-d electrons in high-Tc cuprates is very strong and that a proper treatment of this effect is essential for a consistent description of the electronic structure. In particular, we find that ab-initio calculations for undoped La2CuO4 yield an insulator only if the frequency dependence of the Coulomb interaction is taken into account. We also identify a collective excitation in the screened interaction at 9 eV which is rather localized on the copper site, and which is responsible for a satellite structure at energy -13 eV, located below the p bands.


title: Velocity-aided Attitude Estimation for Accelerated Rigid Bodies

abstract: Two nonlinear observers for velocity-aided attitude estimation, relying on gyrometers, accelerometers, magnetometers, and velocity measured in the body-fixed frame, are proposed. As opposed to state-of-the-art body-fixed velocity-aided attitude observers endowed with local properties, both observers are (almost) globally asymptotically stable, with very simple and flexible tuning. Moreover, the roll and pitch estimates are globally decoupled from magnetometer measurements.


title: Optimal mixture weights in multiple importance sampling

abstract: In multiple importance sampling we combine samples from a finite list of proposal distributions. When those proposal distributions are used to create control variates, it is possible (Owen and Zhou, 2000) to bound the ratio of the resulting variance to that of the unknown best proposal distribution in our list. The minimax regret arises by taking a uniform mixture of proposals, but that is conservative when there are many components. In this paper we optimize the mixture component sampling rates to gain further efficiency. We show that the sampling variance of mixture importance sampling with control variates is jointly convex in the mixture probabilities and control variate regression coefficients. We also give a sequential importance sampling algorithm to estimate the optimal mixture from the sample data.


title: Zeroth Order Phase Transition in a Holographic Superconductor with Single Impurity

abstract: We investigate the single normal impurity effect in a superconductor by the holographic method. When the size of impurity is much smaller than the host superconductor, we can reproduce the Anderson theorem, which states that a conventional s-wave superconductor is robust to a normal (non-magnetic) impurity with small impurity strength. However, by increasing the size of the impurity in a fixed-size host superconductor, we find a decreasing critical temperature $T_c$ of the host superconductor, which agrees with the results in condensed matter literatures. More importantly, the phase transition at the critical impurity strength (or the critical temperature) is of zeroth order.


title: Self-organization of microcircuits in networks of neurons with plastic synapses

abstract: The synaptic connectivity of cortical networks features an overrepresentation of certain wiring motifs compared to simple random-network models. This structure is shaped, in part, by synaptic plasticity that promotes or suppresses connections between neurons depending on their spiking activity. Frequently, theoretical studies focus on how feedforward inputs drive plasticity to create this network structure. We study the complementary scenario of self-organized structure in a recurrent network, with spike timing-dependent plasticity driven by spontaneous dynamics. We develop a self-consistent theory that describes the evolution of network structure by combining fast spiking covariance with a fast-slow theory for synaptic weight dynamics. Through a finite-size expansion of network dynamics, we obtain a low-dimensional set of nonlinear differential equations for the evolution of two-synapse connectivity motifs. With this theory in hand, we explore how the form of the plasticity rule drives the evolution of microcircuits in cortical networks. When potentiation and depression are in approximate balance, synaptic dynamics depend on the frequency of weighted divergent, convergent, and chain motifs. For additive, Hebbian STDP, these motif interactions create instabilities in synaptic dynamics that either promote or suppress the initial network structure. Our work provides a consistent theoretical framework for studying how spiking activity in recurrent networks interacts with synaptic plasticity to determine network structure.


title: Stability and Tunneling Dynamics of a Dark-Bright Soliton Pair in a Harmonic Trap

abstract: We consider a binary repulsive Bose-Einstein condensate in a harmonic trap in one spatial dimension and investigate particular solutions consisting of two dark-bright (DB) solitons. There are two different stationary solutions characterized by the phase difference in the bright component, in-phase and out-of-phase states. We show that above a critical particle number in the bright component, a symmetry breaking bifurcation of the pitchfork type occurs that leads to a new asymmetric solution whereas the parental branch, i.e., the out-of-phase state becomes unstable. These three different states support different small amplitude oscillations, characterized by an almost stationary density of the dark component and a tunneling of the bright component between the two dark solitons. Within a suitable effective double-well picture, these can be understood as the characteristic features of a Bosonic Josephson Junction (BJJ), and we show within a two-mode approach that all characteristic features of the BJJ phase space are recovered. For larger deviations from the stationary states, the simplifying double-well description breaks down due to the feedback of the bright component onto the dark one, causing the solitons to move. In this regime we observe intricate anharmonic and aperiodic dynamics, exhibiting remnants of the BJJ phase space.


title: Edge-state enhanced transport in a 2-dimensional quantum walk

abstract: Quantum walks on translation invariant regular graphs spread quadratically faster than their classical counterparts. The same coherence that gives them this quantum speedup inhibits, or even stops their spread in the presence of disorder. We ask how to create an efficient transport channel from a fixed source site (A) to fixed target site (B) in a disordered 2-dimensional discrete-time quantum walk by cutting some of the links. We show that the somewhat counterintuitive strategy of cutting links along a single line connecting A to B creates such a channel. The efficient transport along the cut is due to topologically protected chiral edge states, which exist even though the bulk Chern number in this system vanishes. We give a realization of the walk as a periodically driven lattice Hamiltonian, and identify the bulk topological invariant responsible for the edge states as the quasienergy winding of this Hamiltonian.


title: Hamilton-Jacobi theory in Cauchy data space

abstract: Recently, M. de Le\'on el al. ([9]) have developed a geometric Hamilton-Jacobi theory for Classical Field Theories in the setting of multisymplectic geometry. Our purpose in the current paper is to establish the corresponding Hamilton-Jacobi theory in the Cauchy data space, and relate both approaches.


title: The Herschel Virgo Cluster Survey. XVIII. Star-forming dwarf galaxies in a cluster environment

abstract: To assess the effects of the cluster environment on the different components of the interstellar medium, we analyse the FIR-submm properties of a sample of star-forming dwarf (SFD) galaxies detected by the Herschel Virgo Cluster Survey (HeViCS). We determine dust masses and dust temperatures by fitting a modified black body (MBB) function to the spectral energy distributions (SEDs). Stellar and gas masses, star formation rates (SFRs), and metallicities are obtained from the analysis of a set of ancillary data. Dust is detected in 49 out of 140 optically identified dwarfs covered by the HeViCS field; considering only dwarfs brighter than $m_B$ = 18 mag, this gives a detection rate of 43%. After evaluating different emissivity indices, we find that the FIR-submm SEDs are best-fit by $\beta$=1.5, with a median dust temperature $T_d$ = 22.4 K. Assuming $\beta$=1.5, 67% of the 23 galaxies detected in all five Herschel bands show emission at 500 $\mu$m in excess of the MBB model. The excess is inversely correlated with SFR and stellar masses. To study the variations in the global properties of our sample due to environmental effects, we compare the Virgo SFDs to other Herschel surveys, such as KINGFISH, the Dwarf Galaxy Survey (DGS), and the HeViCS bright galaxy catalogue (BGC). We explore the relations between stellar mass and HI fraction, specific SFR, dust fraction, gas-to-dust ratio over a wide range of stellar masses. Highly HI-deficient Virgo dwarf galaxies are mostly characterised by quenched star formation activity and lower dust fractions giving hints for dust stripping in cluster dwarfs. However, we find that the fraction of dust removed has to be less than that of the HI component. Since the Virgo SFDs are likely to be crossing the cluster for the first time, a longer timescale might be necessary to strip the more centrally concentrated dust distribution.


title: Privacy-preserving Loyalty Programs

abstract: Loyalty programs are promoted by vendors to incentivize loyalty in buyers. Although such programs have become widespread, they have been criticized by business experts and consumer associations: loyalty results in profiling and hence in loss of privacy of consumers. We propose a protocol for privacy-preserving loyalty programs that allows vendors and consumers to enjoy the benefits of loyalty (returning customers and discounts, respectively), while allowing consumers to stay anonymous and empowering them to decide how much of their profile they reveal to the vendor. The vendor must offer additional reward if he wants to learn more details on the consumer's profile. Our protocol is based on partially blind signatures and generalization techniques, and provides anonymity to consumers and their purchases, while still allowing negotiated consumer profiling.


title: Galois Transformers and Modular Abstract Interpreters

abstract: The design and implementation of static analyzers has become increasingly systematic. Yet for a given language or analysis feature, it often requires tedious and error prone work to implement an analyzer and prove it sound. In short, static analysis features and their proofs of soundness do not compose well, causing a dearth of reuse in both implementation and metatheory. We solve the problem of systematically constructing static analyzers by introducing Galois transformers: monad transformers that transport Galois connection properties. In concert with a monadic interpreter, we define a library of monad transformers that implement building blocks for classic analysis parameters like context, path, and heap (in)sensitivity. Moreover, these can be composed together independent of the language being analyzed. Significantly, a Galois transformer can be proved sound once and for all, making it a reusable analysis component. As new analysis features and abstractions are developed and mixed in, soundness proofs need not be reconstructed, as the composition of a monad transformer stack is sound by virtue of its constituents. Galois transformers provide a viable foundation for reusable and composable metatheory for program analysis. Finally, these Galois transformers shift the level of abstraction in analysis design and implementation to a level where non-specialists have the ability to synthesize sound analyzers over a number of parameters.


title: Search for heavy neutrinos in $K^+\to\mu^+\nu_H$ decays

abstract: Evidence of a heavy neutrino, $\nu_H$, in the $K^+\to\mu^+\nu_H$ decays was sought using the E949 experimental data with an exposure of $1.70\times 10^{12}$ stopped kaons. With the major background from the radiative $K^+\to\mu^+\nu_\mu\gamma$ decay understood and suppressed, upper limits (90% C.L.) on the neutrino mixing matrix element between muon and heavy neutrino, $|U_{\mu H}|^2$, were set at the level of $10^{-7}$ to $10^{-9}$ for the heavy neutrino mass region 175 to 300 MeV/$c^2$.


title: Fast Algorithm for Simulating Lipid Vesicle Deformation I: Spherical Harmonic Approximation

abstract: Lipid vesicles appear ubiquitously in biological systems. Understanding how the mechanical and intermolecular interations deform vesicle membrane is a fundamental question in biophysics. In this article we developed a fast algorithm to compute the surface configurations of lipid vesicles by introducing the surface harmonic functions to approximate the surfaces. This parameterization of the surfaces allows an analytical computation of the membrane curvature energy and its gradient for the efficient minimization of the curvature energy using a nonlinear conjugate gradient method. Our approach drastically reduces the degrees of freedom for approximating the membrane surfaces compared to the previously developed finite element and finite difference methods. Vesicle deformations with a reduced volume larger than 0.65 can be well approximated by using as small as 49 surface harmonic functions. The method thus has a great potential to reduce the computational expense of tracking multiple vesicles which deform for their interaction with external fields.


title: Mechanism of Outflows in Accretion System: Advective Cooling Cannot Balance Viscous Heating?

abstract: Based on no-outflow assumption, we investigate steady state, axisymmetric, optically thin accretion flows in spherical coordinates. By comparing the vertically integrated advective cooling rate with the viscous heating rate, we find that the former is generally less than 30% of the latter, which indicates that the advective cooling itself cannot balance the viscous heating. As a consequence, for radiatively inefficient flows with low accretion rates such as $\dot M \la 10^{-3} \dot M_{Edd}$, where $\dot M_{Edd}$ is the Eddington accretion rate, the viscous heating rate will be larger than the sum of the advective cooling rate and the radiative cooling one. Thus, no thermal equilibrium can be established under the no-outflow assumption. We therefore argue that in such case outflows ought to occur and take away more than 70% of the thermal energy generated by viscous dissipation. Similarly, for optically thick flows with extremely large accretion rates such as $\dot M \ga 10 \dot M_{Edd}$, outflows should also occur due to the limited advection and the low efficiency of radiative cooling. Our results may help to understand the mechanism of outflows found in observations and numerical simulations.


title: Multiwavelength Achromatic Metasurfaces by Dispersive Phase Compensation

abstract: The replacement of bulk refractive optical elements with diffractive planar components enables the miniaturization of optical systems. However, diffractive optics suffers from large chromatic aberrations due to the dispersion of the phase accumulated by light during propagation. We show that this limitation can be overcome with an engineered wavelength-dependent phase shift imparted by a metasurface and demonstrate a design that deflects three wavelengths without dispersion. A planar lens without chromatic aberrations at three wavelengths is also presented. Our design is based on low-loss dielectric resonators which introduce a dense spectrum of optical modes to enable dispersive phase compensation. The suppression of chromatic aberrations in metasurface-based planar photonics will find applications in lightweight collimators for displays, and chromatically-corrected imaging systems.


title: Relatively Complete Counterexamples for Higher-Order Programs

abstract: In this paper, we study the problem of generating inputs to a higher-order program causing it to error. We first study the problem in the setting of PCF, a typed, core functional language and contribute the first relatively complete method for constructing counterexamples for PCF programs. The method is relatively complete in the sense of Hoare logic; completeness is reduced to the completeness of a first-order solver over the base types of PCF. In practice, this means an SMT solver can be used for the effective, automated generation of higher-order counterexamples for a large class of programs. We achieve this result by employing a novel form of symbolic execution for higher-order programs. The remarkable aspect of this symbolic execution is that even though symbolic higher-order inputs and values are considered, the path condition remains a first-order formula. Our handling of symbolic function application enables the reconstruction of higher-order counterexamples from this first-order formula. After establishing our main theoretical results, we sketch how to apply the approach to untyped, higher-order, stateful languages with first-class contracts and show how counterexample generation can be used to detect contract violations in this setting. To validate our approach, we implement a tool generating counterexamples for erroneous modules written in Racket.


title: Parallel Neutrino Triggers using GPUs for an underwater telescope

abstract: Graphics Processing Units are high performance co-processors originally intended to improve the use and the acceleration of computer graphics applications. Because of their performance, researchers have extended their use beyond the computer graphics scope. We have investigated the possibility of implementing online neutrino trigger algorithms in the KM3Net-It experiment using a CPU-GPU system. The results of a neutrino trigger simulation on a NEMO Phase II tower and a KM3-It 14 floors tower are reported.


title: Formal Semantic Annotations for Models Interoperability in a PLM environment

abstract: Nowadays, the need for system interoperability in or across enterprises has become more and more ubiquitous. Lots of research works have been carried out in the information exchange, transformation, discovery and reuse. One of the main challenges in these researches is to overcome the semantic heterogeneity between enterprise applications along the lifecycle of a product. As a possible solution to assist the semantic interoperability, semantic annotation has gained more and more attentions and is widely used in different domains. In this paper, based on the investigation of the context and the related works, we identify some existing drawbacks and propose a formal semantic annotation approach to support the semantics enrichment of models in a PLM environment.


title: Approaches for Synthesis Conjectures in an SMT Solver

abstract: This report describes several approaches for handling synthesis conjectures within an Satisfiability Modulo Theories (SMT) solver. We describe approaches that primarily focus on determining the unsatisfiability of the negated form of synthesis conjectures using new techniques for quantifier instantiation.


title: Finite-horizon optimal multiple switching with signed switching costs

abstract: This paper is concerned with optimal switching over multiple modes in continuous time and on a finite horizon. The performance index includes a running reward, terminal reward and switching costs that can belong to a large class of stochastic processes. Particularly, the switching costs are modelled by right-continuous with left-limits processes that are quasi-left-continuous and can take both positive and negative values. We provide sufficient conditions leading to a well known probabilistic representation of the value function for the switching problem in terms of interconnected Snell envelopes. We also prove the existence of an optimal strategy within a suitable class of admissible controls, defined iteratively in terms of the Snell envelope processes.


title: Causal Inference by Identification of Vector Autoregressive Processes with Hidden Components

abstract: A widely applied approach to causal inference from a non-experimental time series $X$, often referred to as "(linear) Granger causal analysis", is to regress present on past and interpret the regression matrix $\hat{B}$ causally. However, if there is an unmeasured time series $Z$ that influences $X$, then this approach can lead to wrong causal conclusions, i.e., distinct from those one would draw if one had additional information such as $Z$. In this paper we take a different approach: We assume that $X$ together with some hidden $Z$ forms a first order vector autoregressive (VAR) process with transition matrix $A$, and argue why it is more valid to interpret $A$ causally instead of $\hat{B}$. Then we examine under which conditions the most important parts of $A$ are identifiable or almost identifiable from only $X$. Essentially, sufficient conditions are (1) non-Gaussian, independent noise or (2) no influence from $X$ to $Z$. We present two estimation algorithms that are tailored towards conditions (1) and (2), respectively, and evaluate them on synthetic and real-world data. We discuss how to check the model using $X$.


title: The impact of global nuclear mass model uncertainties on $r$-process abundance predictions

abstract: Rapid neutron capture or `$r$-process' nucleosynthesis may be responsible for half the production of heavy elements above iron on the periodic table. Masses are one of the most important nuclear physics ingredients that go into calculations of $r$-process nucleosynthesis as they enter into the calculations of reaction rates, decay rates, branching ratios and Q-values. We explore the impact of uncertainties in three nuclear mass models on $r$-process abundances by performing global monte carlo simulations. We show that root-mean-square (rms) errors of current mass models are large so that current $r$-process predictions are insufficient in predicting features found in solar residuals and in $r$-process enhanced metal poor stars. We conclude that the reduction of global rms errors below $100$ keV will allow for more robust $r$-process predictions.


title: CMB probes on the correlated axion isocurvature perturbation

abstract: We explore the possible cosmological consequence of the gravitational coupling between the inflaton and axion-like fields. In view of the forthcoming cosmic microwave background (CMB) polarization and lensing data, we study the sensitivity of the CMB data on the cross-correlation between the curvature and axion isocurvature perturbations. Through a concrete example, we illustrate the explicit dependence of the scale dependent cross-correlation power spectrum on the axion parameters.


title: Jamming Criticality Revealed by Removing Localized Buckling Excitations

abstract: Recent theoretical advances offer an exact, first-principle theory of jamming criticality in infinite dimension as well as universal scaling relations between critical exponents in all dimensions. For packings of frictionless spheres near the jamming transition, these advances predict that nontrivial power-law exponents characterize the critical distribution of (i) small inter-particle gaps and (ii) weak contact forces, both of which are crucial for mechanical stability. The scaling of the inter-particle gaps is known to be constant in all spatial dimensions $d$ -- including the physically relevant $d=2$ and 3, but the value of the weak force exponent remains the object of debate and confusion. Here, we resolve this ambiguity by numerical simulations. We construct isostatic jammed packings with extremely high accuracy, and introduce a simple criterion to separate the contribution of particles that give rise to localized buckling excitations, i.e., bucklers, from the others. This analysis reveals the remarkable dimensional robustness of mean-field marginality and its associated criticality.


title: Generation of Cold Argon Plasma Jet at the End of Flexible Plastic Tube

abstract: This brief communication reports a new method for generation of cold atmospheric pressure plasma jet at the downstream end of a flexible plastic tube. The device consists of a small chamber where dielectric barrier discharge (DBD) is ignited in Argon. The discharge is driven by a conventional low frequency AC power supply. The exit of DBD reactor is connected to a commercial flexible plastic tube (up to 4 meters long) with a thin floating Cu wire inside. Under certain conditions an Ar plasma jet can be extracted from the downstream tube end and there is no discharge inside the plastic tube. The jet obtained by this method is cold enough to be put in direct contact with human skin without electric shock and can be used for medical treatment and decontamination.


title: Multi-curve HJM modelling for risk management

abstract: We present a HJM approach to the projection of multiple yield curves developed to capture the volatility content of historical term structures for risk management purposes. Since we observe the empirical data at daily frequency and only for a finite number of time-to-maturity buckets, we propose a modelling framework which is inherently discrete. In particular, we show how to approximate the HJM continuous time description of the multi-curve dynamics by a Vector Autoregressive process of order one. The resulting dynamics lends itself to a feasible estimation of the model volatility-correlation structure and market risk-premia. Then, resorting to the Principal Component Analysis we further simplify the dynamics reducing the number of covariance components. Applying the constant volatility version of our model on a sample of curves from the Euro area, we demonstrate its forecasting ability through an out-of-sample test.


title: Nuclear masses near $N=82$ that influence $r$-process abundances

abstract: Nuclear masses are one of the key ingredients of nuclear physics that go into astrophysical simulations of the $r$ process. Nuclear masses effect $r$-process abundances by entering into calculations of Q-values, neutron capture rates, photo-dissociation rates, beta-decay rates, branching ratios and the properties of fission. Most of the thousands of short-lived neutron-rich nuclei which are believed to participate in the $r$ process lack any experimental verification, thus the identification of the most influential nuclei is of paramount importance. We have conducted mass sensitivity studies near the $N=82$ closed shell in the context of a main $r$-process. Our studies take into account how an uncertainty in a single nuclear mass propagates to influence the relevant quantities of neighboring nuclei and finally to $r$-process abundances. We identify influential nuclei in various astrophysical conditions using the FRDM mass model. We show that our conclusions regarding these key nuclei are still retained when a superposition of astrophysical trajectories is considered.


title: Brane structure from a scalar field in general covariant Horava-Lifshitz gravity

abstract: In this paper we have considered the structure of the non-projectable Horava-Melby-Thompson (HMT) gravity to find braneworld scenarios. A relativistic scalar field is considered in the matter sector and we have shown how to reduce the equations of motion to first-order differential equations. In particular, we have studied thick brane solutions of both the dilatonic and Randall-Sundrum types.


title: Addressing Hadronic Uncertainties in Extractions of $\phi_s$

abstract: In light of recent LHC results for the extraction of the $B_s$ mixing phase $\phi_s$, we can already conclude that if New Physics (NP) is present in this observable, it is hiding pretty well. Thus, as our hunt continues, we must be weary not to confuse NP for penguin effects, or vice versa. In this talk the progress made towards addressing hadronic uncertainties in extractions of $\phi_s$ from $B_s\to J/\psi \phi$ is reviewed, and the nature of the scalar $f_0(980)$ state, which plays a dominant role in the extraction of $\phi_s$ from the $B_s\to J/\psi \pi^+\pi^-$ decay, is discussed.


title: Dynamic programming for discrete-time finite horizon optimal switching problems with negative switching costs

abstract: This paper studies a discrete-time optimal switching problem on a finite horizon. The underlying model has a running reward, terminal reward and signed (positive and negative) switching costs. Using the martingale approach to optimal stopping problems, we extend a well known explicit dynamic programming method for computing the value function and the optimal strategy to the case of signed switching costs.


title: Simulation of QCD with N_f=2+1 flavors of non-perturbatively improved Wilson fermions

abstract: We describe a new set of gauge configurations generated within the CLS effort. These ensembles have N_f=2+1 flavors of non-perturbatively improved Wilson fermions in the sea with the Luescher-Weisz action used for the gluons. Open boundary conditions in time are used to address the problem of topological freezing at small lattice spacings and twisted-mass reweighting for improved stability of the simulations. We give the bare parameters at which the ensembles have been generated and how these parameters have been chosen. Details of the algorithmic setup and its performance are presented as well as measurements of the pion and kaon masses alongside the scale parameter t_0.


title: Light Higgs channel of the resonant decay of magnon condensate in superfluid $^3$He-B

abstract: In superfluids the order parameter, which describes spontaneous symmetry breaking, is an analogue of the Higgs field in the Standard Model of particle physics. Oscillations of the field amplitude are massive Higgs bosons, while oscillations of the orientation are massless Nambu-Goldstone bosons. The 125~GeV Higgs boson, discovered at Large Hadron Collider, is light compared to electroweak energy scale, which led to a suggestion of the ``little Higgs'' extension of the Standard Model, in which the light Higgs appears as a NG mode acquiring mass due to violation of a hidden symmetry. Here we show that such light Higgs exists in superfluid $^3$He-B, where one of three Nambu-Goldstone spin-wave modes acquires small mass due to the spin-orbit interaction. Other modes become optical and acoustic magnons. We observe parametric decay of Bose-Einstein condensate of optical magnons to light Higgs modes and decay of optical to acoustic magnons. Formation of a light Higgs from a Nambu-Goldstone mode observed in $^3$He-B opens a possibility that such scenario can be realized in other systems, where violation of some hidden symmetry is possible, including the Standard Model.


title: Qualitative Robustness in Bayesian Inference

abstract: The practical implementation of Bayesian inference requires numerical approximation when closed-form expressions are not available. What types of accuracy (convergence) of the numerical approximations guarantee robustness and what types do not? In particular, is the recursive application of Bayes' rule robust when subsequent data or posteriors are approximated? When the prior is the push forward of a distribution by the map induced by the solution of a PDE, in which norm should that solution be approximated? Motivated by such questions, we investigate the sensitivity of the distribution of posterior distributions (i.e. posterior distribution-valued random variables, randomized through the data) with respect to perturbations of the prior and data generating distributions in the limit when the number of data points grows towards infinity.


title: Intermittency in the relative separations of tracers and of heavy particles in turbulent flows

abstract: Results from Direct Numerical Simulations of particle relative dispersion in three dimensional homogeneous and isotropic turbulence at Reynolds number $Re_\lambda \sim 300$ are presented. We study point-like passive tracers and heavy particles, at Stokes number St = 0, 0.6, 1 and 5. Particles are emitted from localised sources, in bunches of thousands, periodically in time, allowing to reach an unprecedented statistical accuracy, with a total number of events for two-point observables of the order of $10^{11}$. The right tail of the probability density function for tracers develops a clear deviation from Richardson's self-similar prediction, pointing to the intermittent nature of the dispersion process. In our numerical experiment, such deviations are manifest once the probability to measure an event becomes of the order of -or rarer than- one part over one million, hence the crucial importance of a large dataset. The role of finite-Reynolds effects and the related fluctuations when pair separations cross the boundary between viscous and inertial range scales are discussed. An asymptotic prediction based on the multifractal theory for inertial range intermittency and valid for large Reynolds numbers is found to agree with the data better than the Richardson theory. The agreement is improved when considering heavy particles, whose inertia filters out viscous scale fluctuations. By using the exit-time statistics we also show that events associated to pairs experiencing unusually slow inertial range separations have a non self-similar probability distribution function.


title: Plasmonic waveguides with hyperbolic multilayer cladding

abstract: Engineering plasmonic metamaterials with anisotropic optical dispersion enables us to tailor the properties of metamaterial-based waveguides. We investigate plasmonic waveguides with dielectric cores and multilayer metal-dielectric claddings with hyperbolic dispersion. Without using any homogenization, we calculate the resonant eigenmodes of the finite-width cladding layers, and find agreement with the resonant features in the dispersion of the cladded waveguides. We show that at the resonant widths, the propagating modes of the waveguides are coupled to the cladding eigenmodes and hence, are strongly absorbed. By avoiding the resonant widths in the design of the actual waveguides, the strong absorption can be eliminated.


title: Photons that travel in free space slower than the speed of light

abstract: That the speed of light in free space is constant is a cornerstone of modern physics. However, light beams have finite transverse size, which leads to a modification of their wavevectors resulting in a change to their phase and group velocities. We study the group velocity of single photons by measuring a change in their arrival time that results from changing the beam's transverse spatial structure. Using time-correlated photon pairs we show a reduction of the group velocity of photons in both a Bessel beam and photons in a focused Gaussian beam. In both cases, the delay is several microns over a propagation distance of the order of 1 m. Our work highlights that, even in free space, the invariance of the speed of light only applies to plane waves. Introducing spatial structure to an optical beam, even for a single photon, reduces the group velocity of the light by a readily measurable amount.


title: Superradiance on the Reissner-Nordstrom metric

abstract: In this article, we study the superradiance of charged scalar fields on the sub-extremal Reissner-Nordstrom metric, a mechanism by which such fields can extract energy from a static spherically symmetric charged black hole. A geometrical way of measuring the amount of energy extracted is proposed. Then we investigate the question numerically. The toy-model and the numerical methods used in our simulations are presented and the problem of long time measurement of the outgoing energy flux is discussed. We provide a numerical example of a field exhibiting a behaviour analogous to the Penrose process: an incoming wave packet which splits, as it approaches the black hole, into an incoming part with negative energy and an outgoing part with more energy than the initial incoming one. We also show another type of superradiant solution for which the energy extraction is more important. Hyperradiant behaviour is not observed, which is an indication that the Reissner-Nordstrom metric is linearly stable in the sub-extremal case.


title: Symplectic non-squeezing in Hilbert space and discrete Schr\"odinger equations

abstract: We prove a generalization of Gromov's symplectic non-squeezing theorem for the case of Hilbert spaces. Our approach is based on filling almost complex Hilbert spaces by complex discs partially extending Gromov's results on existence of $J$-complex curves. We apply our result to the flow of the discrete nonlinear Schr\"odinger equation.


title: Dynamic Modeling of Cascading Failure in Power Systems

abstract: The modeling of cascading failure in power systems is difficult because of the many different mechanisms involved; no single model captures all of these mechanisms. Understanding the relative importance of these different mechanisms is an important step in choosing which mechanisms need to be modeled for particular types of cascading failure analysis. This work presents a dynamic simulation model of both power networks and protection systems, which can simulate a wider variety of cascading outage mechanisms, relative to existing quasi-steady state (QSS) models. The model allows one to test the impact of different load models and protections on cascading outage sizes. This paper describes each module of the developed dynamic model and demonstrates how different mechanisms interact. In order to test the model we simulated a batch of randomly selected $N-2$ contingencies for several different static load configurations, and found that the distribution of blackout sizes and event lengths from the proposed dynamic simulator correlates well with historical trends. The results also show that load models have significant impacts on the cascading risks. This dynamic model was also compared against a QSS model based on the dc power flow approximations; we find that the two models largely agree, but produce substantially different results for later stages of cascading.


title: Topologically Massive Spin-1 Particles and Spin-Dependent Potentials

abstract: We investigate the role played by particular field representations of an intermediate massive spin-1 boson in the context of spin-dependent interparticle potentials between fermionic sources in the limit of low momentum transfer. The comparison between the well-known case of the Proca field and that of an exchanged spin-1 boson (with gauge-invariant mass) described by a 2-form potential mixed with a 4-vector gauge field is established in order to pursue an analysis of spin- as well as velocity-dependent profiles of the interparticle potentials. We discuss possible applications and derive an upper bound on the product of vector and pseudo-tensor coupling constants.


title: The Einstein-Maxwell Equations, Kaehler Metrics, and Hermitian Geometry

abstract: Any constant-scalar-curvature Kaehler (cscK) metric on a complex surface may be viewed as a solution of the Einstein-Maxwell equations, and this allows one to produce solutions of these equations on any 4-manifold that arises as a compact complex surface with b_1 even. It is shown, however, that not all solutions of the Einstein-Maxwell equations on such manifolds arise in this way; new examples can be constructed by means of conformally Kaehler geometry.


title: Pseudoholomorphic discs and symplectic structures in Hilbert space

abstract: We develop the theory of $J$-holomorphic discs in Hilbert spaces with almost complex structures. As an aplication, we prove a version of Gromov's symplectic non-squeezing theorem for Hilbert spaces. It can be applied to short-time symplectic flows of a wide class of Hamiltonian PDEs.


title: Coordinated Observations with Pulsar Timing Arrays and ISS-Lobster

abstract: Supermassive black hole binaries are the strongest gravitational wave sources in the universe. The systems most likely to be observed with pulsar timing arrays (PTAs) will have particularly high masses ($\gtrsim 10^9 M_\odot$), long periods ($T_{\rm orb} \gtrsim 1$ yr), and be in the local universe ($z \lesssim 1$). These features are also the most favorable for bright electromagnetic counterparts, which should be easily observable with existing ground- and space-based telescopes. Wide-field X-ray observatories such as ISS-Lobster will provide independent candidates that can be used to lower the threshold for PTA detections of resolvable binary sources. The primary challenge lies in correctly identifying and characterizing binary sources with long orbital periods, as opposed to "normal" active galactic nuclei (AGN) hosting single black holes. Here too ISS-Lobster will provide valuable new understanding into the wide range of behaviors seen in AGN by vastly expanding our sample of X-ray light curves from accreting supermassive black holes.


title: Sub-computable Boundedness Randomness

abstract: This paper defines a new notion of bounded computable randomness for certain classes of sub-computable functions which lack a universal machine. In particular, we define such versions of randomness for primitive recursive functions and for PSPACE functions. These new notions are robust in that there are equivalent formulations in terms of (1) Martin-L\"of tests, (2) Kolmogorov complexity, and (3) martingales. We show these notions can be equivalently defined with prefix-free Kolmogorov complexity. We prove that one direction of van Lambalgen's theorem holds for relative computability, but the other direction fails. We discuss statistical properties of these notions of randomness.


title: Quantized magneto-thermoelectric transport in low-dimensional junctions

abstract: Quantization of the magneto-thermoelectric transport is studied when an external d.c. magnetic field is applied to the C/N-knot formed as crossing between a narrow stripe of conducting atomic monolayer C on the one hand and metal stripe N on the other hand. The temperature gradient in C is created by injecting the non-equilibrium electrons, holes and phonons from the heater H thereby directing them toward the C/N-knot. A non-linear coupling between electron states of the C/N-knot counter electrodes causes splitting of the heat flow into several fractions owing to the Lorentz force acting in the C/N-knot vicinity, thereby inducing the magneto-thermoelectric current in N whereas the phonons pass and propagate along C further ahead. The heat flow along C generates a transversal electric current in N showing a series of maximums when dimensions of the Landau orbits and the C/N-knot match each other. It allows observing the interplay between the quantum Hall effect and the spatial quantization.


title: The numerical relativity breakthrough for binary black holes

abstract: The evolution of black-hole binaries in vacuum spacetimes constitutes the two-body problem in general relativity. The solution of this problem in the framework of the Einstein field equations is a substantially more complex exercise than that of the dynamics of two point masses in Newtonian gravity, but it also presents us with a wealth of new exciting physics. Numerical methods are likely the only method to compute the dynamics of black-hole systems in the fully non-linear regime and have been pursued since the 1960s, culminating in dramatic breakthroughs in 2005. Here we review the methodology and the developments that finally gave us a solution of this fundamental problem of Einstein's theory and discuss the breakthrough's implication for the wide range of contemporary black-hole physics.


title: On-off intermittency and amplitude-phase synchronization in Keplerian shear flows

abstract: We study the development of coherent structures in local simulations of the magnetorotational instability in accretion discs in regimes of on-off intermittency. In a previous paper [Chian et al., Phys. Rev. Lett. 104, 254102 (2010)], we have shown that the laminar and bursty states due to the on-off spatiotemporal intermittency in a one-dimensional model of nonlinear waves correspond, respectively, to nonattracting coherent structures with higher and lower degrees of amplitude-phase synchronization. In this paper we extend these results to a three-dimensional model of magnetized Keplerian shear flows. Keeping the kinetic Reynolds number and the magnetic Prandtl number fixed, we investigate two different intermittent regimes by varying the plasma beta parameter. The first regime is characterized by turbulent patterns interrupted by the recurrent emergence of a large-scale coherent structure known as two-channel flow, where the state of the system can be described by a single Fourier mode. The second regime is dominated by the turbulence with sporadic emergence of coherent structures with shapes that are reminiscent of a perturbed channel flow. By computing the Fourier power and phase spectral entropies in three-dimensions, we show that the large-scale coherent structures are characterized by a high degree of amplitude-phase synchronization.

categories: astro-ph.SR physics.data-an physics.flu-dyn physics.plasm-ph


title: Nonthermal correction to black hole spectroscopy

abstract: Area spectrum of black holes have been obtained via various methods such as quasinormal modes, adiabatic invariance and angular momentum. Among those methods, calculations were done by assuming black holes in thermal equilibrium. Nevertheless, black holes in the asymptotically flat space usually have negative specific heat and therefore tend to stay away from thermal equilibrium. Even for those black holes with positive specific heat, temperature may still not be well defined in the process of radiation, due to the back reaction of decreasing mass. Respect to these facts, it is very likely that Hawking radiation is nonthermal and the area spectrum is no longer equidistant. In this note, we would like to illustrate how the area spectrum of black holes is corrected by this nonthermal effect.


title: How to Scale Up Kernel Methods to Be As Good As Deep Neural Nets

abstract: The computational complexity of kernel methods has often been a major barrier for applying them to large-scale learning problems. We argue that this barrier can be effectively overcome. In particular, we develop methods to scale up kernel models to successfully tackle large-scale learning problems that are so far only approachable by deep learning architectures. Based on the seminal work by Rahimi and Recht on approximating kernel functions with features derived from random projections, we advance the state-of-the-art by proposing methods that can efficiently train models with hundreds of millions of parameters, and learn optimal representations from multiple kernels. We conduct extensive empirical studies on problems from image recognition and automatic speech recognition, and show that the performance of our kernel models matches that of well-engineered deep neural nets (DNNs). To the best of our knowledge, this is the first time that a direct comparison between these two methods on large-scale problems is reported. Our kernel methods have several appealing properties: training with convex optimization, cost for training a single model comparable to DNNs, and significantly reduced total cost due to fewer hyperparameters to tune for model selection. Our contrastive study between these two very different but equally competitive models sheds light on fundamental questions such as how to learn good representations.


title: MultiDark simulations: the story of dark matter halo concentrations and density profiles

abstract: Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of dark matter halos in all the simulations exceeds 60 billion. We find that in order to understand the structure of dark matter halos and to make ~1% accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the traditional ratio of the virial radius to the core radius in the NFW profile. For massive halos the averge density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter alpha in the Einasto approximation. Combining results from different redshifts, masses and cosmologies, we show that halos progress through three stages of evolution. (1) They start as rare density peaks that experience very fast and nearly radial infall. This radial infall brings mass closer to the center producing a high concentrated halo. Here, the halo concentration increases with the increasing halo mass and the concentration is defined by the alpha parameter with nearly constant core radius. Later halos slide into (2) the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. Now the concentration does not depend on halo mass. (3) Once the rate of accretion slows down, halos move into the domain of declining concentration-mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. We provide accurate analytical fits to the numerical results for halo density profiles and concentrations.


title: Collaborative Randomized Beamforming for Phased Array Radio Interferometers

abstract: The Square Kilometre Array (SKA) will form the largest radio telescope ever built and such a huge instrument in the desert poses enormous engineering and logistic challenges. Algorithmic and architectural breakthroughs are needed. Data is collected and processed in groups of antennas before transport for central processing. This processing includes beamforming, primarily so as to reduce the amount of data sent. The principal existing technique points to a region of interest independently of the sky model and how the other stations beamform. We propose a new collaborative beamforming algorithm in order to maximize information captured at the stations (thus reducing the amount of data transported). The method increases the diversity in measurements through randomized beam- forming. We demonstrate through numerical simulation the effectiveness of the method. In particular, we show that randomized beamforming can achieve the same image quality while producing 40% less data when compared to the prevailing method matched beamforming.


title: Decorated marked surfaces II: Intersection numbers and dimensions of Homs

abstract: We study the 3-Calabi-Yau categories $\mathcal{D}$ arising from quivers with potential associated to a decorated marked surface $\mathbf{S}_\bigtriangleup$ introduced by the first author. We prove two conjectures in the prequel, that under a bijection between certain objects in $\mathcal{D}$ and certain arcs in $\mathbf{S}_\bigtriangleup$, the dimensions of morphisms between these objects equal the intersection numbers between the corresponding arcs.


title: Resurgence, Operator Product Expansion, and Remarks on Renormalons in Supersymmetric Yang-Mills

abstract: I discuss similarities and differences between the resurgence program in quantum mechanics and the operator product expansion in strongly coupled Yang-Mills theories. In ${\mathcal N}=1$ super-Yang-Mills renormalons possess peculiar features that make them different from renormalons in QCD. Their conventional QCD interpretation does not seem to be applicable in supersymmetric Yang-Mills in a straightforward manner.


title: A convex formulation for hyperspectral image superresolution via subspace-based regularization

abstract: Hyperspectral remote sensing images (HSIs) usually have high spectral resolution and low spatial resolution. Conversely, multispectral images (MSIs) usually have low spectral and high spatial resolutions. The problem of inferring images which combine the high spectral and high spatial resolutions of HSIs and MSIs, respectively, is a data fusion problem that has been the focus of recent active research due to the increasing availability of HSIs and MSIs retrieved from the same geographical area. We formulate this problem as the minimization of a convex objective function containing two quadratic data-fitting terms and an edge-preserving regularizer. The data-fitting terms account for blur, different resolutions, and additive noise. The regularizer, a form of vector Total Variation, promotes piecewise-smooth solutions with discontinuities aligned across the hyperspectral bands. The downsampling operator accounting for the different spatial resolutions, the non-quadratic and non-smooth nature of the regularizer, and the very large size of the HSI to be estimated lead to a hard optimization problem. We deal with these difficulties by exploiting the fact that HSIs generally "live" in a low-dimensional subspace and by tailoring the Split Augmented Lagrangian Shrinkage Algorithm (SALSA), which is an instance of the Alternating Direction Method of Multipliers (ADMM), to this optimization problem, by means of a convenient variable splitting. The spatial blur and the spectral linear operators linked, respectively, with the HSI and MSI acquisition processes are also estimated, and we obtain an effective algorithm that outperforms the state-of-the-art, as illustrated in a series of experiments with simulated and real-life data.


title: A Discriminative CNN Video Representation for Event Detection

abstract: In this paper, we propose a discriminative video representation for event detection over a large scale video dataset when only limited hardware resources are available. The focus of this paper is to effectively leverage deep Convolutional Neural Networks (CNNs) to advance event detection, where only frame level static descriptors can be extracted by the existing CNN toolkit. This paper makes two contributions to the inference of CNN video representation. First, while average pooling and max pooling have long been the standard approaches to aggregating frame level static features, we show that performance can be significantly improved by taking advantage of an appropriate encoding method. Second, we propose using a set of latent concept descriptors as the frame descriptor, which enriches visual information while keeping it computationally affordable. The integration of the two contributions results in a new state-of-the-art performance in event detection over the largest video datasets. Compared to improved Dense Trajectories, which has been recognized as the best video representation for event detection, our new representation improves the Mean Average Precision (mAP) from 27.6% to 36.8% for the TRECVID MEDTest 14 dataset and from 34.0% to 44.6% for the TRECVID MEDTest 13 dataset. This work is the core part of the winning solution of our CMU-Informedia team in TRECVID MED 2014 competition.


title: Decaying light particles in the SHiP experiment: Signal rate estimates for hidden photons

abstract: For the extension of the Standard Model with light hidden photons we present preliminary estimates of the signal rate expected at the recently proposed fixed target SHiP experiment exploiting the CERN SPS beam of 400 GeV protons.


title: Multyphase solutions to the vector Allen-Cahn equation: Crystalline and other complex symmetric structures

abstract: We present a systematic study of entire symmetric solutions $u:R^n\rightarrow R^m$ of the vector Allen-Cahn equation $\Delta u-W_u(u)=0, x \in R^n$, where $W:R^m\rightarrow R$ is smooth, symmetric, nonnegative with a finite number of zeros and $W_u=(\frac{\partial W}{\partial u_1},\ldots,\frac{\partial W}{\partial u_m})^\top$. We introduce a general notion of equivariance with respect to a homomorphism $f:G\rightarrow\Gamma$ ($G,\Gamma$ reflection groups) and prove two abstract results, concerning the cases of $G$ finite and $G$ discrete, for the existence of equivariant solutions. Our approach is variational and based on a mapping property of the parabolic vector Allen-Cahn equation and on a pointwise estimate for vector minimizers.


title: On the number of large triangles in the Brownian triangulation and fragmentation processes

abstract: The Brownian triangulation is a random compact subset of the unit disk introduced by Aldous. For $\epsilon>0$, let $N(\epsilon)$ be the number of triangles whose sizes (measured in different ways) are greater than $\epsilon$ in the Brownian triangulation. We determine the asymptotic behaviour of $N(\epsilon)$ as $\epsilon \to 0$. To obtain this result, a novel concept of "large" dislocations in fragmentations has been proposed. We develop an approach to study the number of large dislocations which is widely applicable to general self-similar fragmentation processes. This technique enables us to study $N(\epsilon)$ because of a bijection between the triangles in the Brownian triangulation and the dislocations of a certain self-similar fragmentation process. Our method also provides a new way to obtain the law of the length of the longest chord in the Brownian triangulation. We further extend our results to the more general class of geodesic stable laminations introduced by Kortchemski.


title: Ab initio equations of state for hydrogen (H-REOS.3) and helium (He-REOS.3) and their implications for the interior of Brown Dwarfs

abstract: We present new equations of state (EOS) for hydrogen and helium covering a wide range of temperatures from 60 K to 10$^7$ K and densities from $10^{-10}$ g/cm$^3$ to $10^3$ g/cm$^3$. They include an extended set of ab initio EOS data for the strongly correlated quantum regime with an accurate connection to data derived from other approaches for the neighboring regions. We compare linear-mixing isotherms based on our EOS tables with available real-mixture data. A first important astrophysical application of this new EOS data is the calculation of interior models for Jupiter and the comparison with recent results. Secondly, mass-radius relations are calculated for Brown Dwarfs which we compare with predictions derived from the widely used EOS of Saumon, Chabrier and van Horn. Furthermore, we calculate interior models for typical Brown Dwarfs with different masses, namely Corot-3b, Gliese-229b and Corot-15b, and the Giant Planet KOI-889b. The predictions for the central pressures and densities differ by up to 10$\%$ dependent on the EOS used. Our EOS tables are made available in the supplemental material of this paper.


title: Optimal Radio Resource Allocation for Hybrid Traffic in Cellular Networks: Centralized and Distributed Architecture

abstract: Optimal resource allocation is of paramount importance in utilizing the scarce radio spectrum efficiently and provisioning quality of service for miscellaneous user applications, generating hybrid data traffic streams in present-day wireless communications systems. A dynamism of the hybrid traffic stemmed from concurrently running mobile applications with temporally varying usage percentages in addition to subscriber priorities impelled from network providers' perspective necessitate resource allocation schemes assigning the spectrum to the applications accordingly and optimally. This manuscript concocts novel centralized and distributed radio resource allocation optimization problems for hybrid traffic-conveying cellular networks communicating users with simultaneously running multiple delay-tolerant and real-time applications modelled as logarithmic and sigmoidal utility functions, volatile application percent usages, and diverse subscriptions. Casting under a utility proportional fairness entail no lost calls for the proposed modi operandi, for which we substantiate the convexity, devise computationally efficient algorithms catering optimal rates to the applications, and prove a mutual mathematical equivalence. Ultimately, the algorithms performance is evaluated via simulations and discussing germane numerical results.


title: Optimal Radio Resource Allocation for Hybrid Traffic in Cellular Networks: Traffic Analysis and Implementation

abstract: In part I of this paper, a distributed and a centralized architecture for an optimal radio resource allocation aware of the traffic delay-tolerance nature, user subscription type, and application usage variations were developed. In the current article, a transmission overhead analysis of the aforementioned distributed and a centralized architectures is investigated and it is proved that the centralized scheme endures a significantly lower transmission overhead than does the distributed approach. Furthermore, the lower bounds of the transmission overhead for both the centralized and the distributed architectures are derived. Moreover, a sensitivity analysis of the resource allocation procedures of the aforesaid centralized and distributed architectures to the changes in the number of users in the system is presented. Besides, a sensitivity analysis of the centralized and distributed approaches to the temporal changes in application usages are investigated. Ultimately, the transmission overhead and sensitivity relevant statements are verified through appropriate simulations. And last but not the least, a real-world implementation of the resource allocation methods developed in Part I is provided.


title: Electronic structure of Ce2RhIn8 2D heavy Fermion system studied by angle resolved photoemission spectroscopy

abstract: We use angle-resolved photoemission spectroscopy to study heavy fermion superconductor Ce2RhIn8. The Fermi surface is rather complicated and consists of several hole and electron pock- ets. We do not observe kz dispersion of Fermi sheets, which is consistent with 2D character of the electronic structure. Comparison of the ARPES data and band structure calculations points to a localized picture of f electrons. Our findings pave the way for understanding the transport and thermodynamical properties of this material.


title: Does asymmetric dark matter always lead to an anti-neutrino signal?

abstract: Under rather generic assumptions, we show that in the asymmetric dark matter (ADM) scenario, the sign of the B-L asymmetry stored in the dark matter sector and the standard model sector are always the same. One particularly striking consequence of this result is that, when the dark matter decays or annihilates in the present universe, the resulting final state always involves an anti-neutrino. As a concrete example of this, we construct a composite ADM model and explore the feasibility of detecting such an anti-neutrino signal in atmospheric neutrino detectors.


title: The Second Order Pole over Split Quaternions

abstract: This is an addition to a series of papers [FL1, FL2, FL3, FL4], where we develop quaternionic analysis from the point of view of representation theory of the conformal Lie group and its Lie algebra. In this paper we develop split quaternionic analogues of certain results from [FL4]. Thus we introduce a space of functions ${\cal D}^h \oplus {\cal D}^a$ with a natural action of the Lie algebra $\mathfrak{gl}(2,\mathbb H_{\mathbb C}) \simeq \mathfrak{sl}(4,\mathbb C)$, decompose ${\cal D}^h \oplus {\cal D}^a$ into irreducible components and find the $\mathfrak{gl}(2,\mathbb H_{\mathbb C})$-equivariant projectors onto each of these irreducible components.


title: The Theory of Surface Enhanced Hyper Raman Scattering

abstract: The Dipole-Quadrupole theory of Surface Enhanced Hyper Raman Scattering (SEHRS), created by the authors is expounded in details. Peculiarities of the behavior of electromagnetic field on rough metal surfaces are considered. It is demonstrated that there is an enhancement of the dipole and quadrupole light-molecule interaction near the places with a large curvature. The expression for the SEHRS cross-section of symmetrical molecules, which consists of several contributions is obtained. Selection rules for the scattering contributions are obtained and a qualitative classification of the contributions after an enhancement degree is performed. Analysis of experimental spectra of pyrazine and phenazine, and also some another molecules is performed too. It is demonstrated a full coincidence of experimental regularities in these spectra with the theory suggested.


title: On the exponential convergence of the Kaczmarz algorithm

abstract: The Kaczmarz algorithm (KA) is a popular method for solving a system of linear equations. In this note we derive a new exponential convergence result for the KA. The key allowing us to establish the new result is to rewrite the KA in such a way that its solution path can be interpreted as the output from a particular dynamical system. The asymptotic stability results of the corresponding dynamical system can then be leveraged to prove exponential convergence of the KA. The new bound is also compared to existing bounds.


title: A new structure exploiting derivation of recursive direct weight optimization

abstract: The recursive direct weight optimization method is used to solve challenging nonlinear system identification problems. This note provides a new derivation and a new interpretation of the method. The key underlying the note is to acknowledge and exploit a certain structure inherent in the problem.


title: Impact of nanostructure configuration on the photovoltaic performance of quantum dot arrays

abstract: In this work, a mesoscopic model based on the non-equilibrium Green's function formalism for a tight-binding-like effective Hamiltonian is used to investigate a selectively contacted quantum dot array designed for operation as a single junction quantum dot solar cell. By establishing a direct relation between nanostructure configuration and optoelectronic properties, the investigation reveals the influence of inter-dot and dot-contact coupling strengths on the rates of charge carrier photogeneration, radiative recombination, and extraction at contacts, and consequently on the ultimate performance of photovoltaic devices with finite quantum dot arrays as the active medium. For long carrier lifetimes, the dominant configuration effects originate in the dependence of the joint density of states on the inter-dot coupling in terms of band width and effective band gap. In the low carrier lifetime regime, where recombination competes with carrier extraction, the extraction efficiency shows a critical dependence on the dot-contact coupling.


title: LAMP: A Locally Adapting Matching Pursuit Framework for Group Sparse Signatures in Ultra-Wide Band Radar Imaging

abstract: It has been found that radar returns of extended targets are not only sparse but also exhibit a tendency to cluster into randomly located, variable sized groups. However, the standard techniques of Compressive Sensing as applied in radar imaging hardly considers the clustering tendency into account while reconstructing the image from the compressed measurements. If the group sparsity is taken into account, it is intuitive that one might obtain better results both in terms of accuracy and time complexity as compared to the conventional recovery techniques like Orthogonal Matching Pursuit (OMP). In order to remedy this, techniques like Block OMP have been used in the existing literature. An alternate approach is via reconstructing the signal by transforming into the Hough Transform Domain where they become point-wise sparse. However, these techniques essentially assume specific size and structure of the groups and are not always effective if the exact characteristics of the groups are not known, prior to reconstruction. In this manuscript, a novel framework that we call locally adapting matching pursuit (LAMP) have been proposed for efficient reconstruction of group sparse signals from compressed measurements without assuming any specific size, location, or structure of the groups. The recovery guarantee of the LAMP and its superiority compared to the existing algorithms has been established with respect to accuracy, time complexity and flexibility in group size. LAMP has been successfully used on a real-world, experimental data set.


title: Cause-of-death estimates for the early and late neonatal periods for 194 countries from 2000-2013

abstract: Objective: Cause-of-death distributions are important for prioritising interventions. We estimated proportions, risks, and numbers of deaths (with uncertainty) for programme-relevant causes of neonatal death for 194 countries for 2000-2013, differentiating between the early (days 0-6) and late (days 7-27) neonatal periods. Methods: For 65 high-quality VR countries, we used the observed early and late neonatal proportional cause distributions. For the remaining 129 countries, we used multinomial logistic models to estimate the early and late proportional cause distributions. We used separate models, with different inputs, for low and high neonatal mortality countries. We applied these cause-specific proportions to neonatal death estimates from the United Nations by country/year to estimate cause-specific risks and numbers of deaths. Findings: Of the 2.76 million neonatal deaths in 2013, 0.99 (uncertainty: 0.70-1.31) million (35.7%) were estimated to be from preterm complications, 0.64 (uncertainty: 0.46-0.84) million (23.4%) from intrapartum-related complications, and 0.43 (0.22-0.66) million (15.6%) from sepsis. Preterm (40.8%) and intrapartum-related (27.0%) complications accounted for the majority of early neonatal deaths while infections caused nearly half of late neonatal deaths. In every region, preterm was the leading cause of neonatal death, with the highest risks in Southern Asia (11.9 per 1000 livebirths) and Sub-Saharan Africa (9.5). Conclusion: The neonatal cause-of-death distribution differs between the early and late periods, and varies with NMR level and over time. To reduce neonatal deaths, this knowledge must be incorporated into policy decisions. The Every Newborn Action Plan provides stimulus for countries to update national strategies and include high-impact interventions to address these causes.


title: Numeric Invariants from Multidimensional Persistence

abstract: We extend the results of Adcock, Carlsson, and Carlsson by constructing numeric invariants from the computation of a multidimensional persistence module as given by Carlsson, Singh, and Zomorodian.


title: Automatic Generation of Alternative Starting Positions for Simple Traditional Board Games

abstract: Simple board games, like Tic-Tac-Toe and CONNECT-4, play an important role not only in the development of mathematical and logical skills, but also in the emotional and social development. In this paper, we address the problem of generating targeted starting positions for such games. This can facilitate new approaches for bringing novice players to mastery, and also leads to discovery of interesting game variants. We present an approach that generates starting states of varying hardness levels for player~$1$ in a two-player board game, given rules of the board game, the desired number of steps required for player~$1$ to win, and the expertise levels of the two players. Our approach leverages symbolic methods and iterative simulation to efficiently search the extremely large state space. We present experimental results that include discovery of states of varying hardness levels for several simple grid-based board games. The presence of such states for standard game variants like $4 \times 4$ Tic-Tac-Toe opens up new games to be played that have never been played as the default start state is heavily biased.


title: A Discrete Logarithm-based Approach to Compute Low-Weight Multiples of Binary Polynomials

abstract: Being able to compute efficiently a low-weight multiple of a given binary polynomial is often a key ingredient of correlation attacks to LFSR-based stream ciphers. The best known general purpose algorithm is based on the generalized birthday problem. We describe an alternative approach which is based on discrete logarithms and has much lower memory complexity requirements with a comparable time complexity.