title: On-Line Estimation of Local Oscillator Noise and Optimisation of Servo Parameters in Atomic Clocks

abstract: For atomic frequency standards in which fluctuations of the local oscillator (LO) frequency are the dominant noise source, we examine the role of the servo algorithm that predicts and corrects these frequency fluctuations. We derive the optimal linear prediction algorithm, showing how to measure the relevant spectral properties of the noise and optimise servo parameters while the standard is running, using only the atomic error signal. We find that, for realistic LO noise spectra, a conventional integrating servo with a properly chosen gain performs nearly as well as the optimal linear predictor. Using simple analytical models and numerical simulations, we establish optimum probe times as a function of clock atom number and of the dominant noise type in the local oscillator. We calculate the resulting LO-dependent scaling of achievable clock stability with atom number for product states as well as for maximally-correlated states.


title: Minimal cut-generating functions are nearly extreme

abstract: We study continuous (strongly) minimal cut generating functions for the model where all variables are integer. We consider both the original Gomory-Johnson setting as well as a recent extension by Cornu\'ejols and Y{\i}ld{\i}z. We show that for any continuous minimal or strongly minimal cut generating function, there exists an extreme cut generating function that approximates the (strongly) minimal function as closely as desired. In other words, the extreme functions are "dense" in the set of continuous (strongly) minimal functions.


title: Imitating Driver Behavior with Generative Adversarial Networks

abstract: The ability to accurately predict and simulate human driving behavior is critical for the development of intelligent transportation systems. Traditional modeling methods have employed simple parametric models and behavioral cloning. This paper adopts a method for overcoming the problem of cascading errors inherent in prior approaches, resulting in realistic behavior that is robust to trajectory perturbations. We extend Generative Adversarial Imitation Learning to the training of recurrent policies, and we demonstrate that our model outperforms rule-based controllers and maximum likelihood models in realistic highway simulations. Our model both reproduces emergent behavior of human drivers, such as lane change rate, while maintaining realistic control over long time horizons.


title: Extreme functions with an arbitrary number of slopes

abstract: For the one dimensional infinite group relaxation, we construct a sequence of extreme valid functions that are piecewise linear and such that for every natural number $k\geq 2$, there is a function in the sequence with $k$ slopes. This settles an open question in this area regarding a universal bound on the number of slopes for extreme functions. The function which is the pointwise limit of this sequence is an extreme valid function that is continuous and has an infinite number of slopes. This provides a new and more refined counterexample to an old conjecture of Gomory and Johnson stating that all extreme functions are piecewise linear. These constructions are extended to obtain functions for the higher dimensional group problems via the sequential-merge operation of Dey and Richard.


title: Multi-dimensional Core-Collapse Supernova Simulations with Neutrino Transport

abstract: We present multi-dimensional core-collapse supernova simulations using the Isotropic Diffusion Source Approximation (IDSA) for the neutrino transport and a modified potential for general relativity in two different supernova codes: FLASH and ELEPHANT. Due to the complexity of the core-collapse supernova explosion mechanism, simulations require not only high-performance computers and the exploitation of GPUs, but also sophisticated approximations to capture the essential microphysics. We demonstrate that the IDSA is an elegant and efficient neutrino radiation transfer scheme, which is portable to multiple hydrodynamics codes and fast enough to investigate long-term evolutions in two and three dimensions. Simulations with a 40 solar mass progenitor are presented in both FLASH (1D and 2D) and ELEPHANT (3D) as an extreme test condition. It is found that the black hole formation time is delayed in multiple dimensions and we argue that the strong standing accretion shock instability before black hole formation will lead to strong gravitational waves.


title: Homogeneity and thermodynamic identities in geometrothermodynamics

abstract: We propose a classification of thermodynamic systems in terms of the homogeneity properties of their fundamental equations. Ordinary systems correspond to homogeneous functions and non-ordinary systems are given by generalized homogeneous functions. This affects the explicit form of the Gibbs-Duhem relation and Euler's identity. We show that these generalized relations can be implemented in the formalism of black hole geometrothermodynamics in order to completely fix the arbitrariness present in Legendre invariant metrics.


title: Placing the spotted T Tauri star LkCa 4 on an HR diagram

abstract: Ages and masses of young stars are often estimated by comparing their luminosities and effective temperatures to pre-main sequence stellar evolution tracks, but magnetic fields and starspots complicate both the observations and evolution. To understand their influence, we study the heavily-spotted weak-lined T-Tauri star LkCa 4 by searching for spectral signatures of radiation originating from the starspot or starspot groups. We introduce a new methodology for constraining both the starspot filling factor and the spot temperature by fitting two-temperature stellar atmosphere models constructed from Phoenix synthetic spectra to a high-resolution near-IR IGRINS spectrum. Clearly discernable spectral features arise from both a hot photospheric component $T_{\mathrm{hot}} \sim4100$ K and to a cool component $T_{\mathrm{cool}} \sim2700-3000$ K, which covers $\sim80\%$ of the visible surface. This mix of hot and cool emission is supported by analyses of the spectral energy distribution, rotational modulation of colors and of TiO band strengths, and features in low-resolution optical/near-IR spectroscopy. Although the revised effective temperature and luminosity make LkCa 4 appear much younger and lower mass than previous estimates from unspotted stellar evolution models, appropriate estimates will require the production and adoption of spotted evolutionary models. Biases from starspots likely afflict most fully convective young stars and contribute to uncertainties in ages and age spreads of open clusters. In some spectral regions starspots act as a featureless veiling continuum owing to high rotational broadening and heavy line-blanketing in cool star spectra. Some evidence is also found for an anti-correlation between the velocities of the warm and cool components.


title: The contribution of bulk Comptonization to the soft X-ray excess in AGN

abstract: Bulk velocities exceed thermal velocities for sufficiently radiation pressure dominated accretion flows. We model the contribution of bulk Comptonization to the soft X-ray excess in AGN. Bulk Comptonization is due to both turbulence and the background shear. We calculate spectra both taking into account and not taking into account bulk velocities using scaled data from radiation magnetohydrodynamic (MHD) shearing box simulations. We characterize our results with temperatures and optical depths to make contact with other warm Comptonization models of the soft excess. We chose our fiducial mass, $M = 2 \times 10^6 M_{\odot}$, and accretion rate, $L/L_{\rm Edd} = 2.5$, to correspond to those fit to the super-Eddington narrow line Seyfert 1 (NLS1) RE1034+396. The temperatures, optical depths, and Compton $y$ parameters we find broadly agree with those fit to RE1034+396. The effect of bulk Comptonization is to shift the Wien tail to higher energy and lower the gas temperature, broadening the spectrum. Observations of the soft excess in NLS1s can constrain the properties of disc turbulence if the bulk Comptonization contribution can be separated out from contributions from other physical effects, such as reflection and absorption.


title: Generalized Gray codes with prescribed ends of small dimensions

abstract: Given pairwise distinct vertices $\{\alpha_i , \beta_i\}^k_{i=1}$ of the $n$-dimensional hypercube $Q_n$ such that the distance of $\alpha_i$ and $\beta_i$ is odd, are there paths $P_i$ between $\alpha_i$ and $\beta_i$ such that $\{V (P_i)\}^k_{i=1}$ partitions $V(Q_n)$? A positive solution for every $n\ge1$ and $k=1$ is known as a Gray code of dimension $n$. In this paper we settle this problem for small values of $n$.


title: Detailed modeling of dust distribution in the disk of HD 142527

abstract: We investigate the dust distribution in the crescent disk around HD 142527 based on the continuum emission at $890 \mathrm{\ \mu m}$ obtained by ALMA Cycle 0. The map is divided into $18$ azimuthal sectors, and the radial intensity profile in each sector is reproduced with a 2D disk model. Our model takes account of scattering and inclination of the disk as well as the azimuthal dependence in intensity. When the dust is assumed to have the conventional composition and maximum size of $1\ \mathrm{mm}$, the northwestern region ($PA=329^{\circ}-29^{\circ}$) cannot be reproduced. This is because the model intensity gets insensitive to the increase in surface density due to heavy self-scattering, reaching its ceiling much lower than the observed intensity. The ceiling depends on the position angle. When the scattering opacity is reduced by a factor of $10$, the intensity distribution is reproduced successfully in all the sectors including those in the northwestern region. The best fit model parameters depend little on the scattering opacity in the southern region where the disk is optically thin. The contrast of dust surface density along $PA$ is derived to be about $40$, much smaller than the value for the cases of conventional opacities ($70-130$). These results strongly suggest that the albedo is lower than considered by some reasons at least in the northwestern region.


title: Chimeras and clusters in networks of hyperbolic chaotic oscillators

abstract: We show that chimera states, where differentiated subsets of synchronized and desynchronized dynamical elements coexist, can emerge in networks of hyperbolic chaotic oscillators subject to global interactions. As local dynamics we employ Lozi maps which possess hyperbolic chaotic attractors. We consider a globally coupled system of these maps and use two statistical quantities to describe its collective behavior: the average fraction of elements belonging to clusters and the average standard deviation of state variables. Chimera states, clusters, complete synchronization, and incoherence are thus characterized on the space of parameters of the system. We find that chimera states are related to the formation of clusters in the system. In addition, we show that chimera states arise for a sufficiently long range of interactions in nonlocally coupled networks of these maps. Our results reveal that, under some circumstances, hyperbolicity does not impede the formation of chimera states in networks of coupled chaotic systems, as it had been previously hypothesized.


title: Speech Map: A Statistical Multimodal Atlas of 4D Tongue Motion During Speech from Tagged and Cine MR Images

abstract: Quantitative measurement of functional and anatomical traits of 4D tongue motion in the course of speech or other lingual behaviors remains a major challenge in scientific research and clinical applications. Here, we introduce a statistical multimodal atlas of 4D tongue motion using healthy subjects, which enables a combined quantitative characterization of tongue motion in a reference anatomical configuration. This atlas framework, termed Speech Map, combines cine- and tagged-MRI in order to provide both the anatomic reference and motion information during speech. Our approach involves a series of steps including (1) construction of a common reference anatomical configuration from cine-MRI, (2) motion estimation from tagged-MRI, (3) transformation of the motion estimations to the reference anatomical configuration, and (4) computation of motion quantities such as Lagrangian strain. Using this framework, the anatomic configuration of the tongue appears motionless, while the motion fields and associated strain measurements change over the time course of speech. In addition, to form a succinct representation of the high-dimensional and complex motion fields, principal component analysis is carried out to characterize the central tendencies and variations of motion fields of our speech tasks. Our proposed method provides a platform to quantitatively and objectively explain the differences and variability of tongue motion by illuminating internal motion and strain that have so far been intractable. The findings are used to understand how tongue function for speech is limited by abnormal internal motion and strain in glossectomy patients.


title: A Complex Quaternion Model for Hyperbolic 3-Space

abstract: In 1900, Macfarlane proposed a hyperbolic variation on Hamilton's quaternions that closely resembles Minkowski spacetime. Viewing this in a modern context, we expand upon Macfarlane's idea and develop a model for real hyperbolic 3-space in which both points and isometries are expressed as complex quaternions, analogous to Hamilton's famous theorem on Euclidean rotations. We use this to give new computational tools for studying isometries of hyperbolic 2- and 3-space. We also give a generalization to other quaternion algebras.


title: Manifestly T-dual formulation of AdS space

abstract: We present a manifestly T-dual formulation of curved spaces such as an AdS space. For group manifolds related by the orthogonal vielbein fields the three form H=dB in the doubled space is universal at least locally. We construct an affine nondegenerate doubled bosonic AdS algebra to define the AdS space with the Ramond-Ramond flux. The non-zero commutator of the left and right momenta leads to that the left momentum is in an AdS space while the right momentum is in a dS space. Dimensional reduction constraints and the physical AdS algebra are shown to preserve all the doubled coordinates.


title: WebSelect: A Research Prototype for Optimizing Ad Exposures based on Network Structure

abstract: This paper describes a Research Prototype, WebSelect, designed to assist online media planners in deciding which websites to target for a media campaign. The salient feature of the prototype is its ability to capture and utilize the overlap information in website traffic and use it for media planning exercises. In addition, the prototype possesses the capability to include targeting parameters like users age group and income along with the varying advertising costs across different websites. The prototype uses a genetic algorithm at the backend to select the final subset to target from the possible website set.


title: Macfarlane Hyperbolic 3-Manfiolds

abstract: We identify and study a class of hyperbolic 3-manifolds (which we call Macfarlane manifolds) whose quaternion algebras admit a geometric interpretation analogous to Hamilton's classical model for Euclidean rotations. We characterize these manifolds arithmetically, and show that infinitely many commensurability classes of them arise in diverse topological and arithmetic settings. We then use this perspective to introduce a new method for computing their Dirichlet domains. We also give similar results for a class of hyperbolic surfaces and explore their occurrence as subsurfaces of Macfarlane manifolds.


title: Theoretical approach to the nanoporous phase diagram of carbon

abstract: Nanoporous carbon has been extensively used in a wide range of applications ranging from water treatment to electrochemical applications, such as in energy storage devices. An effort to relate structural to thermodynamical properties has not been explored from an atomistic approach. In this work we present numerical strategies to produce and study nanoporous carbon structures, using molecular dynamics simulations and a many-body potential. We designed a heating-quenching procedure in a thermodynamic region bounded by the critical, and triple point densities of carbon to study an ensemble of 1750 atomic arrangements produced at different densities, quench rates, and using graphite and diamond unit cells as precursor structures. All these samples were numerically characterized through the calculation of the free volumes, surface areas, radial distribution functions, and structure factors. We found particularly useful the potential energy dependence with $sp^3$ hybridization content, to determine structural phases through clustering methods. Three phases were related to graphite-like, sponge-like, and unstable states. We showed that our results are compatible with available experiments and different theoretical schemes, concluding that the use of Tersoff potential is a reliable choice to produce nanoporous structures with low computational cost.


title: Coherent backscattering of light from ultracold and optically dense atomic ensembles

abstract: We review experimental and theoretical studies of coherent backscattering of near resonant radiation from an ultracold atomic gas in the weak localization regime. Recent accomplishments in high resolution spectroscopy of atomic ensembles based on the coherent backscattering process are discussed. We also propose several new experimental schemes for time-dependent spectroscopy as applied to multiple scattering in the regime of weak localization.


title: A graph cut approach to 3D tree delineation, using integrated airborne LiDAR and hyperspectral imagery

abstract: Recognising individual trees within remotely sensed imagery has important applications in forest ecology and management. Several algorithms for tree delineation have been suggested, mostly based on locating local maxima or inverted basins in raster canopy height models (CHMs) derived from Light Detection And Ranging (LiDAR) data or photographs. However, these algorithms often lead to inaccurate estimates of forest stand characteristics due to the limited information content of raster CHMs. Here we develop a 3D tree delineation method which uses graph cut to delineate trees from the full 3D LiDAR point cloud, and also makes use of any optical imagery available (hyperspectral imagery in our case). First, conventional methods are used to locate local maxima in the CHM and generate an initial map of trees. Second, a graph is built from the LiDAR point cloud, fused with the hyperspectral data. For computational efficiency, the feature space of hyperspectral imagery is reduced using robust PCA. Third, a multi-class normalised cut is applied to the graph, using the initial map of trees to constrain the number of clusters and their locations. Finally, recursive normalised cut is used to subdivide, if necessary, each of the clusters identified by the initial analysis. We call this approach Multiclass Cut followed by Recursive Cut (MCRC). The effectiveness of MCRC was tested using three datasets: i) NewFor, ii) a coniferous forest in the Italian Alps, and iii) a deciduous woodland in the UK. The performance of MCRC was usually superior to that of other delineation methods, and was further improved by including high-resolution optical imagery. Since MCRC delineates the entire LiDAR point cloud in 3D, it allows individual crown characteristics to be measured. By making full use of the data available, graph cut has the potential to considerably improve the accuracy of tree delineation.


title: An Enhanced Threshold Resummation Formalism for Lepton Pair Production and Its Effects in the Determination of Parton Distribution Functions

abstract: We demonstrate that theoretical predictions using current resummation techniques for the lepton pair production (LPP) rapidity and $x_F$ distributions can be inconsistent with data in high rapidity and $x_F$ kinematic regions by observing their effect in global fits of parton distribution functions (PDFs). We present an enhanced resummation technique for the LPP rapidity and $x_F$ distributions that agrees with LPP data. The enhanced resummation method is used in conjunction with threshold resummation in deep inelastic scattering (DIS) to perform two global fits of PDFs using the minimal and Borel prescriptions. The results are analyzed to determine the effect of threshold resummation on global fits of PDFs.


title: Lower Bounds on the Complexity of Solving Two Classes of Non-cooperative Games

abstract: This paper studies the complexity of solving two classes of non-cooperative games in a distributed manner in which the players communicate with a set of system nodes over noisy communication channels. The complexity of solving each game class is defined as the minimum number of iterations required to find a Nash equilibrium (NE) of any game in that class with $\epsilon$ accuracy. First, we consider the class $\mathcal{G}$ of all $N$-player non-cooperative games with a continuous action space that admit at least one NE. Using information-theoretic inequalities, we derive a lower bound on the complexity of solving $\mathcal{G}$ that depends on the Kolmogorov $2\epsilon$-capacity of the constraint set and the total capacity of the communication channels. We also derive a lower bound on the complexity of solving games in $\mathcal{G}$ which depends on the volume and surface area of the constraint set. We next consider the class of all $N$-player non-cooperative games with at least one NE such that the players' utility functions satisfy a certain (differential) constraint. We derive lower bounds on the complexity of solving this game class under both Gaussian and non-Gaussian noise models. Our result in the non-Gaussian case is derived by establishing a connection between the Kullback-Leibler distance and Fisher information.


title: Leveraging Quality Prediction Models for Automatic Writing Feedback

abstract: User-generated, multi-paragraph writing is pervasive and important in many social media platforms (i.e. Amazon reviews, AirBnB host profiles, etc). Ensuring high-quality content is important. Unfortunately, content submitted by users is often not of high quality. Moreover, the characteristics that constitute high quality may even vary between domains in ways that users are unaware of. Automated writing feedback has the potential to immediately point out and suggest improvements during the writing process. Most approaches, however, focus on syntax/phrasing, which is only one characteristic of high-quality content. Existing research develops accurate quality prediction models. We propose combining these models with model explanation techniques to identify writing features that, if changed, will most improve the text quality. To this end, we develop a perturbation-based explanation method for a popular class of models called tree-ensembles. Furthermore, we use a weak-supervision technique to adapt this method to generate feedback for specific text segments in addition to feedback for the entire document. Our user study finds that the perturbation-based approach, when combined with segment-specific feedback, can help improve writing quality on Amazon (review helpfulness) and Airbnb (host profile trustworthiness) by > 14% (3X improvement over recent automated feedback techniques).


title: Fabrication of a centimeter-long cavity on a nanofiber for cavity QED

abstract: We report the fabrication of a 1.2 cm long cavity directly on a nanofiber using femtosecond laser ablation. The cavity modes with finesse value in the range 200-400 can still maintain the transmission between 40-60%, which can enable "strong-coupling" regime of cavity QED for a single atom trapped 200 nm away from the fiber surface. For such cavity modes, we estimate the one-pass intra-cavity transmission to be 99.53%. Other cavity modes, which can enable high cooperativity in the range 3-10, show transmission over 60-85% and are suitable for fiber-based single photon sources and quantum nonlinear optics in the "Purcell" regime.


title: Bayesian inference with Monte Carlo approximation: Measuring regional differentiation in ceramic and glass vessel assemblages in Republican Italy, ca. 200 BCE - 20 CE

abstract: Methods of measuring differentiation in archaeological assemblages have long been based on attribute-level analyses of assemblages. This paper considers a method of comparing assemblages as probability distributions via the Hellinger distance, as calculated through a Dirichlet-categorical model of inference using Monte Carlo methods of approximation. This method has application within practice-theory traditions of archaeology, an approach which seeks to measure and associate different factors that comprise the habitus of society. It is implemented here focusing on the question of regional food consumption habits in Republican Italy in the last two centuries BCE, toward informing a perspective on mass social change.


title: Photoionization Models for the Inner Gaseous Disks of Herbig Be Stars: Evidence Against Magnetospheric Accretion?

abstract: We investigate the physical properties of the inner gaseous disks of the three, hot, Herbig B2e stars, HD 76534, HD 114981 and HD 216629, by modelling CFHT-ESPaDOns spectra using non-LTE radiative transfer codes. We assume that the emission lines are produced in a circumstellar disk heated solely by the photospheric radiation from the central star in order to test if the optical and near-IR emission lines can be reproduced without invoking magnetospheric accretion. The inner gaseous disk density was assumed to follow a simple power-law in the equatorial plane, and we searched for models that could reproduce observed lines of HI (H$\alpha$ and H$\beta$), HeI, CaII and FeII. For the three stars, good matches were found for all emission line profiles individually; however, no density model based on a single power-law was able to reproduce all of the observed emission lines. Among the single power-law models, the one with the gas density varying as $\sim 10^{-10}\,(R_{*}/R)^3\, \rm g\,cm^{-3}$ in the equatorial plane of a 25 R$_{*}$ ($0.78$ AU) disk did the best overall job of representing the optical emission lines of the three stars. This model implies a mass for the H$\alpha$-emitting portion of the inner gaseous disk of $\sim 10^{-9} M_*$. We conclude that the optical emission line spectra of these HBe stars can be qualitatively reproduced by a $\approx\,1$ AU, geometrically thin, circumstellar disk of negligible mass compared to the central star in Keplerian rotation and radiative equilibrium.


title: Characterization of the strong divisibility property for generalized Fibonacci polynomials

abstract: It is known that the greatest common divisor of two Fibonacci numbers is again a Fibonacci number. It is called the strong divisibility property. However, this property does not hold for every second order sequence. In this paper we study the generalized Fibonacci polynomials and classify them in two types depending on their Binet formula. We give a complete characterization for those polynomials that satisfy the strong divisibility property. We also giveformulas to calculate the $\gcd$ of those polynomials that do not satisfy the strong divisibility property. Note. This paper is now published in INTEGERS.


title: A Stochastic Programming Approach for Electric Vehicle Charging Network Design

abstract: Advantages of electric vehicles (EV) include reduction of greenhouse gas and other emissions, energy security, and fuel economy. The societal benefits of large-scale adoption of EVs cannot be realized without adequate deployment of publicly accessible charging stations. We propose a two-stage stochastic programming model to determine the optimal network of charging stations for a community considering uncertainties in arrival and dwell time of vehicles, battery state of charge of arriving vehicles, walkable range and charging preferences of drivers, demand during weekdays and weekends, and rate of adoption of EVs within a community. We conducted studies using sample average approximation (SAA) method which asymptotically converges to an optimal solution for a two-stage stochastic problem, however it is computationally expensive for large-scale instances. Therefore, we developed a heuristic to produce near to optimal solutions quickly for our data instances. We conducted computational experiments using various publicly available data sources, and benefits of the solutions are evaluated both quantitatively and qualitatively for a given community.


title: Phantom solution in a non-linear Israel-Stewart theory

abstract: In this paper we present a phantom solution with a big rip singularity in a non-linear regime of the Israel-Stewart formalism. In this framework it is possible to extend this causal formalism in order to describe accelerated expansion, where assumption of near equilibrium is no longer valid. We assume a flat universe filled with a single viscous fluid ruled by a barotropic EoS, $p=\omega \rho$, which can represent a late time accelerated phase of the cosmic evolution. The solution allows to cross the phantom divide without evoking an exotic matter fluid and the effective EoS parameter is always lesser than $-1$ and constant in time.


title: A Contextual Bandit Approach for Stream-Based Active Learning

abstract: Contextual bandit algorithms -- a class of multi-armed bandit algorithms that exploit the contextual information -- have been shown to be effective in solving sequential decision making problems under uncertainty. A common assumption adopted in the literature is that the realized (ground truth) reward by taking the selected action is observed by the learner at no cost, which, however, is not realistic in many practical scenarios. When observing the ground truth reward is costly, a key challenge for the learner is how to judiciously acquire the ground truth by assessing the benefits and costs in order to balance learning efficiency and learning cost. From the information theoretic perspective, a perhaps even more interesting question is how much efficiency might be lost due to this cost. In this paper, we design a novel contextual bandit-based learning algorithm and endow it with the active learning capability. The key feature of our algorithm is that in addition to sending a query to an annotator for the ground truth, prior information about the ground truth learned by the learner is sent together, thereby reducing the query cost. We prove that by carefully choosing the algorithm parameters, the learning regret of the proposed algorithm achieves the same order as that of conventional contextual bandit algorithms in cost-free scenarios, implying that, surprisingly, cost due to acquiring the ground truth does not increase the learning regret in the long-run. Our analysis shows that prior information about the ground truth plays a critical role in improving the system performance in scenarios where active learning is necessary.


title: Instantaneous Decentralized Poker

abstract: We present efficient protocols for amortized secure multiparty computation with penalties and secure cash distribution, of which poker is a prime example. Our protocols have an initial phase where the parties interact with a cryptocurrency network, that then enables them to interact only among themselves over the course of playing many poker games in which money changes hands. The high efficiency of our protocols is achieved by harnessing the power of stateful contracts. Compared to the limited expressive power of Bitcoin scripts, stateful contracts enable richer forms of interaction between standard secure computation and a cryptocurrency. We formalize the stateful contract model and the security notions that our protocols accomplish, and provide proofs using the simulation paradigm. Moreover, we provide a reference implementation in Ethereum/Solidity for the stateful contracts that our protocols are based on. We also adopt our off-chain cash distribution protocols to the special case of stateful duplex micropayment channels, which are of independent interest. In comparison to Bitcoin based payment channels, our duplex channel implementation is more efficient and has additional features.


title: Regular approximations of spectra of singular discrete linear Hamiltonian systems with one singular endpoint?

abstract: This paper is concerned with regular approximations of spectra of singular discrete linear Hamiltonian systems with one singular endpoint. For any given self-adjoint subspace extension (SSE) of the corresponding minimal subspace, its spectrum can be approximated by eigenvalues of a sequence of induced regular SSEs, generated by the same difference expression on smaller finite intervals. It is shown that every SSE of the minimal subspace has a pure discrete spectrum, and the k-th eigenvalue of any given SSE is exactly the limit of the k-th eigenvalues of the induced regular SSEs; that is, spectral exactness holds, in the limit circle case. Furthermore, error estimates for the approximations of eigenvalues are given in this case. In addition, in the limit point and intermediate cases, spectral inclusive holds.


title: Shock formation for $2D$ quasilinear wave systems featuring multiple speeds: Blowup for the fastest wave, with non-trivial interactions up to the singularity

abstract: We prove a stable shock formation result for a large class of systems of quasilinear wave equations in two spatial dimensions. We give a precise description of the dynamics all the way up to the singularity. Our main theorem applies to systems of two wave equations featuring two distinct wave speeds and various quasilinear and semilinear nonlinearities, while the solutions under study are (non-symmetric) perturbations of simple outgoing plane symmetric waves. The two waves are allowed to interact all the way up to the singularity. Our approach is robust and could be used to prove shock formation results for other related systems with many unknowns and multiple speeds, in various solution regimes, and in higher spatial dimensions. However, a fundamental aspect of our framework is that it applies only to solutions in which the "fastest wave" forms a shock while the remaining solution variables do not. Our approach is based on an extended version of the geometric vectorfield method developed by D. Christodoulou in his study of shock formation for scalar wave equations as well as the framework developed in our recent joint work with J. Luk, in which we proved a shock formation result for a quasilinear wave-transport system featuring a single wave operator. A key new difficulty that we encounter is that the geometric vectorfields that we use to commute the equations are, by necessity, adapted to the wave operator of the (shock-forming) fast wave and therefore exhibit very poor commutation properties with the slow wave operator, much worse than their commutation properties with a transport operator. To overcome this difficulty, we rely on a first-order reformulation of the slow wave equation, which, though somewhat limiting in the precision it affords, allows us to avoid uncontrollable commutator terms.


title: Theory of nonlinear polaritonics: $\chi^{(2)}$ scattering on a $\beta$-SiC surface

abstract: In this article we provide a practical prescription to harness the rigorous microscopic, quantum level descriptions of light-matter systems provided by Hopfield diagonalisation for quantum description of nonlinear scattering. A general frame to describe the practically important second-order optical nonlinearities which underpin sum and difference frequency generation is developed for arbitrarily inhomogeneous dielectric environments. Specific attention is then focussed on planar systems with optical nonlinearity mediated by a polar dielectric $\beta$-SiC halfspace. In this system we calculate the rate of second harmonic generation and the result compared to recent experimental measurements. Furthermore the rate of difference frequency generation of subdiffraction surface phonon polaritons on the $\beta$-SiC halfspace by two plane waves is calculated. The developed theory is easily integrated with commercial finite element solvers, opening the way for calculation of second-order nonlinear scattering coefficients in complex geometries which lack analytical linear solutions.


title: Optimal interdependence enhances robustness of complex systems

abstract: While interdependent systems have usually been associated with increased fragility, we show that strengthening the interdependence between dynamical processes on different networks can make them more robust. By coupling the dynamics of networks that in isolation exhibit catastrophic collapse with extinction of nodal activity, we demonstrate system-wide persistence of activity for an optimal range of interdependence between the networks. This is related to the appearance of attractors of the global dynamics comprising disjoint sets ("islands") of stable activity.


title: Weak Adaptive Submodularity and Group-Based Active Diagnosis with Applications to State Estimation with Persistent Sensor Faults

abstract: In this paper, we consider adaptive decision-making problems for stochastic state estimation with partial observations. First, we introduce the concept of weak adaptive submodularity, a generalization of adaptive submodularity, which has found great success in solving challenging adaptive state estimation problems. Then, for the problem of active diagnosis, i.e., discrete state estimation via active sensing, we show that an adaptive greedy policy has a near-optimal performance guarantee when the reward function possesses this property. We further show that the reward function for group-based active diagnosis, which arises in applications such as medical diagnosis and state estimation with persistent sensor faults, is also weakly adaptive submodular. Finally, in experiments of state estimation for an aircraft electrical system with persistent sensor faults, we observe that an adaptive greedy policy performs equally well as an exhaustive search.


title: On a binary system of Prendiville: The cubic case

abstract: We prove sharp decoupling inequalities for a class of two dimensional non-degenerate surfaces in R^5, introduced by Prendiville. As a consequence, we obtain sharp bounds on the number of integer solutions of the Diophantine systems associated with these surfaces.


title: Two-Dimensional Source Coding by Means of Subblock Enumeration

abstract: A technique of lossless compression via substring enumeration (CSE) attains compression ratios as well as popular lossless compressors for one-dimensional (1D) sources. The CSE utilizes a probabilistic model built from the circular string of an input source for encoding the source.The CSE is applicable to two-dimensional (2D) sources such as images by dealing with a line of pixels of 2D source as a symbol of an extended alphabet. At the initial step of the CSE encoding process, we need to output the number of occurrences of all symbols of the extended alphabet, so that the time complexity increase exponentially when the size of source becomes large. To reduce the time complexity, we propose a new CSE which can encode a 2D source in block-by-block instead of line-by-line. The proposed CSE utilizes the flat torus of an input 2D source as a probabilistic model for encoding the source instead of the circular string of the source. Moreover, we analyze the limit of the average codeword length of the proposed CSE for general sources.


title: Remote Estimation of the Wiener Process over a Channel with Random Delay

abstract: In this paper, we consider a problem of sampling a Wiener process, with samples forwarded to a remote estimator via a channel that consists of a queue with random delay. The estimator reconstructs a real-time estimate of the signal from causally received samples. Motivated by recent research on age-of-information, we study the optimal sampling strategy that minimizes the mean square estimation error subject to a sampling frequency constraint. We prove that the optimal sampling strategy is a threshold policy, and find the optimal threshold. This threshold is determined by the sampling frequency constraint and how much the Wiener process varies during the channel delay. An interesting consequence is that even in the absence of the sampling frequency constraint, the optimal strategy is not zero-wait sampling in which a new sample is taken once the previous sample is delivered; rather, it is optimal to wait for a non-zero amount of time after the previous sample is delivered, and then take the next sample. Further, if the sampling times are independent of the observed Wiener process, the optimal sampling problem reduces to an age-of-information optimization problem that has been recently solved. Our comparisons show that the estimation error of the optimal sampling policy is much smaller than those of age-optimal sampling, zero-wait sampling, and classic uniform sampling.


title: Cache-Aided Heterogeneous Networks: Coverage and Delay Analysis

abstract: This paper characterizes the performance of a generic $K$-tier cache-aided heterogeneous network (CHN), in which the base stations (BSs) across tiers differ in terms of their spatial densities, transmission powers, pathloss exponents, activity probabilities conditioned on the serving link and placement caching strategies. We consider that each user connects to the BS which maximizes its average received power and at the same time caches its file of interest. Modeling the locations of the BSs across different tiers as independent homogeneous Poisson Point processes (HPPPs), we derive closed-form expressions for the coverage probability and local delay experienced by a typical user in receiving each requested file. We show that our results for coverage probability and delay are consistent with those previously obtained in the literature for a single tier system.


title: Area minimizing discs in locally non-compact metric spaces

abstract: We solve the classical problem of Plateau in every metric space which is $1$-complemented in an ultra-completion of itself. This includes all proper metric spaces as well as many locally non-compact metric spaces, in particular, all dual Banach spaces, some non-dual Banach spaces such as $L^1$, all Hadamard spaces, and many more. Our results generalize corresponding results of Lytchak and the second author from the setting of proper metric spaces to that of locally non-compact ones. We furthermore solve the Dirichlet problem in the same class of spaces. The main new ingredient in our proofs is a suitable generalization of the Rellich-Kondrachov compactness theorem, from which we deduce a result about ultra-limits of sequences of Sobolev maps.


title: A model of mesons in finite extra-dimension

abstract: Recently,problem of stability of H-atom has been reported in extra-finite dimension,and found out that it is stable in extra-finite dimension of size,$R\leq\frac{a_0}{4}$,where,$a_0$ is the Bohr radius.Assuming that,the heavy flavoured mesons have also such stability controlled by the scale of coupling constant,we obtain corresponding QCD Bohr radius and it is found to be well within the present theoretical and experimental limit of higher dimension.We then study its consequences in their masses using effective string inspired potential model in higher dimension pursued by us.Within the uncertainty of masses of known Heavy Flavoured mesons the allowed range of extra dimension is $L\leq10^{-16}m$,which is well below the present theoretical and experimental limit,and far above the Planck length $\simeq1.5\times10^{-35}$ m.


title: Koszul cycles and Golod rings

abstract: Let $S$ be the power series ring or the polynomial ring over a field $K$ in the variables $x_1,\ldots,x_n$, and let $R=S/I$, where $I$ is proper ideal which we assume to be graded if $S$ is the polynomial ring. We give an explicit description of the cycles of the Koszul complex whose homology classes generate the Koszul homology of $R=S/I$ with respect to $x_1,\ldots,x_n$. The description is given in terms of the data of the free $S$-resolution of $R$. The result is used to determine classes of Golod ideals, among them proper ordinary powers and proper symbolic powers of monomial ideals. Our theory is also applied to stretched local rings.


title: Partial Bridging of Vaccine Efficacy to New Populations

abstract: Suppose one has data from one or more completed vaccine efficacy trials and wishes to estimate the efficacy in a new setting. Often logistical or ethical considerations make running another efficacy trial impossible. Fortunately, if there is a biomarker that is the primary modifier of efficacy, then the biomarker-conditional efficacy may be identical in the completed trials and the new setting, or at least informative enough to meaningfully bound this quantity. Given a sample of this biomarker from the new population, we might hope we can bridge the results of the completed trials to estimate the vaccine efficacy in this new population. Unfortunately, even knowing the true conditional efficacy in the new population fails to identify the marginal efficacy due to the unknown conditional unvaccinated risk. We define a curve that partially identifies (lower bounds) the marginal efficacy in the new population as a function of the population's marginal unvaccinated risk, under the assumption that one can identify bounds on the conditional unvaccinated risk in the new population. Interpreting the curve only requires identifying plausible regions of the marginal unvaccinated risk in the new population. We present a nonparametric estimator of this curve and develop valid lower confidence bounds that concentrate at a parametric rate. We use vaccine terminology throughout, but the results apply to general binary interventions and bounded outcomes.


title: Gamma-ray observations of Tycho's SNR with VERITAS and Fermi

abstract: High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it is a young, type Ia SNR that is well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-LAT, there have been several theoretical models proposed to explain its broadband emission and high-energy morphology. We report on an update to the gamma-ray measurements of Tycho's SNR with 147 hours of VERITAS and 84 months of Fermi-LAT observations, which represents about a factor of two increase in exposure over previously published data. About half of the VERITAS data benefited from a camera upgrade, which has made it possible to extend the TeV measurements toward lower energies. The TeV spectral index measured by VERITAS is consistent with previous results, but the expanded energy range softens a straight power-law fit. At energies higher than 400 GeV, the power-law index is $2.92 \pm 0.42_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$. It is also softer than the spectral index in the GeV energy range, $2.14 \pm 0.09_{\mathrm{stat}} \pm 0.02_{\mathrm{sys}}$, measured by this study using Fermi--LAT data. The centroid position of the gamma-ray emission is coincident with the center of the remnant, as well as with the centroid measurement of Fermi--LAT above 1 GeV. The results are consistent with an SNR shell origin of the emission, as many models assume. The updated spectrum points to a lower maximum particle energy than has been suggested previously.


title: Variability-Aware Design for Energy Efficient Computational Artificial Intelligence Platform

abstract: Portable computing devices, which include tablets, smart phones and various types of wearable sensors, experienced a rapid development in recent years. One of the most critical limitations for these devices is the power consumption as they use batteries as the power supply. However, the bottleneck of the power saving schemes in both hardware design and software algorithm is the huge variability in power consumption. The variability is caused by a myriad of factors, including the manufacturing process, the ambient environment (temperature, humidity), the aging effects and etc. As the technology node scaled down to 28nm and even lower, the variability becomes more severe. As a result, a platform for variability characterization seems to be very necessary and helpful.


title: A study of halo and relic radio emission in merging clusters using the Murchison Widefield Array

abstract: We have studied radio haloes and relics in nine merging galaxy clusters using the Murchison Widefield Array (MWA). The images used for this study were obtained from the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey which was carried out at 5 frequencies, viz. 88, 118, 154, 188 and 215 MHz. We detect diffuse radio emission in 8 of these clusters. We have estimated the spectra of haloes and relics in these clusters over the frequency range 80-1400 MHz; the first such attempt to estimate their spectra at low frequencies. The spectra follow a power law with a mean value of $\alpha = -1.13\pm0.21$ for haloes and $\alpha = -1.2\pm0.19$ for relics where, $S \propto \nu^{\alpha}$. We reclassify two of the cluster sources as radio galaxies. The low frequency spectra are thus an independent means of confirming the nature of cluster sources. Five of the nine clusters host radio haloes. For the remaining four clusters, we place upper limits on the radio powers of possible haloes in them. These upper limits are a factor of 2-20 below those expected from the $L_{\rm X}-P_{\rm 1.4}$ relation. These limits are the lowest ever obtained and the implications of these limits to the hadronic model of halo emission are discussed.


title: Reasoning about Probabilistic Defense Mechanisms against Remote Attacks

abstract: Despite numerous countermeasures proposed by practitioners and researchers, remote control-flow alteration of programs with memory-safety vulnerabilities continues to be a realistic threat. Guaranteeing that complex software is completely free of memory-safety vulnerabilities is extremely expensive. Probabilistic countermeasures that depend on random secret keys are interesting, because they are an inexpensive way to raise the bar for attackers who aim to exploit memory-safety vulnerabilities. Moreover, some countermeasures even support legacy systems. However, it is unclear how to quantify and compare the effectiveness of different probabilistic countermeasures or combinations of such countermeasures. In this paper we propose a methodology to rigorously derive security bounds for probabilistic countermeasures. We argue that by representing security notions in this setting as events in probabilistic games, similarly as done with cryptographic security definitions, concrete and asymptotic guarantees can be obtained against realistic attackers. These guarantees shed light on the effectiveness of single countermeasures and their composition and allow practitioners to more precisely gauge the risk of an attack.


title: Minimal control power of controlled dense coding and genuine tripartite entanglement

abstract: We investigate minimal control power (MCP) for controlled dense coding defined by the channel capacity. We obtain MCPs for extended three-qubit Greenberger-Horne-Zeilinger (GHZ) states and generalized three-qubit $W$ states. Among those GHZ states, the standard GHZ state is found to maximize the MCP and so does the standard $W$ state among the $W$-type states. We find the lower and upper bounds of the MCP and show for pure states that the lower bound, zero, is achieved if and only if the three-qubit state is biseparable or fully separable. The upper bound is achieved only for the standard GHZ state. Since the MCP is nonzero only when a three-qubit entanglement exists, this quantity may be a good candidate to measure the degree of genuine tripartite entanglement.


title: Proceedings of the 12th Workshop on User Interfaces for Theorem Provers

abstract: The UITP workshop series brings together researchers interested in designing, developing and evaluating user interfaces for automated reasoning tools, such as interactive proof assistants, automated theorem provers, model finders, tools for formal methods, and tools for visualising and manipulating logical formulas and proofs. The twelth edition of UITP took place in Coimbra, Portugal, and was part of the International Joint Conference on Automated Reasoning (IJCAR'16). The workshop consisted of an invited talk, six presentations of submitted papers and lively hands-on session for reasoning tools and their user-interface. These post-proceedings contain four contributed papers accepted for publication after a second round of reviewing after the workshop as well as the invited paper.


title: Molecular Structuring and Percolation Transition in Hydrated Sulfonated Poly(ether ether ketone) Membranes

abstract: The extent of phase separation and water percolation in sulfonated membranes are the key to their performance in fuel cells. Toward this, the effect of hydration on the morphology and transport characteristics of sulfonated poly(ether ether ketone), sPEEK, membrane is investigated using atomistic molecular dynamics (MD) simulation at various hydration levels ({\lambda}: number of water molecules per sulfonate group) between 4 and 15. At the molecular level, the evolution of local morphology is investigated in terms of structural pair correlations and minimum pair distances, and the transport properties are studied in terms of mean squared displacements (MSDs) and diffusion coefficients. The water-sulfur interaction in sPEEK is found to be stronger than that in Nafion, as observed in experiments. As opposed to Nafion, a weaker interaction of hydronium, with sulfonate, than water is observed. The behavior of water in sPEEK membrane is found to remain far from bulk as indicated by its diffusion coefficient. Analysis of simulation data indicate that at low {\lambda}, the largest water cluster forms a narrow connected path of water molecules and hydronium ions. With increasing {\lambda}, larger water domains appear, spanning more than half of the simulation box at {\lambda} = 15. Small isolated clusters are present at all hydration levels, demonstrating the extent of phase separation in sPEEK to be lesser than that in Nafion. Various analyses, both at molecular and collective level, suggest the occurrence of a percolation transition between {\lambda} = 8 and 10, which leads to a connected network of water channels in the membrane, thereby boosting the mobility of hydronium ions.


title: Embedded Spiral Patterns in the Cool Core of the Massive Cluster of Galaxies Abell 1835

abstract: We present the properties of intracluster medium (ICM) in the cool core of the massive cluster of galaxies Abell 1835 obtained with the data by $ Chandra$ $X$-$ray$ $Observatory$. We find distinctive spiral patterns with the radius of 70 kpc (or 18 arcsec) as a whole in the residual image of X-ray surface brightness after the 2-dimensional ellipse model of surface brightness is subtracted. The size is smaller by a factor of 2 -- 4 than that of other clusters known to have a similar pattern. The spiral patterns consist of two arms. One of them appears as positive, and the other does as negative excesses in the residual image. Their X-ray spectra show that the ICM temperatures in the positive- and negative-excess regions are $5.09^{+0.12}_{-0.13}$ keV and $6.52^{+0.18}_{-0.15}$ keV, respectively. In contrast, no significant difference is found in the abundance or pressure, the latter of which suggests that the ICM in the two regions of the spiral patterns is in pressure equilibrium or close. The spatially-resolved X-ray spectroscopy of the central region ($r<40$ arcsec) divided into 92 sub-regions indicates that Abell 1835 is a typical cool core cluster. We also find that the spiral patterns extend from the cool core out to the hotter surrounding ICM. The residual image reveals some lumpy sub-structure in the cool core. The line-of-sight component of the disturbance velocity responsible for the sub-structures is estimated to be lower than 600 km/s. Abell 1835 may be now experiencing an off-axis minor merger.


title: Heat kernel approach to relations between covariant and consistent currents in chiral gauge theories

abstract: We apply the heat kernel method to relations between covariant and consistent currents in anomalous chiral gauge theories. Banerjee et al. have shown that the relation between these currents is expressed by a "functional curl" of the covariant current. Using the heat kernel method, we evaluate the functional curl explicitly in arbitrary even dimensions. We also apply the heat kernel method to evaluate Osabe and Suzuki's results of the difference between covariant and consistent currents in two and four dimensions. Applying the arguments of Banerjee et al. to gravitational anomalies, we investigate the relationship between the covariant and consistent energy-momentum tensors. The relation is found to be expressed by a functional curl of the covariant energy-momentum tensor.


title: Robust mixture modelling using sub-Gaussian stable distribution

abstract: Heavy-tailed distributions are widely used in robust mixture modelling due to possessing thick tails. As a computationally tractable subclass of the stable distributions, sub-Gaussian $\alpha$-stable distribution received much interest in the literature. Here, we introduce a type of expectation maximization algorithm that estimates parameters of a mixture of sub-Gaussian stable distributions. A comparative study, in the presence of some well-known mixture models, is performed to show the robustness and performance of the mixture of sub-Gaussian $\alpha$-stable distributions for modelling, simulated, synthetic, and real data.


title: Two properties of maximal antichains in strict chain product posets

abstract: We present two results on maximal antichains in the strict chain product poset $[t_1+1]\times[t_2+1]\times\ldots\times[t_n+1]$. First, we prove that these maximal antichains are also maximum. Second, we prove that there is a bijection between maximal antichains in the strict chain product poset $[t_1+1]\times[t_2+1]\times\ldots\times[t_n+1]$ and antichains in the non-strict chain product poset $[t_1]\times[t_2]\times\ldots\times[t_n]$.


title: Collective Vertex Classification Using Recursive Neural Network

abstract: Collective classification of vertices is a task of assigning categories to each vertex in a graph based on both vertex attributes and link structure. Nevertheless, some existing approaches do not use the features of neighbouring vertices properly, due to the noise introduced by these features. In this paper, we propose a graph-based recursive neural network framework for collective vertex classification. In this framework, we generate hidden representations from both attributes of vertices and representations of neighbouring vertices via recursive neural networks. Under this framework, we explore two types of recursive neural units, naive recursive neural unit and long short-term memory unit. We have conducted experiments on four real-world network datasets. The experimental results show that our frame- work with long short-term memory model achieves better results and outperforms several competitive baseline methods.


title: Expected number of critical points of random holomorphic sections over complex projective space

abstract: We study the high dimensional asymptotics of the expected number of critical points of a given Morse index of Gaussian random holomorphic sections over complex projective space. We explicitly compute the exponential growth rate of the expected number of critical points of the largest index and of diverging indices at various rates as well as the exponential growth rate for the expected number of critical points (regardless of index). We also compute the distribution of the critical values for the expected number of critical points of smallest index.


title: Space of $C^2$-smooth geometrically continuous isogeometric functions on planar multi-patch geometries: Dimension and numerical experiments

abstract: We study the space of $C^{2}$-smooth isogeometric functions on bilinearly parameterized multi-patch domains $\Omega \subset \mathbb{R}^{2}$, where the graph of each isogeometric function is a multi-patch spline surface of bidegree $(d,d)$, $d \in \{5,6 \}$. The space is fully characterized by the equivalence of the $C^2$-smoothness of an isogeometric function and the $G^2$-smoothness of its graph surface, cf. (Groisser and Peters,2015; Kapl et al.,2015). This is the reason to call its functions $C^{2}$-smooth geometrically continuous isogeometric functions. In particular, the dimension of this $C^{2}$-smooth isogeometric space is investigated. The study is based on the decomposition of the space into three subspaces and is an extension of the work (Kapl and Vitrih, 2017) to the multi-patch case. In addition, we present an algorithm for the construction of a basis, and use the resulting globally $C^{2}$-smooth functions for numerical experiments, such as performing $L^{2}$ approximation and solving triharmonic equation, on bilinear multi-patch domains. The numerical results indicate optimal approximation order.


title: Estimating Time-Varying Effective Connectivity in High-Dimensional fMRI Data Using Regime-Switching Factor Models

abstract: Recent studies on analyzing dynamic brain connectivity rely on sliding-window analysis or time-varying coefficient models which are unable to capture both smooth and abrupt changes simultaneously. Emerging evidence suggests state-related changes in brain connectivity where dependence structure alternates between a finite number of latent states or regimes. Another challenge is inference of full-brain networks with large number of nodes. We employ a Markov-switching dynamic factor model in which the state-driven time-varying connectivity regimes of high-dimensional fMRI data are characterized by lower-dimensional common latent factors, following a regime-switching process. It enables a reliable, data-adaptive estimation of change-points of connectivity regimes and the massive dependencies associated with each regime. We consider the switching VAR to quantity the dynamic effective connectivity. We propose a three-step estimation procedure: (1) extracting the factors using principal component analysis (PCA) and (2) identifying dynamic connectivity states using the factor-based switching vector autoregressive (VAR) models in a state-space formulation using Kalman filter and expectation-maximization (EM) algorithm, and (3) constructing the high-dimensional connectivity metrics for each state based on subspace estimates. Simulation results show that our proposed estimator outperforms the K-means clustering of time-windowed coefficients, providing more accurate estimation of regime dynamics and connectivity metrics in high-dimensional settings. Applications to analyzing resting-state fMRI data identify dynamic changes in brain states during rest, and reveal distinct directed connectivity patterns and modular organization in resting-state networks across different states.


title: Low Frequency Observations of the Subpulse Drifter PSR J0034-0721 with the Murchison Widefield Array

abstract: The phenomenon of subpulse drifting may hold the key to understanding the pulsar emission mechanism. Here, we report on new observations of PSR J0034-0721 (B0031-07), carried out with the Murchison Widefield Array at 185 MHz. We observe three distinct drift modes whose "vertical" drift band separations ($P_3$) and relative abundances are consistent with previous studies at similar and higher frequencies. The driftbands, however, are observed to change their slopes over the course of individual drift modes, which can be interpreted as a continuously changing drift rate. The implied acceleration of the intrinsic carousel rotation cannot easily be explained by plasma models based on ExB drift. Furthermore, we find that methods of classifying the drift modes by means of $P_3$ measurements can sometimes produce erroneous identifications in the presence of a changing drift rate. The "horizontal" separation between driftbands ($P_2$) is found to be larger at later rotation phases within the pulse window, which is inconsistent with the established effects of retardation, aberration, and the motion of the visible point. Longer observations spanning at least ~10,000 pulses are required to determine how the carousel rotation parameters change from one drift sequence to the next.


title: Topological edge state engineering with high-frequency electromagnetic radiation

abstract: We outline here how strong light-matter interaction can be used to induce quantum phase transition between normal and topological phases in two-dimensional topological insulators. We consider the case of a HgTe quantum well, in which band inversion occurs above a critical value of the well thickness, and demonstrate that coupling between electron states and the $E$ field from an off-resonant linearly polarized laser provides a powerful tool to control topological transitions, even for a thickness of the quantum well that is below the critical value. We also show that topological phase properties of the edge states, including their group velocity, can be tuned in a controllable way by changing the intensity of the laser field. These findings open up the possibility for new experimental means with which to investigate topological insulators and shed new light on topological-insulator-based technologies that are under active discussion.


title: On hyperbolic cobweb manifolds

abstract: A compact hyperbolic "cobweb" manifold (hyperbolic space form) of symbol $Cw(6,6,6)$ will be constructed in Fig.1,4,5 as a representant of a presumably infinite series $Cw(2p,2p,2p)$ $(3 \le p \in \bN$ natural numbers). This is a by-product of our investigations \cite{MSz16}. In that work dense ball packings and coverings of hyperbolic space $\HYP$ have been constructed on the base of complete hyperbolic Coxeter orthoschemes $\mathcal{O}=W_{uvw}$ and its extended reflection groups $\bG$ (see diagram in Fig.~3. and picture of fundamental domain in Fig.~2). Now $u=v=w=6 (=2p)$. Thus the maximal ball contained in $Cw(6,6,6)$, moreover its minimal covering bal l (so diameter) can also be determined. The algorithmic procedure provides us with the proof of our statements.


title: Residual Gas & Dust Around Transition Objects and Weak T Tauri Stars

abstract: Residual gas in disks around young stars can spin down stars, circularize the orbits of terrestrial planets, and whisk away the dusty debris that is expected to serve as a signpost of terrestrial planet formation. We have carried out a sensitive search for residual gas and dust in the terrestrial planet region surrounding young stars ranging in age from a few Myr to ~10 Myr in age. Using high resolution 4.7 micron spectra of transition objects and weak T Tauri stars, we searched for weak continuum excesses and CO fundamental emission, after making a careful correction for the stellar contribution to the observed spectrum. We find that the CO emission from transition objects is weaker and located further from the star than CO emission from non-transition T Tauri stars with similar stellar accretion rates. The difference is possibly the result of chemical and/or dynamical effects (i.e., a low CO abundance or close-in low-mass planets). The weak T Tauri stars show no CO fundamental emission down to low flux levels (5 x 10^(-20) - 10^{-18} W/m^2). We illustrate how our results can be used to constrain the residual disk gas content in these systems and discuss their potential implications for star and planet formation.


title: The second knee observed in the local muon density spectra at various zenith angles

abstract: Local muon density spectra (LMDS) at various zenith angles have been reconstructed from the data of two detectors of the Experimental complex NEVOD. The inclined muon bundles at the ground level were detected with the coordinate detector DECOR, and for the near-vertical direction with the calibration telescope system (CTS) of the Cherenkov water detector. In comparison with the earlier DECOR results, the experimental statistics has been increased by 2-3 times for different ranges of zenith angle and muon bundle multiplicity and is now based on about 40,000 h of the setup operation. The live time of measurements with CTS is about 12,000 h. It is found with both setups that the slope of LMDS is increasing above the primary energy of about 1017 eV. The details of the experiment and data analysis are presented.


title: On the dense Preferential Attachment Graph models and their graphon induced counterpart

abstract: Letting $\mathcal{M}$ denote the space of finite measures on $\mathbb{N}$, and $\mu_\lambda\in\mathcal{M}$ denote the Poisson distribution with parameter $\lambda$, the function $W:[0,1]^2\to\mathcal{M}$ given by \[ W(x,y)=\mu_{c\log x\log y} \] is called the PAG graphon with density $c$. It is known that this is the limit, in the multigraph homomorphism sense, of the dense Preferential Attachment Graph (PAG) model with edge density $c$. This graphon can then in turn be used to generate the so-called W-random graphs in a natural way. The aim of this paper is to compare the dense PAG model with the W-random graph model obtained from the corresponding graphon. Motivated by the multigraph limit theory, we investigate the expected jumble norm distance of the two models in terms on the number of vertices $n$. We present a coupling for which the expectation can be bounded from above by $O(\log^2 n\cdot n^{-1/3})$, and provide a universal lower bound that is coupling independent, but with a worse exponent.


title: Octupolar Tensors for Liquid Crystals

abstract: A third-order three-dimensional symmetric traceless tensor, called the \emph{octupolar} tensor, has been introduced to study tetrahedratic nematic phases in liquid crystals. The octupolar \emph{potential}, a scalar-valued function generated on the unit sphere by that tensor, should ideally have four maxima capturing the most probable molecular orientations (on the vertices of a tetrahedron), but it was recently found to possess an equally generic variant with \emph{three} maxima instead of four. It was also shown that the irreducible admissible region for the octupolar tensor in a three-dimensional parameter space is bounded by a dome-shaped surface, beneath which is a \emph{separatrix} surface connecting the two generic octupolar states. The latter surface, which was obtained through numerical continuation, may be physically interpreted as marking a possible \emph{intra-octupolar} transition. In this paper, by using the resultant theory of algebraic geometry and the E-characteristic polynomial of spectral theory of tensors, we give a closed-form, algebraic expression for both the dome-shaped surface and the separatrix surface. This turns the envisaged intra-octupolar transition into a quantitative, possibly observable prediction. Some other properties of octupolar tensors are also studied.


title: Multiplicative partition functions for reverse plane partitions derived from an integrable dynamical system

abstract: A close connection of reverse plane partitions with an integrable dynamical system called the discrete two-dimensional (2D) Toda molecule is clarified. It is shown that a multiplicative partition function for reverse plane partition of arbitrary shape with bounded parts can be obtained from each non-vanishing solution to the discrete 2D Toda molecule. As an example a partition function which generalizes MacMahon's triple product formula as well as Gansner's multi-trace generating function is derived from a specific solution to the dynamical system.


title: Restrictions for the Causal Inferences in an Interferometric System

abstract: Causal discovery algorithms allow for the inference of causal structures from probabilistic relations of random variables. A natural field for the application of this tool is quantum mechanics, where a long-standing debate about the role of causality in the theory has flourished since its early days. In this paper, a causal discovery algorithm is applied in the search for causal models to describe a quantum version of Wheeler's delayed-choice experiment. The outputs explicitly show the restrictions for the introduction of classical concepts in this system. The exclusion of models with two hidden variables is one of them. A consequence of such a constraint is the impossibility to construct a causal model that avoids superluminal causation and assumes an objective view of the wave and particle properties simultaneously.


title: A new air-shower observable to constrain hadronic interaction models

abstract: The energy spectrum of muons at ground level in air showers are studied and a new observable is proposed to constrain hadronic interaction models used in air shower simulations. An asymmetric Gaussian function is proposed to describe the muon ground energy spectrum and its parameters are studied regarding primary particle, energy and hadronic interaction models. Based on two realistic measurements of the muon density at a given distance from the shower axis, a new observable (\rmu) is defined. Considering realistic values of detector resolutions and number of measured events, it is also shown \rmu can be successfully used to constrain low and high hadronic interaction models. The study is focused in the energy range between $10^{17.5}$ and $10^{18.0}$ eV because of the importance of this interval for particle physics and astrophysical models. The constraining power of the new observable is shown to be large within current experimental capabilities.


title: The core collapse supernova rate from 24 years of data of the Large Volume Detector

abstract: The Large Volume Detector (LVD) at INFN Laboratori Nazionali del Gran Sasso, Italy is a 1 kt liquid scintillator neutrino observatory mainly designed to study low energy neutrinos from Gravi- tational Stellar Collapses (GSC) with 100% efficiency over the entire Galaxy. Here we summarize the results of the search for supernova neutrino bursts over the full data set lasting from June 1992 to May 2016 for a total live time of 8211 days. In the lack of a positive observation, either in standalone mode or in coincidence with other experiments, we establish the upper limit to the rate of GSC event in the Milky Way: 0.1/year at 90% c.l.


title: Weak convergence on Wiener space: targeting the first two chaoses

abstract: We consider sequences of random variables living in a finite sum of Wiener chaoses. We find necessary and sufficient conditions for convergence in law to a target variable living in the sum of the first two Wiener chaoses. Our conditions hold notably for sequences of multiple Wiener integrals. Malliavin calculus and in particular the Gamma-operators are used. Our results extend previous findings by Azmoodeh, Peccati and Poly (2014) and are applied to central and non-central convergence situations. Our methods are applied as well to investigate stable convergence. We finally exclude certain classes of random variables as target variables for sequences living in a fixed Wiener chaos.


title: Towards the quality improvement of cross-platform mobile applications

abstract: During last ten years, the number of smartphones and mobile applications has been constantly growing. Android, iOS and Windows Mobile are three mobile platforms that cover almost all smartphones in the world in 2017. Developing a mobile app involves first to choose the platforms the app will run, and then to develop specific solutions (i.e., native apps) for each chosen platform using platform-related toolkits such as AndroidSDK. Across-platform mobile application is an app that runs on two or more mobile platforms. Several frameworks have been proposed to simplify the development of cross-platform mobile applications and to reduce development and maintenance costs.They are called cross-platform mobile app development frameworks.However, to our knowledge, the life-cycle and the quality of cross-platforms mobile applications built using those frameworks have not been studied in depth. Our main goal is to first study the processes of development and maintenance of mobile applications built using cross-platform mobile app development frameworks, focusing particularly on the bug-fixing activity. Then, we aim at defining tools for automated repairing bugs from cross-platform mobile applications.


title: Impact of leptoquarks in semileptonic B decays

abstract: We study the impact of vector leptoquarks on the recent anomalies in semileptonic $B$ meson decays such as $R_K$ and $R_{D^{(*)}}$ etc. We constrain the leptoquark couplings by using the existing data on the branching ratios of $B_s \to l^+ l^-$ and $\bar{B} \to X_s l^+ l^-$ processes, where $l=e,\mu,\tau$. We estimate the branching ratios of $\bar B \to D^{(*)} l \bar \nu_l$ processes using the constrained leptoquark parameter space. We also investigate the possibility of simultaneous explanation of $R_{D^{(*)}}$, $R_K$ anomalies in this model.


title: Density of states of Dirac-Landau levels in a gapped graphene monolayer under strain gradient

abstract: We study a gapped graphene monolayer in a combination of uniform magnetic field and strain-induced uniform pseudomagnetic field. The presence of two fields completely removes the valley degeneracy. The resulting density of states shows a complicated behaviour that can be tuned by adjusting the strength of the fields. We analyze how these features can be observed in the sublattice, valley and full density of states. The analytical expression for the valley DOS is derived.


title: Analysis of Breakdown Probability of Wireless Sensor Networks with Unreliable Relay Nodes

abstract: In the present paper, we derive an upper bound of the average network breakdown probability of packet networks with unreliable relay nodes. We here assume that relay nodes get independently broken with a given node breakdown probability. A survivor graph is the induced subgraph obtained by removing the broken relay nodes and their connecting edges from the original graph. If the survivor network is disconnected, we consider a network breakdown happens. The primal contribution of the paper is to derive an upper bound of the average network breakdown probability, where the expectation is taken over a regular graph ensemble. The proof of the bound is based on a natural one-to-one correspondence between a regular graph and a regular bipartite graph, and also on enumeration of bipartite graphs satisfying certain conditions. This proof argument is inspired by the analysis of weight distribution for low-density parity-check codes. Compared with estimates of the average network breakdown probability obtained by computer experiments, it is observed that the upper bound provides the values which are not only upper bounds but also precise estimates of the network breakdown probability when the node breakdown probability is small.


title: Community detection based on significance optimization in complex networks

abstract: Community structure is an important structural property that extensively exists in various complex networks. In the past decade, much attention has been paid to the design of community-detection methods, but analyzing the behaviors of the methods is also of interest in the theoretical research and real applications. Here, we focus on an important measure for community structure, significance [Sci. Rep. 3 (2013) 2930]. Specifically, we study the effect of various network parameters on this measure in detail, analyze the critical behaviors of it in partition transition, and analytically give the formula of the critical points and the phase diagrams. The results shows that the critical number of communities in partition transition increases dramatically with the difference between inter- and intra-community link densities, and thus significance optimization displays higher resolution in community detection than many other methods, but it is also easily to lead to the excessive splitting of communities. By Louvain algorithm for significance optimization, we confirmed the theoretical results on artificial and real-world networks, and give a series of comparisons with some classical methods.


title: Training Group Orthogonal Neural Networks with Privileged Information

abstract: Learning rich and diverse representations is critical for the performance of deep convolutional neural networks (CNNs). In this paper, we consider how to use privileged information to promote inherent diversity of a single CNN model such that the model can learn better representations and offer stronger generalization ability. To this end, we propose a novel group orthogonal convolutional neural network (GoCNN) that learns untangled representations within each layer by exploiting provided privileged information and enhances representation diversity effectively. We take image classification as an example where image segmentation annotations are used as privileged information during the training process. Experiments on two benchmark datasets -- ImageNet and PASCAL VOC -- clearly demonstrate the strong generalization ability of our proposed GoCNN model. On the ImageNet dataset, GoCNN improves the performance of state-of-the-art ResNet-152 model by absolute value of 1.2% while only uses privileged information of 10% of the training images, confirming effectiveness of GoCNN on utilizing available privileged knowledge to train better CNNs.


title: Digit frequencies and self-affine sets with non-empty interior

abstract: In this paper we study digit frequencies in the setting of expansions in non-integer bases, and self-affine sets with non-empty interior. Within expansions in non-integer bases we show that if $\beta\in(1,1.787\ldots)$ then every $x\in(0,\frac{1}{\beta-1})$ has a simply normal $\beta$-expansion. We also prove that if $\beta\in(1,\frac{1+\sqrt{5}}{2})$ then every $x\in(0,\frac{1}{\beta-1})$ has a $\beta$-expansion for which the digit frequency does not exist, and a $\beta$-expansion with limiting frequency of zeros $p$, where $p$ is any real number sufficiently close to $1/2$. For a class of planar self-affine sets we show that if the horizontal contraction lies in a certain parameter space and the vertical contractions are sufficiently close to $1,$ then every nontrivial vertical fibre contains an interval. Our approach lends itself to explicit calculation and give rise to new examples of self-affine sets with non-empty interior. One particular strength of our approach is that it allows for different rates of contraction in the vertical direction.


title: Properties of AGN coronae in the NuSTAR era II: hybrid plasma

abstract: The corona, a hot cloud of electrons close to the centre of the accretion disc, produces the hard X-ray power-law continuum commonly seen in luminous Active Galactic Nuclei (AGN). The continuum has a high-energy turnover, typically in the range of one to several 100 keV and is suggestive of Comptonization by thermal electrons. We are studying hard X-ray spectra of AGN obtained with NuSTAR after correction for X-ray reflection and under the assumption that coronae are compact, being only a few gravitational radii in size as indicated by reflection and reverberation modelling. Compact coronae raise the possibility that the temperature is limited and indeed controlled by electron-positron pair production, as explored earlier (Paper I). Here we examine hybrid plasmas in which a mixture of thermal and nonthermal particles is present. Pair production from the nonthermal component reduces the temperature leading to a wider temperature range more consistent with observations.


title: Dynamics of a Camphoric Acid boat at the air-water interface

abstract: We report experiments on an agarose gel tablet loaded with camphoric acid (c-boat) set into self-motion by interfacial tension gradients at the air-water interface. We observe three distinct modes of c-boat motion: harmonic mode where the c-boat speed oscillates sinusoidally in time, a steady mode where the c-boat maintains constant speed, and a relaxation oscillation mode where the c-boat maintains near-zero speed between sudden jumps in speed and position at regular time intervals. Whereas all three modes have been separately reported before in different systems, we show they belong to a common description. Through control of the air-water surface tension with Sodium Dodecyl Sulfate (SDS), we experimentally deduce the three self-propulsive modes result from surface tension difference between Camphoric Acid (CA) and the ambient surroundings.


title: A generalized linear model for decomposing cis-regulatory, parent-of-origin, and maternal effects on allele-specific gene expression

abstract: Joint quantification of genetic and epigenetic effects on gene expression is important for understanding the establishment of complex gene regulation systems in living organisms. In particular, genomic imprinting and maternal effects play important roles in the developmental process of mammals and flowering plants. However, the influence of these effects on gene expression are difficult to quantify because they act simultaneously with cis-regulatory mutations. Here we propose a simple method to decompose cis-regulatory (i.e., allelic genotype, AG), genomic imprinting (i.e., parent-of-origin, PO), and maternal (i.e., maternal genotype, MG) effects on allele-specific gene expression using RNA-seq data obtained from reciprocal crosses. We evaluated the efficiency of method using a simulated dataset and applied the method to whole-body Drosophila and mouse trophoblast stem cell (TSC) and liver RNA-seq data. Consistent with previous studies, we found little evidence of PO and MG effects in adult Drosophila samples. In contrast, we identified dozens and hundreds of mouse genes with significant PO and MG effects, respectively. Interestingly, a similar number of genes with significant PO effect were detect in mouse TSCs and livers, whereas more genes with significant MG effect were observed in livers. Further application of this method will clarify how these three effects influence gene expression levels in different tissues and developmental stages, and provide novel insight into the evolution of gene expression regulation.


title: The N/O ratio in early B-type main sequence stars as an indicator of their evolution

abstract: It is shown that, in the case of early B-type main-sequence stars, of the three ratios N/C, C/O, and N/O which are regarded as indicators of stellar evolution, the ratio N/O is more reliable since it seems to be insensitive to overionization of the NII and OII ions. On the other hand, the N/C and C/O ratios, which include carbon, may contain systematic errors for stars with Teff > 18500 K because of neglected overionization of CII ions. The ratio N/O is studied in the atmospheres of 46 early main-sequence B stars. These values of N/O are examined as functions of the effective temperature, age, rotation speed, and mass of the stars. Most early B-stars in the main sequence are found to have [N / O] = 0 , which indicates that N/O varies little during the main sequence stage, and this result is independent of the basic parameters listed above. There are two explanations for the large number of stars with [N / O] = 0 : it is predicted theoretically that for an initial rotation velocity V0 < 100 km/s, N/O varies little toward the end of the main sequence stage ([N/O] < 0.2) and observations show that most early main-sequence B-stars do actually have low initial rotation velocities V0. The few early main-sequence B-stars with higher [N/O] = 0.4-0.8 correspond to models with rotation velocities V0 = 200-300 km/s. This conclusion is consistent with earlier data for stars with the same masses in a later stage of evolution: the AFG-supergiant and bright giant stage.


title: Truncation Dimension for Linear Problems on Multivariate Function Spaces

abstract: The paper considers linear problems on weighted spaces of multivariate functions of many variables. The main questions addressed are: When is it possible to approximate the solution for the original function of very many variables by the solution for the same function; however with all but the first $k$ variables set to zero, so that the corresponding error is small? What is the truncation dimension, i.e., the smallest number $k=k(\varepsilon)$ such that the corresponding error is bounded by a given error demand $\varepsilon$? Surprisingly, $k(\varepsilon)$ could be very small even for weights with a modest speed of convergence to zero.


title: Monotone Martingale Transport Plans and Skorohod Embedding

abstract: We show that the left-monotone martingale coupling is optimal for any given performance function satisfying the martingale version of the Spence-Mirrlees condition, without assuming additional structural conditions on the marginals. We also give a new interpretation of the left monotone coupling in terms of Skorokhod embedding which allows us to give a short proof of uniqueness.


title: Sparse reconstruction of electric fields from radial magnetic data

abstract: Accurate estimates of the horizontal electric field on the Sun's visible surface are important not only for estimating the Poynting flux of magnetic energy into the corona but also for driving time-dependent magnetohydrodynamic models of the corona. In this paper, a method is developed for estimating the horizontal electric field from a sequence of radial-component magnetic field maps. This problem of inverting Faraday's law has no unique solution. Unfortunately, the simplest solution (a divergence-free electric field) is not realistically localized in regions of non-zero magnetic field, as would be expected from Ohm's law. Our new method generates instead a localized solution, using a basis pursuit algorithm to find a sparse solution for the electric field. The method is shown to perform well on test cases where the input magnetic maps are flux balanced, in both Cartesian and spherical geometries. However, we show that if the input maps have a significant imbalance of flux - usually arising from data assimilation - then it is not possible to find a localized, realistic, electric field solution. This is the main obstacle to driving coronal models from time sequences of solar surface magnetic maps.


title: Microenvironmental cooperation promotes early spread and bistability of a Warburg-like phenotype

abstract: We introduce an in silico model for the initial spread of an aberrant phenotype with Warburg-like overflow metabolism within a healthy homeostatic tissue in contact with a nutrient reservoir (the blood), aimed at characterizing the role of the microenvironment for aberrant growth. Accounting for cellular metabolic activity, competition for nutrients, spatial diffusion and their feedbacks on aberrant replication and death rates, we obtain a phase portrait where distinct asymptotic whole-tissue states are found upon varying the tissue-blood turnover rate and the level of blood-borne primary nutrient. Over a broad range of parameters, the spreading dynamics is bistable as random fluctuations can impact the final state of the tissue. Such a behaviour turns out to be linked to the re-cycling of overflow products by non-aberrant cells. Quantitative insight on the overall emerging picture is provided by a spatially homogeneous version of the model.

categories: q-bio.MN cond-mat.dis-nn physics.bio-ph q-bio.TO


title: Scaling Fixed-Field Alternating-Gradient accelerators with reverse bend and spiral edge angle

abstract: A novel scaling type of Fixed-Field Alternating-Gradient (FFAG) accelerator is proposed that solves the major problems of conventional scaling FFAGs. This scaling FFAG accelerator combines reverse bending magnets of the radial sector type and a spiral edge angle of the spiral sector type to ensure sufficient vertical focusing without relying on extreme values of either parameter. This new concept makes it possible to design a scaling FFAG for high energy (above GeV range) applications such as a proton driver for a spallation neutron source and an accelerator driven subcritical reactor.


title: The Autonomic Architecture of the Licas System

abstract: Licas (lightweight internet-based communication for autonomic services) is a distributed framework for building service-based systems. The framework provides a p2p server and more intelligent processing of information through its AI algorithms. Distributed communication includes XML-RPC, REST, HTTP and Web Services. It can now provide a robust platform for building different types of system, where Microservices or SOA would be possible. However, the system may be equally suited for the IoT, as it provides classes to connect with external sources and has an optional Autonomic Manager with a MAPE control loop integrated into the communication process. The system is also mobile-compatible with Android. This paper focuses in particular on the autonomic setup and how that might be used. A novel linking mechanism has been described previously and is considered again, as part of the autonomous framework.


title: Proof of Koml\'os's conjecture on Hamiltonian subsets

abstract: Koml\'os conjectured in 1981 that among all graphs with minimum degree at least $d$, the complete graph $K_{d+1}$ minimises the number of Hamiltonian subsets, where a subset of vertices is Hamiltonian if it contains a spanning cycle. We prove this conjecture when $d$ is sufficiently large. In fact we prove a stronger result: for large $d$, any graph $G$ with average degree at least $d$ contains almost twice as many Hamiltonian subsets as $K_{d+1}$, unless $G$ is isomorphic to $K_{d+1}$ or a certain other graph which we specify.


title: Diffeological Dirac operators and diffeological gluing

abstract: This manuscript attempts to present a way in which the classical construction of the Dirac operator can be carried over to the setting of diffeology. A more specific aim is to describe a procedure for gluing together two usual Dirac operators and to explain in what sense the result is again a Dirac operator. Since versions of cut-and-paste (surgery) operations have already appeared in the context of Atiyah-Singer theory, we specify that our gluing procedure is designed to lead to spaces that are not smooth manifolds in any ordinary sense, and since much attention has been paid in recent years to Dirac operators on spaces with singularities, we also specify that our approach is more of a piecewise-linear nature (although, hopefully, singular spaces in a more analytic sense will enter the picture sooner or later; but this work is not yet about them). Most of it is devoted to the diffeological versions of the components that go into the standard definition of a Dirac operator as the composition of a Clifford connection with Clifford action by sections of the cotangent bundle; a diffeological Dirac operator is then standardly defined.


title: Nucleosynthesis in the Innermost Ejecta of Neutrino-Drive Supernova Explosions in Two Dimensions

abstract: We examine the nucleosynthesis in the innermost, neutrino-processed ejecta (a few $10^{-3}\,M_\odot$) of self-consistent, two-dimensional explosion models of core-collapse supernovae for six progenitor stars with different initial masses. Three models have initial masses near the low-mass end of the supernova range, $8.8\,M_\odot$ (e8.8; electron-capture supernova), $9.6\,M_\odot$ (z9.6), and $8.1\,M_\odot$ (u8.1), with initial metallicities of 1, 0, and $10^{-4}$ times the solar metallicity, respectively. The other three are solar-metallicity models with initial masses of $11.2\,M_\odot$ (s11), $15\,M_\odot$ (s15), and $27\,M_\odot$ (s27). The low-mass models e8.8, z9.6, and u8.1 exhibit high production factors (nucleosynthetic abundances relative to the solar ones) of 100--200 for light trans-iron elements from Zn to Zr. This is associated with appreciable ejection of neutron-rich matter in these models. Remarkably, the nucleosynthetic outcomes for progenitors e8.8 and z9.6 are almost identical, including interesting productions of $^{48}$Ca and $^{60}$Fe, irrespective of their quite different (O-Ne-Mg and Fe) cores prior to collapse. In the more massive models s11, s15, and s27, several proton-rich isotopes of light trans-iron elements, including the $p$-isotope $^{92}$Mo (for s27) are made, up to production factors of $\sim$30. Both electron-capture and core-collapse supernovae near the low-mass end can therefore be dominant contributors to the Galactic inventory of light trans-iron elements from Zn to Zr and probably $^{48}$Ca and live $^{60}$Fe. The innermost ejecta of more massive supernovae may have only sub-dominant contributions to the chemical enrichment of the Galaxy except for $^{92}$Mo.