title: Decentralized Task Allocation in Multi-Robot Systems via Bipartite Graph Matching Augmented with Fuzzy Clustering

abstract: Robotic systems, working together as a team, are becoming valuable players in different real-world applications, from disaster response to warehouse fulfillment services. Centralized solutions for coordinating multi-robot teams often suffer from poor scalability and vulnerability to communication disruptions. This paper develops a decentralized multi-agent task allocation (Dec-MATA) algorithm for multi-robot applications. The task planning problem is posed as a maximum-weighted matching of a bipartite graph, the solution of which using the blossom algorithm allows each robot to autonomously identify the optimal sequence of tasks it should undertake. The graph weights are determined based on a soft clustering process, which also plays a problem decomposition role seeking to reduce the complexity of the individual-agents' task assignment problems. To evaluate the new Dec-MATA algorithm, a series of case studies (of varying complexity) are performed, with tasks being distributed randomly over an observable 2D environment. A centralized approach, based on a state-of-the-art MILP formulation of the multi-Traveling Salesman problem is used for comparative analysis. While getting within 7-28% of the optimal cost obtained by the centralized algorithm, the Dec-MATA algorithm is found to be 1-3 orders of magnitude faster and minimally sensitive to task-to-robot ratios, unlike the centralized algorithm.


title: Quantile contours and allometric modelling for risk classification of abnormal ratios with an application to asymmetric growth-restriction in preterm infants

abstract: We develop an approach to risk classification based on quantile contours and allometric modelling of multivariate anthropometric measurements. We propose the definition of allometric direction tangent to the directional quantile envelope, which divides ratios of measurements into half-spaces. This in turn provides an operational definition of directional quantile that can be used as cutoff for risk assessment. We show the application of the proposed approach using a large dataset from the Vermont Oxford Network containing observations of birthweight (BW) and head circumference (HC) for more than 150,000 preterm infants. Our analysis suggests that disproportionately growth-restricted infants with a larger HC-to-BW ratio are at increased mortality risk as compared to proportionately growth-restricted infants. The role of maternal hypertension is also investigated.


title: A Fully Convolutional Neural Network Approach to End-to-End Speech Enhancement

abstract: This paper will describe a novel approach to the cocktail party problem that relies on a fully convolutional neural network (FCN) architecture. The FCN takes noisy audio data as input and performs nonlinear, filtering operations to produce clean audio data of the target speech at the output. Our method learns a model for one specific speaker, and is then able to extract that speakers voice from babble background noise. Results from experimentation indicate the ability to generalize to new speakers and robustness to new noise environments of varying signal-to-noise ratios. A potential application of this method would be for use in hearing aids. A pre-trained model could be quickly fine tuned for an individuals family members and close friends, and deployed onto a hearing aid to assist listeners in noisy environments.


title: Alpha-rooting color image enhancement method by two-side 2-D quaternion discrete Fourier transform followed by spatial transformation

abstract: In this paper a quaternion approach of enhancement method is proposed in which color in the image is considered as a single entity. This new method is referred as the alpha-rooting method of color image enhancement by the two-dimensional quaternion discrete Fourier transform (2-D QDFT) followed by a spatial transformation. The results of the proposed color image enhancement method are compared with its counterpart channel-by-channel enhancement algorithm by the 2-D DFT. The image enhancements are quantified to the enhancement measure that is based on visual perception referred as the color enhancement measure estimation (CEME). The preliminary experiment results show that the quaternion approach of image enhancement is an effective color image enhancement technique.


title: Twitter Sentiment Analysis via Bi-sense Emoji Embedding and Attention-based LSTM

abstract: Sentiment analysis on large-scale social media data is important to bridge the gaps between social media contents and real world activities including political election prediction, individual and public emotional status monitoring and analysis, and so on. Although textual sentiment analysis has been well studied based on platforms such as Twitter and Instagram, analysis of the role of extensive emoji uses in sentiment analysis remains light. In this paper, we propose a novel scheme for Twitter sentiment analysis with extra attention on emojis. We first learn bi-sense emoji embeddings under positive and negative sentimental tweets individually, and then train a sentiment classifier by attending on these bi-sense emoji embeddings with an attention-based long short-term memory network (LSTM). Our experiments show that the bi-sense embedding is effective for extracting sentiment-aware embeddings of emojis and outperforms the state-of-the-art models. We also visualize the attentions to show that the bi-sense emoji embedding provides better guidance on the attention mechanism to obtain a more robust understanding of the semantics and sentiments.


title: A Novel Color Image Enhancement Method by the Transformation of Color Images to 2-D Grayscale Images

abstract: A novel method of color image enhancement is proposed, in which three or four color channels of the image are transformed to one channel 2-D grayscale image. This paper describes different models of such transformations in the RGB and other color models. Color image enhancement is achieved by enhancing first the transformed grayscale image and, then, transforming back the grayscale image into the colors. The color image enhancement is done on the transformed 2-D grayscale image rather than on the color image. New algorithms of color image enhancement are described in both frequency and time domains. The enhancement by this novel method shows good results. The enhancement of the image is measured with respect to the metric referred to as the Color Enhancement Measure Estimation (CEME).


title: Deep Transfer Learning for Cross-domain Activity Recognition

abstract: Human activity recognition plays an important role in people's daily life. However, it is often expensive and time-consuming to acquire sufficient labeled activity data. To solve this problem, transfer learning leverages the labeled samples from the source domain to annotate the target domain which has few or none labels. Unfortunately, when there are several source domains available, it is difficult to select the right source domains for transfer. The right source domain means that it has the most similar properties with the target domain, thus their similarity is higher, which can facilitate transfer learning. Choosing the right source domain helps the algorithm perform well and prevents the negative transfer. In this paper, we propose an effective Unsupervised Source Selection algorithm for Activity Recognition (USSAR). USSAR is able to select the most similar $K$ source domains from a list of available domains. After this, we propose an effective Transfer Neural Network to perform knowledge transfer for Activity Recognition (TNNAR). TNNAR could capture both the time and spatial relationship between activities while transferring knowledge. Experiments on three public activity recognition datasets demonstrate that: 1) The USSAR algorithm is effective in selecting the best source domains. 2) The TNNAR method can reach high accuracy when performing activity knowledge transfer.


title: Question-Aware Sentence Gating Networks for Question and Answering

abstract: Machine comprehension question answering, which finds an answer to the question given a passage, involves high-level reasoning processes of understanding and tracking the relevant contents across various semantic units such as words, phrases, and sentences in a document. This paper proposes the novel question-aware sentence gating networks that directly incorporate the sentence-level information into word-level encoding processes. To this end, our model first learns question-aware sentence representations and then dynamically combines them with word-level representations, resulting in semantically meaningful word representations for QA tasks. Experimental results demonstrate that our approach consistently improves the accuracy over existing baseline approaches on various QA datasets and bears the wide applicability to other neural network-based QA models.


title: An Efficient End-to-End Neural Model for Handwritten Text Recognition

abstract: Offline handwritten text recognition from images is an important problem for enterprises attempting to digitize large volumes of handmarked scanned documents/reports. Deep recurrent models such as Multi-dimensional LSTMs have been shown to yield superior performance over traditional Hidden Markov Model based approaches that suffer from the Markov assumption and therefore lack the representational power of RNNs. In this paper we introduce a novel approach that combines a deep convolutional network with a recurrent Encoder-Decoder network to map an image to a sequence of characters corresponding to the text present in the image. The entire model is trained end-to-end using Focal Loss, an improvement over the standard Cross-Entropy loss that addresses the class imbalance problem, inherent to text recognition. To enhance the decoding capacity of the model, Beam Search algorithm is employed which searches for the best sequence out of a set of hypotheses based on a joint distribution of individual characters. Our model takes as input a downsampled version of the original image thereby making it both computationally and memory efficient. The experimental results were benchmarked against two publicly available datasets, IAM and RIMES. We surpass the state-of-the-art word level accuracy on the evaluation set of both datasets by 3.5% & 1.1%, respectively.


title: Exploring Combinations of Ontological Features and Keywords for Text Retrieval

abstract: Named entities have been considered and combined with keywords to enhance information retrieval performance. However, there is not yet a formal and complete model that takes into account entity names, classes, and identifiers together. Our work explores various adaptations of the traditional Vector Space Model that combine different ontological features with keywords, and in different ways. It shows better performance of the proposed models as compared to the keyword-based Lucene, and their advantages for both text retrieval and representation of documents and queries.


title: Combining Named Entities with WordNet and Using Query-Oriented Spreading Activation for Semantic Text Search

abstract: Purely keyword-based text search is not satisfactory because named entities and WordNet words are also important elements to define the content of a document or a query in which they occur. Named entities have ontological features, namely, their aliases, classes, and identifiers. Words in WordNet also have ontological features, namely, their synonyms, hypernyms, hyponyms, and senses. Those features of concepts may be hidden from their textual appearance. Besides, there are related concepts that do not appear in a query, but can bring out the meaning of the query if they are added. We propose an ontology-based generalized Vector Space Model to semantic text search. It exploits ontological features of named entities and WordNet words, and develops a query-oriented spreading activation algorithm to expand queries. In addition, it combines and utilizes advantages of different ontologies for semantic annotation and searching. Experiments on a benchmark dataset show that, in terms of the MAP measure, our model is 42.5% better than the purely keyword-based model, and 32.3% and 15.9% respectively better than the ones using only WordNet or named entities. Keywords: semantic search, spreading activation, ontology, named entity, WordNet.


title: Scaling in the eigenvalue fluctuations of the empirical correlation matrices

abstract: The spectra of empirical correlation matrices, constructed from multivariate data, are widely used in many areas of sciences, engineering and social sciences as a tool to understand the information contained in typically large datasets. In the last two decades, random matrix theory-based tools such as the nearest neighbour eigenvalue spacing and eigenvector distributions have been employed to extract the significant modes of variability present in such empirical correlations. In this work, we present an alternative analysis in terms of the recently introduced spacing ratios, which does not require the cumbersome unfolding process. It is shown that the higher order spacing ratio distributions for the Wishart ensemble of random matrices, characterized by the Dyson index $\beta$, is related to the first order spacing ratio distribution with a modified value of co-dimension $\beta'$. This scaling is demonstrated for Wishart ensemble and also for the spectra of empirical correlation matrices drawn from the observed stock market and atmospheric pressure data. Using a combination of analytical and numerics, such scalings in spacing distributions are also discussed.


title: Deep Inelastic Scattering on an Extremal RN-AdS Black Hole

abstract: We consider deep inelastic scattering (DIS) on a large nucleus described as an extremal RN-AdS black hole using the holographic principle. Using the R-current correlators we determine the structure functions as a function Bjorken-x, and map it on a finite but large nucleus with fixed atomic number. The R-ratio of the nuclear structure functions exhibit strong shadowing at low-x.


title: Simulating the Magnetorotational Collapse of Supermassive Stars: Incorporating Gas Pressure Perturbations and Different Rotation Profiles

abstract: Collapsing supermassive stars (SMSs) with masses $M \gtrsim 10^{4-6}M_\odot$ have long been speculated to be the seeds that can grow and become supermassive black holes (SMBHs). We previously performed GRMHD simulations of marginally stable magnetized $\Gamma = 4/3$ polytropes uniformly rotating at the mass-shedding limit to model the direct collapse of SMSs. These configurations are supported entirely by thermal radiation pressure and model SMSs with $M \gtrsim 10^{6}M_\odot$. We found that around $90\%$ of the initial stellar mass forms a spinning black hole (BH) surrounded by a massive, hot, magnetized torus, which eventually launches an incipient jet. Here we perform GRMHD simulations of $\Gamma \gtrsim 4/3$, polytropes to account for the perturbative role of gas pressure in SMSs with $M \lesssim 10^{6}M_\odot$. We also consider different initial stellar rotation profiles. The stars are initially seeded with a dynamically weak dipole magnetic field that is either confined to the stellar interior or extended from its interior into the stellar exterior. We find that the mass of the BH remnant is $90\%-99\%$ of the initial stellar mass, depending sharply on $\Gamma -4/3$ as well as on the initial stellar rotation profile. After $t\sim 250-550M\approx 1-2\times 10^3(M/10^6M_\odot)$s following the BH formation, a jet is launched and it lasts for $\sim 10^4-10^5(M/10^6M_\odot)$s, consistent with the duration of long gamma-ray bursts. Our results suggest that the Blandford-Znajek mechanism powers the jet. They are also in agreement with our proposed universal model that estimates accretion rates and luminosities that characterize magnetized BH-disk remnant systems that launch a jet. This model helps explain why the outgoing luminosities for vastly different BH-disk formation scenarios all reside within a narrow range ($\sim 10^{52 \pm 1} \rm erg/s$), roughly independent of $M$.


title: The physics of multiphase gas flows: fragmentation of a radiatively cooling gas cloud in a hot wind

abstract: Galactic winds exhibit a multiphase structure that consists of hot-diffuse and cold-dense phases. Here we present high-resolution idealised simulations of the interaction of a hot supersonic wind with a cold cloud with the moving-mesh code arepo in setups with and without radiative cooling. We demonstrate that cooling causes clouds with sizes larger than the cooling length to fragment in two- and three-dimensional simulations (2D and 3D). We confirm earlier 2D simulations by McCourt et al. 2018 and highlight differences of the shattering processes of 3D clouds that are exposed to a hot wind. The fragmentation process is quantified with a friends-of-friends analysis of shattered cloudlets and density power spectra. Those show that radiative cooling causes the power spectral index to gradually increase when the initial cloud radius is larger than the cooling length and with increasing time until the cloud is fully dissolved in the hot wind. A resolution of around 1 pc is required to reveal the effect of cooling-induced fragmentation of a 100 pc outflowing cloud. Thus, state-of-the-art cosmological zoom simulations of the circumgalactic medium (CGM) fall short by orders of magnitudes from resolving this fragmentation process. This physics is, however, necessary to reliably model observed column densities and covering fractions of Lyman-$\alpha$ haloes, high-velocity clouds, and broad-line regions of active galactic nuclei.


title: Broadening Dark Matter Searches at the LHC: Mono-X versus Darkonium Channels

abstract: Current searches for dark matter at the LHC focus on mono-X signatures: the production of dark matter in association with a Standard Model (SM) particle. The simplest benchmark introduces a massive spin-1 mediator, the $Z^\prime$ boson, between the dark matter $\chi$ and the SM. Limits derived from mono-X channels are most effective when the mediator can decay into two on-shell dark matter particles: $M_{Z'}\gtrsim 2M_\chi$. We broaden the experimental reach into the complementary region, where the $Z^\prime$ mediator is much lighter than the dark matter. In this scenario the $Z^\prime$ mediates an effective long-range force between the dark matter, thereby facilitating the formation of darkonium bound states, as is common in many dark sector models. The darkonium becomes active when $M_{\chi}>M_{Z'}/\alpha_{\rm eff}$, where $\alpha_{\rm eff}$ is the effective fine-structure constant in the dark sector. Moreover, the darkonium could decay back into SM quarks, without producing missing transverse momentum in the detector. Considering multijet final states, we reinterpret existing searches to constrain the simple $Z^\prime$ benchmark beyond the region probed by mono-X searches. Assuming a baryonic $Z^\prime$ mediator and a Dirac dark matter, direct detection bounds can be loosened by giving a small Majorana mass to the dark matter. We also consider the interplay between mono-X and darkonium channels at future high energy colliders, which is at the frontier of probing the model parameter space.


title: Chromatic transit light curves of disintegrating rocky planets

abstract: Context. Kepler observations have revealed a class of short period exoplanets, of which Kepler-1520 b is the prototype, which have comet-like dust tails thought to be the result of small, rocky planets losing mass. The shape and chromaticity of the transits constrain the properties of the dust particles originating from the planet's surface, offering a unique opportunity to probe the composition and geophysics of rocky exoplanets. Aims. We aim to approximate the average Kepler long-cadence light curve of Kepler-1520 b and investigate how the optical thickness and transit cross-section of a general dust tail can affect the observed wavelength dependence and depth of transit light curves. Methods. We developed a new 3D model that ejects sublimating particles from the planet surface to build up a dust tail, assuming it to be optically thin, and used 3D radiative transfer computations that fully treat scattering using the distribution of hollow spheres (DHS) method, to generate transit light curves between 0.45 and 2.5 $\mu$m. Results. We show that the transit depth is wavelength independent for optically thick tails, potentially explaining why only some observations indicate a wavelength dependence. From the 3D nature of our simulated tails, we show that their transit cross-sections are related to the component of particle ejection velocity perpendicular to the planet's orbital plane and use this to derive a minimum ejection velocity of 1.2 kms$^{-1}$. To fit the average transit depth of Kepler-1520 b of 0.87%, we require a high dust mas-loss rate of 7 $-$ 80 M$_\oplus$ Gyr$^{-1}$ which implies planet lifetimes that may be inconsistent with the observed sample. Therefore, these mass-loss rates should be considered to be upper limits.


title: Heat-Bath Algorithmic Cooling with optimal thermalization strategies

abstract: Heat-Bath Algorithmic Cooling is a set of techniques for producing highly pure quantum systems by utilizing a surrounding heat-bath and unitary interactions. These techniques originally used the thermal environment only to fully thermalize ancillas at the environment temperature. Here we extend HBAC protocols by optimizing over the thermalization strategy. We find, for any $d$-dimensional system in an arbitrary initial state, provably optimal cooling protocols with surprisingly simple structure and exponential convergence to the ground state. Compared to the standard ones, these schemes can use fewer or no ancillas and exploit memory effects to enhance cooling. We verify that the optimal protocols are robusts to various deviations from the ideal scenario. For a single target qubit, the optimal protocol can be well approximated with a Jaynes-Cummings interaction between the system and a single thermal bosonic mode for a wide range of environmental temperatures. This admits an experimental implementation close to the setup of a micromaser, with a performance competitive with leading proposals in the literature. The proposed protocol provides an experimental setup that illustrates how non-Markovianity can be harnessed to improve cooling. On the technical side we 1. introduce a new class of states called maximally active states and discuss their thermodynamic significance in terms of optimal unitary control, 2. introduce a new set of thermodynamic processes, called $\beta$-permutations, whose access is sufficient to simulate a generic thermalization process, 3. show how to use abstract toolbox developed within the resource theory approach to thermodynamics to perform challenging optimizations, while combining it with open quantum system dynamics tools to approximate optimal solutions within physically realistic setups.


title: Direct randomized benchmarking for multi-qubit devices

abstract: Benchmarking methods that can be adapted to multi-qubit systems are essential for assessing the overall or "holistic" performance of nascent quantum processors. The current industry standard is Clifford randomized benchmarking (RB), which measures a single error rate that quantifies overall performance. But scaling Clifford RB to many qubits is surprisingly hard. It has only been performed on 1, 2, and 3 qubits as of this writing. This reflects a fundamental inefficiency in Clifford RB: the $n$-qubit Clifford gates at its core have to be compiled into large circuits over the 1- and 2-qubit gates native to a device. As $n$ grows, the quality of these Clifford gates quickly degrades, making Clifford RB impractical at relatively low $n$. In this Letter, we propose a direct RB protocol that mostly avoids compiling. Instead, it uses random circuits over the native gates in a device, seeded by an initial layer of Clifford-like randomization. We demonstrate this protocol experimentally on 2 -- 5 qubits, using the publicly available IBMQX5. We believe this to be the greatest number of qubits holistically benchmarked, and this was achieved on a freely available device without any special tuning up. Our protocol retains the simplicity and convenient properties of Clifford RB: it estimates an error rate from an exponential decay. But it can be extended to processors with more qubits -- we present simulations on 10+ qubits -- and it reports a more directly informative and flexible error rate than the one reported by Clifford RB. We show how to use this flexibility to measure separate error rates for distinct sets of gates, which includes tasks such as measuring an average CNOT error rate.


title: On the discontinuity of the $\pi_1$-action

abstract: We show the classical $\pi_1$-action on the $n$-th homotopy group can fail to be continuous for any $n$ when the homotopy groups are equipped with the natural quotient topology. In particular, we prove the action $\pi_1(X)\times\pi_n(X)\to\pi_n(X)$ fails to be continuous for a one-point union $X=A\vee \mathbb{H}_n$ where $A$ is an aspherical space such that $\pi_1(A)$ is a topological group and $\mathbb{H}_n$ is the $(n-1)$-connected, n-dimensional Hawaiian earring space $\mathbb{H}_n$ for which $\pi_n(\mathbb{H}_n)$ is a topological abelian group.


title: The GALAH survey: a catalogue of carbon-enhanced stars and CEMP candidates

abstract: Swan bands - characteristic molecular absorption features of the C$_2$ molecule - are a spectroscopic signature of carbon-enhanced stars. They can also be used to identify carbon-enhanced metal-poor (CEMP) stars. The GALAH (GALactic Archaeology with Hermes) is a magnitude-limited survey of stars producing high-resolution, high signal-to-noise spectra. We used 627,708 GALAH spectra to search for carbon-enhanced stars with a supervised and unsupervised classification algorithm, relying on the imprint of the Swan bands. We identified 918 carbon-enhanced stars, including 12 already described in the literature. An unbiased selection function of the GALAH survey allows us to perform a population study of carbon-enhanced stars. Most of them are giants, out of which we find 28 CEMP candidates. A large fraction of our carbon-enhanced stars with repeated observations show variation in radial velocity, hinting that there is a large fraction of variables among them. 32 of the detected stars also show strong Lithium enhancement in their spectra.


title: Prior Convictions: Black-Box Adversarial Attacks with Bandits and Priors

abstract: We study the problem of generating adversarial examples in a black-box setting in which only loss-oracle access to a model is available. We introduce a framework that conceptually unifies much of the existing work on black-box attacks, and we demonstrate that the current state-of-the-art methods are optimal in a natural sense. Despite this optimality, we show how to improve black-box attacks by bringing a new element into the problem: gradient priors. We give a bandit optimization-based algorithm that allows us to seamlessly integrate any such priors, and we explicitly identify and incorporate two examples. The resulting methods use two to four times fewer queries and fail two to five times less often than the current state-of-the-art.


title: Multi-criteria Evolution of Neural Network Topologies: Balancing Experience and Performance in Autonomous Systems

abstract: Majority of Artificial Neural Network (ANN) implementations in autonomous systems use a fixed/user-prescribed network topology, leading to sub-optimal performance and low portability. The existing neuro-evolution of augmenting topology or NEAT paradigm offers a powerful alternative by allowing the network topology and the connection weights to be simultaneously optimized through an evolutionary process. However, most NEAT implementations allow the consideration of only a single objective. There also persists the question of how to tractably introduce topological diversification that mitigates overfitting to training scenarios. To address these gaps, this paper develops a multi-objective neuro-evolution algorithm. While adopting the basic elements of NEAT, important modifications are made to the selection, speciation, and mutation processes. With the backdrop of small-robot path-planning applications, an experience-gain criterion is derived to encapsulate the amount of diverse local environment encountered by the system. This criterion facilitates the evolution of genes that support exploration, thereby seeking to generalize from a smaller set of mission scenarios than possible with performance maximization alone. The effectiveness of the single-objective (optimizing performance) and the multi-objective (optimizing performance and experience-gain) neuro-evolution approaches are evaluated on two different small-robot cases, with ANNs obtained by the multi-objective optimization observed to provide superior performance in unseen scenarios.


title: Anisotropic Purity of Entangled Photons From Cooper Pairs in Heterostructures

abstract: It was theoretically proposed that a forward-biased p-n junction with a superconducting layer (P-N-S), would produce pure, polarization entangled photons due to inherent spin-singlet pairing of electrons. However, any heterostructure interface generically induces Rashba spin-orbit coupling, which in turn generates a mixed singlet-triplet superconducting order parameter. Here we study the effect of triplet pairing on the purity of photons produced through Cooper pair recombination. A unique directional dependence of the state purity is found for a triplet superconductor with fixed $\vec{d}$: pure, entangled photons are produced when the photon polarization axis is parallel to $\vec{d}$. Induced triplet pairing in a singlet superconductor is shown to degrade the state purity, while induced singlet pairing in a triplet superconductor is shown to enhance the production of entangled pairs. These considerations may aid the design of functional devices to produce entangled photons.


title: Non-Riemannian Description of Robinson-Trautman Spacetimes in Brans-Dicke Theory Of Gravity

abstract: The variational field equations of Brans-Dicke scalar-tensor theory of gravitation are given in a non-Riemannian setting in the language of exterior differential forms over 4-dimensional spacetimes. A conformally re-scaled Robinson-Trautman metric together with the Brans-Dicke scalar field are used to characterise algebraically special Robinson-Trautman spacetimes. All the relevant tensors are worked out in a complex null basis and given explicitly in an appendix for future reference. Some special families of solutions are also given and discussed.


title: Visitors to urban greenspace have higher sentiment and lower negativity on Twitter

abstract: With more people living in cities, we are witnessing a decline in exposure to nature. A growing body of research has demonstrated an association between nature contact and improved mood. Here, we used Twitter and the Hedonometer, a world analysis tool, to investigate how sentiment, or the estimated happiness of the words people write, varied before, during, and after visits to San Francisco's urban park system. We found that sentiment was substantially higher during park visits and remained elevated for several hours following the visit. Leveraging differences in vegetative cover across park types, we explored how different types of outdoor public spaces may contribute to subjective well-being. Tweets during visits to Regional Parks, which are greener and have greater vegetative cover, exhibited larger increases in sentiment than tweets during visits to Civic Plazas and Squares. Finally, we analyzed word frequencies to explore several mechanisms theorized to link nature exposure with mental and cognitive benefits. Negation words such as 'no', 'not', and 'don't' decreased in frequency during visits to urban parks. These results can be used by urban planners and public health officials to better target nature contact recommendations for growing urban populations.


title: Hydrodynamics with spacetime-dependent scattering length

abstract: Hydrodynamics provides a concise but powerful description of long-time and long-distance physics of correlated systems out of thermodynamic equilibrium. Here we construct hydrodynamic equations for nonrelativistic particles with a spacetime-dependent scattering length and show that it enters constitutive relations uniquely so as to represent the fluid expansion and contraction in both normal and superfluid phases. As a consequence, we find that a leading dissipative correction to the contact density due to the spacetime-dependent scattering length is proportional to the bulk viscosity ($\zeta_2$ in the superfluid phase). Also, when the scattering length is slowly varied over time in a uniform system, the entropy density is found to be produced even without fluid flows in proportion to the bulk viscosity, which may be useful as a novel probe to measure the bulk viscosity in ultracold-atom experiments.


title: Attention Models in Graphs: A Survey

abstract: Graph-structured data arise naturally in many different application domains. By representing data as graphs, we can capture entities (i.e., nodes) as well as their relationships (i.e., edges) with each other. Many useful insights can be derived from graph-structured data as demonstrated by an ever-growing body of work focused on graph mining. However, in the real-world, graphs can be both large - with many complex patterns - and noisy which can pose a problem for effective graph mining. An effective way to deal with this issue is to incorporate "attention" into graph mining solutions. An attention mechanism allows a method to focus on task-relevant parts of the graph, helping it to make better decisions. In this work, we conduct a comprehensive and focused survey of the literature on the emerging field of graph attention models. We introduce three intuitive taxonomies to group existing work. These are based on problem setting (type of input and output), the type of attention mechanism used, and the task (e.g., graph classification, link prediction, etc.). We motivate our taxonomies through detailed examples and use each to survey competing approaches from a unique standpoint. Finally, we highlight several challenges in the area and discuss promising directions for future work.


title: A Luminous Transient Event in a Sample of WISE-Selected Variable AGN

abstract: Recently Assef et al.(2018) presented two catalogs of AGN candidates over 30,093 deg^2 selected from the Wide-field Infrared Survey Explorer (WISE) observations. From their most reliable sample, Assef et al. (2018) identified 45 AGN candidates with the highest variability levels in the AllWISE catalog, but that are not blazars. Here we present new spectroscopic observations of some of these targets to further constrain their nature. We also study their optical lightcurves using observations from CRTS, and find that only seven show significant optical variability, and that five of those seven are spectroscopically classified as AGN. In one of them, WISEA J094806.56+031801.7 (W0948+0318), we identify a transient event in the CRTS lightcurve. We present a detailed analysis of this transient, and characterize it through its CRTS lightcurve and its multi-wavelength spectral energy distribution obtained from GALEX, Pan-STARRS and WISE observations. We find that the most likely source of the transient is a super-luminous supernova (SLSN) in W0948+0318. We estimate the total radiated energy to be E=1.6\pm 0.3 x 10^52 erg, making it one of the most energetic SLSN observed. Based on the lack of change in mid-IR color throughout and after the transient event, we speculate that the location of the SLSN is within the torus of the AGN. We identify 9 possible analogs to W0948+0318 based on their WISE lightcurves. None show optically detected transients and hence suggest significant dust obscuration. Finally, we estimate a rate of >2x10^-7 yr^-1 per AGN for these transients under the conservative assumption that none of the identified analogs have a common origin with the transient in W0948+0318.


title: Spatial Imaging of Charge Transport in Silicon at Low Temperature

abstract: We present direct imaging measurements of charge transport across a 1 cm$\times$ 1 cm$\times$ 4 mm crystal of high purity silicon ($\sim$20 k$\Omega$cm) at temperatures between 500 mK and and 5 K. We use these data to determine the intervalley scattering rate of electrons as a function of the electric field applied along the $\langle 111 \rangle$ crystal axis, and we present a phenomenological model of intervalley scattering that explains the constant scattering rate seen at low-voltage for cryogenic temperatures. We also demonstrate direct imaging measurements of effective hole mass anisotropy, which is strongly dependent on both temperature and electric field strength. The observed effects can be explained by a warping of the valence bands for carrier energies near the spin-orbit splitting energy in silicon.


title: Deep Learning

abstract: Deep learning (DL) is a high dimensional data reduction technique for constructing high-dimensional predictors in input-output models. DL is a form of machine learning that uses hierarchical layers of latent features. In this article, we review the state-of-the-art of deep learning from a modeling and algorithmic perspective. We provide a list of successful areas of applications in Artificial Intelligence (AI), Image Processing, Robotics and Automation. Deep learning is predictive in its nature rather then inferential and can be viewed as a black-box methodology for high-dimensional function estimation.


title: Radiative neutrino mass via fermion kinetic mixing

abstract: We propose that the radiative generation of the neutrino mass can be achieved by incorporating the kinetic mixing of fermion fields which arises radiatively at one-loop level. As a demonstrative example of the application of the mechanism, we present the particular case of the Standard Model extension by $U(1)_D$ symmetry. As a result, we show how neutrino masses can be generated via a kinetic mixing portal instead of a mass matrix with residual symmetries responsible for the stability of multicomponent dark matter.


title: Einstein gravity from the ${\N=4}$ spinning particle

abstract: We obtain a manifestly background independent BRST quantization of the $\N=4$ supersymmetric spinning particle. We show that nilpotency of the BRST charge $Q$ implies the Einstein equations admitting a cosmological constant of indefinite sign. The physical graviton states are given by the vertex operator, obtained by an infinitesimal variation of $Q$, acting on diffeomorphism ghost states. In addition, the tree-point graviton scattering vertex is correctly reproduced by the worldline computation with these vertex operators.


title: Intermediate scattering function and quantum recoil in non-Markovian quantum diffusion

abstract: Exact expressions are derived for the intermediate scattering function (ISF) of a quantum particle diffusing in a harmonic potential and linearly coupled to a harmonic bath. The results are valid for arbitrary strength and spectral density of the coupling. The general, exact non-Markovian result is expressed in terms of the classical velocity autocorrelation function, which represents an accumulated phase during a scattering event. The imaginary part of the exponent of the ISF is proportional to the accumulated phase, which is an antisymmetric function of the correlation time $t$. The expressions extend previous results given in the quantum Langevin framework where the classical response of the bath was taken as Markovian. For a special case of non-Markovian friction, where the friction kernel decays exponentially in time rather than instantaneously, we provide exact results relating to unconfined quantum diffusion, and identify general features that allow insight to be exported to more complex examples. The accumulated phase as a function of the t has a universal gradient at the origin, depending only on the mass of the diffusing system particle. At large t the accumulated phase reaches a constant limit that depends only on the classical diffusion coefficient and is therefore independent of the detailed memory properties of the friction kernel. Non-Markovian properties of the friction kernel are encoded in the details of how the accumulated phase switches from its $t\rightarrow -\infty$ to its $t\rightarrow -\infty$ limit, subject to the constraint of the universal gradient. When memory effects are significant, the transition from one limit to the other becomes non-monotonic, owing to oscillations in the classical velocity autocorrelation. The result is interpreted in terms of a solvent caging effect, in which slowly fluctuating bath modes create transient wells for the system particle.


title: Knowledge Integration for Disease Characterization: A Breast Cancer Example

abstract: With the rapid advancements in cancer research, the information that is useful for characterizing disease, staging tumors, and creating treatment and survivorship plans has been changing at a pace that creates challenges when physicians try to remain current. One example involves increasing usage of biomarkers when characterizing the pathologic prognostic stage of a breast tumor. We present our semantic technology approach to support cancer characterization and demonstrate it in our end-to-end prototype system that collects the newest breast cancer staging criteria from authoritative oncology manuals to construct an ontology for breast cancer. Using a tool we developed that utilizes this ontology, physician-facing applications can be used to quickly stage a new patient to support identifying risks, treatment options, and monitoring plans based on authoritative and best practice guidelines. Physicians can also re-stage existing patients or patient populations, allowing them to find patients whose stage has changed in a given patient cohort. As new guidelines emerge, using our proposed mechanism, which is grounded by semantic technologies for ingesting new data from staging manuals, we have created an enriched cancer staging ontology that integrates relevant data from several sources with very little human intervention.


title: On graphs with 2 trivial distance ideals

abstract: Distance ideals generalize the Smith normal form of the distance matrix of a graph. The family of graphs with 2 trivial distance ideals contains the family of graphs whose distance matrix has at most 2 invariant factors equal to 1. Here we give an infinite family of forbidden induced subgraphs for the graphs with 2 trivial distance ideals. These are also related with other well known graph classes.


title: Is the Cometary Nucleus Extraction Technique Reliable?

abstract: It depends. Our experiment reveals that, given an optically thin coma, generally, the smaller the signal ratio of nucleus to coma, the less reliable is the cometary nucleus-extraction technique. We strongly suggest the technique only be applied to cases where the nucleus signal occupies $\gtrsim$10% of the total signal wherein the bias is no more than a few percent. Otherwise there is probably no way to debias results from this technique in reality, since its reliability is highly sensitive to entangling complications, including the coma profile, and the point-spread function (PSF).


title: A Stochastic Line Search Method with Convergence Rate Analysis

abstract: For deterministic optimization, line-search methods augment algorithms by providing stability and improved efficiency. We adapt a classical backtracking Armijo line-search to the stochastic optimization setting. While traditional line-search relies on exact computations of the gradient and values of the objective function, our method assumes that these values are available up to some dynamically adjusted accuracy which holds with some sufficiently large, but fixed, probability. We show the expected number of iterations to reach a near stationary point matches the worst-case efficiency of typical first-order methods, while for convex and strongly convex objective, it achieves rates of deterministic gradient descent in function values.


title: On exceptional sets in Manin's Conjecture

abstract: In this survey paper we study Manin's Conjecture from a geometric perspective. The focus of the paper is the recent conjectural description of the exceptional set in Manin's Conjecture due to Lehmann-Sengupta-Tanimoto. After giving an extensive background, we give a precise description of this set and compute it in many examples.


title: Variance propagation for density surface models

abstract: Spatially-explicit estimates of population density, together with appropriate estimates of uncertainty, are required in many management contexts. Density Surface Models (DSMs) are a two-stage approach for estimating spatially-varying density from distance-sampling data. First, detection probabilities -- perhaps depending on covariates -- are estimated based on details of individual encounters; next, local densities are estimated using a GAM, by fitting local encounter rates to location and/or spatially-varying covariates while allowing for the estimated detectabilities. One criticism of DSMs has been that uncertainty from the two stages is not usually propagated correctly into the final variance estimates. We show how to reformulate a DSM so that the uncertainty in detection probability from the distance sampling stage (regardless of its complexity) is captured as an extra random effect in the GAM stage. In effect, we refit an approximation to the detection function model at the same time as fitting the spatial model. This allows straightforward computation of the overall variance via exactly the same software already needed to fit the GAM. A further extension allows for spatial variation in group size, which can be an important covariate for detectability as well as directly affecting abundance. We illustrate these models using point transect survey data of Island Scrub-Jays on Santa Cruz Island, CA and harbour porpoise from the SCANS-II line transect survey of European waters.


title: Lifetime of Majorana qubits in Rashba nanowires with non-uniform chemical potential

abstract: We study the lifetime of topological qubits based on Majorana bound states hosted in a one-dimensional Rashba nanowire (NW) with proximity-induced superconductivity and non-uniform chemical potential needed for manipulation and read-out. If nearby gates tune the chemical potential locally so that part of the NW is in the trivial phase, Andreev bound states (ABSs) can emerge which are localized at the interface between topological and trivial phases with energies significantly less than the gap. The emergence of such subgap states strongly decreases the Majorana qubit lifetime at finite temperatures due to local perturbations that can excite the system into these ABSs. Using Keldysh formalism, we study such excitations caused by fluctuating charges in capacitively coupled gates and calculate the corresponding Majorana lifetimes due to thermal noise, which are shown to be much shorter than those in NWs with uniform chemical potential.


title: Convolutional Neural Networks Analyzed via Inverse Problem Theory and Sparse Representations

abstract: Inverse problems in imaging such as denoising, deblurring, superresolution (SR) have been addressed for many decades. In recent years, convolutional neural networks (CNNs) have been widely used for many inverse problem areas. Although their indisputable success, CNNs are not mathematically validated as to how and what they learn. In this paper, we prove that during training, CNN elements solve for inverse problems which are optimum solutions stored as CNN neuron filters. We discuss the necessity of mutual coherence between CNN layer elements in order for a network to converge to the optimum solution. We prove that required mutual coherence can be provided by the usage of residual learning and skip connections. We have set rules over training sets and depth of networks for better convergence, i.e. performance.


title: The intermediate scattering function for quasi-elastic scattering in the presence of memory friction

abstract: We derive an analytical expression for the intermediate scattering function of a particle on a flat surface obeying the Generalised Langevin Equation, with exponential memory friction. Numerical simulations based on an extended phase space method confirm the analytical results. The simulated trajectories provide qualitative insight into the effect that introducing a finite memory timescale has on the analytical line shapes. The relative amplitude of the long-time exponential tail of the line shape is suppressed, but its decay rate is unchanged, reflecting the fact that the cutoff frequency of the exponential kernel affects short-time correlations but not the diffusion coefficient which is defined in terms of a long-time limit. The exponential sensitivity of the relative amplitudes to the decay time of the chosen memory kernel is a very strong indicator for the prospect of inferring a friction kernel and the physical insights from experimentally measured intermediate scattering functions.


title: Abstractive and Extractive Text Summarization using Document Context Vector and Recurrent Neural Networks

abstract: Sequence to sequence (Seq2Seq) learning has recently been used for abstractive and extractive summarization. In current study, Seq2Seq models have been used for eBay product description summarization. We propose a novel Document-Context based Seq2Seq models using RNNs for abstractive and extractive summarizations. Intuitively, this is similar to humans reading the title, abstract or any other contextual information before reading the document. This gives humans a high-level idea of what the document is about. We use this idea and propose that Seq2Seq models should be started with contextual information at the first time-step of the input to obtain better summaries. In this manner, the output summaries are more document centric, than being generic, overcoming one of the major hurdles of using generative models. We generate document-context from user-behavior and seller provided information. We train and evaluate our models on human-extracted-golden-summaries. The document-contextual Seq2Seq models outperform standard Seq2Seq models. Moreover, generating human extracted summaries is prohibitively expensive to scale, we therefore propose a semi-supervised technique for extracting approximate summaries and using it for training Seq2Seq models at scale. Semi-supervised models are evaluated against human extracted summaries and are found to be of similar efficacy. We provide side by side comparison for abstractive and extractive summarizers (contextual and non-contextual) on same evaluation dataset. Overall, we provide methodologies to use and evaluate the proposed techniques for large document summarization. Furthermore, we found these techniques to be highly effective, which is not the case with existing techniques.


title: $\hat p \cdot \hat A$ vs $\hat x \cdot \hat E$: Gauge invariance in quantum optics

abstract: We compare the predictions of the fundamentally motivated minimal coupling ($\hat{\boldsymbol{p}}\cdot\hat{\boldsymbol{A}}$) and the ubiquitous dipole coupling ($\hat{\boldsymbol{x}}\cdot\hat{\boldsymbol{E}}$) in the light-matter interaction. By studying the light-matter interaction for hydrogen-like atoms we find that the dipole approximation cannot be a-priori justified to analyze the physics of vacuum excitations (a very important phenomenon in relativistic quantum information) since a dominant wavelength is absent in those problems, no matter how small (as compared to any frequency scale) the atom is. Remarkably, we show that the dipole approximation in those regimes can still be valid as long as the interaction time is longer than the light-crossing time of the atoms, which is a very reasonable assumption. We also highlight some of the subtleties that one has to be careful with when working with the explicitly non-gauge invariant nature of the minimal coupling and we compare it with the explicitly gauge invariant dipole coupling.


title: Regularity of the singular set in a two-phase problem for harmonic measure with H\"older data

abstract: In non-variational two-phase free boundary problems for harmonic measure, we examine how the relationship between the interior and exterior harmonic measures of a domain $\Omega \subset \mathbb{R}^n$ influences the geometry of its boundary. This type of free boundary problem was initially studied by Kenig and Toro in 2006 and was further examined in a series of separate and joint investigations by several authors. The focus of the present paper is on the singular set in the free boundary, where the boundary looks infinitesimally like zero sets of homogeneous harmonic polynomials of degree at least 2. We prove that if the Radon-Nikodym derivative of the exterior harmonic measure with respect to the interior harmonic measure has a H\"older continuous logarithm, then the free boundary admits unique geometric blowups at every singular point and the singular set can be covered by countably many $C^{1,\beta}$ submanifolds of dimension at most $n-3$. This result is partly obtained by adapting tools such as Garofalo and Petrosyan's Weiss type monotonicity formula and an epiperimetric inequality for harmonic functions from the variational to the non-variational setting.


title: ScaRR: Scalable Runtime Remote Attestation for Complex Systems

abstract: The introduction of remote attestation (RA) schemes has allowed academia and industry to enhance the security of their systems. The commercial products currently available enable only the validation of static properties, such as applications fingerprint, and do not handle runtime properties, such as control-flow correctness. This limitation pushed researchers towards the identification of new approaches, called runtime RA. However, those mainly work on embedded devices, which share very few common features with complex systems, such as virtual machines in a cloud. A naive deployment of runtime RA schemes for embedded devices on complex systems faces scalability problems, such as the representation of complex control-flows or slow verification phase. In this work, we present ScaRR: the first Scalable Runtime Remote attestation schema for complex systems. Thanks to its novel control-flow model, ScaRR enables the deployment of runtime RA on any application regardless of its complexity, by also achieving good performance. We implemented ScaRR and tested it on the benchmark suite SPEC CPU 2017. We show that ScaRR can validate on average 2M control-flow events per second, definitely outperforming existing solutions.


title: Robust Resilient Signal Reconstruction under Adversarial Attacks

abstract: We consider the problem of signal reconstruction for a system under sparse signal corruption by a malicious agent. The reconstruction problem follows the standard error coding problem that has been studied extensively in the literature. We include a new challenge of robust estimation of the attack support. The problem is then cast as a constrained optimization problem merging promising techniques in the area of deep learning and estimation theory. A pruning algorithm is developed to reduce the ``false positive" uncertainty of data-driven attack localization results, thereby improving the probability of correct signal reconstruction. Sufficient conditions for the correct reconstruction and the associated reconstruction error bounds are obtained for both exact and inexact attack support estimation. Moreover, a simulation of a water distribution system is presented to validate the proposed techniques.


title: Interior X-ray diffraction tomography with low-resolution exterior information

abstract: X-ray diffraction tomography (XDT) resolves spatially-variant XRD profiles within macroscopic objects, and provides improved material contrast compared to the conventional transmission-based computed tomography (CT). However, due to the small diffraction cross-section, XDT suffers from long imaging acquisition time, which could take tens of hours for a full scan using a table-top X-ray tube. In medical and industrial imaging applications, oftentimes only the XRD measurement within a region-of-interest (ROI) is required, which, together with the demand to reduce imaging time and radiation dose to the sample, motivates the development of interior XDT systems that scan and reconstruct only an internal region within the sample. The interior problem does not have a unique solution, and a direct inversion on the truncated projection data often leads to large reconstruction errors in ROI. To reduce the truncation artifacts, conventional attenuation-based interior reconstruction problems rely on a known region or piecewise constant constraint within the ROI, which do not apply to all the samples in XDT. Here we propose a quasi-interior XDT scheme that incorporates a small fraction of projection information from the exterior region to assist interior reconstruction. The low-resolution exterior projection data obviates the requirement for prior knowledge on the object, and allows the ROI reconstruction to be performed with the fast, widely-used filtered back-projection algorithm for easy integration into real-time XDT imaging modules.


title: Magnetotransport in Sr3PbO antiperovskite with three-dimensional massive Dirac electrons

abstract: Novel topological phenomena are anticipated for three-dimensional (3D) Dirac electrons. The magnetotransport properties of cubic ${\rm Sr_{3}PbO}$ antiperovskite, theoretically proposed to be a 3D massive Dirac electron system, are studied. The measurements of Shubnikov-de Haas oscillations and Hall resistivity indicate the presence of a low density ($\sim 1 \times 10^{18}$ ${\rm cm^{-3}}$) of holes with an extremely small cyclotron mass of 0.01-0.06$m_{e}$. The magnetoresistance $\Delta\rho_{xx}(B)$ is linear in magnetic field $B$ with the magnitude independent of temperature. These results are fully consistent with the presence of 3D massive Dirac electrons in ${\rm Sr_{3}PbO}$. The chemical flexibility of the antiperovskites and our findings in the family member, ${\rm Sr_{3}PbO}$, point to their potential as a model system in which to explore exotic topological phases.


title: Cosmological Perturbations in the Early Universe

abstract: One essence in capturing the history of primordial fluctuations that arise during inflation and eventually lead to formation of large scale structures in the universe, relies on quantizing general relativity coupled to a scalar field. This is a system that respects diffeomorphism invariance, the freedom of choosing the coordinate system to work with without changing the physics. Hence the way to handle its quantization requires an understanding of how to quantize diffeomorphisms.Deciding on suitable coordinate choices and making sure that this gauge fixing, which is unavoidable in any calculation, does not effect the end result, is tricky. In this thesis we make full use of the effects of diffeomorphism invariance of the theory on the primordial fluctuations to pursue two different approaches that focus on treating perturbations after gauge fixing. On one hand we work towards developing our understanding of how to handle quantization in terms of Dirac quantization and Becchi, Rouet, Stora, Tyutin (BRST) quantization. On another, we focus on how to generalize the allowed interactions and understand the scales they bring in the era of preheating that follows inflation, with effective field theory (EFT) methods on cosmological backgrounds.


title: Ensemble of Deep Learned Features for Melanoma Classification

abstract: The aim of this work is to propose an ensemble of descriptors for Melanoma Classification, whose performance has been evaluated on validation and test datasets of the melanoma challenge 2018. The system proposed here achieves a strong discriminative power thanks to the combination of multiple descriptors. The proposed system represents a very simple yet effective way of boosting the performance of trained CNNs by composing multiple CNNs into an ensemble and combining scores by sum rule. Several types of ensembles are considered, with different CNN architectures along with different learning parameter sets. Moreover CNN are used as feature extractors: an input image is processed by a trained CNN and the response of a particular layer (usually the classification layer, but also internal layers can be employed) is treated as a descriptor for the image and used for training a set of Support Vector Machines (SVM).


title: On the Cantor-Bendixson rank of the Grigorchuk group and the Gupta-Sidki $3$ group

abstract: We study the Cantor--Bendixson rank of the space of subgroups for members of a general class of finitely generated self-replicating branch groups. In particular, we show for $G$ either the Grigorchuk group or the Gupta--Sidki $3$ group, the Cantor--Bendixson rank of $\mathrm{Sub}(G)$ is $\omega$. For each natural number $n$, we additionally characterize the subgroups of rank $n$ and give a description of subgroups in the perfect kernel.


title: Measurement of the Iron Spectrum in Cosmic Rays by VERITAS

abstract: We present a new measurement of the energy spectrum of iron nuclei in cosmic rays from 20 to 500 TeV. The measurement makes use of a template-based analysis method, which, for the first time, is applied to the energy reconstruction of iron-induced air showers recorded by the VERITAS array of imaging atmospheric Cherenkov telescopes. The event selection makes use of the direct Cherenkov light which is emitted by charged particles before the first interaction, as well as other parameters related to the shape of the recorded air shower images. The measured spectrum is well described by a power law $\frac{\mathrm{d} F}{\mathrm{d} E}=f_0\cdot \left(\frac{E}{E_0}\right)^{-\gamma}$ over the full energy range, with $\gamma = 2.82 \pm 0.30 \mathrm{(stat.)} ^{+0.24}_{-0.27} \mathrm{(syst.)}$ and $f_0 = \left( 4.82 \pm 0.98 \mathrm{(stat.)}^{+2.12}_{-2.70} \mathrm{(syst.)} \right)\cdot 10^{-7}$m$^{-2}$s$^{-1}$sr$^{-1}$TeV$^{-1}$ at $E_0=50$TeV, with no indication of a cutoff or spectral break. The measured differential flux is compatible with previous results, with improved statistical uncertainty at the highest energies.


title: Shared Processor Scheduling of Multiprocessor Jobs

abstract: We study shared processor scheduling of $\textit{multiprocessor}$ weighted jobs where each job can be executed on its private processor and simultaneously on possibly $\textit{many}$ processors shared by all jobs in order to reduce their completion times due to processing time overlap. Each of $m$ shared processors may charge different fee but otherwise the processors are identical. The total weighted overlap of all jobs is to be maximized. This problem is key to subcontractor scheduling in extended enterprises and supply chains, and divisible load scheduling in computing. We prove that, quite surprisingly, $\textit{synchronized}$ schedules that complete each job using shared processors at the same time on its private and shared processors include optimal schedules. We show that optimal $\alpha$-$\textit{private}$ schedules that require each job to use its private processor for at least $\alpha=1/2+1/(4(m+1))$ of the time required by the job guarantee more than an $\alpha$ fraction of the total weighted overlap of the optimal schedules. This gives an $\alpha$-approximation algorithm that runs in strongly polynomial time for the problem, and improves the $1/2$-approximation reported recently in the literature to $5/8$-approximation for a single shared processor problem. The computational complexity of the problem, both single and multi-shared processor, remains open. We show however an LP-based optimal algorithm for $\textit{antithetical}$ instances where for any pair of jobs $j$ and $i$, if the processing time of $j$ is smaller than or equal to the processing time of $i$, then the weight of $j$ is greater than or equal to the weight of $i$.


title: Anharmonic stabilization and lattice heat transport in rocksalt $\beta$-GeTe

abstract: Peierls-Boltzmann transport equation, coupled with third-order anharmonic lattice dynamics calculations, has been widely used to model lattice thermal conductivity ($\kappa_{l}$) in bulk crystals. However, its application to materials with structural phase transition at relatively high temperature is fundamentally challenged by the presence of lattice instabilities (imaginary phonon modes). Additionally, its accuracy suffers from the absence of higher-than-third-order phonon scattering processes, which are important near/above the Debye temperature. In this letter, we present an effective scheme that combines temperature-induced anharmonic phonon renormalization and four-phonon scattering to resolve these two theoretical challenges. We apply this scheme to investigate the lattice dynamics and thermal transport properties of GeTe, which undergoes a second-order ferroelectric phase transition from rhombohedral $\alpha$-GeTe to rocksalt $\beta$-GeTe at about 700~K. Our results on the high-temperature phase $\beta$-GeTe at 800~K confirm the stabilization of $\beta$-GeTe by temperature effects. We find that considering only three-phonon scattering leads to significantly overestimated $\kappa_{l}$ of 3.8~W/mK at 800~K, whereas including four-phonon scattering reduces $\kappa_{l}$ to 1.7~W/mK, a value comparable with experiments. To explore the possibility to further suppress $\kappa_{l}$, we show that alloying $\beta$-GeTe with heavy cations such as Pb and Bi can effectively reduce $\kappa_{l}$ to about 1.0~W/mK, whereas sample size needs to be around 10nm through nanostructuring to achieve a comparable reduction of $\kappa_{l}$.


title: Coarse-Graining Simulation Approaches for Polymer Melts: Range of Potential and Computational Efficiency

abstract: The integral equation coarse-graining (IECG) approach is a promising high-level coarse-graining (CG) method for polymer melts, with variable resolution from soft spheres to multi CG sites, which preserves the structural and thermodynamical consistencies with the related atomistic simulations. When compared to the atomistic description, the procedure of coarse-graining results in smoother free energy surfaces, longer-ranged potentials, a decrease in the number of interaction sites for a given polymer, and more. Because these changes have competing effects on the computational efficiency of the CG model, care needs to be taken when studying the effect of coarse-graining on the computational speed-up in CG molecular dynamics simulations. For instance, treatment of long-range CG interactions requires the selection of cutoff distances that include the attractive part of the effective CG potential and force. In particular, we show how the complex nature of the range and curvature of the effective CG potential, the selection of a suitable CG timestep, the choice of the cutoff distance, the molecular dynamics algorithms, and the smoothness of the CG free energy surface affect the efficiency of IECG simulations. By direct comparison with the atomistic simulations of relatively short chain polymer melts, we find that the overall computational efficiency is highest for the highest level of CG (soft spheres), with an overall improvement of the computational efficiency being about $10^6-10^8$ for various CG levels/resolutions. Therefore, the IECG method can have important applications in molecular dynamics simulations of polymeric systems. Finally, making use of the standard spatial decomposition algorithm, the parallel scalability of the IECG simulations for various levels of CG is presented. Optimal parallel scaling is observed for a reasonably large number of processors.


title: Maunakea Spectroscopic Explorer Advancing from Conceptual Design

abstract: The Maunakea Spectroscopic Explorer (MSE) project has completed its Conceptual Design Phase. This paper is a status report of the MSE project regarding its technical and programmatic progress. The technical status includes its conceptual design and system performance, and highlights findings and recommendations from the System and various subsystems design reviews. The programmatic status includes the project organization and management plan for the Preliminary Design Phase. In addition, this paper provides the latest information related to the permitting process for Maunakea construction.


title: Uma an\'alise comparativa de ferramentas de an\'alise est\'atica para dete\c{c}\~ao de erros de mem\'oria

abstract: --- Portuguese version As falhas de software est\~ao com frequ\^encia associadas a acidentes com graves consequ\^encias econ\'omicas e/ou humanas, pelo que se torna imperioso investir na valida\c{c}\~ao do software, nomeadamente daquele que \'e cr\'itico. Este artigo endere\c{c}a a tem\'atica da qualidade do software atrav\'es de uma an\'alise comparativa da usabilidade e efic\'acia de quatro ferramentas de an\'alise est\'atica de programas em C/C++. Este estudo permitiu compreender o grande potencial e o elevado impacto que as ferramentas de an\'alise est\'atica podem ter na valida\c{c}\~ao e verifica\c{c}\~ao de software. Como resultado complementar, foram identificados novos erros em programas de c\'odigo aberto e com elevada popularidade, que foram reportados. --- English version Software bugs are frequently associated with accidents with serious economical and/or human consequences, being thus imperative the investment in the validation of software, namely of the critical one. This article addresses the topic of software quality by making a comparative analysis of the usability and efficiency of four static analysis tools for C/C++ programs. This study allow to understand the big potential and high impact that these tools may have in the validation and verification of software. As a complementary result, we identified new errors in very popular open source projects, which have been reported.


title: Impact of high-order tidal terms on binary neutron-star waveforms

abstract: GW170817, the milestone gravitational-wave event originated from a binary neutron star merger, has allowed scientific community to place a constraint on the equation of state of neutron stars by extracting the leading-order, tidal-deformability term from the gravitational waveform. Here we incorporate tidal corrections to the gravitational-wave phase at next-to-leading and next-to-next-to-leading order, including the magnetic tidal Love numbers, tail effects, and the spin-tidal couplings recently computed in Tiziano Abdelsalhin [Phys. Rev. D 98, 104046 (2018)]. These effects have not yet been included in the waveform approximants for the analysis of GW170817. We provide a qualitative and quantitative analysis of the impact of these new terms by studying the parameter bias induced on events compatible with GW170817 assuming second-generation (advanced LIGO) and third-generation (Einstein Telescope) ground-based gravitational-wave interferometers. We find that including the tidal-tail term deteriorates the convergence properties of the post-Newtonian expansion in the relevant frequency range. We also find that the effect of magnetic tidal Love numbers could be measurable for an optimal GW170817 event with signal-to-noise ratio $\rho \approx 1750$ detected with the Einstein Telescope. On the same line, spin-tidal couplings may be relevant if mildly high-spin $\chi \gtrsim 0.1$ neutron star binaries exist in nature.


title: Maunakea Spectroscopic Explorer (MSE) Instrumentation Suite

abstract: The Maunakea Spectroscopic Explorer (MSE) is replacement of the existing 3.6-m Canada France Hawaii Telescope into a dedicated wide field highly multiplexed fiber fed spectroscopic facility. MSE is capable of observing over four thousand science targets simultaneously in two resolution modes. The paper describes the unique instrument system capabilities and its components starting from the telescope prime focus and ending at the spectrograph suite. The instrument system components have completed their conceptual designs and they include a Sphinx positioner system, fiber transmission system, low/moderate resolution and high resolution spectrographs and a calibration system. These components will be procured separately and the Project is responsible for their integration and the overall system performance afterward. The paper describes from a system perspective the specific design and interface constraints imposed on the components given the extra interface and integration considerations.


title: Information estimation using nonparametric copulas

abstract: Estimation of mutual information between random variables has become crucial in a range of fields, from physics to neuroscience to finance. Estimating information accurately over a wide range of conditions relies on the development of flexible methods to describe statistical dependencies among variables, without imposing potentially invalid assumptions on the data. Such methods are needed in cases that lack prior knowledge of their statistical properties and that have limited sample numbers. Here we propose a powerful and generally applicable information estimator based on non-parametric copulas. This estimator, called the non-parametric copula-based estimator (NPC), is tailored to take into account detailed stochastic relationships in the data independently of the data's marginal distributions. The NPC estimator can be used both for continuous and discrete numerical variables and thus provides a single framework for the mutual information estimation of both continuous and discrete data. By extensive validation on artificial samples drawn from various statistical distributions, we found that the NPC estimator compares well against commonly used alternatives. Unlike methods not based on copulas, it allows an estimation of information that is robust to changes of the details of the marginal distributions. Unlike parametric copula methods, it remains accurate regardless of the precise form of the interactions between the variables. In addition, the NPC estimator had accurate information estimates even at low sample numbers, in comparison to alternative estimators. The NPC estimator therefore provides a good balance between general applicability to arbitrarily shaped statistical dependencies in the data and shows accurate and robust performance when working with small sample sizes.


title: Optimal scheduling and science delivery of spectra for millions of targets in thousands of fields: the operational concept of the Maunakea spectroscopic explorer (MSE)

abstract: The Maunakea Spectroscopic Explorer (MSE) will each year obtain millions of spectra in the optical to near-infrared, at low (R~3000) to high (R~40000) spectral resolution by observing >3000 spectra per pointing via a highly multiplexed fiber-fed system. Key science programs for MSE include black hole reverberation mapping, stellar population analysis of faint galaxies at high redshift, and sub-km/s velocity accuracy for stellar astrophysics. The architecture of MSE is an assembly of subsystems designed to meet the science requirements and describes what MSE will look like. In this paper we focus on the operations concept of MSE, which describes how to operate a fiber fed, highly multiplexed, dedicated observatory given its architecture and the science requirements. The operations concept details the phases of operations, from selecting proposals within the science community to distributing back millions of spectra to this community. For each phase, the operations concept describes the tools required to support the science community in their analyses and the operations staff in their work. It also highlights the specific needs related to the complexity of MSE with millions of targets to observe, thousands of fibers to position, and different spectral resolution to use. Finally, the operations concept shows how the science requirements on calibration and observing efficiency can be met.


title: More on a Curious Nucleus

abstract: Harold Simmons introduced a pre-nucleus and its associated nucleus that measure the subfitness of a frame in Simmons [2010]. In this paper we continue the study if this pre-nucleus. We answer the questions posed in Simmons [2010].


title: On ideals generated by fold products of linear forms

abstract: Let $\mathbb K$ be a field of characteristic 0. Given $n$ linear forms in $R=\mathbb K[x_1,\ldots,x_k]$, with no two proportional, in one of our main results we show that the ideal $I\subset R$ generated by all $(n-2)$-fold products of these linear forms has linear graded free resolution. This result helps determining a complete set of generators of the symmetric ideal of $I$. Via Sylvester forms we can analyze from a different perspective the generators of the presentation ideal of the Orlik-Terao algebra of the second order; this is the algebra generated by the reciprocals of the products of any two (distinct) of the linear forms considered. We also show that when $k=2$, and when the collection of $n$ linear forms may contain proportional linear forms, for any $1\leq a\leq n$, the ideal generated by $a$-fold products of these linear forms has linear graded free resolution.


title: Canonical threefold singularities with a torus action of complexity one and $k$-empty polytopes

abstract: We classify the canonical threefold singularities that allow an effective two-torus action. This extends classification results of Mori on terminal threefold singularities and of Ishida and Iwashita on toric canonical threefold singularities. Our classification relies on lattice point emptiness of certain polytopal complexes with rational vertices. Scaling the polytopes by the least common multiple $k$ of the respective denominators, we investigate $k$-emptiness of polytopes with integer vertices. We show that two dimensional $k$-empty polytopes either are sporadic or come in series given by Farey sequences. We finally present the Cox ring iteration tree of the classified singularities, where all roots, i.e. all spectra of factorial Cox rings, are generalized compound du Val singularities.


title: signProx: One-Bit Proximal Algorithm for Nonconvex Stochastic Optimization

abstract: Stochastic gradient descent (SGD) is one of the most widely used optimization methods for parallel and distributed processing of large datasets. One of the key limitations of distributed SGD is the need to regularly communicate the gradients between different computation nodes. To reduce this communication bottleneck, recent work has considered a one-bit variant of SGD, where only the sign of each gradient element is used in optimization. In this paper, we extend this idea by proposing a stochastic variant of the proximal-gradient method that also uses one-bit per update element. We prove the theoretical convergence of the method for non-convex optimization under a set of explicit assumptions. Our results indicate that the compressed method can match the convergence rate of the uncompressed one, making the proposed method potentially appealing for distributed processing of large datasets.


title: Explaining Image Classifiers by Counterfactual Generation

abstract: When an image classifier makes a prediction, which parts of the image are relevant and why? We can rephrase this question to ask: which parts of the image, if they were not seen by the classifier, would most change its decision? Producing an answer requires marginalizing over images that could have been seen but weren't. We can sample plausible image in-fills by conditioning a generative model on the rest of the image. We then optimize to find the image regions that most change the classifier's decision after in-fill. Our approach contrasts with ad-hoc in-filling approaches, such as blurring or injecting noise, which generate inputs far from the data distribution, and ignore informative relationships between different parts of the image. Our method produces more compact and relevant saliency maps, with fewer artifacts compared to previous methods.


title: Maximising the sensitivity of next generation multi-object spectroscopy: system budget development and design optimizations for the Maunakea Spectroscopic Explorer

abstract: MSE is an 11.25m telescope with a 1.5 sq.deg. field of view. It can simultaneously obtain 3249 spectra at R=3000 from 360-1800nm, and 1083 spectra at R=40000 in the optical. Absolutely critical to the scientific success of MSE is to efficiently access the faint Universe. Here, we describe the adopted systems engineering methodology to ensure MSE meets the challenging sensitivity requirements, and how these requirements are partitioned across three budgets, relating to the throughput, noise and fiber injection efficiency. We then describe how the sensitivity of MSE as a system was estimated at the end of Conceptual Design Phase, and how this information was used to revisit the system design in order to meet the sensitivity requirements while maintaining the overall architectural concept of the Observatory. Finally, we present the anticipated sensitivity performance of MSE and describe the key science that these capabilities will enable.


title: TCP SYN Cookie Vulnerability

abstract: TCP SYN Cookies were implemented to mitigate against DoS attacks. It ensured that the server did not have to store any information for half-open connections. A SYN cookie contains all information required by the server to know the request is valid. However, the usage of these cookies introduces a vulnerability that allows an attacker to guess the initial sequence number and use that to spoof a connection or plant false logs.


title: Maunakea Spectroscopic Explorer (MSE): Implementing systems engineering methodology for the development of a new facility

abstract: Maunakea Spectroscopic Explorer will be a 10-m class highly multiplexed survey telescope, including a segmented primary mirror and robotic fiber positioners at the prime focus. MSE will replace the Canada France Hawaii Telescope (CFHT) on the summit of Mauna Kea, Hawaii. The multiplexing includes an array of over four thousand fibres feeding banks of spectrographs several tens of meters away. We present an overview of the requirements flow-down for MSE, from Science Requirements Document to Observatory Requirements Document. We have developed the system performance budgets, along with updating the budget architecture of our evolving project. We have also identified the links between subsystems and system budgets (and subsequently science requirements) and included system budget that are unique to MSE as a fiber-fed facility. All of this has led to a set of Observatory Requirements that is fully consistent with the Science Requirements.


title: On An Application of Hayashi's Inequality for Differentiable Functions

abstract: In this work, we show that the proof of the main result in [An Application of Hayashi's Inequality for Differentiable Functions, Computers & Mathematics with Applications, 32 (6) (1996), 95--99, by R.P. Agarwal and S.S. Dragomir] was wrong. A correction of the proof is given. More general related inequalities are also provided.


title: The science calibration challenges of next generation highly multiplexed optical spectroscopy: the case of the Maunakea Spectroscopic Explorer

abstract: MSE is an 11.25m telescope with a 1.5 sq.deg. field of view. It can simultaneously obtain 3249 spectra at R=3000 from 360-1800nm, and 1083 spectra at R=40000 in the optical. The large field of view, large number of targets, as well as the use of more than 4000 optical fibres to transport the light from the focal plane to the spectrographs, means that precise and accurate science calibration is difficult but essential to obtaining the science goals. As a large aperture telescope focusing on the faint Universe, precision sky subtraction and spectrophotometry are especially important. Here, we discuss the science calibration requirements, and the adopted calibration strategy, including operational features and hardware, that will enable the successful scientific exploitation of the vast MSE dataset.


title: Running on empty: Recharge dynamics from animal movement data

abstract: Vital rates such as survival and recruitment have always been important in the study of population and community ecology. At the individual level, physiological processes such as energetics are critical in understanding biomechanics and movement ecology and also scale up to influence food webs and trophic cascades. Although vital rates and population-level characteristics are tied with individual-level animal movement, most statistical models for telemetry data are not equipped to provide inference about these relationships because they lack the explicit, mechanistic connection to physiological dynamics. We present a framework for modeling telemetry data that explicitly includes an aggregated physiological process associated with decision making and movement in heterogeneous environments. Our framework accommodates a wide range of movement and physiological process specifications. We illustrate a specific model formulation in continuous-time to provide direct inference about gains and losses associated with physiological processes based on movement. Our approach can also be extended to accommodate auxiliary data when available. We demonstrate our model to infer mountain lion (in Colorado, USA) and African buffalo (in Kruger National Park, South Africa) recharge dynamics.


title: Sterile Neutrinos in String Derived Models

abstract: The MiniBooNE collaboration recently reported further evidence for the existence of sterile neutrinos, implying substantial mixing with the left-handed active neutrinos and at a comparable mass scale. I argue that while sterile neutrinos may arise naturally in large volume string models, they prove more of a challenge in heterotic-string models that replicate the Grand Unified Theory structure of the Standard Model matter states. Sterile neutrinos in heterotic-string models may imply the existence of an additional Abelian gauge symmetry of order 10-100TeV.


title: Out-of-equilibrium Gubser flow attractors

abstract: We discuss the non-equilibrium attractors of systems undergoing Gubser flow within kinetic theory by means of nonlinear dynamical systems. We obtain the attractors of anisotropic hydrodynamics, Israel-Stewart (IS) and transient fluid (DNMR) theories. These attractors are non-planar and the basin of attraction is three dimensional. We compare the asymptotic attractors of each hydrodynamic model with the one obtained from the exact Gubser solution of the Boltzmann equation within the relaxation time approximation. Anisotropic hydrodynamics matches, up to high numerical accuracy, the attractor of the exact theory while the other hydrodynamic theories fail to do so. Thus, anisotropic hydrodynamics is an effective theory for far-from-equilibrium fluids, which consists of the dissipative (nonperturbative) contributions at any order in the gradient expansion.


title: Expected observing efficiency of the Maunakea Spectroscopic Explorer (MSE)

abstract: The Maunakea Spectroscopic Explorer (MSE) will obtain millions of spectra each year in the optical to near-infrared, at low (R ~ 3,000) to high (R ~ 40,000) spectral resolution by observing >4,000 spectra per pointing via a highly multiplexed fiber-fed system. Key science programs for MSE include black hole reverberation mapping, stellar population analysis of faint galaxies at high redshift, and sub-km/s velocity accuracy for stellar astrophysics. One key metric of the success of MSE will be its survey speed, i.e. how many spectra of good signal-to-noise ratio will MSE be able to obtain every night and every year. This is defined at the higher level by the observing efficiency of the observatory and should be at least 80%, as indicated in the Science Requirements. In this paper we present the observing efficiency budget developed for MSE based on historical data at the Canada-France-Hawaii Telescope and other Maunakea Observatories. We describe the typical sequence of events at night to help us compute the observing efficiency and how we envision to optimize it to meet the science requirements.


title: Two novel approaches to the hadron-quark mixed phase in compact stars

abstract: First-order phase transitions, like the liquid-gas transition, proceed via formation of structures such as bubbles and droplets. In strongly interacting compact star matter, at the crust-core transition, but also at the hadron-quark transition in the core, these structures form different shapes dubbed "pasta phases". We describe two methods to obtain one-parameter families of hybrid equations of state (EoS) which mimic the thermodynamic behavior of pasta phases in between a low-density hadron and a high-density quark matter phase, thus generalizing the Maxwell construction. The first method replaces the behavior of pressure vs. chemical potential in a finite region around the critical %chemical potential pressure of the Maxwell construction by a polynomial interpolation. The second method uses extrapolations of the hadronic and quark matter EoS beyond the Maxwell point to define a mixing of both with weight functions bounded by finite limits around the Maxwell point. We apply both methods to the case of a hybrid EoS with a strong first order transition that entails the formation of a third family of compact stars and the corresponding mass twin phenomenon. We investigate for both models the robustness of this phenomenon against variation of the single parameter, the pressure increment at the critical chemical potential which quantifies the deviation from the Maxwell construction. We also show sets of results for other compact star observables than mass and radius, namely the moment of inertia and the baryon mass.


title: Modeling and budgeting fiber injection efficiency for the Maunakea spectroscopic explorer (MSE)

abstract: The Maunakea Spectroscopic Explorer (MSE) will each year obtain millions of spectra in the optical to near-infrared, at low (R ~ 3,000) to high (R ~ 40,000) spectral resolution by observing >4,000 spectra per pointing via a highly multiplexed fiber-fed system. Key science programs for MSE include black hole reverberation mapping, stellar population analysis of faint galaxies at high redshift, and sub-km/s velocity accuracy for stellar astrophysics. One key metric of the success of MSE will be its survey speed, i.e. how many spectra of good signal-to-noise ratio will MSE be able to obtain every night and every year. The survey speed is directly linked to the allocation efficiency - how many fibers in the focal surface can be allocated to targets - and to the injection efficiency - what fraction of light from a target can enter the fiber at the focal surface. In this paper we focus on the injection efficiency and how to optimize it to increase the signal-to-noise ratio of targets observed in sky dominated conditions. The injection efficiency depends on the size of the fiber and requires highly precise, repeatable and stable positioning of the fiber in the focal surface. We present the allocation budget used for Conceptual Design Review and the modeling that allows to estimate the injection efficiency, which is just one part necessary to meet the science requirements on sensitivities.


title: MSE FiTS: the ultimate multi-fiber optic transmission system

abstract: The Maunakea Spectroscopic Explorer (MSE) is a next-generation observatory, designed to provide highly multiplexed, multi-object spectroscopy over a wide field of view. The observatory will consist of (1) a telescope with an 11.25 m aperture, (2) a 1.5 square-degree science field of view, (3) fibre optic positioning and transmission systems, and (4) a suite of low (R=3000), moderate (R=6000) and high resolution (R=40,000) spectrographs. The Fibre Transmission System (FiTS) consists of 4332 optical fibres, designed to transmit the light from the telescope prime focus to the dedicated spectrographs. The ambitious science goals of MSE require the Fibre Transmission System to deliver performance well beyond the current state of the art for multi-fibre systems, e.g., the sensitivity to observe magnitude 24 objects over a very broad wavelength range (0.37 - 1.8 microns) while achieving relative spectrophotometric accuracy of <3% and radial velocity precision of 20 km/s.


title: Bipedal Hopping: Reduced-order Model Embedding via Optimization-based Control

abstract: This paper presents the design and validation of controlling hopping on the 3D bipedal robot Cassie. A spring-mass model is identified from the kinematics and compliance of the robot. The spring stiffness and damping are encapsulated by the leg length, thus actuating the leg length can create and control hopping behaviors. Trajectory optimization via direct collocation is performed on the spring-mass model to plan jumping and landing motions. The leg length trajectories are utilized as desired outputs to synthesize a control Lyapunov function based quadratic program (CLF-QP). Centroidal angular momentum, taking as an addition output in the CLF-QP, is also stabilized in the jumping phase to prevent whole body rotation in the underactuated flight phase. The solution to the CLF-QP is a nonlinear feedback control law that achieves dynamic jumping behaviors on bipedal robots with compliance. The framework presented in this paper is verified experimentally on the bipedal robot Cassie.


title: Additive and multiplicative effects network models

abstract: Network datasets typically exhibit certain types of statistical dependencies, such as within-dyad correlation, row and column heterogeneity, and third-order dependence patterns such as transitivity and clustering. The first two of these can be well-represented statistically with a social relations model, a type of additive random effects model originally developed for continuous dyadic data. Third-order patterns can be represented with multiplicative random effects models, which are related to matrix decompositions commonly used for matrix-variate data analysis. Additionally, these multiplicative random effects models generalize other popular latent variable network models, such as the stochastic blockmodel and the latent space model. In this article we review a general regression framework for the analysis of network data that combines these two types of random effects and accommodates a variety of network data types, including continuous, binary and ordinal network relations.


title: PGLasso: Microbial Community Detection through Phylogenetic Graphical Lasso

abstract: Due to the recent advances in high-throughput sequencing technologies, it becomes possible to directly analyze microbial communities in the human body and in the environment. Knowledge of how microbes interact with each other and form functional communities can provide a solid foundation to understand microbiome related diseases; this can serve as a key step towards precision medicine. In order to understand how microbes form communities, we propose a two step approach: First, we infer the microbial co-occurrence network by integrating a graph inference algorithm with phylogenetic information obtained directly from metagenomic data. Next, we utilize a network-based community detection algorithm to cluster microbes into functional groups where microbes in each group are highly correlated. We also curate a "gold standard" network based on the microbe-metabolic relationships which are extracted directly from the metagenomic data. Utilizing community detection on the resulting microbial metabolic pathway bipartite graph, the community membership for each microbe can be viewed as the true label when evaluating against other existing methods. Overall, our proposed framework Phylogenetic Graphical Lasso (PGLasso) outperforms existing methods with gains larger than 100% in terms of Adjusted Rand Index (ARI) which is commonly used to quantify the goodness of clusterings.


title: Analysis of a Reaction Diffusion Model for a Reservoir Supported Spread of Infectious Disease

abstract: Motivated by recent outbreaks of the Ebola Virus, we are concerned with the role that a vector reservoir plays in supporting the spatio-temporal spread of a highly lethal disease through a host population. In our context, the reservoir is a species capable of harboring and sustaining the pathogen. We develop models that describe the horizontal spread of the disease among the host population when the host population is in contact with the reservoir and when it is not in contact with the host population. These models are of reaction diffusion type, and they are analyzed, and their long term asymptotic behavior is determined.


title: Radiative Signatures of Plasmoid-Dominated Reconnection in Blazar Jets

abstract: The multi-wavelength spectral and temporal variability observed in blazars set tight constraints on current theoretical emission models. Here, we investigate the relativistic magnetic reconnection process as a source of blazar emission in which quasi-spherical plasmoids, containing relativistic particles and magnetic fields, are associated with the emission sites in blazar jets. By coupling recent two-dimensional particle-in-cell simulations of relativistic reconnection with a time-dependent radiative transfer code, we compute the non-thermal emission from a chain of plasmoids formed during a reconnection event. The derived photon spectra display characteristic features observed in both BL Lac sources and flat spectrum radio quasars, with the distinction made by varying the strength of the external photon fields, the jet magnetization, and the number of pairs per proton contained within. Light curves produced from reconnection events are composed of many fast and powerful flares that appear on excess of a slower evolving envelope produced by the cumulative emission of medium-sized plasmoids. The observed variability is highly dependent upon the orientation of the reconnection layer with respect to the blazar jet axis and to the observer. Our model provides a physically motivated framework for explaining the multi-timescale blazar variability across the entire electromagnetic spectrum.