title: A Highly Collimated Jet from the Red Square Nebula, MWC 922

abstract: Deep, narrow-band \Ha\ and 6584 \AA\ [\Nii ] CCD images of the peculiar, infrared excess B[e] star MWC~922 reveal a collimated, bipolar jet orthogonal to the previously detected extended nebula. The jet consists of a pair of $\sim$0.15 pc segments on either side of MWC~922 separated by gaps. The most distant jet segments disappear $\sim$0.6 pc from the star. The northwest beam points to a faint emission-line feature 1.65 pc from MWC~922 that may be a terminal bow shock where the jet rams the ambient medium. The narrow opening-angle of the jet combined with an estimated internal sound speed of $\sim$10 \kms\ implies a jet speed $\sim$500 \kms . The previously detected nebula extends up to 0.6~pc to the southwest of MWC~922 at right angles to the jet and appears to be an extension of the compact, edge-on disk surrounding the star. It points toward the HII region Messier 16 located $\sim$1\degr\ ($\sim$30 pc in projection) to the southwest. This nebula and jet appear to be externally ionized by the ambient Lyman continuum radiation field and have electron densities of n$_e \sim$ 50 to 100 cm$^{-3}$. The southwest nebula and jet have similar surface brightness in \Ha\ and [\Nii ]. Faint 70 $\mu$m emission traces the southwest ejecta that likely originates from $\sim$50 K dust embedded in the photo-ionized plasma which may shadow the dimmer ejecta northeast of MWC~922. MWC~922 may be a massive member of the Serpens OB1 or OB2 associations surrounding Messier 16 and Sh2-54.


title: Readership Data and Research Impact

abstract: Reading academic publications is a key scholarly activity. Scholars accessing and recording academic publications online are producing new types of readership data. These include publisher, repository, and academic social network download statistics as well as online reference manager records. This chapter discusses the use of download and reference manager data for research evaluation and library collection development. The focus is on the validity and application of readership data as an impact indicator for academic publications across different disciplines. Mendeley is particularly promising in this regard, although all data sources are not subjected to rigorous quality control and can be manipulated.


title: Equal heartbeat intervals and their effects on the nonlinearity of permutation-based time irreversibility in heart rate

abstract: The association of equal heartbeat intervals with cardiac conditions and the effect of the equality on permutation-based time irreversibility are investigated in this paper. We measure the distributions of equal heartbeat intervals under three conditions, namely congestive heart failure (CHF), healthy young and elderly, whose time irreversibility of heartbeats is detected by measuring the probabilistic difference between permutations instead of raw vectors. We demonstrate that heartbeats contain high rates of equal states, particularly the CHF with around 20% equalities, and the distributions of equal values discriminate the heartbeats at very short data length. The CHF have more equal values than the healthy young (p<1.47*10^{-15}) and elderly (p<2.48*10^{-11}), and the healthy young have less equalities than the elderly (p<3.16*10^{-4}). Time irreversibility considering equal values is promising to extract nonlinear behaviors of heartbeats, confirming the decreased nonlinear complexity of the diseased and aging heart rates, while that involving no equality leads to erroneous nonlinearity detection. In our contribution, we highlight the pathological or physiological information contained by the distribution of equal heartbeat intervals that might contribute to develop relevant biomarkers in the area of heart analysis, and demonstrate the effectiveness of equality-based time irreversibility in the nonlinearity detection of heartbeats.


title: Weak Lensing Effect on CMB in the Presence of a Dipole Anisotropy

abstract: We investigate weak lensing effect on cosmic microwave background (CMB) in the presence of dipole anisotropy. The approach of flat-sky approximation is considered. We determine the functions $\sigma_0^2$ and $\sigma_2^2$ that appear in expressions of the lensed CMB power spectrum in the presence of a dipole anisotropy. We determine the correction to B-mode power spectrum which is found to be appreciable at low multipoles ($l$). However, the temperature and E-mode power spectrum are not altered significantly.


title: Observable signatures of dark photons from supernovae

abstract: A dark photon is a well-motivated new particle which, as a component of an associated dark sector, could explain dark matter. One strong limit on dark photons arises from excessive cooling of supernovae. We point out that even at couplings where too few dark photons are produced in supernovae to violate the cooling bound, they can be observed directly through their decays. Supernovae produce dark photons which decay to positrons, giving a signal in the 511 keV annihilation line observed by SPI/INTEGRAL. Further, prompt gamma-ray emission by these decaying dark photons gives a signal for gamma-ray telescopes. Existing GRS observations of SN1987a already constrain this, and a future nearby SN could provide a detection. Finally, dark photon decays from extragalactic SN would produce a diffuse flux of gamma rays observable by detectors such as SMM and HEAO-1. Together these observations can probe dark photon couplings several orders of magnitude beyond current constraints for masses of roughly 1 - 100 MeV.


title: Realizing Hopf Insulators in Dipolar Spin Systems

abstract: The Hopf insulator is a weak topological insulator characterized by an insulating bulk with conducting edge states protected by an integer-valued linking number invariant. The state exists in three-dimensional two-band models. We demonstrate that the Hopf insulator can be naturally realized in lattices of dipolar-interacting spins, where spin exchange plays the role of particle hopping. The long-ranged, anisotropic nature of the dipole-dipole interactions allows for the precise detail required in the momentum-space structure, while different spin orientations ensure the necessary structure of the complex phases of the hoppings. Our model features robust gapless edge states at both smooth edges, as well as sharp edges obeying a certain crystalline symmetry, despite the breakdown of the two-band picture at the latter. In a companion manuscript [2105.10504], we provide a specific experimental blueprint for implementing our proposal using ultracold polar molecules of $^{40}$K$^{87}$Rb.


title: Integrability and Holographic Aspects of Six-Dimensional ${\cal N}=(1,0)$ Superconformal Field Theories

abstract: In the framework of six-dimensional conformal field theories with ${\cal N}=(1,0)$ supersymmetry we develop the map between the holographic description, the field theoretical description and the associated Hanany-Witten set-ups. General expressions that calculate various observables are presented. The study of string solitons singles out a special background of Massive IIA on which we show (by explicitly finding a Lax pair) that the Neveu-Schwarz part of the string sigma model is classically integrable. We study the particular dual conformal field theory and compute some of its observables.


title: Characterisation of multi-level quantum coherence without ideal measurements

abstract: Coherent superpositions are one of the hallmarks of quantum mechanics and are vital for any quantum mechanical device to outperform the classically achievable. Generically, superpositions are verified in interference experiments, but despite their longstanding central role we know very little about how to extract the number of coherently superposed amplitudes from a general interference pattern. A fundamental issue is that performing a phase-sensitive measurement is as challenging as creating a coherent superposition, so that assuming a perfectly implemented measurement for verification of quantum coherence is hard to justify. In order to overcome this issue, we construct a coherence certifier derived from simple statistical properties of an interference pattern, such that any imperfection in the measurement can never over-estimate the number of coherently superposed amplitudes. We numerically test how robust this measure is to under-estimating the coherence in the case of imperfect state preparation or measurement, and find it to be very resilient in both cases.


title: Small-scale chemical abundance analysis in a blue compact dwarf galaxy SBS 1415+437

abstract: We use integral field spectroscopic (IFS) observations from Gemini Multi-Object Spectrograph-North (GMOS-N) to analyse the ionised gas in the principal star-forming region in the blue compact dwarf galaxy SBS 1415+437. The IFS data enable us to map the weak auroral line [O III] $\lambda$4363 at a spatial scale of $\sim$6.5 pc across a region of $\sim$143 $\times$ 143 pc$^2$. This in turn allows us to use the robust direct T$_e$-method to map the ionic and elemental abundances of nitrogen (N) along with the alpha-elements, oxygen (O), neon (Ne), sulphur (S) and argon (Ar). We utilise these abundances to map the relative abundances of N, Ne, S and Ar with respect to O. We segment this predominantly photoionised region of study into elliptical annuli on the basis of the H$\alpha$ flux distribution to study the variation of chemical abundances and their ratios, and find no significant chemical variation. We also perform chemical abundance analysis on the integrated spectra of the region under study and elliptical annuli within it. We find that the inferred abundances are in agreement with the median of the abundances obtained from the chemical abundance maps of the principal star-forming region and the mapped values within annuli.The finding has important implications for direct comparison with high-redshift observations, where spatial resolution is not available, and for a consistent approach to track chemical evolution across cosmic time.


title: DASH: a library of dynamical subhalo evolution

abstract: The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological $N$-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that accurately describe subhalo evolution exist either. We therefore present a library of idealized, high resolution $N$-body simulations of the tidal evolution of individual subhaloes that can be used to calibrate semi-analytical models and to complement cosmological simulations. The simulations focus on minor mergers, i.e., the mass of the subhalo is much smaller than that of the host halo, such that the impact of dynamical friction is negligible. This setup allows the adoption of a fixed analytical potential for modelling the host halo. The dynamical evolution of subhaloes is followed with $N$-body computations. In the library, four parameters, two of which characterize the subhalo orbit with respect to the host halo, and the two concentrations of the host- and subhalo, are varied over the ranges encountered in cosmological simulations. We show several representative examples from the library that illustrate the evolution of the subhalo mass and velocity dispersion profiles. Additionally, we make publicly available a pre-trained non-parametric model of the subhalo mass evolution based on random forest regression. This model is able to interpolate the simulation data at the 0.1\,dex level and provides efficient access to the data for further use in modelling.


title: $n$-Hourglass Weyl fermions in nonsymmorphic materials

abstract: Hourglass-like band structures protected by nonsymmorphic space group symmetries can appear along high-symmetry lines or in high-symmetry surfaces in the Brillouin zone. In this work, from symmetry analysis, we demonstrate that $n$-hourglass-like band structures, a generalization of hourglass-like band structures, which host a number of Weyl points, are enforced along screw-invariant lines in non-magnetic materials with a single $N$-fold screw axis when spin-orbit coupling is finite, where $n$, a non-unity factor of $N$, denotes the degree of the screw-invariant line. The "standard" $n$-hourglass has $n-1$ crossings, which is minimal, and its variants can have more crossings. Purely by symmetry considerations, we find there are minimally two particle and two hole Fermi pockets enclosing Weyl points with opposite monopole charges at proper fillings, which can result in distinct physical effects including the possible formation of topological density waves and the quantum nonlinear Hall effect. We construct a minimal model which respects all the symmetries, and from which we see how the $n$-hourglasses appear when spin-orbit coupling is turned on. The same results are derived from compatibility relations. As exemplary, BiPd, ZnTe under high pressure, and the high-temperature phase of Tl$_3$PbBr$_5$, are shown from first-principles calculations to exhibit $n$-hourglass-like band structures, with $N=2,3,4$, respectively, which confirms our symmetry analysis and minimal model.


title: Hyperinflation generalised: from its attractor mechanism to its tension with the `swampland conjectures'

abstract: In negatively curved field spaces, inflation can be realised even in steep potentials. Hyperinflation invokes the `centrifugal force' of a field orbiting the hyperbolic plane to sustain inflation. We generalise hyperinflation by showing that it can be realised in models with any number of fields ($N_f\geq2$), and in broad classes of potentials that, in particular, don't need to be rotationally symmetric. For example, hyperinflation can follow a period of radial slow-roll inflation that undergoes geometric destabilisation, yet this inflationary phase is not identical to the recently proposed scenario of `side-tracked inflation'. We furthermore provide a detailed proof of the attractor mechanism of (the original and generalised) hyperinflation, and provide a novel set of characteristic, explicit models. We close by discussing the compatibility of hyperinflation with observations and the recently much discussed `swampland conjectures'. Observationally viable models can be realised that satisfy either the `de Sitter conjecture' ($V'/V\gtrsim 1$) or the `distance conjecture' ($\Delta \phi \lesssim 1$), but satisfying both simultaneously brings hyperinflation in some tension with successful reheating after inflation. However, hyperinflation can get much closer to satisfying all of these criteria than standard slow-roll inflation. Furthermore, while the original model is in stark tension with the weak gravity conjecture, generalisations can circumvent this issue.


title: Implications on Spatial Models of Interstellar Gamma-Ray Inverse-Compton Emission from Synchrotron Emission Studies in Radio and Microwaves

abstract: Cosmic rays interacting with gas and photon fields in the Galaxy produce interstellar gamma-ray emission (IGE), which accounts for almost 50% of the photons detected at gamma-ray energies. Models of this IGE have to be very accurate for interpreting the high-quality observations by present gamma-ray telescopes, such as Fermi Large Area Telescope (LAT). Standard models of IGE, used as reference models for analyses of the Fermi LAT data, show spatial discrepancies with respect to the data, underlining the necessity of more realistic models. The same CR electrons that produce the inverse-Compton component of the IGE produce also interstellar synchrotron emission observed in radio and microwave. However, present standard models do not take advantage of results coming from studies of this interstellar synchrotron emission. Accounting for such results, in this work we show how they affect the calculated spatial maps of the large-scale inverse-Compton component of the IGE, which are usually used in studies of Fermi LAT data. It is found that these results significantly affect these spatial model maps even at a 60% level. In particular, propagation models based on synchrotron studies produce a more peaked inverse-Compton emission in the inner Galaxy region with respect to the standard models used to analyze Fermi LAT data. The conclusion is that radio and microwave observations can be included in a multifrequency self-consistent approach for a more accurate modeling of the IGE finalized to a physical comprehensive interpretation of gamma-ray data and its present unexplained features. Model parameters are provided, which supply a more realistic basis for high-energy gamma-ray studies.


title: Investing for Discovery and Sustainability in Astronomy in the 2020s

abstract: As the next decade approaches, it is once again time for the US astronomical community to assess its investment priorities for the coming decade on the ground and in space. This report, created to aid NOAO in its planning for the 2020 Decadal Survey on Astronomy and Astrophysics, reviews the outcome of the previous Decadal Survey (Astro2010); describes the themes that emerged from the 2018 NOAO community planning workshop "NOAO Community Needs for Science in the 2020s"; and based on the above, offers thoughts for the coming review. We find that a balanced set of investments in small- to large-scale initiatives is essential to a sustainable future, based on the experience of previous decades. While large facilities are the "value" investments that are guaranteed to produce compelling science and discoveries, smaller facilities are the "growth stocks" that are likely to deliver the biggest science bang per buck, sometimes with outsize returns. Investments in data-intensive missions also have benefits to society beyond the science they deliver. By training scientists who are well equipped to use their data science skills to solve problems in the public or private sector, astronomy can provide a valuable service to society by contributing to a data-capable workforce.


title: Fast Markov Chain Monte Carlo Algorithms via Lie Groups

abstract: From basic considerations of the Lie group that preserves a target probability measure, we derive the Barker, Metropolis, and ensemble Markov chain Monte Carlo (MCMC) algorithms, as well as variants of waste-recycling Metropolis-Hastings and an altogether new MCMC algorithm. We illustrate these constructions with explicit numerical computations, and we empirically demonstrate on a spin glass that the new algorithm converges more quickly than its siblings.


title: Fate of dynamical phases of a BCS superconductor beyond the dissipationless regimen

abstract: The BCS model of an isolated superconductor initially prepared in a nonequilibrium state, predicts the existence of interesting dynamical phenomena in the time-dependent order parameter as decaying oscillations, persistent oscillations and overdamped dynamics. To make contact with real systems remains an open challenge as one needs to introduce dissipation due to the environment in a self-consistent computation. Here, we reach this goal with the use of the Keldysh formalism to treat the effect of a thermal bath. We show that, contrary to the dissipationless case, all dynamical phases reach the equilibrium order parameter in a characteristic time that depends on the coupling with the bath. Remarkably, as time evolves, the overdamped phase shows a fast crossover where the superconducting order parameter recovers to reach a state with a well-developed long range order that tends towards equilibrium with the damped Higgs mode oscillations. Our results provide a benchmark for the description of the dynamics of real out-of-equilibrium superconductors relevant for quantum technological applications.


title: Multi-stream Network With Temporal Attention For Environmental Sound Classification

abstract: Environmental sound classification systems often do not perform robustly across different sound classification tasks and audio signals of varying temporal structures. We introduce a multi-stream convolutional neural network with temporal attention that addresses these problems. The network relies on three input streams consisting of raw audio and spectral features and utilizes a temporal attention function computed from energy changes over time. Training and classification utilizes decision fusion and data augmentation techniques that incorporate uncertainty. We evaluate this network on three commonly used data sets for environmental sound and audio scene classification and achieve new state-of-the-art performance without any changes in network architecture or front-end preprocessing, thus demonstrating better generalizability.


title: General field theory and weak Euler-Lagrange equation for classical particle-field systems in plasma physics

abstract: A general field theory for classical particle-field systems is developed. Compared with the standard classical field theory, the distinguish feature of a classical particle-field system is that the particles and fields reside on different manifolds. The fields are defined on the 4D space-time, whereas each particle's trajectory is defined on the 1D time-axis. As a consequence, the standard Noether's procedure for deriving local conservation laws in space-time from symmetries is not applicable without modification. To overcome this difficulty, a weak Euler-Lagrange equation for particles is developed on the 4D space-time, which plays a pivotal role in establishing the connections between symmetries and local conservation laws in space-time. Especially, the non-vanishing Euler derivative in the weak Euler-Lagrangian equation generates a new current in the conservation laws. Several examples from plasma physics are studied as special cases of the general field theory. In particular, the relations between the rotational symmetry and angular momentum conservation for the Klimontovich-Poisson system and the Klimontovich-Darwin system are established.


title: Triple-lens Gravitational Microlensing: Critical Curves for Arbitrary Spatial Configuration

abstract: Since the first observation of triple-lens gravitational microlensing in 2006, analyses of six more events have been published by the end of 2018. In three events the lens was a star with two planets; four involved a binary star with a planet. Other possible triple lenses, such as triple stars or stars with a planet with a moon, are yet to be detected. The analysis of triple-lens events is hindered by the lack of understanding of the diversity of their caustics and critical curves. We present a method for identifying the full range of critical curves for a triple lens with a given combination of masses in an arbitrary spatial configuration. We compute their boundaries in parameter space, identify the critical-curve topologies in the partitioned regions, and evaluate their probabilities of occurrence. We demonstrate the analysis on three triple-lens models. For three equal masses the computed boundaries divide the parameter space into 39 regions yielding nine different critical-curve topologies. The other models include a binary star with a planet, and a hierarchical star--planet--moon combination of masses. Both have the same set of 11 topologies, including new ones with doubly nested critical-curve loops. The number of lensing regimes thus depends on the combination of masses -- unlike in the double lens, which has the same three regimes for any mass ratio. The presented approach is suitable for further investigations, such as studies of the changes occurring in nonstatic lens configurations due to orbital motion of the components or other parallax-type effects.


title: Protostellar Outflows at the EarliesT Stages (POETS). II. A possible radio synchrotron jet associated with the EGO G035.02+0.35

abstract: Centimeter continuum observations of protostellar jets have revealed the presence of knots of shocked gas where the flux density decreases with frequency. This spectrum is characteristic of nonthermal synchrotron radiation and implies the presence of both magnetic fields and relativistic electrons in protostellar jets. Here, we report on one of the few detections of nonthermal jet driven by a young massive star in the star-forming region G035.02$+$0.35. We made use of the NSF's Karl G. Jansky Very Large Array (VLA) to observe this region at C, Ku, and K bands with the A- and B-array configurations, and obtained sensitive radio continuum maps down to a rms of 10 $\mu$Jy beam$^{-1}$. These observations allow for a detailed spectral index analysis of the radio continuum emission in the region, which we interpret as a protostellar jet with a number of knots aligned with extended 4.5 $\mu$m emission. Two knots clearly emit nonthermal radiation and are found at similar distances, of approximately 10,000 au, each side of the central young star, from which they expand at velocities of hundreds km s$^{-1}$. We estimate both the mechanical force and the magnetic field associated with the radio jet, and infer a lower limit of $0.4\times10^{-4} $M$_{\odot}$ yr$^{-1}$ km s$^{-1}$ and values in the range $0.7-1.3 $mG, respectively.


title: Deep Neural Linear Bandits: Overcoming Catastrophic Forgetting through Likelihood Matching

abstract: We study the neural-linear bandit model for solving sequential decision-making problems with high dimensional side information. Neural-linear bandits leverage the representation power of deep neural networks and combine it with efficient exploration mechanisms, designed for linear contextual bandits, on top of the last hidden layer. Since the representation is being optimized during learning, information regarding exploration with "old" features is lost. Here, we propose the first limited memory neural-linear bandit that is resilient to this phenomenon, which we term catastrophic forgetting. We evaluate our method on a variety of real-world data sets, including regression, classification, and sentiment analysis, and observe that our algorithm is resilient to catastrophic forgetting and achieves superior performance.


title: Rainbow numbers of $[n]$ for $\sum_{i=1}^{k-1} x_i = x_k$

abstract: Consider the set $\{1,2,\dots,n\} = [n]$ and an equation $eq$. The rainbow number of $[n]$ for $eq$, denoted $\operatorname{rb}([n],eq)$, is the smallest number of colors such that for every exact $\operatorname{rb}([n], eq)$-coloring of $[n]$, there exists a solution to $eq$ with every member of the solution set assigned a distinct color. This paper focuses on linear equations and, in particular, establishes the rainbow number for the equations $\sum_{i=1}^{k-1} x_i = x_k$ for $k=3$ and $k=4$. The paper also establishes a general lower bound for $k \ge 5$.


title: Dilute magnetism in graphene

abstract: The phase diagram of graphene decorated with magnetic adatoms distributed either on a single sublattice, or evenly over the two sublattices, is computed for adatom concentrations as low as $\sim1\%$. Within the framework of the $s$-$d$ interaction, we take into account disorder effects due to the random positioning of the adatoms and/or to the thermal fluctuations in the direction of magnetic moments. Despite the presence of disorder, the magnetic phases are shown to be stable down to the lowest concentration accessed here. This result agrees with several experimental observations where adatom decorated graphene has been shown to have a magnetic response. In particular, the present theory provides a qualitative understanding for the results of Hwang et al. [Sci. Rep. 6, 21460 (2016)], where a ferromagnetic phase has been found below $\sim30\,\text{K}$ for graphene decorated with S-atoms.


title: Approximating power of machine-learning ansatz for quantum many-body states

abstract: An artificial neural network (ANN) with the restricted Boltzmann machine (RBM) architecture was recently proposed as a versatile variational quantum many-body wave function. In this work we provide physical insights into the performance of this ansatz. We uncover the connection between the structure of RBM and perturbation series, which explains the excellent precision achieved by RBM ansazt in certain simple models, demonstrated in the literature. Based on this relation, we improve the numerical algorithm to achieve better performance of RBM in cases where local minima complicate the convergence to the global one. We introduce other classes of variational wave-functions, which are also capable of reproducing the perturbative structure, and show that their performance is comparable to that of RBM. Furthermore, we study the performance of a few-layer RBM for approximating ground states of random, translationally-invariant models in 1d, as well as random matrix-product states (MPS). We find that the error in approximating such states exhibits a broad distribution, and is largely determined by the entanglement properties of the targeted state.


title: Boosting Standard Classification Architectures Through a Ranking Regularizer

abstract: We employ triplet loss as a feature embedding regularizer to boost classification performance. Standard architectures, like ResNet and Inception, are extended to support both losses with minimal hyper-parameter tuning. This promotes generality while fine-tuning pretrained networks. Triplet loss is a powerful surrogate for recently proposed embedding regularizers. Yet, it is avoided due to large batch-size requirement and high computational cost. Through our experiments, we re-assess these assumptions. During inference, our network supports both classification and embedding tasks without any computational overhead. Quantitative evaluation highlights a steady improvement on five fine-grained recognition datasets. Further evaluation on an imbalanced video dataset achieves significant improvement. Triplet loss brings feature embedding characteristics like nearest neighbor to classification models. Code available at \url{http://bit.ly/2LNYEqL}.


title: An Explicit Method for Schrieffer-Wolff Transformation

abstract: Schrieffer-Wolff transformation is a very important transformation in Quantum Many Body physics. Yet, there isn't an explicit method in the literature to calculate the generator of this unitary transformation directly from the hamiltonian. In this paper we present an explicit method to compute the generator of the Schrieffer-Wolff transformation in general and then calculate the same for two very important hamiltonians in condensed matter physics, the Single Impurity Anderson model (SIAM) and its lattice generalization , the Periodic Anderson Model (PAM). Our method can be used to carry out the transformation directly from the hamiltonian and the generator of the transformation can be obtained without any recourse to perturbative expansion or intuitive guesswork. In fact here, the full generator is calculated rather than its expansion to first two orders which is usually the case.


title: Verifiable Round-Robin Scheme for Smart Homes

abstract: Advances in sensing, networking, and actuation technologies have resulted in the IoT wave that is expected to revolutionize all aspects of modern society. This paper focuses on the new challenges of privacy that arise in IoT in the context of smart homes. Specifically, the paper focuses on preventing the user's privacy via inferences through channel and in-home device activities. We propose a method for securely scheduling the devices while decoupling the device and channels activities. The proposed solution avoids any attacks that may reveal the coordinated schedule of the devices, and hence, also, assures that inferences that may compromise individual's privacy are not leaked due to device and channel level activities. Our experiments also validate the proposed approach, and consequently, an adversary cannot infer device and channel activities by just observing the network traffic.


title: A Stable Combinatorial Particle Swarm Optimization for Scalable Feature Selection in Gene Expression Data

abstract: Evolutionary computation (EC) algorithms, such as discrete and multi-objective versions of particle swarm optimization (PSO), have been applied to solve the Feature selection (FS) problem, tackling the combinatorial explosion of search spaces that are peppered with local minima. Furthermore, high-dimensional FS problems such as finding a small set of biomarkers to make a diagnostic call add an additional challenge as such methods ability to pick out the most important features must remain unchanged in decision spaces of increasing dimensions and presence of irrelevant features. We developed a combinatorial PSO algorithm, called COMB-PSO, that scales up to high-dimensional gene expression data while still selecting the smallest subsets of genes that allow reliable classification of samples. In particular, COMB-PSO enhances the encoding, speed of convergence, control of divergence and diversity of the conventional PSO algorithm, balancing exploration and exploitation of the search space. Applying our approach on real gene expression data of different cancers, COMB-PSO finds gene sets of smallest size that allow a reliable classification of the underlying disease classes.


title: Atomistic simulation of nearly defect-free models of amorphous silicon: An information-based approach

abstract: We present an information-based total-energy optimization method to produce nearly defect-free structural models of amorphous silicon. Using geometrical, structural and topological information from disordered tetrahedral networks, we have shown that it is possible to generate structural configurations of amorphous silicon, which are superior than the models obtained from conventional reverse Monte Carlo and molecular-dynamics simulations. The new data-driven hybrid approach presented here is capable of producing atomistic models with structural and electronic properties which are on a par with those obtained from the modified Wooten-Winer-Weaire (WWW) models of amorphous silicon. Structural, electronic and thermodynamic properties of the hybrid models are compared with the best dynamical models obtained from using machine-intelligence-based potentials and efficient classical molecular-dynamics simulations, reported in the recent literature. We have shown that, together with the WWW models, our hybrid models represent one of the best structural models so far produced by total-energy-based Monte Carlo methods in conjunction with experimental diffraction data and a few structural constraints.


title: Learned Belief-Propagation Decoding with Simple Scaling and SNR Adaptation

abstract: We consider the weighted belief-propagation (WBP) decoder recently proposed by Nachmani et al. where different weights are introduced for each Tanner graph edge and optimized using machine learning techniques. Our focus is on simple-scaling models that use the same weights across certain edges to reduce the storage and computational burden. The main contribution is to show that simple scaling with few parameters often achieves the same gain as the full parameterization. Moreover, several training improvements for WBP are proposed. For example, it is shown that minimizing average binary cross-entropy is suboptimal in general in terms of bit error rate (BER) and a new "soft-BER" loss is proposed which can lead to better performance. We also investigate parameter adapter networks (PANs) that learn the relation between the signal-to-noise ratio and the WBP parameters. As an example, for the (32,16) Reed-Muller code with a highly redundant parity-check matrix, training a PAN with soft-BER loss gives near-maximum-likelihood performance assuming simple scaling with only three parameters.


title: Branko Gr\"unbaum in many dimensions

abstract: This survey article collects a few of my favorite open problems of Branko Gr\"{u}nbaum.


title: Petrophysical Property Estimation from Seismic Data Using Recurrent Neural Networks

abstract: Reservoir characterization involves the estimation petrophysical properties from well-log data and seismic data. Estimating such properties is a challenging task due to the non-linearity and heterogeneity of the subsurface. Various attempts have been made to estimate petrophysical properties using machine learning techniques such as feed-forward neural networks and support vector regression (SVR). Recent advances in machine learning have shown promising results for recurrent neural networks (RNN) in modeling complex sequential data such as videos and speech signals. In this work, we propose an algorithm for property estimation from seismic data using recurrent neural networks. An applications of the proposed workflow to estimate density and p-wave impedance using seismic data shows promising results compared to feed-forward neural networks.


title: AutoShuffleNet: Learning Permutation Matrices via an Exact Lipschitz Continuous Penalty in Deep Convolutional Neural Networks

abstract: ShuffleNet is a state-of-the-art light weight convolutional neural network architecture. Its basic operations include group, channel-wise convolution and channel shuffling. However, channel shuffling is manually designed empirically. Mathematically, shuffling is a multiplication by a permutation matrix. In this paper, we propose to automate channel shuffling by learning permutation matrices in network training. We introduce an exact Lipschitz continuous non-convex penalty so that it can be incorporated in the stochastic gradient descent to approximate permutation at high precision. Exact permutations are obtained by simple rounding at the end of training and are used in inference. The resulting network, referred to as AutoShuffleNet, achieved improved classification accuracies on CIFAR-10 and ImageNet data sets. In addition, we found experimentally that the standard convex relaxation of permutation matrices into stochastic matrices leads to poor performance. We prove theoretically the exactness (error bounds) in recovering permutation matrices when our penalty function is zero (very small). We present examples of permutation optimization through graph matching and two-layer neural network models where the loss functions are calculated in closed analytical form. In the examples, convex relaxation failed to capture permutations whereas our penalty succeeded.


title: Automatic Parallel Corpus Creation for Hindi-English News Translation Task

abstract: The parallel corpus for multilingual NLP tasks, deep learning applications like Statistical Machine Translation Systems is very important. The parallel corpus of Hindi-English language pair available for news translation task till date is of very limited size as per the requirement of the systems are concerned. In this work we have developed an automatic parallel corpus generation system prototype, which creates Hindi-English parallel corpus for news translation task. Further to verify the quality of generated parallel corpus we have experimented by taking various performance metrics and the results are quite interesting.


title: RADYNVERSION: Learning to Invert a Solar Flare Atmosphere with Invertible Neural Networks

abstract: During a solar flare, it is believed that reconnection takes place in the corona followed by fast energy transport to the chromosphere. The resulting intense heating strongly disturbs the chromospheric structure, and induces complex radiation hydrodynamic effects. Interpreting the physics of the flaring solar atmosphere is one of the most challenging tasks in solar physics. Here we present a novel deep learning approach, an invertible neural network, to understanding the chromospheric physics of a flaring solar atmosphere via the inversion of observed solar line profiles in H{\alpha} and Ca II {\lambda}8542. Our network is trained using flare simulations from the 1D radiation hydrodynamics code RADYN as the expected atmosphere and line profile. This model is then applied to single pixels from an observation of an M1.1 solar flare taken with SST/CRISP instrument just after the flare onset. The inverted atmospheres obtained from observations provide physical information on the electron number density, temperature and bulk velocity flow of the plasma throughout the solar atmosphere ranging from 0-10 Mm in height. The density and temperature profiles appear consistent with the expected atmospheric response, and the bulk plasma velocity provides the gradients needed to produce the broad spectral lines whilst also predicting the expected chromospheric evaporation from flare heating. We conclude that we have taught our novel algorithm the physics of a solar flare according to RADYN and that this can be confidently used for the analysis of flare data taken in these two wavelengths. This algorithm can also be adapted for a menagerie of inverse problems providing extremely fast ($\sim$10 {\mu}s) inversion samples.


title: High order concentrated non-negative matrix-exponential functions

abstract: Highly concentrated functions play an important role in many research fields including control system analysis and physics, and they turned out to be the key idea behind inverse Laplace transform methods as well. This paper uses the matrix-exponential family of functions to create highly concentrated functions, whose squared coefficient of variation (SCV) is very low. In the field of stochastic modeling, matrix-exponential functions have been used for decades. They have many advantages: they are easy to manipulate, always non-negative, and integrals involving matrix-exponential functions often have closed-form solutions. For the time being there is no symbolic construction available to obtain the most concentrated matrix-exponential functions, and the numerical optimization-based approach has many pitfalls, too. In this paper, we present a numerical optimization-based procedure to construct highly concentrated matrix-exponential functions. To make the objective function explicit and easy to evaluate we introduce and use a new representation called hyper-trigonometric representation. This representation makes it possible to achieve very low SCV.


title: Fair k-Center Clustering for Data Summarization

abstract: In data summarization we want to choose $k$ prototypes in order to summarize a data set. We study a setting where the data set comprises several demographic groups and we are restricted to choose $k_i$ prototypes belonging to group $i$. A common approach to the problem without the fairness constraint is to optimize a centroid-based clustering objective such as $k$-center. A natural extension then is to incorporate the fairness constraint into the clustering problem. Existing algorithms for doing so run in time super-quadratic in the size of the data set, which is in contrast to the standard $k$-center problem being approximable in linear time. In this paper, we resolve this gap by providing a simple approximation algorithm for the $k$-center problem under the fairness constraint with running time linear in the size of the data set and $k$. If the number of demographic groups is small, the approximation guarantee of our algorithm only incurs a constant-factor overhead.


title: Realization of digraphs in Abelian groups and its consequences

abstract: Let $\overrightarrow{G}$ be a directed graph with no component of orderless than~$3$, and let $\Gamma$ be a finite Abelian group such that $|\Gamma|\geq 4|V(\overrightarrow{G})|$ or if $|V(\overrightarrow{G})|$ is large enough with respect to an arbitrarily fixed $\varepsilon>0$ then $|\Gamma|\geq (1+\varepsilon)|V(\overrightarrow{G})|$. We show that there exists an injective mapping $\varphi$ from $V(\overrightarrow{G})$ to the group $\Gamma$ such that $\sum_{x\in V(C)}\varphi(x) = 0$ for every connected component $C$ of $\overrightarrow{G}$, where $0$ is the identity element of $\Gamma$. Moreover we show some applications of this result to group distance magic labelings.


title: Real-time Scene Segmentation Using a Light Deep Neural Network Architecture for Autonomous Robot Navigation on Construction Sites

abstract: Camera-equipped unmanned vehicles (UVs) have received a lot of attention in data collection for construction monitoring applications. To develop an autonomous platform, the UV should be able to process multiple modules (e.g., context-awareness, control, localization, and mapping) on an embedded platform. Pixel-wise semantic segmentation provides a UV with the ability to be contextually aware of its surrounding environment. However, in the case of mobile robotic systems with limited computing resources, the large size of the segmentation model and high memory usage requires high computing resources, which a major challenge for mobile UVs (e.g., a small-scale vehicle with limited payload and space). To overcome this challenge, this paper presents a light and efficient deep neural network architecture to run on an embedded platform in real-time. The proposed model segments navigable space on an image sequence (i.e., a video stream), which is essential for an autonomous vehicle that is based on machine vision. The results demonstrate the performance efficiency of the proposed architecture compared to the existing models and suggest possible improvements that could make the model even more efficient, which is necessary for the future development of the autonomous robotics systems.


title: A comparison theorem for nonsmooth nonlinear operators

abstract: We prove a comparison theorem for super- and sub-solutions with non-vanishing gradients to semilinear PDEs provided a nonlinearity $f$ is $L^p$ function with $p > 1$. The proof is based on a strong maximum principle for solutions of divergence type elliptic equations with VMO leading coefficients and with lower order coefficients from a Kato class. An application to estimation of periodic water waves profiles is given.


title: Dynamical systems and nonlinear transient rheology of the far-from-equilibrium Bjorken flow

abstract: In relativistic kinetic theory, the one-particle distribution function is approximated by an asymptotic perturbative power series in Knudsen number which is divergent. For the Bjorken flow, we expand the distribution function in terms of its moments and study their nonlinear evolution equations. The resulting coupled dynamical system can be solved for each moment consistently using a multi-parameter transseries which makes the constitutive relations inherit the same structure. A new non-perturbative dynamical renormalization scheme is born out of this formalism that goes beyond the linear response theory. We show that there is a Lyapunov function, aka dynamical potential, which is, in general, a function of the moments and time satisfying Lyapunov stability conditions along RG flows connected to the asymptotic hydrodynamic fixed point. As a result, the transport coefficients get dynamically renormalized at every order in the time-dependent perturbative expansion by receiving non-perturbative corrections present in the transseries. The connection between the integration constants and the UV data is discussed using the language of dynamical systems. Furthermore, we show that the first dissipative correction in the Knudsen number to the distribution function is not only determined by the known effective shear viscous term but also a new high energy non-hydrodynamic mode. It is demonstrated that the survival of this new mode is intrinsically related to the nonlinear mode-to-mode coupling with the shear viscous term. Finally, we comment on some possible phenomenological applications of the proposed non-hydrodynamic transport theory.


title: Structure and superconductivity in the binary Re$_{1-x}$Mo$_x$ alloys

abstract: The binary Re$_{1-x}$Mo$_x$ alloys, known to cover the full range of solid solutions, were successfully synthesized and their crystal structures and physical properties investigated via powder x-ray diffraction, electrical resistivity, magnetic susceptibility, and heat capacity. By varying the Re/Mo ratio we explore the full Re$_{1-x}$Mo$_x$ binary phase diagram, in all its four different solid phases: hcp-Mg ($P6_3/mmc$), $\alpha$-Mn ($I\overline{4}3m$), $\beta$-CrFe ($P4_2/mnm$), and bcc-W ($Im\overline{3}m$), of which the second is non-centrosymmetric with the rest being centrosymmetric. All Re$_{1-x}$Mo$_x$ alloys are superconductors, whose critical temperatures exhibit a peculiar phase diagram, characterized by three different superconducting regions. In most alloys the $T_c$ is almost an order of magnitude higher than in pure Re and Mo. Low-temperature electronic specific-heat data evidence a fully-gapped superconducting state, whose enhanced gap magnitude and specific-heat discontinuity suggest a moderately strong electron-phonon coupling across the series. Considering that several $\alpha$-Mn-type Re$T$ alloys ($T$ = transition metal) show time-reversal symmetry breaking (TRSB) in the superconducting state, while TRS is preserved in the isostructural Mg$_{10}$Ir$_{19}$B$_{16}$ or Nb$_{0.5}$Os$_{0.5}$, the Re$_{1-x}$Mo$_x$ alloys represent another suitable system for studying the interplay of space-inversion, gauge, and time-reversal symmetries in future experiments expected to probe TRSB in the Re$T$ family.


title: A BERT Baseline for the Natural Questions

abstract: This technical note describes a new baseline for the Natural Questions. Our model is based on BERT and reduces the gap between the model F1 scores reported in the original dataset paper and the human upper bound by 30% and 50% relative for the long and short answer tasks respectively. This baseline has been submitted to the official NQ leaderboard at ai.google.com/research/NaturalQuestions. Code, preprocessed data and pretrained model are available at https://github.com/google-research/language/tree/master/language/question_answering/bert_joint.


title: Stochastic modeling of phenotypic switching and chemoresistance in cancer cell populations

abstract: Phenotypic heterogeneity in cancer cells is widely observed and is often linked to drug resistance. In several cases, such heterogeneity in drug sensitivity of tumors is driven by stochastic and reversible acquisition of a drug tolerant phenotype by individual cells even in an isogenic population. Accumulating evidence further suggests that cell-fate transitions such as the epithelial to mesenchymal transition (EMT) are associated with drug resistance. In this study, we analyze stochastic models of phenotypic switching to provide a framework for analyzing cell-fate transitions such as EMT as a source of phenotypic variability in drug sensitivity. Motivated by our cell-culture based experimental observations connecting phenotypic switching in EMT and drug resistance, we analyze a coarse-grained model of phenotypic switching between two states in the presence of cytotoxic stress from chemotherapy. We derive analytical results for time-dependent probability distributions that provide insights into the rates of phenotypic switching and characterize initial phenotypic heterogeneity of cancer cells. The results obtained can also shed light on fundamental questions relating to adaptation and selection scenarios in tumor response to cytotoxic therapy.


title: Evolution Boussinesq model with nonmonotone friction and heat flux boundary conditions

abstract: In this paper we prove the existence and regularity of a solution to a two-dimensional system of evolutionary hemivariational inequalities which describes the Boussinesq model with nonmonotone friction and heat flux. We use the time retardation and regularization technique, combined with a regularized Galerkin method, and recent results from the theory of hemivariational inequalities.


title: State-of-the-Art Perovskite Solar Cells Benefit from Photon Recycling at Maximum Power Point

abstract: Photon recycling is required for a solar cell to achieve an open-circuit voltage ($V_{OC}$) and power conversion efficiency (PCE) approaching the Shockley-Queisser theoretical limit. In metal halide perovskite solar cells, the achievable performance gains from photon recycling remain uncertain due to high variability in perovskite material quality and the non-radiative recombination rate ($k_{1}$). In this work, we study state-of-the-art $\textrm{Cs}_{0.05}(\textrm{MA}_{0.17}\textrm{FA}_{0.83})_{0.95}\textrm{Pb}(\textrm{I}_{0.83}\textrm{Br}_{0.17})_{3}$ films and analyze the impact of varying non-radiative recombination rates on photon recycling and device performance. Importantly, we predict the impact of photon recycling at the maximum power point (MPP), demonstrating an absolute PCE increase of up to 2.0% in the radiative limit, primarily due to a 77 mV increase in $V_{MPP}$. Even with finite non-radiative recombination, benefits from photon recycling can be achieved when non-radiative lifetimes and external LED electroluminescence efficiencies measured at open-circuit, $Q_{e}^{LED}(\textrm{V}_{OC})$, exceed 2 $\mu$s and 10%, respectively. This analysis clarifies the opportunity to fully exploit photon recycling to push the real-world performance of perovskite solar cells toward theoretical limits.


title: Towards Real-Time Respiratory Motion Prediction based on Long Short-Term Memory Neural Networks

abstract: Radiation therapy of thoracic and abdominal tumors requires incorporating the respiratory motion into treatments. To precisely account for the patient respiratory motions and predict the respiratory signals, a generalized model for predictions of different types of respiratory motions is desired. The aim of this study is to explore the feasibility of developing a Long Short-Term Memory (LSTM)-based model for the respiratory signal prediction. To achieve that, 1703 sets of Real-Time Position Management data were collected from retrospective studies across three clinical institutions. These datasets were separated as the training, internal validity and external validity groups. 1187 datasets were used for model development and the remaining 516 datasets were used to test the generality power of the model. Furthermore, an exhaustive grid search was implemented to find the optimal hyper-parameters of the LSTM model. The hyper-parameters are the number of LSTM layers, the number of hidden units, the optimizer, the learning rate, the number of epochs, and the length of time lags. Our model achieved superior accuracy over conventional artificial neural network models: with a prediction window of 500ms, the LSTM model achieved an average relative Mean Absolute Error (MAE) of 0.037, an average Root Mean Square Error (RMSE) of 0.048, and a Maximum Error (ME) of 1.687 in the internal validity data, and an average relative MAE of 0.112, an average RMSE of 0.139 and an ME of 1.811 in the external validity data. Compared to the LSTM model trained with default hyper-parameters, the MAE of the optimized model results decreased by 20%, indicating the importance of tuning the hyper-parameters of LSTM models to obtain superior accuracy. This study demonstrates the potential of deep LSTM models for the respiratory signal prediction and illustrates the impacts of major hyper-parameters in LSTM models.


title: Dolha - an Efficient and Exact Data Structure for Streaming Graphs

abstract: A streaming graph is a graph formed by a sequence of incoming edges with time stamps. Unlike static graphs, the streaming graph is highly dynamic and time related. In the real world, the high volume and velocity streaming graphs such as internet traffic data, social network communication data and financial transfer data are bringing challenges to the classic graph data structures. We present a new data structure: double orthogonal list in hash table (Dolha) which is a high speed and high memory efficiency graph structure applicable to streaming graph. Dolha has constant time cost for single edge and near linear space cost that we can contain billions of edges information in memory size and process an incoming edge in nanoseconds. Dolha also has linear time cost for neighborhood queries, which allow it to support most algorithms in graphs without extra cost. We also present a persistent structure based on Dolha that has the ability to handle the sliding window update and time related queries.


title: Sparkling nights and very hot days on WASP-18b: the formation of clouds and the emergence of an ionosphere

abstract: WASP-18b is an utra-hot Jupiter with a temperature difference of upto 2500K between day and night. Such giant planets begin to emerge as planetary laboratory for understanding cloud formation and gas chemistry in well-tested parameter regimes in order to better understand planetary mass loss and for linking observed element ratios to planet formation and evolution. We aim to understand where clouds form, their interaction with the gas phase chemistry through depletion and enrichment, the ionisation of the atmospheric gas and the possible emergence of an ionosphere on ultra-hot Jupiters. We utilize 1D profiles from a 3D atmosphere simulations for WASP-18b as input for kinetic cloud formation and gas-phase chemical equilibrium calculations. We solve our kinetic cloud formation model for these 1D profiles that sample the atmosphere of WASP-18b at 16 different locations along the equator and in the mid-latitudes and derive consistently the gas-phase composition. The dayside of WASP-18b emerges as completely cloud-free due to the very high atmospheric temperatures. In contrast, the nightside is covered in geometrically extended and chemically heterogeneous clouds with disperse particle size distributions. The atmospheric C/O increases to $>0.7$ and the enrichment of the atmospheric gas with cloud particles is $\rho_{\rm d}/\rho_{\rm gas}>10^{-3}$. The clouds that form at the limbs appear located farther inside the atmosphere and they are the least extended. Not all day-night terminator regions form clouds. The gas-phase is dominated by H$_2$, CO, SiO, H$_2$O, H$_2$S, CH$_4$, SiS. In addition, the dayside has a substantial degree of ionisation due to ions like Na$^+$, K$^+$, Ca$^+$, Fe$^+$. Al$^+$ and Ti$^+$ are the most abundant of their element classes. We find that WASP-18b, as one example for ultra-hot Jupiters, develops an ionosphere on the dayside.


title: Gibbs posterior convergence and the thermodynamic formalism

abstract: In this paper we consider a Bayesian framework for making inferences about dynamical systems from ergodic observations. The proposed Bayesian procedure is based on the Gibbs posterior, a decision theoretic generalization of standard Bayesian inference. We place a prior over a model class consisting of a parametrized family of Gibbs measures on a mixing shift of finite type. This model class generalizes (hidden) Markov chain models by allowing for long range dependencies, including Markov chains of arbitrarily large orders. We characterize the asymptotic behavior of the Gibbs posterior distribution on the parameter space as the number of observations tends to infinity. In particular, we define a limiting variational problem over the space of joinings of the model system with the observed system, and we show that the Gibbs posterior distributions concentrate around the solution set of this variational problem. In the case of properly specified models our convergence results may be used to establish posterior consistency. This work establishes tight connections between Gibbs posterior inference and the thermodynamic formalism, which may inspire new proof techniques in the study of Bayesian posterior consistency for dependent processes.


title: 3d "chiral" Kutasov-Schwimmer duality

abstract: We propose a "chiral" version of the Kutasov-Schwimmer duality in a 3d $\mathcal{N}=2$ $SU(N)$ gauge theory with $F$ fundamental matters, $\bar{F}$ anti-fundamental matters and an adjoint matter $X$ with a tree-level superpotential $W= \mathrm{tr} \, X^{k+1}$. The theory exhibits a rich structure of the baryonic and (dressed) Coulomb branch operators. At first sight, the duality seems bad due to the mismatch of the anti-baryonic branch in the moduli space of vacua. The duality well works by realizing that the anti-baryonic operators are identified with some of the dressed Coulomb branch coordinates under the proposed duality. This generalizes the $SU(N)$ ``chiral'' duality with (anti-)fundamental matters, which we previously proposed.


title: A class of dynamic frictional contact problems governed by a system of hemivariational inequalities in thermoviscoelasticity

abstract: In this paper we prove the existence and uniqueness of the weak solution for a dynamic thermoviscoelastic problem which describes frictional contact between a body and a foundation. We employ the nonlinear constitutive viscoelastic law with a long-term memory, which include the thermal effects and consider the general nonmonotone and multivalued subdifferential boundary conditions for the contact, friction and heat flux. The model consists of the system of the hemivariational inequality of for the displacement and the parabolic hemivariational inequality for the temperature. The existence of solutions is proved by using recent results from the theory of hemivariational inequalities and a fixed point argument.


title: Ablation Studies in Artificial Neural Networks

abstract: Ablation studies have been widely used in the field of neuroscience to tackle complex biological systems such as the extensively studied Drosophila central nervous system, the vertebrate brain and more interestingly and most delicately, the human brain. In the past, these kinds of studies were utilized to uncover structure and organization in the brain, i.e. a mapping of features inherent to external stimuli onto different areas of the neocortex. considering the growth in size and complexity of state-of-the-art artificial neural networks (ANNs) and the corresponding growth in complexity of the tasks that are tackled by these networks, the question arises whether ablation studies may be used to investigate these networks for a similar organization of their inner representations. In this paper, we address this question and performed two ablation studies in two fundamentally different ANNs to investigate their inner representations of two well-known benchmark datasets from the computer vision domain. We found that features distinct to the local and global structure of the data are selectively represented in specific parts of the network. Furthermore, some of these representations are redundant, awarding the network a certain robustness to structural damages. We further determined the importance of specific parts of the network for the classification task solely based on the weight structure of single units. Finally, we examined the ability of damaged networks to recover from the consequences of ablations by means of recovery training. We argue that ablations studies are a feasible method to investigate knowledge representations in ANNs and are especially helpful to examine a networks robustness to structural damages, a feature of ANNs that will become increasingly important for future safety-critical applications.


title: Local cohomology of binomial edge ideals and their generic ideals

abstract: We provide a Hochster type formula for the local cohomology modules of binomial edge ideals. As a consequence we obtain a simple criterion for the Cohen-Macaulayness of these ideals and we describe their Castelnuovo-Mumford regularity and their Hilbert series. We also prove a conjecture of Conca, De Negri and Gorla relating the graded components of the local cohomology modules of binomial edge ideals and their generic initial ideals.


title: Approximation by sequence of operators including Dunkl Appell polynomials

abstract: In this article, we give a sequence of operators for producing an approximation result. The relation between the rate of approximation of sequence operators including Dunkl variant of exponential function with first and second-order modulus of continuity are shown. A specific application of sequence of operators which include Gould-Hopper type polynomials is constructed.


title: A Hot Ultraviolet Flare on the M Dwarf Star GJ 674

abstract: As part of the Mega MUSCLES Hubble Space Telescope (HST) Treasury program, we obtained time-series ultraviolet spectroscopy of the M2.5V star, GJ~674. During the FUV monitoring observations, the target exhibited several small flares and one large flare (E_FUV = 10^{30.75} ergs) that persisted over the entirety of a HST orbit and had an equivalent duration >30,000 sec, comparable to the highest relative amplitude event previously recorded in the FUV. The flare spectrum exhibited enhanced line emission from chromospheric, transition region, and coronal transitions and a blue FUV continuum with an unprecedented color temperature of T_c ~ 40,000+/-10,000 K. In this paper, we compare the flare FUV continuum emission with parameterizations of radiative hydrodynamic model atmospheres of M star flares. We find that the observed flare continuum can be reproduced using flare models but only with the ad hoc addition of hot, dense emitting component. This observation demonstrates that flares with hot FUV continuum temperatures and significant EUV/FUV energy deposition will continue to be of importance to exoplanet atmospheric chemistry and heating even as the host M dwarfs age beyond their most active evolutionary phases.


title: Krickeberg mixing for Z extensions of Gibbs Markov semiflows

abstract: We obtain Krickeberg mixing for a class of Z extensions of Gibbs Markov semiflows with roof function and displacement function not in L2, where previous methods have not been employed. This is done via a 'smooth tail' estimate for the isomorphic suspension flow.


title: Learning Independently-Obtainable Reward Functions

abstract: We present a novel method for learning a set of disentangled reward functions that sum to the original environment reward and are constrained to be independently obtainable. We define independent obtainability in terms of value functions with respect to obtaining one learned reward while pursuing another learned reward. Empirically, we illustrate that our method can learn meaningful reward decompositions in a variety of domains and that these decompositions exhibit some form of generalization performance when the environment's reward is modified. Theoretically, we derive results about the effect of maximizing our method's objective on the resulting reward functions and their corresponding optimal policies.


title: Adaptive Optimal Control of Linear Periodic Systems: An Off-Policy Value Iteration Approach

abstract: This paper studies the infinite-horizon adaptive optimal control of continuous-time linear periodic (CTLP) systems. A novel value iteration (VI) based off-policy ADP algorithm is proposed for a general class of CTLP systems, so that approximate optimal solutions can be obtained directly from the collected data, without the exact knowledge of system dynamics. Under mild conditions, the proofs on uniform convergence of the proposed algorithm to the optimal solutions are given for both the model-based and model-free cases. The VI-based ADP algorithm is able to find suboptimal controllers without assuming the knowledge of an initial stabilizing controller. Application to the optimal control of a triple inverted pendulum subjected to a periodically varying load demonstrates the feasibility and effectiveness of the proposed method.


title: Decoupling feature extraction from policy learning: assessing benefits of state representation learning in goal based robotics

abstract: Scaling end-to-end reinforcement learning to control real robots from vision presents a series of challenges, in particular in terms of sample efficiency. Against end-to-end learning, state representation learning can help learn a compact, efficient and relevant representation of states that speeds up policy learning, reducing the number of samples needed, and that is easier to interpret. We evaluate several state representation learning methods on goal based robotics tasks and propose a new unsupervised model that stacks representations and combines strengths of several of these approaches. This method encodes all the relevant features, performs on par or better than end-to-end learning with better sample efficiency, and is robust to hyper-parameters change.


title: Learning agile and dynamic motor skills for legged robots

abstract: Legged robots pose one of the greatest challenges in robotics. Dynamic and agile maneuvers of animals cannot be imitated by existing methods that are crafted by humans. A compelling alternative is reinforcement learning, which requires minimal craftsmanship and promotes the natural evolution of a control policy. However, so far, reinforcement learning research for legged robots is mainly limited to simulation, and only few and comparably simple examples have been deployed on real systems. The primary reason is that training with real robots, particularly with dynamically balancing systems, is complicated and expensive. In the present work, we introduce a method for training a neural network policy in simulation and transferring it to a state-of-the-art legged system, thereby leveraging fast, automated, and cost-effective data generation schemes. The approach is applied to the ANYmal robot, a sophisticated medium-dog-sized quadrupedal system. Using policies trained in simulation, the quadrupedal machine achieves locomotion skills that go beyond what had been achieved with prior methods: ANYmal is capable of precisely and energy-efficiently following high-level body velocity commands, running faster than before, and recovering from falling even in complex configurations.


title: An Identity for Vertically Aligned Entries in Pascal's Triangle

abstract: The classic way to write down Pascal's triangle leads to entries in alternating rows being vertically aligned. In this paper, we prove a linear dependence on vertically aligned entries in Pascal's triangle. Furthermore, we give an application of this dependence to morphisms between hyperelliptic curves.


title: The Assistive Multi-Armed Bandit

abstract: Learning preferences implicit in the choices humans make is a well studied problem in both economics and computer science. However, most work makes the assumption that humans are acting (noisily) optimally with respect to their preferences. Such approaches can fail when people are themselves learning about what they want. In this work, we introduce the assistive multi-armed bandit, where a robot assists a human playing a bandit task to maximize cumulative reward. In this problem, the human does not know the reward function but can learn it through the rewards received from arm pulls; the robot only observes which arms the human pulls but not the reward associated with each pull. We offer sufficient and necessary conditions for successfully assisting the human in this framework. Surprisingly, better human performance in isolation does not necessarily lead to better performance when assisted by the robot: a human policy can do better by effectively communicating its observed rewards to the robot. We conduct proof-of-concept experiments that support these results. We see this work as contributing towards a theory behind algorithms for human-robot interaction.


title: Small ball probability for the condition number of random matrices

abstract: Let $A$ be an $n\times n$ random matrix with i.i.d. entries of zero mean, unit variance and a bounded subgaussian moment. We show that the condition number $s_{\max}(A)/s_{\min}(A)$ satisfies the small ball probability estimate $${\mathbb P}\big\{s_{\max}(A)/s_{\min}(A)\leq n/t\big\}\leq 2\exp(-c t^2),\quad t\geq 1,$$ where $c>0$ may only depend on the subgaussian moment. Although the estimate can be obtained as a combination of known results and techniques, it was not noticed in the literature before. As a key step of the proof, we apply estimates for the singular values of $A$, ${\mathbb P}\big\{s_{n-k+1}(A)\leq ck/\sqrt{n}\big\}\leq 2 \exp(-c k^2), \quad 1\leq k\leq n,$ obtained (under some additional assumptions) by Nguyen.


title: Edges control clustering in levitated granular matter

abstract: The properties of small clusters depend dramatically on the interactions between their constituent particles. However, it remains challenging to design and tune the interactions between macroscopic particles, such as in a granular material. Here, we use acoustic levitation to trap macroscopic grains and induce forces between them. Our main results show that particles levitated in an acoustic field prefer to make contact along sharp edges. The radius of curvature of the edges directly controls the magnitude of these forces. These highly directional interactions, combined with local contact forces, give rise to a diverse array of cluster shapes. Our results open up new possibilities for the design of specific forces between macroscopic particles, directing their assembly, and actuating their motion.


title: On the Monoidal Center of Deligne's Category Rep(S_t)

abstract: We explicitly compute a monoidal subcategory of the monoidal center of Deligne's interpolation category Rep(S_t), for t not necessarily a natural number, and we show that this subcategory is a ribbon category. For t=n, a natural number, there exists a functor onto the braided monoidal category of modules over the Drinfeld double of S_n which is essentially surjective and full. Hence the new ribbon category interpolates the categories of crossed modules over the symmetric groups. As an application, we obtain invariants of framed ribbon links which are polynomials in the interpolating variable t. These polynomials interpolate untwisted Dijkgraaf-Witten invariants of the symmetric groups.


title: Is Pretraining Necessary for Hyperspectral Image Classification?

abstract: We address two questions for training a convolutional neural network (CNN) for hyperspectral image classification: i) is it possible to build a pre-trained network? and ii) is the pre-training effective in furthering the performance? To answer the first question, we have devised an approach that pre-trains a network on multiple source datasets that differ in their hyperspectral characteristics and fine-tunes on a target dataset. This approach effectively resolves the architectural issue that arises when transferring meaningful information between the source and the target networks. To answer the second question, we carried out several ablation experiments. Based on the experimental results, a network trained from scratch performs as good as a network fine-tuned from a pre-trained network. However, we observed that pre-training the network has its own advantage in achieving better performances when deeper networks are required.


title: Projected Stein Variational Newton: A Fast and Scalable Bayesian Inference Method in High Dimensions

abstract: We propose a fast and scalable variational method for Bayesian inference in high-dimensional parameter space, which we call projected Stein variational Newton (pSVN) method. We exploit the intrinsic low-dimensional geometric structure of the posterior distribution in the high-dimensional parameter space via its Hessian (of the log posterior) operator and perform a parallel update of the parameter samples projected into a low-dimensional subspace by an SVN method. The subspace is adaptively constructed using the eigenvectors of the averaged Hessian at the current samples. We demonstrate fast convergence of the proposed method and its scalability with respect to the number of parameters, samples, and processor cores.


title: Multimodal phase velocity-frequency dispersion images using different MASW transformation techniques

abstract: Three different transformation techniques, namely, (i)w-c, (ii) w-k and (iii) tau-p, has been employed for generating multimodal dispersion images on the basis of multi-channel analysis of surface waves (MASW) data recorded in distance-time domain; here w= circular frequency, c = phase velocity, tau = time intercept, p = phase slowness (1/c) and k = wavenumber. All the three methods have been first clearly described. The results from these three different transforms have been examined by using synthetic as well as field data obtained from field tests using 48 geophones. The effect of sensor spread length (X) and geophone numbers (M) on multimodal dispersion images were examined. The solutions from these three transforms were found to match generally well with each other. The w-c transform has been noted to provide the most clarity since it does not require either high sampling rate as normally needed for the tau-p method or inclusion of the zero padding of the data in a distance domain for the w-k approach. The paper will be useful since it not only describes the methods, but it brings out simultaneously their merits in implementation and for generating the dispersion images.


title: Intelligent Sensor Based Bayesian Neural Network for Combined Parameters and States Estimation of a Brushed DC Motor

abstract: The objective of this paper is to develop an Artificial Neural Network (ANN) model to estimate simultaneously, parameters and state of a brushed DC machine. The proposed ANN estimator is novel in the sense that his estimates simultaneously temperature, speed and rotor resistance based only on the measurement of the voltage and current inputs. Many types of ANN estimators have been designed by a lot of researchers during the last two decades. Each type is designed for a specific application. The thermal behavior of the motor is very slow, which leads to large amounts of data sets. The standard ANN use often Multi-Layer Perceptron (MLP) with Levenberg-Marquardt Backpropagation (LMBP), among the limits of LMBP in the case of large number of data, so the use of MLP based on LMBP is no longer valid in our case. As solution, we propose the use of Cascade-Forward Neural Network (CFNN) based Bayesian Regulation backpropagation (BRBP). To test our estimator robustness a random white-Gaussian noise has been added to the sets. The proposed estimator is in our viewpoint accurate and robust.


title: An Identity for Second Order Sequences Obeying the Same Recurrence Relation

abstract: We derive an identity connecting any two second-order linear recurrence sequences having the same recurrence relation but whose initial terms may be different. Binomial and ordinary summation identities arising from the identity are developed. Illustrative examples are drawn from Fibonacci, Fibonacci-Lucas, Pell, Pell-Lucas, Jacobsthal and Jacobsthal-Lucas sequences and their generalizations. Our new results subsume previously known identities.


title: New nonasymptotic convergence rates of stochastic proximal pointalgorithm for convex optimization problems

abstract: Large sectors of the recent optimization literature focused in the last decade on the development of optimal stochastic first order schemes for constrained convex models under progressively relaxed assumptions. Stochastic proximal point is an iterative scheme born from the adaptation of proximal point algorithm to noisy stochastic optimization, with a resulting iteration related to stochastic alternating projections. Inspired by the scalability of alternating projection methods, we start from the (linear) regularity assumption, typically used in convex feasiblity problems to guarantee the linear convergence of stochastic alternating projection methods, and analyze a general weak linear regularity condition which facilitates convergence rate boosts in stochastic proximal point schemes. Our applications include many non-strongly convex functions classes often used in machine learning and statistics. Moreover, under weak linear regularity assumption we guarantee $\mathcal{O}\left(\frac{1}{k}\right)$ convergence rate for SPP, in terms of the distance to the optimal set, using only projections onto a simple component set. Linear convergence is obtained for interpolation setting, when the optimal set of the expected cost is included into the optimal sets of each functional component.


title: Localized bright luminescence of indirect excitons and trions in MoSe$_2$/WSe$_2$ van der Waals heterostructure

abstract: Indirect excitons (IX) in semiconductor heterostructures are bosons, which can cool below the temperature of quantum degeneracy and can be effectively controlled by voltage and light. IX quantum Bose gases and IX devices were explored in GaAs heterostructures where an IX range of existence is limited to low temperatures due to low IX binding energies. IXs in van der Waals transition-metal dichalcogenide (TMD) heterostructures are characterized by large binding energies giving the opportunity for exploring excitonic quantum gases and for creating excitonic devices at high temperatures. TMD heterostructures also offer a new platform for studying single-exciton phenomena and few-particle complexes. In this work, we present studies of IXs in MoSe$_2$/WSe$_2$ heterostructures and report on two IX luminescence lines whose energy splitting and temperature dependence identify them as neutral and charged IXs. The experimentally found binding energy of the indirect charged excitons, i.e. indirect trions, is close to the calculated binding energy of 28 meV for negative indirect trions in TMD heterostructures [Deilmann, Thygesen, Nano Lett. 18, 1460 (2018)]. We also report on the realization of IXs with a luminescence linewidth reaching 4~meV at low temperatures. An enhancement of IX luminescence intensity and the narrow linewidth are observed in localized spots.


title: Fairness risk measures

abstract: Ensuring that classifiers are non-discriminatory or fair with respect to a sensitive feature (e.g., race or gender) is a topical problem. Progress in this task requires fixing a definition of fairness, and there have been several proposals in this regard over the past few years. Several of these, however, assume either binary sensitive features (thus precluding categorical or real-valued sensitive groups), or result in non-convex objectives (thus adversely affecting the optimisation landscape). In this paper, we propose a new definition of fairness that generalises some existing proposals, while allowing for generic sensitive features and resulting in a convex objective. The key idea is to enforce that the expected losses (or risks) across each subgroup induced by the sensitive feature are commensurate. We show how this relates to the rich literature on risk measures from mathematical finance. As a special case, this leads to a new convex fairness-aware objective based on minimising the conditional value at risk (CVaR).


title: HRDBMS: Combining the Best of Modern and Traditional Relational Databases

abstract: HRDBMS is a novel distributed relational database that uses a hybrid model combining the best of traditional distributed relational databases and Big Data analytics platforms such as Hive. This allows HRDBMS to leverage years worth of research regarding query optimization, while also taking advantage of the scalability of Big Data platforms. The system uses an execution framework that is tailored for relational processing, thus addressing some of the performance challenges of running SQL on top of platforms such as MapReduce and Spark. These include excessive materialization of intermediate results, lack of a global cost-based optimization, unnecessary sorting, lack of index support, no statistics, no support for DML and ACID, and excessive communication caused by the rigid communication patterns enforced by these platforms.


title: Parylene Based Memristive Devices with Multilevel Resistive Switching for Neuromorphic Applications

abstract: In this paper, the resistive switching and neuromorphic behavior of memristive devices based on parylene, a polymer both low-cost and safe for the human body, is comprehensively studied. The Metal/Parylene/ITO sandwich structures were prepared by means of the standard gas phase surface polymerization method with different top active metal electrodes (Ag, Al, Cu or Ti of about 500 nm thickness). These organic memristive devices exhibit excellent performance: low switching voltage (down to 1 V), large OFF/ON resistance ratio (about 10^3), retention (> 10^4 s) and high multilevel resistance switching (at least 16 stable resistive states in the case of Cu electrodes). We have experimentally shown that parylene-based memristive elements can be trained by a biologically inspired spike-timing-dependent plasticity (STDP) mechanism. The obtained results have been used to implement a simple neuromorphic network model of classical conditioning. The described advantages allow considering parylene-based organic memristors as prospective devices for hardware realization of spiking artificial neuron networks capable of supervised and unsupervised learning and suitable for biomedical applications.


title: Guarantees for Spectral Clustering with Fairness Constraints

abstract: Given the widespread popularity of spectral clustering (SC) for partitioning graph data, we study a version of constrained SC in which we try to incorporate the fairness notion proposed by Chierichetti et al. (2017). According to this notion, a clustering is fair if every demographic group is approximately proportionally represented in each cluster. To this end, we develop variants of both normalized and unnormalized constrained SC and show that they help find fairer clusterings on both synthetic and real data. We also provide a rigorous theoretical analysis of our algorithms on a natural variant of the stochastic block model, where $h$ groups have strong inter-group connectivity, but also exhibit a "natural" clustering structure which is fair. We prove that our algorithms can recover this fair clustering with high probability.


title: SAGA with Arbitrary Sampling

abstract: We study the problem of minimizing the average of a very large number of smooth functions, which is of key importance in training supervised learning models. One of the most celebrated methods in this context is the SAGA algorithm. Despite years of research on the topic, a general-purpose version of SAGA---one that would include arbitrary importance sampling and minibatching schemes---does not exist. We remedy this situation and propose a general and flexible variant of SAGA following the {\em arbitrary sampling} paradigm. We perform an iteration complexity analysis of the method, largely possible due to the construction of new stochastic Lyapunov functions. We establish linear convergence rates in the smooth and strongly convex regime, and under a quadratic functional growth condition (i.e., in a regime not assuming strong convexity). Our rates match those of the primal-dual method Quartz for which an arbitrary sampling analysis is available, which makes a significant step towards closing the gap in our understanding of complexity of primal and dual methods for finite sum problems.


title: Python script used as simulator for the teaching of electric field in electromagnetism course

abstract: We present this work like software tool developed in Python, based on a methodology to obtain the electric field produced by n charges. The tool was developed and implemented in courses of electromagnetism and laboratory in three institutions of higher education. The aim for this work is to incorporate information and communication technologies (ICTs) at the university, in accordance with the programs promoted by the Colombian Ministry of Education. We wanted to connect the students with sensitives experiences of the physical phenomena that allow them to improve their experience of learning of subjects traditionally studied through the board course.


title: Two-dimensional electron gas in a metal/amorphous oxide interface with spin orbit interaction

abstract: The formation of novel two-dimensional electron gas (2DEG) with high mobility in metal/amorphous interfaces has motivated an ongoing debate regarding the formation and novel characteristics of these 2DEGs. Here we report an optical study, based on infrared spectroscopic ellipsometry, of nonmagnetic metal and amorphous semiconducting oxide (Cu/Bi$_2$O$_3$) interfaces that confirms the formation of a 2DEG with spin orbit coupling (SOC). The 2DEG optical response was simulated with a uniaxial diagonal dielectric tensor within a sub-nanometer thin layer, where its $x$ and $z$ components lineshapes resolved in both free-electron and peak-like contributions, resulted very similar to theoretical predictions [M. Xie et al., Phys. Rev. B $\bf{89}$, 245417 (2014)] of a two dimensional electron gas confined in the normal direction of a perovskite interface. In particular, the small but finite conducting character of the $z$ component provides a unambiguous signature of the presence of the 2DEG in the Cu/Bi$_2$O$_3$ system. Although the original constituent materials do not possess spin-orbit coupling (SOC), the resulting interfacial hybridization of such states induce electronic asymmetric wave functions. This work demonstrates the detection of 2DEG in amorphous crystals allowing to study its challenging interfacial phenomena such as SOC and interface-bulk coupling, overcoming an experimental impediment which has hold back for decades important advancements for the understanding of 2DEGs in amorphous materials.


title: Exotic Bohmian arrival times of spin-1/2 particles I-An analytical treatment

abstract: It is well known that orthodox quantum mechanics does not make unambiguous predictions for the statistics in arrival time (or time-of-flight) experiments. Bohmian mechanics (or de Broglie-Bohm theory) offers a distinct conceptual advantage in this regard, owing to the well defined concepts of point particles and trajectories embedded in this theory. We revisit a recently proposed experiment [S. Das and D. D\"urr, Sci. Rep. (2019)], the numerical analysis of which revealed a striking spin dependence in the (Bohmian) time-of-arrival distributions of a spin-1/2 particle. We present here a mathematically tractable variant of the same experiment, where the predicted effects can be established rigorously. We also obtain some new results that can be compared with experiment.


title: A hidden constraint on the Hamiltonian formulation of relativistic worldlines

abstract: Gauge theories with general covariance are particularly reluctant to quantization. We discuss the example of the Hamiltonian formulation of the relativistic point particle that, despite its apparent simplicity, is of crucial importance since a number of point particle systems can be cast into this form on a higher dimensional Rindler background, as recently pointed out by Hojman. It is shown that this system can be equipped with a hidden local, symmetry generating, constraint which on the one hand does not bother the classical evolution and on the other hand simplifies the realization of the path integral quantization. Even though the positive impact of the hidden symmetry is more evident in the Lagrangian version of the theory, it is still present through the suggested Hamiltonian constraint.


title: Longitudinal Double-Spin Asymmetries for Dijet Production at Intermediate Pseudorapidity in Polarized $pp$ Collisions at $\sqrt{s}$ = 200 GeV

abstract: One of the primary goals of the RHIC spin program is to determine the spin-dependent gluon distribution, $\Delta g(x)$, of the proton. The measurements of the 2009 longitudinal double-helicity asymmetry, $A_{LL}$, for mid-rapidity inclusive jet and $\pi^{0}$ production place strong constraints on $\Delta g(x)$ and, for the first time, find evidence for non-zero gluon polarization values for partonic momentum fraction $x$ greater than 0.05. In contrast to inclusive jets, dijet correlation measurements provide access to partonic kinematics at leading order, and thus give better constraints on the behavior of $\Delta g(x)$ as a function of gluon momentum fraction. Furthermore, dijet measurements at higher rapidity probe the lower $x$ values where $\Delta G$ is poorly constrained.\\ In these proceedings, we present the first measurement of $A_{LL}$ for dijets with at least one jet reconstructed within the pseudorapidity range 0.8 $< \eta <$ 1.8 at STAR. The dijets were measured in polarized proton+proton collisions at a center-of-mass energy $\sqrt{s}$ = 200 GeV. Values of $A_{LL}$ are determined for several distinct event topologies, defined by the jet pseudorapidities, and span a range of parton momentum fraction $x$ down to $x \sim 0.01$. The measured asymmetries are found to be consistent with the predictions of global analyses that incorporate the results of previous RHIC measurements. They will provide new constraints on $\Delta g(x)$ in this poorly constrained region when included in future global analyses.


title: Properties of Interaction Networks in Compressed Two and Three Dimensional Particulate Systems

abstract: We consider two (2D) and three (3D) dimensional granular systems exposed to compression, and ask what is the influence of the number of physical dimensions on the properties of the interaction networks that spontaneously form as these systems evolve. The study is carried out based on discrete element simulations of frictional disks in 2D and spheres in 3D. The main finding is that both the number of physical dimensions and the type of particle-particle interaction influence significantly the properties of interaction networks. These networks play an important role in bridging the microscale (particle size) and macroscale (system size), thus both aspects (the interaction model and dimensionality) are carefully considered. Our work uses a combination of tools and techniques, including percolation study, statistical analysis, as well as algebraic topology-based techniques. In many instances we find that different techniques and measures provide complementary information that, when combined, allow for gaining better insight into the properties of interaction networks.


title: ALMA Autocorrelation Spectroscopy of Comets: The HCN/H^13CN ratio in C/2012 S1 (ISON)

abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful tool for high-resolution mapping of comets, but the main interferometer (comprised of 50x12-m antennas) is insensitive to the largest coma scales due to a lack of very short baselines. In this work, we present a new technique employing ALMA autocorrelation data (obtained simultaneously with the interferometric observations), effectively treating the entire 12-m array as a collection of single-dish telescopes. Using combined autocorrelation spectra from 28 active antennas, we recovered extended HCN coma emission from comet C/2012 S1 (ISON), resulting in a fourteen-fold increase in detected line brightness compared with the interferometer. This resulted in the first detection of rotational emission from H^13CN in this comet. Using a detailed coma radiative transfer model accounting for optical depth and non-LTE excitation effects, we obtained an H^12CN/H^13CN ratio of 88+-18, which matches the terrestrial value of 89, consistent with a lack of isotopic fractionation in HCN during comet formation in the protosolar accretion disk. The possibility of future discoveries in extended sources using autocorrelation spectroscopy from the main ALMA array is thus demonstrated.


title: Investigation of the Mode-Switching Phenomenon in Pulsar B0329+54 Through Polarimetric Analysis

abstract: The phenomenon of profile mode switching in pulsars, where the stable average pulse profile changes to another stable state on the timescale of a pulsar's period, remains poorly understood. We sought to understand how pulsars undergo profile mode switching through a comparative analysis of the polarization and geometry of the two different profile modes of PSR B0329+54. The polarization behavior and fitted parameters of the rotation-vector model remain constant between modes, and the emission height remains constant as well. These similarities lend support to a model of pair production in the surface plasma that would cause a change in the available electrons and therefore the differential emission intensity.


title: Urban sensing as a random search process

abstract: We study a new random search process: the \textit{taxi-drive}. The motivation for this process comes from urban sensing, in which sensors are mounted on moving vehicles such as taxis, allowing urban environments to be opportunistically monitored. Inspired by the movements of real taxis, the taxi-drive is composed of both random and regular parts; passengers are brought to randomly chosen locations via deterministic (i.e. shortest paths) routes. We show through a numerical study that this hybrid motion endows the taxi-drive with advantageous spreading properties. In particular, on certain graph topologies it offers reduced cover times compared to persistent random walks.


title: Fourier-series expansion of the dark-energy equation of state

abstract: The dark energy component of the universe still remains as a mystery, however, several papers based on observational data have shown that its equation of state may have an oscillatory behaviour. In this paper, we provide a general description for the dark-energy equation-of-state $w(z)$ in the form of Fourier series. This description generalises some previous dynamical dark energy models and is in agreement with the $w(z)$ reconstructions. We make use of a modified version of a simple and fast Markov Chain Monte Carlo code to constraint the model parameters. For the analysis we use data from supernovae type-Ia , baryon acoustic oscillations, $H(z)$ measurements and cosmic microwave background. We provide a comparison of the proposed model with $\Lambda$CDM, $w$CDM and the standard Taylor approximation. The Fourier series expansion of $w(z)$ is preferred from $\Lambda$CDM at more than $3\sigma$ significance level based on the improvement in the fit alone. We use the Akaike criteria to perform the model comparison and found that, even though there are extra parameters, there is a slight preference of the Fourier series compared with the $\Lambda$CDM model. The preferred shape of $w(z)$ found here puts in jeopardy the single scalar field models, as they as they cannot reproduce the crossing the phantom divide line $w=-1$.