title: The Art of DNA Strings: Sixteen Years of DNA Coding Theory

abstract: The idea of computing with DNA was given by Tom Head in 1987, however in 1994 in a seminal paper, the actual successful experiment for DNA computing was performed by Adleman. The heart of the DNA computing is the DNA hybridization, however, it is also the source of errors. Thus the success of the DNA computing depends on the error control techniques. The classical coding theory techniques have provided foundation for the current information and communication technology (ICT). Thus it is natural to expect that coding theory will be the foundational subject for the DNA computing paradigm. For the successful experiments with DNA computing usually we design DNA strings which are sufficiently dissimilar. This leads to the construction of a large set of DNA strings which satisfy certain combinatorial and thermodynamic constraints. Over the last 16 years, many approaches such as combinatorial, algebraic, computational have been used to construct such DNA strings. In this work, we survey this interesting area of DNA coding theory by providing key ideas of the area and current known results.


title: Automated 5-year Mortality Prediction using Deep Learning and Radiomics Features from Chest Computed Tomography

abstract: We propose new methods for the prediction of 5-year mortality in elderly individuals using chest computed tomography (CT). The methods consist of a classifier that performs this prediction using a set of features extracted from the CT image and segmentation maps of multiple anatomic structures. We explore two approaches: 1) a unified framework based on deep learning, where features and classifier are automatically learned in a single optimisation process; and 2) a multi-stage framework based on the design and selection/extraction of hand-crafted radiomics features, followed by the classifier learning process. Experimental results, based on a dataset of 48 annotated chest CTs, show that the deep learning model produces a mean 5-year mortality prediction accuracy of 68.5%, while radiomics produces a mean accuracy that varies between 56% to 66% (depending on the feature selection/extraction method and classifier). The successful development of the proposed models has the potential to make a profound impact in preventive and personalised healthcare.


title: Well-posedness for mean-field evolutions arising in superconductivity

abstract: We establish the existence of a global solution for a new family of fluid-like equations, which are obtained in a joint work with Serfaty in certain regimes as the mean-field evolution of the supercurrent density in a (2D section of a) type-II superconductor with pinning and with imposed electric current. We also consider general vortex-sheet initial data, and investigate the uniqueness and regularity properties of the solution. For some choice of parameters, the equation under investigation coincides with the so-called lake equation from 2D shallow water fluid dynamics, and our analysis then leads to a new existence result for rough initial data.


title: Does Your DNS Recursion Really Time Out as Intended? A Timeout Vulnerability of DNS Recursive Servers

abstract: Parallelization is featured by DNS recursive servers to do time-consuming recursions on behalf on clients. As common DNS configurations, recursive servers should allow a reasonable timeout for each recursion which may take as long as several seconds. However, it is proposed in this paper that recursion parallelization may be exploited by attackers to compromise the recursion timeout mechanism for the purpose of DoS or DDoS attacks. Attackers can have recursive servers drop early existing recursions in service by saturating recursion parallelization. The key of the proposed attack model is to reliably prolong service times for any attacking queries. As means of prolong service times, serval techniques are proposed to effectively avoiding cache hit and prolonging overall latency of external DNS lookups respectively. The impacts of saturated recursion parallelization on timeout are analytically provided. The testing on BIND servers demonstrates that with carefully crafted queries, an attacker can use a low or moderate level of query load to successfully overwhelm a target recursive server from serving the legitimate clients.


title: Constraining the dark energy equation of state using Bayes theorem and the Kullback-Leibler divergence

abstract: Data-driven model-independent reconstructions of the dark energy equation of state $w(z)$ are presented using Planck 2015 era CMB, BAO, SNIa and Lyman-$\alpha$ data. These reconstructions identify the $w(z)$ behaviour supported by the data and show a bifurcation of the equation of state posterior in the range $1.5{<}z{<}3$. Although the concordance $\Lambda$CDM model is consistent with the data at all redshifts in one of the bifurcated spaces, in the other a supernegative equation of state (also known as `phantom dark energy') is identified within the $1.5 \sigma$ confidence intervals of the posterior distribution. To identify the power of different datasets in constraining the dark energy equation of state, we use a novel formulation of the Kullback--Leibler divergence. This formalism quantifies the information the data add when moving from priors to posteriors for each possible dataset combination. The SNIa and BAO datasets are shown to provide much more constraining power in comparison to the Lyman-$\alpha$ datasets. Further, SNIa and BAO constrain most strongly around redshift range $0.1-0.5$, whilst the Lyman-$\alpha$ data constrains weakly over a broader range. We do not attribute the supernegative favouring to any particular dataset, and note that the $\Lambda$CDM model was favoured at more than $2$ log-units in Bayes factors over all the models tested despite the weakly preferred $w(z)$ structure in the data.


title: A cluster realization of $U_q(\mathfrak{sl_n})$ from quantum character varieties

abstract: We construct an injective algebra homomorphism of the quantum group $U_q(\mathfrak{sl}_{n+1})$ into a quantum cluster algebra $\mathbf{L}_n$ associated to the moduli space of framed $PGL_{n+1}$-local systems on a marked punctured disk. We obtain a description of the coproduct of $U_q(\mathfrak{sl}_{n+1})$ in terms of the corresponding quantum cluster algebra associated to the marked twice punctured disk, and express the action of the $R$-matrix in terms of a mapping class group element corresponding to the half-Dehn twist rotating one puncture about the other. As a consequence, we realize the algebra automorphism of $U_q(\mathfrak{sl}_{n+1})^{\otimes 2}$ given by conjugation by the $R$-matrix as an explicit sequence of cluster mutations, and derive a refined factorization of the $R$-matrix into quantum dilogarithms of cluster monomials.


title: Connections between centrifugal, stratorotational and radiative instabilities in viscous Taylor-Couette flow

abstract: The `Rayleigh line' $\mu=\eta^2$, where $\mu=\Omega_o/\Omega_i$ and $\eta=r_i/r_o$ are respectively the rotation and radius ratios between inner (subscript `i') and outer (subscript `o') cylinders, is regarded as marking the limit of centrifugal instability (CI) in unstratified inviscid Taylor-Couette flow, for both axisymmetric and non-axisymmetric modes. Non-axisymmetric stratorotational instability (SRI) is known to set in for anticyclonic rotation ratios beyond that line, i.e. $\eta^2<\mu<1$ for axially stably-stratified Taylor-Couette flow, but the competition between CI and SRI in the range $\mu<\eta^2$ has not yet been addressed. In this paper, we establish continuous connections between the two instabilities at finite Reynolds number Re, as previously suggested by M. Le Bars & P. Le Gal, Phys. Rev. Lett. 99, 064502 (2007), making them indistinguishable at onset. Both instabilities are also continuously connected to the radiative instability at finite Re. These results demonstrate the complex impact viscosity has on the linear stability properties of this flow. Several other qualitative differences with inviscid theory were found, among which the instability of a non-axisymmetric mode localized at the outer cylinder without stratification, and the instability of a mode propagating against the inner cylinder rotation with stratification. The combination of viscosity and stratification can also lead to a `collision' between (axisymmetric) Taylor vortex branches, causing the axisymmetric oscillatory state already observed in past experiments. Perhaps more surprising is the instability of a centrifugal-like helical mode beyond the Rayleigh line, caused by the joint effects of stratification and viscosity. The threshold $\mu=\eta^2$ seems to remain, however, an impassable instability limit for axisymmetric modes, regardless of stratification, viscosity, and even disturbance amplitude.


title: Noise Models in Feature-based Stereo Visual Odometry

abstract: Feature-based visual structure and motion reconstruction pipelines, common in visual odometry and large-scale reconstruction from photos, use the location of corresponding features in different images to determine the 3D structure of the scene, as well as the camera parameters associated with each image. The noise model, which defines the likelihood of the location of each feature in each image, is a key factor in the accuracy of such pipelines, alongside optimization strategy. Many different noise models have been proposed in the literature; in this paper we investigate the performance of several. We evaluate these models specifically w.r.t. stereo visual odometry, as this task is both simple (camera intrinsics are constant and known; geometry can be initialized reliably) and has datasets with ground truth readily available (KITTI Odometry and New Tsukuba Stereo Dataset). Our evaluation shows that noise models which are more adaptable to the varying nature of noise generally perform better.


title: A new analytical approach to consistency and overfitting in regularized empirical risk minimization

abstract: This work considers the problem of binary classification: given training data $x_1, \dots, x_n$ from a certain population, together with associated labels $y_1,\dots, y_n \in \left\{0,1 \right\}$, determine the best label for an element $x$ not among the training data. More specifically, this work considers a variant of the regularized empirical risk functional which is defined intrinsically to the observed data and does not depend on the underlying population. Tools from modern analysis are used to obtain a concise proof of asymptotic consistency as regularization parameters are taken to zero at rates related to the size of the sample. These analytical tools give a new framework for understanding overfitting and underfitting, and rigorously connect the notion of overfitting with a loss of compactness.


title: The Sivers asymmetry in $J/\Psi$ and lepton pair production at COMPASS

abstract: The abundant production of lepton pairs via $J/\Psi$ creation at COMPASS, $\pi ^\pm \, p^\uparrow \to J/\Psi \, X \to \ell^+ \ell^- X$, allows a measurement of the transverse Single Spin Asymmetry generated by the Sivers effect. The crucial issue of the sign change of the Sivers function in lepton pair production, with respect to Semi Inclusive Deep Inelastic Scattering processes, can be solved. Predictions for the expected magnitude of the Single Spin Asymmetry, which turns out to be large, are given.


title: Analysis of Chromosome 20 - A Study

abstract: Since the arrival of next-generation sequencing technologies the amount of genetic sequencing data has increased dramatically. This has has fueled an increase in human genetics research. At the same time, with the recent advent of technologies in processing large data sets, lot of these technologies are proving valuable and efficient in analyzing these huge datasets. In this paper we use some of these technologies to analyze genetic sequencing data of 1000 Genomes Project,produce and evaluate a framework to process the sequencing data thereof and look into structural variations with respect to population groups.


title: Wavefront Sensing in Space from the PICTURE-B Sounding Rocket

abstract: A NASA sounding rocket for high contrast imaging with a visible nulling coronagraph, the Planet Imaging Coronagraphic Technology Using a Reconfigurable Experimental Base (PICTURE-B) payload has made two suborbital attempts to observe the warm dust disk inferred around Epsilon Eridani. We present results from the November 2015 launch demonstrating active wavefront sensing in space with a piezoelectric mirror stage and a micromachine deformable mirror along with precision pointing and lightweight optics in space.


title: Obstructing Visibilities with One Obstacle

abstract: Obstacle representations of graphs have been investigated quite intensely over the last few years. We focus on graphs that can be represented by a single obstacle. Given a (topologically open) polygon $C$ and a finite set $P$ of points in general position in the complement of $C$, the visibility graph $G_C(P)$ has a vertex for each point in $P$ and an edge $pq$ for any two points $p$ and $q$ in $P$ that can see each other, that is, $\overline{pq} \cap C=\emptyset$. We draw $G_C(P)$ straight-line. Given a graph $G$, we want to compute an obstacle representation of $G$, that is, an obstacle $C$ and a set of points $P$ such that $G=G_C(P)$. The complexity of this problem is open, even for the case that the points are exactly the vertices of a simple polygon and the obstacle is the complement of the polygon-the simple-polygon visibility graph problem. There are two types of obstacles; an inside obstacle lies in a bounded component of the complement of the visibility drawing, whereas an outside obstacle lies in the unbounded component. We show that the class of graphs with an inside-obstacle representation is incomparable with the class of graphs that have an outside-obstacle representation. We further show that any graph with at most seven vertices or circumference at most 6 has an outside-obstacle representation, which does not hold for a specific graph with 8 vertices and circumference 8. Finally, we consider the outside-obstacle graph sandwich problem: given graphs $G$ and $H$ on the same vertex set, is there a graph $K$ such that $G \subseteq K \subseteq H$ and $K$ has an outside-obstacle representation? We show that this problem is NP-hard even for co-bipartite graphs. With slight modifications, our proof also shows that the inside-obstacle graph sandwich problem, the single-obstacle graph sandwich problem, and the simple-polygon visibility graph sandwich problem are all NP-hard.


title: Meaningful Models: Utilizing Conceptual Structure to Improve Machine Learning Interpretability

abstract: The last decade has seen huge progress in the development of advanced machine learning models; however, those models are powerless unless human users can interpret them. Here we show how the mind's construction of concepts and meaning can be used to create more interpretable machine learning models. By proposing a novel method of classifying concepts, in terms of 'form' and 'function', we elucidate the nature of meaning and offer proposals to improve model understandability. As machine learning begins to permeate daily life, interpretable models may serve as a bridge between domain-expert authors and non-expert users.


title: Solution of deformed Einstein equations and quantum black holes

abstract: Recently one and two-parameter deformed Einstein equations have been studied for extremal quantum black holes which have been proposed to obey deformed statistics by Strominger. In this study, we give a deeper insight to the deformed Einstein equations and consider the solutions of these equations for the extremal quantum black holes. We then represent the implications of the solutions, such that the deformation parameters lead the charged black holes to have a smaller mass than the usual Reissner-Nordstr\"om black holes. This reduction in mass of a usual black hole can be considered as a transition from classical to quantum black hole regime.


title: The System of Molecular Clouds in the Gould Belt

abstract: Based on high-latitude molecular clouds with highly accurate distance estimates taken from the literature, we have redetermined the parameters of their spatial orientation. This system can be approximated by a 350x235x140 pc ellipsoid inclined by the angle i=17+/-2 degrees to the Galactic plane with the longitude of the ascending node l{\Omega}=337+/-1 degrees. Based on the radial velocities of the clouds, we have found their group velocity relative to the Sun to be (u0,v0,w0) = (10.6,18.2,6.8)+/-(0.9,1.7,1.5) km/s. The trajectory of the center of the molecular cloud system in the past in a time interval of ~60 Myr has been constructed. Using data on masers associated with low-mass protostars, we have calculated the space velocities of the molecular complexes in Orion, Taurus, Perseus, and Ophiuchus. Their motion in the past is shown to be not random.


title: Clues on the Majorana scale from scalar resonances at the LHC

abstract: In order to address the observation of the neutrino oscillations and the metastability of the Standard Model, we extend the fermion sector with two right-handed (i.e. sterile) neutrinos, and the scalar sector of the SM with a real scalar, the Hill field. The latter takes the role of a Majoron and generates the Majorana masses for the neutrino sector, such that the particle spectrum features two CP-even scalars $h_1$ and $h_2$, and also two heavy, mass degenerate neutrinos. When the $h_1$ is identified with the scalar resonance at ~125 GeV and the condition is imposed that the $h_1$ self coupling and its running vanish at the Planck scale, the scalar mixing and the vacuum expectation value of the Hill field are fixed by the $h_2$ mass. The $h_2$ can be searched for at the LHC, and it has prospects of being discovered for the target integrated luminosities of the HL-LHC and the Future Circular hadron Collider (FCC-hh) when its mass is on the weak scale. The knowledge of the $h_2$ mass and its decay properties can yield an insight into its coupling to the heavy neutrinos, and thus also on the heavy neutrino mass scale. This yields an interesting connection between potentially detectable heavy scalars in high-energy proton collisions and the mass scale of the heavy neutrinos that is testable at the LHC and at future colliders.


title: Excited-state quantum phase transition in the Rabi model

abstract: The Rabi model, a two-level atom coupled to a harmonic oscillator, can undergo a second-order quantum phase transition (QPT) [M. -J. Hwang et al, Phys. Rev. Lett. 115, 180404 (2015)]. Here we show that the Rabi QPT accompanies critical behavior in the higher energy excited states, i.e., the excited-state QPT (ESQPT). We derive analytic expressions for the semiclassical density of states, which shows a logarithmic divergence at a critical energy eigenvalue in the broken symmetry (superradiant) phase. Moreover, we find that the logarithmic singularities in the density of states leads to singularities in the relevant observables in the system such as photon number and atomic polarization. We corroborate our analytical semiclassical prediction of the ESQPT in the Rabi model with its numerically exact quantum mechanical solution.


title: Frustrated double and single ionization in a two-electron triatomic molecule H$^+_3$

abstract: Using a semi-classical model, we study the formation of highly excited neutral fragments during the fragmentation of $\mathrm{H_3^+}$, a two-electron triatomic molecule, driven by an intense near-IR laser field. To do so, we first formulate a microcanonical distribution for arbitrary one-electron triatomic molecules. We then study frustrated double and single ionization in strongly-driven $\mathrm{H_3^+}$ and compute the kinetic energy release of the nuclei for these two processes. Moreover, we investigate the dependence of frustrated ionization on the strength of the laser field as well as on the geometry of the initial molecular state.


title: Liquid-grain mixing suppresses droplet spreading and splashing during impact

abstract: Would a raindrop impacting on a coarse beach behave differently from that impacting on a desert of fine sand? We study this question by a series of model experiments, where the packing density of the granular target, the wettability of individual grains, the grain size, the impacting liquid, and the impact speed are varied. We find that by increasing the grain size and/or the wettability of individual grains the maximum droplet spreading undergoes a transition from a capillary regime towards a viscous regime, and splashing is suppressed. The liquid-grain mixing is discovered to be the underlying mechanism. An effective viscosity is defined accordingly to quantitatively explain the observations.


title: Quantile Graphical Models: Prediction and Conditional Independence with Applications to Systemic Risk

abstract: We propose two types of Quantile Graphical Models (QGMs) --- Conditional Independence Quantile Graphical Models (CIQGMs) and Prediction Quantile Graphical Models (PQGMs). CIQGMs characterize the conditional independence of distributions by evaluating the distributional dependence structure at each quantile index. As such, CIQGMs can be used for validation of the graph structure in the causal graphical models (\cite{pearl2009causality, robins1986new, heckman2015causal}). One main advantage of these models is that we can apply them to large collections of variables driven by non-Gaussian and non-separable shocks. PQGMs characterize the statistical dependencies through the graphs of the best linear predictors under asymmetric loss functions. PQGMs make weaker assumptions than CIQGMs as they allow for misspecification. Because of QGMs' ability to handle large collections of variables and focus on specific parts of the distributions, we could apply them to quantify tail interdependence. The resulting tail risk network can be used for measuring systemic risk contributions that help make inroads in understanding international financial contagion and dependence structures of returns under downside market movements. We develop estimation and inference methods for QGMs focusing on the high-dimensional case, where the number of variables in the graph is large compared to the number of observations. For CIQGMs, these methods and results include valid simultaneous choices of penalty functions, uniform rates of convergence, and confidence regions that are simultaneously valid. We also derive analogous results for PQGMs, which include new results for penalized quantile regressions in high-dimensional settings to handle misspecification, many controls, and a continuum of additional conditioning events.


title: Pauli Blockade in a Few-Hole PMOS Double Quantum Dot limited by Spin-Orbit Interaction

abstract: We report on hole compact double quantum dots fabricated using conventional CMOS technology. We provide evidence of Pauli spin blockade in the few hole regime which is relevant to spin qubit implementations. A current dip is observed around zero magnetic field, in agreement with the expected behavior for the case of strong spin-orbit. We deduce an intradot spin relaxation rate $\approx$120\,kHz for the first holes, an important step towards a robust hole spin-orbit qubit.


title: Evolution of the dust-to-metals ratio in high-redshift galaxies probed by GRB-DLAs

abstract: Context. Several issues regarding the nature of dust at high redshift remain unresolved: its composition, its production and growth mechanisms, and its effect on background sources. Aims. We provide a more accurate relation between dust depletion levels and dust-to-metals ratio (DTM), and to use the DTM to investigate the origin and evolution of dust in the high-redshift Universe via Gamma-ray burst damped Lyman-alpha absorbers (GRB-DLAs). Methods. We use absorption-line measured metal column densities for a total of 19 GRB-DLAs, including five new GRB afterglow spectra from VLT/X-shooter. We use the latest linear models to calculate the dust depletion strength factor in each DLA. Using these values we calculate total dust and metal column densities to determine a DTM. We explore the evolution of DTM with metallicity, and compare it to previous trends in DTM measured with different methods. Results. We find significant dust depletion in 16 of our 19 GRB-DLAs, yet 18 of the 19 have a DTM significantly lower than the Milky Way. We find that DTM is positively correlated with metallicity, which supports a dominant ISM grain-growth mode of dust formation. We find a substantial discrepancy between the dust content measured from depletion and that derived from the total V-band extinction, $A_V$ , measured by fitting the afterglow SED. We advise against using a measurement from one method to estimate that from the other until the discrepancy can be resolved.


title: Degenerating Hermitian metrics and spectral geometry of the canonical bundle

abstract: Let $(X,h)$ be a compact and irreducible Hermitian complex space of complex dimension $m$. In this paper we are interested in the Dolbeault operator acting on the space of $L^2$ sections of the canonical bundle of $reg(X)$, the regular part of $X$. More precisely let $\overline{\mathfrak{d}}_{m,0}:L^2\Omega^{m,0}(reg(X),h)\rightarrow L^2\Omega^{m,1}(reg(X),h)$ be an arbitrarily fixed closed extension of $\overline{\partial}_{m,0}:L^2\Omega^{m,0}(reg(X),h)\rightarrow L^2\Omega^{m,1}(reg(X),h)$ where the domain of the latter operator is $\Omega_c^{m,0}(reg(X))$. We establish various properties such as closed range of $\overline{\mathfrak{d}}_{m,0}$, compactness of the inclusion $\mathcal{D}(\overline{\mathfrak{d}}_{m,0})\hookrightarrow L^2\Omega^{m,0}(reg(X),h)$ where $\mathcal{D}(\overline{\mathfrak{d}}_{m,0})$, the domain of $\overline{\mathfrak{d}}_{m,0}$, is endowed with the corresponding graph norm, and discreteness of the spectrum of the associated Hodge-Kodaira Laplacian $\overline{\mathfrak{d}}_{m,0}^*\circ \overline{\mathfrak{d}}_{m,0}$ with an estimate for the growth of its eigenvalues. Several corollaries such as trace class property for the heat operator associated to $\overline{\mathfrak{d}}_{m,0}^*\circ \overline{\mathfrak{d}}_{m,0}$, with an estimate for its trace, are derived. Finally in the last part we provide several applications to the Hodge-Kodaira Laplacian in the setting of both compact irreducible Hermitian complex spaces with isolated singularities and complex projective surfaces.


title: Implications of gamma-ray observations on proton models of UHECR

abstract: The origin of ultra high energy cosmic rays (UHECR) is still unknown. However, great progress has been achieved in past years due to the good quality and large statistics in experimental data collected by the current observatories. The data of the Pierre Auger Observatory show that the composition of the UHECRs becomes progressively lighter starting from $10^{17}$ eV up to $\sim 10^{18.3}$ eV and then, beyond that energy, it becomes increasingly heavier. These analyses are subject to important systematic uncertainties due to the use of hadronic interaction models that extrapolate lower energy accelerator data to the highest energies. Although proton models of UHECRs are disfavored by these results, they cannot be completely ruled out. It is well known that the energy spectra of gamma rays and neutrinos, produced during propagation of these very energetic particles through the intergalactic medium, are a useful tool to constrain the spectrum models. In particular, it has recently been shown that the neutrino upper limits obtained by IceCube challenge the proton models at 95% CL. In this work we study the constraints imposed by the extragalactic gamma-ray background, measured by Fermi-LAT, on proton models of UHECRs. In particular, we make use of the extragalactic gamma-ray background flux, integrated from 50 GeV to 2 TeV, that originates in point sources, which has recently been obtained by the Fermi-LAT collaboration, in combination with the neutrino upper limits, to constrain the emission of UHECRs at high redshits ($z>1$), in the context of the proton models.


title: High-Performance Tensor Contraction without Transposition

abstract: Tensor computations--in particular tensor contraction (TC)--are important kernels in many scientific computing applications. Due to the fundamental similarity of TC to matrix multiplication (MM) and to the availability of optimized implementations such as the BLAS, tensor operations have traditionally been implemented in terms of BLAS operations, incurring both a performance and a storage overhead. Instead, we implement TC using the flexible BLIS framework, which allows for transposition (reshaping) of the tensor to be fused with internal partitioning and packing operations, requiring no explicit transposition operations or additional workspace. This implementation, TBLIS, achieves performance approaching that of MM, and in some cases considerably higher than that of traditional TC. Our implementation supports multithreading using an approach identical to that used for MM in BLIS, with similar performance characteristics. The complexity of managing tensor-to-matrix transformations is also handled automatically in our approach, greatly simplifying its use in scientific applications.


title: Bayes factor consistency

abstract: Good large sample performance is typically a minimum requirement of any model selection criterion. This article focuses on the consistency property of the Bayes factor, a commonly used model comparison tool, which has experienced a recent surge of attention in the literature. We thoroughly review existing results. As there exists such a wide variety of settings to be considered, e.g. parametric vs. nonparametric, nested vs. non-nested, etc., we adopt the view that a unified framework has didactic value. Using the basic marginal likelihood identity of Chib (1995), we study Bayes factor asymptotics by decomposing the natural logarithm of the ratio of marginal likelihoods into three components. These are, respectively, log ratios of likelihoods, prior densities, and posterior densities. This yields an interpretation of the log ratio of posteriors as a penalty term, and emphasizes that to understand Bayes factor consistency, the prior support conditions driving posterior consistency in each respective model under comparison should be contrasted in terms of the rates of posterior contraction they imply.


title: Neutrino oscillations at DUNE with improved energy reconstruction

abstract: We study the physics reach of the long-baseline oscillation analysis of the DUNE experiment when realistic simulations are used to estimate its neutrino energy reconstruction capabilities. Our studies indicate that significant improvements in energy resolution compared to what is customarily assumed are plausible. This improved energy resolution can increase the sensitivity to leptonic CP violation in two ways. On the one hand, the CP-violating term in the oscillation probability has a characteristic energy dependence that can be better reproduced. On the other hand, the second oscillation maximum, especially sensitive to $\delta_{CP}$, is better reconstructed. These effects lead to a significant improvement in the fraction of values of $\delta_{CP}$ for which a $5 \sigma$ discovery of leptonic CP-violation would be possible. The precision of the $\delta_{CP}$ measurement could also be greatly enhanced, with a reduction of the maximum uncertainties from $26^\circ$ to $18^\circ$ for a 300~MW$\cdot$kt$\cdot$yr exposure. We therefore believe that this potential gain in physics reach merits further investigations of the detector performance achievable in DUNE.


title: Testing Lorentz symmetry with Lunar Laser Ranging

abstract: Lorentz symmetry violations can be parametrized by an effective field theory framework that contains both general relativity and the standard model of particle physics called the standard-model extension (SME). We present new constraints on pure gravity SME coefficients obtained by analyzing lunar laser ranging (LLR) observations. We use a new numerical lunar ephemeris computed in the SME framework and we perform a LLR data analysis using a set of 20721 normal points covering the period of August, 1969 to December, 2013. We emphasize that linear combination of SME coefficients to which LLR data are sensitive and not the same as those fitted in previous postfit residuals analysis using LLR observations and based on theoretical grounds. We found no evidence for Lorentz violation at the level of $10^{-8}$ for $\bar{s}^{TX}$, $10^{-12}$ for $\bar{s}^{XY}$ and $\bar{s}^{XZ}$, $10^{-11}$ for $\bar{s}^{XX}-\bar{s}^{YY}$ and $\bar{s}^{XX}+\bar{s}^{YY}-2\bar{s}^{ZZ}-4.5\bar{s}^{YZ}$ and $10^{-9}$ for $\bar{s}^{TY}+0.43\bar{s}^{TZ}$. We improve previous constraints on SME coefficient by a factor up to 5 and 800 compared to postfit residuals analysis of respectively binary pulsars and LLR observations.


title: Studies on charge production from Cs2Te photocathodes in the PITZ L-band normal conducting radio frequency photo injector

abstract: This paper discusses the behavior of electron bunch charge produced in an L-band normal conducting radio frequency cavity (RF gun) from Cs2Te photocathodes illuminated with ps-long UV laser pulses when the laser transverse distribution consists of a flat-top core with Gaussian-like decaying halo. The produced charge shows a linear dependence at low laser pulse energies as expected in the quantum efficiency limited emission regime, while its dependence on laser pulse energy is observed to be much weaker for higher values, due to space charge limited emission. However, direct plug-in of experimental parameters into the space charge tracking code ASTRA yields lower output charge in the space charge limited regime compared to measured values. The rate of increase of the produced charge at high laser pulse energies close to the space charge limited emission regime seems to be proportional to the amount of halo present in the radial laser profile since the charge from the core has saturated already. By utilizing core + halo particle distributions based on measured radial laser profiles, ASTRA simulations and semi-analytical emission models reproduce the behavior of the measured charge for a wide range of RF gun and laser operational parameters within the measurement uncertainties.


title: Charge-exchange, ionization and excitation in low-energy Li$^{+}-$ Ar, K$^{+}-$ Ar, and Na$^{+}-$He collisions

abstract: Absolute cross sections are measured for charge-exchange, ionization, and excitation within the same experimental setup for the Li$^{+}-$Ar, K$^{+}-$ Ar, and Na$^{+}-$ He collisions in the ion energy range $0.5-10$ keV. Results of our measurements along with existing experimental data and the schematic correlation diagrams are used to analyze and determine the mechanisms for these processes. The experimental results show that the charge-exchange processes are realized with high probabilities and electrons are predominately captured in ground states. The cross section ratio for charge exchange, ionization and excitation processes roughly attains the value $10:2:1$, respectively. The contributions of various partial inelastic channels to the total ionization cross sections are estimated and a primary mechanism for the process is defined. The energy-loss spectrum, in addition, is applied to estimate the relative contribution of different inelastic channels and to determine the mechanisms for the ionization and for some excitation processes of Ar resonance lines for the K$^{+}-$Ar collision system. The excitation function for the helium, as well as for the sodium doublet lines for the Na$^{+}-$He collision system, reveals some unexpected features and a mechanism to explain this observation is suggested.


title: Characterizing Corson and Valdivia compact spaces

abstract: We give a new characterization of Valdivia compact spaces: A compact space is Valdivia if and only if it has a dense commutatively monotonically retractable subspace. This result solves Problem 5.12 from \cite{sal-rey}. Besides, we introduce the notion of full $c$-skeleton and prove that a compact space is Corson if and only if it has a full $c$-skeleton.


title: CLEX: Yet Another Supercomputer Architecture?

abstract: We propose the CLEX supercomputer topology and routing scheme. We prove that CLEX can utilize a constant fraction of the total bandwidth for point-to-point communication, at delays proportional to the sum of the number of intermediate hops and the maximum physical distance between any two nodes. Moreover, % applying an asymmetric bandwidth assignment to the links, all-to-all communication can be realized $(1+o(1))$-optimally both with regard to bandwidth and delays. This is achieved at node degrees of $n^{\varepsilon}$, for an arbitrary small constant $\varepsilon\in (0,1]$. In contrast, these results are impossible in any network featuring constant or polylogarithmic node degrees. Through simulation, we assess the benefits of an implementation of the proposed communication strategy. Our results indicate that, for a million processors, CLEX can increase bandwidth utilization and reduce average routing path length by at least factors $10$ respectively $5$ in comparison to a torus network. Furthermore, the CLEX communication scheme features several other properties, such as deadlock-freedom, inherent fault-tolerance, and canonical partition into smaller subsystems.


title: Review of lattice results concerning low-energy particle physics

abstract: We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.


title: Lie bialgebra structures on 2-step nilpotent graph algebras

abstract: We generalize a result on the Heisenberg Lie algebra that gives restrictions to possible Lie bialgebra cobrackets on 2-step nilpotent algebras with some additional properties. For the class of 2-step nilpotent Lie algebras coming from graphs, we describe these extra properties in a very easy graph-combinatorial way. We exhibit applications for $\mathfrak f_n$, the free 2-step nilpotent Lie algebra.


title: A time-stepping DPG scheme for the heat equation

abstract: We introduce and analyze a discontinuous Petrov-Galerkin method with optimal test functions for the heat equation. The scheme is based on the backward Euler time stepping and uses an ultra-weak variational formulation at each time step. We prove the stability of the method for the field variables (the original unknown and its gradient weighted by the square root of the time step) and derive a C\'ea-type error estimate. For low-order approximation spaces this implies certain convergence orders when time steps are not too small in comparison with mesh sizes. Some numerical experiments are reported to support our theoretical results.


title: Observation of the quantum paradox of separation of a single photon from one of its properties

abstract: We report an experimental realization of the quantum paradox of the separation of a single photon from one of its properties (the so-called "quantum Cheshire cat"). We use a modified Sagnac interferometer with displaced paths to produce appropriately pre- and postselected states of heralded single photons. Weak measurements of photon presence and circular polarization are performed in each arm of the interferometer by introducing weak absorbers and small polarization rotations and analyzing changes in the postselected signal. The absorber is found to have an appreciable effect only in one arm of the interferometer, while the polarization rotation significantly affects the signal only when performed in the other arm. We carry out both sequential and simultaneous weak measurements and find good agreement between measured and predicted weak values. In the language of Aharonov et al. and in the sense of the ensemble averages described by weak values, the experiment establishes the separation of a particle from one its properties during the passage through the interferometer.


title: Performance analysis for the new G-2 experiment

abstract: The new g-2 experiment at Fermilab aims to measure the muon anomalous magnetic moment by a fourfold improvement in precision compared to the BNL experiment. Achieving this goal requires the delivery of highly polarized 3.094 GeV/c muons with a narrow +-0.5% {\Delta}p/p acceptance to the storage ring. In this study, we describe a muon capture and transport scheme that should meet this requirement. First, we present the conceptual design of our proposed scheme wherein we describe its basic features. Then, we detail our numerical model and present a complete end-to-end simulation of all g-2 beamlines.


title: Public Naming of Planets and Planetary Satellites

abstract: While one of the IAU's missions is to "serve as the internationally recognized authority for assigning designations to celestial bodies and surface features on them", the participation of the public in the naming of celestial objects has been a little-known, but decade-long tradition of the IAU. While reiterating its opposition to having the public pay to give a name to an exoplanet, the IAU Executive Committee nonetheless recognized the right of organizations to invite public, international exoplanet naming or voting campaigns. To this end, clear selection rules were to be defined by the IAU, inviting mutual collaboration, the goal being to sanction the campaign and officially approve the resulting names, for the sake of boosting the public's interest in astronomy and at the same time reaffirm the authority of the IAU. In no way were these names supposed to supersede the designations in use by professional astronomers. Since the field of exoplanet research was still "virgin" for naming but is now becoming scientifically mature (with hundreds of well-characterized, confirmed objects) and so appealing to the public's interest and imagination, IAU decided to organize a naming contest as the leader of a new project with a potentially high worldwide impact, the first since the International Year of Astronomy in 2009, called "NameExoWorlds".


title: The Representation Dimension of a Selfinjective Algebra of Wild Tilted Type

abstract: We prove that the representation dimension of a selfinjective algebra of wild tilted type is equal to three, and give an explicit construction of an Auslander generator of its module category. We also show that if a connected selfinjective algebra admits an acyclic generalised standard Auslander-Reiten component then its representation dimension is equal to three.


title: Multi-wavelength Study of Transition Region Penumbral Subarcsecond Bright Dots Using IRIS and NST

abstract: Using high-resolution transition region (TR) observations taken by the Interface Region Imaging Spectrograph (IRIS) mission, Tian et al. (2014b) revealed numerous short-lived subarcsecond bright dots (BDs) above sunspots (mostly located in the penumbrae), which indicate yet unexplained small-scale energy releases. Moreover, whether these subarcsecond TR brightenings have any signature in the lower atmosphere and how they are formed are still not fully resolved. This paper presents a multi-wavelength study of the TR penumbral BDs using a coordinated observation of a near disk-center sunspot with IRIS and the 1.6 m New Solar Telescope (NST) at the Big Bear Solar Observatory. NST provides high-resolution chromospheric and photospheric observations with narrow-band H-alpha imaging spectroscopy and broad-band TiO images, respectively, complementary to IRIS TR observations. A total of 2692 TR penumbral BDs are identified from a 37-minute time series of IRIS 1400 A slitjaw images. Their locations tend to be associated more with downflowing and darker fibrils in the chromosphere, and weakly associated with bright penumbral features in the photosphere. However, temporal evolution analyses of the BDs show that there is no consistent and convincing brightening response in the chromosphere. These results are compatible with a formation mechanism of the TR penumbral BDs by falling plasma from coronal heights along more vertical and dense magnetic loops. The BDs may also be produced by small-scale impulsive magnetic reconnection taking place sufficiently high in the atmosphere that has no energy release in the chromosphere.


title: Optimal Tree Hash Modes: the Case of Trees Having their Leaves at All the Levels

abstract: A recent work shows how we can optimize a tree based mode of operation for a hash function where the sizes of input message blocks and digest are the same, subject to the constraint that the involved tree structure has all its leaves at the same depth. In this work, we show that we can further optimize the running time of such a mode by using a tree having leaves at all its levels. We make the assumption that the input message block has a size a multiple of that of the digest and denote by $d$ the ratio block size over digest size. The running time is evaluated in terms of number of operations performed by the hash function, i.e. the number of calls to its underlying function. It turns out that a digest can be computed in $\lceil \log_{d+1} (l/2) \rceil+2$ evaluations of the underlying function using $\lceil l/2 \rceil$ processors, where $l$ is the number of blocks of the message. Other results of interest are discussed, such as the optimization of the parallel running time for a tree of restricted height.


title: Mechanical stability study for Integrable Optics Test Accelerator at Fermilab

abstract: The Integrable Optics Test Accelerator (IOTA) is proposed for operation at Fermilab. The goal of IOTA is to create practical nonlinear accelerator focusing systems with a large frequency spread and stable particle motion. The IOTA is a 40 m circumference, 150 MeV (e-), 2.5 MeV (p+) diagnostic test ring. A heavy low frequency steel floor girder is proposed as the primary tier for IOTA device component support. Two design lengths; (8) 3.96 m and (2) 3.1 m long girders with identical cross section completely encompass the ring. This study focuses on the 3.96 m length girder and the development of a working prototype. Hydrostatic Level Sensor (HLS), temperature, metrology and fast motion measurements characterize the anticipated mechanical stability of the IOTA ring.


title: Full Monte-Carlo description of the Moscow State University Extensive Air Shower experiment

abstract: The Moscow State University Extensive Air Shower (EAS-MSU) array studied high-energy cosmic rays with primary energies ~(1-500) PeV in the Northern hemisphere. The EAS-MSU data are being revisited following recently found indications to an excess of muonless showers, which may be interpreted as the first observation of cosmic gamma rays at ~100 PeV. In this paper, we present a complete Monte-Carlo model of the surface detector which results in a good agreement between data and simulations. The model allows us to study the performance of the detector and will be used to obtain physical results in further studies.


title: A WENO-based Method of Line Transpose Approach for Vlasov Simulations

abstract: In this paper, a high order implicit Method of Line Transpose (MOL$^T$ ) method based on a weighted essentially non-oscillatory (WENO) methodology is developed for one-dimensional linear transport equations and further applied to the Vlasov-Poisson (VP) simulations via dimensional splitting. In the MOL$^T$ framework, the time variable is first discretized by a diagonally implicit strong-stability-preserving Runge-Kutta method, resulting in a boundary value problem (BVP) at the discrete time levels. Then an integral formulation coupled with a high order WENO methodology is employed to solve the BVP. As a result, the proposed scheme is high order accurate in both space and time and free of oscillations even though the solution is discontinuous or has sharp gradients. Moreover, the scheme is able to take larger time step evolution compared with an explicit MOL WENO scheme with the same order of accuracy. The desired positivity-preserving (PP) property of the scheme is further attained by incorporating a newly proposed high order PP limiter. We perform numerical experiments on several benchmarks including linear advection, solid body rotation problem; and on the Landau damping, two-stream instabilities, bump-on-tail, and plasma sheath by solving the VP system. The efficacy and efficiency of the proposed scheme is numerically verified.


title: Local integrands for two-loop QCD amplitudes

abstract: In this talk we review the recent computation of the five- and six-gluon two-loop amplitudes in Yang-Mills theory using local integrands which make the infrared pole structure manifest. We make some remarks on the connection with BCJ relations and the all-multiplicity structure.


title: Hamburger Sternwarte plate archives: Historic long-term variability study of active galaxies based on digitized photographic plates

abstract: At the Hamburger Sternwarte an effort was started in 2010 with the aim of digitizing its more than 45000 photographic plates and films stored in its plate archives. At the time of writing, more than 31000 plates have already been made available on the Internet for researchers, historians, as well as for the interested public. The digitization process and the Internet presentation of the plates and accompanying hand written material (plate envelopes, logbooks, observer notes) are presented here. To fully exploit the unique photometric and astrometric data, stored on the plates, further processing steps are required including registering the plate to celestial coordinates, masking of the plates, and a calibration of the photo-emulsion darkening curve. To demonstrate the correct functioning of these procedures, historical light curves of two bright BL Lac type active galactic nuclei are extracted. The resulting light curve of the blazar 1ES 1215+303 exhibits a large decrease in the magnitude from $14.25^{+0.07}_{-0.12}$ to $15.94^{+0.09}_{-0.13}$ in about 300 days, which proves the variability in the optical region. Furthermore, we compare the measured magnitudes for the quasar 3C~273 with contemporaneous measurements.


title: Notes on S-folds and $\mathcal{N}=3$ theories

abstract: We consider D3 branes in presence of an S-fold plane. The latter is a non perturbative object, arising from the combined projection of an S-duality twist and a discrete orbifold of the R-symmetry group. This construction naively gives rise to 4d $\mathcal{N}=3$ SCFTs. Nevertheless it has been observed that in some cases supersymmetry is enhanced to $\mathcal{N}=4$. In this paper we study the explicit counting of degrees of freedom arising from vector multiplets associated to strings suspended between the D3 branes probing the S-fold. We propose that, for trivial discrete torsion, there is no vector multiplet associated to $(1,0)$ strings stretched between a brane and its image. We then focus on the case of rank 2 $\mathcal{N}=3$ theory that enhances to $SU(3)$ $\mathcal{N}=4$ SYM, explicitly spelling out the isomorphism between the BPS-spectrum of the manifestly $\mathcal{N}=3$ theory and that of three D3 branes in flat spacetime. Subsequently, we consider 3-pronged strings in these setups and show how wall-crossing in the S-fold background implies wall crossing in the flat geometry. This can be considered a consistency check of the \emph{conjectured} SUSY enhancement. We also find that the above isomorphism implies that a $(1,0)$ string, suspended between a brane and its image in the S-fold, corresponds to a 3-string junction in the flat geometry. This is in agreement with our claim on the absence of a vector multiplet associated to such $(1,0)$ strings. This is because the 3-string junction in flat geometry gives rise to a $1/4$-th BPS multiplet of the $\mathcal{N}=4$ algebra. Such multiplets always include particles with spin $>1$ as opposed to a vector multiplet which is restricted by the requirement that the spins must be $\leq 1$.


title: HOMFLY-PT homology for general link diagrams and braidlike isotopy

abstract: Khovanov and Rozansky's categorification of the HOMFLY-PT polynomial is invariant under braidlike isotopies for any link diagram and Markov moves for braid closures. To define HOMFLY-PT homology, they required a link to be presented as a braid closure, because they did not prove invariance under the other oriented Reidemeister moves. In this text we prove that the Reidemeister IIb move fails in HOMFLY-PT homology by using virtual crossing filtrations of the author and Rozansky. The decategorification of HOMFLY-PT homology for general link diagrams gives a deformed version of the HOMFLY-PT polynomial, $P^{b}(D)$, which can be used to detect nonbraidlike isotopies. Finally, we will use $P^{b}(D)$ to prove that HOMFLY-PT homology is not an invariant of virtual links, even when virtual links are presented as virtual braid closures.


title: A multilevel framework for sparse optimization with application to inverse covariance estimation and logistic regression

abstract: Solving l1 regularized optimization problems is common in the fields of computational biology, signal processing and machine learning. Such l1 regularization is utilized to find sparse minimizers of convex functions. A well-known example is the LASSO problem, where the l1 norm regularizes a quadratic function. A multilevel framework is presented for solving such l1 regularized sparse optimization problems efficiently. We take advantage of the expected sparseness of the solution, and create a hierarchy of problems of similar type, which is traversed in order to accelerate the optimization process. This framework is applied for solving two problems: (1) the sparse inverse covariance estimation problem, and (2) l1-regularized logistic regression. In the first problem, the inverse of an unknown covariance matrix of a multivariate normal distribution is estimated, under the assumption that it is sparse. To this end, an l1 regularized log-determinant optimization problem needs to be solved. This task is challenging especially for large-scale datasets, due to time and memory limitations. In the second problem, the l1-regularization is added to the logistic regression classification objective to reduce overfitting to the data and obtain a sparse model. Numerical experiments demonstrate the efficiency of the multilevel framework in accelerating existing iterative solvers for both of these problems.


title: Influences of transversely-isotropic rheology and translational diffusion on the stability of active suspensions

abstract: Suspensions of self-motile, elongated particles are a topic of significant current interest, exemplifying a form of `active matter'. Examples include self-propelling bacteria, algae and sperm, and artificial swimmers. Ericksen's model of a transversely-isotropic fluid [J. L. Ericksen, Colloid Polym. Sci. 173(2):117-122 (1960)] treats suspensions of non-motile particles as a continuum with an evolving preferred direction; this model describes fibrous materials as diverse as extracellular matrix, textile tufts and plant cell walls. Director-dependent effects are incorporated through a modified stress tensor with four viscosity-like parameters. By making fundamental connections with recent models for active suspensions, we propose a modification to Ericksen's model, mainly the inclusion of self motility; this can be considered the simplest description of an oriented suspension including transversely-isotropic effects. Motivated by the fact that transversely-isotropic fluids exhibit modified flow stability, we conduct a linear stability analysis of two distinct cases, aligned and isotropic suspensions of elongated active particles. Novel aspects include the anisotropic rheology and translational diffusion. In general anisotropic effects increase the instability of small perturbations, whilst translational diffusion stabilises a range of wave-directions and, in some cases, a finite range of wave-numbers, thus emphasising that both anisotropy and translational diffusion can have important effects in these systems.


title: Phosphorene as an Anode Material for High Performance Lithium-Ion Battery: First Principle Study and Experimental Measurement

abstract: The prospects of phosphorene as an anode material for high performance Li-ion battery was systematically investigated from the first principle calculations and experimental measurements. The diffusion energy barriers of a Li atom moving along various orientations on phosphorene layer were calculated from the Li adsorption energy landscape. It was found that the diffusion mobility of a Li atom along the zigzag direction in the valley of phosphorene could be about 7 to 11 orders of magnitude faster than that along the other directions, indicating its ultrafast and anisotropic diffusivity. The lithium insertion in phosphorene was studied considering various LinP16 configurations (n=1~16). It was found that phosphorene could accommodate up to one Li per P atom (i.e., Li16P16), and the predicted theoretical value of the Li capacity for a single layered phosphorene can reach about 865 mAh/g. Our experimental measurement on the Li capacity for a network of a few layered phosphorene can reach a reversible stable value of ~ 453 mAh/g even after 50 cycles. In particular, it was found that, even at the high Li concentration (e.g., x = 1 in LixP), there was no Li clustering and the structure of phosphorene is reversible during the lithium intercalation. Our results clearly show that phosphorene has promise as a novel anode material for high performance Li-ion batteries.


title: The Evolution of Sex through the Baldwin Effect

abstract: This paper suggests that the fundamental haploid-diploid cycle of eukaryotic sex exploits a rudimentary form of the Baldwin effect. With this explanation for the basic cycle, the other associated phenomena can be explained as evolution tuning the amount and frequency of learning experienced by an organism. Using the well-known NK model of fitness landscapes it is shown that varying landscape ruggedness varies the benefit of the haploid-diploid cycle, whether based upon endomitosis or syngamy. The utility of pre-meiotic doubling and recombination during the cycle are also shown to vary with landscape ruggedness. This view is suggested as underpinning, rather than contradicting, many existing explanations for sex.


title: Quantum smoothing for classical mixtures

abstract: In quantum mechanics, wave functions and density matrices represent our knowledge about a quantum system and give probabilities for the outcomes of measurements. If the combined dynamics and measurements on a system lead to a density matrix $\rho(t)$ with only diagonal elements in a given basis $\{|n\rangle\}$, it may be treated as a classical mixture, i.e., a system which randomly occupies the basis states $|n\rangle$ with probabilities $\rho_{nn}(t)$. Fully equivalent to so-called smoothing in classical probability theory, subsequent probing of the occupation of the states $|n\rangle$ improves our ability to retrodict what was the outcome of a projective state measurement at time $t$. Here, we show with experiments on a superconducting qubit that the smoothed probabilities do not, in the same way as the diagonal elements of $\rho$, permit a classical mixture interpretation of the state of the system at the past time $t$.


title: Probing Intergalactic Magnetic Fields with Simulations of Electromagnetic Cascades

abstract: We determine the effect of intergalactic magnetic fields on the distribution of high energy gamma rays by performing three-dimensional Monte Carlo simulations of the development of gamma-ray-induced electromagnetic cascades in the magnetized intergalactic medium. We employ the so-called 'Large Sphere Observer' method to efficiently simulate blazar gamma ray halos. We study magnetic fields with a Batchelor spectrum and with maximal left- and right-handed helicities. We also consider the case of sources whose jets are tilted with respect to the line of sight. We verify the formation of extended gamma ray halos around the source direction, and observe spiral-like patterns if the magnetic field is helical. We apply the $Q$-statistics to the simulated halos to extract their spiral nature and also propose an alternative method, the $S$-statistics. Both methods provide a quantative way to infer the helicity of the intervening magnetic fields from the morphology of individual blazar halos for magnetic field strengths $B \gtrsim 10^{-15}\,{\rm G}$ and magnetic coherence lengths $L_{\rm c} \gtrsim 100\,{\rm Mpc}$. We show that the $S$-statistics has a better performance than the $Q$-statistics when assessing magnetic helicity from the simulated halos.


title: Formal Definition of QoE Metrics

abstract: This technical report formally defines the QoE metrics which are introduced and discussed in the article "QoE Beyond the MOS: An In-Depth Look at QoE via Better Metrics and their Relation to MOS" by Tobias Ho{\ss}feld, Poul E. Heegaard, Martin Varela, Sebastian M\"oller, accepted for publication in the Springer journal "Quality and User Experience". Matlab scripts for computing the QoE metrics for given data sets are available in GitHub.


title: HATS-19b, HATS-20b, HATS-21b: Three Transiting Hot-Saturns Discovered by the HATSouth Survey

abstract: We report the discovery by the HATSouth exoplanet survey of three hot-Saturn transiting exoplanets: HATS-19b, HATS-20b, and HATS-21b. The planet host HATS-19 is a slightly evolved V = 13.0 G0 star with [Fe/H] = 0.240, a mass of 1.303 Msun, and a radius of 1.75 Rsun. HATS-19b is in an eccentric orbit (e = 0.30) around this star with an orbital period of 4.5697 days and has a mass of 0.427 Mjup and a highly inflated radius of 1.66 Rjup. The planet HATS-20b has a Saturn-like mass and radius of 0.273 Mjup and 0.776 Rjup respectively. It orbits the V = 13.8 G9V star HATS-20 (Ms = 0.910 Msun; Rs = 0.892 Rsun) with a period of 3.7993 days. Finally, HATS-21 is a V = 12.2 G4V star with [Fe/H] = 0.300, a mass of 1.080 Msun, and a radius of 1.021 Rsun. Its accompanying planet HATS-21b has a 3.5544-day orbital period, a mass of 0.332 Mjup, and a moderately inflated radius of 1.123 Rjup. With the addition of these three very different planets to the growing sample of hot-Saturns, we re-examine the relations between the observed giant planet radii, stellar irradiation, and host metallicity. We find a significant positive correlation between planet equilibrium temperature and radius, and a weak negative correlation between host metallicity and radius. To assess the relative influence of various physical parameters on observed planet radii, we train and fit models using Random Forest regression. We find that for hot-Saturns (0.1 < Mp < 0.5 Mjup), the planetary mass and equilibrium temperature play dominant roles in determining radii. For hot-Jupiters (0.5 < Mp < 2.0 Mjup), the most important parameter is equilibrium temperature alone. Finally, for irradiated higher-mass planets (Mp > 2.0 Mjup), we find that equilibrium temperature dominates in influence, with smaller contributions from planet mass and host metallicity.


title: Rotationally-Driven Fragmentation for the Formation of the Binary Protostellar System L1551 IRS 5

abstract: Either bulk rotation or local turbulence is widely invoked to drive fragmentation in collapsing cores so as to produce multiple star systems. Even when the two mechanisms predict different manners in which the stellar spins and orbits are aligned, subsequent internal or external interactions can drive multiple systems towards or away from alignment thus masking their formation process. Here, we demonstrate that the geometrical and dynamical relationship between the binary system and its surrounding bulk envelope provide the crucial distinction between fragmentation models. We find that the circumstellar disks of the binary protostellar system L1551 IRS 5 are closely parallel not just with each other but also with their surrounding flattened envelope. Measurements of the relative proper motion of the binary components spanning nearly 30 yr indicate an orbital motion in the same sense as the envelope rotation. Eliminating orbital solutions whereby the circumstellar disks would be tidally truncated to sizes smaller than are observed, the remaining solutions favor a circular or low-eccentricity orbit tilted by up to $\sim$25$^\circ$ from the circumstellar disks. Turbulence-driven fragmentation can generate local angular momentum to produce a coplanar binary system, but which bears no particular relationship with its surrounding envelope. Instead, the observed properties conform with predictions for rotationally-driven fragmentation. If the fragments were produced at different heights or on opposite sides of the midplane in the flattened central region of a rotating core, the resulting protostars would then exhibit circumstellar disks parallel with the surrounding envelope but tilted from the orbital plane as is observed.


title: Finite-energy pseudoholomorphic planes with multiple asymptotic limits

abstract: It's known from from work of Hofer, Wysocki, and Zehnder [1996] and Bourgeois [2002] that in a contact manifold equipped with either a nondegenerate or Morse-Bott contact form, a finite-energy pseudoholomorphic curve will be asymptotic at each of its non removable punctures to a single periodic orbit of the Reeb vector field and that the convergence is exponential. We provide examples here to show that this need not be the case if the contact form is degenerate. More specifically, we show that on any contact manifold $(M, \xi)$ with cooriented contact structure one can choose a contact form $\lambda$ with $\ker\lambda=\xi$ and a compatible complex structure $J$ on $\xi$ so that for the associated $\mathbb{R}$-invariant almost complex structure $\tilde J$ on $\mathbb{R}\times M$ there exist families of embedded finite-energy $\tilde J$-holomorphic cylinders and planes having embedded tori as limit sets.


title: Permutation Invariant Training of Deep Models for Speaker-Independent Multi-talker Speech Separation

abstract: We propose a novel deep learning model, which supports permutation invariant training (PIT), for speaker independent multi-talker speech separation, commonly known as the cocktail-party problem. Different from most of the prior arts that treat speech separation as a multi-class regression problem and the deep clustering technique that considers it a segmentation (or clustering) problem, our model optimizes for the separation regression error, ignoring the order of mixing sources. This strategy cleverly solves the long-lasting label permutation problem that has prevented progress on deep learning based techniques for speech separation. Experiments on the equal-energy mixing setup of a Danish corpus confirms the effectiveness of PIT. We believe improvements built upon PIT can eventually solve the cocktail-party problem and enable real-world adoption of, e.g., automatic meeting transcription and multi-party human-computer interaction, where overlapping speech is common.


title: Evaluating linear response in active systems with no perturbing field

abstract: We present a method for the evaluation of time-dependent linear response functions for systems of active Ornstein-Uhlenbeck particles from unperturbed simulations. The method is inspired by the Malliavin weights sampling method proposed by Warren and Allen [Phys. Rev. Lett. 109, 250601 (2012)] for systems of Brownian particles. We illustrate our method by evaluating a linear response function for a single active particle in an external harmonic potential. As an application, we calculate the time-dependent mobility function and an effective temperature, defined through the Einstein relation between the self-diffusion and mobility coefficients, for a system of active particles interacting via a screened-Coulomb potential. We find that this effective temperature decreases with increasing persistence time of the self-propulsion. Initially, for not too large persistence times, it changes rather slowly, but then it decreases markedly when the persistence length of the self-propelled motion becomes comparable with the particle size.


title: The field and Killing spinor equations of M-theory and type IIA/IIB supergravity in coordinate-free notation

abstract: We review the actions of the supergravity theory in eleven dimensions as well as the type IIA and IIB supergravities in ten dimensions and derive the bosonic equations of motion in a coordinate-free notation. We also consider the existence of supersymmetries and the associated generalized Killing spinor equations. The aim of this note is to serve as a formulary and make the equations of supergravity more easily accessible to mathematicians.


title: Amenability of coarse spaces and K-algebras

abstract: In this article we analyze the notions of amenability and paradoxical decomposition from an algebraic perspective. We consider this dichotomy for locally finite extended metric spaces and for general algebras over commutative fields. In the context of algebras we also study the relation of amenability with proper infiniteness. We apply our general analysis to two important classes of algebras: the unital Leavitt path algebras and the translation algebras on locally finite extended metric spaces. In particular, we show that the amenability of a metric space is equivalent to the algebraic amenability of the corresponding translation algebra.


title: Proactive Location-Based Scheduling of Delay-Constrained Traffic Over Fading Channels

abstract: In this paper, proactive resource allocation based on user location for point-to-point communication over fading channels is introduced, whereby the source must transmit a packet when the user requests it within a deadline of a single time slot. We introduce a prediction model in which the source predicts the request arrival $T_p$ slots ahead, where $T_p$ denotes the prediction window (PW) size. The source allocates energy to transmit some bits proactively for each time slot of the PW with the objective of reducing the transmission energy over the non-predictive case. The requests are predicted based on the user location utilizing the prior statistics about the user requests at each location. We also assume that the prediction is not perfect. We propose proactive scheduling policies to minimize the expected energy consumption required to transmit the requested packets under two different assumptions on the channel state information at the source. In the first scenario, offline scheduling, we assume the channel states are known a-priori at the source at the beginning of the PW. In the second scenario, online scheduling, it is assumed that the source has causal knowledge of the channel state. Numerical results are presented showing the gains achieved by using proactive scheduling policies compared with classical (reactive) networks. Simulation results also show that increasing the PW size leads to a significant reduction in the consumed transmission energy even with imperfect prediction.


title: Beyond LCDM: Bits join the Dark Side

abstract: We are encouraged to look beyond LambdaCDM as there are no satisfactory explanations for dark energy or dark matter. An evidence based approach applies two foundational principles of information theory to show dark energy and dark matter may be Holographic Dark Information Energy, HDIE. HDIE explains many effects attributed separately to Lambda and CDM. HDIE mimics Lambda with sufficient energy and an equation of state parameter, w= -1.03+-0.05 at redshifts z<1.35 to account for accelerating expansion. HDIE is clumped around galaxies at densities that distort space-time, explaining many CDM attributed effects. The present ratio of HDIE:baryons ~2.15, required for observed expansion history, is equivalent to a dark matter fraction ~68%, consistent with several galaxy surveys. The HDIE/baryon model is based largely on proven physics, provides a common explanation for dark energy and dark matter, and solves the cosmological coincidence problem. At earlier times, z > ~1.35, HDIE was phantom, w = -1.82+-0.08, enabling the model to be falsified. HDIE fits Planck dark energy wo-wa plots at least as well as Lambda, and is consistent with other results that suggest dark energy was phantom at earlier times. A new w-parameterisation is proposed, as the usual CPL parameterisation is unsuitable for distinguishing between HDIE/baryon and LambdaCDM models.


title: Machine-based Multimodal Pain Assessment Tool for Infants: A Review

abstract: Bedside caregivers assess infants' pain at constant intervals by observing specific behavioral and physiological signs of pain. This standard has two main limitations. The first limitation is the intermittent assessment of pain, which might lead to missing pain when the infants are left unattended. Second, it is inconsistent since it depends on the observer's subjective judgment and differs between observers. The intermittent and inconsistent assessment can induce poor treatment and, therefore, cause serious long-term consequences. To mitigate these limitations, the current standard can be augmented by an automated system that monitors infants continuously and provides quantitative and consistent assessment of pain. Several automated methods have been introduced to assess infants' pain automatically based on analysis of behavioral or physiological pain indicators. This paper comprehensively reviews the automated approaches (i.e., approaches to feature extraction) for analyzing infants' pain and the current efforts in automatic pain recognition. In addition, it reviews the databases available to the research community and discusses the current limitations of the automated pain assessment.


title: Short-Pulse Amplification by Strongly-Coupled Stimulated Brillouin Scattering

abstract: We examine the feasibility of strongly-coupled stimulated Brillouin scattering as a mechanism for the plasma-based amplification of sub-picosecond pulses. In particular, we use fluid theory and particle-in-cell simulations to compare the relative advantages of Raman and Brillouin amplification over a broad range of achievable parameters.


title: On SPDE and backward filtering equations for SDE systems (direct approach)

abstract: A direct approach to linear backward filtering equations for SDE systems is proposed. This preprint is a corrected version of the paper 1995 in the LMS Lecture Notes combined with another paper by the author on the direct approach to linear SPDEs for SDEs.


title: Massive scalar field evolution in de Sitter

abstract: The behaviour of a massive, non-interacting and non-minimally coupled quantised scalar field in an expanding de Sitter background is investigated by solving the field evolution for an arbitrary initial state. In this approach there is no need to choose a vacuum in order to provide a definition for particle states. We conclude that the expanding de Sitter space is a stable equilibrium configuration under small perturbations of the initial conditions. Depending on the initial state, the energy density can approach its asymptotic value from above or below, the latter of which implies a violation of the weak energy condition. The backreaction of the quantum corrections can therefore lead to a phase of super-acceleration also in the non-interacting massive case.


title: Understanding the pseudocapacitance of RuO2 from joint density functional theory

abstract: Pseudocapacitors have been experimentally studied for many years in electric energy storage. However, first principles understanding of the pseudocapacitive behavior is still not satisfactory due to the complexity involved in modeling electrochemistry. In this paper, we applied a novel simulation technique called Joint Density Functional Theory (JDFT) to simulate the pseudocapacitive behavior of RuO2, a prototypical material, in a model electrolyte. We obtained from JDFT a capacitive curve which showed a redox peak position comparable to that in the experimental cyclic voltammetry (CV) curve. We found that the experimental turning point from double-layer to pseudocapacitive charge storage at low scan rates could be explained by the hydrogen adsorption at low coverage. As the electrode voltage becomes more negative, H coverage increases and causes the surface structure change, leading to bended OH bonds at the on-top oxygen atoms and large capacitance. This H coverage-dependent capacitance can explain the high pseudocapacitance of hydrous RuO2. Our work here provides a first principles understanding of the pseudocapacitance for RuO2 in particular and for transition-metal oxides in general.


title: Two problems from the Polishchuk and Positselski book on Quadratic algebras

abstract: In the book 'Quadratic algebras' by Polishchuk and Positselski [23] algebras with a small number of generators (n=2,3) are considered. For some number r of relations possible Hilbert series are listed, and those appearing as series of Koszul algebras are specified. The first case, where it was not possible to do, namely the case of three generators n=3 and six relations r=6 is formulated as an open problem. We give here a complete answer to this question, namely for quadratic algebras with dim A_1=dim A_2=3, we list all possible Hilbert series, and find out which of them can come from Koszul algebras, and which can not. As a consequence of this classification, we found an algebra, which serves as a counterexample to another problem from the same book [23] (Chapter 7, Sec. 1, Conjecture 2), saying that Koszul algebra of finite global homological dimension d has dim A_1 >= d. Namely, the 3-generated algebra A given by relations xx+yx=xz=zy=0 is Koszul and its Koszul dual algebra A^! has Hilbert series of degree 4: H_{A^!}(t)= 1+3t+3t^2+2t^3+t^4, hence A has global homological dimension 4.


title: Late Time Emission of Prompt Fission Gamma Rays

abstract: The emission of prompt fission $\gamma$ rays within a few nanoseconds to a few microseconds following the scission point is studied in the Hauser-Feshbach formalism applied to the deexcitation of primary excited fission fragments. Neutron and $\gamma$-ray evaporations from fully accelerated fission fragments are calculated in competition at each stage of the decay, and the role of isomers in the fission products, before $\beta$-decay, is analyzed. The time evolution of the average total $\gamma$-ray energy, average total $\gamma$-ray multiplicity, and fragment-specific $\gamma$-ray spectra, is presented in the case of neutron-induced fission reactions of $^{235}$U and $^{239}$Pu, as well as spontaneous fission of $^{252}$Cf. The production of specific isomeric states is calculated and compared to available experimental data. About 7% of all prompt fission $\gamma$ rays are predicted to be emitted between 10 nsec and 5 $\mu$sec following fission, in the case of $^{235}$U and $^{239}$Pu $(n_{\rm th},f)$ reactions, and up to 3% in the case of $^{252}$Cf spontaneous fission. The cumulative average total $\gamma$-ray energy increases by 2 to 5% in the same time interval. Finally, those results are shown to be robust against significant changes in the model input parameters.


title: Convergence properties of $\eta\to 3\pi$ decays in chiral perturbation theory

abstract: Theoretical efforts to describe and explain the $\eta\to 3\pi$ decays reach far back in time. Even today, the convergence of the decay widths and some of the Dalitz plot parameters seems problematic in low energy QCD. In the framework of resummed CHPT, we explore the question of compatibility of experimental data with a reasonable convergence of a carefully defined chiral series, where NNLO remainders are assumed to be small. By treating the uncertainties in the higher orders statistically, we numerically generate a large set of theoretical predictions, which are then confronted with experimental information. In the case of the decay widths, the experimental values can be reconstructed for a reasonable range of the free parameters and thus no tension is observed, in spite of what some of the traditional calculations suggest. The Dalitz plot parameters $a$ and $d$ can be described very well too. When the parameters $b$ and $\alpha$ are concerned, we find a mild tension for the whole range of the free parameters, at less than 2$\sigma$ C.L. This can be interpreted in two ways - either some of the higher order corrections are indeed unexpectedly large or there is a specific configuration of the remainders, which is, however, not completely improbable. Also, the distribution of the theoretical uncertainties is found to be significantly non-gaussian, so the consistency cannot be simply judged by the 1$\sigma$ error bars.


title: Orbit Portraits in non-Autonomous Iteration

abstract: We extend the definition of an orbit portrait to the context of non-autonomous iteration, both for the combinatorial version involving collections of angles and for the dynamic version involving external rays where combinatorial portraits can be realized by the dynamics associated with sequences of polynomials with suitably uniformly bounded degrees and coefficients. We show that, in the case of sequences of polynomials of constant degree, the portraits which arise are eventually periodic which is somewhat similar to the classical theory of polynomial iteration. However, if the degrees of the polynomials in the sequence are allowed to vary, one can obtain portraits with complementary arcs of irrational length which are fundamentally different from the classical ones.


title: Elastic Graphs

abstract: An elastic graph is a graph with an elasticity associated to each edge. It may be viewed as a network made out of ideal rubber bands. If the rubber bands are stretched on a target space there is an elastic energy. We characterize when a homotopy class of maps from one elastic graph to another is loosening, i.e., decreases this elastic energy for all possible targets. This fits into a more general framework of energies for maps between graphs.


title: Role of spin-transfer torques on synchronization and resonance phenomena in stochastic magnetic oscillators

abstract: A theoretical study on how synchronization and resonance-like phenomena in superparamagnetic tunnel junctions can be driven by spin-transfer torques is presented. We examine the magnetization of a superparamagnetic free layer that reverses randomly between two well-defined orientations due to thermal fluctuations, acting as a stochastic oscillator. When subject to an external ac forcing this system can present stochastic resonance and noise-enhanced synchronization. We focus on the roles of the mutually perpendicular damping-like and field-like torques, showing that the response of the system is very different at low and high-frequencies. We also demonstrate that the field-like torque can increase the efficiency of the current-driven forcing, specially at sub-threshold electric currents. These results can be useful for possible low-power, more energy efficient, applications.


title: Beyond the standard gauging: gauge symmetries of Dirac Sigma Models

abstract: In this paper we study the general conditions that have to be met for a gauged extension of a two-dimensional bosonic sigma-model to exist. In an inversion of the usual approach of identifying a global symmetry and then promoting it to a local one, we focus directly on the gauge symmetries of the theory. This allows for action functionals which are gauge invariant for rather general background fields in the sense that their invariance conditions are milder than the usual case. In particular, the vector fields that control the gauging need not be Killing. The relaxation of isometry for the background fields is controlled by two connections on a Lie algebroid L in which the gauge fields take values, in a generalization of the common Lie-algebraic picture. Here we show that these connections can always be determined when L is a Dirac structure in the H-twisted Courant algebroid. This also leads us to a derivation of the general form for the gauge symmetries of a wide class of two-dimensional topological field theories called Dirac sigma-models, which interpolate between the G/G Wess-Zumino-Witten model and the (Wess-Zumino-term twisted) Poisson sigma model.


title: NuSTAR Observations of the Black Hole GS 1354-645: Evidence of Rapid Black Hole Spin

abstract: We present the results of a NuSTAR study of the dynamically confirmed stellar-mass black hole GS 1354-645. The source was observed during its 2015 "hard" state outburst; we concentrate on spectra from two relatively bright phases. In the higher-flux observation, the broadband NuSTAR spectra reveal a clear, strong disk reflection spectrum, blurred by a degree that requires a black hole spin of a = cJ/GM^2 > 0.98 (1 sigma statistical limits only). The fits also require a high inclination: theta = 75(2) degrees. Strong "dips" are sometimes observed in the X-ray light curves of sources viewed at such an angle; these are absent, perhaps indicating that dips correspond to flared disk structures that only manifest at higher accretion rates. In the lower-flux observation, there is evidence of radial truncation of the thin accretion disk. We discuss these results in the context of spin in stellar-mass black holes, and inner accretion flow geometries at moderate accretion rates.


title: Today a Duo, But Once a Trio? The Double White Dwarf HS 2220$+$2146 May Be A Post-Blue Straggler Binary

abstract: For sufficiently wide orbital separations {\it a}, the two members of a stellar binary evolve independently. This implies that in a wide double white dwarf (DWD), the more massive WD should always be produced first, when its more massive progenitor ends its main-sequence life, and should therefore be older and cooler than its companion. The bound, wide DWD HS 2220$+$2146 ($a\approx500$ AU) does not conform to this picture: the more massive WD is the younger, hotter of the pair. We show that this discrepancy is unlikely to be due to past mass-transfer phases or to the presence of an unresolved companion. Instead, we propose that HS 2220$+$2146 formed through a new wide DWD evolutionary channel involving the merger of the inner binary in a hierarchical triple system. The resulting blue straggler and its wide companion then evolved independently, forming the WD pair seen today. Although we cannot rule out other scenarios, the most likely formation channel has the inner binary merging while both stars are still on the main sequence. This provides us with the tantalizing possibility that Kozai-Lidov oscillations may have played a role in the inner binary's merger. {\it Gaia} may uncover hundreds more wide DWDs, leading to the identification of other systems like HS 2220$+$2146. There are already indications that other WD systems may have been formed through different, but related, hierarchical triple evolutionary scenarios. Characterizing these populations may allow for thorough testing of the efficiency with which KL oscillations induce stellar mergers.


title: Convergence Rate of Frank-Wolfe for Non-Convex Objectives

abstract: We give a simple proof that the Frank-Wolfe algorithm obtains a stationary point at a rate of $O(1/\sqrt{t})$ on non-convex objectives with a Lipschitz continuous gradient. Our analysis is affine invariant and is the first, to the best of our knowledge, giving a similar rate to what was already proven for projected gradient methods (though on slightly different measures of stationarity).


title: Distributed-memory Hierarchical Interpolative Factorization

abstract: The hierarchical interpolative factorization (HIF) offers an efficient way for solving or preconditioning elliptic partial differential equations. By exploiting locality and low-rank properties of the operators, the HIF achieves quasi-linear complexity for factorizing the discrete positive definite elliptic operator and linear complexity for solving the associated linear system. In this paper, the distributed-memory HIF (DHIF) is introduced as a parallel and distributed-memory implementation of the HIF. The DHIF organizes the processes in a hierarchical structure and keep the communication as local as possible. The computation complexity is $O\left(\frac{N\log N}{P}\right)$ and $O\left(\frac{N}{P}\right)$ for constructing and applying the DHIF, respectively, where $N$ is the size of the problem and $P$ is the number of processes. The communication complexity is $O\left(\sqrt{P}\log^3 P\right)\alpha + O\left(\frac{N^{2/3}}{\sqrt{P}}\right)\beta$ where $\alpha$ is the latency and $\beta$ is the inverse bandwidth. Extensive numerical examples are performed on the NERSC Edison system with up to 8192 processes. The numerical results agree with the complexity analysis and demonstrate the efficiency and scalability of the DHIF.


title: Colorful simplicial depth, Minkowski sums, and generalized Gale transforms

abstract: The colorful simplicial depth of a collection of d+1 finite sets of points in Euclidean d-space is the number of choices of a point from each set such that the origin is contained in their convex hull. We use methods from combinatorial topology to prove a tight upper bound on the colorful simplicial depth. This implies a conjecture of Deza et al. (2006). Furthermore, we introduce colorful Gale transforms as a bridge between colorful configurations and Minkowski sums. Our colorful upper bound then yields a tight upper bound on the number of totally mixed facets of certain Minkowski sums of simplices. This resolves a conjecture of Burton (2003) in the theory of normal surfaces.


title: Gauged B-L Number and Neutron--Antineutron Oscillation: Long-range Forces Mediated by Baryophotons

abstract: Transformation of neutron to antineutron is a small effect that has not yet been experimentally observed. %\cite{Phillips:2014fgb}. In principle, it can occur with free neutrons in the vacuum or with bound neutrons inside the nuclear environment different for neutrons and antineutrons and for that reason in the latter case it is heavily suppressed. Free neutron transformation also can be suppressed if environmental vector field exists destinguishing neutron from antineutron. We consider here the case of a vector field coupled to $B-L$ charge of the particles ($B-L$ photons) and study a possibility of this to lead to the observable suppression of neutron to antineutron transformation. The suppression effect however can be removed by applying external magnetic field. If the neutron--antineutron oscillation will be discovered in free neutron oscillation experiments, this will imply limits on $B-L$ photon coupling constant and interaction radius few order of magnitudes stronger than present limits form the tests of the equivalence principle. If $n-\bar n$ oscillation will be discovered via nuclear instability, but not in free neutron oscillations in corresponding level, this would indicate to the presence of fifth-forces mediated by such baryophotons.


title: Painlev\'e analysis for two 1D parabolic-parabolic models of chemotaxis; some travelling wave solutions

abstract: In this paper we study the Painlev\'e analysis for two models of chemotaxis. We find that in some cases the reductions of these models in terms of travelling wave variable allow exact analytical solutions.


title: Quantum solvable models with nonlocal one point interactions

abstract: Within the framework of quantum mechanics working with one-dimensional, manifestly non-Hermitian Hamiltonians $H=T+V$ the traditional class of the exactly solvable models with local point interactions $V=V(x)$ is generalized. The consequences of the use of the nonlocal point interactions such that $(V f)(x) = \int K(x,s) f(s) ds$ are discussed using the suitably adapted formalism of boundary triplets.


title: A Comparison of Preconditioned Krylov Subspace Methods for Large-Scale Nonsymmetric Linear Systems

abstract: Preconditioned Krylov subspace (KSP) methods are widely used for solving large-scale sparse linear systems arising from numerical solutions of partial differential equations (PDEs). These linear systems are often nonsymmetric due to the nature of the PDEs, boundary or jump conditions, or discretization methods. While implementations of preconditioned KSP methods are usually readily available, it is unclear to users which methods are the best for different classes of problems. In this work, we present a comparison of some KSP methods, including GMRES, TFQMR, BiCGSTAB, and QMRCGSTAB, coupled with three classes of preconditioners, namely Gauss-Seidel, incomplete LU factorization (including ILUT, ILUTP, and multilevel ILU), and algebraic multigrid (including BoomerAMG and ML). Theoretically, we compare the mathematical formulations and operation counts of these methods. Empirically, we compare the convergence and serial performance for a range of benchmark problems from numerical PDEs in 2D and 3D with up to millions of unknowns and also assess the asymptotic complexity of the methods as the number of unknowns increases. Our results show that GMRES tends to deliver better performance when coupled with an effective multigrid preconditioner, but it is less competitive with an ineffective preconditioner due to restarts. BoomerAMG with proper choice of coarsening and interpolation techniques typically converges faster than ML, but both may fail for ill-conditioned or saddle-point problems while multilevel ILU tends to succeed. We also show that right preconditioning is more desirable. This study helps establish some practical guidelines for choosing preconditioned KSP methods and motivates the development of more effective preconditioners.


title: Minimizing statistical and systematic bias in transverse momentum correlations for relativistic heavy-ion collisions

abstract: Two-particle correlation measurements and analysis are an important component of the relativistic heavy-ion physics program. In particular, particle pair-number correlations on two-dimensional transverse momentum ($p_t$) allow unique access to soft, semi-hard and hard-scattering processes in these collisions. Precise measurements of this type of correlation are essential for understanding the dynamics in heavy-ion collisions. However, transverse momentum correlation measurements are especially vulnerable to statistical and systematic biases. In this paper the origins of these large bias effects are explained and mathematical correlation forms are derived from mean-$p_t$ fluctuation quantities in the literature. Monte Carlo simulations are then used to determine the conditions, e.g. multiplicity and collision centrality bin widths, where each correlation form is minimally biased. The ranges of applicability for each correlation quantity are compared. Several are found to reproduce the assumed input correlations with reasonable fidelity over a wide range of conditions encountered in practical analysis of data.


title: To fission or not to fission

abstract: The fission-fragments mass-yield of 236U is obtained by an approximate solution of the eigenvalue problem of the collective Hamiltonian that describes the dynamics of the fission process whose degrees of freedom are: the fission (elongation), the neck and the mass-asymmetry mode. The macroscopic-microscopic method is used to evaluate the potential energy surface. The macroscopic energy part is calculated using the liquid drop model and the microscopic corrections are obtained using the Woods-Saxon single-particle levels. The four dimensional modified Cassini ovals shape parametrization is used to describe the shape of the fissioning nucleus. The mass tensor is taken within the cranking-type approximation. The final fragment mass distribution is obtained by weighting the adiabatic density distribution in the collective space with the neck-dependent fission probability. The neck degree of freedom is found to play a significant role in determining that final fragment mass distribution.


title: STAM: A Framework for Spatio-Temporal Affordance Maps

abstract: Affordances have been introduced in literature as action opportunities that objects offer, and used in robotics to semantically represent their interconnection. However, when considering an environment instead of an object, the problem becomes more complex due to the dynamism of its state. To tackle this issue, we introduce the concept of Spatio-Temporal Affordances (STA) and Spatio-Temporal Affordance Map (STAM). Using this formalism, we encode action semantics related to the environment to improve task execution capabilities of an autonomous robot. We experimentally validate our approach to support the execution of robot tasks by showing that affordances encode accurate semantics of the environment.


title: A note on some inequalities used in channel polarization and polar coding

abstract: We give a unified treatment of some inequalities that are used in the proofs of channel polarization theorems involving a binary-input discrete memoryless channel.