title: Unimodality-Constrained Matrix Factorization for Non-Parametric Source Localization

abstract: Herein, the problem of simultaneous localization of multiple sources given a number of energy samples at different locations is examined. The strategies do not require knowledge of the signal propagation models, nor do they exploit the spatial signatures of the source. A non-parametric source localization framework based on a matrix observation model is developed. It is shown that the source location can be estimated by localizing the peaks of a pair of location signature vectors extracted from the incomplete energy observation matrix. A robust peak localization algorithm is developed and shown to decrease the source localization mean squared error (MSE) faster than O(1/M^1.5) with M samples, when there is no measurement noise. To extract the source signature vectors from a matrix with mixed energy from multiple sources, a unimodality-constrained matrix factorization (UMF) problem is formulated, and two rotation techniques are developed to solve the UMF efficiently. Our numerical experiments demonstrate that the proposed scheme achieves similar performance as the kernel regression baseline using only 1/5 energy measurement samples in detecting a single source, and the performance gain is more significant in the cases of detecting multiple sources.


title: Comparison of Methods for modelling Coronal Magnetic Fields

abstract: Four different approximate approaches used to model the stressing of coronal magnetic fields due to an imposed photospheric motion are compared with each other and the results from a full time-dependent magnetohydrodynamic (MHD) code. The assumptions used for each of the approximate methods are tested by considering large photospheric footpoint displacements. We consider a simple model problem, comparing the full nonlinear magnetohydrodynamic evolution, determined with the Lare2D numerical code, with four approximate approaches. Two of these, magneto-frictional relaxation and a quasi-1D Grad-Shafranov approach, assume sequences of equilibria, whilst the other two methods, a second-order linearisation of the MHD equations and Reduced MHD, are time-dependent. The relaxation method is very accurate compared to full MHD for force-free equilibria for all footpoint displacements but has significant errors when the plasma $\beta_0$ is of order unity. The 1D approach gives an extremely accurate description of the equilibria away from the photospheric boundary layers, and agrees well with Lare2D for all parameter values tested. The linearised MHD equations correctly predict the existence of photospheric boundary layers that are present in the full MHD results. As soon as the footpoint displacement becomes a significant fraction of the loop length, the RMHD method fails to model the sequences of equilibria correctly. The full numerical solution is interesting in its own right, and care must be taken for low $\beta_0$ plasmas if the viscosity is too large.


title: SquishedNets: Squishing SqueezeNet further for edge device scenarios via deep evolutionary synthesis

abstract: While deep neural networks have been shown in recent years to outperform other machine learning methods in a wide range of applications, one of the biggest challenges with enabling deep neural networks for widespread deployment on edge devices such as mobile and other consumer devices is high computational and memory requirements. Recently, there has been greater exploration into small deep neural network architectures that are more suitable for edge devices, with one of the most popular architectures being SqueezeNet, with an incredibly small model size of 4.8MB. Taking further advantage of the notion that many applications of machine learning on edge devices are often characterized by a low number of target classes, this study explores the utility of combining architectural modifications and an evolutionary synthesis strategy for synthesizing even smaller deep neural architectures based on the more recent SqueezeNet v1.1 macroarchitecture for applications with fewer target classes. In particular, architectural modifications are first made to SqueezeNet v1.1 to accommodate for a 10-class ImageNet-10 dataset, and then an evolutionary synthesis strategy is leveraged to synthesize more efficient deep neural networks based on this modified macroarchitecture. The resulting SquishedNets possess model sizes ranging from 2.4MB to 0.95MB (~5.17X smaller than SqueezeNet v1.1, or 253X smaller than AlexNet). Furthermore, the SquishedNets are still able to achieve accuracies ranging from 81.2% to 77%, and able to process at speeds of 156 images/sec to as much as 256 images/sec on a Nvidia Jetson TX1 embedded chip. These preliminary results show that a combination of architectural modifications and an evolutionary synthesis strategy can be a useful tool for producing very small deep neural network architectures that are well-suited for edge device scenarios.


title: Regular 1-harmonic flow

abstract: We consider the 1-harmonic flow of maps from a bounded domain into a submanifold of a Euclidean space, i.e. the gradient flow of the total variation functional restricted to maps taking values in the manifold. We restrict ourselves to Lipschitz initial data. We prove uniqueness and, in the case of a convex domain, local existence of solutions to the flow equations. If the target manifold has non-positive sectional curvature or in the case that the datum is small, solutions are shown to exist globally and to become constant in finite time. We also consider the case where the domain is a compact Riemannian manifold without boundary, solving the homotopy problem for 1-harmonic maps under some assumptions.


title: Bidirectional Conditional Generative Adversarial Networks

abstract: Conditional Generative Adversarial Networks (cGANs) are generative models that can produce data samples ($x$) conditioned on both latent variables ($z$) and known auxiliary information ($c$). We propose the Bidirectional cGAN (BiCoGAN), which effectively disentangles $z$ and $c$ in the generation process and provides an encoder that learns inverse mappings from $x$ to both $z$ and $c$, trained jointly with the generator and the discriminator. We present crucial techniques for training BiCoGANs, which involve an extrinsic factor loss along with an associated dynamically-tuned importance weight. As compared to other encoder-based cGANs, BiCoGANs encode $c$ more accurately, and utilize $z$ and $c$ more effectively and in a more disentangled way to generate samples.


title: Real-time brain machine interaction via social robot gesture control

abstract: Brain-Machine Interaction (BMI) system motivates interesting and promising results in forward/feedback control consistent with human intention. It holds great promise for advancements in patient care and applications to neurorehabilitation. Here, we propose a novel neurofeedback-based BCI robotic platform using a personalized social robot in order to assist patients having cognitive deficits through bilateral rehabilitation and mental training. For initial testing of the platform, electroencephalography (EEG) brainwaves of a human user were collected in real time during tasks of imaginary movements. First, the brainwaves associated with imagined body kinematics parameters were decoded to control a cursor on a computer screen in training protocol. Then, the experienced subject was able to interact with a social robot via our real-time BMI robotic platform. Corresponding to subject's imagery performance, he/she received specific gesture movements and eye color changes as neural-based feedback from the robot. This hands-free neurofeedback interaction not only can be used for mind control of a social robot's movements, but also sets the stage for application to enhancing and recovering mental abilities such as attention via training in humans by providing real-time neurofeedback from a social robot.


title: Quantum simulation of the spin-boson model with a microwave circuit

abstract: We consider superconducting circuits for the purpose of simulating the spin-boson model. The spin-boson model consists of a single two-level system coupled to bosonic modes. In most cases, the model is considered in a limit where the bosonic modes are sufficiently dense to form a continuous spectral bath. A very well known case is the ohmic bath, where the density of states grows linearly with the frequency. In the limit of weak coupling or large temperature, this problem can be solved numerically. If the coupling is strong, the bosonic modes can become sufficiently excited to make a classical simulation impossible. Here, we discuss how a quantum simulation of this problem can be performed by coupling a superconducting qubit to a set of microwave resonators. We demonstrate a possible implementation of a continuous spectral bath with individual bath resonators coupling strongly to the qubit. Applying a microwave drive scheme potentially allows us to access the strong-coupling regime of the spin-boson model. We discuss how the resulting spin relaxation dynamics with different initialization conditions can be probed by standard qubit-readout techniques from circuit quantum electrodynamics.


title: Measuring the Sterile Neutrino CP Phase at DUNE and T2HK

abstract: The CP phases associated with the sterile neutrino cannot be measured in the dedicated short-baseline experiments being built to test the sterile neutrino hypothesis. On the other hand, these phases can be measured in long-baseline experiments, even though the main goal of these experiments is not to test or measure sterile neutrino parameters. In particular, the sterile neutrino phase $\delta_{24}$ affects the charged-current electron appearance data in long-baseline experiment. In this paper we show for the first time how well the sterile neutrino phase $\delta_{24}$ can be measured by the next-generation long-baseline experiments DUNE, T2HK (and T2HKK). We also show the expected precision with which this sterile phase can be measured by combining the DUNE data with data from T2HK or T2HKK. We also present the sensitivity of these experiments to the sterile mixing angles, both by themselves, as well as when DUNE is combined with T2HK or T2HKK.


title: Better Agnostic Clustering Via Relaxed Tensor Norms

abstract: We develop a new family of convex relaxations for $k$-means clustering based on sum-of-squares norms, a relaxation of the injective tensor norm that is efficiently computable using the Sum-of-Squares algorithm. We give an algorithm based on this relaxation that recovers a faithful approximation to the true means in the given data whenever the low-degree moments of the points in each cluster have bounded sum-of-squares norms. We then prove a sharp upper bound on the sum-of-squares norms for moment tensors of any distribution that satisfies the \emph{Poincare inequality}. The Poincare inequality is a central inequality in probability theory, and a large class of distributions satisfy it including Gaussians, product distributions, strongly log-concave distributions, and any sum or uniformly continuous transformation of such distributions. As an immediate corollary, for any $\gamma > 0$, we obtain an efficient algorithm for learning the means of a mixture of $k$ arbitrary \Poincare distributions in $\mathbb{R}^d$ in time $d^{O(1/\gamma)}$ so long as the means have separation $\Omega(k^{\gamma})$. This in particular yields an algorithm for learning Gaussian mixtures with separation $\Omega(k^{\gamma})$, thus partially resolving an open problem of Regev and Vijayaraghavan \citet{regev2017learning}. Our algorithm works even in the outlier-robust setting where an $\epsilon$ fraction of arbitrary outliers are added to the data, as long as the fraction of outliers is smaller than the smallest cluster. We, therefore, obtain results in the strong agnostic setting where, in addition to not knowing the distribution family, the data itself may be arbitrarily corrupted.


title: Lensing convergence in galaxy clustering in LambdaCDM and beyond

abstract: We study the impact of neglecting lensing magnification in galaxy clustering analyses for future galaxy surveys, considering the LambdaCDM model and two extensions: massive neutrinos and modifications of General Relativity. Our study focuses on the biases on the constraints and on the estimation of the cosmological parameters. We perform a comprehensive investigation of these two effects for the upcoming photometric and spectroscopic galaxy surveys Euclid and SKA for different redshift binning configurations. We also provide a fitting formula for the magnification bias of SKA. Our results show that the information present in the lensing contribution does improve the constraints on the modified gravity parameters whereas the lensing constraining power is negligible for the LambdaCDM parameters. For photometric surveys the estimation is biased for all the parameters if lensing is not taken into account. This effect is particularly significant for the modified gravity parameters. Conversely for spectroscopic surveys the bias is below one sigma for all the parameters. Our findings show the importance of including lensing in galaxy clustering analyses for testing General Relativity and to constrain the parameters which describe its modifications.


title: Complete super-sample lensing covariance in the response approach

abstract: We derive the complete super-sample covariance (SSC) of the matter and weak lensing convergence power spectra using the power spectrum response formalism to accurately describe the coupling of super- to sub-survey modes. The SSC term is completely characterized by the survey window function, the nonlinear matter power spectrum and the full first-order nonlinear power spectrum response function, which describes the response to super-survey density and tidal field perturbations. Generalized separate universe simulations can efficiently measure these responses in the nonlinear regime of structure formation, which is necessary for lensing applications. We derive the lensing SSC formulae for two cases: one under the Limber and flat-sky approximations, and a more general one that goes beyond the Limber approximation in the super-survey mode and is valid for curved sky applications. Quantitatively, we find that for sky fractions $f_{\rm sky} \approx 0.3$ and a single source redshift at $z_S=1$, the use of the flat-sky and Limber approximation underestimates the total SSC contribution by $\approx 10\%$. The contribution from super-survey tidal fields to the lensing SSC, which has not been included in cosmological analyses so far, is shown to represent about $5\%$ of the total lensing covariance on multipoles $\ell_1,\ell_2 \gtrsim 300$. The SSC is the dominant off-diagonal contribution to the total lensing covariance, making it appropriate to include these tidal terms and beyond flat-sky/Limber corrections in cosmic shear analyses.


title: Glitch Classification and Clustering for LIGO with Deep Transfer Learning

abstract: The detection of gravitational waves with LIGO and Virgo requires a detailed understanding of the response of these instruments in the presence of environmental and instrumental noise. Of particular interest is the study of anomalous non-Gaussian noise transients known as glitches, since their high occurrence rate in LIGO/Virgo data can obscure or even mimic true gravitational wave signals. Therefore, successfully identifying and excising glitches is of utmost importance to detect and characterize gravitational waves. In this article, we present the first application of Deep Learning combined with Transfer Learning for glitch classification, using real data from LIGO's first discovery campaign labeled by Gravity Spy, showing that knowledge from pre-trained models for real-world object recognition can be transferred for classifying spectrograms of glitches. We demonstrate that this method enables the optimal use of very deep convolutional neural networks for glitch classification given small unbalanced training datasets, significantly reduces the training time, and achieves state-of-the-art accuracy above 98.8%. Once trained via transfer learning, we show that the networks can be truncated and used as feature extractors for unsupervised clustering to automatically group together new classes of glitches and anomalies. This novel capability is of critical importance to identify and remove new types of glitches which will occur as the LIGO/Virgo detectors gradually attain design sensitivity.


title: Mapping the Tidal Disruption of the Hercules Dwarf: A Wide-Field DECam Imaging Search for RR Lyrae Stars

abstract: We investigate the hypothesized tidal disruption of the Hercules ultra-faint dwarf galaxy (UFD). Previous tidal disruption studies of the Hercules UFD have been hindered by the high degree of foreground contamination in the direction of the dwarf. We bypass this issue by using RR Lyrae stars, which are standard candles with a very low field-volume density at the distance of Hercules. We use wide-field imaging from the Dark Energy Camera on CTIO to identify candidate RR Lyrae stars, supplemented with observations taken in coordination with the Beijing-Arizona Sky Survey on the Bok Telescope. Combining color, magnitude, and light-curve information, we identify three new RR Lyrae stars associated with Hercules. All three of these new RR Lyrae stars lie outside its published tidal radius. When considered with the nine RR Lyrae stars already known within the tidal radius, these results suggest that a substantial fraction of Hercules' stellar content has been stripped. With this degree of tidal disruption, Hercules is an interesting case between a visibly disrupted dwarf (such as the Sagittarius dwarf spheroidal galaxy) and one in dynamic equilibrium. The degree of disruption also shows that we must be more careful with the ways we determine object membership when estimating dwarf masses in the future. One of the three discovered RR Lyrae stars sits along the minor axis of Hercules, but over two tidal radii away. This type of debris is consistent with recent models that suggest Hercules' orbit is aligned with its minor axis.


title: Theory and Practice of Electron Diffraction from Single Atoms and Extended Objects using an Electron Microscope Pixel Array Detector

abstract: What does the diffraction pattern from a single atom look like? How does it differ from the scattering from long range potential? With the development of new high-dynamic range pixel array detectors to measure the complete momentum distribution, these questions have immediate relevance for designing and understanding momentum-resolved imaging modes. We explore the asymptotic limits of long range and short range potentials. We use a simple quantum mechanical model to explain the general and asymptotic limits for the probability distribution in both and real and reciprocal space. Features in the scattering potential much larger than the probe size cause the bright-field disk to deflect uniformly, while features much smaller than the probe size, instead of a deflection cause a redistribution of intensity within the bright-field disk. Because long range and short range features are encoded differently in the diffraction pattern, it is possible to separate their contributions in differential phase contrast (DPC) or Center-of-Mass (CoM) imaging. The shape profiles for atomic resolution CoM imaging are dominated by the shape of the probe gradient and not the highly-singular atomic potentials or their local fields. Instead, only the peak height shows an atomic-number sensitivity, whose precise dependence is determined by the convergence angle. At lower convergence angles, the contrast oscillates with increasing atomic number, similar to bright field imaging. The range of collection angles impacts DPC and CoM imaging differently, with CoM being more sensitive to the upper cutoff limit, while DPC is more sensitive to the lower cutoff.


title: Global current circuit structure in a resistive pulsar magnetosphere model

abstract: Pulsar magnetospheres have strong magnetic fields and large amounts of plasma. The structures of these magnetospheres are studied using force-free electrodynamics. To understand pulsar magnetospheres, discussions must include their outer region. However, force-free electrodynamics is limited in it does not handle dissipation. Therefore, a resistive pulsar magnetic field model is needed. To break the ideal magnetohydrodynamic (MHD) condition $E \cdot B = 0$, Ohm's law is used. In this work, I introduce resistivity depending upon the distance from the star and obtain a self-consistent steady state by time integration. Poloidal current circuits form in the magnetosphere while the toroidal magnetic-field region expands beyond the light cylinder and the Poynting flux radiation appears. High electric resistivity causes a large space scale poloidal current circuit and the magnetosphere radiates a larger Poynting flux than the linear increase outside of the light cylinder radius. The formed poloidal current circuit has width, which grows with the electric conductivity. This result contributes to a more concrete dissipative pulsar magnetosphere model.


title: KIC 8462852: Potential repeat of the Kepler day 1540 dip in August 2017

abstract: We report 33 V-band observations by the Hereford Arizona Observatory (HAO) of the enigmatic star KIC 8462852 during the two week period 3-17 August 2017. We find a striking resemblance of these observations to the Kepler day 1540 dip with HAO observations tracking the Kepler light curve (adjusted for egress symmetry). A possible explanation of this potential repeat transit is a brown dwarf and extensive ring system in a 1601-day eccentric orbit. We suggest this object may be detectable through radial velocity observations in October and November 2017, with an amplitude of ~ 1-2 kms-1.


title: Determination of carbohydrates in infant milk powders by ultra-performance liquid chromatography with evaporative light scattering detector and BEH HILIC column

abstract: The possibility of carbohydrate separation in BEH HILIC (Ethylene Bridged Hybride, Hydrophilic Interaction Liquid Chromatography) column was studied by ultra-performance liquid chromatography (UPLC) with evaporative light scattering detector (ELSD) and mobile phase containing amine compounds as modifiers. The chromatography conditions and ELSD parameters were optimized to separate five typical carbohydrates and applied to analysis of four infant milk powders. The linear ranges of carbohydrate determination were 20-300mg/L for fructose and glucose, 20-250mg/L for sucrose and lactose, and 35-180mg/L for fructo-oligosaccharide. The LODs were 16.4mg/L for fructose and glucose, 17.3mg/L for sucrose, 20.0mg/L for lactose, and 46.7mg/L for fructo-oligosaccharide. Relative standard deviations (RSDs) ranged between 3.45-4.23%, 1.46-4.17%, 4.14-5.60%, 1.39-4.09%, and 2.49-3.61% for fructose, glucose, sucrose, lactose, and fructo-oilgosaccharide, respectively and recoveries ranged between 95.0 and 105.4%


title: XSAT of Linear CNF Formulas

abstract: Open questions with respect to the computational complexity of linear CNF formulas in connection with regularity and uniformity are addressed. In particular it is proven that any l-regular monotone CNF formula is XSAT-unsatisfiable if its number of clauses m is not a multiple of l. For exact linear formulas one finds surprisingly that l-regularity implies k-uniformity, with m = 1 + k(l-1)) and allowed k-values obey k(k-1) = 0 (mod l). Then the computational complexity of the class of monotone exact linear and l-regular CNF formulas with respect to XSAT can be determined: XSAT-satisfiability is either trivial, if m is not a multiple of l, or it can be decided in sub-exponential time, namely O(exp(n^^1/2)). Sub-exponential time behaviour for the wider class of regular and uniform linear CNF formulas can be shown for certain subclasses.


title: Analytical evaluation of relativistic molecular integrals. I. Auxiliary functions

abstract: The auxiliary functions provide efficient computation of integrals arising at the self-consistent field (SCF) level for molecules using Slater-type bases. This applies both in relativistic and non-relativistic electronic structure theory. The relativistic molecular auxiliary functions derived in our previous paper [Phys. Rev. E 91, 023303 (2015)] are discussed here in detail. Two solution methods are proposed in the present study. The ill-conditioned binomial series representation formulae first, are replaced by convergent series representation for incomplete beta functions then, they are improved by inserting an extra parameter used to extend the domain of convergence. Highly accurate results can be achieved for integrals by the procedures discussed in the present study which also places no restrictions on quantum numbers in all ranges of orbital parameters. The difficulty of obtaining analytical relations associated with using non-integer Slater-type orbitals which are non-analytic in the sense of complex analysis at r=0 is therefore, eliminated.


title: Virtual Adversarial Ladder Networks For Semi-supervised Learning

abstract: Semi-supervised learning (SSL) partially circumvents the high cost of labeling data by augmenting a small labeled dataset with a large and relatively cheap unlabeled dataset drawn from the same distribution. This paper offers a novel interpretation of two deep learning-based SSL approaches, ladder networks and virtual adversarial training (VAT), as applying distributional smoothing to their respective latent spaces. We propose a class of models that fuse these approaches. We achieve near-supervised accuracy with high consistency on the MNIST dataset using just 5 labels per class: our best model, ladder with layer-wise virtual adversarial noise (LVAN-LW), achieves 1.42% +/- 0.12 average error rate on the MNIST test set, in comparison with 1.62% +/- 0.65 reported for the ladder network. On adversarial examples generated with L2-normalized fast gradient method, LVAN-LW trained with 5 examples per class achieves average error rate 2.4% +/- 0.3 compared to 68.6% +/- 6.5 for the ladder network and 9.9% +/- 7.5 for VAT.


title: MaMaDroid: Detecting Android Malware by Building Markov Chains of Behavioral Models (Extended Version)

abstract: As Android has become increasingly popular, so has malware targeting it, thus pushing the research community to propose different detection techniques. However, the constant evolution of the Android ecosystem, and of malware itself, makes it hard to design robust tools that can operate for long periods of time without the need for modifications or costly re-training. Aiming to address this issue, we set to detect malware from a behavioral point of view, modeled as the sequence of abstracted API calls. We introduce MaMaDroid, a static-analysis based system that abstracts the API calls performed by an app to their class, package, or family, and builds a model from their sequences obtained from the call graph of an app as Markov chains. This ensures that the model is more resilient to API changes and the features set is of manageable size. We evaluate MaMaDroid using a dataset of 8.5K benign and 35.5K malicious apps collected over a period of six years, showing that it effectively detects malware (with up to 0.99 F-measure) and keeps its detection capabilities for long periods of time (up to 0.87 F-measure two years after training). We also show that MaMaDroid remarkably outperforms DroidAPIMiner, a state-of-the-art detection system that relies on the frequency of (raw) API calls. Aiming to assess whether MaMaDroid's effectiveness mainly stems from the API abstraction or from the sequencing modeling, we also evaluate a variant of it that uses frequency (instead of sequences), of abstracted API calls. We find that it is not as accurate, failing to capture maliciousness when trained on malware samples that include API calls that are equally or more frequently used by benign apps.


title: Implementing the Deep Q-Network

abstract: The Deep Q-Network proposed by Mnih et al. [2015] has become a benchmark and building point for much deep reinforcement learning research. However, replicating results for complex systems is often challenging since original scientific publications are not always able to describe in detail every important parameter setting and software engineering solution. In this paper, we present results from our work reproducing the results of the DQN paper. We highlight key areas in the implementation that were not covered in great detail in the original paper to make it easier for researchers to replicate these results, including termination conditions and gradient descent algorithms. Finally, we discuss methods for improving the computational performance and provide our own implementation that is designed to work with a range of domains, and not just the original Arcade Learning Environment [Bellemare et al., 2013].


title: Teaching a Machine to Read Maps with Deep Reinforcement Learning

abstract: The ability to use a 2D map to navigate a complex 3D environment is quite remarkable, and even difficult for many humans. Localization and navigation is also an important problem in domains such as robotics, and has recently become a focus of the deep reinforcement learning community. In this paper we teach a reinforcement learning agent to read a map in order to find the shortest way out of a random maze it has never seen before. Our system combines several state-of-the-art methods such as A3C and incorporates novel elements such as a recurrent localization cell. Our agent learns to localize itself based on 3D first person images and an approximate orientation angle. The agent generalizes well to bigger mazes, showing that it learned useful localization and navigation capabilities.


title: E-PUR: An Energy-Efficient Processing Unit for Recurrent Neural Networks

abstract: Recurrent Neural Networks (RNNs) are a key technology for emerging applications such as automatic speech recognition, machine translation or image description. Long Short Term Memory (LSTM) networks are the most successful RNN implementation, as they can learn long term dependencies to achieve high accuracy. Unfortunately, the recurrent nature of LSTM networks significantly constrains the amount of parallelism and, hence, multicore CPUs and many-core GPUs exhibit poor efficiency for RNN inference. In this paper, we present E-PUR, an energy-efficient processing unit tailored to the requirements of LSTM computation. The main goal of E-PUR is to support large recurrent neural networks for low-power mobile devices. E-PUR provides an efficient hardware implementation of LSTM networks that is flexible to support diverse applications. One of its main novelties is a technique that we call Maximizing Weight Locality (MWL), which improves the temporal locality of the memory accesses for fetching the synaptic weights, reducing the memory requirements by a large extent. Our experimental results show that E-PUR achieves real-time performance for different LSTM networks, while reducing energy consumption by orders of magnitude with respect to general-purpose processors and GPUs, and it requires a very small chip area. Compared to a modern mobile SoC, an NVIDIA Tegra X1, E-PUR provides an average energy reduction of 92x.


title: Simulating and interpreting deep observations in the Hubble Ultra Deep Field with the JWST/NIRSpec low-resolution 'prism'

abstract: The James Webb Space Telescope (JWST) will enable the detection of optical emission lines in galaxies spanning a broad range of luminosities out to redshifts z>10. Measurements of key galaxy properties, such as star formation rate and metallicity, through these observations will provide unique insight into, e.g., the role of feedback from stars and active galactic nuclei (AGNs) in regulating galaxy evolution, the co-evolution of AGNs and host galaxies, the physical origin of the 'main sequence' of star-forming galaxies and the contribution by star-forming galaxies to cosmic reionization. We present an original framework to simulate and analyse observations performed with the Near Infrared Spectrograph (NIRSpec) on board JWST. We use the BEAGLE tool (BayEsian Analysis of GaLaxy sEds) to build a semi-empirical catalogue of galaxy spectra based on photometric spectral energy distributions (SEDs) of dropout galaxies in the Hubble Ultra Deep Field (HUDF). We demonstrate that the resulting catalogue of galaxy spectra satisfies different types of observational constraints on high redshift galaxies, and use it as input to simulate NIRSpec/prism (R~100) observations. We show that a single 'deep' (~100 ks) NIRSpec/prism pointing in the HUDF will enable S/N>3 detections of multiple optical emission lines in ~30 (~60) galaxies at z>6 (z~4-6) down to m_F160W<30 AB mag. Such observations will allow measurements of galaxy star formation rates, ionization parameters and gas-phase metallicities within factors of 1.5, mass-to-light ratios within a factor of 2, galaxy ages within a factor of 3 and V-band attenuation optical depths with a precision of 0.3.


title: Measuring the Local Diffusion Coefficient with H.E.S.S. Observations of Very High-Energy Electrons

abstract: The HAWC Collaboration has recently reported the detection of bright and spatially extended multi-TeV gamma-ray emission from Geminga, Monogem, and a handful of other nearby, middle-aged pulsars. The angular profile of the emission observed from these pulsars is surprising, in that it implies that cosmic-ray diffusion is significantly inhibited within ~25 pc of these objects, compared to the expectations of standard Galactic diffusion models. This raises the important question of whether the diffusion coefficient in the local interstellar medium is also low, or whether it is instead better fit by the mean Galactic value. Here, we utilize recent observations of the cosmic-ray electron spectrum (extending up to ~20 TeV) by the H.E.S.S. Collaboration to show that the local diffusion coefficient cannot be as low as it is in the regions surrounding Geminga and Monogem. Instead, we conclude that cosmic rays efficiently diffuse through the bulk of the local interstellar medium. Among other implications, this further supports the conclusion that pulsars significantly contribute to the observed positron excess.


title: Thermalized axion inflation: natural and monomial inflation with small $r$

abstract: A safe way to reheat the universe, in models of natural and quadratic inflation, is through shift symmetric couplings between the inflaton $\phi$ and the Standard Model (SM), since they do not generate loop corrections to the potential $V(\phi)$. We consider such a coupling to SM gauge fields, of the form $\phi F\tilde{F}/f$, with sub-Planckian $f$. In this case gauge fields can be exponentially produced already {\it during inflation} and thermalize via interactions with charged particles, as pointed out in previous work. This can lead to a plasma of temperature $T$ during inflation and the thermal masses $gT$ of the gauge bosons can equilibrate the system. In addition, inflaton perturbations $\delta \phi$ can also have a thermal spectrum if they have sufficiently large cross sections with the plasma. In this case inflationary predictions are strongly modified: (1) scalar perturbations are thermal, and so enhanced over the vacuum, leading to a generic way to {\it suppress} the tensor-to-scalar ratio $r$; (2) the spectral index is $n_s-1=\eta-4\epsilon$. After presenting the relevant conditions for thermalization, we show that thermalized natural and monomial models of inflation agree with present observations and have $r\approx 10^{-3} - 10^{-2}$, which is within reach of next generation CMB experiments.


title: Consistent Searches for SMEFT Effects in Non-Resonant Dijet Events

abstract: We investigate the bounds which can be placed on generic new-physics contributions to dijet production at the LHC using the framework of the Standard Model Effective Field Theory, deriving the first consistently-treated EFT bounds from non-resonant high-energy data. We recast an analysis searching for quark compositeness, equivalent to treating the SM with one higher-dimensional operator as a complete UV model. In order to reach consistent, model-independent EFT conclusions, it is necessary to truncate the EFT effects consistently at order $1/\Lambda^2$ and to include the possibility of multiple operators simultaneously contributing to the observables, neither of which has been done in previous searches of this nature. Furthermore, it is important to give consistent error estimates for the theoretical predictions of the signal model, particularly in the region of phase space where the probed energy is approaching the cutoff scale of the EFT. There are two linear combinations of operators which contribute to dijet production in the SMEFT with distinct angular behavior; we identify those linear combinations and determine the ability of LHC searches to constrain them simultaneously. Consistently treating the EFT generically leads to weakened bounds on new-physics parameters. These constraints will be a useful input to future global analyses in the SMEFT framework, and the techniques used here to consistently search for EFT effects are directly applicable to other off-resonance signals.


title: Unmixing the Galactic Halo with RR Lyrae tagging

abstract: We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously.


title: Quantum equivalence of $f(R)$ gravity and scalar-tensor theories

abstract: We investigate whether the classical equivalence of $f(R)$ gravity and its formulation as scalar-tensor theory still holds at the quantum level. We explicitly compare the corresponding one-loop divergences and find that the equivalence is broken by off-shell quantum corrections, but recovered on-shell.


title: The emergence of the galactic stellar mass function from a non-universal IMF in clusters

abstract: We investigate how a single generation galactic mass function (SGMF) depends on the existence of variations in the initial stellar mass functions (IMF) of stellar clusters. We show that cluster-to-cluster variations of the IMF lead to a multicomponent SGMF where each component in a given mass range can be described by a distinct power-law function. We also show that a dispersion of $\approx 0.3$ M$_{\odot}$ in the characteristic mass of the IMF, as observed for young Galactic clusters, leads to a low mass slope of the SGMF that matches the observed Galactic stellar mass function even when the IMFs in the low mass end of individual clusters are much steeper.


title: Uncovering the identities of compact objects in high-mass X-ray binaries and gamma-ray binaries by astrometric measurements

abstract: We develop a method for identifying a compact object in binary systems with astrometric measurements and apply it to some binaries. Compact objects in some high-mass X-ray binaries and gamma-ray binaries are unknown, which is responsible for the fact that emission mechanisms in such systems have not yet confirmed. The accurate estimate of the mass of the compact object allows us to identify the compact object in such systems. Astrometric measurements are expected to enable us to estimate the masses of the compact objects in the binary systems via a determination of a binary orbit. We aim to evaluate the possibility of the identification of the compact objects for some binary systems. We then calculate probabilities that the compact object is correctly identified with astrometric observation (= confidence level) by taking into account a dependence of the orbital shape on orbital parameters and distributions of masses of white dwarfs, neutron stars, and black holes. We find that the astrometric measurements with the precision of 70 $\mu$as for $\gamma$~Cas allows us to identify the compact object at 99% confidence level if the compact object is a white dwarf with 0.6 M$_{\odot}$. In addition, we can identify the compact object with the precision of 10 $\mu$as at 97% or larger confidence level for LS~I~+61$^{\circ}$~303 and 99% or larger for HESS~J0632+057. These results imply that the astrometric measurements with the 10-$\mu$as precision level can realize the identification of compact objects for $\gamma$~Cas, LS~I~+61$^{\circ}$~303, and HESS~J0632+057.


title: ALMA and VLA Observations of the HD 141569 System

abstract: We present VLA 9 mm (33 GHz) observations of the HD 141569 system from semester 16A. The observations achieve a resolution of 0.25 arcsec ($\sim28$ au) and a sensitivity of $4.7~\mu \rm Jy~beam^{-1}$. We find (1) a $52\pm 5~\mu$Jy point source at the location of HD 141569A that shows potential variability, (2) the detected flux is contained within the SED-inferred central clearing of the disc meaning the spectral index of the dust disc is steeper than previously inferred, and (3) the M dwarf companions are also detected and variable. Previous lower-resolution VLA observations (semester 14A) found a higher flux density, interpreted as solely dust emission. When combined with ALMA observations, the VLA 14A observations suggested the spectral index and grain size distribution of HD 141569's disc was shallow and an outlier among debris systems. Using archival ALMA observations of HD 141569 at 0.87 mm and 2.9 mm we find a dust spectral index of $\alpha_{\rm mm} = 1.81\pm 0.20$. The VLA 16A flux corresponds to a brightness temperature of $\sim5\times10^{6}$ K, suggesting strong non-disc emission is affecting the inferred grain properties. The VLA 16A flux density of the M2V companion HD 141569B is $149\pm9~\mu$Jy, corresponding to a brightness temperature of $\sim2\times10^{8}$ K and suggesting significant stellar variability when compared to the VLA14A observations, which are smaller by a factor of $\sim6$.


title: Spacetimes dressed with stealth electromagnetic fields

abstract: Stealth field configurations by definition have a vanishing energy-momentum tensor, thus do not contribute to the gravitational field equations. While only trivial fields can be stealth in Maxwell's electrodynamics, nontrivial stealth fields appear in some nonlinear models of electromagnetism. We find the necessary and sufficient conditions for the electromagnetic fields to be stealth and analyse which models admit such configurations. Furthermore, we present some concrete exact solutions, featuring a class of black holes dressed with the stealth electromagnetic hair, closely related to force-free solutions. Stealth hair does not alter the generalized Smarr formula, but may contribute to the Komar charges.


title: Neutrino Oscillation Measurements Computed in Quantum Field Theory

abstract: We perform a calculation in quantum field theory of neutrino oscillation probabilities, where we include simultaneously the source, detector, and neutrino fields in the Hamiltonian. Within the appropriate limits associated with current neutrino oscillation experiments, we recover the standard oscillation formula. On the other hand, we find that the dominant contributions to the amplitude are associated with different neutrino mass eigenstates being emitted at different times, such that they arrive at the detector at the same time. This is contrary to the neutrino wave packet picture, where they are emitted simultaneously and separate as they travel to the detector. This has direct consequences regarding the mechanisms that lead to a damping of neutrino oscillations for very long baselines. Our analysis also provides a pedagogical example of a measurement process in quantum mechanics.


title: Quasar lenses and pairs in the VST-ATLAS and Gaia

abstract: We report on discovery results from a quasar lens search in the ATLAS public footprint, extending quasar lens searches to a regime without $u-$band or fiber-spectroscopic information, using a combination of data mining techniques on multi-band catalog magnitudes and image-cutout modelling. Spectroscopic follow-up campaigns, conducted at the 2.6m Nordic Optical Telescope (La Palma) and 3.6m New Technology Telescope (La Silla) in 2016, yielded seven pairs of quasars exhibiting the same lines at the same redshift and monotonic flux-ratios with wavelength (hereafter NIQs, Nearly Identical Quasar pairs). The quasar redshifts range between $\approx1.2$ and $\approx 2.7;$ contaminants are typically pairs of bright blue stars, quasar-star alignments along the line of sight, and narrow-line galaxies at $0.3<z<0.7.$ Magellan data of A0140-1152 (01$^h$40$^m$03.0$^s$-11$^d$52$^m$19.0$^s$, $z_{s}=1.807$) confirm it as a lens with deflector at $z_{l}=0.277$ and Einstein radius $\theta_{\rm E}=(0.73\pm0.02)^\ase$. We show the use of spatial resolution from the Gaia mission to select lenses and list additional systems from a WISE-Gaia-ATLAS search, yielding three additional lenses (02$^h$35$^m$27.4$^s$-24$^d$33$^m$13.2$^s$, 02$^h$59$^m$33.$^s$-23$^d$38$^m$01.8$^s$, 01$^h$46$^m$32.9$^s$-11$^d$33$^m$39.0$^s$). The overall sample consists of 11 lenses/NIQs, plus three lenses known before 2016, over the ATLAS-DR3 footprint ($\approx3500$~deg$^2$). Finally, we discuss future prospects for objective classification of pair/NIQ/contaminant spectra.


title: A Magellanic origin for the Virgo substructure

abstract: Iorio et al. (2018) mapped out the Milky Way halo using a sample of RR Lyrae stars drawn from a cross-match of Gaia with 2MASS. We investigate the significant residual in their model which we constrain to lie at Galactocentric radii $12<R<27\;\mathrm{kpc}$ and extend over $2600\;\mathrm{deg}^2$ of the sky. A counterpart of this structure exists in both the Catalina Real Time Survey and the sample of RR Lyrae variables identified in Pan-STARRS by Hernitschek et al. (2016), demonstrating that this structure is not caused by the spatial inhomogeneity of Gaia. The structure is likely the Virgo Stellar Stream and/or Virgo Over-Density. We show the structure is aligned with the Magellanic Stream and suggest that it is either debris from a disrupted dwarf galaxy that was a member of the Vast Polar Structure or that it is SMC debris from a tidal interaction of the SMC and LMC $3\;\mathrm{Gyr}$ ago. If the latter then the sub-structure in Virgo may have a Magellanic origin.


title: Light Resonances and the Low-$q^2$ Bin of $R_{K^*}$

abstract: LHCb has reported hints of lepton-flavor universality violation in the rare decays $B \to K^{(*)} \ell^+\ell^-$, both in high- and low-$q^2$ bins. Although the high-$q^2$ hint may be explained by new short-ranged interactions, the low-$q^2$ one cannot. We thus explore the possibility that the latter is explained by a new light resonance. We find that LHCb's central value of $R_{K^*}$ in the low-$q^2$ bin is achievable in a restricted parameter space of new-physics scenarios in which the new, light resonance decays preferentially to electrons and has a mass within approximately $10$ MeV of the di-muon threshold. Interestingly, such an explanation can have a kinematic origin and does not require a source of lepton-flavor universality violation. A model-independent prediction is a narrow peak in the differential $B \to K^* e^+e^-$ rate close to the di-muon threshold. If such a peak is observed, other observables, such as the differential $B \to K e^+e^-$ rate and $R_K$, may be employed to distinguish between models. However, if a low-mass resonance is not observed and the low-$q^2$ anomaly increases in significance, then the case for an experimental origin of the lepton-flavor universality violating anomalies would be strengthened. To further explore this, we also point out that, in analogy to $J/\psi$ decays, $e^+e^-$ and $\mu^+\mu^-$ decays of $\phi$ mesons can be used as a cross check of lepton-flavor universality by LHCb with $5$ fb$^{-1}$ of integrated luminosity.


title: Detection of HC$_5$N and HC$_7$N Isotopologues in TMC-1 with the Green Bank Telescope

abstract: We report the first interstellar detection of DC$_7$N and six $^{13}$C-bearing isotopologues of HC$_7$N toward the dark cloud TMC-1 through observations with the Green Bank Telescope, and confirm the recent detection of HC$_5$$^{15}$N. For the average of the $^{13}$C isotopomers, DC$_7$N, and HC$_5$$^{15}$N, we derive column densities of 1.9(2)$\times$10$^{11}$, 2.5(9)$\times$10$^{11}$, and 1.5(4)$\times$10$^{11}$ cm$^{-2}$, respectively. The resulting isotopic ratios are consistent with previous values derived from similar species in the source, and we discuss the implications for the formation chemistry of the observed cyanopolyynes. Within our uncertainties, no significant $^{13}$C isotopomer variation is found for HC$_7$N, limiting the significance CN could have in its production. The results further show that, for all observed isotopes, HC$_5$N may be isotopically depleted relative to HC$_3$N and HC$_7$N, suggesting that reactions starting from smaller cyanopolyynes may not be efficient to form HC$_{n}$N. This leads to the conclusion that the dominant production route may be the reaction between hydrocarbon ions and nitrogen atoms.


title: The globular cluster NGC7492 and the Sagittarius tidal stream: together but unmixed

abstract: We have derived from VIMOS spectroscopy the radial velocities for a sample of 71 stars selected from CFHT/Megacam photometry around the Galactic globular cluster NGC7492. In the resulting velocity distribution, it is possible to distinguish two relevant non-Galactic kinematic components along the same line of sight: a group of stars at $\langle{v_{\rm r}}\rangle \sim 125$km s$^{-1}$ which is compatible with the velocity of the old leading arm of the Sagittarius tidal stream, and a larger number of objects at $\langle{v_{\rm r}}\rangle \sim -110$km s$^{-1}$ that might be identified as members of the trailing wrap of the same stream. The systemic velocity of NGC7492 set at $v_{\rm r} \sim -177$km s$^{-1}$ differs significantly from that of both components, thus our results confirm that this cluster is not one of the globular clusters deposited by the Sagittarius dwarf spheroidal in the Galactic halo, even if it is immersed in the stream. A group of stars with $<v_{\rm r}> \sim -180$km s$^{-1}$ might be comprised of cluster members along one of the tidal tails of NGC7492.


title: Lense-Thirring precession in ULXs as a possible means to constrain the neutron star equation-of-state

abstract: The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in the disc (implying surface dipole fields that are sub-magnetar in strength), then the large scale-height flow can precess through the Lense-Thirring effect which could provide an explanation for the observed super-orbital periods. By connecting the details of the Lense-Thirring effect with the observed pulsar spin period, we are able to infer the moment-of-inertia and therefore equation-of-state of the neutron star without relying on the inclination of, or distance to the system. We apply our technique to the case of NGC 7793 P13 and demonstrate that stronger magnetic fields imply stiffer equations of state. We discuss the caveats and uncertainties, many of which can be addressed through forthcoming radiative magnetohydrodynamic (RMHD) simulations and their connection to observation.


title: Topological Defects in a Deformed Gauge Theory

abstract: In this paper, we will analyse the topological defects in a deformation of a non-abelian gauge theory using the Polyakov variables. The gauge theory will be deformed by the existence of a minimum measurable length scale in the background spacetime. We will construct the Polyakov loops for this deformed non-abelian gauge theory, and use these deformed loop space variables for obtaining a deformed loop space curvature. It will be demonstrated that this curvature will vanish if the deformed Bianchi identities are satisfied. However, it is possible that the original Bianchi identities are satisfied, but the deformed Bianchi identities are violated at the leading order in the deformation parameter, due to some topological defects. Thus, topological defects could be produced purely from a deformation of the background geometry.


title: Elliptic Genus Derivation of 4d Holomorphic Blocks

abstract: We study elliptic vortices on $\mathbb{C}\times T^2$ by considering the 2d quiver gauge theory describing their moduli spaces. The elliptic genus of these moduli spaces is the elliptic version of vortex partition function of the 4d theory. We focus on two examples: the first is a $\mathcal{N}=1$, $\mathrm{U}(N)$ gauge theory with fundamental and anti-fundamental matter; the second is a $\mathcal{N}=2$, $\mathrm{U}(N)$ gauge theory with matter in the fundamental representation. The results are instances of 4d "holomorphic blocks" into which partition functions on more complicated surfaces factorize. They can also be interpreted as free-field representations of elliptic Virasoro algebrae.


title: Robust entanglement renormalization on a noisy quantum computer

abstract: A method to study strongly interacting quantum many-body systems at and away from criticality is proposed. The method is based on a MERA-like tensor network that can be efficiently and reliably contracted on a noisy quantum computer using a number of qubits that is much smaller than the system size. We prove that the outcome of the contraction is stable to noise and that the estimated energy upper bounds the ground state energy. The stability, which we numerically substantiate, follows from the positivity of operator scaling dimensions under renormalization group flow. The variational upper bound follows from a particular assignment of physical qubits to different locations of the tensor network plus the assumption that the noise model is local. We postulate a scaling law for how well the tensor network can approximate ground states of lattice regulated conformal field theories in d spatial dimensions and provide evidence for the postulate. Under this postulate, a $O(\log^{d}(1/\delta))$-qubit quantum computer can prepare a valid quantum-mechanical state with energy density $\delta$ above the ground state. In the presence of noise, $\delta = O(\epsilon \log^{d+1}(1/\epsilon))$ can be achieved, where $\epsilon$ is the noise strength.


title: The Combinatorics of Higher Derivatives of Implicit Functions

abstract: We prove a closed formula for the derivative, of any order, of a implicit function, in terms of some binomial building blocks, and explain the combinatorics behind the coefficients appearing in the formula.


title: Light scalars on cosmological backgrounds

abstract: We study the behaviour of a light quartically self-interacting scalar field $\phi$ on curved backgrounds that may be described with the cosmological equation state parameter $w$. At leading order in the non-perturbative 2PI expansion we find a general formula for the variance $\langle\hat{\phi}^2\rangle$ and show for several previously unexplored cases, including matter domination and kination, that the curvature of space can induce a significant excitation of the field. We discuss how the generation of a non-zero variance for $w\neq-1$ can be understood as a process of self-regulation of the infrared divergences very similarly to what is known to occur in de Sitter space. To conclude, the appearance of an effective mass due to self-interaction is generic for a light scalar in curved space and can have important implications for reheating, vacuum stability and dark matter generation.


title: Determining the Elemental and Isotopic Composition of the preSolar Nebula from Genesis Data Analysis: The Case of Oxygen

abstract: We compare element and isotopic fractionations measured in solar wind samples collected by NASA's Genesis mission with those predicted from models incorporating both the ponderomotive force in the chromosphere and conservation of the first adiabatic invariant in the low corona. Generally good agreement is found, suggesting that these factors are consistent with the process of solar wind fractionation. Based on bulk wind measurements, we also consider in more detail the isotopic and elemental abundances of O. We find mild support for an O abundance in the range 8.75 - 8.83, with a value as low as 8.69 disfavored. A stronger conclusion must await solar wind regime specific measurements from the Genesis samples.


title: The dynamical matter density in the solar neighbourhood inferred from Gaia DR1

abstract: We determine the total dynamical density in the solar neighbourhood using the Tycho-Gaia Astrometric Solution (TGAS) catalogue. Astrometric measurements of proper motion and parallax of stars inform us of both the stellar number density distribution and the velocity distribution of stars close to the plane. Assuming equilibrium, these distributions are interrelated through the local dynamical density. For the first time, we do a full joint fit of the velocity and stellar number density distributions while accounting for the astrometric error of individual stars, in the framework of Bayesian Hierarchical Modelling. We use a sample of stars whose distance extends to approximately 160 pc from the Sun. We find a local matter density of $\rho_0 = 0.119^{+0.015}_{-0.012}~M_\odot$pc$^{-3}$, where the result is presented as the median to the posterior distribution, plus/minus its 16th and 84th percentiles. We find the Sun's position above the Galactic plane to be $z_\odot = 15.3^{+2.24}_{-2.16}$ pc, and the Sun's velocity perpendicular to the Galactic plane to be $W_\odot = 7.19^{+0.18}_{-0.18}$ km/s.


title: Light-cone spreading of perturbations and the butterfly effect in a classical spin chain

abstract: We find that localised perturbations in a chaotic classical many-body system-- the classical Heisenberg We find that the effects of a localised perturbation in a chaotic classical many-body system--the classical Heisenberg chain at infinite temperature--spread ballistically with a finite speed even when the local spin dynamics is diffusive. We study two complementary aspects of this butterfly effect: the rapid growth of the perturbation, and its simultaneous ballistic (light-cone) spread, as characterised by the Lyapunov exponents and the butterfly speed respectively. We connect this to recent studies of the out-of-time-ordered commutators (OTOC), which have been proposed as an indicator of chaos in a quantum system. We provide a straightforward identification of the OTOC with a natural correlator in our system and demonstrate that many of its interesting qualitative features are present in the classical system. Finally, by analysing the scaling forms, we relate the growth, spread and propagation of the perturbation with the growth of one-dimensional interfaces described by the Kardar-Parisi-Zhang (KPZ) equation.


title: A matrix analysis approach to discrete comparison principles for nonmonotone PDE

abstract: We consider a linear algebra approach to establishing a discrete comparison principle for a nonmonotone class of quasilinear elliptic partial differential equations. In the absence of a lower order term, we require local conditions on the mesh to establish the comparison principle and uniqueness of the piecewise linear finite element solution. We consider the assembled matrix corresponding to the linearized problem satisfied by the difference of two solutions to the nonlinear problem. Monotonicity of the assembled matrix establishes a maximum principle for the linear problem and a comparison principle for the nonlinear problem. The matrix analysis approach to the discrete comparison principle yields sharper constants and more relaxed mesh conditions than does the argument by contradiction used in previous work.


title: Quantum Quench dynamics in Non-local Luttinger Model: Rigorous Results

abstract: We investigate, in the Luttinger model with fixed box potential, the time evolution of an inhomogeneous state prepared as a localized fermion added to the noninteracting ground state. We proved that, if the state is evolved with the interacting Hamiltonian, the averaged density has two peaks moving in opposite directions, with a constant but renormalized velocity. We also proved that a dynamical `Landau quasi-particle weight' appears in the oscillating part of the averaged density, asymptotically vanishing with large time. The results are proved with the Mattis-Lieb diagonalization method. A simpler proof with the exact Bosonization formulas is also provided.

categories: math-ph cond-mat.stat-mech math.MP math.OA


title: Driving forbidden transitions in the fluxonium artificial atom

abstract: Atomic systems display a rich variety of quantum dynamics due to the different possible symmetries obeyed by the atoms. These symmetries result in selection rules that have been essential for the quantum control of atomic systems. Superconducting artificial atoms are mainly governed by parity symmetry. Its corresponding selection rule limits the types of quantum systems that can be built using electromagnetic circuits at their optimal coherence operation points ("sweet spots"). Here, we use third-order nonlinear coupling between the artificial atom and its readout resonator to drive transitions forbidden by the parity selection rule for linear coupling to microwave radiation. A Lambda-type system emerges from these newly accessible transitions, implemented here in the fluxonium artificial atom coupled to its "antenna" resonator. We demonstrate coherent manipulation of the fluxonium artificial atom at its sweet spot by stimulated Raman transitions. This type of transition enables the creation of new quantum operations, such as the control and readout of physically protected artificial atoms.

categories: quant-ph cond-mat.mes-hall cond-mat.supr-con physics.atom-ph


title: Continuum of quantum fluctuations in a three-dimensional $S\!=\!1$ Heisenberg magnet

abstract: Conventional crystalline magnets are characterized by symmetry breaking and normal modes of excitation called magnons with quantized angular momentum $\hbar$. Neutron scattering correspondingly features extra magnetic Bragg diffraction at low temperatures and dispersive inelastic scattering associated with single magnon creation and annihilation. Exceptions are anticipated in so-called quantum spin liquids as exemplified by the one-dimensional spin-1/2 chain which has no magnetic order and where magnons accordingly fractionalize into spinons with angular momentum $\hbar/2$. This is spectacularly revealed by a continuum of inelastic neutron scattering associated with two-spinon processes and the absence of magnetic Bragg diffraction. Here, we report evidence for these same key features of a quantum spin liquid in the three-dimensional Heisenberg antiferromagnet NaCaNi$_2$F$_7$. Through specific heat and neutron scattering measurements, Monte Carlo simulations, and analytic approximations to the equal time correlations, we show that NaCaNi$_2$F$_7$ is an almost ideal realization of the spin-1 antiferromagnetic Heisenberg model on a pyrochlore lattice with weak connectivity and frustrated interactions. Magnetic Bragg diffraction is absent and 90\% of the spectral weight forms a continuum of magnetic scattering not dissimilar to that of the spin-1/2 chain but with low energy pinch points indicating NaCaNi$_2$F$_7$ is in a Coulomb phase. The residual entropy and diffuse elastic scattering points to an exotic state of matter driven by frustration, quantum fluctuations and weak exchange disorder.


title: Robust Environmental Mapping by Mobile Sensor Networks

abstract: Constructing a spatial map of environmental parameters is a crucial step to preventing hazardous chemical leakages, forest fires, or while estimating a spatially distributed physical quantities such as terrain elevation. Although prior methods can do such mapping tasks efficiently via dispatching a group of autonomous agents, they are unable to ensure satisfactory convergence to the underlying ground truth distribution in a decentralized manner when any of the agents fail. Since the types of agents utilized to perform such mapping are typically inexpensive and prone to failure, this results in poor overall mapping performance in real-world applications, which can in certain cases endanger human safety. This paper presents a Bayesian approach for robust spatial mapping of environmental parameters by deploying a group of mobile robots capable of ad-hoc communication equipped with short-range sensors in the presence of hardware failures. Our approach first utilizes a variant of the Voronoi diagram to partition the region to be mapped into disjoint regions that are each associated with at least one robot. These robots are then deployed in a decentralized manner to maximize the likelihood that at least one robot detects every target in their associated region despite a non-zero probability of failure. A suite of simulation results is presented to demonstrate the effectiveness and robustness of the proposed method when compared to existing techniques.


title: Optimistic Robust Optimization With Applications To Machine Learning

abstract: Robust Optimization has traditionally taken a pessimistic, or worst-case viewpoint of uncertainty which is motivated by a desire to find sets of optimal policies that maintain feasibility under a variety of operating conditions. In this paper, we explore an optimistic, or best-case view of uncertainty and show that it can be a fruitful approach. We show that these techniques can be used to address a wide variety of problems. First, we apply our methods in the context of robust linear programming, providing a method for reducing conservatism in intuitive ways that encode economically realistic modeling assumptions. Second, we look at problems in machine learning and find that this approach is strongly connected to the existing literature. Specifically, we provide a new interpretation for popular sparsity inducing non-convex regularization schemes. Additionally, we show that successful approaches for dealing with outliers and noise can be interpreted as optimistic robust optimization problems. Although many of the problems resulting from our approach are non-convex, we find that DCA or DCA-like optimization approaches can be intuitive and efficient.


title: The Accelerating Pace of Star Formation

abstract: We study the temporal and spatial distribution of star formation rates in four well-studied star-forming regions in local molecular clouds(MCs): Taurus, Perseus, $\rho$ Ophiuchi, and Orion A. Using published mass and age estimates for young stellar objects in each system, we show that the rate of star formation over the last 10 Myrs has been accelerating and is (roughly) consistent with a $t^2$ power law. This is in line with previous studies of the star formation history of molecular clouds and with recent theoretical studies. We further study the clustering of star formation in the Orion Nebula Cluster(ONC). We examine the distribution of young stellar objects as a function of their age by computing an effective half-light radius for these young stars subdivided into age bins. We show that the distribution of young stellar objects is broadly consistent with the star formation being entirely localized within the central region. We also find a slow radial expansion of the newly formed stars at a velocity of $v=0.17\,{\rm km\,s}^{-1}$, which is roughly the sound speed of the cold molecular gas. This strongly suggests the dense structures that form stars persist much longer than the local dynamical time. We argue that this structure is quasi-static in nature and is likely the result of the density profile approaching an attractor solution as suggested by recent analytic and numerical analysis.


title: Self-Similarity Based Time Warping

abstract: In this work, we explore the problem of aligning two time-ordered point clouds which are spatially transformed and re-parameterized versions of each other. This has a diverse array of applications such as cross modal time series synchronization (e.g. MOCAP to video) and alignment of discretized curves in images. Most other works that address this problem attempt to jointly uncover a spatial alignment and correspondences between the two point clouds, or to derive local invariants to spatial transformations such as curvature before computing correspondences. By contrast, we sidestep spatial alignment completely by using self-similarity matrices (SSMs) as a proxy to the time-ordered point clouds, since self-similarity matrices are blind to isometries and respect global geometry. Our algorithm, dubbed "Isometry Blind Dynamic Time Warping" (IBDTW), is simple and general, and we show that its associated dissimilarity measure lower bounds the L1 Gromov-Hausdorff distance between the two point sets when restricted to warping paths. We also present a local, partial alignment extension of IBDTW based on the Smith Waterman algorithm. This eliminates the need for tedious manual cropping of time series, which is ordinarily necessary for global alignment algorithms to function properly.


title: Pomeron Physics at the LHC

abstract: We present current and ongoing research aimed at identifying Pomeron effects at the LHC in both the weak and strongly coupled regimes of QCD.


title: Follower, Predecessor, and Extender Entropies

abstract: Using follower/predecessor/extender set sequences, we define quantities which we call the follower/predecessor/extender entropies, which can be associated to any shift space. We analyze the behavior of these quantities under conjugacies and factor maps, most notably showing that extender entropy is a conjugacy invariant and that having follower entropy zero is a conjugacy invariant. We give some applications, including examples of shift spaces with equal entropy which can be distinguished by extender entropy, and examples of shift spaces which can be shown to not be isomorphic to their inverse by using follower/predecessor entropy.


title: Non-Gaussian Autoregressive Processes with Tukey g-and-h Transformations

abstract: When performing a time series analysis of continuous data, for example from climate or environmental problems, the assumption that the process is Gaussian is often violated. Therefore, we introduce two non-Gaussian autoregressive time series models that are able to fit skewed and heavy-tailed time series data. Our two models are based on the Tukey g-and-h transformation. We discuss parameter estimation, order selection, and forecasting procedures for our models and examine their performances in a simulation study. We demonstrate the usefulness of our models by applying them to two sets of wind speed data.


title: Temperature dependent relaxation of dipole-exchange magnons in yttrium iron garnet films

abstract: Low energy consumption enabled by charge-free information transport, which is free from ohmic heating, and the ability to process phase-encoded data by nanometer-sized interference devices at GHz and THz frequencies are just a few benefits of spin-wave-based technologies. Moreover, when approaching cryogenic temperatures, quantum phenomena in spin-wave systems pave the path towards quantum information processing. In view of these applications, the lifetime of magnons$-$spin-wave quanta$-$is of high relevance for the fields of magnonics, magnon spintronics and quantum computing. Here, the relaxation behavior of parametrically excited magnons having wavenumbers from zero up to $6\cdot 10^5 \mathrm{rad~cm}^{-1}$ was experimentally investigated in the temperature range from 20 K to 340 K in single crystal yttrium iron garnet (YIG) films epitaxially grown on gallium gadolinium garnet (GGG) substrates as well as in a bulk YIG crystal$-$the magnonic materials featuring the lowest magnetic damping known so far. As opposed to the bulk YIG crystal in YIG films we have found a significant increase in the magnon relaxation rate below 150 K$-$up to 10.5 times the reference value at 340 K$-$in the entire range of probed wavenumbers. This increase is associated with rare-earth impurities contaminating the YIG samples with a slight contribution caused by coupling of spin waves to the spin system of the paramagnetic GGG substrate at the lowest temperatures.


title: Treatment Effect Quantification for Time-to-event Endpoints -- Estimands, Analysis Strategies, and beyond

abstract: A draft addendum to ICH E9 has been released for public consultation in August 2017. The addendum focuses on two topics particularly relevant for randomized confirmatory clinical trials: estimands and sensitivity analyses. The need to amend ICH E9 grew out of the realization of a lack of alignment between the objectives of a clinical trial stated in the protocol and the accompanying quantification of the "treatment effect" reported in a regulatory submission. We embed time-to-event endpoints in the estimand framework, and discuss how the four estimand attributes described in the addendum apply to time-to-event endpoints. We point out that if the proportional hazards assumption is not met, the estimand targeted by the most prevalent methods used to analyze time-to-event endpoints, logrank test and Cox regression, depends on the censoring distribution. We discuss for a large randomized clinical trial how the analyses for the primary and secondary endpoints as well as the sensitivity analyses actually performed in the trial can be seen in the context of the addendum. To the best of our knowledge, this is the first attempt to do so for a trial with a time-to-event endpoint. Questions that remain open with the addendum for time-to-event endpoints and beyond are formulated, and recommendations for planning of future trials are given. We hope that this will provide a contribution to developing a common framework based on the final version of the addendum that can be applied to design, protocols, statistical analysis plans, and clinical study reports in the future.


title: Miyachi's Theorem for the Quaternion Fourier Transform

abstract: The quaternion Fourier transform (QFT) satisfies some uncertainty principles similar to the Euclidean Fourier transform. In this paper, we establish Miyachi's theorem for this transform.


title: Dropping Activation Outputs with Localized First-layer Deep Network for Enhancing User Privacy and Data Security

abstract: Deep learning methods can play a crucial role in anomaly detection, prediction, and supporting decision making for applications like personal health-care, pervasive body sensing, etc. However, current architecture of deep networks suffers the privacy issue that users need to give out their data to the model (typically hosted in a server or a cluster on Cloud) for training or prediction. This problem is getting more severe for those sensitive health-care or medical data (e.g fMRI or body sensors measures like EEG signals). In addition to this, there is also a security risk of leaking these data during the data transmission from user to the model (especially when it's through Internet). Targeting at these issues, in this paper we proposed a new architecture for deep network in which users don't reveal their original data to the model. In our method, feed-forward propagation and data encryption are combined into one process: we migrate the first layer of deep network to users' local devices, and apply the activation functions locally, and then use "dropping activation output" method to make the output non-invertible. The resulting approach is able to make model prediction without accessing users' sensitive raw data. Experiment conducted in this paper showed that our approach achieves the desirable privacy protection requirement, and demonstrated several advantages over the traditional approach with encryption / decryption


title: Medium-mass hypernuclei and the nucleon-isospin dependence of the three-body hyperon-nucleon-nucleon force

abstract: We report quantum Monte Carlo calculations of single-$\Lambda$ hypernuclei for $A<50$ based on phenomenological two- and three-body hyperon-nucleon forces. We present results for the $\Lambda$ separation energy in different hyperon orbits, showing that the accuracy of theoretical predictions exceeds that of currently available experimental data, especially for medium-mass hypernuclei. We show the results of a sensitivity study that indicates the possibility to investigate the nucleon-isospin dependence of the three-body hyperon-nucleon-nucleon force in the medium-mass region of the hypernuclear chart, where new spectroscopy studies are currently planned. The importance of such a dependence for the description of the physics of hypernuclei, and the consequences for the prediction of neutron star properties are discussed.


title: Gemini IFU, VLA and HST observations of the OH Megamaser galaxy IRASF23199+0123: the hidden monster and its outflow

abstract: We present Gemini Multi-Object Spectrograph (GMOS) Integral field Unit (IFU), Very Large Array (VLA) and Hubble Space Telescope (HST) observations of the OH Megamaser (OHM) galaxy IRASF23199+0123. Our observations show that this system is an interacting pair, with two OHM sources associated to the eastern (IRAS23199E) member. The two members of the pair present somewhat extended radio emission at 3 and 20~cm, with flux peaks at each nucleus. The GMOS-IFU observations cover the inner $\sim$6kpc of IRAS23199E at a spatial resolution of 2.3~kpc. The GMOS-IFU flux distributions in H$\alpha$ and [NII]$\lambda$6583 are similar to that of an HST [NII]+H$\alpha$ narrow-band image, being more extended along the northeast-southwest direction, as also observed in the continuum HST F814W image. The GMOS-IFU H$\alpha$ flux map of IRAS23199E shows three extranuclear knots attributed to star-forming complexes. We have discovered a Seyfert 1 nucleus in this galaxy, as its nuclear spectrum shows an unresolved broad (FWHM$\approx$2170 kms$^{-1}$) double-peaked H$\alpha$ component, from which we derive a black hole mass of M$_{BH}$= 3.8$^{+0.3}_{-0.2}\times 10^{6}$M$_{\odot}$. The gas kinematics shows low velocity dispersions ($\sigma$) and low [NII]/H$\alpha$ ratios for the star-forming complexes and higher $\sigma$ and [NII]/H$\alpha$ surrounding the radio emission region, supporting interaction between the radio-plasma and ambient gas. The two OH masers detected in IRASF23199E are observed in the vicinity of these enhanced $\sigma$ regions, supporting their association with the active nucleus and its interaction with the surrounding gas. The gas velocity field can be partially reproduced by rotation in a disk, with residuals along the north-south direction being tentatively attributed to emission from the front walls of a bipolar outflow.


title: Phase Synchronization on Spacially Embeded Duplex Networks with Total Cost Constraint

abstract: Synchronization on multiplex networks have attracted increasing attention in the past few years. We investigate collective behaviors of Kuramoto oscillators on single layer and duplex spacial networks with total cost restriction, which was introduced by Li et. al [Li G., Reis S. D., Moreira A. A., Havlin S., Stanley H. E. and Jr A. J., {\it Phys. Rev. Lett.} 104, 018701 (2010)] and termed as the Li network afterwards. In the Li network model, with the increase of its spacial exponent, the network's structure will vary from the random type to the small-world one, and finally to the regular lattice.We first explore how the spacial exponent influences the synchronizability of Kuramoto oscillators on single layer Li networks and find that the closer the Li network is to a regular lattice, the more difficult for it to evolve into synchronization. Then we investigate synchronizability of duplex Li networks and find that the existence of inter-layer interaction can greatly enhance inter-layer and global synchronizability. When the inter-layer coupling strength is larger than a certain critical value, whatever the intra-layer coupling strength is, the inter-layer synchronization will always occur. Furthermore, on single layer Li networks, nodes with larger degrees more easily reach global synchronization, while on duplex Li networks, this phenomenon becomes much less obvious. Finally, we study the impact of inter-link density on global synchronization and obtain that sparse inter-links can lead to the emergence of global synchronization for duplex Li networks just as dense inter-links do. In a word, inter-layer interaction plays a vital role in determining synchronizability for duplex spacial networks with total cost constraint.


title: On k-neighbor separated permutations

abstract: Two permutations of $[n]=\{1,2 \ldots n\}$ are \textit{$k$-neighbor separated} if there are two elements that are neighbors in one of the permutations and that are separated by exactly $k-2$ other elements in the other permutation. Let the maximal number of pairwise $k$-neighbor separated permutations of $[n]$ be denoted by $P(n,k)$. In a previous paper, the authors have determined $P(n,3)$ for every $n$, answering a question of K\"orner, Messuti and Simonyi affirmatively. In this paper we prove that for every fixed positive integer $\ell $, $$P(n,2^\ell+1) = 2^{n-o(n)}. $$ We conjecture that for every fixed even $k$, $P(n,k)=2^{n-o(n)}$. We also show that this conjecture is asymptotically true in the following sense $$\lim_{k \rightarrow \infty} \lim_{n \rightarrow \infty} \sqrt[n]{P(n,k)}=2.$$ Finally, we show that for even $n$, $P(n,n)= 3n/2$.


title: Glivenko-Cantelli Theory, Ornstein-Weiss quasi-tilings, and uniform Ergodic Theorems for distribution-valued fields over amenable groups

abstract: We consider random fields indexed by finite subsets of an amenable discrete group, taking values in the Banach-space of bounded right-continuous functions. The field is assumed to be equivariant, local, coordinate-wise monotone, and almost additive, with finite range dependence. Using the theory of quasi-tilings we prove an uniform ergodic theorem, more precisely, that averages along a Foelner sequence converge uniformly to a limiting function. Moreover we give explicit error estimates for the approximation in the sup norm.


title: Radiative corrections to masses and couplings in Universal Extra Dimensions

abstract: Models with an orbifolded universal extra dimension receive important loop-induced corrections to the masses and couplings of Kaluza-Klein (KK) particles. The dominant contributions stem from so-called boundary terms which violate KK number. Previously, only the parts of these boundary terms proportional to $\ln(\Lambda R)$ have been computed, where $R$ is the radius of the extra dimension and $\Lambda$ is cut-off scale. However, for typical values of $\Lambda R \sim 10 \cdots 50$, the logarithms are not particularly large and non-logarithmic contributions may be numerically important. In this paper, these remaining finite terms are computed and their phenomenological impact is discussed. It is shown that the finite terms have a significant impact on the KK mass spectrum. Furthermore, one finds new KK-number violating interactions that do not depend on $\ln(\Lambda R)$ but nevertheless are non-zero. These lead to new production and decay channels for level-2 KK particles at colliders.


title: Subgroup Identification and Interpretation with Bayesian Nonparametric Models in Health Care Claims Data

abstract: Inpatient care is a large share of total health care spending, making analysis of inpatient utilization patterns an important part of understanding what drives health care spending growth. Common features of inpatient utilization measures include zero inflation, over-dispersion, and skewness, all of which complicate statistical modeling. Mixture modeling is a popular approach that can accommodate these features of health care utilization data. In this work, we add a nonparametric clustering component to such models. Our fully Bayesian model framework allows for an unknown number of mixing components, so that the data determine the number of mixture components. When we apply the modeling framework to data on hospital lengths of stay for patients with lung cancer, we find distinct subgroups of patients with differences in means and variances of hospital days, health and treatment covariates, and relationships between covariates and length of stay.


title: Cluster tilting subcategories and torsion pairs in Igusa--Todorov cluster categories of Dynkin type $A_{ \infty }$

abstract: We give a combinatorial classification of cluster tilting subcategories and torsion pairs in Igusa--Todorov cluster categories of Dynkin type $A_{ \infty }$.


title: Dissipativity of system abstractions obtained using approximate input-output simulation

abstract: This work focuses on the invariance of important properties between continuous and discrete models of systems which can be useful in the control design of large-scale systems and their software implementations. In particular, this paper discusses the relationships between the QSR dissipativity of a continuous state dynamical system and of its abstractions obtained through approximate input-output simulation relations. First, conditions to guarantee the dissipativity of the continuous system from its abstractions are provided. The reverse problem of determining the Q, S and R dissipativity matrices of the abstract system from that of the continuous system is also considered. Results characterizing the change in the dissipativity matrices are provided when the system abstraction is obtained. Since, under certain conditions, QSR dissipative systems are known to be stable, the results of this paper can be used to construct stable system abstractions as well. In the second part of this paper, we analyze the dissipativity of the approximate feedback composition of a continuous dynamical system and a discrete controller. We present illustrative examples to demonstrate the results of this paper.


title: Online Resource Inference in Network Utility Maximization Problems

abstract: The amount of transmitted data in computer networks is expected to grow considerably in the future, putting more and more pressure on the network infrastructures. In order to guarantee a good service, it then becomes fundamental to use the network resources efficiently. Network Utility Maximization (NUM) provides a framework to optimize the rate allocation when network resources are limited. Unfortunately, in the scenario where the amount of available resources is not known a priori, classical NUM solving methods do not offer a viable solution. To overcome this limitation we design an overlay rate allocation scheme that attempts to infer the actual amount of available network resources while coordinating the users rate allocation. Due to the general and complex model assumed for the congestion measurements, a passive learning of the available resources would not lead to satisfying performance. The coordination scheme must then perform active learning in order to speed up the resources estimation and quickly increase the system performance. By adopting an optimal learning formulation we are able to balance the tradeoff between an accurate estimation, and an effective resources exploitation in order to maximize the long term quality of the service delivered to the users.


title: Development of Solar Flares and Features of the Fine Structure of Solar Radio Emission

abstract: The reason for the occurrence of different elements of the fine structure of solar radio bursts in the decimeter and centimeter wavelength ranges has been determined based on all available data from terrestrial and satellite observations. In some phenomena, fast pulsations, a zebra structure, fiber bursts, and spikes have been observed almost simultaneously. Two phenomena have been selected to show that the pulsations of radio emission are caused by particles accelerated in the magnetic reconnection region and that the zebra structure is excited in a source, such as a magnetic trap for fast particles. The complex combination of unusual fiber bursts, zebra structure, and spikes in the phenomenon on December 1, 2004, is associated with a single source, a magnetic island formed after a coronal mass ejection.


title: Path properties of the solution to the stochastic heat equation with L\'evy noise

abstract: We consider sample path properties of the solution to the stochastic heat equation, in $\mathbb{R}^d$ or bounded domains of $\mathbb{R}^d$, driven by a L\'evy space-time white noise. When viewed as a stochastic process in time with values in an infinite-dimensional space, the solution is shown to have a c\`adl\`ag modification in fractional Sobolev spaces of index less than $-\frac d 2$. Concerning the partial regularity of the solution in time or space when the other variable is fixed, we determine critical values for the Blumenthal-Getoor index of the L\'evy noise such that noises with a smaller index entail continuous sample paths, while L\'evy noises with a larger index entail sample paths that are unbounded on any non-empty open subset. Our results apply to additive as well as multiplicative L\'evy noises, and to light- as well as heavy-tailed jumps.


title: Decay of homogeneous two dimensional quantum turbulence

abstract: We numerically simulate the free decay of two-dimensional quantum turbulence in a large, homogeneous Bose-Einstein condensate. The large number of vortices, the uniformity of the density profile and the absence of boundaries (where vortices can drift out of the condensate) isolate the annihilation of vortex-antivortex pairs as the only mechanism which reduces the number of vortices, $\Nv$, during the turbulence decay. The results clearly reveal that vortex annihilation is a four-vortex process, confirming the decay law $\Nv \sim t^{-1/3}$ where $t$ is time, which was inferred from experiments with relatively few vortices in small harmonically trapped condensates.


title: Solar-wind predictions for the Parker Solar Probe orbit

abstract: The scope of this study is to model the solar-wind environment for the Parker Solar Probe's unprecedented distances down to 9.86 Rs in its mission phase during 2018-2025. The study is performed within the CGAUSS project which is the German contribution to the PSP mission as part of the WISPR imager on PSP. We present an empirical solar-wind model for the inner heliosphere which is derived from OMNI and Helios data. The sunspot number (SSN) and its predictions are used to derive dependencies of solar-wind parameters on solar activity and to forecast them for the PSP mission. The frequency distributions for the solar-wind key parameters magnetic field strength, proton velocity, density, and temperature, are represented by lognormal functions, considering the velocity distribution's bi-componental shape. Functional relations to the SSN are compiled using OMNI data and based on data from both Helios probes, the parameters' frequency distributions are fitted with respect to solar distance. Thus, an empirical solar-wind model for the inner heliosphere is derived, accounting for solar activity and solar distance. The inclusion of SSN predictions and the extrapolation down to PSP's perihelion region enables us to estimate the solar-wind environment for PSP's planned trajectory during its mission duration. This empirical model yields estimated solar-wind values for PSP's 1st perihelion in 2018 at 0.16 au: 87 nT, 340 km s-1, 214 cm-3 and 503000 K. The estimates for PSP's first closest perihelion, occurring in 2024 at 0.046 au, are 943 nT, 290 km s-1, 2951 cm-3, and 1930000 K. Since the modeled velocity and temperature values below approximately 20 Rs appear overestimated in comparison with existing observations, this suggests that PSP will directly measure solar-wind acceleration and heating processes below 20 Rs as planned.


title: 1I/2017 U1 (Oumuamua) Might Be A Cometary Nucleus

abstract: In this work we find evidence that the object is of cometary origin.


title: Temperature-reflection I: field theory, ensembles, and interactions

abstract: In this paper, we revisit the claim that many partition functions are invariant under reflecting temperatures to negative values (T-reflection). The goal of this paper is to demarcate which partition functions should be invariant under T-reflection, and why. Our main claim is that finite-temperature path integrals for quantum field theories (QFTs) should be T-reflection invariant. Because multi-particle partition functions are equal to Euclidean path integrals for QFTs, we expect them to be T-reflection invariant. Single-particle partition functions though are often not invariant under T-reflection. Several exactly solvable systems are non-invariant under naive T-reflection, but are likely invariant under an extended T-reflection. We give example systems that are T-reflection invariant but are (1) non-unitary, (2) chiral, (3) interacting, (4) non-supersymmetric, or (5) non-conformal, and (6) argue that T-reflection is unrelated to time-reversal. Finally, we study the interplay between T-reflection and perturbation theory in the anharmonic harmonic oscillator in quantum mechanics and in Yang-Mills in four-dimensions. This is the first in a series of papers on temperature-reflections.


title: Topological Susceptibility in $N_f=2$ QCD at Finite Temperature

abstract: We study the topological charge in $N_f=2$ QCD at finite temperature using M\"obius domain-wall fermions. The susceptibility $\chi_t$ of the topological charge defined either by the index of overlap Dirac operator or a gluonic operator is investigated at several values of temperature $T\ (>T_c)$ varying the quark mass. A strong suppression of the susceptibility is observed below a certain value of the quark mass. The relation with the restoration of $U_A(1)$ is discussed.


title: The broad band SEDs of four 'hypervariable' AGN

abstract: We present an optical to X-ray spectral analysis of four 'hypervariable' AGN (HVAs) discovered by comparing Pan-STARRS data to that from SDSS over a 10 year baseline (Lawrence et al 2016). There is some evidence that these objects are X-ray loud for their corresponding UV luminosities, but given that we measured them in a historic high state, it is not clear whether to take the high-state or low-state as typical of the properties of these HVAs. We estimate black hole masses based on MgII and H$\alpha$ emission line profiles, and either the high and low state luminosities, finding mass ranges $\log(M_{\rm BH}/M_{\odot}) = 8.2-8.8$ and $\log(M_{\rm BH}/M_{\odot}) = 7.9-8.3$ respectively. We then fit energy conserving models to the SEDs, obtaining strong constraints on the bolometric luminosity and $\alpha_{\rm OX}$. We compare the SED properties with a larger, X-ray selected AGN sample for both of these scenarios, and observe distinct groupings in spectral shape versus luminosity parameter space. In general, the SED properties are closer to normal if we assume that the low-state is representative. This supports the idea that the large slow outbursts may be due to extrinsic effects (for example microlensing) as opposed to accretion rate changes, but a larger sample of HVAs is needed to be confident of this conclusion.


title: Transition density estimates for diagonal systems of SDEs driven by cylindrical $\alpha$-stable processes

abstract: We consider the system of stochastic differential equation $dX_t = A(X_{t-}) \, dZ_t$, $ X_0 = x$, driven by cylindrical $\alpha$-stable process $Z_t$ in $\mathbb{R}^d$. We assume that $A(x) = (a_{ij}(x))$ is diagonal and $a_{ii}(x)$ are bounded away from zero, from infinity and H\"older continuous. We construct transition density $p^A(t,x,y)$ of the process $X_t$ and show sharp two-sided estimates of this density. We also prove H\"older and gradient estimates of $x \to p^A(t,x,y)$. Our approach is based on the method developed by Chen and Zhang.


title: Post-prior equivalence for transfer reactions with complex potentials

abstract: In this paper, we address the problem of the post-prior equivalence in the calculation of inclusive breakup and transfer cross sections. For that, we employ the model proposed by Ichimura, Austern, and Vincent [Phys. Rev. C 32, 431 (1985)], conveniently generalized to include the part of the cross section corresponding the transfer to bound states. We pay particular attention to the case in which the unobserved particle is left in a bound state of the residual nucleus, in which case the theory prescribes the use of a complex potential, responsible for the spreading width of the populated single-particle states. We see that the introduction of this complex potential gives rise to an additional term in the prior cross section formula, not present in the usual case of real binding potentials. The equivalence is numerically tested for reaction induced by deuterons.