title: A Fog Computing Framework for Autonomous Driving Assist: Architecture, Experiments, and Challenges

abstract: Autonomous driving is expected to provide a range of far-reaching economic, environmental and safety benefits. In this study, we propose a fog computing based framework to assist autonomous driving. Our framework relies on overhead views from cameras and data streams from vehicle sensors to create a network of distributed digital twins, called an edge twin, on fog machines. The edge twin will be continuously updated with the locations of both autonomous and human-piloted vehicles on the road segments. The vehicle locations will be harvested from overhead cameras as well as location feeds from the vehicles themselves. Although the edge twin can make fair road space allocations from a global viewpoint, there is a communication cost (delay) in reaching it from the cameras and vehicular sensors. To address this, we introduce a machine learning forecaster as a part of the edge twin which is responsible for predicting the future location of vehicles. Lastly, we introduce a box algorithm that will use the forecasted values to create a hazard map for the road segment which would be used by the framework to suggest safe manoeuvres for the autonomous vehicles such as lane changes and accelerations. We present the complete fog computing framework for autonomous driving assist and evaluate key portions of the proposed framework using simulations based on a real-world dataset of vehicle position traces on a highway


title: Latent Function Decomposition for Forecasting Li-ion Battery Cells Capacity: A Multi-Output Convolved Gaussian Process Approach

abstract: A latent function decomposition method is proposed for forecasting the capacity of lithium-ion battery cells. The method uses the Multi-Output Gaussian Process, a generative machine learning framework for multi-task and transfer learning. The MCGP decomposes the available capacity trends from multiple battery cells into latent functions. The latent functions are then convolved over kernel smoothers to reconstruct and/or forecast capacity trends of the battery cells. Besides the high prediction accuracy the proposed method possesses, it provides uncertainty information for the predictions and captures nontrivial cross-correlations between capacity trends of different battery cells. These two merits make the proposed MCGP a very reliable and practical solution for applications that use battery cell packs. The MCGP is derived and compared to benchmark methods on an experimental lithium-ion battery cells data. The results show the effectiveness of the proposed method.


title: Signal Processing on PV Time-Series Data: Robust Degradation Analysis without Physical Models

abstract: A novel unsupervised machine learning approach for analyzing time-series data is applied to the topic of photovoltaic (PV) system degradation rate estimation, sometimes referred to as energy-yield degradation analysis. This approach only requires a measured power signal as an input--no irradiance data, temperature data, or system configuration information. We present results on a data set that was previously analyzed and presented by NREL using RdTools, validating the accuracy of the new approach and showing increased robustness to data anomalies while reducing the data requirements to carry out the analysis.


title: A GN-model closed-form formula considering coherency terms in the Link function and covering all possible islands in 2-D GN integration

abstract: Efficient evaluation of nonlinearity in modern coherent optical telecommunication networks is necessary for design, optimization and management. The original GN model formula can not be used in real time applications due to long processing time needed for numerical 2-D integration, particularly when facing with full C band or C+L band. There are closed form approximations of the GN model which has been mathematically derived based on some assumptions which fail in low dispersion regime. In this work, we derive a GN closed form formula without any restrictive assumption capable of handling arbitrary WDM combs and arbitrary link structures in closed form which can be used with a very good accuracy with respect to original GN formula for networks employing both high and low dispersion fibers.


title: A Stochastic Model for Uncontrolled Charging of Electric Vehicles Using Cluster Analysis

abstract: This paper proposes a probabilistic model for uncontrolled charging of electric vehicles (EVs). EV charging will add significant load to power systems in the coming years and, due to the convenience of charging at home, this is likely to occur in residential distribution systems. Estimating the size and shape of the load will allow necessary reinforcements to be identified. Models predicting EV charging are usually based on data from travel surveys, or from small trials. Travel surveys are recorded by hand and typically describe conventional vehicles, but represent a much larger and more diverse sample of the population. The model here utilizes both sources: trial data to parameterize the model, and survey data as the model input. Clustering is used to identify modes of vehicle use, thus reducing vehicle use to a single parameter -- which can be incorporated into the model without adding significant computational burden. Two case studies are included: one investigating the aggregated charging of 50 vehicles, and one predicting the increase in after diversity maximum demand for different regions of the UK.


title: Galactic Gas Flows from Halo to Disk: Tomography and Kinematics at the Milky Way's Disk-Halo Interface

abstract: We present a novel absorption line survey using 54 blue horizontal branch stars (BHBs) in the Milky Way halo as background sources for detecting gas flows at the disk-halo interface. Distance measurements to high-latitude ($b$ > 60{\deg}) background stars at 3.1-13.4 kpc, combined with unprecedented spatial sampling and spectral resolution, allow us to examine the 3-dimensional spatial distribution and kinematics of gas flows near the disk. We detect absorption signatures of extraplanar CaII and NaI in Keck HIRES spectra and find that their column densities exhibit no trend with distance to the background sources, indicating that these clouds lie within 3.1 kpc of the disk. We calculate covering fractions of $f_{CaII}$ = 63%, $f_{NaI}$ = 26%, and $f_{HI}$ = 52%, consistent with a picture of the CGM that includes multi-phase clouds containing small clumps of cool gas within hotter, more diffuse gas. Our measurements constrain the scale of any substructure within these cool clouds to <0.5 kpc. CaII and NaI absorption features exhibit an intermediate-velocity (IV) component inflowing at velocities of -75 km/s < v < -25 km/s relative to the local standard of rest, consistent with previously-studied HI structures in this region. We report the new detection of an inflow velocity gradient ${\Delta}v_z$ ~ 6-9 km/s/kpc across the Galactic plane. These findings place constraints on the physical and kinematic properties of CGM gas flows through the disk-halo interface, and support a galactic fountain model in which cold gas rains back onto the disk.


title: From Gauged Linear Sigma Models to Geometric Representation of $\mathbb{WCP}(N,\tilde{N})$ in 2D

abstract: In this paper two issues are addressed. First, we discuss renormalization properties of a class of gauged linear sigma models (GLSM) which reduce to $\mathbb{WCP}(N,\tilde{N})$ non-linear sigma models (NLSM) in the low-energy limit. Sometimes they are referred to as the Hanany-Tong models. If supersymmetry is ${\cal N} =(2,2)$ the ultraviolet-divergent logarithm in LGSM appears, in the renormalization of the Fayet-Iliopoulos parameter, and is exhausted by a single tadpole graph. This is not the case in the daughter NLSMs. As a result, the one-loop renormalizations are different in GLSMs and their daughter NLSMs We explain this difference and identify its source. In particular, we show why at $N=\tilde N$ there is no UV logarithms in the parent GLSM, while they do appear on the corresponding NLSM does not vanish. In the second part of the paper we discuss the same problem for a class of ${\cal N} =(0,2)$ GLSMs considered previously. In this case renormalization is not limited to one loop; all-orders exact $\beta$ functions for GLSMs are known. We discuss divergent loops at one and two-loop levels.


title: Dependence of the static quark free energy on $\mu_B$ and the crossover temperature of $N_f = 2+1$ QCD

abstract: We study the dependence of the static quark free energy on the baryon chemical potential for $N_f = 2+1$ QCD with physical quark masses, in a range of temperature spanning from 120~MeV up to 1~GeV and adopting a stout staggered discretization with two different values of the Euclidean temporal extension, $N_t = 6$ and $N_t = 8$. In order to deal with the sign problem, we exploit both Taylor expansion and analytic continuation, obtaining consistent results. We show that the dependence of the free energy on $\mu_B$ is sensitive to the location of the chiral crossover, in particular the $\mu_B$-susceptibility, i.e. the linear term in $\mu_B^2$ in the Taylor expansion of the free energy, has a peak around 150 MeV. We also discuss the behavior expected in the high temperature regime based on perturbation theory, and obtain a good quantitative agreement with numerical results.


title: On spectral spread of generalized distance matrix of a graph

abstract: For a simple connected graph $G$, let $D(G)$, $Tr(G)$, $D^{L}(G)$ and $D^{Q}(G)$, respectively be the distance matrix, the diagonal matrix of the vertex transmissions, distance Laplacian matrix and the distance signless Laplacian matrix of a graph $G$. The convex linear combinations $D_{\alpha}(G)$ of $Tr(G)$ and $D(G)$ is defined as $D_{\alpha}(G)=\alpha Tr(G)+(1-\alpha)D(G)$, $0\leq \alpha\leq 1$. As $D_{0}(G)=D(G), ~~~ 2D_{\frac{1}{2}}(G)=D^{Q}(G), ~~~ D_{1}(G)=Tr(G)$ and $D_{\alpha}(G)-D_{\beta}(G)=(\alpha-\beta)D^{L}(G)$, this matrix reduces to merging the distance spectral, distance Laplacian spectral and distance signless Laplacian spectral theories. Let $\partial_{1}(G)\geq \partial_{2}(G)\geq \dots \geq \partial_{n}(G)$ be the eigenvalues of $D_{\alpha}(G)$ and let $D_{\alpha}S(G)=\partial_{1}(G)-\partial_{n}(G)$ be the generalized distance spectral spread of the graph $G$. In this paper, we obtain some bounds for the generalized distance spectral spread $D_{\alpha}(G)$. We also obtain relation between the generalized distance spectral spread $D_{\alpha}(G)$ and the distance spectral spread $S_{D}(G)$. Further, we obtain the lower bounds for $D_{\alpha}S(G)$ of bipartite graphs involving different graph parameters and we characterize the extremal graphs for some cases. We also obtain lower bounds for $D_{\alpha}S(G)$ in terms of clique number and independence number of the graph $G$ and characterize the extremal graphs for some cases.


title: Ellis wormhole without a phantom scalar field

abstract: In this paper, we present an exact solution for $(3+1)$-dimensional Einstein-scalar-Gauss-Bonnet theory (EsGB) in electrovacuum. The solution is characterized by only one parameter, $Q$, which in general can be associated with the electromagnetic field and the scalar field. We show that the solution corresponds to a charged wormhole with throat at the region $r = |Q|$, and is also supported by a real scalar field having a positive kinetic term. We show that the solution belongs to the most general class of solutions known as Ellis wormholes but without the need for `exotic matter' or a phantom scalar field.


title: Flat Littlewood Polynomials Exist

abstract: We show that there exist absolute constants $\Delta > \delta > 0$ such that, for all $n \geqslant 2$, there exists a polynomial $P$ of degree $n$, with $\pm 1$ coefficients, such that $$\delta\sqrt{n} \leqslant |P(z)| \leqslant \Delta\sqrt{n}$$ for all $z\in\mathbb{C}$ with $|z|=1$. This confirms a conjecture of Littlewood from 1966.


title: Inhomogeneous broadening of optical transitions observed in photoluminescence and modulated reflectance of polar and non-polar InGaN quantum wells

abstract: In this work the broadening of interband transitions in InGaN/GaN quantum wells (QWs) resulting from structural inhomogeneities is analyzed. The role of polarization-induced electric field in the mechanism behind the inhomogeneous broadening observed in photoluminescence (PL) and electromodulated reflectance (ER) spectra of InGaN QWs dedicated to green/blue lasers is explained. Spectra of both polar and non-polar QWs are simulated within the random QW model distinguishing contributions of individual transitions taking into account QW inhomogeneities (QW width and indium content fluctuations). On this basis we interpret the ER and PL spectra measured for a polar multiple QW InGaN/GaN structure. The built-in electric field shifts the emission wavelength to red, and enhances the broadening of optical transitions. It is clearly shown that for polar QWs the Stokes shift can be easily overestimated if PL spectra are compared with ER spectra since the intensity of the fundamental transition observed in ER spectra significantly decreases with the increase in QW width. In this way ER signal related to excited states start to dominate. This effect is strongly enhanced by QW inhomogeneities.


title: An Actor-Critic-Attention Mechanism for Deep Reinforcement Learning in Multi-view Environments

abstract: In reinforcement learning algorithms, leveraging multiple views of the environment can improve the learning of complicated policies. In multi-view environments, due to the fact that the views may frequently suffer from partial observability, their level of importance are often different. In this paper, we propose a deep reinforcement learning method and an attention mechanism in a multi-view environment. Each view can provide various representative information about the environment. Through our attention mechanism, our method generates a single feature representation of environment given its multiple views. It learns a policy to dynamically attend to each view based on its importance in the decision-making process. Through experiments, we show that our method outperforms its state-of-the-art baselines on TORCS racing car simulator and three other complex 3D environments with obstacles. We also provide experimental results to evaluate the performance of our method on noisy conditions and partial observation settings.


title: Arena: a toolkit for Multi-Agent Reinforcement Learning

abstract: We introduce Arena, a toolkit for multi-agent reinforcement learning (MARL) research. In MARL, it usually requires customizing observations, rewards and actions for each agent, changing cooperative-competitive agent-interaction, and playing with/against a third-party agent, etc. We provide a novel modular design, called Interface, for manipulating such routines in essentially two ways: 1) Different interfaces can be concatenated and combined, which extends the OpenAI Gym Wrappers concept to MARL scenarios. 2) During MARL training or testing, interfaces can be embedded in either wrapped OpenAI Gym compatible Environments or raw environment compatible Agents. We offer off-the-shelf interfaces for several popular MARL platforms, including StarCraft II, Pommerman, ViZDoom, Soccer, etc. The interfaces effectively support self-play RL and cooperative-competitive hybrid MARL. Also, Arena can be conveniently extended to your own favorite MARL platform.


title: Heuristic solutions to robust variants of the minimum-cost integer flow problem

abstract: This paper deals with robust optimization applied to network flows. Two robust variants of the minimum-cost integer flow problem are considered. Thereby, uncertainty in problem formulation is limited to arc unit costs and expressed by a finite set of explicitly given scenarios. It is shown that both problem variants are NP-hard. To solve the considered variants, several heuristics based on local search or evolutionary computing are proposed. The heuristics are experimentally evaluated on appropriate problem instances.


title: On-the-fly particle metrology in hollow-core photonic crystal fibre

abstract: Efficient monitoring of airborne particulate matter (PM), especially particles with aerodynamic diameter less than 2.5 um (PM2.5), is crucial for improving public health. Reliable information on the concentration, size distribution and chemical characteristics of PMs is key to evaluating air pollution and identifying its sources. Standard methods for PM2.5 characterization require sample collection from the atmosphere and post-analysis using sophisticated equipment in a laboratory environment, and are normally very time-consuming. Although optical methods based on analysis of scattering of free-space laser beams or evanescent fields are in principle suitable for real-time particle counting and sizing, lack of knowledge of the refractive index in these methods not only leads to inevitable sizing ambiguity but also prevents identification of the particle material. In the case of evanescent wave detection, the system lifetime is strongly limited by adhesion of particles to the surfaces. Here we report a novel technique for airborne particle metrology based on hollow-core photonic crystal fibre. It offers in situ particle counting, sizing and refractive index measurement with effectively unlimited device lifetime, and relies on optical forces that automatically capture airborne particles in front of the hollow core and propel them into the fibre. The resulting transmission drop, together with the time-of-flight of the particles passing through the fibre, provide unambiguous mapping of particle size and refractive index with high accuracy. The technique represents a considerable advance over currently available real-time particle metrology systems, and can be directly applied to monitoring air pollution in the open atmosphere as well as precise particle characterization in a local environment such as a closed room or a reaction vessel.


title: Characterizing Attacks on Deep Reinforcement Learning

abstract: Recent studies show that Deep Reinforcement Learning (DRL) models are vulnerable to adversarial attacks, which attack DRL models by adding small perturbations to the observations. However, some attacks assume full availability of the victim model, and some require a huge amount of computation, making them less feasible for real world applications. In this work, we make further explorations of the vulnerabilities of DRL by studying other aspects of attacks on DRL using realistic and efficient attacks. First, we adapt and propose efficient black-box attacks when we do not have access to DRL model parameters. Second, to address the high computational demands of existing attacks, we introduce efficient online sequential attacks that exploit temporal consistency across consecutive steps. Third, we explore the possibility of an attacker perturbing other aspects in the DRL setting, such as the environment dynamics. Finally, to account for imperfections in how an attacker would inject perturbations in the physical world, we devise a method for generating a robust physical perturbations to be printed. The attack is evaluated on a real-world robot under various conditions. We conduct extensive experiments both in simulation such as Atari games, robotics and autonomous driving, and on real-world robotics, to compare the effectiveness of the proposed attacks with baseline approaches. To the best of our knowledge, we are the first to apply adversarial attacks on DRL systems to physical robots.


title: Model Adaptation via Model Interpolation and Boosting for Web Search Ranking

abstract: This paper explores two classes of model adaptation methods for Web search ranking: Model Interpolation and error-driven learning approaches based on a boosting algorithm. The results show that model interpolation, though simple, achieves the best results on all the open test sets where the test data is very different from the training data. The tree-based boosting algorithm achieves the best performance on most of the closed test sets where the test data and the training data are similar, but its performance drops significantly on the open test sets due to the instability of trees. Several methods are explored to improve the robustness of the algorithm, with limited success.


title: Learning Probabilities: Towards a Logic of Statistical Learning

abstract: We propose a new model for forming beliefs and learning about unknown probabilities (such as the probability of picking a red marble from a bag with an unknown distribution of coloured marbles). The most widespread model for such situations of 'radical uncertainty' is in terms of imprecise probabilities, i.e. representing the agent's knowledge as a set of probability measures. We add to this model a plausibility map, associating to each measure a plausibility number, as a way to go beyond what is known with certainty and represent the agent's beliefs about probability. There are a number of standard examples: Shannon Entropy, Centre of Mass etc. We then consider learning of two types of information: (1) learning by repeated sampling from the unknown distribution (e.g. picking marbles from the bag); and (2) learning higher-order information about the distribution (in the shape of linear inequalities, e.g. we are told there are more red marbles than green marbles). The first changes only the plausibility map (via a 'plausibilistic' version of Bayes' Rule), but leaves the given set of measures unchanged; the second shrinks the set of measures, without changing their plausibility. Beliefs are defined as in Belief Revision Theory, in terms of truth in the most plausible worlds. But our belief change does not comply with standard AGM axioms, since the revision induced by (1) is of a non-AGM type. This is essential, as it allows our agents to learn the true probability: we prove that the beliefs obtained by repeated sampling converge almost surely to the correct belief (in the true probability). We end by sketching the contours of a dynamic doxastic logic for statistical learning.


title: What Are Observables in Hamiltonian Einstein-Maxwell Theory?

abstract: Is change missing in Hamiltonian Einstein-Maxwell theory? Given the most common definition of observables (having weakly vanishing Poisson bracket with each first-class constraint), observables are constants of the motion and nonlocal. Unfortunately this definition also implies that the observables for massive electromagnetism with gauge freedom (Stueckelberg) are inequivalent to those of massive electromagnetism without gauge freedom (Proca). The alternative Pons-Salisbury-Sundermeyer definition of observables, aiming for Hamiltonian-Lagrangian equivalence, uses the gauge generator G, a tuned sum of first-class constraints, rather than each first-class constraint separately, and implies equivalent observables for equivalent massive electromagnetisms. For General Relativity, G generates 4-dimensional Lie derivatives for solutions. The Lie derivative compares different space-time points with the same coordinate value in different coordinate systems, like 1 a.m. summer time vs. 1 a.m. standard time, so a vanishing Lie derivative implies constancy rather than covariance. Requiring equivalent observables for equivalent formulations of massive gravity confirms that $G$ must generate the $4$-dimensional Lie derivative (not $0$) for observables. These separate results indicate that observables are invariant under internal gauge symmetries but covariant under external gauge symmetries, but can this bifurcated definition work for mixed theories such as Einstein-Maxwell theory? Pons, Salisbury and Shepley have studied $G$ for Einstein-Yang-Mills. For Einstein-Maxwell, both $F_{\mu\nu}$ and $g_{\mu\nu}$ are invariant under electromagnetic gauge transformations and covariant (changing by a Lie derivative) under $4$-dimensional coordinate transformations. Using the bifurcated definition, these quantities count as observables, as one would expect on non-Hamiltonian grounds.


title: Design of one-year mortality forecast at hospital admission based: a machine learning approach

abstract: Background: Palliative care is referred to a set of programs for patients that suffer life-limiting illnesses. These programs aim to guarantee a minimum level of quality of life (QoL) for the last stage of life. They are currently based on clinical evaluation of risk of one-year mortality. Objectives: The main objective of this work is to develop and validate machine-learning based models to predict the exitus of a patient within the next year using data gathered at hospital admission. Methods: Five machine learning techniques were applied in our study to develop machine-learning predictive models: Support Vector Machines, K-neighbors Classifier, Gradient Boosting Classifier, Random Forest and Multilayer Perceptron. All models were trained and evaluated using the retrospective dataset. The evaluation was performed with five metrics computed by a resampling strategy: Accuracy, the area under the ROC curve, Specificity, Sensitivity, and the Balanced Error Rate. Results: All models for forecasting one-year mortality achieved an AUC ROC from 0.858 to 0.911. Specifically, Gradient Boosting Classifier was the best model, producing an AUC ROC of 0.911 (CI 95%, 0.911 to 0.912), a sensitivity of 0.858 (CI 95%, 0.856 to 0.86) and a specificity of 0.807 (CI 95%, 0.806 to 0808) and a BER of 0.168 (CI 95%, 0.167 to 0.169). Conclusions: The analysis of common information at hospital admission combined with machine learning techniques produced models with competitive discriminative power. Our models reach the best results reported in state of the art. These results demonstrate that they can be used as an accurate data-driven palliative care criteria inclusion.


title: Deep Reinforcement Learning for Clinical Decision Support: A Brief Survey

abstract: Owe to the recent advancements in Artificial Intelligence especially deep learning, many data-driven decision support systems have been implemented to facilitate medical doctors in delivering personalized care. We focus on the deep reinforcement learning (DRL) models in this paper. DRL models have demonstrated human-level or even superior performance in the tasks of computer vision and game playings, such as Go and Atari game. However, the adoption of deep reinforcement learning techniques in clinical decision optimization is still rare. We present the first survey that summarizes reinforcement learning algorithms with Deep Neural Networks (DNN) on clinical decision support. We also discuss some case studies, where different DRL algorithms were applied to address various clinical challenges. We further compare and contrast the advantages and limitations of various DRL algorithms and present a preliminary guide on how to choose the appropriate DRL algorithm for particular clinical applications.


title: The evolution of CNO isotopes: the impact of massive stellar rotators

abstract: Chemical abundances and abundance ratios measured in galaxies provide precious information about the mechanisms, modes and time scales of the assembly of cosmic structures. Yet, the nucleogenesis and chemical evolution of elements heavier than helium are dictated mostly by the physics of the stars and the shape of the stellar mass spectrum. In particular, estimates of CNO isotopic abundances in the hot, dusty media of high-redshift starburst galaxies offer a unique glimpse into the shape of the stellar initial mass function (IMF) in extreme environments that can not be accessed with direct observations (star counts). Underlying uncertainties in stellar evolution and nucleosynthesis theory, however, may hurt our chances of getting a firm grasp of the IMF in these galaxies. In this work, we adopt new yields for massive stars, covering different initial rotational velocities. First, we implement the new yield set in a well-tested chemical evolution model for the Milky Way. The calibrated model is then adapted to the specific case of a prototype submillimeter galaxy (SMG). We show that, if the formation of fast-rotating stars is favoured in the turbulent medium of violently star-forming galaxies irrespective of metallicity, the IMF needs to be skewed towards high-mass stars in order to explain the CNO isotope ratios observed in SMGs. If, instead, stellar rotation becomes negligible beyond a given metallicity threshold, as is the case for our own Galaxy, there is no need to invoke a top-heavy IMF in starbursts.


title: Multiple block sizes and overlapping blocks for multivariate time series extremes

abstract: Block maxima methods constitute a fundamental part of the statistical toolbox in extreme value analysis. However, most of the corresponding theory is derived under the simplifying assumption that block maxima are independent observations from a genuine extreme value distribution. In practice however, block sizes are finite and observations from different blocks are dependent. Theory respecting the latter complications is not well developed, and, in the multivariate case, has only recently been established for disjoint blocks of a single block size. We show that using overlapping blocks instead of disjoint blocks leads to a uniform improvement in the asymptotic variance of the multivariate empirical distribution function of rescaled block maxima and any smooth functionals thereof (such as the empirical copula), without any sacrifice in the asymptotic bias. We further derive functional central limit theorems for multivariate empirical distribution functions and empirical copulas that are uniform in the block size parameter, which seems to be the first result of this kind for estimators based on block maxima in general. The theory allows for various aggregation schemes over multiple block sizes, leading to substantial improvements over the single block length case and opens the door to further methodology developments. In particular, we consider bias correction procedures that can improve the convergence rates of extreme-value estimators and shed some new light on estimation of the second-order parameter when the main purpose is bias correction.


title: Context-Aware Convolutional Neural Network for Grading of Colorectal Cancer Histology Images

abstract: Digital histology images are amenable to the application of convolutional neural network (CNN) for analysis due to the sheer size of pixel data present in them. CNNs are generally used for representation learning from small image patches (e.g. 224x224) extracted from digital histology images due to computational and memory constraints. However, this approach does not incorporate high-resolution contextual information in histology images. We propose a novel way to incorporate larger context by a context-aware neural network based on images with a dimension of 1,792x1,792 pixels. The proposed framework first encodes the local representation of a histology image into high dimensional features then aggregates the features by considering their spatial organization to make a final prediction. The proposed method is evaluated for colorectal cancer grading and breast cancer classification. A comprehensive analysis of some variants of the proposed method is presented. Our method outperformed the traditional patch-based approaches, problem-specific methods, and existing context-based methods quantitatively by a margin of 3.61%. Code and dataset related information is available at this link: https://tia-lab.github.io/Context-Aware-CNN


title: Cyanoacetylene in the outflow/hot molecular core G331.512-0.103

abstract: Using APEX-1 and APEX-2 observations, we have detected and studied the rotational lines of the HC$_3$N molecule (cyanoacetylene) in the powerful outflow/hot molecular core G331.512-0.103. We identified thirty-one rotational lines at $J$ levels between 24 and 39; seventeen of them in the ground vibrational state $v$=0 (9 lines corresponding to the main C isotopologue and 8 lines corresponding to the $^{13}$C isotopologues), and fourteen in the lowest vibrationally excited state $v_7$=1. Using LTE-based population diagrams for the beam-diluted $v$=0 transitions, we determined $T_{\rm exc}$=85$\pm$4 K and $N$(HC$_3$N)=(6.9$\pm$0.8)$\times$10$^{14}$ cm$^{-2}$, while for the beam-diluted $v_7$=1 transitions we obtained $T_{\rm exc}$=89$\pm$10 K and $N$(HC$_3$N)=2$\pm$1$\times$10$^{15}$ cm$^{-2}$. Non-LTE calculations using H$_2$ collision rates indicate that the HC$_3$N emission is in good agreement with LTE-based results. From the non-LTE method we estimated $T_{\rm kin}$ $\simeq$90~K, $n$(H$_2$)$\simeq$2$\times$10$^7$~cm$^{-3}$ for a central core of 6 arcsec in size. A vibrational temperature in the range from 130~K to 145~K was also determined, values which are very likely lower limits. Our results suggest that rotational transitions are thermalized, while IR radiative pumping processes are probably more efficient than collisions in exciting the molecule to the vibrationally excited state $v_7$=1. Abundance ratios derived under LTE conditions for the $^{13}$C isotopologues suggest that the main formation pathway of HC$_3$N is ${\rm C}_2{\rm H}_2 + {\rm CN} \rightarrow {\rm HC}_3{\rm N} + {\rm H}$.


title: EXOFASTv2: A public, generalized, publication-quality exoplanet modeling code

abstract: We present the next generation public exoplanet fitting software, EXOFASTv2. It is capable of fitting an arbitrary number of planets, radial velocity data sets, astrometric data sets, and/or transits observed with any combination of wavelengths. We model the star simultaneously in the fit and provide several state-of-the-art ways to constrain its properties, including taking advantage of the now-ubiquitous all-sky catalog photometry and Gaia parallaxes. EXOFASTv2 can model the star by itself, too. Multi-planet systems are modeled self-consistently with the same underlying stellar mass that defines their semi-major axes through Kepler's law and the planetary period. Transit timing, duration, and depth variations can be modeled with a simple command line option. We explain our methodology and rationale as well as provide an improved version of the core transit model that is both 25\% faster and more accurate. We highlight several potential pitfalls in exoplanet modeling, including the handling of eccentricity in transit-only fits, that the standard exoplanet convention for $\omega$ uses a left-handed coordinate system, contrary to most modern textbooks, how to avoid an important degeneracy when allowing negative companion masses, and a widely unappreciated, potential 10-minute ambiguity in the reported transit times. EXOFASTv2 is available at https://github.com/jdeast/EXOFASTv2 . The code is written in IDL, and includes an executable that can be run freely and legally without an IDL license or any knowledge of the language. Extensive documentation and tutorials are included in the distribution for a variety of example fits. Advanced amateurs and undergrads have successfully performed sophisticated global fits of complex planetary systems with EXOFASTv2. It is therefore a powerful tool for education and outreach as well as the broader professional community.


title: The Southern Stellar Stream Spectroscopic Survey (${S}^5$): Overview, Target Selection, Data Reduction, Validation, and Early Science

abstract: We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}.


title: Fibration structure in toric hypersurface Calabi-Yau threefolds

abstract: We find through a systematic analysis that all but 29,223 of the 473.8 million 4D reflexive polytopes found by Kreuzer and Skarke have a 2D reflexive subpolytope. Such a subpolytope is generally associated with the presence of an elliptic or genus one fibration in the corresponding birational equivalence class of Calabi-Yau threefolds. This extends the growing body of evidence that most Calabi-Yau threefolds have an elliptically fibered phase.


title: Fossil group origins X. Velocity segregation in fossil systems

abstract: We want to study how the velocity segregation and the radial profile of the velocity dispersion depend on the prominence of the brightest cluster galaxies (BCGs). We divide a sample of 102 clusters and groups of galaxies into four bins of magnitude gap between the two brightest cluster members. We then compute the velocity segregation in bins of absolute and relative magnitudes. Moreover, for each bin of magnitude gap we compute the radial profile of the velocity dispersion. When using absolute magnitudes, the segregation in velocity is limited to the two brightest bins and no significant difference is found for different magnitude gaps. However, when we use relative magnitudes, a trend appears in the brightest bin: the larger the magnitude gap, the larger the velocity segregation. We also show that this trend is mainly due to the presence, in the brightest bin, of satellite galaxies in systems with small magnitude gaps: in fact, if we study separately central galaxies and satellites, this trend is mitigated and central galaxies are more segregated than satellites for any magnitude gap. A similar result is found in the radial velocity dispersion profiles: a trend is visible in central regions (where the BCGs dominate) but, if we analyse the profile using satellites alone, the trend disappears. In the latter case, the shape of the velocity dispersion profile in the centre of systems with different magnitude gaps show three types of behaviours: systems with the smallest magnitude gaps have an almost flat profile from the centre to the external regions; systems with the largest magnitude gaps show a monothonical growth from the low values of the central part to the flat ones in the external regions; finally, systems with $1.0 < \Delta m_{12} \le 1.5$ show a profile that peaks in the centres and then decreases towards the external regions. We suggest that two mechanisms could be respons....


title: Limit on the LMC mass from a census of its satellites

abstract: We study the orbits of ultra-faint dwarf galaxies in the combined presence of the Milky Way and LMC and we find 6 dwarfs which were likely accreted with the LMC (Car 2, Car 3, Hor 1, Hyi 1, Phe 2, Ret 2), in addition to the SMC, representing strong evidence of dwarf galaxy group infall. This procedure depends on the gravitational pull of the LMC, thus allowing us to place a lower bound on the Cloud's mass of $M_{\rm LMC} > 1.24\times10^{11} M_\odot$. This mass estimate is validated by applying the technique to a cosmological zoom-in simulation of a Milky Way-like galaxy with an LMC analogue where we find that while this lower bound may be overestimated, it will improve in the future with smaller observational errors. We apply this technique to dwarf galaxies lacking radial velocities and find that Eri 3 has a broad range of radial velocities for which it has a significant chance ($> 0.4$) of having being bound to the Cloud. We study the non-Magellanic classical satellites and find that Fornax has an appreciable probability of being an LMC satellite if the LMC is sufficiently massive. In addition, we explore how the orbits of the Milky Way satellites change in the presence of the LMC and find a significant change for several objects. Finally, we find that the LMC satellites are slightly smaller than the Milky Way satellites at a fixed luminosity, possibly due to the different tidal environments they have experienced.


title: Superconducting order of $\mathrm{Sr}_2\mathrm{RuO}_4$ from a three-dimensional microscopic model

abstract: We compute and compare even- and odd-parity superconducting order parameters of strontium ruthenate ($\mathrm{Sr}_2\mathrm{RuO}_4$) in the limit of weak interactions, resulting from a fully microscopic three-dimensional model including spin-orbit coupling. We find that odd-parity helical and even-parity $d$-wave order are favored for smaller and larger values of the Hund's coupling parameter $J$, respectively. Both orders are found compatible with specific heat data and the recently-reported nuclear magnetic resonance (NMR) Knight shift drop [A. Pustogow et al. Nature 574, 72 (2019)]. The chiral $p$-wave order, numerically very competitive with helical order, sharply conflicts with the NMR experiment.


title: Pseudo-Evolution of Galaxies in LCDM Cosmology

abstract: Our knowledge about galaxy evolution comes from transforming observed galaxy properties at different redshifts to co-moving physical scales. This transformation depends on using a cosmological model. Here the effects of unintentional mixing of two different cosmological models on the size evolution of galaxies is studied. As a gedanken experiment, a galaxy of fixed proper size and luminosity is moved across different redshifts. The apparent size of this galaxy is then interpreted with a cosmological model presumed by the observer, which is different compared to the cosmology exhibited by the Universe. In such a case, a spurious size evolution of the galaxy is observed. A galaxy behaving according to the R_h=ct and Neumann's cosmology, when interpreted with the LCDM cosmological model, shows an increase in size by a factor of 1.1 and 1.3 from z=7.5 to z approx. 0, respectively. The apparent size of a galaxy in a static Euclidean cosmology, when interpreted in the LCDM model, shows a factor of 23.8 increase in size between z=7.5 to z approx. 0. This is in close agreement with the observational data with a size increase of a factor of 6.8 between z=3.2 to z approx. 0. Furthermore, using the apparent size data, it is shown that the difference between the derived proper sizes in R_h=ct, Neumann's and LCDM cosmological models are minimal.


title: Barycentric interpolation on Riemannian and semi-Riemannian spaces

abstract: Interpolation of data represented in curvilinear coordinates and possibly having some non-trivial, typically Riemannian or semi-Riemannian geometry is an ubiquitous task in all of physics. In this work we present a covariant generalization of the barycentric coordinates and the barycentric interpolation method for Riemannian and semi-Riemannian spaces of arbitrary dimension. We show that our new method preserves the linear accuracy property of barycentric interpolation in a coordinate-invariant sense. In addition, we show how the method can be used to interpolate constrained quantities so that the given constraint is automatically respected. We showcase the method with two astrophysics related examples situated in the curved Kerr spacetime. The first problem is interpolating a locally constant vector field, in which case curvature effects are expected to be maximally important. The second example is a General Relativistic Magnetohydrodynamics simulation of a turbulent accretion flow around a black hole, wherein high intrinsic variability is expected to be at least as important as curvature effects.

categories: astro-ph.IM astro-ph.HE cs.NA math.NA physics.data-an


title: Proper Motions of Stellar Streams Discovered in the Dark Energy Survey

abstract: We cross-match high-precision astrometric data from Gaia DR2 with accurate multi-band photometry from the Dark Energy Survey (DES) DR1 to confidently measure proper motions for nine stellar streams in the DES footprint: Aliqa Uma, ATLAS, Chenab, Elqui, Indus, Jhelum, Phoenix, Tucana III, and Turranburra. We determine low-confidence proper motion measurements for four additional stellar streams: Ravi, Wambelong, Willka Yaku, and Turbio. We find evidence for a misalignment between stream tracks and the systemic proper motion of streams that may suggest a systematic gravitational influence from the Large Magellanic Cloud. These proper motions, when combined with radial velocity measurements, will allow for detailed orbit modeling which can be used to constrain properties of the LMC and its affect on nearby streams, as well as global properties of the Milky Way's gravitational potential.


title: Light composite fermions from holography

abstract: Motivated by Beyond the Standard Model theories of composite fermions or top partners, we propose a holographic mechanism that generates light baryonic states in a strongly coupled gauge theory. The starting point are the fermionic fluctuations of massive probe branes embedded into AdS$_5 \times S^5$. We first consider the D3/probe D7-brane system. We derive in detail the fermionic fluctuation equations and show the supersymmetric degeneracy of the mesinos with the mesons. Here we view the fermionic mesinos as potential realizations of composite fermions or top partners. We then add higher dimension operators and study their impact on the mesino spectrum. In particular we show that the ground state mesino mass can be pushed to an arbitrarily light value by a suitable choice of the coupling of the higher dimension operator, $g$. No matter the value of $g$, the masses of higher excited states never fall below the mass of the ground state at $g=0$. We also present similar results for the supersymmetric D3/D3 and D3/D5 systems.


title: Emergence and stability of spin-valley entangled quantum liquids in moir\'e heterostructures

abstract: Twisting moir\'e heterostructures to the flatband regime allows for the formation of strongly correlated quantum states, since the dramatic reduction of the bandwidth can cause the residual electronic interactions to set the principal energy scale. An effective description for such correlated moir\'e heterostructures, derived in the strong-coupling limit at integer filling, generically leads to spin-valley Heisenberg models. Here we explore the emergence and stability of spin liquid behavior in an SU(2)$^{\mathrm{spin}}\otimes$SU(2)$^{\mathrm{valley}}$ Heisenberg model upon inclusion of Hund's-induced and longer-ranged exchange couplings, employing a pseudofermion functional renormalization group approach. We consider two lattice geometries, triangular and honeycomb (relevant to different moir\'e heterostructures), and find, for both cases, an extended parameter regime surrounding the SU(4) symmetric point where no long-range order occurs, indicating a stable realm of quantum spin liquid behavior. For large Hund's coupling, we identify the adjacent magnetic orders, with both antiferromagnetic and ferromagnetic ground states emerging in the separate spin and valley degrees of freedom. For both lattice geometries the inclusion of longer-ranged exchange couplings is found to have both stabilizing and destabilizing effects on the spin liquid regime depending on the sign of the additional couplings.


title: QCD Resummation in Hard Diffractive Dijet Production at the Electron-Ion Collider

abstract: Diffractive dijet production at the electron-ion collider (EIC) has been proposed to study the gluon Wigner distribution at small-$x$. We investigate the soft gluon radiation associated with the final state jets and an all order resummation formula is derived. We show that the soft gluon resummation plays an important role to describe E791 data on $\pi$-induced diffractive dijet production at Fermilab. Predictions for the EIC are presented, and we emphasize that the soft gluon resummation is an important aspect to explore the nucleon/nucleus tomography through these processes.


title: Free-form Grale reconstruction of Abell 2744: robustness of uncertainties against changes in lensing data

abstract: Abell 2744, a massive Hubble Frontier Fields merging galaxy cluster with many multiple images in the core has been the subject of many lens inversions using different methods. While most existing studies compare various inversion methods, we focus on a comparison of reconstructions that use different input lensing data. Since the quantity and quality of lensing data is constantly improving, it makes sense to ask if the estimated uncertainties are robust against changes in the data. We address this question using free-form Grale, which takes only image information as input, and nothing pertaining to cluster galaxies. We reconstruct Abell 2744 using two sets of strong lensing data from the Hubble Frontier Fields community. Our first and second reconstructions use 55 and 91 images, respectively, and only 10 of the 91 images have the same positions and redshifts as in the first reconstruction. Comparison of the two mass maps shows that Grale uncertainties are robust against these changes, as well as small modifications in the inversion routine. Additionally, applying the methods used in Sebesta et al. (2016) for MACS J0416, we conclude that, in a statistical sense, light follows mass in Abell 2744, with brighter galaxies clustering stronger with the recovered mass than the fainter ones. We also show that the faintest galaxies are anti-correlated with mass, which is likely the result of light contamination from bright galaxies, and lensing magnification bias acting on galaxies background to the cluster.


title: Elementary Excitations in Fractional Quantum Hall Effect from Classical Constraints

abstract: Classical constraints on the reduced density matrix of quantum fluids in a single Landau level, termed as local exclusion conditions (LECs) [B. Yang, arXiv:1901.00047], have recently been shown to characterize the ground state of many FQH phases. In this work, we extend the LEC construction to build the elementary excitations, namely quasiholes and quasielectrons, of these FQH phases. In particular, we elucidate the quasihole counting, categorize various types of quasielectrons, and construct their microscopic wave functions. Our extensive numerical calculations indicate that the undressed quasielectron excitations of the Laughlin state obtained from LECs are topologically equivalent to those obtained from the composite fermion theory. Intriguingly, the LEC construction unveils interesting connections between different FQH phases and offers a novel perspective on exotic states such as the Gaffnian and the Fibonacci state.


title: Classification of Non-local Actions: Area versus Volume Entanglement Entropy

abstract: From the partition function for two classes of classically non-local actions containing the fractional Laplacian, we show that as long as there exists a suitable (non-local) Hilbert-space transform the underlying action can be mapped onto a purely local theory. In all such cases the partition function is equivalent to that of a local theory and an area law for the entanglement entropy obtains. When such a reduction fails, the entanglement entropy deviates strongly from an area law and can in some cases scale as the volume. As these two criteria are coincident, we conjecture that they are equivalent and provide the ultimate test for locality of Gaussian theories rather than a simple inspection of the explicit operator content.


title: IsoNN: Isomorphic Neural Network for Graph Representation Learning and Classification

abstract: Deep learning models have achieved huge success in numerous fields, such as computer vision and natural language processing. However, unlike such fields, it is hard to apply traditional deep learning models on the graph data due to the 'node-orderless' property. Normally, adjacency matrices will cast an artificial and random node-order on the graphs, which renders the performance of deep models on graph classification tasks extremely erratic, and the representations learned by such models lack clear interpretability. To eliminate the unnecessary node-order constraint, we propose a novel model named Isomorphic Neural Network (IsoNN), which learns the graph representation by extracting its isomorphic features via the graph matching between input graph and templates. IsoNN has two main components: graph isomorphic feature extraction component and classification component. The graph isomorphic feature extraction component utilizes a set of subgraph templates as the kernel variables to learn the possible subgraph patterns existing in the input graph and then computes the isomorphic features. A set of permutation matrices is used in the component to break the node-order brought by the matrix representation. Three fully-connected layers are used as the classification component in IsoNN. Extensive experiments are conducted on benchmark datasets, the experimental results can demonstrate the effectiveness of ISONN, especially compared with both classic and state-of-the-art graph classification methods.


title: Astrobites as a Community-led Model for Education, Science Communication, and Accessibility in Astrophysics

abstract: Support for early career astronomers who are just beginning to explore astronomy research is imperative to increase retention of diverse practitioners in the field. Since 2010, Astrobites has played an instrumental role in engaging members of the community -- particularly undergraduate and graduate students -- in research. In this white paper, the Astrobites collaboration outlines our multi-faceted online education platform that both eases the transition into astronomy research and promotes inclusive professional development opportunities. We additionally offer recommendations for how the astronomy community can reduce barriers to entry to astronomy research in the coming decade.


title: The BAHAMAS project: Effects of a running scalar spectral index on large-scale structure

abstract: Recent analyses of the cosmic microwave background (CMB) and the Lyman-alpha forest indicate a mild preference for a deviation from a power law primordial matter power spectrum (a so-called negative `running'). We use an extension to the BAHAMAS suite of cosmological hydrodynamic simulations to explore the effects that a running scalar spectral index has on large-scale structure (LSS), using Planck CMB constraints to initialize the simulations. We focus on 5 key statistics: i) the non-linear matter power spectrum ii) the halo mass function; iii) the halo two-point auto correlation function; iv) total mass halo density profiles; and v) the halo concentration-mass relation. In terms of the matter power spectrum, we find that a running scalar spectral index affects all k-scales examined in this study, with a negative (positive) running leading to an amplification (suppression) of power. These effects should be easily detectable with upcoming surveys such as LSST and Euclid. In the mass range sampled, a positive running leads to an increase in the mass of galaxy groups and clusters, with the favoured negative running leading to a decrease in mass of lower-mass (M <~ 10^13 M_solar) halos, but an increase for the most massive (M >~ 10^13 M_solar) halos. Changes in the mass are generally confined to 5-10% which, while not insignificant, cannot by itself reconcile the claimed tension between the primary CMB and cluster number counts. We find that running does not significantly affect the shapes of density profiles of matched halos, changing only their amplitude. Finally, we demonstrate that the observed effects on LSS due to a running scalar spectral index are separable from those of baryonic effects to typically a few percent precision.


title: Non-Standard Interactions in Radiative Neutrino Mass Models

abstract: Models of radiative Majorana neutrino masses require new scalars and/or fermions to induce lepton number violating interactions. We show that these new particles also generate observable neutrino nonstandard interactions (NSI) with matter. We classify radiative models as type-I or II, with type-I models containing at least one Standard Model (SM) particle inside the loop diagram generating neutrino mass, and type-II models having no SM particle inside the loop. While type-II radiative models do not generate tree-level NSI, popular models which fall under the type-I category are shown, somewhat surprisingly, to generate observable NSI at tree-level, while being consistent with direct and indirect constraints from colliders, electroweak precision data and charged-lepton flavor violation (cLFV). We survey such models where neutrino masses arise at one, two and three loops. In the prototypical Zee model which generates neutrino masses via one-loop diagrams involving charged scalars, we find that diagonal NSI can be as large as ($8\%, 3.8\%, 9.3\%$) for ($\varepsilon_{ee},\varepsilon_{\mu \mu}, \varepsilon_{\tau \tau}$), while off-diagonal NSI can be at most ($10^{-3}\%, 0.56\%, 0.34\%$) for ($\varepsilon_{e\mu},\varepsilon_{e \tau}, \varepsilon_{\mu \tau}$). In one-loop neutrino mass models using leptoquarks (LQs), $(\varepsilon_{\mu\mu},\, \varepsilon_{\tau\tau})$ can be as large as $(21.6\%,\,51.7\%)$, while $\varepsilon_{ee}$ and $(\varepsilon_{e\mu},\, \varepsilon_{e\tau},\,\varepsilon_{\mu\tau})$ can at most be 0.6\%. The most stringent constraints on the diagonal NSI are found to come from neutrino oscillation and scattering experiments, while the off-diagonal NSI are mostly constrained by low-energy processes, such as atomic parity violation and cLFV. While our analysis is focused on radiative neutrino mass models, it essentially covers all NSI possibilities with heavy mediators.


title: Global alignment of solution-based single-wall carbon nanotube films via machine-vision controlled filtration

abstract: Over the past decade, substantial progress has been made in the chemical processing (chiral enrichment, length sorting, handedness selectivity, and filling substance) of single-wall carbon nanotubes (SWCNTs). Recently, it was shown that large, horizontally-aligned films can be created out of post-processed SWCNT solutions. Here, we use machine-vision automation and parallelization to simultaneously produce globally-aligned SWCNT films using pressure-driven filtration. Feedback control enables filtration to occur with a constant flow rate that not only improves the nematic ordering of the SWCNT films, but also provides the ability to align a wide range of SWCNT types and on a variety of nanoporous membranes using the same filtration parameters. Using polarized optical spectroscopic techniques, we show that meniscus combing produces a two-dimensional radial SWCNT alignment on one side of the film. After we flatten the meniscus through silanation, spatially-resolved nematicity maps on both sides of the SWCNT film reveal global alignment across the entire structure. From experiments changing ionic strength and membrane tribocharging, we provide evidence that the SWCNT alignment mechanism stems from an interplay of intertube interactions and ordered membrane charging. This work opens up the possibility of creating globally-aligned SWCNT film structures for a new-generation of nanotube electronics and optical control elements.


title: deepCR: Cosmic Ray Rejection with Deep Learning

abstract: Cosmic ray (CR) identification and replacement are critical components of imaging and spectroscopic reduction pipelines involving solid-state detectors. We present deepCR, a deep learning based framework for CR identification and subsequent image inpainting based on the predicted CR mask. To demonstrate the effectiveness of this framework, we train and evaluate models on Hubble Space Telescope ACS/WFC images of sparse extragalactic fields, globular clusters, and resolved galaxies. We demonstrate that at a false positive rate of 0.5%, deepCR achieves close to 100% detection rates in both extragalactic and globular cluster fields, and 91% in resolved galaxy fields, which is a significant improvement over the current state-of-the-art method LACosmic. Compared to a multicore CPU implementation of LACosmic, deepCR CR mask predictions run up to 6.5 times faster on CPU and 90 times faster on a single GPU. For image inpainting, the mean squared errors of deepCR predictions are 20 times lower in globular cluster fields, 5 times lower in resolved galaxy fields, and 2.5 times lower in extragalactic fields, compared to the best performing non-neural technique tested. We present our framework and the trained models as an open-source Python project, with a simple-to-use API. To facilitate reproducibility of the results we also provide a benchmarking codebase.


title: Higher curvature Bianchi identities, generalised geometry and $L_{\infty}$ algebras

abstract: The Bianchi identities for bosonic fluxes in supergravity can receive higher derivative quantum and string corrections, the most well known being that of Heterotic theory $d H = \tfrac{1}{4}\alpha'(\text{tr } F^2 - \text{tr } R^2)$. Less studied are the modifications at order $R^4$ that may arise, for example, in the Bianchi identity for the seven-form flux of M theory compactifications. We argue that such corrections appear to be incompatible with the exceptional generalised geometry description of the lower order supergravity, and seem to imply a gauge algebra for the bosonic potentials that cannot be written in terms of an (exceptional) Courant bracket. However, we show that this algebra retains the form of an $L_{\infty}$ gauge field theory, which terminates at a level ten multibracket for the case involving just the seven-form flux.


title: Monomialization of a quasianalytic morphism

abstract: We prove a monomialization theorem for mappings in general classes of infinitely differentiable functions that are called quasianalytic. Examples include Denjoy-Carleman classes, the class of $\cC^\infty$ functions definable in a polynomially bounded $o$-minimal structure, as well as the classes of real- or complex analytic functions, and algebraic functions over any field of characteristic zero. The monomialization theorem asserts that a mapping in a quasianalytic class can be transformed to a mapping whose components are monomials with respect to suitable local coordinates, by sequences of simple modifications of the source and target -- local blowings-up and power substitutions in the real cases, in general, and local blowings-up alone in the algebraic or analytic cases. Monomialization is a version of resolution of singularities for a mapping. We show that it is not possible, in general, to monomialize by global blowings-up, even in the real-analytic case.


title: Pure Lovelock black hole in the dimension, $d=3N+1$, is stable

abstract: In this paper we show that pure Lovelock static Schwarzschild's analogue black hole in dimensions $d>3N+1$, where $N$ is the degree of Lovelock polynomial action, is stable even though pure Gauss-Bonnet $N=2$ black hole is unstable in dimension $d<7$. We also discuss and compare quasinormal modes for pure Lovelock and the corresponding Einstein black hole in the same dimension. We find that perturbations decay with characteristic time which is weakly dimensional dependent as it depends only on the gravitational potential of the background solution, while frequency of oscillations however depend on the dimension. Also we show that spectrum of perturbations is not isospectral except in $d=4$.


title: Reservoir Computing Models for Patient-Adaptable ECG Monitoring in Wearable Devices

abstract: The reservoir computing paradigm is employed to classify heartbeat anomalies online based on electrocardiogram signals. Inspired by the principles of information processing in the brain, reservoir computing provides a framework to design, train, and analyze recurrent neural networks (RNNs) for processing time-dependent information. Due to its computational efficiency and the fact that training amounts to a simple linear regression, this supervised learning algorithm has been variously considered as a strategy to implement useful computations not only on digital computers but also on emerging unconventional hardware platforms such as neuromorphic microchips. Here, this biological-inspired learning framework is exploited to devise an accurate patient-adaptive model that has the potential to be integrated into wearable cardiac events monitoring devices. The proposed patient-customized model was trained and tested on ECG recordings selected from the MIT-BIH arrhythmia database. Restrictive inclusion criteria were used to conduct the study only on ECGs including, at least, two classes of heartbeats with highly unequal number of instances. The results of extensive simulations showed this model not only provides accurate, cheap and fast patient-customized heartbeat classifier but also circumvents the problem of "imbalanced classes" when the readout weights are trained using weighted ridge-regression.


title: Improving Brain Magnetic Resonance Image MRI Segmentation via a Novel Algorithm based on Genetic and Regional Growth

abstract: Background: Regarding the importance of right diagnosis in medical applications, various methods have been exploited for processing medical images solar. The method of segmentation is used to analyze anal to miscall structures in medical imaging. Objective: This study describes a new method for brain Magnetic Resonance Image (MRI) segmentation via a novel algorithm based on genetic and regional growth. Methods: Among medical imaging methods, brains MRI segmentation is important due to the high contrast of non-intrusive soft tissue and high spatial resolution. Size variations of brain tissues are often accompanied by various diseases such as Alzheimers disease. As our knowledge about the relationship between various brain diseases and deviation of brain anatomy increases, MRI segmentation is exploited as the first step in early diagnosis. In this paper, the regional growth method and auto-mate selection of initial points by genetic algorithm are used to introduce a new method for MRI segmentation. Primary pixels and similarity criterion are automatically by genetic algorithms to maximize the accuracy and validity in image segmentation. Results: By using genetic algorithms and defining the fixed function of image segmentation, the initial points for the algorithm were found. The proposed algorithms are applied to the images and results are manually selected by regional growth in which the initial points were compared. The results showed that the proposed algorithm could reduce segmentation error effectively. Conclusion: The study concluded that the proposed algorithm could reduce segmentation error effectively and help us to diagnose brain diseases.


title: The X-ray-Resolved Supernova Remnant S8 in the Dwarf Irregular Galaxy IC 1613

abstract: We conducted an observation of the nearby irregular galaxy IC 1613 with the Chandra X-ray Observatory using the S3 chip of the ACIS with an effective exposure time of 49.9 ksec. The observation primarily targeted the extensive bubble and star formation region in the northeast quadrant of the galaxy. The only known supernova remnant (SNR) in IC 1613, S8, is also the galaxy's most luminous X-ray source (L_X (0.3-8 keV) ~5-6x10^{36} erg/s). We resolve the SNR with Chandra and compare its nearly circular X-ray morphology with Halpha and radio images. We assign an upper limit on the luminosity of any possible associated compact central object of ~4x10^{35}$ erg/s (0.3-8 keV band) and conclude that we would detect a Crab-like pulsar but not a Cas A-like object. We infer an age for S8 of ~3400-5600 years and compare it to other SNRs in the Local Group. We suggest that S8 is a young composite SNR.


title: Robust and optimal sparse regression for nonlinear PDE models

abstract: This paper investigates how models of spatiotemporal dynamics in the form of nonlinear partial differential equations can be identified directly from noisy data using a combination of sparse regression and weak formulation. Using the 4th-order Kuramoto-Sivashinsky equation for illustration, we show how this approach can be optimized in the limits of low and high noise, achieving accuracy that is orders of magnitude better than what existing techniques allow. In particular, we derive the scaling relation between the accuracy of the model, the parameters of the weak formulation, and the properties of the data, such as its spatial and temporal resolution and the level of noise.


title: High-Contrast Testbeds for Future Space-Based Direct Imaging Exoplanet Missions

abstract: Instrumentation techniques in the field of direct imaging of exoplanets have greatly advanced over the last two decades. Two of the four NASA-commissioned large concept studies involve a high-contrast instrument for the imaging and spectral characterization of exo-Earths from space: LUVOIR and HabEx. This whitepaper describes the status of 8 optical testbeds in the US and France currently in operation to experimentally validate the necessary technologies to image exo-Earths from space. They explore two complementary axes of research: (i) coronagraph designs and manufacturing and (ii) active wavefront correction methods and technologies. Several instrument architectures are currently being analyzed in parallel to provide more degrees of freedom for designing the future coronagraphic instruments. The necessary level of performance has already been demonstrated in-laboratory for clear off-axis telescopes (HabEx-like) and important efforts are currently in development to reproduce this accomplishment on segmented and/or on-axis telescopes (LUVOIR-like) over the next two years.


title: Some Results on Tighter Bayesian Lower Bounds on the Mean-Square Error

abstract: In random parameter estimation, Bayesian lower bounds (BLBs) for the mean-square error have been noticed to not be tight in a number of cases, even when the sample size, or the signal-to-noise ratio, grow to infinity. In this paper, we study alternative forms of BLBs obtained from a covariance inequality, where the inner product is based on the \textit{a posteriori} instead of the joint probability density function. We hence obtain a family of BLBs, which is shown to form a counterpart at least as tight as the well-known Weiss-Weinstein family of BLBs, and we extend it to the general case of vector parameter estimation. Conditions for equality between these two families are provided. Focusing on the Bayesian Cram\'er-Rao bound (BCRB), a definition of efficiency is proposed relatively to its tighter form, and efficient estimators are described for various types of common estimation problems, e.g., scalar, exponential family model parameter estimation. Finally, an example is provided, for which the classical BCRB is known to not be tight, while we show its tighter form is, based on formal proofs of asymptotic efficiency of Bayesian estimators. This analysis is finally corroborated by numerical results.


title: Symplectic fillings of asymptotically dynamically convex manifolds I

abstract: We consider exact fillings with vanishing first Chern class of asymptotically dynamically convex (ADC) manifolds. We construct two structure maps on the positive symplectic cohomology and prove that they are independent of the filling for ADC manifolds. The invariance of the structure maps implies that the vanishing of symplectic cohomology and the existence of symplectic dilations are properties independent of the filling for ADC manifolds. Using them, various topological applications on symplectic fillings are obtained, including the uniqueness of diffeomorphism types of fillings for many contact manifolds. We use the structure maps to define the first symplectic obstructions to Weinstein fillability. In particular, we show that for all dimension $4k+3, k\ge 1$, there exist infinitely many contact manifolds that are exactly fillable, almost Weinstein fillable but not Weinstein fillable. The invariance of the structure maps generalizes to strong fillings with vanishing first Chern class. We show that any strong filling with vanishing first Chern class of a class of manifolds, including $(S^{2n-1},\xi_{std}), \partial(T^*L \times \mathbb{C}^n)$ with $L$ simply connected, must be exact and have unique diffeomorphism type. As an application of the proof, we show that the existence of symplectic dilation implies uniruledness. In particular any affine exotic $\mathbb{C}^n$ with non-negative log Kodaira dimension is a symplectic exotic $\mathbb{C}^{n}$.


title: Universal Person Re-Identification

abstract: Most state-of-the-art person re-identification (re-id) methods depend on supervised model learning with a large set of cross-view identity labelled training data. Even worse, such trained models are limited to only the same-domain deployment with significantly degraded cross-domain generalization capability, i.e. "domain specific". To solve this limitation, there are a number of recent unsupervised domain adaptation and unsupervised learning methods that leverage unlabelled target domain training data. However, these methods need to train a separate model for each target domain as supervised learning methods. This conventional "{\em train once, run once}" pattern is unscalable to a large number of target domains typically encountered in real-world deployments. We address this problem by presenting a "train once, run everywhere" pattern industry-scale systems are desperate for. We formulate a "universal model learning' approach enabling domain-generic person re-id using only limited training data of a "{\em single}" seed domain. Specifically, we train a universal re-id deep model to discriminate between a set of transformed person identity classes. Each of such classes is formed by applying a variety of random appearance transformations to the images of that class, where the transformations simulate the camera viewing conditions of any domains for making the model training domain generic. Extensive evaluations show the superiority of our method for universal person re-id over a wide variety of state-of-the-art unsupervised domain adaptation and unsupervised learning re-id methods on five standard benchmarks: Market-1501, DukeMTMC, CUHK03, MSMT17, and VIPeR.


title: The effect of valley, spin and band nesting on the electronic properties of gated quantum dots in a single layer of transition metal dichalcogenides (TMDCs)

abstract: We present here results of atomistic theory of electrons confined by metallic gates in a single layer of transition metal dichalcogenides. The electronic states are described by the tight-binding model and computed using a computational box including up to million atoms with periodic boundary conditions and parabolic confining potential due to external gates embedded in it. With this methodology applied to MoS2, we find a twofold degenerate energy spectrum of electrons confined in the two non-equivalent K-valleys by the metallic gates as well as six-fold degenerate spectrum associated with Q-valleys. We compare the electron spectrum with the energy levels of electrons confined in GaAs/GaAlAs and in self-assembled quantum dots. We discuss the role of spin splitting and topological moments on the K and Q valley electronic states in quantum dots with sizes comparable to experiment.


title: Taming Koepke's Zoo II: Register Machines

abstract: We study the computational strength of resetting $\alpha$-register machines, a model of transfinite computability introduced by P. Koepke in \cite{K1}. Specifically, we prove the following strengthening of a result from \cite{C}: For an exponentially closed ordinal $\alpha$, we have $L_{\alpha}\models$ZF$^{-}$ if and only if COMP$^{\text{ITRM}}_{\alpha}=L_{\alpha+1}\cap\mathfrak{P}(\alpha)$, i.e. if and only if the set of $\alpha$-ITRM-computable subsets of $\alpha$ coincides with the set of subsets of $\alpha$ in $L_{\alpha+1}$. Moreover, we show that, if $\alpha$ is exponentially closed and $L_{\alpha}\not\models$ZF$^{-}$, then COMP$^{\text{ITRM}}_{\alpha}=L_{\beta(\alpha)}\cap\mathfrak{P}(\alpha)$, where $\beta(\alpha)$ is the supremum of the $\alpha$-ITRM-clockable ordinals, which coincides with the supremum of the $\alpha$-ITRM-computable ordinals. We also determine the set of subsets of $\alpha$ computable by an $\alpha$-ITRM with time bounded below $\delta$ when $\delta>\alpha$ is an exponentially closed ordinal smaller than the supremum of the $\alpha$-ITRM-clockable ordinals. Moreover, we obtain some sufficient and necessary conditions on ordinals $\alpha$ for which the $\alpha$-wITRM-clockable ordinals are bounded by $\alpha$.


title: Unsupervised Learning Technique to Obtain the Coordinates of Wi-Fi Access Points

abstract: Given that the accuracy of range-based positioning techniques generally increases with the number of available anchor nodes, it is important to secure more of these nodes. To this end, this paper studies an unsupervised learning technique to obtain the coordinates of unknown nodes that coexist with anchor nodes. As users use the location services in an area of interests, the proposed method automatically discovers unknown nodes and estimates their coordinates. In addition, this method learns an appropriate calibration curve to correct the distortion of raw distance measurements. As such, the positioning accuracy can be greatly improved using more anchor nodes and well-calibrated distance measurements. The performance of the proposed method was verified using commercial Wi-Fi devices in a practical indoor environment. The experiment results show that the coordinates of unknown nodes and the calibration curve are simultaneously determined without any ground truth data.


title: Vertical electrolyte transistor operating at very low voltage and high current density

abstract: In this work it is reported a vertical electrolyte transistor (VET) whose structure is based on stacked layers as described below: bottom contact (source or drain) - channel - permeable intermediate electrode (drain or source) - ion gel (electrolyte gate dielectric) - gate top contact. This VET depicts versatility to work as Electrolyte-Gated Vertical Organic Field Effect Transistor (Electrolyte-Gated VOFET) or Vertical Organic Electrochemical Transistor (VOECT) as never reported before. The distinction of these operation modes is regarding to the transistor transconductance that occurs due to induced charge carriers or ionic current, respectively. Both modes of operation show that this VET is able to work at very low voltage range and drive a high current density. These observed features make VETs a good candidate for applications in iontronic devices, bio-sensors and/or very low power optoelectronic circuits.


title: Distributed Quantum Fiber Magnetometry

abstract: Nitrogen-vacancy (NV) quantum magnetometers offer exceptional sensitivity and long-term stability. However, their use to date in distributed sensing applications, including remote detection of ferrous metals, geophysics, and biosensing, has been limited due to the need to combine optical, RF, and magnetic excitations into a single system. Existing approaches have yielded localized devices but not distributed capabilities. In this study, we report on a continuous system-in-a-fiber architecture that enables distributed magnetic sensing over extended lengths. Key to this realization is a thermally drawn fiber that has hundreds of embedded photodiodes connected in parallel and a hollow optical waveguide that contains a fluid with NV diamonds. This fiber is placed in a larger coaxial cable to deliver the required RF excitation. We realize this distributed quantum sensor in a water-immersible 90-meter-long cable with 102 detection sites. A sensitivity of 63 nT Hz-1/2 per site, limited by laser shot noise, was established along a 90 m test section. This fiber architecture opens new possibilities as a robust and scalable platform for distributed quantum sensing technologies.

categories: physics.app-ph physics.ins-det physics.optics quant-ph


title: Ultra-high-speed color imaging with single-pixel detectors under low light level

abstract: Single-pixel imaging is suitable for low light level scenarios because a bucket detector is employed to maximally collect the light from an object. However, one of the challenges is its slow imaging speed, mainly due to the slow light modulation technique. We here demonstrate 1.4MHz video imaging based on computational ghost imaging with a RGB LED array having a full-range frame rate up to 100MHz. With this method, the motion of a high speed propeller is observed. Moreover, by exploiting single-photon detectors to increase the detection efficiency, this method is developed for ultra-high-speed imaging under low light level.


title: ${\cal N}{=}4$ Supersymmetric $d=1$ Sigma Models on Group Manifolds

abstract: We construct manifestly ${\cal N}=4$ supersymmetric off-shell superfield actions for the HKT $d=1$ sigma models on the group manifolds U(2) and SU(3), using the harmonic $d=1$ approach. The underlying $({\bf 4, 4, 0})$ and $({\bf 4, 4, 0})\oplus({\bf 4, 4, 0})$ multiplets are described, respectively, by one and two harmonic analytic superfields $q^+$ satisfying the appropriate nonlinear harmonic constraints. The invariant actions in both cases are bilinear in the superfields. We present the corresponding superfield realizations of the U(2) and SU(3) isometries and show that in fact they are enlarged to the products U(2)$\times$SU(2) and SU(3)$\times$U(2). We prove the corresponding invariances at both the superfield and component levels and present the bosonic $d=1$ sigma model actions, as integral over $t$ in the U(2) case and over $t$ and ${\rm SU}(2)$ harmonics in the SU(3) case. In the U(2) case we also give a detailed comparison with the general harmonic approach to HKT models and establish a correspondence with a particular action of the off-shell nonlinear multiplet $({\bf 3, 4, 1})$. A possible way of generalizing U(2) model to the matrix U($2n$) case is suggested.


title: Resonant inelastic x-ray scattering of magnetic excitations under pressure

abstract: Resonant inelastic x-ray scattering (RIXS) is an extremely valuable tool for the study of elementary, including magnetic, excitations in matter. Latest developments of this technique mostly aimed at improving the energy resolution and performing polarization analysis of the scattered radiation, with a great impact on the interpretation and applicability of RIXS. Instead, this article focuses on the sample environment and presents a setup for high-pressure low-temperature RIXS measurements of low-energy excitations. The feasibility of these experiments is proved by probing the magnetic excitations of the bilayer iridate Sr$_3$Ir$_2$O$_7$ at pressures up to 12 GPa.


title: Vadere: An open-source simulation framework to promote interdisciplinary understanding

abstract: Pedestrian dynamics is an interdisciplinary field of research. Psychologists, sociologists, traffic engineers, physicists, mathematicians and computer scientists all strive to understand the dynamics of a moving crowd. In principle, computer simulations offer means to further this understanding. Yet, unlike for many classic dynamical systems in physics, there is no universally accepted locomotion model for crowd dynamics. On the contrary, a multitude of approaches, with very different characteristics, compete. Often only the experts in one special model type are able to assess the consequences these characteristics have on a simulation study. Therefore, scientists from all disciplines who wish to use simulations to analyze pedestrian dynamics need a tool to compare competing approaches. Developers, too, would profit from an easy way to get insight into an alternative modeling ansatz. Vadere meets this interdisciplinary demand by offering an open-source simulation framework that is lightweight in its approach and in its user interface while offering pre-implemented versions of the most widely spread models.


title: On the link between the Maxwell and linearised Einstein equations on Schwarzschild

abstract: In this short note we shall demonstrate that given a smooth solution $\gamma$ to the linearised Einstein equations on Schwarzschild which is supported on the $l\geq 2$ spherical harmonics and expressed relative to a transverse and traceless gauge then one can construct from it a smooth solution to the sourced Maxwell equations expressed relative to a generalised Lorentz gauge. Here the Maxwell current is constructed from those gauge-invariant combinations of the components of $\gamma$ which are determined by solutions to the Regge--Wheeler and Zerilli equations. The result thus provides an elegant link between the spin 1 and spin 2 equations on Schwarzschild.


title: Factor Analysis for High-Dimensional Time Series with Change Point

abstract: We consider change-point latent factor models for high-dimensional time series, where a structural break may exist in the underlying factor structure. In particular, we propose consistent estimators for factor loading spaces before and after the change point, and the problem of estimating the change-point location is also considered. Compared with existing results on change-point factor analysis of high-dimensional time series, a distinguished feature of the current paper is that our results allow strong cross-sectional dependence in the noise process. To accommodate the unknown degree of cross-sectional dependence strength, we propose to use self-normalization to pivotalize the change-point test statistic. Numerical experiments including a Monte Carlo simulation study and a real data application are presented to illustrate the proposed methods.


title: An end-to-end (deep) neural network applied to raw EEG, fNIRs and body motion data for data fusion and BCI classification task without any pre-/post-processing

abstract: Brain computer interfaces (BCI) using EEG, fNIRS and body motion (MoCap) data are getting more attention due to the fact that fNIRS and MoCap are not prone to movement artifacts similar to other brain imaging techniques such as EEG. Advancements in deep learning (neural networks) would allow the use of raw data for efficient feature extraction without any pre-/post-processing. In this work, we are performing human activity recognition (BCI classification task) for 5 activity classes using an end-to-end (deep) neural network (NN) (from input all the way to the output) on raw fNIRS, EEG and MoCap data. Our core contribution is focused on applying an end-to-end NN model without any pre-/post-processing on the data. The entire NN model is being trained using backpropagation algorithm. Our end-to-end model is composed of a four-layered MLP: input layer, two hidden layers (using fully connected (dense) layer, batch normalization and leaky-RELU as non-linearity and activation function), and output layer using softmax. We have reached minimum 90\% accuracy on the test dataset for the classification task on 10 subjects data and 5 classes of activity.


title: Development and Applicability of Online Passivity Enforced Wide-Band Multi-Port Equivalents For Hybrid Transient Simulation

abstract: This paper presents a method for developing single and multi-port frequency dependent network equivalent (FDNE) based on a passivity enforced online recursive least squares identification algorithm, which identifies the input admittance matrix in z-domain. Furthermore, with the proposed architecture, a real-time hybrid model of the reduced power system is developed that integrate transient stability analysis and FDNE. Main advantages of the proposed architecture are, it identifies the FDNE even with unknown network parameters in the frequency range of interest, and yet can be implemented directly due to discrete formulation while maintaining desired accuracy, stability, and passivity conditions. The accuracy and characteristics of the proposed method are verified by implementing on two-area, IEEE 39 and 68 bus power system models.


title: Interaction between water and carbon nanostructures: How good are current density functional approximations?

abstract: Due to their current and future technological applications, including realisation of water filters and desalination membranes, water adsorption on graphitic sp$^{2}$-bonded carbon is of overwhelming interest. However, these systems are notoriously challenging to model, even for electronic structure methods such as density functional theory (DFT), because of the crucial role played by London dispersion forces and non-covalent interactions in general. Recent efforts have established reference quality interactions of several carbon nanostructures interacting with water. Here, we compile a new benchmark set (dubbed \textbf{WaC18}), which includes a single water molecule interacting with a broad range of carbon structures, and various bulk (3D) and two dimensional (2D) ice polymorphs. The performance of 28 approaches, including semi-local exchange-correlation functionals, non-local (Fock) exchange contributions, and long-range van der Waals (vdW) treatments, are tested by computing the deviations from the reference interaction energies. The calculated mean absolute deviations on the WaC18 set depends crucially on the DFT approach, ranging from 135 meV for LDA to 12 meV for PBE0-D4. We find that modern vdW corrections to DFT significantly improve over their precursors. Within the 28 tested approaches, we identify the best performing within the functional classes of: generalized gradient approximated (GGA), meta-GGA, vdW-DF, and hybrid DF, which are BLYP-D4, TPSS-D4, rev-vdW-DF2, and PBE0-D4, respectively.


title: Flattened stellar systems based on distribution functions depending on actions

abstract: We address an issue that arises when self-consistently flattened dynamical stellar systems are constructed by adopting a distribution function (DF) that depends on the action integrals. The velocity distribution at points on the symmetry axis is controlled by the the dependence of the DF on just one action, while at points off the symmetry axis two actions are involved. Consequently, the physical requirement that the velocity distribution evolves continuously in the neighbourhood of the symmetry axis restricts the functional forms of acceptable DFs. An algorithm for conforming to this restriction is presented and used to construct a variety of flattened models.


title: Maximizing Stylistic Control and Semantic Accuracy in NLG: Personality Variation and Discourse Contrast

abstract: Neural generation methods for task-oriented dialogue typically generate from a meaning representation that is populated using a database of domain information, such as a table of data describing a restaurant. While earlier work focused solely on the semantic fidelity of outputs, recent work has started to explore methods for controlling the style of the generated text while simultaneously achieving semantic accuracy. Here we experiment with two stylistic benchmark tasks, generating language that exhibits variation in personality, and generating discourse contrast. We report a huge performance improvement in both stylistic control and semantic accuracy over the state of the art on both of these benchmarks. We test several different models and show that putting stylistic conditioning in the decoder and eliminating the semantic re-ranker used in earlier models results in more than 15 points higher BLEU for Personality, with a reduction of semantic error to near zero. We also report an improvement from .75 to .81 in controlling contrast and a reduction in semantic error from 16% to 2%.


title: Topological and subsystem codes on low-degree graphs with flag qubits

abstract: In this work we introduce two code families, which we call the heavy hexagon code and heavy square code. Both code families are implemented by assigning physical data and ancilla qubits to both vertices and edges of low degree graphs. Such a layout is particularly suitable for superconducting qubit architectures to minimize frequency collisions and crosstalk. In some cases, frequency collisions can be reduced by several orders of magnitude. The heavy hexagon code is a hybrid surface/Bacon-Shor code mapped onto a (heavy) hexagonal lattice whereas the heavy square code is the surface code mapped onto a (heavy) square lattice. In both cases, the lattice includes all the ancilla qubits required for fault-tolerant error-correction. Naively, the limited qubit connectivity might be thought to limit the error-correcting capability of the code to less than its full distance. Therefore, essential to our construction is the use of flag qubits. We modify minimum weight perfect matching decoding to efficiently and scalably incorporate information from measurements of the flag qubits and correct up to the full code distance while respecting the limited connectivity. Simulations show that high threshold values for both codes can be obtained using our decoding protocol. Further, our decoding scheme can be adapted to other topological code families.


title: Path-accelerated molecular dynamics: Parallel-in-time integration using path integrals

abstract: Massively parallel computer architectures create new opportunities for the performance of long-timescale molecular dynamics (MD) simulations. Here, we introduce the path-accelerated molecular dynamics (PAMD) method that takes advantage of distributed computing to reduce the wall-clock time of MD simulation via parallelization with respect to MD timesteps. The marginal distribution for the time evolution of a system is expressed in terms of a path integral, enabling the use of path sampling techniques to numerically integrate MD trajectories. By parallelizing the evaluation of the path action with respect to time and by initializing the path configurations from a non-equilibrium distribution, the algorithm enables significant speedups in terms of the length of MD trajectories that can be integrated in a given amount of wall-clock time. The method is demonstrated for Brownian dynamics, although it is generalizable to other stochastic equations of motion including open systems. We apply the method to two simple systems, a harmonic oscillator and a Lennard-Jones liquid, and we show that in comparison to the conventional Euler integration scheme for Brownian dynamics, the new method can reduce the wall-clock time for integrating trajectories of a given length by more than three orders of magnitude in the former system and more than two in the latter. This new method for parallelizing MD in the dimension of time can be trivially combined with algorithms for parallelizing the MD force evaluation to achieve further speedup.


title: Random Hamiltonians with Arbitrary Point Interactions

abstract: We consider disordered Hamiltonians given by the Laplace operator subject to arbitrary random self-adjoint singular perturbations supported on random discrete subsets of the real line. Under minimal assumptions on the type of disorder, we prove the following dichotomy: Either every realization of the random operator has purely absolutely continuous spectrum or spectral and exponential dynamical localization hold. In particular, we establish Anderson localization for Schr\"odinger operators with Bernoulli-type random singular potential and singular density.


title: Between the deterministic and non-deterministic query complexity

abstract: We consider problems that can be solved by asking certain queries. The deterministic query complexity $D(P)$ of a problem $P$ is the smallest number of queries needed to ask in order to find the solution (in the worst case), while the non-deterministic query complexity $D_0(P)$ is the smallest number of queries needed to ask, in case we know the solution, to prove that it is indeed the solution (in the worst case). Equivalently, $D(P)$ is the largest number of queries needed to find the solution in case an Adversary is answering the queries, while $D_0(P)$ is the largest number of queries needed to find the solution in case an Adversary chooses the input. We define a series of quantities between these two values, $D_k(P)$ is the largest number of queries needed to find the solution in case an Adversary chooses the input, and answers the queries, but he can change the input at most $k$ times. We give bounds on $D_k(P)$ for various problems $P$.


title: Computational p-Willmore Flow with Conformal Penalty

abstract: The unsigned p-Willmore functional introduced in \cite{mondino2011} generalizes important geometric functionals which measure the area and Willmore energy of immersed surfaces. Presently, techniques from \cite{dziuk2008} are adapted to compute the first variation of this functional as a weak-form system of equations, which are subsequently used to develop a model for the p-Willmore flow of closed surfaces in $\mathbb{R}^3$. This model is amenable to constraints on surface area and enclosed volume, and is shown to decrease the p-Willmore energy monotonically over time. In addition, a penalty-based regularization procedure is formulated to prevent artificial mesh degeneration along the flow; inspired by a conformality condition derived in \cite{kamberov1996}, this procedure encourages angle-preservation in a closed and oriented surface immersion as it evolves. Following this, a finite-element discretization of both systems is discussed, and an application to mesh editing is presented.