title: A criterion for bubble merging in liquid metal: computational and experimental study

abstract: An innovative model is presented for merging of bubbles inside a liquid metal. The proposed model is based on forming a thin film (narrow channel) between merging bubbles during growth. Rupturing of the film occurs when an oscillation in velocity and pressure arises inside the channel followed by merging of the bubbles. The proposed model based on lattice Boltzmann Method is capable of simulating merging bubbles in micro, meso, and macro-scales with no limitation on the number of bubbles. Experimental studies reveal a good consistency between modeling results and real conditions.


title: Deep searches for broadband extended gravitational-wave emission bursts by heterogeneous computing

abstract: We present a heterogeneous search algorithm for broadband extended gravitational-wave emission (BEGE), expected from gamma-ray bursts and energetic core-collapse supernovae. It searches the $(f,\dot{f})$-plane for long duration bursts by inner engines slowly exhausting their energy reservoir by matched filtering on a {\em Graphics Processor Unit} (GPU) over a template bank of millions of one-second duration chirps. Parseval's Theorem is used to predict the standard deviation $\sigma$ of filter output, taking advantage of near-Gaussian noise in LIGO S6 data over 350-2000 Hz. Tails exceeding a mulitple of $\sigma$ are communicated back to a {\em Central Processing Unit} (CPU). This algorithm attains about 65\% efficiency overall, normalized to the Fast Fourier Transform (FFT). At about one million correlations per second over data segments of 16 s duration $(N=2^{16}$ samples), better than real-time analysis is achieved on a cluster of about a dozen GPUs. We demonstrate its application to the capture of high frequency hardware LIGO injections. This algorithm serves as a starting point for deep all-sky searches in both archive data and real-time analysis in current observational runs.


title: A novel X-FEM based fast computational method for crack propagation

abstract: This study suggests a fast computational method for crack propagation, which is based on the extended finite element method (X-FEM). It is well known that the X-FEM might be the most popular numerical method for crack propagation. However, with the increase of complexity of the given problem, the size of FE model and the number of iterative steps are increased correspondingly. To improve the efficiency of X-FEM, an efficient computational method termed decomposed updating reanalysis (DUR) method is suggested. For most of X-FEM simulation procedures, the change of each iterative step is small and it will only lead a local change of stiffness matrix. Therefore, the DUR method is proposed to predict the modified response by only calculating the changed part of equilibrium equations. Compared with other fast computational methods, the distinctive characteristic of the proposed method is to update the modified stiffness matrix with a local updating strategy, which only the changed part of stiffness matrix needs to be updated. To verify the performance of the DUR method, several typical numerical examples have been analyzed and the results demonstrate that this method is a highly efficient method with high accuracy.


title: Identification of Probabilities

abstract: Within psychology, neuroscience and artificial intelligence, there has been increasing interest in the proposal that the brain builds probabilistic models of sensory and linguistic input: that is, to infer a probabilistic model from a sample. The practical problems of such inference are substantial: the brain has limited data and restricted computational resources. But there is a more fundamental question: is the problem of inferring a probabilistic model from a sample possible even in principle? We explore this question and find some surprisingly positive and general results. First, for a broad class of probability distributions characterised by computability restrictions, we specify a learning algorithm that will almost surely identify a probability distribution in the limit given a finite i.i.d. sample of sufficient but unknown length. This is similarly shown to hold for sequences generated by a broad class of Markov chains, subject to computability assumptions. The technical tool is the strong law of large numbers. Second, for a large class of dependent sequences, we specify an algorithm which identifies in the limit a computable measure for which the sequence is typical, in the sense of Martin-Lof (there may be more than one such measure). The technical tool is the theory of Kolmogorov complexity. We analyse the associated predictions in both cases. We also briefly consider special cases, including language learning, and wider theoretical implications for psychology.


title: Turbulence Momentum Transport and Prediction of the Reynolds Stress in Canonical Flows

abstract: We present a unique method for solving for the Reynolds stress in turbulent canonical flows, based on the momentum balance for a control volume moving at the local mean velocity. A differential transform converts this momentum balance to a solvable form. Comparisons with experimental and computational data in simple geometries show quite good agreements. An alternate picture for the turbulence momentum transport is offered, as verified with data, where the turbulence momentum is transported by the mean velocity while being dissipated by viscosity. The net momentum transport is the Reynolds stress. This turbulence momentum balance is verified using DNS and experimental data. Implications of this work are that the Reynolds stress can be written explicitly in terms of basic turbulence parameters in a simple form, derived from pure fluid physics, and that potential exists for applications of this concept to more complex turbulent flows.


title: Multiphase Aluminum A356 Foam Formation Process Simulation Using Lattice Boltzmann Method

abstract: Shan-Chen model is a numerical scheme to simulate multiphase fluid flows using Lattice Boltzmann approach. The original Shan-Chen model suffers from inability to accurately predict behavior of air bubbles interacting in a non-aqueous fluid. In the present study, we extended the Shan-Chen model to take the effect of the attraction-repulsion barriers among bubbles in to account. The proposed model corrects the interaction and coalescence criterion of the original Shan-Chen scheme in order to have a more accurate simulation of bubbles morphology in a metal foam. The model is based on forming a thin film (narrow channel) between merging bubbles during growth. Rupturing of the film occurs when an oscillation in velocity and pressure arises inside the channel followed by merging of the bubbles. Comparing numerical results obtained from proposed model with mettallorgraphy images for aluminum A356 demonstrated a good consistency in mean bubble size and bubbles distribution


title: Virtual signatures of dark sectors in Higgs couplings

abstract: Where collider searches for resonant invisible particles loose steam, dark sectors might leave their trace as virtual effects in precision observables. Here we explore this option in the framework of Higgs portal models, where a sector of dark fermions interacts with the standard model through a strong renormalizable coupling to the Higgs boson. We show that precise measurements of Higgs-gauge and triple Higgs interactions can probe dark fermions up to the TeV scale through virtual corrections. Observation prospects at the LHC and future lepton colliders are discussed for the so-called singlet-doublet model of Majorana fermions, a generalization of the bino-higgsino scenario in supersymmetry. We advocate a two-fold search strategy for dark sectors through direct and indirect observables.


title: Galaxies with Shells in the Illustris Simulation: Metallicity Signatures

abstract: Stellar shells are low surface brightness arcs of overdense stellar regions, extending to large galactocentric distances. In a companion study, we identified 39 shell galaxies in a sample of 220 massive ellipticals ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_\odot$) from the Illustris cosmological simulation. We used stellar history catalogs to trace the history of each individual star particle inside the shell substructures, and we found that shells in high-mass galaxies form through mergers with massive satellites (stellar mass ratios $\mu_{\mathrm{stars}}\gtrsim1:10$). Using the same sample of shell galaxies, the current study extends the stellar history catalogs in order to investigate the metallicity of stellar shells around massive galaxies. Our results indicate that outer shells are often times more metal-rich than the surrounding stellar material in a galaxy's halo. For a galaxy with two different satellites forming $z=0$ shells, we find a significant difference in the metallicity of the shells produced by each progenitor. We also find that shell galaxies have higher mass-weighted logarithmic metallicities ([Z/H]) at $2$-$4\,\mathrm{R}_{\mathrm{eff}}$ compared to galaxies without shells. Our results indicate that observations comparing the metallicities of stars in tidal features, such as shells, to the average metallicities in the stellar halo can provide information about the assembly histories of galaxies.


title: Spatially resolved spectroscopy across stellar surfaces. I. Using exoplanet transits to analyze 3-D stellar atmospheres

abstract: CONTEXT: High-precision stellar analyses require hydrodynamic modeling to interpret chemical abundances or oscillation modes. Exoplanet atmosphere studies require stellar background spectra to be known along the transit path while detection of Earth analogs require stellar microvariability to be understood. Hydrodynamic 3-D models can be computed for widely different stars but have been tested in detail only for the Sun with its resolved surface features. Model predictions include spectral line shapes, asymmetries, and wavelength shifts, and their center-to-limb changes across stellar disks. AIMS: To observe high-resolution spectral line profiles across spatially highly resolved stellar surfaces, which are free from the effects of spatial smearing and rotational broadening present in full-disk spectra, enabling comparisons to synthetic profiles from 3-D models. METHODS: During exoplanet transits, successive stellar surface portions become hidden and differential spectroscopy between various transit phases provides spectra of small surface segments temporarily hidden behind the planet. Planets cover no more than about 1% of any main-sequence star, enabling high spatial resolution but demanding very precise observations. Realistically measurable quantities are identified through simulated observations of synthetic spectral lines. RESULTS: In normal stars, line profile ratios between various transit phases may vary by some 0.5%, requiring S/N ratios of 5,000 or more for meaningful spectral reconstruction. While not yet realistic for individual spectral lines, this is achievable for cool stars by averaging over numerous lines with similar parameters. CONCLUSIONS: For bright host stars of large transiting planets, spatially resolved spectroscopy is currently practical. More observable targets are likely to be found in the near future by ongoing photometric searches.


title: Large Synoptic Survey Telescope Galaxies Science Roadmap

abstract: The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.


title: Spatially resolved spectroscopy across stellar surfaces. II. High-resolution spectra across HD209458 (G0V)

abstract: CONTEXT: High-resolution spectroscopy across spatially resolved stellar surfaces aims at obtaining spectral-line profiles that are free from rotational broadening; the gradual changes of these profiles from disk center toward the stellar limb reveal properties of atmospheric fine structure, which are possible to model with 3-D hydrodynamics. AIMS: Previous such studies have only been carried out for the Sun but are now extended to other stars. In this work, profiles of photospheric spectral lines are retrieved across the disk of the planet-hosting star HD209458 (G0V). METHODS: During exoplanet transit, stellar surface portions successively become hidden and differential spectroscopy provides spectra of small surface segments temporarily hidden behind the planet. The method was elaborated in Paper I, with observable signatures quantitatively predicted from hydrodynamic simulations. RESULTS: From observations of HD209458 with spectral resolution R=80,000, photospheric FeI line profiles are obtained at several center-to-limb positions, reaching adequately high S/N after averaging over numerous similar lines CONCLUSIONS: Retrieved line profiles are compared to synthetic line profiles. Hydrodynamic 3-D models predict, and current observations confirm, that photospheric absorption lines become broader and shallower toward the stellar limb, reflecting that horizontal velocities in stellar granulation are greater than vertical velocities. Additional types of 3-D signatures will become observable with the highest resolution spectrometers at large telescopes.


title: Hydrogen-Poor Superluminous Supernovae from the Pan-STARRS1 Medium Deep Survey

abstract: We present light curves and classification spectra of 17 hydrogen-poor superluminous supernovae (SLSNe) from the Pan-STARRS1 Medium Deep Survey (PS1 MDS). Our sample contains all objects from the PS1 MDS sample with spectroscopic classification that are similar to either of the prototypes SN2005ap or SN2007bi, without an explicit limit on luminosity. With a redshift range $0.3 < z < 1.6$, PS1MDS is the first SLSN sample primarily probing the high-redshift population; our multi-filter PS1 light curves probe the rest-frame UV emission, and hence the peak of the spectral energy distribution. We measure the temperature evolution and construct bolometric light curves, and find peak luminosities of $(0.5-5) \times 10^{44}$ erg s$^{-1}$ and lower limits on the total radiated energies of $(0.3-2) \times 10^{51}$ erg. The light curve shapes are diverse, with both rise- and decline times spanning a factor of $\sim 5$, and several examples of double-peaked light curves. When correcting for the flux-limited nature of our survey, we find a median peak luminosity at 4000 {\AA} of $M_{\rm 4000} = -21.1$ mag, and a spread of $\sigma = 0.7$ mag.


title: Exponentially Slow Heating in Short and Long-range Interacting Floquet Systems

abstract: We analyze the dynamics of periodically-driven (Floquet) Hamiltonians with short- and long-range interactions, finding clear evidence for a thermalization time, $\tau^*$, that increases exponentially with the drive frequency. We observe this behavior, both in systems with short-ranged interactions, where our results are consistent with rigorous bounds, and in systems with long-range interactions, where such bounds do not exist at present. Using a combination of heating and entanglement dynamics, we explicitly extract the effective energy scale controlling the rate of thermalization. Finally, we demonstrate that for times shorter than $\tau^*$, the dynamics of the system is well-approximated by evolution under a time-independent Hamiltonian $D_{\mathrm{eff}}$, for both short- and long-range interacting systems.

categories: quant-ph cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el


title: Core-powered mass loss and the radius distribution of small exoplanets

abstract: Recent observations identify a valley in the radius distribution of small exoplanets, with planets in the range $1.5-2.0\,{\rm R}_{\oplus}$ significantly less common than somewhat smaller or larger planets. This valley may suggest a bimodal population of rocky planets that are either engulfed by massive gas envelopes that significantly enlarge their radius, or do not have detectable atmospheres at all. One explanation of such a bimodal distribution is atmospheric erosion by high-energy stellar photons. We investigate an alternative mechanism: the luminosity of the cooling rocky core, which can completely erode light envelopes while preserving heavy ones, produces a deficit of intermediate sized planets. We evolve planetary populations that are derived from observations using a simple analytical prescription, accounting self-consistently for envelope accretion, cooling and mass loss, and demonstrate that core-powered mass loss naturally reproduces the observed radius distribution, regardless of the high-energy incident flux. Observations of planets around different stellar types may distinguish between photoevaporation, which is powered by the high-energy tail of the stellar radiation, and core-powered mass loss, which depends on the bolometric flux through the planet's equilibrium temperature that sets both its cooling and mass-loss rates.


title: The structure of young embedded protostellar discs

abstract: Young protostellar discs provide the initial conditions for planet formation. The properties of these discs may be different from those of late-phase (T Tauri) discs due to continuing infall from the envelope and protostellar variability resulting from irregular gas accretion. We use a set of hydrodynamic simulations to determine the structure of discs forming in collapsing molecular clouds. We examine how radiative feedback from the host protostar affects the disc properties by examining three regimes: without radiative feedback, with continuous radiative feedback and with episodic feedback, similar to FU Ori-type outbursts. We find that the radial surface density and temperature profiles vary significantly as the disc accretes gas from the infalling envelope. These profiles are sensitive to the presence of spiral structure, induced by gravitational instabilities, and the radiative feedback provided by the protostar, especially in the case when the feedback is episodic. We also investigate whether mass estimates from position-velocity (PV) diagrams are accurate for early-phase discs. We find that the protostellar system mass (i.e. the mass of the protostar and its disc) is underestimated by up to 20%, due to the impact of an enhanced radial pressure gradient on the gas. The mass of early-phase discs is a significant fraction of the mass of the protostar, so position-velocity diagrams cannot accurately provide the mass of the protostar alone. The enhanced radial pressure gradient expected in young discs may lead to an increased rate of dust depletion due to gas drag, and therefore to a reduced dust-to-gas ratio.


title: Light curves of hydrogen-poor Superluminous Supernovae from the Palomar Transient Factory

abstract: We investigate the light-curve properties of a sample of 26 spectroscopically confirmed hydrogen-poor superluminous supernovae (SLSNe-I) in the Palomar Transient Factory (PTF) survey. These events are brighter than SNe Ib/c and SNe Ic-BL, on average, by about 4 and 2~mag, respectively. The peak absolute magnitudes of SLSNe-I in rest-frame $g$ band span $-22\lesssim M_g \lesssim-20$~mag, and these peaks are not powered by radioactive $^{56}$Ni, unless strong asymmetries are at play. The rise timescales are longer for SLSNe than for normal SNe Ib/c, by roughly 10 days, for events with similar decay times. Thus, SLSNe-I can be considered as a separate population based on photometric properties. After peak, SLSNe-I decay with a wide range of slopes, with no obvious gap between rapidly declining and slowly declining events. The latter events show more irregularities (bumps) in the light curves at all times. At late times, the SLSN-I light curves slow down and cluster around the $^{56}$Co radioactive decay rate. Powering the late-time light curves with radioactive decay would require between 1 and 10${\rm M}_\odot$ of Ni masses. Alternatively, a simple magnetar model can reasonably fit the majority of SLSNe-I light curves, with four exceptions, and can mimic the radioactive decay of $^{56}$Co, up to $\sim400$ days from explosion. The resulting spin values do not correlate with the host-galaxy metallicities. Finally, the analysis of our sample cannot strengthen the case for using SLSNe-I for cosmology.


title: Signal Formation Processes in Micromegas Detectors and Quality Control for large size Detector Construction for the ATLAS New Small Wheel

abstract: The Micromegas technology is one of the most successful MPGD concepts and widely utilized in nuclear and particle physics experiments. Twenty years of research and development rendered the technology sufficiently mature to be selected as precision tracking detector for the New Small Wheel (NSW) upgrade of the ATLAS Muon spectrometer. This will be the first large scale application of Micromegas in one of the major LHC experiments. The studies on signal formation in Micromegas focuses on the microscopic signal electron loss mechanisms and the amplification processes in electron gas interaction. Based on a model of detector parameter dependencies, these processes are scrutinized in an iterating comparison between exper- imental results, theory prediction of the macroscopic observables and process simulation on the microscopic level. Utilizing the specialized detectors developed in the scope of this thesis as well as refined simulation algorithms, an unprecedented level of accuracy in the process description is reached, deepening the understanding of the fundamental process in gaseous detectors. The second part is dedicated to the challenges arising with the large scale Micro- megas production for the ATLAS NSW. A selection of technological choices, are discussed alongside a final report on two production related tasks addressing the detectors core components: For the industrial production of resistive anode PCBs a detailed quality control (QC) and quality assurance (QA) scheme as well as the therefore required testing tools have been developed. In parallel the study on micromesh parameter optimization and production feasibility resulted in the selection of the proposed mesh by the NSW community. The successful completion of both tasks were important milestones towards the construction of large size Micromegas detectors clearing the path for NSW series production.


title: Traffic flow optimization using a quantum annealer

abstract: Quantum annealing algorithms belong to the class of meta-heuristic tools, applicable for solving binary optimization problems. Hardware implementations of quantum annealing, such as the quantum processing units (QPUs) produced by D-Wave Systems, have been subject to multiple analyses in research, with the aim of characterizing the technology's usefulness for optimization and sampling tasks. In this paper, we present a real-world application that uses quantum technologies. Specifically, we show how to map certain parts of the real-world traffic flow optimization problem to be suitable for quantum annealing. We show that time-critical optimization tasks, such as continuous redistribution of position data for cars in dense road networks, are suitable candidates for quantum applications. Due to the limited size and connectivity of current-generation D-Wave QPUs, we use a hybrid quantum and classical approach to solve the traffic flow problem.


title: Expansion of Kes 73, a Shell Supernova Remnant Containing a Magnetar

abstract: Of the 30 or so Galactic magnetars, about 8 are in supernova remnants (SNRs). One of the most extreme magnetars, 1E 1841-045, is at the center of the SNR Kes 73 (G27.4+0.0), whose age is uncertain. We measure its expansion using three Chandra observations over 15 yr, obtaining a mean rate of 0.023% +/- 0.002% per yr. For a distance of 8.5 kpc, we obtain a shell velocity of 1100 km/s and infer a blast-wave speed of 1400 km/s. For Sedov expansion into a uniform medium, this gives an age of 1800 yr. Derived emission measures imply an ambient density of about 2 cm$^{-3}$ and an upper limit on the swept-up mass of about 70 solar masses, with lower limits of tens of solar masses, confirming that Kes 73 is in an advanced evolutionary stage. Our spectral analysis shows no evidence for enhanced abundances as would be expected from a massive progenitor. Our derived total energy is $1.9 \times 10^{51}$ erg, giving a very conservative lower limit to the magnetar's initial period of about 3 ms, unless its energy was lost by non-electromagnetic means. We see no evidence of a wind-blown bubble as would be produced by a massive progenitor, or any evidence that the progenitor of Kes 73/1E 1841-045 was anything but a normal red supergiant producing a Type IIP supernova, though a short-lived stripped-envelope progenitor cannot be absolutely excluded. Kes 73's magnetar thus joins SGR 1900+14 as magnetars resulting from relatively low-mass progenitors.


title: Probing the formation history of the nuclear star cluster at the Galactic Centre with millisecond pulsars

abstract: The origin of the Nuclear Star Cluster in the centre of our Galaxy is still unknown. One possibility is that it formed after the disruption of stellar clusters that spiralled into the Galactic Centre due to dynamical friction. We trace the formation of the Nuclear Star Cluster around the central black hole, using state-of-the-art N-body simulations, and follow the dynamics of the neutron stars born in the clusters. We then estimate the number of Millisecond Pulsars (MSPs) that are released in the Nuclear Star Cluster, during its formation. The assembly and tidal dismemberment of globular clusters lead to a population of MSPs distributed over a radius of about 20 pc, with a peak near 3 pc. No clustering is found on the sub-parsec scale. We simulate the detectability of this population with future radio telescopes like the MeerKAT radio telescope and SKA1, and find that about of order ten MSPs can be observed over this large volume, with a paucity of MSPs within the central parsec. This helps discriminating this scenario from the in-situ formation model for the Nuclear Star Cluster that would predict an over abundance of MSPs closer to the black hole. We then discuss the potential contribution of our MSP population to the gamma-ray excess at the Galactic Centre.


title: Validation of Enhanced Emotion Enabled Cognitive Agent Using Virtual Overlay Multi-Agent System Approach

abstract: Making roads safer by avoiding road collisions is one of the main reasons for inventing Autonomous vehicles (AVs). In this context, designing agent-based collision avoidance components of AVs which truly represent human cognition and emotions look is a more feasible approach as agents can replace human drivers. However, to the best of our knowledge, very few human emotion and cognition-inspired agent-based studies have previously been conducted in this domain. Furthermore, these agent-based solutions have not been validated using any key validation technique. Keeping in view this lack of validation practices, we have selected state-of-the-art Emotion Enabled Cognitive Agent (EEC_Agent), which was proposed to avoid lateral collisions between semi-AVs. The architecture of EEC_Agent has been revised using Exploratory Agent Based Modeling (EABM) level of the Cognitive Agent Based Computing (CABC) framework and real-time fear emotion generation mechanism using the Ortony, Clore & Collins (OCC) model has also been introduced. Then the proposed fear generation mechanism has been validated using the Validated Agent Based Modeling level of CABC framework using a Virtual Overlay MultiAgent System (VOMAS). Extensive simulation and practical experiments demonstrate that the Enhanced EEC_Agent exhibits the capability to feel different levels of fear, according to different traffic situations and also needs a smaller Stopping Sight Distance (SSD) and Overtaking Sight Distance (OSD) as compared to human drivers.


title: The s-process nucleosynthesis: impact of the uncertainties in the nuclear physics determined by Monte Carlo variations

abstract: We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars and massive stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. We found beta-decay rate uncertainties affect only a few nuclei near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in the neutron capture rates. We suggest a list of uncertain rates as candidates for improved measurement by future experiments.


title: Probing the Sea Quark Content of the Proton with One-Particle-Inclusive Processes

abstract: We investigate the feasibility of constraining parton distribution functions in the proton through a comparison with data on semi-inclusive deep-inelastic lepton-nucleon scattering. Specifically, we reweight replicas of these distributions according to how well they reproduce recent, very precise charged kaon multiplicity measurements and analyze how this procedure optimizes the determination of the sea quark densities and improves their uncertainties. The results can help to shed new light on the long standing question on the size of the flavor and charge symmetry breaking among quarks of radiative origin. An iterative method is proposed and adopted to account for the inevitable correlation with what is assumed about the parton-to-hadron fragmentation functions in the reweighting procedure. It is shown how the fragmentation functions can be optimized simultaneously in each step of the iteration. As a first case study, we implement this method to analyze kaon production data.


title: Missing Link Identifcation Using SIS Epidemic Traces

abstract: The study of SIS epidemics on networks has stressed the role of the network topology on the spreading process. However, accurate models of SIS epidemics rely on the complete knowledge of the network topology, which is often not available. This paper tackles the problem of inferring the network topology from observed infection time traces, especially where the network topology is partially known or known with some uncertainty. We propose a Bayesian method to infer the posterior probability of uncertain links in the network, and we derive closed form equations for these probabilities. We also propose a numerical approach based on a Gibbs sampling when the number of uncertain links is large such that using the closed form equations becomes impractical. Numerical results show the capability of the proposed approach to assign high probability to existing links and low probability to non-existing links of the network when the SIS traces are sufficiently long.


title: Localization of Electrical Flows

abstract: We show that in any graph, the average length of a flow path in an electrical flow between the endpoints of a random edge is $O(\log^2 n)$. This is a consequence of a more general result which shows that the spectral norm of the entrywise absolute value of the transfer impedance matrix of a graph is $O(\log^2 n)$. This result implies a simple oblivious routing scheme based on electrical flows in the case of transitive graphs.


title: Constructing Types in Differentially Closed Fields that are Analysable in the Constants

abstract: Analysability of finite $U$-rank types are explored both in general and in the theory $\mathrm{DCF}_0$. The well-known fact that the equation $\delta(\mathrm{log}\delta x)=0$ is analysable in but not almost internal to the constants is generalized to show that $\underbrace{\mathrm{log}\delta...\mathrm{log}\delta}_n x=0$ is not analysable in the constants in $(n-1)$-steps. The notion of a \emph{canonical analysis} is introduced -- namely an analysis that is of minimal length and interalgebraic with every other analysis of that length. Not every analysable type admits a canonical analysis. Using properties of reductions and coreductions in theories with the canonical base property, it is constructed, for any sequence of positive integers $(n_1,...,n_\ell)$, a type in $\mathrm{DCF}_0$ that admits a canonical analysis with the property that the $i$th step has $U$-rank $n_i$.


title: Proportional Mean Residual Life Model with Censored Survival Data under Case-cohort Design

abstract: Proportional mean residual life model is studied for analysing survival data from the case-cohort design. To simultaneously estimate the regression parameters and the baseline mean residual life function, weighted estimating equations based on an inverse selection probability are proposed. The resulting regression coefficients estimates are shown to be consistent and asymptotic normal with easily estimated variance-covariance. Simulation studies show that the proposed estimators perform very well. An application to a real dataset from the South Welsh nickel refiners study is also given to illustrate the methodology.


title: Eruptive Event Generator Based on the Gibson-Low Magnetic Configuration

abstract: Coronal Mass Ejections (CMEs), a kind of energetic solar eruptions, are an integral subject of space weather research. Numerical magnetohydrodynamic (MHD) modeling, which requires powerful computational resources, is one of the primary means of studying the phenomenon. With increasing accessibility of such resources, grows the demand for user-friendly tools that would facilitate the process of simulating CMEs for scientific and operational purposes. The Eruptive Event Generator based on Gibson-Low flux rope (EEGGL), a new publicly available computational model presented in this paper, is an effort to meet this demand. EEGGL allows one to compute the parameters of a model flux rope driving a CME via an intuitive graphical user interface (GUI). We provide a brief overview of the physical principles behind EEGGL and its functionality. Ways towards future improvements of the tool are outlined.


title: Game theory models for communication between agents: a review

abstract: In the real world, agents or entities are in a continuous state of interactions. These inter- actions lead to various types of complexity dynamics. One key difficulty in the study of complex agent interactions is the difficulty of modeling agent communication on the basis of rewards. Game theory offers a perspective of analysis and modeling these interactions. Previously, while a large amount of literature is available on game theory, most of it is from specific domains and does not cater for the concepts from an agent- based perspective. Here in this paper, we present a comprehensive multidisciplinary state-of-the-art review and taxonomy of game theory models of complex interactions between agents.


title: Matrix biorthogonal polynomials, associated polynomials and functions of the second kind

abstract: In this work the interplay between matrix biorthogonal polynomials with respect to a matrix of linear functionals, the $k$-th associated matrix polynomials and the second kind matrix functions, is studied in terms of quasideterminants. A sort of Poincar\'{e}'s theorem for the ratio of two consecutive matrix functions solutions of a linear difference equation is also presented. Some new formulas connecting these families of matrix functions are given.


title: Ratio asymptotic for bi-orthogonal matrix polynomials with unbounded recurrence coefficients

abstract: In this work is presented a study on matrix biorthogonal polynomials sequences that satisfy a nonsymmetric recurrence relation with unbounded coefficients. The ratio asymptotic for this family of matrix biorthogonal polynomials is derived in quite general assumptions. It is considered some illustrative examples.


title: Treatment of Reactive Routing Protocols Using Second Chance Based On Malicious Behavior of Nodes in MANETS

abstract: Mobile nodes of various routing protocols in Mobile Ad hoc Networks follow different strategies in transmission and receiving of data. Security, packet delivery and routing overhead are important concerns for any protocol during designing them. The presence and absence of malicious nodes in the network affect a lot on the performance of the protocol. This research is mainly focused on the study of the threats, attacks and reasons for malicious behavior of nodes in the network for reactive routing protocols in MANETS. DSR and AODV are the two reactive routing protocols that were considered for the study to propose a second chance strategy to be given to the nodes considering the reason for malicious behavior to improve the packet delivery ratio and reduce the routing overhead in the network. A simulative study has been conducted using Ad hoc Simulator (ASIM) considering the DSR and AODV routing protocols in the presence of malicious nodes and in the absence of malicious nodes which showed that the packet delivery ratio is low and routing overhead is high in the absence of malicious nodes. The second chance strategy proposed considers the reasons for malicious behavior and helps the node to get reintegrated in the network to improve the packet delivery ratio and reduce the routing overhead.


title: Speech-driven Animation with Meaningful Behaviors

abstract: Conversational agents (CAs) play an important role in human computer interaction. Creating believable movements for CAs is challenging, since the movements have to be meaningful and natural, reflecting the coupling between gestures and speech. Studies in the past have mainly relied on rule-based or data-driven approaches. Rule-based methods focus on creating meaningful behaviors conveying the underlying message, but the gestures cannot be easily synchronized with speech. Data-driven approaches, especially speech-driven models, can capture the relationship between speech and gestures. However, they create behaviors disregarding the meaning of the message. This study proposes to bridge the gap between these two approaches overcoming their limitations. The approach builds a dynamic Bayesian network (DBN), where a discrete variable is added to constrain the behaviors on the underlying constraint. The study implements and evaluates the approach with two constraints: discourse functions and prototypical behaviors. By constraining on the discourse functions (e.g., questions), the model learns the characteristic behaviors associated with a given discourse class learning the rules from the data. By constraining on prototypical behaviors (e.g., head nods), the approach can be embedded in a rule-based system as a behavior realizer creating trajectories that are timely synchronized with speech. The study proposes a DBN structure and a training approach that (1) models the cause-effect relationship between the constraint and the gestures, (2) initializes the state configuration models increasing the range of the generated behaviors, and (3) captures the differences in the behaviors across constraints by enforcing sparse transitions between shared and exclusive states per constraint. Objective and subjective evaluations demonstrate the benefits of the proposed approach over an unconstrained model.


title: Localizing Moments in Video with Natural Language

abstract: We consider retrieving a specific temporal segment, or moment, from a video given a natural language text description. Methods designed to retrieve whole video clips with natural language determine what occurs in a video but not when. To address this issue, we propose the Moment Context Network (MCN) which effectively localizes natural language queries in videos by integrating local and global video features over time. A key obstacle to training our MCN model is that current video datasets do not include pairs of localized video segments and referring expressions, or text descriptions which uniquely identify a corresponding moment. Therefore, we collect the Distinct Describable Moments (DiDeMo) dataset which consists of over 10,000 unedited, personal videos in diverse visual settings with pairs of localized video segments and referring expressions. We demonstrate that MCN outperforms several baseline methods and believe that our initial results together with the release of DiDeMo will inspire further research on localizing video moments with natural language.


title: Cut, Paste and Learn: Surprisingly Easy Synthesis for Instance Detection

abstract: A major impediment in rapidly deploying object detection models for instance detection is the lack of large annotated datasets. For example, finding a large labeled dataset containing instances in a particular kitchen is unlikely. Each new environment with new instances requires expensive data collection and annotation. In this paper, we propose a simple approach to generate large annotated instance datasets with minimal effort. Our key insight is that ensuring only patch-level realism provides enough training signal for current object detector models. We automatically `cut' object instances and `paste' them on random backgrounds. A naive way to do this results in pixel artifacts which result in poor performance for trained models. We show how to make detectors ignore these artifacts during training and generate data that gives competitive performance on real data. Our method outperforms existing synthesis approaches and when combined with real images improves relative performance by more than 21% on benchmark datasets. In a cross-domain setting, our synthetic data combined with just 10% real data outperforms models trained on all real data.


title: Ensuring patients privacy in a cryptographic-based-electronic health records using bio-cryptography

abstract: Several recent works have proposed and implemented cryptography as a means to preserve privacy and security of patients health data. Nevertheless, the weakest point of electronic health record (EHR) systems that relied on these cryptographic schemes is key management. Thus, this paper presents the development of privacy and security system for cryptography-based-EHR by taking advantage of the uniqueness of fingerprint and iris characteristic features to secure cryptographic keys in a bio-cryptography framework. The results of the system evaluation showed significant improvements in terms of time efficiency of this approach to cryptographic-based-EHR. Both the fuzzy vault and fuzzy commitment demonstrated false acceptance rate (FAR) of 0%, which reduces the likelihood of imposters gaining successful access to the keys protecting patients protected health information. This result also justifies the feasibility of implementing fuzzy key binding scheme in real applications, especially fuzzy vault which demonstrated a better performance during key reconstruction.


title: The modal logic of set-theoretic potentialism and the potentialist maximality principles

abstract: We analyze the precise modal commitments of several natural varieties of set-theoretic potentialism, using tools we develop for a general model-theoretic account of potentialism, building on those of Hamkins, Leibman and L\"owe, including the use of buttons, switches, dials and ratchets. Among the potentialist conceptions we consider are: rank potentialism (true in all larger $V_\beta$); Grothendieck-Zermelo potentialism (true in all larger $V_\kappa$ for inaccessible cardinals $\kappa$); transitive-set potentialism (true in all larger transitive sets); forcing potentialism (true in all forcing extensions); countable-transitive-model potentialism (true in all larger countable transitive models of ZFC); countable-model potentialism (true in all larger countable models of ZFC); and others. In each case, we identify lower bounds for the modal validities, which are generally either S4.2 or S4.3, and an upper bound of S5, proving in each case that these bounds are optimal. The validity of S5 in a world is a potentialist maximality principle, an interesting set-theoretic principle of its own. The results can be viewed as providing an analysis of the modal commitments of the various set-theoretic multiverse conceptions corresponding to each potentialist account.


title: Existence of locally maximally entangled quantum states via geometric invariant theory

abstract: We study a question which has natural interpretations in both quantum mechanics and in geometry. Let $V_1,..., V_n$ be complex vector spaces of dimension $d_1,...,d_n$ and let $G= SL_{d_1} \times \dots \times SL_{d_n}$. Geometrically, we ask given $(d_1,...,d_n)$, when is the geometric invariant theory quotient $\mathbb{P}(V_1 \otimes \dots \otimes V_n)// G$ non-empty? This is equivalent to the quantum mechanical question of whether the multipart quantum system with Hilbert space $V_1\otimes \dots \otimes V_n$ has a locally maximally entangled state, i.e. a state such that the density matrix for each elementary subsystem is a multiple of the identity. We show that the answer to this question is yes if and only if $R(d_1,...,d_n)\geqslant 0$ where \[ R(d_1,...,d_n) = \prod_i d_i +\sum_{k=1}^n (-1)^k \sum_{1\leq i_1<\dotsb <i_k\leq n} (\gcd(d_{i_1},\dotsc ,d_{i_k}) )^{2}. \] We also provide a simple recursive algorithm which determines the answer to the question, and we compute the dimension of the resulting quotient in the non-empty cases.


title: Matrix rigidity and the Croot-Lev-Pach lemma

abstract: Matrix rigidity is a notion put forth by Valiant as a means for proving arithmetic circuit lower bounds. A matrix is rigid if it is far, in Hamming distance, from any low rank matrix. Despite decades of efforts, no explicit matrix rigid enough to carry out Valiant's plan has been found. Recently, Alman and Williams showed, contrary to common belief, that the $2^n \times 2^n$ Hadamard matrix could not be used for Valiant's program as it is not sufficiently rigid. In this note we observe a similar `non rigidity' phenomena for any $q^n \times q^n$ matrix $M$ of the form $M(x,y) = f(x+y)$, where $f:F_q^n \to F_q$ is any function and $F_q$ is a fixed finite field of $q$ elements ($n$ goes to infinity). The theorem follows almost immediately from a recent lemma of Croot, Lev and Pach which is also the main ingredient in the recent solution of the cap-set problem.


title: Dynamics on networks. Case of Heterogeneous Opinion Status Model

abstract: Here we developed a new conceptual, stochastic Heterogeneous Opinion-Status model (HOpS model), which is adaptive network model. The HOpS model admits to identify the main attributes of dynamics on networks and to study analytically the relation between topological network properties and processes taking place on a network. Another key point of the HOpS model is the possibility to study network dynamics via the novel parameter of heterogeneity. We show that not only clear topological network properties, such as node degree, but also, the nodes' status distribution (the factor of network heterogeneity) play an important role in so-called opinion spreading and information diffusion on a network. This model can be potentially used for studying the co-evolution of globally aggregated or averaged key observables of the earth system. These include natural variables such as atmospheric, oceanic and land carbon stocks, as well as socio-economic quantities such as global human population, economic production or wellbeing.


title: 3D-PRNN: Generating Shape Primitives with Recurrent Neural Networks

abstract: The success of various applications including robotics, digital content creation, and visualization demand a structured and abstract representation of the 3D world from limited sensor data. Inspired by the nature of human perception of 3D shapes as a collection of simple parts, we explore such an abstract shape representation based on primitives. Given a single depth image of an object, we present 3D-PRNN, a generative recurrent neural network that synthesizes multiple plausible shapes composed of a set of primitives. Our generative model encodes symmetry characteristics of common man-made objects, preserves long-range structural coherence, and describes objects of varying complexity with a compact representation. We also propose a method based on Gaussian Fields to generate a large scale dataset of primitive-based shape representations to train our network. We evaluate our approach on a wide range of examples and show that it outperforms nearest-neighbor based shape retrieval methods and is on-par with voxel-based generative models while using a significantly reduced parameter space.


title: Boundary value problems in dimensions seven, four and three related to exceptional holonomy

abstract: In the paper we introduce a boundary value problem for a G_{2} structure on a 7-manifold with boundary, with prescribed 3-form on the boundary. We make some general observations about this problem and then study in more detail reductions to 4 and 3 dimensions by imposing symmetry. The 3-dimensional reduction admits a dual variational formulation in terms of solutions of the real Monge-Ampere equation


title: BDCI: Behavioral Driven Conflict Identification

abstract: Source Code Management (SCM) systems support software evolution by providing features, such as version control, branching, and conflict detection. Despite the presence of these features, support to parallel software development is often limited. SCM systems can only address a subset of the conflicts that might be introduced by developers when concurrently working on multiple parallel branches. In fact, SCM systems can detect textual conflicts, which are generated by the concurrent modification of the same program locations, but they are unable to detect higher-order conflicts, which are generated by the concurrent modification of different program locations that generate program misbehaviors once merged. Higher-order conflicts are painful to detect and expensive to fix because they might be originated by the interference of apparently unrelated changes. In this paper we present Behavioral Driven Conflict Identification (BDCI), a novel approach to conflict detection. BDCI moves the analysis of conflicts from the source code level to the level of program behavior by generating and comparing behavioral models. The analysis based on behavioral models can reveal interfering changes as soon as they are introduced in the SCM system, even if they do not introduce any textual conflict. To evaluate the effectiveness and the cost of the proposed approach, we developed BDCIf , a specific instance of BDCI dedicated to the detection of higher-order conflicts related to the functional behavior of a program. The evidence collected by analyzing multiple versions of Git and Redis suggests that BDCIf can effectively detect higher-order conflicts and report how changes might interfere.


title: Quarter-Flux Hofstadter Lattice in Qubit-Compatible Microwave Cavity Array

abstract: Topological- and strongly-correlated- materials are exciting frontiers in condensed matter physics, married prominently in studies of the fractional quantum hall effect [1]. There is an active effort to develop synthetic materials where the microscopic dynamics and ordering arising from the interplay of topology and interaction may be directly explored. In this work we demonstrate a novel architecture for exploration of topological matter constructed from tunnel-coupled, time-reversalbroken microwave cavities that are both low loss and compatible with Josephson junction-mediated interactions [2]. Following our proposed protocol [3] we implement a square lattice Hofstadter model at a quarter flux per plaquette ({\alpha} = 1/4), with time-reversal symmetry broken through the chiral Wannier-orbital of resonators coupled to Yttrium-Iron-Garnet spheres. We demonstrate site-resolved spectroscopy of the lattice, time-resolved dynamics of its edge channels, and a direct measurement of the dispersion of the edge channels. Finally, we demonstrate the flexibility of the approach by erecting a tunnel barrier investigating dynamics across it. With the introduction of Josephson-junctions to mediate interactions between photons, this platform is poised to explore strongly correlated topological quantum science for the first time in a synthetic system.


title: A NuSTAR observation of the low mass X-ray binary GX 349+2 throughout the Z-track

abstract: Although the most luminous class of neutron star low mass X-ray binaries, known as Z sources, have been well studied, their behavior is not fully understood. In particular, what causes these sources to trace out the characteristic Z-shaped pattern on color-color or hardness-intensity diagrams is not well known. By studying the physical properties of the different spectral states of these sources, we may better understand such variability. With that goal in mind, we present a recent NuSTAR observation of the Z source GX 349+2, which spans approximately 2 days, and covers all its spectral states. By creating a hardness-intensity diagram we were able to extract five spectra and trace the change in spectral parameters throughout the Z-track. GX 349+2 shows a strong, broad Fe K$\alpha$ line in all states, regardless of the continuum model used. Through modeling of the reflection spectrum and Fe K$\alpha$ line we find that in most states the inner disk radius is consistent with remaining unchanged at an average radius of 17.5 $R_g$ or 36.4 km for a canonical 1.4 $M_\odot$ neutron star. During the brightest flaring branch, however, the inner disk radius from reflection is not well constrained.


title: Strong-pinning regimes by spherical inclusions in anisotropic type-II superconductors

abstract: The current-carrying capacity of type-II superconductors is decisively determined by how well material defect structures can immobilize vortex lines. In order to gain deeper insights into the fundamental pinning mechanisms, we have explored the case of vortex trapping by randomly distributed spherical inclusions using large-scale simulations of the time-dependent Ginzburg-Landau equations. We find that for a small density of particles having diameters of two coherence lengths, the vortex lattice preserves its structure and the critical current $j_c$ decays with the magnetic field following a power-law $B^{-\alpha}$ with $\alpha \approx 0.66$, which is consistent with predictions of strong-pinning theory. For a higher density of particles and/or larger inclusions, the lattice becomes progressively more disordered and the exponent smoothly decreases down to $\alpha \approx 0.3$. At high magnetic fields, all inclusions capture a vortex and the critical current decays faster than $B^{-1}$ as would be expected by theory. In the case of larger inclusions with a diameter of four coherence length, the magnetic-field dependence of the critical current is strongly affected by the ability of inclusions to capture multiple vortex lines. We found that at small densities, the fraction of inclusions trapping two vortex lines rapidly grows within narrow field range leading to a peak in $j_c(B)$-dependence within this range. With increasing inclusion density, this peak transforms into a plateau, which then smooths out. Using the insights gained from simulations, we determine the limits of applicability of strong-pinning theory and provide different routes to describe vortex pinning beyond those bounds.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.mtrl-sci physics.app-ph physics.comp-ph


title: Better Together: Joint Reasoning for Non-rigid 3D Reconstruction with Specularities and Shading

abstract: We demonstrate the use of shape-from-shading (SfS) to improve both the quality and the robustness of 3D reconstruction of dynamic objects captured by a single camera. Unlike previous approaches that made use of SfS as a post-processing step, we offer a principled integrated approach that solves dynamic object tracking and reconstruction and SfS as a single unified cost function. Moving beyond Lambertian S f S , we propose a general approach that models both specularities and shading while simultaneously tracking and reconstructing general dynamic objects. Solving these problems jointly prevents the kinds of tracking failures which can not be recovered from by pipeline approaches. We show state-of-the-art results both qualitatively and quantitatively.


title: Possible Beyond the Standard Model Physics Motivated by Muonic Puzzles

abstract: Recent measurements of the proton radius using the Lamb shift in muonic hydrogen are troublingly discrepant with values extracted from hydrogen spectroscopy and electron-proton scattering experiments. This discrepancy, which differs by more than five standard deviations, may be a signal of new physics caused by a violation of lepton universality. Another candidate for a new physics signal is the muon anomalous magnetic moment. The measurement at BNL differs from the standard model prediction by at least three standard deviations. Motivated by these two puzzles, first we use polarized lepton-nucleon elastic scattering to search for a new scalar boson, and furthermore we suggest new measurements of the nucleon form factors. Next, we display a method to analyze the beam dump experiments without using approximation on phase space, and we use it to constrain all possible new spin-0 and spin-1 particles, and a variety of other measurements to study the possibility of the new physics. Finally, assuming a new scalar boson can solve the two puzzles simultaneously, we present a general model-independent analysis and constrain the existence of the new physics.


title: On a Vizing-type integer domination conjecture

abstract: Given a simple graph $G$, a dominating set in $G$ is a set of vertices $S$ such that every vertex not in $S$ has a neighbor in $S$. Denote the domination number, which is the size of any minimum dominating set of $G$, by $\gamma(G)$. For any integer $k\ge 1$, a function $f : V (G) \rightarrow \{0, 1, . . ., k\}$ is called a \emph{$\{k\}$-dominating function} if the sum of its function values over any closed neighborhood is at least $k$. The weight of a $\{k\}$-dominating function is the sum of its values over all the vertices. The $\{k\}$-domination number of $G$, $\gamma_{\{k\}}(G)$, is defined to be the minimum weight taken over all $\{k\}$-domination functions. Bre\v{s}ar, Henning, and Klav\v{z}ar (On integer domination in graphs and Vizing-like problems. \emph{Taiwanese J. Math.} {10(5)} (2006) pp. 1317--1328) asked whether there exists an integer $k\ge 2$ so that $\gamma_{\{k\}}(G\square H)\ge \gamma(G)\gamma(H)$. In this note we use the Roman $\{2\}$-domination number, $\gamma_{R2}$ of Chellali, Haynes, Hedetniemi, and McRae, (Roman $\{2\}$-domination. \emph{Discrete Applied Mathematics} {204} (2016) pp. 22-28.) to prove that if $G$ is a claw-free graph and $H$ is an arbitrary graph, then $\gamma_{\{2\}}(G\square H)\ge \gamma_{R2}(G\square H)\ge \gamma(G)\gamma(H)$, which also implies the conjecture for all $k\ge 2$.


title: A Simple PTAS for the Dual Bin Packing Problem and Advice Complexity of Its Online Version

abstract: Recently, Renault (2016) studied the dual bin packing problem in the per-request advice model of online algorithms. He showed that given $O(1/\epsilon)$ advice bits for each input item allows approximating the dual bin packing problem online to within a factor of $1+\epsilon$. Renault asked about the advice complexity of dual bin packing in the tape-advice model of online algorithms. We make progress on this question. Let $s$ be the maximum bit size of an input item weight. We present a conceptually simple online algorithm that with total advice $O\left(\frac{s + \log n}{\epsilon^2}\right)$ approximates the dual bin packing to within a $1+\epsilon$ factor. To this end, we describe and analyze a simple offline PTAS for the dual bin packing problem. Although a PTAS for a more general problem was known prior to our work (Kellerer 1999, Chekuri and Khanna 2006), our PTAS is arguably simpler to state and analyze. As a result, we could easily adapt our PTAS to obtain the advice-complexity result. We also consider whether the dependence on $s$ is necessary in our algorithm. We show that if $s$ is unrestricted then for small enough $\epsilon > 0$ obtaining a $1+\epsilon$ approximation to the dual bin packing requires $\Omega_\epsilon(n)$ bits of advice. To establish this lower bound we analyze an online reduction that preserves the advice complexity and approximation ratio from the binary separation problem due to Boyar et al. (2016). We define two natural advice complexity classes that capture the distinction similar to the Turing machine world distinction between pseudo polynomial time algorithms and polynomial time algorithms. Our results on the dual bin packing problem imply the separation of the two classes in the advice complexity world.


title: Exploring Features for Predicting Policy Citations

abstract: In this study we performed an initial investigation and evaluation of altmetrics and their relationship with public policy citation of research papers. We examined methods for using altmetrics and other data to predict whether a research paper is cited in public policy and applied receiver operating characteristic curve on various feature groups in order to evaluate their potential usefulness. From the methods we tested, classifying based on tweet count provided the best results, achieving an area under the ROC curve of 0.91.


title: HTM-MAT: An online prediction software toolbox based on cortical machine learning algorithm

abstract: HTM-MAT is a MATLAB based toolbox for implementing cortical learning algorithms (CLA) including related cortical-like algorithms that possesses spatiotemporal properties. CLA is a suite of predictive machine learning algorithms developed by Numenta Inc. and is based on the hierarchical temporal memory (HTM). This paper presents an implementation of HTM-MAT with several illustrative examples including several toy datasets and compared with two sequence learning applications employing state-of-the-art algorithms - the recurrentjs based on the Long Short-Term Memory (LSTM) algorithm and OS-ELM which is based on an online sequential version of the Extreme Learning Machine. The performance of HTM-MAT using two historical benchmark datasets and one real world dataset is also compared with one of the existing sequence learning applications, the OS-ELM. The results indicate that HTM-MAT predictions are indeed competitive and can outperform OS-ELM in sequential prediction tasks.


title: Towards N-mode parametric electromechanical resonances

abstract: The ubiquity of parametric resonance is continually evident in the repeated experimental observations of this phenomenon in multiple physical systems. The elementary case of 2 mode parametric resonance of order 1 involves the excitation of a spectral tone of a parametrically driven mode at a sub-harmonic frequency of the higher directly driven mode. Historically, such examples of parametric resonance have been predominantly researched in a system of micro- and nanoelectromechanical resonators. Here, in this paper, we break this convention by showcasing a collection of experimental signatures in support of the concept of "N-mode parametric resonance" using a number of elementary microelectromechanical devices. Specifically, we present observations of 2, 3, (2+3) and (3+3) mode parametric resonances demonstrating co-existence of different regimes within the same device. In addition, we also present observations of intrinsic "Four-Wave Mixing" of parametric excitations. This paper presents contributions towards the existence proof for such multimode parametric resonances which can also be exploited for engineering benefit within the field of "micro and nanoelectromechanical resonators". The experimental results further point towards the possibility of the ultimate observation of N-mode parametric resonance in such physical system.


title: Homogenization of Maxwell's equations in layered system beyond static approximation

abstract: The propagation of electromagnetic waves through disordered layered system is considered in the paradigm of Maxwell's equations homogenization. In spite of the impossibility to describe the system in terms of effective dielectric permittivity and/or magnetic permeability the unified way to describe the propagation and Anderson localization of electromagnetic waves is proposed in terms of the introduced effective wave vector (effective refractive index). It is demonstrated that both real and imaginary parts of the effective wave vector (contrary to effective dielectric permittivity and/or magnetic permeability) are the self-averaging quantities at any frequency. The introduced effective wave vector is analytical function of frequency; corresponding Kramers-Kronig relation generalizes the Jones-Herbert-Thouless formula.


title: Proliferation of Twinning in HCP Metals: Application to Magnesium

abstract: Plastic deformation of metallic alloys usually takes place through slip, but occasionally involves twinning. In particular, twinning is important in hexagonal close packed materials where the easy slip systems are insufficient to accommodate arbitrary deformations. While deformation by slip mechanisms is reasonably well understood, less remains known about deformation by twinning. Indeed, the identification of relevant twinning modes remains an art. In this paper, we develop an universal framework combining fundamental kinematic definition of twins with large scale atomistic calculations to predict twinning modes of crystalline materials. We apply this framework to magnesium where there are two accepted twin modes -- tension and compression, but a number of anomalous observations. Surprisingly, our framework shows that there are a very large number of twinning modes that are important in the deformation process of magnesium consistent with the anomalous observations. Thus, in contrast to the traditional view where plastic deformation is kinematically partitioned between a few modes, our result argues that the physics of deformation in HCP materials is governed by an energetic and kinetic competition between numerous possibilities. Consequently, our findings suggest that the commonly used models of deformation physics need to be revisited in order to take into account a broader and richer variety of twin modes, and potentially points to new avenues of improving the mechanical properties.


title: Accelerated Image Reconstruction for Nonlinear Diffractive Imaging

abstract: The problem of reconstructing an object from the measurements of the light it scatters is common in numerous imaging applications. While the most popular formulations of the problem are based on linearizing the object-light relationship, there is an increased interest in considering nonlinear formulations that can account for multiple light scattering. In this paper, we propose an image reconstruction method, called CISOR, for nonlinear diffractive imaging, based on a nonconvex optimization formulation with total variation (TV) regularization. The nonconvex solver used in CISOR is our new variant of fast iterative shrinkage/thresholding algorithm (FISTA). We provide fast and memory-efficient implementation of the new FISTA variant and prove that it reliably converges for our nonconvex optimization problem. In addition, we systematically compare our method with other state-of-the-art methods on simulated as well as experimentally measured data in both 2D and 3D settings.


title: Waveform and Spectrum Management for Unmanned Aerial Systems Beyond 2025

abstract: The application domains of civilian unmanned aerial systems (UASs) include agriculture, exploration, transportation, and entertainment. The expected growth of the UAS industry brings along new challenges: Unmanned aerial vehicle (UAV) flight control signaling requires low throughput, but extremely high reliability, whereas the data rate for payload data can be significant. This paper develops UAV number projections and concludes that small and micro UAVs will dominate the US airspace with accelerated growth between 2028 and 2032. We analyze the orthogonal frequency division multiplexing (OFDM) waveform because it can provide the much needed flexibility, spectral efficiency, and, potentially, reliability and derive suitable OFDM waveform parameters as a function of UAV flight characteristics. OFDM also lends itself to agile spectrum access. Based on our UAV growth predictions, we conclude that dynamic spectrum access is needed and discuss the applicability of spectrum sharing techniques for future UAS communications.


title: A Two Factor Forward Curve Model with Stochastic Volatility for Commodity Prices

abstract: We describe a model for evolving commodity forward prices that incorporates three important dynamics which appear in many commodity markets: mean reversion in spot prices and the resulting Samuelson effect on volatility term structure, decorrelation of moves in different points on the forward curve, and implied volatility skew and smile. This model is a "forward curve model" - it describes the stochastic evolution of forward prices - rather than a "spot model" that models the evolution of the spot commodity price. Two Brownian motions drive moves across the forward curve, with a third Heston-like stochastic volatility process scaling instantaneous volatilities of all forward prices. In addition to an efficient numerical scheme for calculating European vanilla and early-exercise option prices, we describe an algorithm for Monte Carlo-based pricing of more generic derivative payoffs which involves an efficient approximation for the risk neutral drift that avoids having to simulate drifts for every forward settlement date required for pricing.


title: An Effective Training Method For Deep Convolutional Neural Network

abstract: In this paper, we propose the nonlinearity generation method to speed up and stabilize the training of deep convolutional neural networks. The proposed method modifies a family of activation functions as nonlinearity generators (NGs). NGs make the activation functions linear symmetric for their inputs to lower model capacity, and automatically introduce nonlinearity to enhance the capacity of the model during training. The proposed method can be considered an unusual form of regularization: the model parameters are obtained by training a relatively low-capacity model, that is relatively easy to optimize at the beginning, with only a few iterations, and these parameters are reused for the initialization of a higher-capacity model. We derive the upper and lower bounds of variance of the weight variation, and show that the initial symmetric structure of NGs helps stabilize training. We evaluate the proposed method on different frameworks of convolutional neural networks over two object recognition benchmark tasks (CIFAR-10 and CIFAR-100). Experimental results showed that the proposed method allows us to (1) speed up the convergence of training, (2) allow for less careful weight initialization, (3) improve or at least maintain the performance of the model at negligible extra computational cost, and (4) easily train a very deep model.


title: Boundary harmonic coordinates on manifolds with boundary in low regularity

abstract: In this paper, we prove the existence of $H^2$-regular coordinates on Riemannian $3$-manifolds with boundary, assuming only $L^2$-bounds on the Ricci curvature, $L^4$-bounds on the second fundamental form of the boundary, and a positive lower bound on the volume radius. The proof follows by extending the theory of Cheeger-Gromov convergence to include manifolds with boundary in the above low regularity setting. The main tools are boundary harmonic coordinates together with elliptic estimates and a geometric trace estimate, and a rigidity argument using manifold doubling. Assuming higher regularity of the Ricci curvature, we also prove corresponding higher regularity estimates for the coordinates.


title: Rigid-layer Raman-active modes in $N$-layer Transition Metal Dichalcogenides: interlayer force constants and hyperspectral Raman imaging

abstract: We report a comparative study of rigid layer Raman-active modes in $N$-layer transition metal dichalcogenides. Trigonal prismatic (2Hc, such as MoSe$_2$, MoTe$_2$, WS$_2$, WSe$_2$) and distorted octahedral (1T', such as ReS$_2$ and ReSe$_2$) phases are considered. The Raman-active in-plane interlayer shear modes and out-of-plane interlayer breathing modes appear as well-defined features with wavenumbers in the range 0-40~cm$^{-1}$. These rigid layer modes are well-described by an elementary linear chain model from which the interlayer force constants are readily extracted. Remarkably, these force constants are all found to be of the same order of magnitude. Finally, we show that the prominent interlayer shear and breathing mode features allow high-precision hyperspectral Raman imaging of $N-$layer domains within a given transition metal dichalcogenide flake.


title: Are fibres in molecular cloud filaments real objects?

abstract: We analyse the morphology and kinematics of dense filamentary structures produced in a numerical simulation of a star--forming cloud of $1.4 \times 10^4 \, \Msun$ evolving under their own self--gravity in a magnetized medium. This study is motivated by recent observations of velocity--coherent substructures ("fibres") in star-forming filaments. We find such "fibres" ubiquitously in our simulated filament. We found that a fibre in one projection is not necessarily a fibre in another projection, and thus, caution should be taken into account when considering them as real objects. We found that only the densest parts of the filament ($\sim$30\% of the densest volume, which contains $\sim$70\% of the mass) belong to fibres in 2 projections. Moreover, it is quite common that they are formed by separated density enhancements superposed along the line of sight. Observations of fibres can yield insight into the level of turbulent substructure driven by gravity, but care should be taken in interpreting the results given the problem of line of sight superposition. We also studied the morphology and kinematics of the 3D skeleton (spine), finding that subfilaments accrete structured material mainly along the magnetic field lines, which are preferentially perpendicular to the skeleton. The magnetic field is at the same time dragged by the velocity field due to the gravitational collapse.


title: Intrinsic3D: High-Quality 3D Reconstruction by Joint Appearance and Geometry Optimization with Spatially-Varying Lighting

abstract: We introduce a novel method to obtain high-quality 3D reconstructions from consumer RGB-D sensors. Our core idea is to simultaneously optimize for geometry encoded in a signed distance field (SDF), textures from automatically-selected keyframes, and their camera poses along with material and scene lighting. To this end, we propose a joint surface reconstruction approach that is based on Shape-from-Shading (SfS) techniques and utilizes the estimation of spatially-varying spherical harmonics (SVSH) from subvolumes of the reconstructed scene. Through extensive examples and evaluations, we demonstrate that our method dramatically increases the level of detail in the reconstructed scene geometry and contributes highly to consistent surface texture recovery.


title: Switching probability of all-perpendicular spin valve nanopillars

abstract: In all-perpendicular spin valve nanopillars the probability density of the free-layer magnetization is independent of the azimuthal angle and its evolution equation simplifies considerably compared to the general, nonaxisymmetric geometry. Expansion of the time-dependent probability density to Legendre polynomials enables analytical integration of the evolution equation and yields a compact expression for the practically relevant switching probability. This approach is valid when the free layer behaves as a single-domain magnetic particle and it can be readily applied to fitting experimental data.


title: On the expected number of internal equilibria in random evolutionary games with correlated payoff matrix

abstract: The analysis of equilibrium points in random games has been of great interest in evolutionary game theory, with important implications for understanding of complexity in a dynamical system, such as its behavioural, cultural or biological diversity. The analysis so far has focused on random games of independent payoff entries. In this paper, we overcome this restrictive assumption by considering multi-player two-strategy evolutionary games where the payoff matrix entries are correlated random variables. Using techniques from the random polynomial theory we establish a closed formula for the mean numbers of internal (stable) equilibria. We then characterise the asymptotic behaviour of this important quantity for large group sizes and study the effect of the correlation. Our results show that decreasing the correlation among payoffs (namely, of a strategist for different group compositions) leads to larger mean numbers of (stable) equilibrium points, suggesting that the system or population behavioural diversity can be promoted by increasing independence of the payoff entries. Numerical results are provided to support the obtained analytical results.

categories: math.AP math.DS math.PR q-bio.PE q-bio.QM


title: On Resource Pooling and Separation for LRU Caching

abstract: Caching systems using the Least Recently Used (LRU) principle have now become ubiquitous. A fundamental question for these systems is whether the cache space should be pooled together or divided to serve multiple flows of data item requests in order to minimize the miss probabilities. In this paper, we show that there is no straight yes or no answer to this question, depending on complex combinations of critical factors, including, e.g., request rates, overlapped data items across different request flows, data item popularities and their sizes. Specifically, we characterize the asymptotic miss probabilities for multiple competing request flows under resource pooling and separation for LRU caching when the cache size is large. Analytically, we show that it is asymptotically optimal to jointly serve multiple flows if their data item sizes and popularity distributions are similar and their arrival rates do not differ significantly; the self-organizing property of LRU caching automatically optimizes the resource allocation among them asymptotically. Otherwise, separating these flows could be better, e.g., when data sizes vary significantly. We also quantify critical points beyond which resource pooling is better than separation for each of the flows when the overlapped data items exceed certain levels. Technically, we generalize existing results on the asymptotic miss probability of LRU caching for a broad class of heavy-tailed distributions and extend them to multiple competing flows with varying data item sizes, which also validates the Che approximation under certain conditions. These results provide new insights on improving the performance of caching systems.


title: Stroboscopic qubit measurement with squeezed illumination

abstract: Microwave squeezing represents the ultimate sensitivity frontier for superconducting qubit measurement. However, observation of enhancement has remained elusive, in part because integration with conventional dispersive readout pollutes the signal channel with antisqueezed vacuum. Here we induce a stroboscopic light-matter coupling with superior squeezing compatibility, and observe an increase in the room-temperature signal-to-noise ratio of 24%. Squeezing the orthogonal phase controls measurement backaction, slowing dephasing by a factor of 1.8. This protocol enables the practical use of microwave squeezing for qubit state measurement.


title: Simple root flows for Hitchin representations

abstract: We study simple root flows and Liouville currents for Hitchin representations. We show that the Liouville current is associated to the measure of maximal entropy for a simple root flow, derive a Liouville volume rigidity result, and construct a Liouville pressure metric on the Hitchin component.


title: Query-guided Regression Network with Context Policy for Phrase Grounding

abstract: Given a textual description of an image, phrase grounding localizes objects in the image referred by query phrases in the description. State-of-the-art methods address the problem by ranking a set of proposals based on the relevance to each query, which are limited by the performance of independent proposal generation systems and ignore useful cues from context in the description. In this paper, we adopt a spatial regression method to break the performance limit, and introduce reinforcement learning techniques to further leverage semantic context information. We propose a novel Query-guided Regression network with Context policy (QRC Net) which jointly learns a Proposal Generation Network (PGN), a Query-guided Regression Network (QRN) and a Context Policy Network (CPN). Experiments show QRC Net provides a significant improvement in accuracy on two popular datasets: Flickr30K Entities and Referit Game, with 14.25% and 17.14% increase over the state-of-the-arts respectively.


title: A network approach to topic models

abstract: One of the main computational and scientific challenges in the modern age is to extract useful information from unstructured texts. Topic models are one popular machine-learning approach which infers the latent topical structure of a collection of documents. Despite their success --- in particular of its most widely used variant called Latent Dirichlet Allocation (LDA) --- and numerous applications in sociology, history, and linguistics, topic models are known to suffer from severe conceptual and practical problems, e.g. a lack of justification for the Bayesian priors, discrepancies with statistical properties of real texts, and the inability to properly choose the number of topics. Here we obtain a fresh view on the problem of identifying topical structures by relating it to the problem of finding communities in complex networks. This is achieved by representing text corpora as bipartite networks of documents and words. By adapting existing community-detection methods -- using a stochastic block model (SBM) with non-parametric priors -- we obtain a more versatile and principled framework for topic modeling (e.g., it automatically detects the number of topics and hierarchically clusters both the words and documents). The analysis of artificial and real corpora demonstrates that our SBM approach leads to better topic models than LDA in terms of statistical model selection. More importantly, our work shows how to formally relate methods from community detection and topic modeling, opening the possibility of cross-fertilization between these two fields.

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


title: On optimal periodic dividend strategies for L\'evy risk processes

abstract: In this paper, we revisit the optimal periodic dividend problem, in which dividend payments can only be made at the jump times of an independent Poisson process. In the dual (spectrally positive L\'evy) model, recent results have shown the optimality of a periodic barrier strategy, which pays dividends at Poissonian dividend-decision times, if and only if the surplus is above some level. In this paper, we show the optimality of this strategy for a spectrally negative L\'evy process whose dual has a completely monotone L\'evy density. The optimal strategies and value functions are concisely written in terms of the scale functions. Numerical results are also provided.


title: Scoped Extension Methods in Dynamically-Typed Languages

abstract: Context. An extension method is a method declared in a package other than the package of its host class. Thanks to extension methods, developers can adapt to their needs classes they do not own: adding methods to core classes is a typical use case. This is particularly useful for adapting software and therefore to increase reusability. Inquiry. In most dynamically-typed languages, extension methods are globally visible. Because any developer can define extension methods for any class, naming conflicts occur: if two developers define an extension method with the same signature in the same class, only one extension method is visible and overwrites the other. Similarly, if two developers each define an extension method with the same name in a class hierarchy, one overrides the other. To avoid such "accidental overrides", some dynamically-typed languages limit the visibility of an extension method to a particular scope. However, their semantics have not been fully described and compared. In addition, these solutions typically rely on a dedicated and slow method lookup algorithm to resolve conflicts at runtime. Approach. In this article, we present a formalization of the underlying models of Ruby refinements, Groovy categories, Classboxes, and Method Shelters that are scoping extension method solutions in dynamically-typed languages. Knowledge. Our formal framework allows us to objectively compare and analyze the shortcomings of the studied solutions and other different approaches such as MultiJava. In addition, language designers can use our formal framework to determine which mechanism has less risk of "accidental overrides". Grounding. Our comparison and analysis of existing solutions is grounded because it is based on denotational semantics formalizations. Importance. Extension methods are widely used in programming languages that support them, especially dynamically-typed languages such as Pharo, Ruby or Python. However, without a carefully designed mechanism, this feature can cause insidious hidden bugs or can be voluntarily used to gain access to protected operations, violate encapsulation or break fundamental invariants.


title: On the Effect of Semantically Enriched Context Models on Software Modularization

abstract: Many of the existing approaches for program comprehension rely on the linguistic information found in source code, such as identifier names and comments. Semantic clustering is one such technique for modularization of the system that relies on the informal semantics of the program, encoded in the vocabulary used in the source code. Treating the source code as a collection of tokens loses the semantic information embedded within the identifiers. We try to overcome this problem by introducing context models for source code identifiers to obtain a semantic kernel, which can be used for both deriving the topics that run through the system as well as their clustering. In the first model, we abstract an identifier to its type representation and build on this notion of context to construct contextual vector representation of the source code. The second notion of context is defined based on the flow of data between identifiers to represent a module as a dependency graph where the nodes correspond to identifiers and the edges represent the data dependencies between pairs of identifiers. We have applied our approach to 10 medium-sized open source Java projects, and show that by introducing contexts for identifiers, the quality of the modularization of the software systems is improved. Both of the context models give results that are superior to the plain vector representation of documents. In some cases, the authoritativeness of decompositions is improved by 67%. Furthermore, a more detailed evaluation of our approach on JEdit, an open source editor, demonstrates that inferred topics through performing topic analysis on the contextual representations are more meaningful compared to the plain representation of the documents. The proposed approach in introducing a context model for source code identifiers paves the way for building tools that support developers in program comprehension tasks such as application and domain concept location, software modularization and topic analysis.


title: Predicting the Law Area and Decisions of French Supreme Court Cases

abstract: In this paper, we investigate the application of text classification methods to predict the law area and the decision of cases judged by the French Supreme Court. We also investigate the influence of the time period in which a ruling was made over the textual form of the case description and the extent to which it is necessary to mask the judge's motivation for a ruling to emulate a real-world test scenario. We report results of 96% f1 score in predicting a case ruling, 90% f1 score in predicting the law area of a case, and 75.9% f1 score in estimating the time span when a ruling has been issued using a linear Support Vector Machine (SVM) classifier trained on lexical features.


title: Deep Metric Learning with Angular Loss

abstract: The modern image search system requires semantic understanding of image, and a key yet under-addressed problem is to learn a good metric for measuring the similarity between images. While deep metric learning has yielded impressive performance gains by extracting high level abstractions from image data, a proper objective loss function becomes the central issue to boost the performance. In this paper, we propose a novel angular loss, which takes angle relationship into account, for learning better similarity metric. Whereas previous metric learning methods focus on optimizing the similarity (contrastive loss) or relative similarity (triplet loss) of image pairs, our proposed method aims at constraining the angle at the negative point of triplet triangles. Several favorable properties are observed when compared with conventional methods. First, scale invariance is introduced, improving the robustness of objective against feature variance. Second, a third-order geometric constraint is inherently imposed, capturing additional local structure of triplet triangles than contrastive loss or triplet loss. Third, better convergence has been demonstrated by experiments on three publicly available datasets.


title: Tunable dimensional crossover and magnetocrystalline anisotropy in Fe$_2$P-based alloys

abstract: Electronic structure calculations are used to examine the magnetic properties of Fe$_2$P-based alloys and the mechanisms through which the Curie temperature and magnetocrystalline anisotropy can be optimized for specific applications. It is found that at elevated temperatures the magnetic interaction in pure Fe$_2$P develops a pronounced two-dimensional character due to the suppression of the magnetization in one of the sublattices, but the interlayer coupling is very sensitive to band filling and structural distortions. This feature suggests a natural explanation of the observed sharp enhancement of the Curie temperature by alloying with multiple elements, such as Co, Ni, Si, and B. The magnetocrystalline anisotropy is also tunable by electron doping, reaching a maximum near the electron count of pure Fe$_2$P. These findings enable the optimization of the alloy content, suggesting co-alloying of Fe$_2$P with Co (or Ni) and Si as a strategy for maximizing the magnetocrystalline anisotropy at and above room temperature.


title: Exactly Solvable Model for Two Dimensional Topological Superconductor

abstract: In this paper, we present an exactly solvable model for two dimensional topological superconductor with helical Majorana edge modes protected by time reversal symmetry. Our construction is based on the idea of decorated domain walls and makes use of the Kasteleyn orientation on a two dimensional lattice, which was used for the construction of the symmetry protected fermion phase with $Z_2$ symmetry in Ref. \onlinecite{Tarantino2016,Ware2016}. By decorating the time reversal domain walls with spinful Majorana chains, we are able to construct a commuting projector Hamiltonian with zero correlation length ground state wave function that realizes a strongly interacting version of the two dimensional topological superconductor. From our construction, it can be seen that the $T^2=-1$ transformation rule for the fermions is crucial for the existence of such a nontrivial phase; with $T^2=1$, our construction does not work.


title: Electron energy-loss spectroscopy of quasi-two-dimensional crystals: Beyond the energy-loss functions formalism

abstract: A consistent theory of electron energy-loss spectroscopy (EELS) includes two indispensable elements: (i) electronic response of the target system and (ii) quantum kinematics of probing electrons. While for the bulk materials and their surfaces separating these two aspects and focusing on the former is the usual satisfactory practice (the energy-loss functions formalism), we show that, for quasi-2D crystals, the interplay of the system's electronic response and the details of the probe's motion affects EEL spectra dramatically, and it must be taken into account for the reliable interpretation of the experiment. To this end, we come up with a unified theory which, on the same footing, treats both the long- and short-range scattering, within both the transmission and reflection experimental setups. Our calculations performed for graphene reveal a phenomenon of a strong coupling between the $\pi+\sigma$ plasmon excitation and elastic {\it scattering resonances}. Freed from the conventions of the energy-loss functions formalism, our theory serves as a consistent and systematic means of the understanding of the EELS of quasi-2D materials.


title: Exponentiated Generalized Pareto Distribution: Properties and applications towards Extreme Value Theory

abstract: The Generalized Pareto Distribution (GPD) plays a central role in modelling heavy tail phenomena in many applications. Applying the GPD to actual datasets however is a non-trivial task. One common way suggested in the literature to investigate the tail behaviour is to take logarithm to the original dataset in order to reduce the sample variability. Inspired by this, we propose and study the Exponentiated Generalized Pareto Distribution (exGPD), which is created via log-transform of the GPD variable. After introducing the exGPD we derive various distributional quantities, including the moment generating function, tail risk measures. As an application we also develop a plot as an alternative to the Hill plot to identify the tail index of heavy tailed datasets, based on the moment matching for the exGPD. Various numerical analyses with both simulated and actual datasets show that the proposed plot works well.


title: Magnetic field-dependent dynamics and field-driven metal-to-insulator transition of the half-filled Hubbard model: A DMFT+DMRG study

abstract: We study the magnetic field driven metal-to-insulator transition in half-filled Hubbard model on the Bethe lattice, using the dynamical mean-field theory by solving the quantum impurity problem with density-matrix renormalization group algorithm. The method enables us to obtain a high-resolution spectral densities in the presence of a magnetic field. It is found that the Kondo resonance at the Fermi level splits at relatively high magnetic field: the spin-up and spin-down components move away from the Fermi level and finally form a spin polarized band insulator. By calculating the magnetization and spin susceptibility, we clarify that an applied magnetic field drives a transition from a paramagnetic metallic phase to a band insulting phase. In the weak interaction regime, the nature of the transition is continuous and captured by the Stoner's description, while in the strong interaction regime the transition is very likely to be metamagnetic, evidenced by the hysteresis curve. Furthermore, we determine the phase boundary by tracking the kink in the magnetic susceptibility, and the step-like change of the entanglement entropy and the entanglement gap closing. Interestingly, the phase boundary determined from these two different ways are largely consistent with each other.


title: Abstract Hidden Markov Models: a monadic account of quantitative information flow

abstract: Hidden Markov Models, HMM's, are mathematical models of Markov processes with state that is hidden, but from which information can leak. They are typically represented as 3-way joint-probability distributions. We use HMM's as denotations of probabilistic hidden-state sequential programs: for that, we recast them as `abstract' HMM's, computations in the Giry monad $\mathbb{D}$, and we equip them with a partial order of increasing security. However to encode the monadic type with hiding over some state $\mathcal{X}$ we use $\mathbb{D}\mathcal{X}\to \mathbb{D}^2\mathcal{X}$ rather than the conventional $\mathcal{X}{\to}\mathbb{D}\mathcal{X}$ that suffices for Markov models whose state is not hidden. We illustrate the $\mathbb{D}\mathcal{X}\to \mathbb{D}^2\mathcal{X}$ construction with a small Haskell prototype. We then present uncertainty measures as a generalisation of the extant diversity of probabilistic entropies, with characteristic analytic properties for them, and show how the new entropies interact with the order of increasing security. Furthermore, we give a `backwards' uncertainty-transformer semantics for HMM's that is dual to the `forwards' abstract HMM's - it is an analogue of the duality between forwards, relational semantics and backwards, predicate-transformer semantics for imperative programs with demonic choice. Finally, we argue that, from this new denotational-semantic viewpoint, one can see that the Dalenius desideratum for statistical databases is actually an issue in compositionality. We propose a means for taking it into account.


title: New concept of connection in signed graphs

abstract: In a signed graph each edge has a sign, $+1$ or $-1$. We introduce in the present paper a new definition of connection in a signed graph by the existence of both positive and negative chains between vertices. We prove some results and properties of this definition, such as sign components, sign articulation vertices, and sign isthmi, and we compare them to corresponding graph and signed-graphic matroid properties. We apply our results to signed graphs without positive cycles. For signed graphs in which every edge is negative our properties become parity properties.