title: An Optimal Policy for Patient Laboratory Tests in Intensive Care Units

abstract: Laboratory testing is an integral tool in the management of patient care in hospitals, particularly in intensive care units (ICUs). There exists an inherent trade-off in the selection and timing of lab tests between considerations of the expected utility in clinical decision-making of a given test at a specific time, and the associated cost or risk it poses to the patient. In this work, we introduce a framework that learns policies for ordering lab tests which optimizes for this trade-off. Our approach uses batch off-policy reinforcement learning with a composite reward function based on clinical imperatives, applied to data that include examples of clinicians ordering labs for patients. To this end, we develop and extend principles of Pareto optimality to improve the selection of actions based on multiple reward function components while respecting typical procedural considerations and prioritization of clinical goals in the ICU. Our experiments show that we can estimate a policy that reduces the frequency of lab tests and optimizes timing to minimize information redundancy. We also find that the estimated policies typically suggest ordering lab tests well ahead of critical onsets--such as mechanical ventilation or dialysis--that depend on the lab results. We evaluate our approach by quantifying how these policies may initiate earlier onset of treatment.


title: Dzyaloshinskii-Moriya interaction at disordered interfaces from ab initio theory: robustness against intermixing and tunability through dusting

abstract: The Dzyaloshinskii-Moriya interaction (DMI), which is essential for the stabilization of topologically non-trivial chiral magnetic textures such as skyrmions, is particularly strong in heterostructures of ultra-thin magnetic materials and heavy elements. We explore by density-functional theory calculations the possibility to modify the magnetic properties at Co/Pt interfaces with chemical disorder. In these systems, we find a particular robustness of the DMI against intermixing. Upon dusting the interface with a third element (all $4d$ transition metals and B, Cu, Au and Bi), a strong reduction of the DMI is predicted. This opens up possibilities to tune the DMI through the degrees of intermixing and dusting.


title: A variational approach to nucleation simulation

abstract: We study by computer simulation the nucleation of a supersaturated Lennard-Jones vapor into the liquid phase. The large free energy barriers to transition make the time scale of this process impossible to study by ordinary molecular dynamics simulations.Therefore we use a recently developed enhanced sampling method [Valsson and Parrinello, Phys. Rev. Lett. 113, 090601 (2014)] based on the variational determination of a bias potential. We differ from previous applications of this method in that the bias is constructed on the basis of the physical model provided by the classical theory of nucleation. We examine the technical problems associated with this approach. Our results are very satisfactory and will pave the way for calculating the nucleation rates in many systems.


title: Wave propagation in infinite nonlinear acoustic metamaterial beam by considering the third harmonic generation

abstract: Nonlinear acoustic metamaterial (NAM) initiates new fields for controlling elastic waves. In this work, the flexural wave propagation in the half-infinite NAM beam consisting of periodic Duffing resonators is reported by considering the third harmonic generation (THG). Different analytical methods are proposed to describe the wave propagation in the equivalent homogenous medium. Then their effectiveness and accuracy are demonstrated in comparison with the finite element methods. We unveil analytically and numerically extensive physical properties of the strongly nonlinear AM, including the nonlinear resonance in a cell, the effective density, nonlinear locally resonant (NLR) bandgap, propagations and couplings of the fundamental and the third harmonics. These characteristics are highly interrelated, which facilitates the prediction of functionalities. In the near field, the identical bifurcation frequency of these features acts as the start frequency of the NLR bandgap for fundamental waves, whose width is narrower for a stronger nonlinearity. While in the far field, the NLR bandgap characterizes a distance-amplitude-dependent behavior leading to a self-adaptive bandwidth. Moreover, the transmission in the passband of the infinite NAM is different from the chaotic band effect of finite NAMs, and it is influenced by the shifted NLR gap. Our work will promote future studies and constructions of NAMs with novel properties.


title: Maintaining both availability and integrity of communications: Challenges and guidelines for data security and privacy during disasters and crises

abstract: Communications play a vital role in the response to disasters and crises. However, existing communications infrastructure is often impaired, destroyed or overwhelmed during such events. This leads to the use of substitute communications solutions including analog two-way radio or unsecured internet access. Often provided by unknown third parties, these solutions may have less sophisticated security characteristics than is desirable. While substitute communications are often invaluable, care is required to minimize the risk to NGOs and individuals stemming from the use of communications channels with reduced or unknown security properties. This is particularly true if private information is involved, including the location and disposition of individuals and first responders. In this work we enumerate the principal risks and challenges that may arise, and provide practical guidelines for mitigating them during crises. We take plausible threats from contemporary disaster and crisis events into account and discuss the security and privacy features of state-of-the-art communications mechanisms.


title: Conducting a Large-scale Field Test of a Smartphone-based Communication Network for Emergency Response

abstract: Smartphone-based communication networks form a basis for services in emergency response scenarios, where communication infrastructure is impaired or overloaded. Still, their design and evaluation are largely based on simulations that rely on generic mobility models and weak assumptions regarding user behavior. For a realistic assessment, scenario-specific models are essential. To this end, we conducted a large-scale field test of a set of emergency services that relied solely on ad hoc communication. Over the course of one day, we gathered data from smartphones distributed to 125 participants in a scripted disaster event. In this paper, we present the scenario, measurement methodology, and a first analysis of the data. Our work provides the first trace combining user interaction, mobility, and additional sensor readings of a large-scale emergency response scenario, facilitating future research.


title: Learning ReLU Networks on Linearly Separable Data: Algorithm, Optimality, and Generalization

abstract: Neural networks with REctified Linear Unit (ReLU) activation functions (a.k.a. ReLU networks) have achieved great empirical success in various domains. Nonetheless, existing results for learning ReLU networks either pose assumptions on the underlying data distribution being e.g. Gaussian, or require the network size and/or training size to be sufficiently large. In this context, the problem of learning a two-layer ReLU network is approached in a binary classification setting, where the data are linearly separable and a hinge loss criterion is adopted. Leveraging the power of random noise perturbation, this paper presents a novel stochastic gradient descent (SGD) algorithm, which can \emph{provably} train any single-hidden-layer ReLU network to attain global optimality, despite the presence of infinitely many bad local minima, maxima, and saddle points in general. This result is the first of its kind, requiring no assumptions on the data distribution, training/network size, or initialization. Convergence of the resultant iterative algorithm to a global minimum is analyzed by establishing both an upper bound and a lower bound on the number of non-zero updates to be performed. Moreover, generalization guarantees are developed for ReLU networks trained with the novel SGD leveraging classic compression bounds. These guarantees highlight a key difference (at least in the worst case) between reliably learning a ReLU network as well as a leaky ReLU network in terms of sample complexity. Numerical tests using both synthetic data and real images validate the effectiveness of the algorithm and the practical merits of the theory.


title: Photoacoustic Sensing of Trapped Fluids in Nanoporous Thin Films: Device Engineering and Sensing Scheme

abstract: Accessing fluid infiltration in nanogranular coatings is an outstanding challenge, of relevance for applications ranging from nanomedicine to catalysis. A sensing platform, allowing to quantify the amount of fluid infiltrated in a nanogranular ultrathin coating, with thickness in the 10 to 40 nm range, is here proposed and theoretically investigated by multiscale modelling. The scheme relies on impulsive photoacoustic excitation of hypersonic mechanical breathing modes in engineered gas-phase synthesised nanogranular metallic ultathin films and time-resolved acousto-optical read-out of the breathing modes frequency shift upon liquid infiltration. A superior sensitivity, exceeding 26x103 cm^2/g, is predicted upon equivalent areal mass loading of a few ng/mm^2. The capability of the present scheme to discriminate among different infiltration patterns is discussed. The platform is an ideal tool to investigate nano fluidics in granular materials and naturally serves as a distributed nanogetter coating, integrating fluid sensing capabilities. The proposed scheme is readily extendable to other nanoscale and mesoscale porous materials.


title: Measuring carbon nanotube vibrations using a single-electron transistor as a fast linear amplifier

abstract: We demonstrate sensitive and fast electrical measurements of a carbon nanotube mechanical resonator. The nanotube is configured as a single-electron transistor, whose conductance is a sensitive transducer for its own displacement. Using an impedance-matching circuit followed by a cryogenic amplifier, the vibrations can be monitored in real time. The sensitivity of this continuous displacement measurement approaches within a factor 470 of the standard quantum limit.


title: On Mechanical Behavior of Elastomeric Networks: Effects of Random Porous Microstructure

abstract: An assumption in micromechanical analysis of polymers is that the constitutive polymeric media is non-porous. Non-porosity of media, however, is merely a simplifying assumption. In this paper, we neglect this assumption and studied polymer networks with a different porosity volume fraction. A random morphology description function is used to model the porosity of the network and nonlinear finite element analyses are conducted to perform structural analysis of porous polymer networks. The results show that the porosity effect is significant in mechanical behavior of polymer networks and may increase the maximum Von-Mises stress drastically.


title: Evolutionary algorithms converge towards evolved biological photonic structures

abstract: Nature features a plethora of extraordinary photonic architectures that have been optimized through natural evolution. While numerical optimization is increasingly and successfully used in photonics, it has yet to replicate any of these complex naturally occurring structures. Using evolutionary algorithms directly inspired by natural evolution, we have retrieved emblematic natural photonic structures, indicating how such regular structures might have spontaneously emerged in nature and to which precise optical or fabrication constraints they respond. Comparisons between algorithms show that recombination between individuals inspired by sexual reproduction confers a clear advantage in this context of modular problems and suggest further ways to improve the algorithms. Such an in silico evolution can also suggest original and elegant solutions to practical problems, as illustrated by the design of counter-intuitive anti-reflective coating for solar cells.


title: Algorithms and image formation in orbital tomography

abstract: Orbital tomography has recently been established as a technique to reconstruct molecular orbitals directly from photoemission data using iterative phase retrieval algorithms. In this work, we present a detailed description of steps for processing of the photoemission data followed by an improved iterative phase retrieval procedure and the interpretation of reconstructed two-dimensional orbital distributions. We address the issue of background subtraction by suggesting a signal restoration routine based on the maximization of mutual information algorithm and solve the problem of finding the geometrical center in the reconstruction by using a tight-centered object support in a two-step phase retrieval procedure. The proposed image processing and improved phase retrieval procedures are used to reconstruct the highest occupied molecular orbital of pentacene on Ag(110), using photoemission data only. The results of the reconstruction agree well with the density functional theory simulation, modified to comply with the experimental conditions. By comparison with photoelectron holography, we show that the reconstructed two-dimensional orbital distribution can be interpreted as a superposition of the in-focus orbital distribution evaluated at the z = 0 plane and out-of-focus distributions evaluated at other z = const planes. Three-dimensional molecular orbital distributions could thus be reconstructed directly from two-dimensional photoemission data, provided the axial resolution of the imaging system is high enough.


title: The Frequency Response Function of the Creep Compliance

abstract: Motivated from the need to convert time-dependent rheometry data into complex frequency response functions, this paper studies the frequency response function of the creep compliance that is coined the complex creep function. While for any physically realizable viscoelastic model the Fourier transform of the creep compliance diverges in the classical sense, the paper shows that the complex creep function, in spite of exhibiting strong singularities, it can be constructed with the calculus of generalized functions. The mathematical expressions of the real and imaginary parts of the Fourier transform of the creep compliance of simple rheological networks derived in this paper are shown to be Hilbert pairs; therefore, returning back in the time domain a causal creep compliance. The paper proceeds by showing how a measured creep compliance of a solid-like or a fluid-like viscoelastic material can be decomposed into elementary functions with parameters that can be identified from best fit of experimental data. The proposed technique allows for a direct determination of the parameters of the corresponding viscoelastic models and leads to dependable expressions of their complex-frequency response functions.


title: FPGA based demodulation of laser induced fluorescence in plasmas

abstract: We present a field programmable gate array (FPGA)-based system that counts photons from laser-induced fluorescence (LIF) on a laboratory plasma. This is accomplished with FPGA-based up/down counters that demodulate the data, giving a background-subtracted LIF signal stream that is updated with a new point as each laser amplitude modulation cycle completes. We demonstrate using the FPGA to modulate a laser at 1 MHz and demodulate the resulting LIF data stream. This data stream is used to calculate an LIF-based measurement sampled at 1 MHz of a plasma ion fluctuation spectrum.


title: Differential Colorimetry Measurements of Fluctuation Growth in Nanofilms Exposed to Large Surface Thermal Gradients

abstract: Slender liquid nanofilms exposed to large surface thermal gradients are known to undergo thickness fluctuations which rapidly self-organize into arrays of nanoprotrusions with a separation distance of tens of microns. We previously reported good agreement between measurements of the characteristic spacing and the wavelength of the most unstable mode predicted by a linear stability analysis based on a long wavelength thermocapillary model. Here we focus on differential colorimetry measurements to quantify early time out-of-plane growth of protrusions for peak heights spanning 20 to 200 nm. Analysis of peak heights based on shape reconstruction reveals robust exponential growth. Good quantitative agreement of the growth rates with the thermocapillary model is obtained using a single fit constant to account for material parameters of nanofilms that could not be measured directly. These findings lend further support to the conjecture that the array protrusions uncovered almost two decades ago stem from a linear instability whose growth rate is controlled by thermocapillary forces counterbalanced by capillary forces.


title: A Hassle-Free Machine Learning Method for Cohort Selection of Clinical Trials

abstract: Traditional text classification techniques in clinical domain have heavily relied on the manually extracted textual cues. This paper proposes a generally supervised machine learning method that is equally hassle-free and does not use clinical knowledge. The employed methods were simple to implement, fast to run and yet effective. This paper proposes a novel named entity recognition (NER) based an ensemble system capable of learning the keyword features in the document. Instead of merely considering the whole sentence/paragraph for analysis, the NER based keyword features can stress the important clinic relevant phases more. In addition, to capture the semantic information in the documents, the FastText features originating from the document level FastText classification results are exploited.


title: Searches for new heavy quarks in ATLAS

abstract: A search for new heavy quarks focusing on recent vector-like quark searches with the ATLAS detector at the CERN Large Hadron Collider is presented. Two recent searches targeting the pair production of type vector-like quarks are described. The first search is sensitive to vector-like up-type quark ($T$) decays to a $t$ quark and either a Standard Model Higgs boson or a $Z$ boson. The second search is primarily sensitive to $T$ decays to $W$ boson and a $b$ quark. Additionally, the results can be interpreted for alternative VLQ decays.


title: SAT.STFR.FRQ (UWA) Detail Design Report (LOW)

abstract: The Square Kilometre Array (SKA) project is an international effort to build the world's most sensitive radio telescope operating in the 50 MHz to 14 GHz frequency range. Construction of the SKA is divided into phases, with the first phase (SKA1) accounting for the first 10% of the telescope's receiving capacity. During SKA1, a Low-Frequency Aperture Array (LFAA) comprising over a hundred thousand individual dipole antenna elements will be constructed in Western Australia (SKA1-LOW), while an array of 197 parabolic-receptor antennas, incorporating the 64 receptors of MeerKAT, will be constructed in South Africa (SKA1-MID). Radio telescope arrays, such as the SKA, require phase-coherent reference signals to be transmitted to each antenna site in the array. In the case of the SKA, these reference signals are generated at a central site and transmitted to the antenna sites via fibre-optic cables up to 175 km in length. Environmental perturbations affect the optical path length of the fibre and act to degrade the phase stability of the reference signals received at the antennas, which has the ultimate effect of reducing the fidelity and dynamic range of the data . Given the combination of long fibre distances and relatively high frequencies of the transmitted reference signals, the SKA needs to employ actively-stabilised frequency transfer technologies to suppress the fibre-optic link noise in order to maintain phase-coherence across the array.


title: Plus-one generated and next to free arrangements of hyperplanes

abstract: We introduce a new class of arrangements of hyperplanes, called (strictly) plus-one generated arrangements, from algebraic point of view. Plus-one generatedness is close to freeness, i.e., plus-one generated arrangements have their logarithmic derivation modules generated by dimension plus one elements, with relations containing one linear form coefficient. We show that strictly plus-one generated arrangements can be obtained if we delete a hyperplane from free arrangements. We show a relative freeness criterion in terms of plus-one generatedness. In particular, for plane arrangements, we show that a free arrangement is in fact surrounded by free or strictly plus-one generated arrangements. We also give several applications.


title: Probabilistic forecasting of heterogeneous consumer transaction-sales time series

abstract: We present new Bayesian methodology for consumer sales forecasting. With a focus on multi-step ahead forecasting of daily sales of many supermarket items, we adapt dynamic count mixture models to forecast individual customer transactions, and introduce novel dynamic binary cascade models for predicting counts of items per transaction. These transactions-sales models can incorporate time-varying trend, seasonal, price, promotion, random effects and other outlet-specific predictors for individual items. Sequential Bayesian analysis involves fast, parallel filtering on sets of decoupled items and is adaptable across items that may exhibit widely varying characteristics. A multi-scale approach enables information sharing across items with related patterns over time to improve prediction while maintaining scalability to many items. A motivating case study in many-item, multi-period, multi-step ahead supermarket sales forecasting provides examples that demonstrate improved forecast accuracy in multiple metrics, and illustrates the benefits of full probabilistic models for forecast accuracy evaluation and comparison. Keywords: Bayesian forecasting; decouple/recouple; dynamic binary cascade; forecast calibration; intermittent demand; multi-scale forecasting; predicting rare events; sales per transaction; supermarket sales forecasting


title: Improving Generalization via Scalable Neighborhood Component Analysis

abstract: Current major approaches to visual recognition follow an end-to-end formulation that classifies an input image into one of the pre-determined set of semantic categories. Parametric softmax classifiers are a common choice for such a closed world with fixed categories, especially when big labeled data is available during training. However, this becomes problematic for open-set scenarios where new categories are encountered with very few examples for learning a generalizable parametric classifier. We adopt a non-parametric approach for visual recognition by optimizing feature embeddings instead of parametric classifiers. We use a deep neural network to learn the visual feature that preserves the neighborhood structure in the semantic space, based on the Neighborhood Component Analysis (NCA) criterion. Limited by its computational bottlenecks, we devise a mechanism to use augmented memory to scale NCA for large datasets and very deep networks. Our experiments deliver not only remarkable performance on ImageNet classification for such a simple non-parametric method, but most importantly a more generalizable feature representation for sub-category discovery and few-shot recognition.


title: On a different expression of the classical limit of quantum mechanics

abstract: The principle of correspondence (or classical limit) is essential in quantum mechanics. Yet, how and why quantum phenomena vanish at the macroscopic scale are issues still open to debate. Here, quantum mechanical predictions for Greenberger-Horne-Zeilinger states of qubits are shown to be easier to reproduce with a classical model as the number of particles increases, even in the absence of loopholes or conspiratorial mechanisms of any kind. It is conjectured that this result may lead to the simplest way to express the principle of correspondence.


title: Measuring Market Performance with Stochastic Demand: Price of Anarchy and Price of Uncertainty

abstract: Globally operating suppliers face the rising challenge of wholesale pricing under scarce data about retail demand, in contrast to better informed, locally operating retailers. At the same time, as local businesses proliferate, markets congest and retail competition increases. To capture these strategic considerations, we employ the classic Cournot model and extend it to a two-stage supply chain with an upstream supplier who operates under demand uncertainty and multiple downstream retailers who compete over quantity. The supplier's belief about retail demand is modeled via a continuous probability distribution function F. If F has the decreasing generalized mean residual life property, then the supplier's optimal pricing policy exists and is the unique fixed point of the mean residual life function. We evaluate the realized Price of Uncertainty and show that there exist demand levels for which market performs better when the supplier prices under demand uncertainty. In general, performance worsens for lower values of realized demand. We examine the effects of increasing competition on supply chain efficiency via the realized Price of Anarchy and complement our findings with numerical results.


title: Pixel Objectness: Learning to Segment Generic Objects Automatically in Images and Videos

abstract: We propose an end-to-end learning framework for segmenting generic objects in both images and videos. Given a novel image or video, our approach produces a pixel-level mask for all "object-like" regions---even for object categories never seen during training. We formulate the task as a structured prediction problem of assigning an object/background label to each pixel, implemented using a deep fully convolutional network. When applied to a video, our model further incorporates a motion stream, and the network learns to combine both appearance and motion and attempts to extract all prominent objects whether they are moving or not. Beyond the core model, a second contribution of our approach is how it leverages varying strengths of training annotations. Pixel-level annotations are quite difficult to obtain, yet crucial for training a deep network approach for segmentation. Thus we propose ways to exploit weakly labeled data for learning dense foreground segmentation. For images, we show the value in mixing object category examples with image-level labels together with relatively few images with boundary-level annotations. For video, we show how to bootstrap weakly annotated videos together with the network trained for image segmentation. Through experiments on multiple challenging image and video segmentation benchmarks, our method offers consistently strong results and improves the state-of-the-art for fully automatic segmentation of generic (unseen) objects. In addition, we demonstrate how our approach benefits image retrieval and image retargeting, both of which flourish when given our high-quality foreground maps. Code, models, and videos are at:http://vision.cs.utexas.edu/projects/pixelobjectness/


title: Profinite groups with an automorphism whose fixed points are right Engel

abstract: An element $g$ of a group $G$ is said to be right Engel if for every $x\in G$ there is a number $n=n(g,x)$ such that $[g,{}_{n}x]=1$. We prove that if a profinite group $G$ admits a coprime automorphism $\varphi$ of prime order such that every fixed point of $\varphi$ is a right Engel element, then $G$ is locally nilpotent.


title: On numerical study of attractors of ODEs

abstract: Due to errors of numerical methods, qualitative behavior of the numerically simulated solution may be completely different from qualitative behavior of the real trajectory. Therefor in study of attractors the more useful idea is to simulate a bundle of trajectories for a (relatively) short time (see Dellnitz etc.) We suggest a method of modelling Milnor's attractor of dissipative dynamical system. It is based on numerical simulation of bundle of trajectories, emanating from some set, and locating regions with high density of trajectories. We present 2D and 3D examples of application of our method.


title: Heavy Neutral Leptons and displaced vertices at LHC

abstract: Heavy neutral leptons are present in many well-motivated beyond the Standard Model theories, sometimes being accessible at present colliders. Depending on their masses and couplings they could be long-lived and lead to events with displaced vertices, and thus to promising signatures due to low Standard Model background. We revisit the potential of the LHC to discover this kind of new particles via searches for displaced vertices, which can probe masses of few GeV for mixings currently allowed by experimental constraints. We also discuss the importance of considering all the possible production channels, including the production in association with light neutrinos, in order to fully explore this region of the parameter space.


title: Neural Machine Translation Inspired Binary Code Similarity Comparison beyond Function Pairs

abstract: Binary code analysis allows analyzing binary code without having access to the corresponding source code. A binary, after disassembly, is expressed in an assembly language. This inspires us to approach binary analysis by leveraging ideas and techniques from Natural Language Processing (NLP), a rich area focused on processing text of various natural languages. We notice that binary code analysis and NLP share a lot of analogical topics, such as semantics extraction, summarization, and classification. This work utilizes these ideas to address two important code similarity comparison problems. (I) Given a pair of basic blocks for different instruction set architectures (ISAs), determining whether their semantics is similar or not; and (II) given a piece of code of interest, determining if it is contained in another piece of assembly code for a different ISA. The solutions to these two problems have many applications, such as cross-architecture vulnerability discovery and code plagiarism detection. We implement a prototype system INNEREYE and perform a comprehensive evaluation. A comparison between our approach and existing approaches to Problem I shows that our system outperforms them in terms of accuracy, efficiency and scalability. And the case studies utilizing the system demonstrate that our solution to Problem II is effective. Moreover, this research showcases how to apply ideas and techniques from NLP to large-scale binary code analysis.


title: Analysis and Classification of Critical Success Factors for Business Conversation Manager (BCM) System Implementation in Service-Oriented Enterprises

abstract: Business conversation manager is a system which makes organization's communications possible on computer system platforms and also within defined processes and relationships. This research reviews the literature of BCM systems and introduces and classifies the most important factors for successful implementation by investigating various references. The aim is to identify critical success factors in order to avoid high risk probabilities and facilitate the implementation and finally getting better outcomes. Next, factors are extracted from experts' points of view and they depicted as a model to show critical factors.


title: Probing heavy ion collisions using quark and gluon jet substructure with machine learning

abstract: Understanding the inner working of the quark-gluon plasma requires complete and precise jet substructure studies in heavy ion collisions. In this proceeding we discuss the use of quark and gluon jets as independent probes, and how their classification allows us to uncover regions of QCD phase space sensitive to medium dynamics. We introduce the telescoping deconstruction (TD) framework to capture complete jet information and show that TD observables reveal fundamental properties of quark and gluon jets and their modifications in the medium. We draw connections to soft-drop subjet distributions and illuminate medium-induced jet modifications using Lund diagrams. The classification is also studied using a physics-motivated, multivariate analysis of jet substructure observables. Moreover, we apply image-recognition techniques by training a deep convolutional neural network on jet images to benchmark classification performances. We find that the quark gluon discrimination performance worsens in \textsc{Jewel}-simulated heavy ion collisions due to significant soft radiation affecting soft jet substructures. This work suggests a systematic framework for jet studies and facilitates direct comparisons between theoretical calculations and measurements in heavy ion collisions.


title: On the exact quantum scale invariance of three-dimensional reduced QED theories

abstract: An effective quantum field theory description of graphene in the ultra-relativistic regime is given by reduced QED aka. pseudo QED aka. mixed-dimensional QED. It has been speculated in the literature that reduced QED constitutes an example of a specific class of hard-to-find theories: an interacting CFT in more than two dimensions. This speculation was based on two-loop perturbation theory. Here, we give a proof of this feature, namely the exact vanishing of the b-function, thereby showing that reduced QED can effectively be considered as an interacting (boundary) CFT, underpinning recent work in this area. The argument, valid for both two- and four-component spinors, also naturally extends to an exactly marginal deformation of reduced QED, thence resulting in a non-supersymmetric conformal manifold. The latter corresponds to boundary layer fermions between two different dielectric half-spaces.


title: Regime-Switching Temperature Dynamics Model for Weather Derivatives

abstract: Weather is a key production factor in agricultural crop production and at the same time the most significant and least controllable source of peril in agriculture. These effects of weather on agricultural crop production have triggered a widespread support for weather derivatives as a means of mitigating the risk associated with climate change on agriculture. However, these products are faced with basis risk as a result of poor design and modelling of the underlying weather variable (temperature). In order to circumvent these problems, a novel time-varying mean-reversion L\'evy regime-switching model is used to model the dynamics of the deseasonalized temperature dynamics. Using plots and test statistics, it is observed that the residuals of the deseasonalized temperature data are not normally distributed. To model the non-normality in the residuals, we propose using the hyperbolic distribution to capture the semi-heavy tails and skewness in the empirical distributions of the residuals for the shifted regime. The proposed regime-switching model has a mean-reverting heteroskedastic process in the base regime and a L\'evy process in the shifted regime. By using the Expectation-Maximization algorithm, the parameters of the proposed model are estimated. The proposed model is flexible as it modelled the deseasonalized temperature data accurately.


title: Testing Cosmic Censorship Conjecture for Extremal and Near-extremal (2+1)-dimensional MTZ Black Holes

abstract: We test the validity of the weak cosmic censorship conjecture for the ($2+1$)-dimensional charged anti-de Sitter black hole solution, which was derived by Martinez, Teitelboim, and Zanelli (MTZ). We first construct a thought experiment by throwing test charged particles on an extremal MTZ black hole. We derive that extremal ($2+1$) dimensional black holes can be overcharged by test particles, unlike their analogues in $4$ and higher dimensions. Nearly-extremal black holes can also be overcharged, by a judicious choice of energy and charge for the test particles in the case when we ignore the second order effects. Contrary to this, nearly-extremal black holes cannot be overcharged by the second order perturbation, obeying the weak cosmic censorship conjecture.


title: Computing Equilibria in Atomic Splittable Polymatroid Congestion Games with Convex Costs

abstract: In this paper, we compute $\epsilon$-approximate Nash equilibria in atomic splittable polymatroid congestion games with convex Lipschitz continuous cost functions. The main approach relies on computing a pure Nash equilibrium for an associated integrally-splittable congestion game, where players can only split their demand in integral multiples of a common packet size. It is known that one can compute pure Nash equilibria for integrally-splittable congestion games within a running time that is pseudo-polynomial in the aggregated demand of the players. As the main contribution of this paper, we decide for every $\epsilon>0$, a packet size $k_{\epsilon}$ and prove that the associated $k_{\epsilon}$-splittable Nash equilibrium is an $\epsilon$-approximate Nash equilibrium for the original game. We further show that our result applies to multimarket oligopolies with decreasing, concave Lipschitz continuous price functions and quadratic production costs: there is a polynomial time transformation to atomic splittable polymatroid congestion games implying that we can compute $\epsilon$-approximate Cournot-Nash equilibria within pseudo-polynomial time.


title: An iterative domain decomposition, spectral finite element method on non-conforming meshes suitable for high frequency Helmholtz problems

abstract: The purpose of this research is to describe an efficient iterative method suitable for obtaining high accuracy solutions to high frequency time-harmonic scattering problems. The method allows for both refinement of local polynomial degree and non-conforming mesh refinement, including multiple hanging nodes per edge. Rather than globally assemble the finite element system, we describe an iterative domain decomposition method which can use element-wise fast solvers for elements of large degree. Since continuity between elements is enforced through moment equations, the resulting constraint equations are hierarchical. We show that, for high frequency problems, a subset of these constraints should be directly enforced, providing the coarse space in the dual-primal domain decomposition method. The subset of constraints is chosen based on a dispersion criterion involving mesh size and wavenumber. By increasing the size of the coarse space according to this criterion, the number of iterations in the domain decomposition method depends only weakly on the wavenumber. We demonstrate this convergence behaviour on standard domain decomposition test problems and conclude the paper with application of the method to electromagnetic problems in two dimensions. These examples include beam steering by lenses and photonic crystal waveguides, as well as radar cross section computation for dielectric, perfect electric conductor, and electromagnetic cloak scatterers.


title: Pseudo-Hermitian position and momentum operators, Hermitian Hamiltonian, and deformed oscillators

abstract: The recently introduced by us two- and three-parameter ($p,q$)- and ($p,q,\mu$)-deformed extensions of the Heisenberg algebra were explored under the condition of their direct link with the respective (nonstandard) deformed quantum oscillator algebras. In this paper we explore certain Hermitian Hamiltonian build in terms of non-Hermitian position and momentum operators obeying definite $\eta(N)$-pseudo-Hermiticity properties. A generalized nonlinear (with the coefficients depending on the excitation number operator $N$) one-mode Bogolyubov transformation is developed as main tool for the corresponding study. Its application enables to obtain the spectrum of "almost free" (but essentially nonlinear) Hamiltonian.


title: On the requirement for mixing-heating to utilize bubble cosmic rays to heat the intracluster medium

abstract: I conduct simple analytical estimates and conclude that mixing by vortices is a more efficient process to transfer the cosmic ray energy of jet-inflated bubbles to the intracluster medium (ICM) than streaming of cosmic rays along magnetic field lines. Jets and the bubbles they inflate transfer heat to the ambient gas in cooling flows in cluster of galaxies and in galaxies. The internal energy of the jet-inflated bubbles is dominate by very hot thermal gas and/or cosmic rays. Cosmic rays that stream along magnetic field lines that connect the bubbles with the ICM heat the ICM as their energy is dissipated there. I find that about half of the cosmic ray energy is dissipated in the bubbles themselves. I also find that the ICM volume that the cosmic ray streaming process heats is only about five times as large as the volume of the bubbles. The outcome of heating by streaming only is that the cosmic rays form a larger bubble filled with very hot thermal gas. Therefore, there is a requirement for a more efficient process to transfer the internal energy of the bubbles to the ICM. I suggest that this process is heating by mixing that operates very well for both cosmic rays and the very hot thermal gas inside the bubbles. This leaves mixing-heating to be the dominant heating process of cooling flows.


title: Evolutionary Models for 15 Galactic Supernova Remnants with New Distances

abstract: Recent studies using 21 cm HI line and $^{13}$CO line observations in the inner part of the Galaxy have resulted in new distances for 30 Galactic supernova remnants (SNRs). 15 of those remnants have observed X-ray spectra, for which shocked-gas temperatures and emission measures are measured. Here we apply spherically symmetric SNR evolution models to these 15 remnants to obtain estimates for ages, explosion energies, circum-stellar medium densities and profiles (uniform or wind-type).From the distribution of ages we obtain a supernova birth rate and estimate incompleteness. The energies and densities can be well fit with log-normal distributions. The distribution of explosion energies is very similar to that of SNRs in the Large Magellanic Cloud (LMC), suggesting SN explosions in the LMC and in the Galaxy are very similar. The density distribution has higher mean density for Galactic SNRs than for LMC SNRs, by a factor $\sim$2.5.


title: $q$-Analogues of some series for powers of $\pi$

abstract: We obtain $q$-analogues of several series for powers of $\pi$. For example, the identity $$\sum_{k=0}^\infty\frac{(-1)^k}{(2k+1)^3}=\frac{\pi^3}{32}$$ has the following $q$-analogue: \begin{equation*} \sum_{k=0}^\infty(-1)^k\frac{q^{2k}(1+q^{2k+1})}{(1-q^{2k+1})^3}=\frac{(q^2;q^4)_{\infty}^2(q^4;q^4)_{\infty}^6} {(q;q^2)_{\infty}^4}, \end{equation*} where $q$ is any complex number with $|q|<1$. We also give $q$-analogues of four new series for powers of $\pi$ found by the second author.


title: Effect of Cavitating-bubble Generators on the Dynamics of Unsteady Cloud Cavitation

abstract: The present paper is devoted to the 3D numerical study of unsteady cloud cavitation around the CAV2003 benchmark hydrofoil using a passive cavitation controller method called cavitation-bubble generator (CGs). We investigate this passive controller on the hydrofoil surface to affect the whole processes of vaporization, bubble generation and bubble implosion. Firstly we simulate the unsteady cavitating flow around the hydrofoil without CGs using a Partially- averaged Navier Stokes (PANS) method. The turbulence approach is coupled with a mass transfer model which implements to an open source finite-volume package. Then we tested different sizes and angles of CGs on the hydrofoil surface. The results of an appropriate design of the CGs shows a significant decrease of the maximum length of the attached cavity and a large decrease of wall-pressure fluctuation intensity on the hydrofoil surface.


title: Intrinsic Point Defects in Ultrathin 1T-PtSe2 Layers

abstract: Among two-dimensional (2D) transition metal dichalcogenides (TMDs), platinum diselenide (PtSe2) stands at a unique place in the sense that it undergoes a phase transition from type-II Dirac semimetal to indirect-gap semiconductor as thickness decreases. Defects in 2D TMDs are ubiquitous and play crucial roles in understanding and tuning electronic, optical, and magnetic properties. Here we investigate intrinsic point defects in ultrathin 1T-PtSe2 layers grown on mica through the chemical vapor transport (CVT) method, using scanning tunneling microscopy and spectroscopy (STM/STS) and first-principles calculations. We observed five types of distinct defects from STM topography images and obtained the local density of states of the defects. By combining the STM results with the first-principles calculations, we identified the types and characteristics of these defects, which are Pt vacancies at the topmost and next monolayers, Se vacancies in the topmost monolayer, and Se antisites at Pt sites within the topmost monolayer. Our study shows that the Se antisite defects are the most abundant with the lowest formation energy in a Se-rich growth condition, in contrast to cases of 2D molybdenum disulfide (MoS2) family. Our findings would provide critical insight into tuning of carrier mobility, charge carrier relaxation, and electron-hole recombination rates by defect engineering or varying growth condition in few-layer 1T-PtSe2 and other related 2D materials.


title: Is the enhancement of type II radio bursts during CME interactions related to the associated solar energetic particle event?

abstract: We investigated 64 pairs of interacting-CME events identified by the simultaneous observations of SOHO and STEREO spacecraft from 2010 January to 2014 August, to examine the relationship between the large SEP events in the energy of 25-60MeV and the properties of the interacting CMEs. We found that during CME interactions the large SEP events in this study were all generated by CMEs with the presence of enhanced type II radio bursts, which also have wider longitudinal distributions comparing to events with the absence of type II radio burst or its enhancement (almost associated with small SEP events). It seems that the signature of type II radio bursts enhancement is a good discriminator between large SEP and small or none SEP event producers during CME interactions. The type II radio burst enhancement is more likely to be generated by CME interactions, with the main CME having larger speed (v), angular width (WD), mass (m) and kinetic energy (E_k), that taking over the preceding CMEs which also have higher v, WD, m and E_k, than those preceding CMEs in CME pairs missing the type II radio bursts or enhancements. Our analysis also revealed that the intensities of associated SEP events correlate negatively with the intersection height of the two CMEs. Most of type-II-enhanced events and SEP events are coincidentally and almost always made by the fast and wide main CMEs that sweeping fully over a relatively slower and narrower preceding CMEs. We suggest that a fast CME with enough energy completely overtaking a relatively narrower preceding CME, especially in low height, can drive more energetic shock signified by the enhanced type II radio bursts. The shock may accelerate ambient particles and lead to large SEP event more easily.


title: X, Y, Z Search at Belle II

abstract: Search for exotics has increased importance since the observation of the X(3872), 13 years ago, announced by the Belle Collaboration. The observation of pentaquark states by LHCb, and the Z-charged states observed at Belle and BES III have raised even more the attention to the field. Presently several states are observed that do not fit potential models, and looking for them in different production mechanisms and search for their decay modes it is important, as well as to do precise measurement of their mass, width, lineshape. We shortly report in this note about the plan in searching for exotics at Belle II at KEK (Tsukuba, Japan), that just ended the Phase-II running period, and show the first re-discovery results using 5 pb$^{-1}$ integrated luminosity.


title: Linearization Instability in Gravity Theories

abstract: In a nonlinear theory, such as gravity, physically relevant solutions are usually hard to find. Therefore, starting from a background exact solution with symmetries, one uses the perturbation theory, which albeit approximately, provides a lot of information regarding a physical solution. But even this approximate information comes with a price: the basic premise of a perturbative solution is that it should be improvable. Namely, by going to higher order perturbation theory, one should be able to improve and better approximate the physical problem or the solution. While this is often the case in many theories and many background solutions, there are important cases where the linear perturbation theory simply fails for various reasons. This issue is well known in the context of general relativity through the works that started in the early 1970s, but it has only been recently studied in modified gravity theories. This thesis is devoted to the study of linearization instability in generic gravity theories where there are spurious solutions to the linearized equations which do not come from the linearization of possible exact solutions. For this purpose we discuss the Taub charges, the ADT charges and the quadratic constraints on the linearized solutions. We give the three dimensional chiral gravity and the $D$ dimensional critical gravity as explicit examples and give a detailed ADM analysis of the topologically massive gravity with a cosmological constant.


title: Asynchronous Sequential Inertial Iterations for Common Fixed Points Problems with an Application to Linear Systems

abstract: The common fixed points problem requires finding a point in the intersection of fixed points sets of a finite collection of operators. Quickly solving problems of this sort is of great practical importance for engineering and scientific tasks (e.g., for computed tomography). Iterative methods for solving these problems often employ a Krasnosel'ski\u{\i}-Mann type iteration. We present an Asynchronous Sequential Inertial (ASI) algorithmic framework in a Hilbert space to solve common fixed points problems with a collection of nonexpansive operators. Our scheme allows use of out-of-date iterates when generating updates, thereby enabling processing nodes to work simultaneously and without synchronization. This method also includes inertial type extrapolation terms to increase the speed of convergence. In particular, we extend the application of the recent ARock algorithm" [Peng, Z. et al, SIAM J. on Scientific Computing 38, A2851-2879, (2016)] in the context of convex feasibility problems. Convergence of the ASI algorithm is proven with no assumption on the distribution of delays, except that they be uniformly bounded. Discussion is provided along with a computational example showing the performance of the ASI algorithm applied in conjunction with a diagonally relaxed orthogonal projections (DROP) algorithm for estimating solutions to large linear systems.


title: Mentor Protege Matchmaking Framework for IT Entrepreneurs: A Qualitative Study

abstract: This study proposes a framework for matchmaking a mentor and mentee by exploring intentions of entrepreneurs towards mentoring. The framework based on Theory of Planned Behaviour, Institutional Theory and Social Exchange Theory were applied to explore factors that influenced intentions of IT-Entrepreneurs in Malaysia towards mentoring. It was applied at three different phases to understand how selection of mentor changes from one phase to another, and identify type of mentoring support that entrepreneurs need during each phase. Qualitative data were collected from entrepreneurs through semi-structured interviews and analysed using NVivo10. Findings showed that entrepreneurs need different mentoring support and skills when they go through each phase of their entrepreneurial career. The conception phase is the most important phase because entrepreneurs are still in the process of developing their product and want a mentor to help them to evaluate the viability of their product. Entrepreneurs at this phase are also more concerned about the surface level characteristics of mentor such as age, gender, race and language than the other phases. At the start-up phase, entrepreneurs want a mentor who can provide them networking support and help them in decision making. Besides that, the IT entrepreneurs think that mentoring services provided in the incubators are very generic and not helpful. On the contrary, entrepreneurs at the growth phase want a mentor whom they see as a role model and is more experienced and successful. Entrepreneurs at this phase do not care much about factors such as gender, race and religion of the mentor.


title: Dynamic programming for optimal stopping via pseudo-regression

abstract: We introduce new variants of classical regression-based algorithms for optimal stopping problems based on computation of regression coefficients by Monte Carlo approximation of the corresponding $L^2$ inner products instead of the least-squares error functional. Coupled with new proposals for simulation of the underlying samples, we call the approach "pseudo regression". A detailed convergence analysis is provided and it is shown that the approach asymptotically leads to less computational cost for a pre-specified error tolerance, hence to lower complexity. The method is justified by numerical examples.


title: Irreducible binary cubics and the generalized superelliptic equation over number fields

abstract: For a large class (heuristically most) of irreducible binary cubic forms $F(x,y) \in \mathbb Z[x,y]$, Bennett and Dahmen proved that the generalized superelliptic equation $F(x,y)=z^l$ has at most finitely many solutions in $x,y \in \mathbb Z$ coprime, $z \in \mathbb Z$ and exponent $l \in \mathbb Z_{\geq 4} $. Their proof uses, among other ingredients, modularity of certain mod $l$ Galois representations and Ribet's level lowering theorem. The aim of this paper is to treat the same problem for binary cubics with coefficients in $\mathcal O_K$, the ring of integers of an arbitrary number field $K$, using by now well-documented modularity conjectures.


title: Momentum chromospin classification of the fundamental fermions and bosons

abstract: We present a momentum chromospin scheme, the $\mathsf{KSC}$-{\it framework}, for parametrizing the elementary particles using the underlying connection between their internal symmetries and their electromagnetic field boundaries. The linear momentum (or $\mathsf{K}$) chromospin branch of the framework deals with electric charge information, including the emergence of quantum chromodynamics color from fractional electric charges and its effect on quark confinement, while the angular momentum (or $\mathsf{S}$) chromospin branch of the framework develops the theory of spin to address the issue of three fermion generations.


title: CosmoFlow: Using Deep Learning to Learn the Universe at Scale

abstract: Deep learning is a promising tool to determine the physical model that describes our universe. To handle the considerable computational cost of this problem, we present CosmoFlow: a highly scalable deep learning application built on top of the TensorFlow framework. CosmoFlow uses efficient implementations of 3D convolution and pooling primitives, together with improvements in threading for many element-wise operations, to improve training performance on Intel(C) Xeon Phi(TM) processors. We also utilize the Cray PE Machine Learning Plugin for efficient scaling to multiple nodes. We demonstrate fully synchronous data-parallel training on 8192 nodes of Cori with 77% parallel efficiency, achieving 3.5 Pflop/s sustained performance. To our knowledge, this is the first large-scale science application of the TensorFlow framework at supercomputer scale with fully-synchronous training. These enhancements enable us to process large 3D dark matter distribution and predict the cosmological parameters $\Omega_M$, $\sigma_8$ and n$_s$ with unprecedented accuracy.

categories: astro-ph.CO astro-ph.IM cs.LG physics.comp-ph


title: Red and Reddened: Ultraviolet through Near-Infrared Observations of Type Ia Supernova 2017erp

abstract: We present space-based ultraviolet/optical photometry and spectroscopy with the Swift Ultra-Violet/Optical Telescope and Hubble Space Telescope, respectively, along with ground-based optical photometry and spectroscopy and near-infrared spectroscopy of supernova SN2017erp. The optical light curves and spectra are consistent with a normal Type Ia supernova (SN Ia). Compared to previous photometric samples in the near-ultraviolet (NUV), SN2017erp has colors similar to the NUV-red category after correcting for Milky Way and host dust reddening. We find the difference between SN2017erp and the NUV-blue SN2011fe is not consistent with dust reddening alone but is similar to the SALT color law, derived from rest-frame UV photometry of higher redshift SNe Ia. This chromatic difference is dominated by the intrinsic differences in the UV and only a small contribution from the expected dust reddening. Differentiating the two can have important consequences for determining cosmological distances with rest-frame UV photometry. This spectroscopic series is important for analyzing SNe Ia with intrinsically redder NUV colors. We also show model comparisons suggesting that metallicity could be the physical difference between NUV-blue and NUV-red SNe Ia, with emission peaks from reverse fluorescence near 3000 Angstroms implying a factor of ten higher metallicity in the upper layers of SN2017erp compared to SN~2011fe. Metallicity estimates are very model dependent however, and there are multiple effects in the UV. Further models and UV spectra of SNe Ia are needed to explore the diversity of SNe Ia which show seemingly independent differences in the near-UV peaks and mid-UV flux levels.


title: Analyzing Inverse Problems with Invertible Neural Networks

abstract: In many tasks, in particular in natural science, the goal is to determine hidden system parameters from a set of measurements. Often, the forward process from parameter- to measurement-space is a well-defined function, whereas the inverse problem is ambiguous: one measurement may map to multiple different sets of parameters. In this setting, the posterior parameter distribution, conditioned on an input measurement, has to be determined. We argue that a particular class of neural networks is well suited for this task -- so-called Invertible Neural Networks (INNs). Although INNs are not new, they have, so far, received little attention in literature. While classical neural networks attempt to solve the ambiguous inverse problem directly, INNs are able to learn it jointly with the well-defined forward process, using additional latent output variables to capture the information otherwise lost. Given a specific measurement and sampled latent variables, the inverse pass of the INN provides a full distribution over parameter space. We verify experimentally, on artificial data and real-world problems from astrophysics and medicine, that INNs are a powerful analysis tool to find multi-modalities in parameter space, to uncover parameter correlations, and to identify unrecoverable parameters.


title: First-order N\'eel-cVBS transition in a model square lattice $S=1$ antiferromagnet

abstract: We study the N\'eel to four-fold columnar valence bond solid (cVBS) quantum phase transition in a sign free $S=1$ square lattice model. This is the same kind of transition that for $S=1/2$ has been argued to realize the prototypical deconfined critical point. Extensive numerical simulations of the square lattice $S=1/2$ N\'eel-VBS transition have found consistency with the DCP scenario with no direct evidence for first order behavior. In contrast to the $S=1/2$ case, in our quantum Monte Carlo simulations for the $S=1$ model, we present unambiguous evidence for a direct conventional first-order quantum phase transition. Classic signs for a first order transition demonstrating co-existence including double peaked histograms and switching behavior are observed. The sharp contrast from the $S=1/2$ case is remarkable, and is a striking demonstration of the role of the size of the quantum spin in the phase diagram of two dimensional lattice models.


title: A model of second-order arithmetic satisfying AC but not DC

abstract: We show that there is a $\beta$-model of second-order arithmetic in which the choice scheme holds, but the dependent choice scheme fails for a $\Pi^1_2$-assertion, confirming a conjecture of Stephen Simpson. We obtain as a corollary that the Reflection Principle, stating that every formula reflects to a transitive set, can fail in models of ${\rm ZFC}^-$. This work is a rediscovery by the first two authors of a result obtained by the third author.


title: Artificial Intelligence for High-Throughput Discovery of Topological Insulators: the Example of Alloyed Tetradymites

abstract: Significant advances have been made in predicting new topological materials using high-throughput empirical descriptors or symmetry-based indicators. To date, these approaches have been applied to materials in existing databases, and are severely limited to systems with well-defined symmetries, leaving a much larger materials space unexplored. Using tetradymites as a prototypical class of examples, we uncover a novel two-dimensional descriptor by applying an artificial intelligence (AI) based approach for fast and reliable identification of the topological characters of a drastically expanded range of materials, without prior determination of their specific symmetries and detailed band structures. By leveraging this descriptor that contains only the atomic number and electronegativity of the constituent species, we have readily scanned a huge number of alloys in the tetradymite family. Strikingly, nearly half of which are identified to be topological insulators, revealing a much larger territory of the topological materials world. The present work also attests the increasingly important role of such AI-based approaches in modern materials discovery.


title: Optimal probabilities and controls for reflecting diffusion processes

abstract: A solution to the optimal problem for determining vector fields which maximize (resp. minimize) the transition probabilities from one location to another for a class of reflecting diffusion processes is obtained in the present paper. The approach is based on a representation for the transition probability density functions. The optimal transition probabilities under the constraint that the drift vector field is bounded by a constant are studied in terms of the HJB equation. In dimension one, the optimal reflecting diffusion processes and the bang-bang diffusion processes are considered. We demonstrate by simulations that, even in this special case, the optimal diffusion processes exhibit an interesting feature of phase transitions. We also solve an optimal stochastic control problem for a class of stochastic control problems involving diffusion processes with reflection.


title: On the removal of the trace mode in lattice ${\cal N}=4$ super Yang-Mills theory

abstract: Twisted and orbifold formulations of lattice ${\cal N}=4$ super Yang-Mills theory which possess an exact supersymmetry require a $U(N)=SU(N)\otimes U(1)$ gauge group. In the naive continuum limit, the $U(1)$ modes trivially decouple and play no role in the theory. However, at non-zero lattice spacing they couple to the $SU(N)$ modes and can drive instabilities in the lattice theory. For example, it is well known that the lattice $U(1)$ theory undergoes a phase transition at strong coupling to a chirally broken phase. An improved action that suppresses the fluctuations in the $U(1)$ sector was proposed in arXiv:1505.03135 . Here, we explore a more aggressive approach to the problem by adding a term to the action which can entirely suppress the $U(1)$ mode. The penalty is that the new term breaks the $\mathcal{Q}$-exact lattice supersymmetry. However, we argue that the term is $1/N^2$ suppressed and the existence of a supersymmetric fixed point in the planar limit ensures that any SUSY-violating terms induced in the action possess couplings that also vanish in this limit. We present numerical results on supersymmetric Ward identities consistent with this conclusion.


title: Adversarial Neural Networks for Cross-lingual Sequence Tagging

abstract: We study cross-lingual sequence tagging with little or no labeled data in the target language. Adversarial training has previously been shown to be effective for training cross-lingual sentence classifiers. However, it is not clear if language-agnostic representations enforced by an adversarial language discriminator will also enable effective transfer for token-level prediction tasks. Therefore, we experiment with different types of adversarial training on two tasks: dependency parsing and sentence compression. We show that adversarial training consistently leads to improved cross-lingual performance on each task compared to a conventionally trained baseline.


title: Uniqueness and Lipschitz stability in Electrical Impedance Tomography with finitely many electrodes

abstract: For the linearized reconstruction problem in Electrical Impedance Tomography (EIT) with the Complete Electrode Model (CEM), Lechleiter and Rieder (2008 Inverse Problems 24 065009) have shown that a piecewise polynomial conductivity on a fixed partition is uniquely determined if enough electrodes are being used. We extend their result to the full non-linear case and show that measurements on a sufficiently high number of electrodes uniquely determine a conductivity in any finite-dimensional subset of piecewise-analytic functions. We also prove Lipschitz stability, and derive analogue results for the continuum model, where finitely many measurements determine a finite-dimensional Galerkin projection of the Neumann-to-Dirichlet operator on a boundary part.


title: Stream Reasoning on Expressive Logics

abstract: Data streams occur widely in various real world applications. The research on streaming data mainly focuses on the data management, query evaluation and optimization on these data, however the work on reasoning procedures for streaming knowledge bases on both the assertional and terminological levels is very limited. Typically reasoning services on large knowledge bases are very expensive, and need to be applied continuously when the data is received as a stream. Hence new techniques for optimizing this continuous process is needed for developing efficient reasoners on streaming data. In this paper, we survey the related research on reasoning on expressive logics that can be applied to this setting, and point to further research directions in this area.


title: Simulating Markov random fields with a conclique-based Gibbs sampler

abstract: For spatial and network data, we consider models formed from a Markov random field (MRF) structure and the specification of a conditional distribution for each observation. Fast simulation from such MRF models is often an important consideration, particularly when repeated generation of large numbers of data sets is required. However, a standard Gibbs strategy for simulating from MRF models involves single-site updates, performed with the conditional univariate distribution of each observation in a sequential manner, whereby a complete Gibbs iteration may become computationally involved even for moderate samples. As an alternative, we describe a general way to simulate from MRF models using Gibbs sampling with "concliques" (i.e., groups of non-neighboring observations). Compared to standard Gibbs sampling, this simulation scheme can be much faster by reducing Gibbs steps and independently updating all observations per conclique at once. The speed improvement depends on the number of concliques relative to the sample size for simulation, and order-of-magnitude speed increases are possible with many MRF models (e.g., having appropriately bounded neighborhoods). We detail the simulation method, establish its validity, and assess its computational performance through numerical studies, where speed advantages are shown for several spatial and network examples.


title: Optimization and projection of coated structures with orthotropic infill material

abstract: The purpose of this work is two-fold. First, we introduce an efficient homogenization-based approach to perform topology optimization of coated structures with orthotropic infill material. By making use of the relaxed design space, we can obtain designs with complex microstructures on a relatively coarse mesh. In the second part of this work, a method is presented to project the homogenization-based designs on a fine but realizable scale. A novel method to adaptively refine the lattice structure is presented to allow for a regular spacing of the infill. Numerical experiments show excellent behavior of the projected designs, which perform almost identical to the homogenization-based designs. Furthermore, a reduction in computational cost of at least an order of magnitude is achieved, compared to a related approach in which the infill is optimized using a density-based approach.


title: A Novel Method for Determining DOA from Far-Field TDOA or FDOA

abstract: Passive source localization is often performed using time difference of arrival (TDOA) measurements, frequency difference of arrival (FDOA) measurements, direction of arrival (DOA) measurements, or a combination of all of these. For a source in the far-field, DOA can be extracted from the TDOA and FDOA measurements due to simplifications that arise in the far-field approximation. This paper presents this relationship and the corresponding DOA estimation method. Utilizing TDOA and FDOA measurements for computation of signal DOA requires only a linear solve, which makes the corresponding source localization technique very efficient. Additionally, the method provides an inherent de-noising of receiver measurements, since they are being projected onto the range of the receiver differencing matrix.


title: Anomalous dimension of subleading-power ${N}$-jet operators II

abstract: We continue the investigation of the anomalous dimension of subleading-power $N$-jet operators. In this paper, we focus on the operators with fermion number one in each collinear direction, corresponding to quark (antiquark) initiated jets in QCD. We investigate the renormalization effects induced by the soft loop and compute the one-loop mixing of time-ordered products involving power-suppressed SCET Lagrangian insertions into $N$-jet currents through soft loops. We discuss fermion number conservation in collinear directions and provide explicit results for the collinear anomalous dimension matrix of the currents. The Feynman rules for the power-suppressed SCET interactions in the position-space formalism are collected in an appendix.


title: Derivative Corrections to the Trapezoidal Rule

abstract: Extensions to the trapezoidal rule using derivative information are studied for periodic integrands and integrals along the entire real line. Integrands which are analytic within a half plane or within a strip containing the path of integration are considered. Derivative-free error bounds are obtained. Alternative approaches to including derivative information are discussed.


title: Spanning Properties of Theta-Theta-6

abstract: We show that, unlike the Yao-Yao graph $YY_6$, the Theta-Theta graph $\Theta\Theta_6$ defined by six cones is a spanner for sets of points in convex position. We also show that, for sets of points in non-convex position, the spanning ratio of $\Theta\Theta_6$ is unbounded.


title: Imagining the Unseen: Learning a Distribution over Incomplete Images with Dense Latent Trees

abstract: Images are composed as a hierarchy of object parts. We use this insight to create a generative graphical model that defines a hierarchical distribution over image parts. Typically, this leads to intractable inference due to loops in the graph. We propose an alternative model structure, the Dense Latent Tree (DLT), which avoids loops and allows for efficient exact inference, while maintaining a dense connectivity between parts of the hierarchy. The usefulness of DLTs is shown for the example task of image completion on partially observed MNIST and Fashion-MNIST data. We verify having successfully learned a hierarchical model of images by visualising its latent states.


title: Neutron star -- axion star collisions in the light of multi-messenger astronomy

abstract: Axions are increasingly favoured as a candidate particle for the dark matter in galaxies, since they satisfy the observational requirements for cold dark matter and are theoretically well motivated. Fluctuations in the axion field give rise to stable localised overdensities known as axion stars, which, for the most massive, compact cases, are potential neutron star mimickers. In principle, there are no fundamental arguments against the multi-messenger observations of GW170817/GRB170817A/AT2017gfo arising from the merger of a neutron star with a neutron star mimicker, rather than from a binary neutron star. To constrain this possibility and better understand the astrophysical signatures of a neutron star--axion star (NSAS) merger, we present in this work a detailed example case of a NSAS merger based on full 3D numerical relativity simulations, and give an overview of the many potential observables - ranging from gravitational waves, to optical and near-infrared electromagnetic signals, radio flares, fast radio bursts, gamma ray bursts, and neutrino emission. We discuss the individual channels and estimate to which distances current and future observatories might be able to detect such a NSAS merger. Such signals could constrain the unknown axion mass and its couplings to standard baryonic matter,thus enhancing our understanding of the dark matter sector of the Universe.


title: A penalty scheme for monotone systems with interconnected obstacles: convergence and error estimates

abstract: We present a novel penalty approach for a class of quasi-variational inequalities (QVIs) involving monotone systems and interconnected obstacles. We show that for any given positive switching cost, the solutions of the penalized equations converge monotonically to those of the QVIs. We estimate the penalization errors and are able to deduce that the optimal switching regions are constructed exactly. We further demonstrate that as the switching cost tends to zero, the QVI degenerates into an equation of HJB type, which is approximated by the penalized equation at the same order (up to a log factor) as that for positive switching cost. Numerical experiments on optimal switching problems are presented to illustrate the theoretical results and to demonstrate the effectiveness of the method.


title: Screening potential topological insulators in half-Heusler compounds via compressed-sensing

abstract: Ternary half-Heusler compounds with widely tunable electronic structures, present a new platform to discover topological insulators. Due to time-consuming computations and synthesis procedures, the identification of new topological insulators is however a rough task. Here, we adopt a compressed-sensing approach to rapidly screen potential topological insulators in half-Heusler family, which is realized via a two-dimensional descriptor that only depends on the fundamental properties of the constituent atoms. Beyond the finite training data, the proposed descriptor is employed to screen many new half-Heusler compounds, including those with integer and fractional stoichiometry, and a larger number of possible topological insulators are predicted.


title: Proton Acceleration in a Laser-induced Relativistic Electron Vortex

abstract: We show that when a solid plasma foil with a density gradient on the front surface is irradiated by an intense laser pulse at a grazing angle, around 80 degrees, a relativistic electron vortex is excited in the near-critical-density layer after the laser pulse depletion. The vortex structure and dynamics are studied using particle-in-cell simulations. Due to the asymmetry introduced by nonuniform background density, the vortex drifts at a constant velocity, typically 0.2 to 0.3 times the speed of light. The strong magnetic field inside the vortex leads to significant charge separation; in the corresponding electric field initially stationary protons can be captured and accelerated to twice the velocity of the vortex (100-200 MeV). A representative scenario - with laser intensity of 10^21 W/cm^2 -is discussed: two dimensional simulations suggest that a quasi-monoenergetic proton beam can be obtained with a mean energy 140 MeV and an energy spread of about 10%. We derive an analytical estimate for the vortex velocity in terms of laser and plasma parameters, demonstrating that the maximum proton energy can be controlled by the incidence angle of the laser and the plasma density gradient.


title: Quantifying the Influences on Probabilistic Wind Power Forecasts

abstract: In recent years, probabilistic forecasts techniques were proposed in research as well as in applications to integrate volatile renewable energy resources into the electrical grid. These techniques allow decision makers to take the uncertainty of the prediction into account and, therefore, to devise optimal decisions, e.g., related to costs and risks in the electrical grid. However, it was yet not studied how the input, such as numerical weather predictions, affects the model output of forecasting models in detail. Therefore, we examine the potential influences with techniques from the field of sensitivity analysis on three different black-box models to obtain insights into differences and similarities of these probabilistic models. The analysis shows a considerable number of potential influences in those models depending on, e.g., the predicted probability and the type of model. These effects motivate the need to take various influences into account when models are tested, analyzed, or compared. Nevertheless, results of the sensitivity analysis will allow us to select a model with advantages in the practical application.


title: The Active Assembly of the Virgo Cluster: Indications for Recent Group Infall From Early-Type Dwarf Galaxies

abstract: Virgo is a dynamically young galaxy cluster with substructure in its spatial and kinematic distribution. Here, we simultaneously study the phase-space distribution and the main characteristics of Virgo's galaxies, particularly its most abundant galaxy population - the early-type dwarfs - to understand their environmental transformation histories. Aside from known correlations with morphological types - like the larger average clustercentric distance of late-type galaxies - we find an intriguing behavior of early types with magnitudes -17 >= M_r >= -18. They show a large velocity spread and an asymmetric phase-space distribution, similar to the late-type galaxies and different from the early types just one magnitude brighter/fainter. Furthermore, we find a close phase-space aggregation of early-type dwarfs at large clustercentric distance and high relative velocity. Nearly all of them show signatures of disk components and their colors imply stellar ages that are younger than the population average. They are not located closely together but spread azimuthally around the cluster center. We show that this is expected from simulations of an infalling galaxy group that slowly gets dispersed after its first pericentric passage. We thus conclude that these galaxies are recent arrivals, and that the peculiar phase-space distribution of early-type dwarfs is evidence for the ongoing growth of this galaxy population. Studying galaxies based on their phase space correlations is a unique way to compare the properties of recent and older cluster members, and to understand which environment most influenced their present-day characteristics.


title: A Survey on Methods and Theories of Quantized Neural Networks

abstract: Deep neural networks are the state-of-the-art methods for many real-world tasks, such as computer vision, natural language processing and speech recognition. For all its popularity, deep neural networks are also criticized for consuming a lot of memory and draining battery life of devices during training and inference. This makes it hard to deploy these models on mobile or embedded devices which have tight resource constraints. Quantization is recognized as one of the most effective approaches to satisfy the extreme memory requirements that deep neural network models demand. Instead of adopting 32-bit floating point format to represent weights, quantized representations store weights using more compact formats such as integers or even binary numbers. Despite a possible degradation in predictive performance, quantization provides a potential solution to greatly reduce the model size and the energy consumption. In this survey, we give a thorough review of different aspects of quantized neural networks. Current challenges and trends of quantized neural networks are also discussed.


title: Estimating the size of a hidden finite set: large-sample behavior of estimators

abstract: A finite set is "hidden" if its elements are not directly enumerable or if its size cannot be ascertained via a deterministic query. In public health, epidemiology, demography, ecology and intelligence analysis, researchers have developed a wide variety of indirect statistical approaches, under different models for sampling and observation, for estimating the size of a hidden set. Some methods make use of random sampling with known or estimable sampling probabilities, and others make structural assumptions about relationships (e.g. ordering or network information) between the elements that comprise the hidden set. In this review, we describe models and methods for learning about the size of a hidden finite set, with special attention to asymptotic properties of estimators. We study the properties of these methods under two asymptotic regimes, "infill" in which the number of fixed-size samples increases, but the population size remains constant, and "outfill" in which the sample size and population size grow together. Statistical properties under these two regimes can be dramatically different.


title: Treepedia 2.0: Applying Deep Learning for Large-scale Quantification of Urban Tree Cover

abstract: Recent advances in deep learning have made it possible to quantify urban metrics at fine resolution, and over large extents using street-level images. Here, we focus on measuring urban tree cover using Google Street View (GSV) images. First, we provide a small-scale labelled validation dataset and propose standard metrics to compare the performance of automated estimations of street tree cover using GSV. We apply state-of-the-art deep learning models, and compare their performance to a previously established benchmark of an unsupervised method. Our training procedure for deep learning models is novel; we utilize the abundance of openly available and similarly labelled street-level image datasets to pre-train our model. We then perform additional training on a small training dataset consisting of GSV images. We find that deep learning models significantly outperform the unsupervised benchmark method. Our semantic segmentation model increased mean intersection-over-union (IoU) from 44.10% to 60.42% relative to the unsupervised method and our end-to-end model decreased Mean Absolute Error from 10.04% to 4.67%. We also employ a recently developed method called gradient-weighted class activation map (Grad-CAM) to interpret the features learned by the end-to-end model. This technique confirms that the end-to-end model has accurately learned to identify tree cover area as key features for predicting percentage tree cover. Our paper provides an example of applying advanced deep learning techniques on a large-scale, geo-tagged and image-based dataset to efficiently estimate important urban metrics. The results demonstrate that deep learning models are highly accurate, can be interpretable, and can also be efficient in terms of data-labelling effort and computational resources.


title: Entanglement of neutral-atom qubits with long ground-Rydberg coherence times

abstract: We report results of a ground-state entanglement protocol for a pair of Cs atoms separated by 6~$\mu$m, combining the Rydberg blockade mechanism with a two-photon Raman transitions to prepare the $\vert\Psi^+\rangle=(\vert 10\rangle+\vert 01\rangle)/\sqrt{2}$ Bell state with a loss-corrected fidelity of 0.81(5), equal to the best demonstrated fidelity for atoms trapped in optical tweezers but without the requirement for dynamically adjustable interatomic spacing. Qubit state coherence is also critical for quantum information applications, and we characterise both ground-state and ground-Rydberg dephasing rates using Ramsey spectroscopy. We demonstrate transverse dephasing times $T_2^*=10(1)$~ms and $T_2'=0.14(1)$~s for the qubit levels and achieve long ground-Rydberg coherence times of $T_2^*=17(2)~\mu$s as required for implementing high-fidelity multi-qubit gate sequences where a control atom remains in the Rydberg state while applying local operations on neighbouring target qubits.


title: Micromechanical theory of strain-stiffening of biopolymer networks

abstract: Filamentous bio-materials such as fibrin or collagen networks exhibit an enormous stiffening of their elastic moduli upon large deformations. This pronounced nonlinear behavior stems from a significant separation between the stiffnesses scales associated with bending vs. stretching the material's constituent elements. Here we study a simple model of such materials - floppy networks of hinged rigid bars embedded in an elastic matrix - in which the effective ratio of bending to stretching stiffnesses vanishes identically. We introduce a theoretical framework and build upon it to construct a numerical method with which the model's micro- and macro-mechanics can be carefully studied. Our model, numerical method and theoretical framework allow us to robustly observe and fully understand the critical properties of the athermal strain-stiffening transition that underlies the nonlinear mechanical response of a broad class of biomaterials.


title: Addressing Johnson graphs, complete multipartite graphs, odd cycles and other graphs

abstract: Graham and Pollak showed that the vertices of any graph $G$ can be addressed with $N$-tuples of three symbols, such that the distance between any two vertices may be easily determined from their addresses. An addressing is optimal if its length $N$ is minimum possible. In this paper, we determine an addressing of length $k(n-k)$ for the Johnson graphs $J(n,k)$ and we show that our addressing is optimal when $k=1$ or when $k=2, n=4,5,6$, but not when $n=6$ and $k=3$. We study the addressing problem as well as a variation of it in which the alphabet used has more than three symbols, for other graphs such as complete multipartite graphs and odd cycles. We also present computations describing the distribution of the minimum length of addressings for connected graphs with up to $10$ vertices. Motivated by these computations we settle a problem of Graham, showing that most graphs on $n$ vertices have an addressing of length at most $n-(2-o(1))\log_2 n$.


title: Finding Minimal Cost Herbrand Models with Branch-Cut-and-Price

abstract: Given (1) a set of clauses $T$ in some first-order language $\cal L$ and (2) a cost function $c : B_{{\cal L}} \rightarrow \mathbb{R}_{+}$, mapping each ground atom in the Herbrand base $B_{{\cal L}}$ to a non-negative real, then the problem of finding a minimal cost Herbrand model is to either find a Herbrand model $\cal I$ of $T$ which is guaranteed to minimise the sum of the costs of true ground atoms, or establish that there is no Herbrand model for $T$. A branch-cut-and-price integer programming (IP) approach to solving this problem is presented. Since the number of ground instantiations of clauses and the size of the Herbrand base are both infinite in general, we add the corresponding IP constraints and IP variables `on the fly' via `cutting' and `pricing' respectively. In the special case of a finite Herbrand base we show that adding all IP variables and constraints from the outset can be advantageous, showing that a challenging Markov logic network MAP problem can be solved in this way if encoded appropriately.


title: An Overview and a Benchmark of Active Learning for Outlier Detection with One-Class Classifiers

abstract: Active learning methods increase classification quality by means of user feedback. An important subcategory is active learning for outlier detection with one-class classifiers. While various methods in this category exist, selecting one for a given application scenario is difficult. This is because existing methods rely on different assumptions, have different objectives, and often are tailored to a specific use case. All this calls for a comprehensive comparison, the topic of this article. This article starts with a categorization of the various methods. We then propose ways to evaluate active learning results. Next, we run extensive experiments to compare existing methods, for a broad variety of scenarios. Based on our results, we formulate guidelines on how to select active learning methods for outlier detection with one-class classifiers.


title: Parallel Statistical and Machine Learning Methods for Estimation of Physical Load

abstract: Several statistical and machine learning methods are proposed to estimate the type and intensity of physical load and accumulated fatigue . They are based on the statistical analysis of accumulated and moving window data subsets with construction of a kurtosis-skewness diagram. This approach was applied to the data gathered by the wearable heart monitor for various types and levels of physical activities, and for people with various physical conditions. The different levels of physical activities, loads, and fitness can be distinguished from the kurtosis-skewness diagram, and their evolution can be monitored. Several metrics for estimation of the instant effect and accumulated effect (physical fatigue) of physical loads were proposed. The data and results presented allow to extend application of these methods for modeling and characterization of complex human activity patterns, for example, to estimate the actual and accumulated physical load and fatigue, model the potential dangerous development, and give cautions and advice in real time.


title: Cache Telepathy: Leveraging Shared Resource Attacks to Learn DNN Architectures

abstract: Deep Neural Networks (DNNs) are fast becoming ubiquitous for their ability to attain good accuracy in various machine learning tasks. A DNN's architecture (i.e., its hyper-parameters) broadly determines the DNN's accuracy and performance, and is often confidential. Attacking a DNN in the cloud to obtain its architecture can potentially provide major commercial value. Further, attaining a DNN's architecture facilitates other, existing DNN attacks. This paper presents Cache Telepathy: a fast and accurate mechanism to steal a DNN's architecture using the cache side channel. Our attack is based on the insight that DNN inference relies heavily on tiled GEMM (Generalized Matrix Multiply), and that DNN architecture parameters determine the number of GEMM calls and the dimensions of the matrices used in the GEMM functions. Such information can be leaked through the cache side channel. This paper uses Prime+Probe and Flush+Reload to attack VGG and ResNet DNNs running OpenBLAS and Intel MKL libraries. Our attack is effective in helping obtain the architectures by very substantially reducing the search space of target DNN architectures. For example, for VGG using OpenBLAS, it reduces the search space from more than $10^{35}$ architectures to just 16.


title: The orientation of galaxy pairs with filamentary structures: dependence on morphology

abstract: Aims. With the aim of performing an analysis of the orientations of galaxy pair systems with respect to the underlying large-scale structure, we study the alignment between the axis connecting the pair galaxies and the host cosmic filament where the pair resides. In addition, we analyze the dependence of the amplitude of the alignment on the morphology of pair members as well as filament properties. Methods. We build a galaxy pair catalog requiring $r_p < 100\kpc$ and $\Delta V < 500 \kms $ within redshift $z<0.1$ from the Sloan Digital Sky Survey (SDSS). We divided the galaxy pair catalog taking into account the morphological classification by defining three pair categories composed by elliptical-elliptical (E-E), elliptical-spiral (E-S) and spiral-spiral (S-S) galaxies. We use a previously defined catalog of filaments obtained from SDSS and we select pairs located closer than $1\mpc$ from the filament spine, which are considered as members of filaments. For these pairs, we calculate the relative angle between the axis connecting each galaxy, and the direction defined by the spine of the parent filament. Results. We find a statistically significant alignment signal between the pair axes and the spine of the host filaments consistent with a relative excess of $ \sim$ 15\% aligned pairs. We obtain that pairs composed by elliptical galaxies exhibit a stronger alignment, showing a higher alignment signal for pairs closer than 200 $\kpc$ to the filament spine. In addition, we find that the aligned pairs are associated with luminous host filaments populated with a high fraction of elliptical galaxies. The findings of this work show that large scale structures play a fundamental role in driving galactic anisotropic accretion as induced by galaxy pairs exhibiting a preferred alignment along the filament direction.


title: Some lower bounds for solutions of Schrodinger evolutions

abstract: We present some lower bounds for regular solutions of Schr\"odinger equations with bounded and time dependent complex potentials. Assuming that the solution has some positive mass at time zero within a ball of certain radius, we prove that this mass can be observed if one looks at the solution and its gradient in space-time regions outside of that ball.


title: Periodic, quasi-periodic, fractal, Kolakoski and random binary polymers: Energy structure and carrier transport

abstract: We study periodic, quasi-periodic (Thue-Morse, Fibonacci, Period Doubling, Rudin-Shapiro), fractal (Cantor, generalized Cantor), Kolakoski and random binary sequences using a tight-binding wire model, where a site is a monomer (e.g., in DNA, a base pair). We use B-DNA as our prototype system. All sequences have purines, guanine (G) or adenine (A) on the same strand, i.e., our prototype binary alphabet is (G,A). Our aim is to examine the influence of sequence intricacy and magnitude of parameters on energy structure, localization and charge transport. We study quantities such as autocorrelation function, eigenspectra, density of states, Lyapunov exponents, transmission coefficients and current-voltage curves. We show that the degree of sequence intricacy and the presence of correlations decisively affect the aforementioned physical properties. Periodic segments have enhanced transport properties. Specifically, in homogeneous sequences transport efficiency is maximum. There are several deterministic aperiodic sequences that can support significant currents, depending on the Fermi level of the leads. Random sequences is the less efficient category.


title: Discrete versus continuous domain models for disease mapping

abstract: The main goal of disease mapping is to estimate disease risk and identify high-risk areas. Such analyses are hampered by the limited geographical resolution of the available data. Typically the available data are counts per spatial unit and the common approach is the Besag--York--Molli{\'e} (BYM) model. When precise geocodes are available, it is more natural to use Log-Gaussian Cox processes (LGCPs). In a simulation study mimicking childhood leukaemia incidence using actual residential locations of all children in the canton of Z\"urich, Switzerland, we compare the ability of these models to recover risk surfaces and identify high-risk areas. We then apply both approaches to actual data on childhood leukaemia incidence in the canton of Z\"urich during 1985-2015. We found that LGCPs outperform BYM models in almost all scenarios considered. Our findings suggest that there are important gains to be made from the use of LGCPs in spatial epidemiology.


title: Early formation of carbon monoxide in the Centaurus A supernova SN 2016adj

abstract: We present near-infrared spectroscopy of the NGC 5128 supernova SN 2016adj in the first 2 months following discovery. We report the detection of first overtone carbon monoxide emission at $\sim58.2$ d after discovery, one of the earliest detections of CO in an erupting supernova. We model the CO emission to derive the CO mass, temperature and velocity, assuming both pure $^{12}$CO and a composition that includes $^{13}$CO; the case for the latter is the isotopic analyses of meteoritic grains, which suggest that core collapse supernovae can synthesise significant amounts of $^{13}$C. Our models show that, while the CO data are adequately explained by pure $^{12}$CO, they do not preclude the presence of $^{13}$CO, to a limit of $^{12}$C/$^{13}$C $> 3$, the first constraint on the $^{12}$C/$^{13}$C ratio determined from near-infrared observations. We estimate the reddening to the object, and the effective temperature from the energy distribution at outburst. We discuss whether the ejecta of SN 2016adj may be carbon-rich, what the infrared data tell us about the classification of this supernova, and what implications the early formation of CO in supernovae may have for CO formation in supernovae in general.