title: Solving differential equations with neural networks: Applications to the calculation of cosmological phase transitions

abstract: Starting from the observation that artificial neural networks are uniquely suited to solving optimisation problems, and most physics problems can be cast as an optimisation task, we introduce a novel way of finding a numerical solution to wide classes of differential equations. We find our approach to be very flexible and stable without relying on trial solutions, and applicable to ordinary, partial and coupled differential equations. We apply our method to the calculation of tunnelling profiles for cosmological phase transitions, which is a problem of relevance for baryogenesis and stochastic gravitational wave spectra. Comparing our solutions with publicly available codes which use numerical methods optimised for the calculation of tunnelling profiles, we find our approach to provide at least as accurate results as these dedicated differential equation solvers, and for some parameter choices even more accurate and reliable solutions. In particular, we compare the neural network approach with two publicly available profile solvers, \texttt{CosmoTransitions} and \texttt{BubbleProfiler}, and give explicit examples where the neural network approach finds the correct solution while dedicated solvers do not. We point out that this approach of using artificial neural networks to solve equations is viable for any problem that can be cast into the form $\mathcal{F}(\vec{x})=0$, and is thus applicable to various other problems in perturbative and non-perturbative quantum field theory.


title: On $SU(3)_{F}$ Breaking through Final State Interactions and CP Asymmetries in $D\to P P$ Decays

abstract: We analyse $D$ decays to two pseudoscalars $(\pi,K)$ assuming the dominant source of $\mathrm{SU}(3)_F$ breaking lies in final state interactions. We obtain an excellent agreement with experimental data and are able to predict CP violation in several channels based on current data on branching ratios and $\Delta \mathrm{A_{CP}}$. We also make predictions for $\delta_{K\pi}$ and the branching fraction for the decay $D_s^+\to K^+K_L$.


title: Curvature Invariants and Lower Dimensional Black Hole Horizons

abstract: It is known that the event horizon of a black hole can often be identified from the zeroes of some curvature invariants. The situation in lower dimensions has not been thoroughly clarified. In this work we investigate both (2+1)- and (1+1)-dimensional black hole horizons of static, stationary and dynamical black holes, identified with the zeroes of scalar polynomial and Cartan curvature invariants, with the purpose of discriminating the different roles played by the Weyl and Riemann curvature tensors. The situations and applicability of the methods are found to be quite different from that in 4-dimensional spacetime. The suitable Cartan invariants employed for detecting the horizon can be interpreted as a local extremum of the tidal force suggesting that the horizon of a black hole is a genuine special hypersurface within the full manifold, contrary to the usual claim that there is nothing special at the horizon, which is said to be a consequence of the equivalence principle.


title: Interest-Related Item Similarity Model Based on Multimodal Data for Top-N Recommendation

abstract: Nowadays, the recommendation systems are applied in the fields of e-commerce, video websites, social networking sites, etc., which bring great convenience to people's daily lives. The types of the information are diversified and abundant in recommendation systems, therefore the proportion of unstructured multimodal data like text, image and video is increasing. However, due to the representation gap between different modalities, it is intractable to effectively use unstructured multimodal data to improve the efficiency of recommendation systems. In this paper, we propose an end-to-end Multimodal Interest-Related Item Similarity model (Multimodal IRIS) to provide recommendations based on multimodal data source. Specifically, the Multimodal IRIS model consists of three modules, i.e., multimodal feature learning module, the Interest-Related Network (IRN) module and item similarity recommendation module. The multimodal feature learning module adds knowledge sharing unit among different modalities. Then IRN learn the interest relevance between target item and different historical items respectively. At last, the multimodal data feature learning, IRN and item similarity recommendation modules are unified into an integrated system to achieve performance enhancements and to accommodate the addition or absence of different modal data. Extensive experiments on real-world datasets show that, by dealing with the multimodal data which people may pay more attention to when selecting items, the proposed Multimodal IRIS significantly improves accuracy and interpretability on top-N recommendation task over the state-of-the-art methods.


title: Discovering Archetypes to Interpret Evolution of Individual Behavior

abstract: In this paper, we aim to discover archetypical patterns of individual evolution in large social networks. In our work, an archetype comprises of $\textit{progressive stages}$ of distinct behavior. We introduce a novel Gaussian Hidden Markov Model (G-HMM) Cluster to identify archetypes of evolutionary patterns. G-HMMs allow for: near limitless behavioral variation; imposing constraints on how individuals can evolve; different evolutionary rates; and are parsimonious. Our experiments with Academic and StackExchange dataset discover insightful archetypes. We identify four archetypes for researchers: $\textit{Steady}$, $\textit{Diverse, Evolving and Diffuse}$. We observe clear differences in the evolution of male and female researchers within the same archetype. Specifically, women and men differ within an archetype (e.g. Diverse) in how they start, how they transition and the time spent in mid-career. We also found that the differences in grant income are better explained by the differences in archetype than by differences in gender. For StackOverflow, discovered archetypes could be labeled as $\textit{Experts, Seekers, Enthusiasts, and Facilitators}$. We have strong quantitative results with competing baselines for activity prediction and perplexity. For future session prediction, the proposed G-HMM cluster model improves by an average of $32\%$ for different Stack Exchanges and $24\%$ for Academic dataset. Our model also exhibits lower perplexity than the baselines.


title: Tractography and machine learning: Current state and open challenges

abstract: Supervised machine learning (ML) algorithms have recently been proposed as an alternative to traditional tractography methods in order to address some of their weaknesses. They can be path-based and local-model-free, and easily incorporate anatomical priors to make contextual and non-local decisions that should help the tracking process. ML-based techniques have thus shown promising reconstructions of larger spatial extent of existing white matter bundles, promising reconstructions of less false positives, and promising robustness to known position and shape biases of current tractography techniques. But as of today, none of these ML-based methods have shown conclusive performances or have been adopted as a de facto solution to tractography. One reason for this might be the lack of well-defined and extensive frameworks to train, evaluate, and compare these methods. In this paper, we describe several datasets and evaluation tools that contain useful features for ML algorithms, along with the various methods proposed in the recent years. We then discuss the strategies that are used to evaluate and compare those methods, as well as their shortcomings. Finally, we describe the particular needs of ML tractography methods and discuss tangible solutions for future works.


title: The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s

abstract: The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.


title: Reinforcement Learning to Optimize Long-term User Engagement in Recommender Systems

abstract: Recommender systems play a crucial role in our daily lives. Feed streaming mechanism has been widely used in the recommender system, especially on the mobile Apps. The feed streaming setting provides users the interactive manner of recommendation in never-ending feeds. In such an interactive manner, a good recommender system should pay more attention to user stickiness, which is far beyond classical instant metrics, and typically measured by {\bf long-term user engagement}. Directly optimizing the long-term user engagement is a non-trivial problem, as the learning target is usually not available for conventional supervised learning methods. Though reinforcement learning~(RL) naturally fits the problem of maximizing the long term rewards, applying RL to optimize long-term user engagement is still facing challenges: user behaviors are versatile and difficult to model, which typically consists of both instant feedback~(e.g. clicks, ordering) and delayed feedback~(e.g. dwell time, revisit); in addition, performing effective off-policy learning is still immature, especially when combining bootstrapping and function approximation. To address these issues, in this work, we introduce a reinforcement learning framework --- FeedRec to optimize the long-term user engagement. FeedRec includes two components: 1)~a Q-Network which designed in hierarchical LSTM takes charge of modeling complex user behaviors, and 2)~an S-Network, which simulates the environment, assists the Q-Network and voids the instability of convergence in policy learning. Extensive experiments on synthetic data and a real-world large scale data show that FeedRec effectively optimizes the long-term user engagement and outperforms state-of-the-arts.


title: The High Time Resolution Universe Pulsar Survey -- XV: completion of the intermediate latitude survey with the discovery and timing of 25 further pulsars

abstract: We report on the latest six pulsars discovered through our standard pipeline in the intermediate-latitude region (|b| < 15 deg) of the Parkes High Time Resolution Universe Survey (HTRU). We also present timing solutions for the new discoveries and for 19 further pulsars for which only discovery parameters were previously published. Highlights of the presented sample include the isolated millisecond pulsar J1826-2415, the long-period binary pulsar J1837-0822 in a mildly eccentric 98-day orbit with a > 0.27 M_sun companion, and the nulling pulsar J1638-4233, detected only 10% of the time. Other interesting objects are PSR J1757-1500, exhibiting sporadic mode changes, and PSR J1635-2616 showing one glitch over 6 years. The new discoveries bring the total count of HTRU intermediate-latitude pulsars to 113, 25% of which are recycled pulsars. This is the higest ratio of recycled over ordinary pulsars discoveries of all recent pulsar surveys in this region of the sky. Among HTRU recycled pulsars, four are isolated objects. Comparing the characteristics of Galactic fully-recycled isolated MSPs with those of eclipsing binaries ('spiders'), from which the former are believed to have formed, we highlight a discrepancy in their spatial distribution. This may reflect a difference in the natal kick, hence, possibly, a different formation path. On the other hand, however, isolated fully-recycled MSPs spin periods are, on average, longer than those of spiders, in line with what one would expect, from simple magnetic-dipole spin-down, if the former were indeed evolved from the latter.


title: Direct determination of mode-projected electron-phonon coupling in the time-domain

abstract: Ultrafast spectroscopies have become an important tool for elucidating the microscopic description and dynamical properties of quantum materials. In particular, by tracking the dynamics of non-thermal electrons, a material's dominant scattering processes -- and thus the many-body interactions between electrons and collective excitations -- can be revealed. Here we present a new method for extracting the electron-phonon coupling strength in the time domain, by means of time and angle-resolved photoemission spectroscopy (TR-ARPES). This method is demonstrated in graphite, where we investigate the dynamics of photo-injected electrons at the K point, detecting quantized energy-loss processes that correspond to the emission of strongly-coupled optical phonons. We show that the observed characteristic timescale for spectral-weight-transfer mediated by phonon-scattering processes allows for the direct quantitative extraction of electron-phonon matrix elements, for specific modes, and with unprecedented sensitivity.


title: Computational framework for applying electrical impedance tomography to head imaging

abstract: This work introduces a computational framework for applying absolute electrical impedance tomography to head imaging without accurate information on the head shape or the electrode positions. A library of fifty heads is employed to build a principal component model for the typical variations in the shape of the human head, which leads to a relatively accurate parametrization for head shapes with only a few free parameters. The estimation of these shape parameters and the electrode positions is incorporated in a regularized Newton-type output least squares reconstruction algorithm. The presented numerical experiments demonstrate that strong enough variations in the internal conductivity of a human head can be detected by absolute electrical impedance tomography even if the geometric information on the measurement configuration is incomplete to an extent that is to be expected in practice.


title: A Consistent Reduced Network for HCN Chemistry in Early Earth and Titan Atmospheres: Quantum Calculations of Reaction Rate Coefficients

abstract: HCN is a key ingredient for synthesizing biomolecules such as nucleobases and amino acids. We calculate 42 reaction rate coefficients directly involved with or in competition with the production of HCN in the early Earth or Titan atmospheres. These reactions are driven by methane and nitrogen radicals produced via UV photodissociation or lightning. For every reaction in this network, we calculate rate coefficients at 298 K using canonical variational transition state theory (CVT) paired with computational quantum chemistry simulations at the BHandHLYP/augcc-pVDZ level of theory. We also calculate the temperature dependence of the rate coefficients for the reactions that have barriers from 50 to 400 K. We present 15 new reaction rate coefficients with no previously known value; 93% of our calculated coefficients are within an order of magnitude of the nearest experimental or recommended values. Above 320 K, the rate coefficient for the new reaction H2CN -> HCN + H dominates. Contrary to experiments, we find the HCN reaction pathway, N + CH3 -> HCN + H2, to be inefficient and suggest that the experimental rate coefficient actually corresponds to an indirect pathway, through the H2CN intermediate. We present CVT using energies computed with density functional theory as a feasible and accurate method for calculating a large network of rate coefficients of small-molecule reactions.


title: Fully Convolutional Networks for Text Classification

abstract: In this work I propose a new way of using fully convolutional networks for classification while allowing for input of any size. I additionally propose two modifications on the idea of attention and the benefits and detriments of using the modifications. Finally, I show suboptimal results on the ITAmoji 2018 tweet to emoji task and provide a discussion about why that might be the case as well as a proposed fix to further improve results.


title: Sub-micron spatial resolution in far-field Raman imaging via positivity constrained super-resolution

abstract: Raman microscopy is a valuable tool for detecting physical and chemical properties of a sample material. When probing nanomaterials or nanocomposites the spatial resolution of Raman microscopy is not always adequate as it is limited by the optical diffraction limit. Numerical post-processing with super-resolution algorithms provides a means to enhance resolution and can be straightforwardly applied. The aim of this work is to present interior-point least squares (IPLS) as a powerful tool for super-resolution in Raman imaging through constrained optimisation. IPLS's potential for super-resolution is illustrated on numerically generated test images. Its resolving power is demonstrated on Raman spectroscopic data of a polymer nanowire sample. Comparison to AFM data of the same sample substantiates that the presented method is a promising technique for analysing nanomaterial samples.


title: A Scalable Platform for Distributed Object Tracking across a Many-camera Network

abstract: Advances in deep neural networks (DNN) and computer vision (CV) algorithms have made it feasible to extract meaningful insights from large-scale deployments of urban cameras. Tracking an object of interest across the camera network in near real-time is a canonical problem. However, current tracking platforms have two key limitations: 1) They are monolithic, proprietary and lack the ability to rapidly incorporate sophisticated tracking models; and 2) They are less responsive to dynamism across wide-area computing resources that include edge, fog and cloud abstractions. We address these gaps using Anveshak, a runtime platform for composing and coordinating distributed tracking applications. It provides a domain-specific dataflow programming model to intuitively compose a tracking application, supporting contemporary CV advances like query fusion and re-identification, and enabling dynamic scoping of the camera network's search space to avoid wasted computation. We also offer tunable batching and data-dropping strategies for dataflow blocks deployed on distributed resources to respond to network and compute variability. These balance the tracking accuracy, its real-time performance and the active camera-set size. We illustrate the concise expressiveness of the programming model for $4$ tracking applications. Our detailed experiments for a network of 1000 camera-feeds on modest resources exhibit the tunable scalability, performance and quality trade-offs enabled by our dynamic tracking, batching and dropping strategies.


title: A Probabilistic framework for Quantum Clustering

abstract: Quantum Clustering is a powerful method to detect clusters in data with mixed density. However, it is very sensitive to a length parameter that is inherent to the Schr\"odinger equation. In addition, linking data points into clusters requires local estimates of covariance that are also controlled by length parameters. This raises the question of how to adjust the control parameters of the Schr\"odinger equation for optimal clustering. We propose a probabilistic framework that provides an objective function for the goodness-of-fit to the data, enabling the control parameters to be optimised within a Bayesian framework. This naturally yields probabilities of cluster membership and data partitions with specific numbers of clusters. The proposed framework is tested on real and synthetic data sets, assessing its validity by measuring concordance with known data structure by means of the Jaccard score (JS). This work also proposes an objective way to measure performance in unsupervised learning that correlates very well with JS.


title: Quantum correlations in optomechanical crystals

abstract: The field of optomechanics provides us with several examples of quantum photon-phonon interface. In this paper, we theoretically investigate the generation and manipulation of quantum correlations in a microfabricated optomechanical array. We consider a system consisting of localized photonic and phononic modes interacting locally via radiation pressure at each lattice site with the possibility of hopping of photons and phonons between neighboring sites. We show that such an interaction can correlate various modes of a driven coupled optomechanical array with well-chosen system parameters. Moreover, in the linearized regime of Gaussian fluctuations, the quantum correlations not only survive in the presence of thermal noise, but may also be generated thermally. We find that these optomechanical arrays provide a suitable platform for quantum simulation of various many-body systems.


title: Phonon-assisted tunnelling in a double quantum dot molecule immersed in a cavity

abstract: The effects caused by phonon-assisted tunnelling (PhAT) in a double quantum dot (QD) molecule immersed in a cavity were studied under the quantum Markovian master equation formalism in order to account for dissipation phenomena. We explain how for higher PhAT rates, a stronger interaction between a QD and the cavity at off-resonance takes place through the resonant interaction of another QD and the cavity, where the QDs interact through tunnelling. A closer look at the system's allowed optical transitions as a function of PhAT rates showed coalescence of the transition energies, leading to a loss of spectroscopic resolution in the photoluminescence (PL) spectrum. This coalescence occurs at several exceptional points, suggesting a quantum dynamical phase transition.


title: Adversarially Approximated Autoencoder for Image Generation and Manipulation

abstract: Regularized autoencoders learn the latent codes, a structure with the regularization under the distribution, which enables them the capability to infer the latent codes given observations and generate new samples given the codes. However, they are sometimes ambiguous as they tend to produce reconstructions that are not necessarily faithful reproduction of the inputs. The main reason is to enforce the learned latent code distribution to match a prior distribution while the true distribution remains unknown. To improve the reconstruction quality and learn the latent space a manifold structure, this work present a novel approach using the adversarially approximated autoencoder (AAAE) to investigate the latent codes with adversarial approximation. Instead of regularizing the latent codes by penalizing on the distance between the distributions of the model and the target, AAAE learns the autoencoder flexibly and approximates the latent space with a simpler generator. The ratio is estimated using generative adversarial network (GAN) to enforce the similarity of the distributions. Additionally, the image space is regularized with an additional adversarial regularizer. The proposed approach unifies two deep generative models for both latent space inference and diverse generation. The learning scheme is realized without regularization on the latent codes, which also encourages faithful reconstruction. Extensive validation experiments on four real-world datasets demonstrate the superior performance of AAAE. In comparison to the state-of-the-art approaches, AAAE generates samples with better quality and shares the properties of regularized autoencoder with a nice latent manifold structure.


title: Improving Catheter Segmentation & Localization in 3D Cardiac Ultrasound Using Direction-Fused FCN

abstract: Fast and accurate catheter detection in cardiac catheterization using harmless 3D ultrasound (US) can improve the efficiency and outcome of the intervention. However, the low image quality of US requires extra training for sonographers to localize the catheter. In this paper, we propose a catheter detection method based on a pre-trained VGG network, which exploits 3D information through re-organized cross-sections to segment the catheter by a shared fully convolutional network (FCN), which is called a Direction-Fused FCN (DF-FCN). Based on the segmented image of DF-FCN, the catheter can be localized by model fitting. Our experiments show that the proposed method can successfully detect an ablation catheter in a challenging ex-vivo 3D US dataset, which was collected on the porcine heart. Extensive analysis shows that the proposed method achieves a Dice score of 57.7%, which offers at least an 11.8 % improvement when compared to state-of-the-art instrument detection methods. Due to the improved segmentation performance by the DF-FCN, the catheter can be localized with an error of only 1.4 mm.


title: Activity-Enhanced Self-Assembly of a Colloidal Kagome Lattice

abstract: Here, we describe a method for the enhanced self-assembly of triblock Janus colloids targeted to form a kagome lattice. Using computer simulations, we demonstrate that the formation of this elusive structure can be significantly improved by self-propelling or activating the colloids along the axis connecting their hydrophobic hemispheres. The process by which metastable aggregates are destabilized and transformed into the favored kagome lattice is quite general, and we argue this active approach provides a systematic pathway to improving the self-assembly of a large number of colloidal structures.


title: The Strong Maximum Principle for Schr\"{o}dinger operators on fractals

abstract: We prove a strong maximum principle for Schr\"odinger operators defined on a class of fractal sets and their blowups without boundary. Our primary interest is in weaker regularity conditions than have previously appeared in the literature; in particular we permit both the fractal Laplacian and the potential to be Radon measures on the fractal. As a consequence of our results, we establish a Harnack inequality for solutions of these operators.


title: Nitrogen-containing Anions and Tholin Growth in Titan's Ionosphere: Implications for Cassini CAPS-ELS Observations

abstract: The Cassini Plasma Spectrometer (CAPS) Electron Spectrometer (ELS) instrument onboard Cassini revealed an unexpected abundance of negative ions above 950 km in Titan's ionosphere. \textit{In situ} measurements indicated the presence of negatively charged particles with mass-over-charge ratios up to 13,800 \textit{u/q}. At present, only a handful of anions have been characterized by photochemical models, consisting mainly of C$_n$H$^-$ carbon chain and C$_{n-1}$N$^-$ cyano compounds ($n=2-6$); their formation occurring essentially through proton abstraction from their parent neutral molecules. However, numerous other species have yet to be detected and identified. Considering the efficient anion growth leading to compounds of thousands of \textit{u/q}, it is necessary to better characterize the first light species. Here, we present new negative ion measurements with masses up to 200 \textit{u/q} obtained in an \ce{N2}:\ce{CH4} dusty plasma discharge reproducing analogous conditions to Titan's ionosphere. We perform a comparison with high altitude CAPS-ELS measurements near the top of Titan's ionosphere from the T18 encounter. The main observed peaks are in agreement with the observations. However, a number of other species (\textit{e.g.} \ce{CNN-}, \ce{CHNN-}) previously not considered suggests an abundance of N-bearing compounds, containing two or three nitrogen atoms, consistent with certain adjacent doubly-bonded nitrogen atoms found in tholins. These results suggest that an N-rich incorporation into tholins may follow mechanisms including anion chemistry, further highlighting the important role of negative ions in Titan's aerosol growth.


title: Can Intelligent Hyperparameter Selection Improve Resistance to Adversarial Examples?

abstract: Convolutional Neural Networks and Deep Learning classification systems in general have been shown to be vulnerable to attack by specially crafted data samples that appear to belong to one class but are instead classified as another, commonly known as adversarial examples. A variety of attack strategies have been proposed to craft these samples; however, there is no standard model that is used to compare the success of each type of attack. Furthermore, there is no literature currently available that evaluates how common hyperparameters and optimization strategies may impact a model's ability to resist these samples. This research bridges that lack of awareness and provides a means for the selection of training and model parameters in future research on evasion attacks against convolutional neural networks. The findings of this work indicate that the selection of model hyperparameters does impact the ability of a model to resist attack, although they alone cannot prevent the existence of adversarial examples.


title: The maximum advantage of quantum illumination

abstract: Discriminating between quantum states is a fundamental problem in quantum information protocols. The optimum approach saturates the Helstrom bound, which quantifies the unavoidable error probability of mistaking one state for another. Computing the error probability directly requires complete knowledge and diagonalization of the density matrices describing these states. Both of these fundamental requirements become impractically difficult to obtain as the dimension of the states grow large. In this article, we analyze quantum illumination as a quantum channel discrimination protocol and circumvent these issues by using the normalized Hilbert-Schmidt inner product as a measure of distinguishability. Using this measure, we show that the greatest advantage gained by quantum illumination over conventional illumination occurs when one uses a Bell state.


title: Observation of $B^0_{(s)} \to J/\psi p \overline{p}$ decays and precision measurements of the $B^0_{(s)}$ masses

abstract: The first observation of the decays $B^0_{(s)} \to J/\psi p \overline{p}$ is reported, using proton-proton collision data corresponding to an integrated luminosity of $5.2{\text{\,fb}}^{-1}$, collected with the \mbox{LHCb}\xspace detector. These decays are suppressed due to limited available phase space, as well as due to Okubo-Zweig-Iizuka or Cabibbo suppression. The measured branching fractions are \begin{align*} \mathcal{B}(B^0 \to J/\psi p \overline{p}) &= (4.51\pm 0.40\; \text{(stat)} \pm 0.44\; \text{(syst)}) \times 10^{-7}, \\ \mathcal{B}(B^0_s \to J/\psi p \overline{p}) &= (3.58\pm 0.19\; \text{(stat)} \pm 0.39\; \text{(syst)}) \times 10^{-6}. \end{align*} For the $B^0_s$ meson, the result is much higher than the expected value of $ {\cal O} (10^{-9})$. The small available phase space in these decays also allows for the most precise single measurement of both the $B^0$ mass as ${5279.74 \pm 0.30\; ({\rm stat})\pm 0.10\; ({\rm syst})}$ MeV, and the $B^0_s$ mass as ${5366.85 \pm 0.19\; ({\rm stat})\pm 0.13\; ({\rm syst})}$ MeV.


title: Photoionization Emission Models for the Cyg X-3 X-ray Spectrum

abstract: We present model fits to the X-ray line spectrum of the well known High Mass X-ray binary Cyg X-3. The primary observational dataset is a spectrum taken with the $Chandra$ X-ray Observatory High Energy Transmission Grating (HETG) in 2006, though we compare it to all the other observations of this source taken so far by this instrument. We show that the density must be $\geq 10^{12}$ cm$^{-3}$ in the region responsible for most of the emission. We discuss the influence of the dust scattering halo on the broad band spectrum and we argue that dust scattering and extinction is not the most likely origin for the narrow featureseen near the Si K edge. We identify the features of a wind in the profiles of the strong resonance lines and we show that the wind is more apparent in the lines from the lighter elements. We argue that this wind is most likely associated with the companion star. We show that the intensities of most lines can be fitted, crudely, by a single component photoionized model. However, the iron K lines do not fit with this model. We show that the iron K line variability as a function of orbital phase is different from the lower energy lines, which indicates that the lines arise in physically distinct regions. We discuss the interpretation of these results in the context of what is known about the system and similar sys


title: The Perils of Exploration under Competition: A Computational Modeling Approach

abstract: We empirically study the interplay between exploration and competition. Systems that learn from interactions with users often engage in exploration: making potentially suboptimal decisions in order to acquire new information for future decisions. However, when multiple systems are competing for the same market of users, exploration may hurt a system's reputation in the near term, with adverse competitive effects. In particular, a system may enter a "death spiral", when the short-term reputation cost decreases the number of users for the system to learn from, which degrades its performance relative to competition and further decreases its market share. We ask whether better exploration algorithms are incentivized under competition. We run extensive numerical experiments in a stylized duopoly model in which two firms deploy multi-armed bandit algorithms and compete for myopic users. We find that duopoly and monopoly tend to favor a primitive "greedy algorithm" that does not explore and leads to low consumer welfare, whereas a temporary monopoly (a duopoly with an early entrant) may incentivize better bandit algorithms and lead to higher consumer welfare. Our findings shed light on the first-mover advantage in the digital economy by exploring the role that data can play as a barrier to entry in online markets.


title: On 3d dipolar Bose-Einstein condensates involving quantum fluctuations and three-body interactions

abstract: We study the following nonlocal mixed order Gross-Pitaevskii equation $$i\,\partial_t \psi=-\frac{1}{2}\,\Delta \psi+V_{ext}\,\psi+\lambda_1\,|\psi|^2\,\psi+\lambda_2\,(K*|\psi|^2)\,\psi+\lambda_3\,|\psi|^{p-2}\,\psi,$$ where $K$ is the classical dipole-dipole interaction kernel, $\lambda_3>0$ and $p\in(4,6]$; the case $p=6$ being energy critical. For $p=5$ the equation is considered currently as the state-of-the-art model for describing the dynamics of dipolar Bose-Einstein condensates (Lee-Huang-Yang corrected dipolar GPE). We prove existence and nonexistence of standing waves in different parameter regimes; for $p\neq 6$ we prove global well-posedness and small data scattering.


title: Matroids and their Dressians

abstract: We study Dressians of matroids using the initial matroids of Dress and Wenzel. These correspond to cells in regular matroid subdivisions of matroid polytopes. An efficient algorithm for computing Dressians is presented, and its implementation is applied to a range of interesting matroids. We give counterexamples to a few plausible statements about matroid subdivisions.


title: Antipodal Hadwiger numbers of finite-dimensional Banach spaces

abstract: Let $X$ be a finite-dimensional Banach space; we introduce and investigate a natural generalization of the concepts of Hadwiger number $H(X)$ and strict Hadwiger number $H'(X)$. More precisely, we define the antipodal Hadwiger number $H_\alpha(X)$ as the largest cardinality of a subset $S \subseteq S_X$, such that $\forall x \neq y \in S \,\,\, \exists f \in B_{X^*}$ with \[1 \le f(x)-f(y) \,\,\, \textrm{and} \,\,\, f(y) \le f(z) \le f(x) \,\,\, \textrm{for} \,\,\, z \in S.\] The strict antipodal Hadwiger number $H'_\alpha(X)$ is defined analogously. We prove that $H'_\alpha(X)=4$ for every Minkowski plane and estimate (or in some cases compute) the numbers $H_\alpha(X)$ and $H'_\alpha(X)$, where $X=\ell_p^n, 1 < p \le +\infty$ and $n \ge 2$. We also show that the number $H'_\alpha(X)$ grows exponentially in $\dim X$.


title: Variability Abstraction and Refinement for Game-based Lifted Model Checking of full CTL (Extended Version)

abstract: Variability models allow effective building of many custom model variants for various configurations. Lifted model checking for a variability model is capable of verifying all its variants simultaneously in a single run by exploiting the similarities between the variants. The computational cost of lifted model checking still greatly depends on the number of variants (the size of configuration space), which is often huge. One of the most promising approaches to fighting the configuration space explosion problem in lifted model checking are variability abstractions. In this work, we define a novel game-based approach for variability-specific abstraction and refinement for lifted model checking of the full CTL, interpreted over 3-valued semantics. We propose a direct algorithm for solving a 3-valued (abstract) lifted model checking game. In case the result of model checking an abstract variability model is indefinite, we suggest a new notion of refinement, which eliminates indefinite results. This provides an iterative incremental variability-specific abstraction and refinement framework, where refinement is applied only where indefinite results exist and definite results from previous iterations are reused.


title: Fredholm determinant solutions of the Painlev\'e II hierarchy and gap probabilities of determinantal point processes

abstract: We study Fredholm determinants of a class of integral operators, whose kernels can be expressed as double contour integrals of a special type. Such Fredholm determinants appear in various random matrix and statistical physics models. We show that the logarithmic derivatives of the Fredholm determinants are directly related to solutions of the Painlev\'e II hierarchy. This confirms and generalizes a recent conjecture by Le Doussal, Majumdar, and Schehr. In addition, we obtain asymptotics at $\pm\infty$ for the Painlev\'e transcendents and large gap asymptotics for the corresponding point processes.


title: Low-Density Hydrodynamic Optical-Field-Ionized Plasma Channels Generated With An Axicon Lens

abstract: We demonstrate optical guiding of high-intensity laser pulses in long, low density hydrodynamic optical-field-ionized (HOFI) plasma channels. An axicon lens is used to generate HOFI plasma channels with on-axis electron densities as low as $n_e(0) = 1.5\times 10^{17}\, \mathrm{cm}^{-3}$ and matched spot sizes in the range $ 20 \mu \mathrm{m} \lesssim W_M \lesssim 40 \mu \mathrm{m}$. Control of these channel parameters via adjustment of the initial cell pressure and the delay after the arrival of the channel-forming pulse is demonstrated. For laser pulses with a peak axial intensity of $4 \times 10^{17}\, \mathrm{W\,cm}^{-2}$, highly reproducible, high-quality guiding over more than 14 Rayleigh ranges is achieved at a pulse repetition rate of 5 Hz, limited by the available channel-forming laser and vacuum pumping system. Plasma channels of this type would seem to be well suited to multi-GeV laser wakefield accelerators operating in the quasi-linear regime.


title: There are at most $2^{d+1}-2$ neighbourly simplices in dimension $d$

abstract: A combinatorial theorem on families of disjoint sub-boxes of a discrete cube, which implies that there at most $2^{d+1}-2$ neighbourly simplices in $\mathbb R^d$, is presented.


title: A theoretical model for realistic local climates

abstract: We write a nonlinear model that predicts the climate (temperature and humidity) on the surface of a small region on Earth, perform numerical investigations using the model, and compare the results to real climate on a variety of regions on Earth. It the parameters are chosen keeping into consideration the climatic K\"oppen zone to which the region belongs, the numerical model accurately reproduces the real climate. The model takes into account the doubly-periodic forcing of the solar radiation (annual and daily), the laws of irradiance, the fact that the Earth has land and oceans with different thermic inertia, and the humidity of the air due to evaporation. This enables us to reproduce remarkable features of Earth's climate such as lag of seasons, lag of noons, and asymmetric evolution of daily temperatures. The model can easily be adapted to planets with non-terrestrial astronomic parameters. We conclude this article with an investigation of an Earth with eccentricity higher than real.


title: Measurement of $B^+$, $B^0$ and $\Lambda_b^0$ production in $p\mkern 1mu\mathrm{Pb}$ collisions at $\sqrt{s_\mathrm{NN}}=8.16\,{\rm TeV}$

abstract: The production of $B^+$, $B^0$ and $\Lambda_b^0$ hadrons is studied in proton-lead collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}=8.16\,{\rm TeV}$ recorded with the LHCb detector at the LHC. The measurement uses a dataset corresponding to an integrated luminosity of $12.2\pm0.3\,\mathrm{nb}^{-1}$ for the case where the proton beam is projected into the LHCb detector (corresponding to measuring hadron production at positive rapidity) and $18.6\pm0.5\,\mathrm{nb}^{-1}$ for the lead beam projected into the LHCb detector (corresponding to measuring hadron production at negative rapidity). Nuclear effects are probed through double-differential cross-sections, forward-to-backward cross-section ratios and nuclear modification factors of the beauty hadrons. The double-differential cross-sections are measured as a function of the beauty-hadron transverse momentum and rapidity in the nucleon-nucleon centre-of-mass frame. Forward-to-backward cross-section ratios and nuclear modification factors indicate a significant nuclear suppression at positive rapidity. The ratio of $\Lambda_b^0$ over $B^0$ production cross-sections is reported and is consistent with the corresponding measurement in $pp$~collisions.


title: Dalton's and Amagat's laws fail in gas mixtures with shock propagation

abstract: As a shock wave propagates through a gas mixture, pressure, temperature, and density increase across the shock front. Rankine-Hugoniot (R-H) relations quantify these changes, correlating post-shock quantities with upstream conditions (pre-shock) and incident shock Mach number [1-5]. These equations describe a calorically perfect gas, but deliver a good approximation for real gases, provided the upstream conditions are well-characterized with a thermodynamic mixing model. Two classic thermodynamic models of gas mixtures are Dalton's law of partial pressures and Amagat's law of partial volumes [6]. Here we show that neither thermodynamic model can accurately predict the post-shock quantities of interest (temperature and pressure), on time scales much longer than those associated with the shock front passage, due to their implicit assumptions about behavior on the molecular level, including mixing reversibility. We found that in non-reacting binary mixtures of sulfur hexafluoride (SF6) and helium (He), kinetic molecular theory (KMT) can be used to quantify the discrepancies found between theoretical and experimental values for post-shock pressure and temperature. Our results demonstrate the complexity of analyzing shock wave interaction with two highly disparate gases, while also providing starting points for future theoretical and experimental work and validation of numerical simulations.


title: Exponentially-Modified Gaussian Mixture Model: Applications in Spectroscopy

abstract: We propose a novel exponentially-modified Gaussian (EMG) mixture residual model. The EMG mixture is well suited to model residuals that are contaminated by a distribution with positive support. This is in contrast to commonly used robust residual models, like the Huber loss or $\ell_1$, which assume a symmetric contaminating distribution and are otherwise asymptotically biased. We propose an expectation-maximization algorithm to optimize an arbitrary model with respect to the EMG mixture. We apply the approach to linear regression and probabilistic matrix factorization (PMF). We compare against other residual models, including quantile regression. Our numerical experiments demonstrate the strengths of the EMG mixture on both tasks. The PMF model arises from considering spectroscopic data. In particular, we demonstrate the effectiveness of PMF in conjunction with the EMG mixture model on synthetic data and two real-world applications: X-ray diffraction and Raman spectroscopy. We show how our approach is effective in inferring background signals and systematic errors in data arising from these experimental settings, dramatically outperforming existing approaches and revealing the data's physically meaningful components.


title: S-wave baryon nonleptonic decay amplitude in large-N_c chiral perturbation theory

abstract: The s-wave decay amplitude in the nonleptonic decay of baryons is analyzed within heavy baryon chiral perturbation theory in the large-N_c limit at one-loop order, where N_c is the number of color charges. Loop graphs with octet and decuplet intermediate states are systematically incorporated into the analysis and the effects of the decuplet-octet mass difference are accounted for. There are large-N_c cancellations between different one-loop graphs as a consequence of the large-N_c spin-flavor symmetry of QCD baryons. The predictions of large-N_c baryon chiral perturbation theory are in very good agreement both with the expectations from the 1/N_c expansion and with the experimental data.


title: Effective and Infinite-Rank Superrigidity in the Context of Representation Stability

abstract: We discuss certain effective improvements on superrigidity for $SL_n(\mathbb{Z})$ for finite $n>2$. Using these ideas we then use superrigidity to prove a representation stability theorem about pointwise finite dimensional $VIC(\mathbb{Z})$-modules, which itself can be viewed as a superrigidity theorem for $VIC(\mathbb{Z})$ and $GL_\infty(\mathbb{Z})$


title: Continuous Reachability for Unordered Data Petri nets is in PTime

abstract: Unordered data Petri nets (UDPN) are an extension of classical Petri nets with tokens that carry data from an infinite domain and where transitions may check equality and disequality of tokens. UDPN are well-structured, so the coverability and termination problems are decidable, but with higher complexity than for Petri nets. On the other hand, the problem of reachability for UDPN is surprisingly complex, and its decidability status remains open. In this paper, we consider the continuous reachability problem for UDPN, which can be seen as an over-approximation of the reachability problem. Our main result is a characterization of continuous reachability for UDPN and polynomial time algorithm for solving it. This is a consequence of a combinatorial argument, which shows that if continuous reachability holds then there exists a run using only polynomially many data values.


title: CrossQ: Batch Normalization in Deep Reinforcement Learning for Greater Sample Efficiency and Simplicity

abstract: Sample efficiency is a crucial problem in deep reinforcement learning. Recent algorithms, such as REDQ and DroQ, found a way to improve the sample efficiency by increasing the update-to-data (UTD) ratio to 20 gradient update steps on the critic per environment sample. However, this comes at the expense of a greatly increased computational cost. To reduce this computational burden, we introduce CrossQ: A lightweight algorithm for continuous control tasks that makes careful use of Batch Normalization and removes target networks to surpass the current state-of-the-art in sample efficiency while maintaining a low UTD ratio of 1. Notably, CrossQ does not rely on advanced bias-reduction schemes used in current methods. CrossQ's contributions are threefold: (1) it matches or surpasses current state-of-the-art methods in terms of sample efficiency, (2) it substantially reduces the computational cost compared to REDQ and DroQ, (3) it is easy to implement, requiring just a few lines of code on top of SAC.


title: A New Twist on the Electroclinic Critical Point: Type I and Type II Smectic $C^*$ Systems

abstract: This analysis of the electroclinic effect in ferroelectric liquid crystals with a first order Smectic-$A^*$--Smectic-$C^*$ (Sm-$A^*$--Sm-$C^*$) transition, shows they can be either Type I or Type II. In temperature--field parameter space Type I systems exhibit a macroscopically achiral (in which the Sm-$C^*$ helical superstructure is expelled) low-tilt (LT) Sm-$C$--high-tilt (HT) Sm-$C$ critical point, which terminates a LT Sm-$C$--HT Sm-$C$ first order boundary. This boundary extends to an achiral-chiral triple point where the achiral LT Sm-$C$ and HT Sm-$C$ phases coexist with the chiral Sm-$C^*$ phase. In Type II systems the critical point, triple point, and first order boundary are replaced by a Sm-$C^*$ region, between LT and HT achiral Sm-$C$ phases, at low and high fields respectively. When the field is ramped up, the Type II system displays a reentrant Sm-$C$--Sm-$C^*$-Sm-$C$ phase sequence. A discontinuity in the tilt of the optical axis at each of the two phase transitions means the Type II system is tristable. The system is Type I or Type II depending on the ratio of two length scales, one is the zero-field Sm-$C^*$ helical pitch, the other depends on the latent heat at the zero-field first order Sm-$A^*$--Sm-$C^*$ transition. A system could be experimentally tuned by varying enantiomeric excess, between Type I and Type II behavior. We also show this Type I vs Type II behavior is the Ising universality class analog of Type I vs Type II behavior in XY universality class systems. Lastly, we make a complete mapping of the phase boundaries in temperature--field--enantiomeric excess parameter space (not just near the critical point) which shows a variety of interesting features, including a multicritical point, tricritical points and a doubly reentrant Sm-$C$--Sm-$C^*$-Sm-$C$--Sm-$C^*$ phase sequence.


title: Learning for DC-OPF: Classifying active sets using neural nets

abstract: The optimal power flow is an optimization problem used in power systems operational planning to maximize economic efficiency while satisfying demand and maintaining safety margins. Due to uncertainty and variability in renewable energy generation and demand, the optimal solution needs to be updated in response to observed uncertainty realizations or near real-time forecast updates. To address the challenge of computing such frequent real-time updates to the optimal solution, recent literature has proposed the use of machine learning to learn the mapping between the uncertainty realization and the optimal solution. Further, learning the active set of constraints at optimality, as opposed to directly learning the optimal solution, has been shown to significantly simplify the machine learning task, and the learnt model can be used to predict optimal solutions in real-time. In this paper, we propose the use of classification algorithms to learn the mapping between the uncertainty realization and the active set of constraints at optimality, thus further enhancing the computational efficiency of the real-time prediction. We employ neural net classifiers for this task and demonstrate the excellent performance of this approach on a number of systems in the IEEE PES PGLib-OPF benchmark library.


title: Coupled nonlinear delay systems as deep convolutional neural networks

abstract: Neural networks are currently transforming the field of computer algorithms, yet their emulation on current computing substrates is highly inefficient. Reservoir computing was successfully implemented on a large variety of substrates and gave new insight in overcoming this implementation bottleneck. Despite its success, the approach lags behind the state of the art in deep learning. We therefore extend time-delay reservoirs to deep networks and demonstrate that these conceptually correspond to deep convolutional neural networks. Convolution is intrinsically realized on a substrate level by generic drive-response properties of dynamical systems. The resulting novelty is avoiding vector-matrix products between layers, which cause low efficiency in today's substrates. Compared to singleton time-delay reservoirs, our deep network achieves accuracy improvements by at least an order of magnitude in Mackey-Glass and Lorenz timeseries prediction.


title: Waveform-sequencing for scanning tunneling microscopy based pump-probe spectroscopy and pulsed-ESR

abstract: We present here the straightforward implementation of pump-probe methods into existing scanning tunneling microscopy (STM) systems. Our method uses the waveform-sequencing abilities of a standard arbitrary waveform generator (AWG) and a simple mechanical relay switch that either connects the regular STM control electronics or the AWG to the STM system. Our upgrade further enables pulsed-ESR excitation for advanced STM based spin-resonance experiments. We demonstrate the technical implementation, signal detection using a lock-in amplifier, and cross-correlation measurements of DC/DC and DC/RF pulses highlighting our ~5 ns time resolution, here limited by the speed of the available electronics. Our setup is highly versatile and can be extended to suit other needs of STM based investigations such as required in diverse mapping schemes or the coherent manipulation of qubits.


title: Partial Identification in Matching Models for the Marriage Market

abstract: We study partial identification of the preference parameters in the one-to-one matching model with perfectly transferable utilities. We do so without imposing parametric distributional assumptions on the unobserved heterogeneity and with data on one large market. We provide a tractable characterisation of the identified set under various classes of nonparametric distributional assumptions on the unobserved heterogeneity. Using our methodology, we re-examine some of the relevant questions in the empirical literature on the marriage market, which have been previously studied under the Logit assumption. Our results reveal that many findings in the aforementioned literature are primarily driven by such parametric restrictions.


title: GeoGAN: A Conditional GAN with Reconstruction and Style Loss to Generate Standard Layer of Maps from Satellite Images

abstract: Automatically generating maps from satellite images is an important task. There is a body of literature which tries to address this challenge. We created a more expansive survey of the task by experimenting with different models and adding new loss functions to improve results. We created a database of pairs of satellite images and the corresponding map of the area. Our model translates the satellite image to the corresponding standard layer map image using three main model architectures: (i) a conditional Generative Adversarial Network (GAN) which compresses the images down to a learned embedding, (ii) a generator which is trained as a normalizing flow (RealNVP) model, and (iii) a conditional GAN where the generator translates via a series of convolutions to the standard layer of a map and the discriminator input is the concatenation of the real/generated map and the satellite image. Model (iii) was by far the most promising of three models. To improve the results we also added a reconstruction loss and style transfer loss in addition to the GAN losses. The third model architecture produced the best quality of sampled images. In contrast to the other generative model where evaluation of the model is a challenging problem. since we have access to the real map for a given satellite image, we are able to assign a quantitative metric to the quality of the generated images in addition to inspecting them visually. While we are continuing to work on increasing the accuracy of the model, one challenge has been the coarse resolution of the data which upper-bounds the quality of the results of our model. Nevertheless, as will be seen in the results, the generated map is more accurate in the features it produces since the generator architecture demands a pixel-wise image translation/pixel-wise coloring. A video presentation summarizing this paper is available at: https://youtu.be/Ur0flOX-Ji0


title: Solving Complex Quadratic Systems with Full-Rank Random Matrices

abstract: We tackle the problem of recovering a complex signal $\boldsymbol x\in\mathbb{C}^n$ from quadratic measurements of the form $y_i=\boldsymbol x^*\boldsymbol A_i\boldsymbol x$, where $\boldsymbol A_i$ is a full-rank, complex random measurement matrix whose entries are generated from a rotation-invariant sub-Gaussian distribution. We formulate it as the minimization of a nonconvex loss. This problem is related to the well understood phase retrieval problem where the measurement matrix is a rank-1 positive semidefinite matrix. Here we study the general full-rank case which models a number of key applications such as molecular geometry recovery from distance distributions and compound measurements in phaseless diffractive imaging. Most prior works either address the rank-1 case or focus on real measurements. The several papers that address the full-rank complex case adopt the computationally-demanding semidefinite relaxation approach. In this paper we prove that the general class of problems with rotation-invariant sub-Gaussian measurement models can be efficiently solved with high probability via the standard framework comprising a spectral initialization followed by iterative Wirtinger flow updates on a nonconvex loss. Numerical experiments on simulated data corroborate our theoretical analysis.


title: Decentralized Privacy-preserving Timed Execution in Blockchain-based Smart Contract Platforms

abstract: In the age of Big Data, enabling task scheduling while protecting users' privacy is critical for various decentralized applications in blockchain-based smart contract platforms. Such a privacy-preserving task scheduler requires the task input data to be secretly maintained until a prescribed task execution time and be automatically recorded into the blockchain to enabling the execution of the task at the execution time, even if the user goes offline. While straight-forward centralized approaches provide a basic solution to the problem, unfortunately they are limited to a single point of trust and involve a single point of control. This paper presents decentralized techniques for supporting privacy-preserving task scheduling using smart contracts in Ethereum blockchain networks. We design a privacy-preserving task scheduling protocol that is managed by a manager smart contract. The protocol requires a user to schedule a task by deploying a proxy smart contract maintaining the non-sensitive information of the task while creating decentralized secret trust and selecting trustees from the network to maintain the sensitive information of the task. With security techniques including secret sharing and layered encryption as well as security deposit paid by trustees as economic deterrence, the protocol can protect the sensitive information against possible attacks including some trustees destroying the sensitive information (drop attack) or secretly releasing the sensitive information before the execution time (release-ahead attack). We demonstrate the attack-resilience of the proposed protocol through rigorous analysis.Our implementation and experimental evaluation on the Ethereum official test network demonstrate the low monetary cost and the low time overhead associated with the proposed approach.


title: Elastic and diffractive scattering in the LHC era

abstract: We review elastic and diffractive scattering of protons (called also "forward physics") with emphasis on the LHC data, especially those deviating from the expectations based on extrapolations from earlier measurements at the ISR, Fermilab and thus triggering searches for new ideas, models and theories. We list these new data and provide a brief introduction of available theoretical approaches, mainly those based on analyticity, crossing symmetry and unitarity, particularly the Regge-pole model realizing these concepts. Fits to the data are presented and tensions between theoretical predictions and the data that may indicate the way to further progress are in the focus of our paper.


title: Limits on Mode Coherence Due to a Non-static Convection Zone

abstract: The standard theory of pulsations deals with the frequencies and growth rates of infinitesimal perturbations in a stellar model. Modes which are calculated to be linearly driven should increase their amplitudes exponentially with time; the fact that nearly constant amplitudes are usually observed is evidence that nonlinear mechanisms inhibit the growth of finite amplitude pulsations. Models predict that the mass of DAV convection zones is very sensitive to temperature (i.e., $M_{\text{CZ}} \propto T_{\text{eff}}^{-90}$) leading to the possibility that even "small amplitude" pulsators may experience significant nonlinear effects. In particular, the outer turning point of finite-amplitude g-mode pulsations can vary with the local surface temperature, producing a reflected wave that is slightly out of phase with that required for a standing wave. This can lead to a lack of coherence of the mode and a reduction in its global amplitude. We compute the size of this effect for specific examples and discuss the results in the context of Kepler and K2 observations.


title: Dualizing Le Cam's method for functional estimation, with applications to estimating the unseens

abstract: Le Cam's method (or the two-point method) is a commonly used tool for obtaining statistical lower bound and especially popular for functional estimation problems. This work aims to explain and give conditions for the tightness of Le Cam's lower bound in functional estimation from the perspective of convex duality. Under a variety of settings it is shown that the maximization problem that searches for the best two-point lower bound, upon dualizing, becomes a minimization problem that optimizes the bias-variance tradeoff among a family of estimators. For estimating linear functionals of a distribution our work strengthens prior results of Donoho-Liu \cite{DL91} (for quadratic loss) by dropping the H\"olderian assumption on the modulus of continuity. For exponential families our results extend those of Juditsky-Nemirovski \cite{JN09} by characterizing the minimax risk for the quadratic loss under weaker assumptions on the exponential family. We also provide an extension to the high-dimensional setting for estimating separable functionals. Notably, coupled with tools from complex analysis, this method is particularly effective for characterizing the ``elbow effect'' -- the phase transition from parametric to nonparametric rates. As the main application we derive sharp minimax rates in the Distinct elements problem (given a fraction $p$ of colored balls from an urn containing $d$ balls, the optimal error of estimating the number of distinct colors is $\tilde \Theta(d^{-\frac{1}{2}\min\{\frac{p}{1-p},1\}})$) and the Fisher's species problem (given $n$ iid observations from an unknown distribution, the optimal prediction error of the number of unseen symbols in the next (unobserved) $r \cdot n$ observations is $\tilde \Theta(n^{-\min\{\frac{1}{r+1},\frac{1}{2}\}})$).


title: The Nowicki Conjecture for free metabelian Lie algebras

abstract: Let $K[X_d]=K[x_1,\ldots,x_d]$ be the polynomial algebra in $d$ variables over a field $K$ of characteristic 0. The classical theorem of Weitzenb\"ock from 1932 states that for linear locally nilpotent derivations $\delta$ (known as Weitzenb\"ock derivations) the algebra of constants $K[X_{d}]^{\delta}$ is finitely generated. When the Weitzenb\"ock derivation $\delta$ acts on the polynomial algebra $K[X_d,Y_d]$ in $2d$ variables by $\delta(y_i)=x_i$, $\delta(x_i)=0$, $i=1,\ldots,d$, Nowicki conjectured that $K[X_d,Y_d]^{\delta}$ is generated by $X_d$ and $x_iy_j-y_ix_j$ for all $1\leq i<j\leq d$. There are several proofs based on different ideas confirming this conjecture. Considering arbitrary Weitzenb\"ock derivations of the free $d$-generated metabelian Lie algebra $F_d$, with few trivial exceptions, the algebra $F_d^{\delta}$ is not finitely generated. However, the vector subspace $(F_d')^{\delta}$ of the commutator ideal $F_d'$ of $F_d$ is finitely generated as a $K[X_d]^{\delta}$-module. In this paper we study an analogue of the Nowicki conjecture in the Lie algebra setting and give an explicit set of generators of the $K[X_d,Y_d]^{\delta}$-module $(F_{2d}')^{\delta}$.


title: Nuclear Transparency in Monte Carlo Neutrino Event Generators

abstract: Hadron cascade model is an essential part of Monte Carlo neutrino event generators that governs final state interactions of knocked-out nucleons and produced pions. It is shown that such model enriched with physically motivated modifications of nucleon-nucleon cross section and incorporation of nuclear correlation effects is able to reproduce experimental nuclear transparency data. Uncertainty of nucleon final state interactions effects is estimated and applied to recent neutrino-nucleus cross section measurements that include an outgoing proton in the experimental signal. Conclusions are drawn on a perspective of identification of events that originate from two-body current mechanism.


title: Superposition principle and schemes for Measure Differential Equations

abstract: Measure Differential Equations (MDE) describe the evolution of probability measures driven by probability velocity fields, i.e. probability measures on the tangent bundle. They are, on one side, a measure-theoretic generalization of ordinary differential equations; on the other side, they allow to describe concentration and diffusion phenomena typical of kinetic equations. In this paper, we analyze some properties of this class of differential equations, especially highlighting their link with nonlocal continuity equations. We prove a representation result in the spirit of the Superposition Principle by Ambrosio-Gigli-Savar\'e, and we provide alternative schemes converging to a solution of the MDE, with a particular view to uniqueness/non-uniqueness phenomena.


title: Calculation and verification of neutron irradiation damage with differential cross sections

abstract: The Displacement per Atom (DPA) rate is conventionally computed with DPA cross sections in reactor applications. The method of direct calculation with energy-angular distributions given in the Center of Mass (CM) frame is proposed and recommended in the present work. The methods for refining and verifying the calculations of DPA cross sections are proposed: (i) Gauss-Legendre-Quadrature-based Piecewise Integration (GLQPI) for ensuring the numeric convergence of integral over emission angle due to the discontinuity of integrand; (ii) verification of the convergence for trapezoidal integration over the secondary energy; (iii) interpolation of double-differential cross sections. For 56Fe of JEFF-3.1.1, the current numeric integration over emission angle is shown not convergent, whereas the direct trapezoidal over the secondary energy and the direct interpolation of energy-angle-integrated damage are shown accurate. On the other hand, it is shown that the DPA cross sections are overestimated if isotropic angular distributions are assumed. However, the DPA cross section is not sensitive to the high-order Legendre polynomials because the former is an angle-integrated quantity. Numerical results of neutron elastic scattering show that 2 orders of Legendre polynomials can give the DPA rates of 56Fe within 0.5% overestimation for fission reactors, while 4 orders are required for fusion reactors. For neutron inelastic scatterings-induced DPA, the first order Legendre polynomial is sufficient for both fission and fusion reactors.


title: A Revision of the Bernoulli Equation as a Controller of the Fick's Diffusion Equation in Drug Delivery Modeling

abstract: Mathematical equations can be used as effectual tools in drug delivery systems modeling and are also highly helpful to have a theoretical understanding of controlled drug release and diffusion mechanisms. In this study we aim to present a mathematical combination between the Bernoulli equation and the Fick's equation as a diffusion controller in drug delivery systems. For this propose we have revised the Bernoulli equation as an additional, controller and complementary method of the Fick's diffusion equation to detect the optimal delivery direction to control the diffusion divergence of the drug carrier in vascular systems during the transportation process in biological tissues. Therefore, by utilizing the Bernoulli equation we could determine the real direction by the route function f.


title: The Shapley Taylor Interaction Index

abstract: The attribution problem, that is the problem of attributing a model's prediction to its base features, is well-studied. We extend the notion of attribution to also apply to feature interactions. The Shapley value is a commonly used method to attribute a model's prediction to its base features. We propose a generalization of the Shapley value called Shapley-Taylor index that attributes the model's prediction to interactions of subsets of features up to some size k. The method is analogous to how the truncated Taylor Series decomposes the function value at a certain point using its derivatives at a different point. In fact, we show that the Shapley Taylor index is equal to the Taylor Series of the multilinear extension of the set-theoretic behavior of the model. We axiomatize this method using the standard Shapley axioms -- linearity, dummy, symmetry and efficiency -- and an additional axiom that we call the interaction distribution axiom. This new axiom explicitly characterizes how interactions are distributed for a class of functions that model pure interaction. We contrast the Shapley-Taylor index against the previously proposed Shapley Interaction index (cf. [9]) from the cooperative game theory literature. We also apply the Shapley Taylor index to three models and identify interesting qualitative insights.


title: Decentralized Release of Self-emerging Data using Smart Contracts

abstract: In the age of Big Data, releasing protected sensitive data at a future point in time is critical for various applications. Such self-emerging data release requires the data to be protected until a prescribed data release time and be automatically released to the recipient at the release time, even if the data sender goes offline. While straight-forward centralized approaches provide a basic solution to the problem, unfortunately they are limited to a single point of trust and involve a single point of control. This paper presents decentralized techniques for supporting self-emerging data using smart contracts in Ethereum blockchain networks. We design a credible and enforceable smart contract for supporting self-emerging data release. The smart contract employs a set of Ethereum peers to jointly follow the proposed timed-release service protocol allowing the participating peers to earn the remuneration paid by the service users.We model the problem as an extensive-form game with imperfect information to protect against possible adversarial attacks including some peers destroying the private data (drop attack) or secretly releasing the private data before the release time (release-ahead attack). We demonstrate the efficacy and attack-resilience of the proposed techniques through rigorous analysis and experimental evaluation. Our implementation and experimental evaluation on the Ethereum official test network demonstrate the low monetary cost and the low time overhead associated with the proposed approach and validate its guaranteed security properties.


title: Quick and Easy Time Series Generation with Established Image-based GANs

abstract: In the recent years Generative Adversarial Networks (GANs) have demonstrated significant progress in generating authentic looking data. In this work we introduce our simple method to exploit the advancements in well established image-based GANs to synthesise single channel time series data. We implement Wasserstein GANs (WGANs) with gradient penalty due to their stability in training to synthesise three different types of data; sinusoidal data, photoplethysmograph (PPG) data and electrocardiograph (ECG) data. The length of the returned time series data is limited only by the image resolution, we use an image size of 64x64 pixels which yields 4096 data points. We present both visual and quantitative evidence that our novel method can successfully generate time series data using image-based GANs.


title: WaveletAE: A Wavelet-enhanced Autoencoder for Wind Turbine Blade Icing Detection

abstract: Wind power, as an alternative to burning fossil fuels, is abundant and inexhaustible. To fully utilize wind power, wind farms are usually located in areas of high altitude and facing serious ice conditions, which can lead to serious consequences. Quick detection of blade ice accretion is crucial for the maintenance of wind farms. Unlike traditional methods of installing expensive physical detectors on wind blades, data-driven approaches are increasingly popular for inspecting the wind turbine failures. In this work, we propose a wavelet enhanced autoencoder model (WaveletAE) to identify wind turbine dysfunction by analyzing the multivariate time series monitored by the SCADA system. WaveletAE is enhanced with wavelet detail coefficients to enforce the autoencoder to capture information from multiple scales, and the CNN-LSTM architecture is applied to learn channel-wise and temporal-wise relations. The empirical study shows that the proposed model outperforms other state-of-the-art time series anomaly detection methods for real-world blade icing detection.


title: Counting square-tiled surfaces with prescribed real and imaginary foliations and connections to Mirzakhani's asymptotics for simple closed hyperbolic geodesics

abstract: We show that the number of square-tiled surfaces of genus $g$, with $n$ marked points, with one or both of its horizontal and vertical foliations belonging to fixed mapping class group orbits, and having at most $L$ squares, is asymptotic to $L^{6g-6+2n}$ times a product of constants appearing in Mirzakhani's count of simple closed hyperbolic geodesics. Many of the results in this paper reflect recent discoveries of Delecroix, Goujard, Zograf, and Zorich, but the approach considered here is very different from theirs. We follow conceptual and geometric methods inspired by Mirzakhani's work.


title: Classification with unknown class-conditional label noise on non-compact feature spaces

abstract: We investigate the problem of classification in the presence of unknown class-conditional label noise in which the labels observed by the learner have been corrupted with some unknown class dependent probability. In order to obtain finite sample rates, previous approaches to classification with unknown class-conditional label noise have required that the regression function is close to its extrema on sets of large measure. We shall consider this problem in the setting of non-compact metric spaces, where the regression function need not attain its extrema. In this setting we determine the minimax optimal learning rates (up to logarithmic factors). The rate displays interesting threshold behaviour: When the regression function approaches its extrema at a sufficient rate, the optimal learning rates are of the same order as those obtained in the label-noise free setting. If the regression function approaches its extrema more gradually then classification performance necessarily degrades. In addition, we present an adaptive algorithm which attains these rates without prior knowledge of either the distributional parameters or the local density. This identifies for the first time a scenario in which finite sample rates are achievable in the label noise setting, but they differ from the optimal rates without label noise.


title: Laser cooling of molecules

abstract: Recently, laser cooling methods have been extended from atoms to molecules. The complex rotational and vibrational energy level structure of molecules makes laser cooling difficult, but these difficulties have been overcome and molecules have now been cooled to a few microkelvin and trapped for several seconds. This opens many possibilities for applications in quantum science and technology, controlled chemistry, and tests of fundamental physics. This article explains how molecules can be decelerated, cooled and trapped using laser light, reviews the progress made in recent years, and outlines some future applications.


title: Environmentally-friendly GR(1) Synthesis

abstract: Many problems in reactive synthesis are stated using two formulas ---an environment assumption and a system guarantee--- and ask for an implementation that satisfies the guarantee in environments that satisfy their assumption. Reactive synthesis tools often produce strategies that formally satisfy such specifications by actively preventing an environment assumption from holding. While formally correct, such strategies do not capture the intention of the designer. We introduce an additional requirement in reactive synthesis, non-conflictingness, which asks that a system strategy should always allow the environment to fulfill its liveness requirements. We give an algorithm for solving GR(1) synthesis that produces non-conflicting strategies. Our algorithm is given by a 4-nested fixed point in the $\mu$-calculus, in contrast to the usual 3-nested fixed point for GR(1). Our algorithm ensures that, in every environment that satisfies its assumptions on its own, traces of the resulting implementation satisfy both the assumptions and the guarantees. In addition, the asymptotic complexity of our algorithm is the same as that of the usual GR(1) solution. We have implemented our algorithm and show how its performance compares to the usual GR(1) synthesis algorithm.


title: Using Key Player Analysis as a Method for Examining the Role of Community Animators in Technology Adoption

abstract: This paper examines the role of community animators in technology adoption. Community animators are individuals that actively build social networks and broker ties between nodes in those networks. The present study observes technology adoption patterns through data collected from a mobile application at a local arts festival. A social network was constructed through photo-sharing and interaction within the app. Given this data, we propose the use of key player analysis to identify community animators. In addition, we use a graph invariant (i.e., fragmentation in the network) to describe the role and impact of key players on the full network of interactions. Our results contribute to literature on technology adoption in usability studies by proposing a method to quantify and identify the theoretical concept of community animators. We further analyze the types of community animators to be found in early adoption of technology: the early adopters themselves, and the initiating developers.


title: Super-Resolution Simulation for Real-Time Prediction of Urban Micrometeorology

abstract: We propose a super-resolution (SR) simulation system that consists of a physics-based meteorological simulation and an SR method based on a deep convolutional neural network (CNN). The CNN is trained using pairs of high-resolution (HR) and low-resolution (LR) images created from meteorological simulation results for different resolutions so that it can map LR simulation images to HR ones. The proposed SR simulation system, which performs LR simulations, can provide HR prediction results in much shorter operating cycles than those required for corresponding HR simulation prediction system. We apply the SR simulation system to urban micrometeorology, which is strongly affected by buildings and human activity. Urban micrometeorology simulations that need to resolve urban buildings are computationally costly and thus cannot be used for operational real-time predictions even when run on supercomputers. We performed HR micrometeorology simulations on a supercomputer to obtain datasets for training the CNN in the SR method. It is shown that the proposed SR method can be used with a spatial scaling factor of 4 and that it outperforms conventional interpolation methods by a large margin. It is also shown that the proposed SR simulation system has the potential to be used for operational urban micrometeorology predictions.


title: Verifiably Safe Off-Model Reinforcement Learning

abstract: The desire to use reinforcement learning in safety-critical settings has inspired a recent interest in formal methods for learning algorithms. Existing formal methods for learning and optimization primarily consider the problem of constrained learning or constrained optimization. Given a single correct model and associated safety constraint, these approaches guarantee efficient learning while provably avoiding behaviors outside the safety constraint. Acting well given an accurate environmental model is an important pre-requisite for safe learning, but is ultimately insufficient for systems that operate in complex heterogeneous environments. This paper introduces verification-preserving model updates, the first approach toward obtaining formal safety guarantees for reinforcement learning in settings where multiple environmental models must be taken into account. Through a combination of design-time model updates and runtime model falsification, we provide a first approach toward obtaining formal safety proofs for autonomous systems acting in heterogeneous environments.


title: Quantum Measurements and Contextuality

abstract: In quantum physics the term `contextual' can be used in more than one way. One usage, here called `Bell contextual' since the idea goes back to Bell, is that if $A$, $B$ and $C$ are three quantum observables, with $A$ compatible (i.e., commuting) with $B$ and also with $C$, whereas $B$ and $C$ are incompatible, a measurement of $A$ might yield a different result (indicating that quantum mechanics is contextual) depending upon whether $A$ is measured along with $B$ (the $\{A,B\}$ context) or with $C$ (the $\{A,C\}$ context). An analysis of what projective quantum measurements measure shows that quantum theory is Bell noncontextual: the outcome of a particular $A$ measurement when $A$ is measured along with $B$ would have been exactly the same if $A$ had, instead, been measured along with $C$. A different definition, here called `globally (non)contextual' refers to whether or not there is ('noncontextual') or is not ('contextual') a single joint probability distribution that simultaneously assigns probabilities in a consistent manner to the outcomes of measurements of a certain collection of observables, not all of which are compatible. A simple example shows that such a joint probability distribution can exist even in a situation where the measurement probabilities cannot refer to properties of a quantum system, and hence lack physical significance, even though mathematically well-defined. It is noted that the quantum sample space, a projective decomposition of the identity, required for interpreting measurements of incompatible properties in different runs of an experiment using different types of apparatus has a tensor product structure, a fact sometimes overlooked.


title: A quantum compiler for qudits of prime dimension greater than 3

abstract: Prevailing proposals for the first generation of quantum computers make use of 2-level systems, or qubits, as the fundamental unit of quantum information. However, recent innovations in quantum error correction and magic state distillation protocols demonstrate that there are advantages of using d-level quantum systems, known as \emph{qudits}, over the qubit analogues. When designing a quantum architecture, it is crucial to consider protocols for compilation, the optimal conversion of high-level instructions used by programmers into low-level instructions interpreted by the machine. In this work, we present a general purpose automated compiler for multiqudit exact synthesis based on previous work on qubits that uses an algebraic representation of quantum circuits called phase polynomials. We assume Clifford gates are low-cost and aim to minimise the number of M gates in a Clifford+M circuit, where M is the qudit analog for the qubit T or pi/8 phase gate. A surprising result that showcases our compiler's capabilities is that we found a unitary implementation of the CCZ or Toffoli gate that uses 4 M gates, which compares to 7 T gates for the qubit analogue.


title: Distributed Processes and Scalability in Sub-networks of Large-Scale Networks

abstract: Performance of standard processes over large distributed networks typically scales with the size of the network. For example, in planar topologies where nodes communicate with their natural neighbors, the scaling factor is $O(n)$, where $n$ is the number of nodes. As the size of the network increases, this makes global convergence over the entire network less practical. On the other hand, for several applications, such as load balancing or detection of local events, global convergence may not be necessary or even relevant. We introduce simple distributed iterative processes which limit the scope of the computational task to a smaller subnetwork of the entire network. This is achieved using one additional local parameter which controls the size of the subnetwork. We establish termination and convergence rate of such processes in theory, in experiment, in comparison to the well understood behavior of Markov processes, and for a variety of network topologies and initial conditions.


title: Calculating temperature-dependent properties of Nd$_2$Fe$_{14}$B permanent magnets by atomistic spin model simulations

abstract: Temperature-dependent magnetic properties of Nd$_2$Fe$_{14}$B permanent magnets, i.e., saturation magnetization $M_\text{s}(T)$, effective magnetic anisotropy constants $K_i^\text{eff}(T)$ ($i=1,2,3$), domain wall width $\delta_w(T)$, and exchange stiffness constant $A_\text{e}(T)$, are calculated by using \textit{ab-initio} informed atomistic spin model simulations. We construct the atomistic spin model Hamiltonian for Nd$_2$Fe$_{14}$B by using the Heisenberg exchange of Fe$-$Fe and Fe$-$Nd atomic pairs, the uniaxial single-ion anisotropy of Fe atoms, and the crystal-field energy of Nd ions which is approximately expanded into an energy formula featured by second, fourth, and sixth-order phenomenological anisotropy constants. After applying a temperature rescaling strategy, we show that the calculated Curie temperature, spin-reorientation phenomenon, $M_\text{s}(T)$, $\delta_w(T)$, and $K_i^\text{eff}(T)$ agree well with the experimental results. $A_\text{e}(T)$ is estimated through a general continuum description of the domain wall profile by mapping atomistic magnetic moments to the macroscopic magnetization. $A_\text{e}$ is found to decrease more slowly than $K_1^\text{eff}$ with increasing temperature, and approximately scale with normalized magnetization as $A_\text{e}(T) \sim m^{1.2}$. This work provokes a scale bridge between \textit{ab-initio} calculations and temperature-dependent micromagnetic simulations of Nd-Fe-B permanent magnets.


title: Mapping the Interstellar Reddening and Extinction towards Baade's Window Using Minimum Light Colors of ab-type RR Lyrae Stars. Revelations from the De-reddened Color-Magnitude Diagrams

abstract: We have obtained repeated images of 6 fields towards the Galactic bulge in 5 passbands (u, g, r, i, z) with the DECam imager on the Blanco 4m telescope at CTIO. From over 1.6 billion individual photometric measurements in the field centered on Baade's window, we have detected 4877 putative variable stars. 474 of these have been confirmed as fundamental mode RR Lyrae stars, whose colors at minimum light yield line-of-sight reddening determinations as well as a reddening law towards the Galactic Bulge which differs significantly from the standard R_V = 3.1 formulation. Assuming that the stellar mix is invariant over the 3 square-degree field, we are able to derive a line-of-sight reddening map with sub-arcminute resolution, enabling us to obtain de-reddened and extinction corrected color-magnitude diagrams (CMD's) of this bulge field using up to 2.5 million well-measured stars. The corrected CMD's show unprecedented detail and expose sparsely populated sequences: e.g., delineation of the very wide red giant branch, structure within the red giant clump, the full extent of the horizontal branch, and a surprising bright feature which is likely due to stars with ages younger than 1 Gyr. We use the RR Lyrae stars to trace the spatial structure of the ancient stars, and find an exponential decline in density with Galactocentric distance. We discuss ways in which our data products can be used to explore the age and metallicity properties of the bulge, and how our larger list of all variables is useful for learning to interpret future LSST alerts.


title: Finding Nearest Neighbors in graphs locally

abstract: Many distributed learning techniques have been motivated by the increasing size of datasets and their inability to fit into main memory on a single machine. We propose an algorithm that finds the nearest neighbor in a graph locally without the need of visiting the whole graph. Our algorithm is distributed which further encourage scalability. We prove the convergence of the algorithm


title: High-speed domain wall racetracks in a magnetic insulator

abstract: Recent reports of current-induced switching of ferrimagnetic oxides coupled to a heavy metal layer have opened realistic prospects for implementing magnetic insulators into electrically addressable spintronic devices. However, key aspects such as the configuration and dynamics of magnetic domain walls driven by electrical currents in insulating oxides remain unexplored. Here, we investigate the internal structure of the domain walls in Tm3Fe5O12 (TmIG) and TmIG/Pt bilayers and demonstrate their efficient manipulation by spin-orbit torques with velocities of up to 400 m s$^{-1}$ and minimal current threshold for domain wall flow of 5 x 10$^{6}$ A cm$^{-2}$. Domain wall racetracks embedded in TmIG are defined by the deposition of Pt current lines, which allow us to control the domain propagation and magnetization switching in selected regions of an extended magnetic layer. Scanning nitrogen-vacancy magnetometry reveals that the domain walls of thin TmIG films are N\'eel walls with left-handed chirality, with the domain wall magnetization rotating towards an intermediate N\'eel-Bloch configuration upon deposition of Pt. These results indicate the presence of a sizable interfacial Dzyaloshinskii-Moriya interaction in TmIG, which leads to novel possibilities to control the formation of chiral spin textures in magnetic insulators. Ultimately, domain wall racetracks provide an efficient scheme to pattern the magnetic landscape of TmIG in a fast and reversible way


title: Achieving Large Sum Rate and Good Fairness in MISO Broadcast Communication

abstract: A tradeoff between sum rate and fairness for MISO broadcast communication employing dirty paper coding or zero-forcing dirty paper coding at physical layer is investigated in this paper. The tradeoff is based on a new design objective termed "tri-stage" approach as well as a new L1-based fairness measure that is much more robust than the well-known Jain's index for comparing fairness levels achieved by various design objectives at a much finer resolution in high SNR regime. The newly proposed tri-stage design also introduces a new concept of statistical power allocation that randomly allocates powers to users based on an optimal probability distribution derived from the tradeoff between sum rate and fairness. Simulation results show that the proposed approach can simultaneously achieve a larger sum rate and better fairness than the reputable proportional fairness criterion. A performance upper bound is also given in the paper to show that the excellent performance of the proposed approach at moderate and high SNR regimes as well as some potential for further improvement in low SNR regime.


title: Suitability and robustness of triangular nanostructured targets for proton acceleration

abstract: Ion acceleration in the MeV range can be routinely achieved with table-top laser technology. One of the current challenges is to improve the energy coupling from the laser to the proton beam without increasing the laser peak power. Introducing nanostructures at the front target surface was shown to be beneficial for an efficient transfer of energy to the electrons. In this manuscript, we study by using full-scale three-dimensional particle-in-cell simulations and finite laser pulses, the process when a proposed optimal target with triangular nanostructure (previously found to allow 97% laser energy absorption) is used . We demonstrate that the absorbed laser energy does not depend on the dimensionality in the range of parameters presented. We also present an analytical model for laser absorption that includes deviations from the ideal conditions. This is supported by a numerical parameter study that establishes the tolerance with respect to the nanostructure size, use of different ion species, existence of preplasma, etc. We found that altering the target thickness or using different ions does not affect the absorption, but it does affect the energy redistribution among the different plasma species. The optimal configuration ($h = 1~\lambda,~ w = 0.7~ \lambda$) is robust with respect to the target fabrication errors. However, high contrast laser pulses are required, because a pre-plasma layer with a thickness on the order of 0.5 lambda is enough to lower the laser absorption by more than a 10% in a non-optimal scenario.


title: Quantum Simulation of Resonant Transitions for Solving the Eigen-problem of an Effective Water Hamiltonian

abstract: It is difficult to calculate the energy levels and eigenstates of a large physical system on a classical computer because of the exponentially growing size of the Hilbert space. In this work, we experimentally demonstrate a quantum algorithm which could solve this problem via simulated resonant transitions. Using a four-qubit quantum simulator in which two qubits are used as ancillas for control and measurement, we obtain the energy spectrum of a 2-qubit low-energy effective Hamiltonian of the water molecule. The simulated transitions allow the state of the quantum simulator to transform and access large regions of the Hilbert space, including states that have no overlap with the initial state. Furthermore, we make use of this algorithm to efficiently prepare specific eigenstates on the simulator according to the measured eigenenergies.


title: Shear rate projection schemes for non-Newtonian fluids

abstract: The operator splitting approach applied to the Navier-Stokes equations, gave rise to various numerical methods for the simulations of the dynamics of fluids. The separate work of Chorin and Temam on this subject gave birth to the so-called projection methods. The basic projection schemes, either the incremental or non-incremental variant (see [1]) induces an artificial Neumann boundary condition on the pressure. By getting rid of this boundary condition on the pressure, the so-call rotational incremental pressure-correction scheme as proposed by Timmermans et al. [2] for Newtonian fluids with constant viscosity gives a consistent equation for the pressure. In this work we propose a family of projection methods for generalized Newtonian fluids based on an extension of the rotational projection scheme. Called shear rate projections, these methods produces consistent pressure when applied to generalized Newtonian fluids. Accuracy of the methods will be illustrated using a manufactured solution. Numerical experiments for the flow past a cylinder, with a Carreau rheological model, will also be presented.


title: Active Perception in Adversarial Scenarios using Maximum Entropy Deep Reinforcement Learning

abstract: We pose an active perception problem where an autonomous agent actively interacts with a second agent with potentially adversarial behaviors. Given the uncertainty in the intent of the other agent, the objective is to collect further evidence to help discriminate potential threats. The main technical challenges are the partial observability of the agent intent, the adversary modeling, and the corresponding uncertainty modeling. Note that an adversary agent may act to mislead the autonomous agent by using a deceptive strategy that is learned from past experiences. We propose an approach that combines belief space planning, generative adversary modeling, and maximum entropy reinforcement learning to obtain a stochastic belief space policy. By accounting for various adversarial behaviors in the simulation framework and minimizing the predictability of the autonomous agent's action, the resulting policy is more robust to unmodeled adversarial strategies. This improved robustness is empirically shown against an adversary that adapts to and exploits the autonomous agent's policy when compared with a standard Chance-Constraint Partially Observable Markov Decision Process robust approach.


title: Degree of irrationality of very general abelian surfaces

abstract: The degree of irrationality of a projective variety $X$ is defined to be the smallest degree rational dominant map to a projective space of the same dimension. For abelian surfaces, Yoshihara computed this invariant in specific cases, while Stapleton gave a sublinear upper bound for very general polarized abelian surfaces $(A, L)$ of degree $d$. Somewhat surprisingly, we show that the degree of irrationality of a very general polarized abelian surface is uniformly bounded above by $4$, independently of the degree of the polarization. This result disproves part of a conjecture of Bastianelli, De Poi, Ein, Lazarsfeld and Ullery.