title: Weyl's polarization theorem in positive characteristic

abstract: Let $V$ be an $n$-dimensional algebraic representation over an algebraically closed field $K$ of a group $G$. For $m > 0$, we study the invariant rings $K[V^{ m}]^G$ for the diagonal action of $G$ on $V^m$. In characteristic zero, a theorem of Weyl tells us that we can obtain all the invariants in $K[V^m]^G$ by the process of polarization and restitution from $K[V^n]^G$. In particular, this means that if $K[V^n]^G$ is generated in degree $\leq d$, then so is $K[V^m]^G$ no matter how large $m$ is. There are several explicit counterexamples to Weyl's theorem in positive characteristic. However, when $G$ is a (connected) reductive affine group scheme over $\mathbb{Z}$ and $V^*$ is a good $G$-module, we show that Weyl's theorem holds in sufficiently large characteristic. As a special case, we consider the ring of invariants $R(n,m)$ for the left-right action of ${\rm SL}_n \times {\rm SL}_n$ on $m$-tuples of $n \times n$ matrices. In this case, we show that the invariants of degree $\leq n^6$ suffice to generate $R(n,m)$ if the characteristic is larger than $2n^6 + n^2$.


title: Learning local, quenched disorder in plasticity and other crackling noise phenomena

abstract: When far from equilibrium, many-body systems display behavior that strongly depends on the initial conditions. A characteristic such example is the phenomenon of plasticity of crystalline and amorphous materials that strongly depends on the material history. In plasticity modeling, the history is captured by a quenched, local and disordered flow stress distribution. While it is this disorder that causes avalanches that are commonly observed during nanoscale plastic deformation, the functional form and scaling properties have remained elusive. In this paper, a generic formalism is developed for deriving local disorder distributions from field-response (e.g. stress/strain) timeseries in models of crackling noise. We demonstrate the efficiency of the method in the hysteretic random-field Ising model and also, models of elastic interface depinning that have been used to model crystalline and amorphous plasticity. We show that the capacity to resolve the quenched disorder distribution improves with the temporal resolution and number of samples.

categories: cond-mat.mtrl-sci cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech nlin.CG


title: SDSS-IV MaNGA: The Spectroscopic Discovery of Strongly Lensed Galaxies

abstract: We present a catalogue of 38 spectroscopically detected strong galaxy-galaxy gravitational lens candidates identified in the Sloan Digital Sky Survey IV (SDSS-IV). We were able to simulate narrow-band images for 8 of them demonstrating evidence of multiple images. Two of our systems are compound lens candidates, each with 2 background source-planes. One of these compound systems shows clear lensing features in the narrow-band image. Our sample is based on 2812 galaxies observed by the Mapping Nearby Galaxies at APO (MaNGA) integral field unit (IFU). This Spectroscopic Identification of Lensing Objects (SILO) survey extends the methodology of the Sloan Lens ACS Survey (SLACS) and BOSS Emission-Line Survey (BELLS) to lower redshift and multiple IFU spectra. We searched ~ 1.5 million spectra, of which 3065 contained multiple high signal-to-noise background emission-lines or a resolved [OII] doublet, that are included in this catalogue. Upon manual inspection, we discovered regions with multiple spectra containing background emission-lines at the same redshift, providing evidence of a common source-plane geometry which was not possible in previous SLACS and BELLS discovery programs. We estimate more than half of our candidates have an Einstein radius > 1.7", which is significantly greater than seen in SLACS and BELLS. These larger Einstein radii produce more extended images of the background galaxy increasing the probability that a background emission-line will enter one of the IFU spectroscopic fibres, making detection more likely.


abstract: We show that the possible Cantor-Bendixson ranks of countable SFTs are exactly the finite ordinals and ordinals of the form $\lambda + 3$, where $\lambda$ is a computable ordinal. This result was claimed by the author in his PhD dissertation, but the proof contains an error, which is fixed in this note.


title: A Simple proof of Johnson-Lindenstrauss extension theorem

abstract: Johnson and Lindenstrauss proved that any Lipschitz mapping from an $n$-point subset of a metric space into Hilbert space can be extended to the whole space, while increasing the Lipschitz constant by a factor of $O(\sqrt{\log n})$. We present a simplification of their argument that avoids dimension reduction and the Kirszbraun theorem.


title: On Generation of Adversarial Examples using Convex Programming

abstract: It has been observed that deep learning architectures tend to make erroneous decisions with high reliability for particularly designed adversarial instances. In this work, we show that the perturbation analysis of these architectures provides a framework for generating adversarial instances by convex programming which, for classification tasks, is able to recover variants of existing non-adaptive adversarial methods. The proposed framework can be used for the design of adversarial noise under various desirable constraints and different types of networks. Moreover, this framework is capable of explaining various existing adversarial methods and can be used to derive new algorithms as well. We make use of these results to obtain novel algorithms. The experiments show the competitive performance of the obtained solutions, in terms of fooling ratio, when benchmarked with well-known adversarial methods.


title: On The Uplink Throughput of Zero-Forcing in Cell-Free Massive MIMO with Coarse Quantization

abstract: The recently proposed Cell-Free massive MIMO architecture is studied for the uplink. In contrast to most previous works, joint detection is performed using global CSI. Therefore, we study strategies for transferring CSI to the CPU taking into account the fronthaul capacity which limits CSI quantization. Two strategies for pilot-based CSI acquisition are considered: estimate-and-quantize and quantize-and-estimate. These are analysed using the Bussgang decomposition. For a given quantization constraint for the data and CSI the achievable rate per user with Zero-Forcing is determined. Numerical results show that quantize-and-estimate (the simpler strategy) is similar to or better than estimate-and-quantize, especially for 1-bit resolution.


title: Thermoelectric performance of topological boundary modes

abstract: We investigate quantum transport and thermoelectrical properties of a finite-size Su-Schrieffer-Heeger model, a paradigmatic model for a one-dimensional topological insulator, which displays topologically protected edge states. By coupling the model to two fermionic reservoirs at its ends, we can explore the non-equilibrium dynamics of the system. Investigating the energy-resolved transmission, the current and the noise, we find that these observables can be used to detect the topologically non-trivial phase. With specific parameters and asymmetric reservoir coupling strengths, we show that we can dissipatively prepare the edge states as stationary states of a non-equilibrium configuration. In addition, we point out that the edge states can be exploited to design a refrigerator driven by chemical work or a heat engine driven by a thermal gradient, respectively. These thermal devices do not require asymmetric couplings and are topologically protected against symmetry-preserving perturbations. Their maximum efficiencies significantly exceed that of a single quantum dot device at comparable coupling strengths.


title: Random Access Schemes in Wireless Systems With Correlated User Activity

abstract: Traditional random access schemes are designed based on the aggregate process of user activation, which is created on the basis of independent activations of the users. However, in Machine-Type Communications (MTC), some users are likely to exhibit a high degree of correlation, e.g. because they observe the same physical phenomenon. This paves the way to devise access schemes that combine scheduling and random access, which is the topic of this work. The underlying idea is to schedule highly correlated users in such a way that their transmissions are less likely to result in a collision. To this end, we propose two greedy allocation algorithms. Both attempt to maximize the throughput using only pairwise correlations, but they rely on different assumptions about the higher-order dependencies. We show that both algorithms achieve higher throughput compared to the traditional random access schemes, suggesting that user correlation can be utilized effectively in access protocols for MTC.


title: The Trade-off between Privacy and Fidelity via Ehrhart Theory

abstract: As an increasing amount of data is gathered nowadays and stored in databases (DBs), the question arises of how to protect the privacy of individual records in a DB even while providing accurate answers to queries on the DB. Differential Privacy (DP) has gained acceptance as a framework to quantify vulnerability of algorithms to privacy breaches. We consider the problem of how to sanitize an entire DB via a DP mechanism, on which unlimited further querying is performed. While protecting privacy, it is important that the sanitized DB still provide accurate responses to queries. The central contribution of this work is to characterize the amount of information preserved in an optimal DP DB sanitizing mechanism (DSM). We precisely characterize the utility-privacy trade-off of mechanisms that sanitize DBs in the asymptotic regime of large DBs. We study this in an information-theoretic framework by modeling a generic distribution on the data, and a measure of fidelity between the histograms of the original and sanitized DBs. We consider the popular $\mathbb{L}_{1}-$distortion metric that leads to the formulation as a linear program (LP). This optimization problem is prohibitive in complexity with the number of constraints growing exponentially in the parameters of the problem. Leveraging tools from discrete geometry, analytic combinatorics, and duality theorems of optimization, we fully characterize the optimal solution in terms of a power series whose coefficients are the number of integer points on a multidimensional convex polytope studied by Ehrhart in 1967. Employing Ehrhart theory, we determine a simple closed form computable expression for the asymptotic growth of the optimal privacy-fidelity trade-off to infinite precision. At the heart of the findings is a deep connection between the minimum expected distortion and the Ehrhart series of an integral convex polytope.


title: Supersymmetric integrable theories from no particle production

abstract: We consider a theory of scalar superfields in two dimensions with arbitrary superpotential. By imposing no particle production in tree level scattering, we constrain the form of the admissible interactions, recovering a supersymmetric extension of the sinh-Gordon model.


title: Detecting Adversarial Examples - A Lesson from Multimedia Forensics

abstract: Adversarial classification is the task of performing robust classification in the presence of a strategic attacker. Originating from information hiding and multimedia forensics, adversarial classification recently received a lot of attention in a broader security context. In the domain of machine learning-based image classification, adversarial classification can be interpreted as detecting so-called adversarial examples, which are slightly altered versions of benign images. They are specifically crafted to be misclassified with a very high probability by the classifier under attack. Neural networks, which dominate among modern image classifiers, have been shown to be especially vulnerable to these adversarial examples. However, detecting subtle changes in digital images has always been the goal of multimedia forensics and steganalysis. In this paper, we highlight the parallels between these two fields and secure machine learning. Furthermore, we adapt a linear filter, similar to early steganalysis methods, to detect adversarial examples that are generated with the projected gradient descent (PGD) method, the state-of-the-art algorithm for this task. We test our method on the MNIST database and show for several parameter combinations of PGD that our method can reliably detect adversarial examples. Additionally, the combination of adversarial re-training and our detection method effectively reduces the attack surface of attacks against neural networks. Thus, we conclude that adversarial examples for image classification possibly do not withstand detection methods from steganalysis, and future work should explore the effectiveness of known techniques from multimedia forensics in other adversarial settings.


title: Divergence-Optimal Fixed-to-Fixed Length Distribution Matching With Shell Mapping

abstract: Distribution matching (DM) transforms independent and Bernoulli(1/2) distributed bits into a sequence of output symbols with a desired distribution. A fixed-to-fixed length, invertible DM architecture based on shell mapping is presented. It is shown that shell mapping for distribution matching (SMDM) is the optimum DM for the informational divergence metric and that finding energy optimal sequences is a special case of divergence minimization. Additionally, it is shown how to find the required shell mapping weight function to approximate arbitrary output distributions. SMDM is combined with probabilistic amplitude shaping (PAS) to operate close to the Shannon limit. SMDM exhibits excellent performance for short blocklengths as required by ultra-reliable low-latency (URLLC) applications. SMDM outperforms constant composition DM (CCDM) by 0.6 dB when used with 64-QAM at a spectral efficiency of 3 bits/channel use and a 5G low-density parity-check code with a short blocklength of 192 bits


title: Ho\v{r}ava-Lifshitz cosmology in light of new data

abstract: We present new observational constraints on Lorentz violating Ho\v{r}ava-Lifshitz cosmological scenarios using an updated cosmological data set from Cosmic Microwave Background (Planck CMB), expansion rates of elliptical and lenticular galaxies, JLA compilation (Joint Light-Curve Analysis) data for Type Ia supernovae (SneIa), Baryon Acoustic Oscillations (BAO) and priors on the Hubble parameter with an alternative parametrisation of the equations. Unlike in other approaches we consider the curvature parameter $\Omega_k$ as a free parameter in the analysis we considered the parameters $\Omega_k$ and $\Delta N_\nu$ as completely free, which helped to place new, updated bounds on several of the theory parameters. Remarkably, the detailed balance scenario exhibits positive spatial curvature to more than 3$\sigma$, whereas for further theory generalizations we found evidence for positive spatial curvature at 1$\sigma$. This could create circumstantial evidence from observations and could be used to single out distinct formulations and scenarios.


title: A Dynamic Jamming Game for Real-Time Status Updates

abstract: We study timely status updates of a real-time system in an adversarial setting. The system samples a physical process, and sends the samples from the source (e.g., a sensor) to the destination (e.g, a control center) through a channel. For real-time monitoring/control tasks, it is crucial for the system to update the status of the physical process "timely". We measure the timeliness of status updates by the time elapsed since the latest update at the destination was generated at the source, and define the time elapsed as age of information, or age in short. To sabotage the system, an attacker aims to maximize the age by jamming the channel and hence causing delay in status updates. The system aims to minimize the age by judiciously choosing when to sample and send the updates. We model the ongoing repeated interaction between the attacker and the system as a dynamic game. In each stage game, the attacker chooses the jamming time according to the jamming time distribution, and the system responds by delaying the sampling according to the sampling policy. We prove that there exists a unique stationary equilibrium in the game, and provide a complete analytical characterization of the equilibrium. Our results shed lights on how the attacker sabotages the system and how the system should defend against the attacker.


title: Evidence for cascaded third harmonic generation in non-centrosymmetric gold nanoantennas

abstract: The optimization of nonlinear optical processes at the nanoscale is a crucial step for the development of nanoscale photon sources for quantum-optical networks. The development of innovative plasmonic nanoantenna designs and hybrid nanostructures to enhance optical nonlinearities in very small volumes represents one of the most promising routes. In such systems, the upconversion of photons can be achieved with high efficiencies via third-order processes, such as third harmonic generation (THG), thanks to the resonantly-enhanced volume currents. Conversely, second-order processes, such as second harmonic generation (SHG), are often inhibited by the symmetry of metal lattices and of common nanoantenna geometries. SHG and THG processes in plasmonic nanostructures are generally treated independently, since they both represent a small perturbation in the light-matter interaction mechanisms. In this work, we demonstrate that this paradigm does not hold in general, by providing evidence of a cascaded process in THG, which is fueled by SHG and sizably contributes to the overall yield. We address this mechanism by unveiling an anomalous fingerprint in the polarization state of the nonlinear emission from non-centrosymmetric gold nanoantennas and point out that such cascaded processes may also appear for structures that exhibit only moderate SHG yields - signifying its general relevance in plasmon-enhanced nonlinear optics. The presence of this peculiar mechanism in THG from plasmonic nanoantennas at telecommunication wavelengths allows gaining further insight on the physics of plasmon-enhanced nonlinear optical processes. This could be crucial in the realization of nanoscale elements for photon conversion and manipulation operating at room-temperature.


title: A Ranking Method For Selection Of $\Eta$ Mesons In High Multiplicity Events

abstract: The selection of $\eta$ mesons with a high efficiency and a high purity can be important in the formation of statistically significant invariant mass spectra in the reconstruction of short-lived particles such as $\eta' \rightarrow \pi^{+} \pi^{-} \eta$. In this study, a cut-based standard method and a Ranking method to reduce combinatorial background in the reconstruction of $\eta \rightarrow \gamma \gamma$ decays in high multiplicity hadronic events are presented. By using recorded ALEPH data and fully simulated events, the performances of the methods are compared. Results show that the Ranking method yields significant improvements in the purity of the selected $\eta$ meson relative to the standard method.


title: Signature of non-isotropic distribution of stellar rotation inclination angles in the Praesepe cluster

abstract: The distribution of the stellar rotation axes of 113 main sequence stars in the open cluster Praesepe are examined by using current photometric rotation periods, spectroscopic rotation velocities, and estimated stellar radii. Three different samples of stellar rotation data on spotted stars from the Galactic field and two independent samples of planetary hosts are used as control samples to support the consistency of the analysis. Considering the high completeness of the Praesepe sample and the behavior of the control samples, we find that the main sequence F - K stars in this cluster are susceptible to rotational axis alignment. Using a cone model, the most likely inclination angle is 76+/-14 degrees with a half opening angle of 47+/-24 degrees. Non-isotropic distribution of the inclination angles is preferred over the isotropic distribution, except if the rotation velocities used in this work are systematically overestimated. We found no indication of this being the case on the basis of the currently available data.


title: Stable and Consistent Membership at Scale with Rapid

abstract: We present the design and evaluation of Rapid, a distributed membership service. At Rapid's core is a scheme for multi-process cut detection (CD) that revolves around two key insights: (i) it suspects a failure of a process only after alerts arrive from multiple sources, and (ii) when a group of processes experience problems, it detects failures of the entire group, rather than conclude about each process individually. Implementing these insights translates into a simple membership algorithm with low communication overhead. We present evidence that our strategy suffices to drive unanimous detection almost-everywhere, even when complex network conditions arise, such as one-way reachability problems, firewall misconfigurations, and high packet loss. Furthermore, we present both empirical evidence and analyses that proves that the almost-everywhere detection happens with high probability. To complete the design, Rapid contains a leaderless consensus protocol that converts multi-process cut detections into a view-change decision. The resulting membership service works both in fully decentralized as well as logically centralized modes. We present an evaluation of Rapid in moderately scalable cloud settings. Rapid bootstraps 2000 node clusters 2-5.8x faster than prevailing tools such as Memberlist and ZooKeeper, remains stable in face of complex failure scenarios, and provides strong consistency guarantees. It is easy to integrate Rapid into existing distributed applications, of which we demonstrate two.


title: Approximate Randomized Benchmarking for Finite Groups

abstract: We investigate randomized benchmarking in a general setting with quantum gates that form a representation, not necessarily an irreducible one, of a finite group. We derive an estimate for the average fidelity, to which experimental data may then be calibrated. Furthermore, we establish that randomized benchmarking can be achieved by the sole implementation of quantum gates that generate the group as well as one additional arbitrary group element. In this case, we need to assume that the noise is close to being covariant. This yields a more practical approach to randomized benchmarking. Moreover, we show that randomized benchmarking is stable with respect to approximate Haar sampling for the sequences of gates. This opens up the possibility of using Markov chain Monte Carlo methods to obtain the random sequences of gates more efficiently. We demonstrate these results numerically using the well-studied example of the Clifford group as well as the group of monomial unitary matrices. For the latter, we focus on the subgroup with nonzero entries consisting of n-th roots of unity, which contains T gates.


title: Virtual Network Embedding Approximations: Leveraging Randomized Rounding

abstract: The Virtual Network Embedding Problem (VNEP) captures the essence of many resource allocation problems of today's infrastructure providers, which offer their physical computation and networking resources to customers. Customers request resources in the form of Virtual Networks, i.e. as a directed graph which specifies computational requirements at the nodes and communication requirements on the edges. An embedding of a Virtual Network on the shared physical infrastructure is the joint mapping of (virtual) nodes to physical servers together with the mapping of (virtual) edges onto paths in the physical network connecting the respective servers. This work initiates the study of approximation algorithms for the VNEP. Concretely, we study the offline setting with admission control: given multiple request graphs the task is to embed the most profitable subset while not exceeding resource capacities. Our approximation is based on the randomized rounding of Linear Programming (LP) solutions. Interestingly, we uncover that the standard LP formulation for the VNEP exhibits an inherent structural deficit when considering general virtual network topologies: its solutions cannot be decomposed into valid embeddings. In turn, focusing on the class of cactus request graphs, we devise a novel LP formulation, whose solutions can be decomposed into convex combinations of valid embedding. Proving performance guarantees of our rounding scheme, we obtain the first approximation algorithm for the VNEP in the resource augmentation model. We propose two types of rounding heuristics and evaluate their performance in an extensive computational study. Our results indicate that randomized rounding can yield good solutions (even without augmentations). Specifically, heuristic rounding achieves 73.8% of the baseline's profit, while not exceeding capacities.


title: Hourly-Similarity Based Solar Forecasting Using Multi-Model Machine Learning Blending

abstract: With the increasing penetration of solar power into power systems, forecasting becomes critical in power system operations. In this paper, an hourly-similarity (HS) based method is developed for 1-hour-ahead (1HA) global horizontal irradiance (GHI) forecasting. This developed method utilizes diurnal patterns, statistical distinctions between different hours, and hourly similarities in solar data to improve the forecasting accuracy. The HS-based method is built by training multiple two-layer multi-model forecasting framework (MMFF) models independently with the same-hour subsets. The final optimal model is a combination of MMFF models with the best-performed blending algorithm at every hour. At the forecasting stage, the most suitable model is selected to perform the forecasting subtask of a certain hour. The HS-based method is validated by 1-year data with six solar features collected by the National Renewable Energy Laboratory (NREL). Results show that the HS-based method outperforms the non-HS (all-in-one) method significantly with the same MMFF architecture, wherein the optimal HS- based method outperforms the best all-in-one method by 10.94% and 7.74% based on the normalized mean absolute error and normalized root mean square error, respectively.


title: The usability of the optical parametric amplification of light for high-angular-resolution imaging and fast astrometry

abstract: High-angular-resolution imaging is crucial for many applications in modern astronomy and astrophysics. The fundamental diffraction limit constrains the resolving power of both ground-based and spaceborne telescopes. The recent idea of a quantum telescope based on the optical parametric amplification (OPA) of light aims to bypass this limit for the imaging of extended sources by an order of magnitude or more. We present an updated scheme of an OPA-based device and a more accurate model of the signal amplification by such a device. The semiclassical model that we present predicts that the noise in such a system will form so-called light speckles as a result of light interference in the optical path. Based on this model, we analysed the efficiency of OPA in increasing the angular resolution of the imaging of extended targets and the precise localization of a distant point source. According to our new model, OPA offers a gain in resolved imaging in comparison to classical optics. For a given time-span, we found that OPA can be more efficient in localizing a single distant point source than classical telescopes.


title: A Gemini view of the galaxy cluster RXC J1504-0248: insights on the nature of the central gaseous filaments

abstract: We revisit the galaxy cluster RXC J1504-0248, a remarkable example of a structure with a strong cool core in a near redshift ($z = 0.216$). We performed a combined analysis using photometric and spectroscopic data obtained at Gemini South Telescope. We estimated the cluster mass through gravitational lensing, obtaining $M_{200} = 5.3\pm0.4 \times 10^{14}$ $h_{70}^{-1}$ M$_\odot$ within $R_{200} = 1.56 \pm 0.04$ $h^{-1}_{70}$ Mpc, in agreement with a virial mass estimate. This cluster presents a prominent filamentary structure associated to its BCG, located mainly along its major axis and aligned with the X-ray emission. A combined study of three emission line diagnostic diagrams has shown that the filament emission falls in the so-called transition region of these diagrams. Consequently, several ionizing sources should be playing an meaningful role. We have argued that old stars, often invoked to explain LINER emission, should not be the major source of ionization. We have noticed that most of the filamentary emission has line ratios consistent with the shock excitation limits obtained from shock models. We also found that line fluxes are related to gas velocities (here estimated from line widths) by power-laws with slopes in the range expected from shock models. These models also show, however, that only ~10% of H$\alpha$ luminosity can be explained by shocks. We conclude that shocks probably associated to the cooling of the intracluster gas in a filamentary structure may indeed be contributing to the filament nebular emission, but can not be the major source of ionizing photons.


title: Examples of cyclic polynomially bounded operators that are not similar to contractions, II

abstract: The question if polynomially bounded operator is similar to a contraction was posed by Halmos and was answered in the negative by Pisier. His counterexample is an operator of infinite multiplicity, while all its restrictions on invariant subspaces of finite multiplicity are similar to contractions. In [G16], cyclic polynomially bounded operators which are not similar to contractions was constructed. The construction was based on a perturbation of the sequence of finite dimensional operators which is uniformly polynomially bounded, but is not uniformly completely polynomially bounded, constructed by Pisier. In this paper, a cyclic polynomially bounded operator $T_0$ such that $T_0$ is not similar to a contraction and $\omega_a(T_0)=\mathbb O$, is constructed. Here $\omega_a(z)=\exp(a\frac{z+1}{z-1})$, $z\in\mathbb D$, $a>0$, and $\mathbb D$ is the open unit disk. To obtain such $T_0$, a slight modification of the construction from [G16] is needed.


title: A local depth measure for general data

abstract: We introduce the Integrated Dual Local Depth which is a local depth measure for data in a Banach space based on the use of one-dimensional projections. The properties of a depth measure are analyzed under this setting and a proper definition of local symmetry is given. Moreover, strong consistency results for the local depth and also for the local depth regions are attained. Finally, applications to descriptive data analysis and classification are analyzed, making the special focus on multivariate functional data, where we obtain very promising results.


title: The Electric Two-echelon Vehicle Routing Problem

abstract: Two-echelon distribution systems are attractive from an economical standpoint and help to keep large vehicles out of city centers. Large trucks can be used to deliver goods to intermediate facilities in accessible locations, whereas smaller vehicles allow to reach the final customers. Due to their reduced size and emissions, companies consider using an electric fleet of terrestrian or aerial vehicles for last mile deliveries. Route planning in multi-tier logistics leads to notoriously difficult problems. This difficulty is accrued in the presence of an electric fleet, since each vehicle operates on a smaller range, and may require visits to charging stations. To study these challenges, we introduce the Electric Two-echelon Vehicle Routing Problem as a prototypical problem. We propose a large neighbourhood search metaheuristic as well as an exact mathematical programming algorithm, which uses decomposition techniques to enumerate promising first-level solutions, in conjunction with bounding functions and route enumeration for the second-level routes. These algorithms produce optimal or near-optimal solutions for the problem, and allow us to evaluate the impact of several defining features of optimized battery-powered distribution networks. We created representative E2EVRP benchmark instances to simulate realistic metropolitan areas. In particular, we observe that the detour miles due to recharging decrease proportionally to $1/\rho^x$ with $x \approx 5/4$ as a function of the charging stations density $\rho$; e.g., in a scenario where the density of charging stations is doubled, recharging detours are reduced by 58\%. Finally, we evaluate the trade-off between battery capacity and detour miles. This estimate is critical for strategic fleet-acquisition decisions, in a context where large batteries are generally more costly and less environment-friendly.


title: Bounds on Dark Matter annihilations from 21 cm data

abstract: The observation of an absorption feature in the 21 cm spectrum at redshift $z\approx 17$ implies bounds on Dark Matter annihilations for a broad range of masses, given that significant heating of the intergalactic medium would have erased such feature. The resulting bounds on the DM annihilation cross sections are comparable to the strongest ones from all other observables.


title: Importance of non-local electron correlations in BaNiS$_{2}$ semimetal from quantum oscillations studies

abstract: By means of Shubnikov-de-Haas and de-Haas-van-Alphen oscillations, and ab initio calculations, we have studied the Fermi surface of high-quality BaNiS$_2$ single crystals, with mean free path $l \sim 400 ~\text{\AA}$. The angle and temperature dependence of quantum oscillations indicates a quasi-two-dimensional Fermi surface, made of an electron-like tube centred at $\Gamma$, and of 4 hole-like cones, generated by Dirac bands, weakly dispersive in the out-of-plane direction. Ab initio electronic structure calculations, in the density functional theory framework, show that the inclusion of screened exchange is necessary to account for the experimental Fermi pockets. Therefore, the choice of the functional becomes crucial. A modified HSE hybrid functional with 7% of exact exchange outperforms both GGA and GGA+U density functionals, signalling the importance of non-local screened-exchange interactions in BaNiS$_2$, and, more generally, in $3d$ compensated semimetals.


title: Mechanics of magma chamber with the implication of the effect of CO2 fluxing

abstract: As the magma ascends from its depth of generation to the surface, it is often stored in a series of chambers along the way. The rheological contrast between the viscous magma in the magmatic chambers and the surrounding rocks disturbed the stress field which can give rise to various modes of rock failure at magmatic pressures less than the lithostatic stress, leading to an eruption. Different modes of mechanical failure of the chamber walls are considered depending on the geometry and the sign the relative pressure. Relaxation of viscous stress around magmatic chambers, which is important on the time scale of weeks to months is considered in the analysis of stability with application to both large and extra-large magmatic chambers such as Yellowstone. The effects of a strong deep CO2 flux in Yellowstone are considered in detail. The analysis shows that variations in the flow rate around the observed mean value of 40 kg/m2/yr in the hydrothermally active areas can change the composition of the magma for several hundred thousand years, and cause periodic uplift and subsidence of the caldera surface with a period of several decades.


title: Measurable selector in Kadison's carpenter's theorem

abstract: We show the existence of a measurable selector in Carpenter's Theorem due to Kadison. This solves a problem posed by Jasper and the first author. As an application we obtain a characterization of all possible spectral functions of shift-invariant subspaces of $L^2(\mathbb R^d)$ and Carpenter's Theorem for type I$_\infty$ von Neumann algebras.


title: Holographic second laws of black hole thermodynamics

abstract: Recently, it has been shown that for out-of-equilibrium systems, there are additional constraints on thermodynamical evolution besides the ordinary second law. These form a new family of second laws of thermodynamics, which are equivalent to the monotonicity of quantum R\'enyi divergences. In black hole thermodynamics, the usual second law is manifest as the area increase theorem. Hence one may ask if these additional laws imply new restrictions for gravitational dynamics, such as for out-of-equilibrium black holes? Inspired by this question, we study these constraints within the AdS/CFT correspondence. First, we show that the R\'enyi divergence can be computed via a Euclidean path integral for a certain class of excited CFT states. Applying this construction to the boundary CFT, the R\'enyi divergence is evaluated as the renormalized action for a particular bulk solution of a minimally coupled gravity-scalar system. Further, within this framework, we show that there exist transitions which are allowed by the traditional second law, but forbidden by the additional thermodynamical constraints. We speculate on the implications of our findings.


title: Smooth orthogonal projections on Riemannian manifold

abstract: We construct a decomposition of the identity operator on a Riemannian manifold $M$ as a sum of smooth orthogonal projections subordinate to an open cover of $M$. This extends a decomposition of the real line by smooth orthogonal projection due to Coifman, Meyer and Auscher, Weiss, Wickerhauser, and a similar decomposition when $M$ is the sphere by the first two authors.


title: The Lottery Ticket Hypothesis: Finding Sparse, Trainable Neural Networks

abstract: Neural network pruning techniques can reduce the parameter counts of trained networks by over 90%, decreasing storage requirements and improving computational performance of inference without compromising accuracy. However, contemporary experience is that the sparse architectures produced by pruning are difficult to train from the start, which would similarly improve training performance. We find that a standard pruning technique naturally uncovers subnetworks whose initializations made them capable of training effectively. Based on these results, we articulate the "lottery ticket hypothesis:" dense, randomly-initialized, feed-forward networks contain subnetworks ("winning tickets") that - when trained in isolation - reach test accuracy comparable to the original network in a similar number of iterations. The winning tickets we find have won the initialization lottery: their connections have initial weights that make training particularly effective. We present an algorithm to identify winning tickets and a series of experiments that support the lottery ticket hypothesis and the importance of these fortuitous initializations. We consistently find winning tickets that are less than 10-20% of the size of several fully-connected and convolutional feed-forward architectures for MNIST and CIFAR10. Above this size, the winning tickets that we find learn faster than the original network and reach higher test accuracy.


title: Spin systems from loop soups

abstract: We study spin systems defined by the winding of a random walk loop soup. For a particular choice of loop soup intensity, we show that the corresponding spin system is reflection-positive and is dual, in the Kramers-Wannier sense, to the spin system $\text{sgn}(\varphi)$ where $\varphi$ is a discrete Gaussian free field. In general, we show that the spin correlation functions have conformally covariant scaling limits corresponding to the one-parameter family of functions studied by Camia, Gandolfi and Kleban (Nuclear Physics B, 902, 2016) and defined in terms of the winding of the Brownian loop soup. These functions have properties consistent with the behavior of correlation functions of conformal primaries in a conformal field theory. Here, we prove that they do correspond to correlation functions of continuum fields (random generalized functions) for values of the intensity of the Brownian loop soup that are not too large.


title: Restrictions of aspherical arrangements

abstract: In this note we present examples of $K(\pi,1)$-arrangements which admit a restriction which fails to be $K(\pi,1)$. This shows that asphericity is not hereditary among hyperplane arrangements.


title: On the Numerical Robustness of the Streaming Instability: Particle Concentration and Gas Dynamics in Protoplanetary Disks

abstract: The Streaming Instability (SI) is a mechanism to concentrate solids in protoplanetary disks. Nonlinear particle clumping from the SI can trigger gravitational collapse into planetesimals. To better understand the numerical robustness of the SI, we perform a suite of vertically-stratified 3D simulations with fixed physical parameters known to produce strong clumping. We vary the numerical implementation, namely the computational domain size and the vertical boundary conditions (vBCs), comparing newly-implemented outflow vBCs to the previously-used periodic and reflecting vBCs. We find strong particle clumping by the SI is mostly independent of the vBCs. However, peak particle densities are higher in larger simulation domains due to a larger particle mass reservoir. We report SI-triggered zonal flows, i.e., azimuthally-banded radial variations of gas pressure. These structures have low amplitudes, insufficient to halt particle radial drift, confirming that particle trapping in gas pressure maxima is not the mechanism of the SI. We find that outflow vBCs produce artificially large gas outflow rates at vertical boundaries. However, the outflow vBCs reduce artificial reflections at vertical boundaries, allowing more particle sedimentation, and showing less temporal variation and better convergence with box size. The radial spacing of dense particle filaments is $\sim0.15$ gas scale heights ($H$) for all vBCs, which sets the feeding zone for planetesimal growth in self-gravitating simulations. Our results validate the use of the outflow vBCs in SI simulations, even with vertical boundaries close ($\leq 0.4H$) to the disk midplane. Overall, our study demonstrates the numerical robustness of nonlinear particle clumping by the SI.


title: Precision and Recall for Time Series

abstract: Classical anomaly detection is principally concerned with point-based anomalies, those anomalies that occur at a single point in time. Yet, many real-world anomalies are range-based, meaning they occur over a period of time. Motivated by this observation, we present a new mathematical model to evaluate the accuracy of time series classification algorithms. Our model expands the well-known Precision and Recall metrics to measure ranges, while simultaneously enabling customization support for domain-specific preferences.


title: Signature Calculation of Area Hermitian Form on Some Spaces of Polygons

abstract: This chapter is motivated by the paper by Thurston on triangulations of the sphere and singular flat metrics on the sphere. Thurston locally parametrized the moduli space of singular flat metrics on the sphere with prescribed positive curvature data by the complex hyperbolic space of appropriate dimension. This work can be considered as a generalization of signature calculation of the Hermitian form that he made in his paper. The moduli space of singular flat metrics having unit area on the sphere with prescribed curvature data can be locally parametrized by certain spaces of polygons. This can be done by cutting singular flat spheres through length minimizing geodesics from a fixed singular point to the others. In that case the space of polygons is a complex vector space of dimension $n-1$ when there are $n+1$ singular points. Also there is natural area Hermitian form of signature $(1,n-2)$ on this vector space. In this chapter we calculate the signature of the area Hermitian form on some spaces of polygons which locally parametrize the moduli space of singular flat metrics having unit area on the sphere with one singular point of negative curvature. The formula we obtain depends only on the sum of the curvatures of the singular points having positive curvature. This paper will appear in the book "In the tradition of Thurston, Vol. II", Springer, 2022.


title: X-ray emission from black-hole and neutron-star binaries

abstract: In this chapter, I present the main X-ray observational characteristics of black-hole binaries and low magnetic field neutron-star binaries, concentrating on what can be considered similarities or differences, with particular emphasis on their fast-timing behaviour.


title: Deep Auxiliary Learning for Visual Localization and Odometry

abstract: Localization is an indispensable component of a robot's autonomy stack that enables it to determine where it is in the environment, essentially making it a precursor for any action execution or planning. Although convolutional neural networks have shown promising results for visual localization, they are still grossly outperformed by state-of-the-art local feature-based techniques. In this work, we propose VLocNet, a new convolutional neural network architecture for 6-DoF global pose regression and odometry estimation from consecutive monocular images. Our multitask model incorporates hard parameter sharing, thus being compact and enabling real-time inference, in addition to being end-to-end trainable. We propose a novel loss function that utilizes auxiliary learning to leverage relative pose information during training, thereby constraining the search space to obtain consistent pose estimates. We evaluate our proposed VLocNet on indoor as well as outdoor datasets and show that even our single task model exceeds the performance of state-of-the-art deep architectures for global localization, while achieving competitive performance for visual odometry estimation. Furthermore, we present extensive experimental evaluations utilizing our proposed Geometric Consistency Loss that show the effectiveness of multitask learning and demonstrate that our model is the first deep learning technique to be on par with, and in some cases outperforms state-of-the-art SIFT-based approaches.


title: Pollen Patterns Form from Modulated Phases

abstract: Pollen grains are known for their impressive variety of species-specific, microscale surface patterning. Despite having similar biological developmental steps, pollen grain surface features are remarkably geometrically varied. Previous work suggests that a physical process may drive this pattern formation and that the observed diversity of patterns can be explained by viewing pollen pattern development as a phase transition to a spatially modulated phase. Several studies have shown that the polysaccharide material of plant cell walls undergoes phase separation in the absence of cross-linking stabilizers of the mixed phase. Here we show experimental evidence of a change in density of the extracellular polysaccharide material (primexine) during pollen cell development leads to a spatially modulated phase. The spatial pattern of this phase-separated primexine is also mechanically coupled to the undulation of the pollen cell membrane. The resulting patterned pools of denser primexine form the negative template of the ultimate sites of sporopollenin deposition, leading to the final micropattern observed in the mature pollen. We then present a general physical model of pattern formation via modulated phases. Using analytical and numerical techniques, we find that most of the pollen micropatterns observed in biological evolution could result from a physical process of modulated phases. However, an analysis of the relative rates of transitions from states that are equilibrated to or from states that are not equilibrated suggests that while equilibrium states of this process have occurred throughout evolutionary history, there has been no particular evolutionary selection for distinctly patterned equilibrated states.


title: Converting non-relativistic dark matter to radiation

abstract: Dark matter in the cosmological concordance model is parameterised by a single number, describing the covariantly conserved energy density of a non-relativistic fluid. Here we test this assumption in a model-independent and conservative way by considering the possibility that, at any point during the cosmological evolution, dark matter may be converted into a non-interacting form of radiation. This scenario encompasses, but is more general than, the cases where dark matter decays or annihilates into these states. We show that observations of the cosmic microwave background allow to strongly constrain this scenario for any conversion time after big bang nucleosynthesis. We discuss in detail, both from a Bayesian and frequentist point of view, in which sense adding large-scale structure observations may even provide a certain preference for a conversion of dark matter to radiation at late times. Finally we apply our general results to a specific particle physics realisation of such a scenario, featuring late kinetic decoupling and Sommerfeld-enhanced dark matter annihilation. We identify a small part of parameter space that both mitigates the tension between cosmic microwave and large-scale structure data and allows for velocity-dependent dark matter self-interactions strong enough to address the small-scale problems of structure formation.


title: The groomed and ungroomed jet mass distribution for inclusive jet production at the LHC

abstract: We study jet mass distributions measured in the single inclusive jet production in proton-proton collisions $pp\to \text{jet}+X$ at the LHC. We consider both standard ungroomed jets as well as soft drop groomed jets. Within the Soft Collinear Effective Theory (SCET), we establish QCD factorization theorems for both cases and we study their relation. The developed framework allows for the joint resummation of several classes of logarithmic corrections to all orders in the strong coupling constant. For the ungroomed case, we resum logarithms in the jet radius parameter and in the small jet mass. For the groomed case, we resum in addition the logarithms in the soft threshold parameter $z_{\text{cut}}$ which is introduced by the soft drop grooming algorithm. In this way, we are able to reliably determine the absolute normalization of the groomed jet mass distribution in proton-proton collisions. All logarithmic corrections are resummed to the next-to-leading logarithmic accuracy. We present numerical results and compare with the available data from the LHC. For both the groomed and ungroomed jet mass distributions we find very good agreement after including non-perturbative corrections.


title: Aspherical Supernovae: Effects on Early Light Curves

abstract: Early light from core-collapse supernovae, now detectable in high-cadence surveys, holds clues to a star and its environment just before it explodes. However, effects that alter the early light have not been fully explored. We highlight the possibility of non-radial flows at the time of shock breakout. These develop in sufficiently non-spherical explosions if the progenitor is not too diffuse. When they do develop, non-radial flows limit ejecta speeds and cause ejecta-ejecta collisions. We explore these phenomena and their observational implications, using global, axisymmetric, non-relativistic FLASH simulations of simplified polytropic progenitors, which we scale to representative stars. We develop a method to track photon production within the ejecta, enabling us to estimate band-dependent light curves from adiabatic simulations. Immediate breakout emission becomes hidden as an oblique flow develops. Non-spherical effects lead the shock-heated ejecta to release a more constant luminosity at a higher, evolving color temperature at early times, effectively mixing breakout light with the early light curve. Collisions between non-radial ejecta thermalize a small fraction of the explosion energy; we address emission from these collisions in a subsequent paper.


title: The Hyperbolic Higgs

abstract: We introduce the Hyperbolic Higgs, a novel solution to the little hierarchy problem that features Standard Model neutral scalar top partners. At one-loop order, the protection from ultraviolet sensitivity is due to an accidental non-compact symmetry of the Higgs potential that emerges in the infrared. Once the general features of the effective description are detailed, a completion that relies on a five dimensional supersymmetric framework is provided. Novel phenomenology is compared and contrasted with the Twin Higgs scenario.


title: Inner Super-Earths, Outer Gas Giants: How Pebble Isolation and Migration Feedback Keep Jupiters Cold

abstract: The majority of gas giants (planets of masses $\gtrsim10^2 M_\oplus$) are found to reside at distances beyond $\sim1$ au from their host stars. Within 1 au, the planetary population is dominated by super-Earths of $2-20 M_\oplus$. We show that this dichotomy between inner super-Earths and outer gas giants can be naturally explained should they form in nearly inviscid disks. In laminar disks, a planet can more easily repel disk gas away from its orbit. The feedback torque from the pile-up of gas inside the planet's orbit slows down and eventually halts migration. A pressure bump outside the planet's orbit traps pebbles and solids, starving the core. Gas giants are born cold and stay cold: more massive cores are preferentially formed at larger distances, and they barely migrate under disk feedback. We demonstrate this using 2D hydrodynamical simulations of disk-planet interaction lasting up to $10^5$ years: we track planet migration and pebble accretion until both come to an end by disk feedback. Whether cores undergo runaway gas accretion to become gas giants or not is determined by computing 1D gas accretion models. Our simulations show that in an inviscid minimum mass solar nebula, gas giants do not form inside $\sim$0.5 au, nor can they migrate there while the disk is present. We also explore the dependence on disk mass, and find that gas giants form further out in less massive disks.


title: Unsqueezing of standing waves due to inflationary domain structure

abstract: The so-called "trans-Planckian" problem of inflation may be evaded by positing that modes come into existence only when they became "cis-Planckian" by virtue of expansion. However, this would imply that for any mode a new random realization would have to be drawn every $N$ wavelengths, with $N$ typically of order 1000 (but it could be larger or smaller). Such a re-drawing of realizations leads to a heteroskodastic distribution if the region under observation contains several such independent domains. This has no effect on the sampled power spectrum for a scale-invariant raw spectrum, but at very small scales it leads to a spectral index bias towards scale-invariance and smooths oscillations in the spectrum. The domain structure would also "unsqueeze" some of the propagating waves, i.e., dismantle their standing wave character. By describing standing waves as travelling waves of the same amplitude moving in opposite directions we determine the observational effects of unsqueezing. We find that it would erase the Doppler peaks in the CMB, but only on very small angular scales, where the primordial signal may not be readily accessible. The standing waves in a primordial gravitational wave background would also be turned into travelling waves. This unsqueezing of the gravitational wave background may constitute a detectable phenomenon.


title: Neutrino Discovery Limit of Dark Matter Direct Detection Experiments in the Presence of Non-Standard Interactions

abstract: The detection of coherent neutrino-nucleus scattering by the COHERENT collaboration has set on quantitative grounds the existence of an irreducible neutrino background in direct detection searches of Weakly Interacting Massive Dark Matter candidates. This background leads to an ultimate discovery limit for these experiments: a minimum Dark Matter interaction cross section below which events produced by the coherent neutrino scattering will mimic the Dark Matter signal, the so-called \emph{neutrino floor}. In this work we study the modification of such neutrino floor induced by non-standard neutrino interactions within their presently allowed values by the global analysis of oscillation and COHERENT data. By using the full likelihood information of such global analysis we consistently account for the correlated effects of non-standard neutrino interactions both in the neutrino propagation in matter and in its interaction in the detector. We quantify their impact on the neutrino floor for five future experiments: DARWIN (Xe), ARGO (Ar), Super-CDMS HV (Ge and Si) and CRESST phase III (CaWO$_4$). Quantitatively, we find that non-standard neutrino interactions allowed at the $3\sigma$ level can result in an increase of the neutrino floor of up to a factor $\sim 5$ with respect to the Standard Model expectations and impact the expected sensitivities of the ARGO, CRESST phase III and DARWIN experiments.


title: Singlet Scalar Top Partners from Accidental Supersymmetry

abstract: We present a model wherein the Higgs mass is protected from the quadratic one-loop top quark corrections by scalar particles that are complete singlets under the Standard Model (SM) gauge group. While bearing some similarity to Folded Supersymmetry, the construction is purely four dimensional and enjoys more parametric freedom, allowing electroweak symmetry breaking to occur easily. The cancelation of the top loop quadratic divergence is ensured by a $Z_3$ symmetry that relates the SM top sector and two hidden top sectors, each charged under its own hidden color group. In addition to the singlet scalars, the hidden sectors contain electroweak-charged supermultiplets below the TeV scale, which provide the main access to this model at colliders. The phenomenology presents both differences and similarities with respect to other realizations of neutral naturalness. Generally, the glueballs of hidden color have longer decay lengths. The production of hidden sector particles results in quirk or squirk bound states, which later annihilate. We survey the possible signatures and corresponding experimental constraints.


title: SPH to Grid: a new integral conserving method

abstract: Analysing data from Smoothed Particle Hydrodynamics (SPH) simulations is about understanding global fluid properties rather than individual fluid elements. Therefore, in order to properly understand the outcome of such simulations it is crucial to transition from a particle to a grid based picture. In this paper we briefly summarise different methods of calculating a representative volume discretisation from SPH data and propose an improved version of commonly used techniques. We present a possibility to generate accurate 2D data directly without the CPU time and memory consuming detour over a 3D grid. We lay out the importance of an accurate algorithm to conserve integral fluid properties and to properly treat small scale structures using a typical galaxy simulation snapshot. For demonstration purposes we additionally calculate velocity power spectra and as expected find the main differences on small scales. Finally we propose two new multi-purpose analysis packages which utilise the new algorithms: Pygad and SPHMapper.


title: Spectroscopic decomposition of NGC 3521: unveiling the properties of the bulge and disc

abstract: We study the kinematics and the stellar populations of the bulge and disc of the spiral galaxy NGC 3521. At each position in the field of view, we separate the contributions of the bulge and the disc from the total observed spectrum and study their kinematics, age, and metallicities independently. Their properties are clearly distinct: the bulge rotates more slowly, has a higher velocity dispersion, and is less luminous than the disc. We identify three main populations of stars in NGC 3521: old ($\geq7$ Gyr), intermediate ($\approx$ 3 Gyr), and young ($\leq$1 Gyr). The mass and light of NGC 3521 are dominated by the intermediate stellar population. The youngest population contributes mostly to the disc component and its contribution increases with radius. We also study the luminosity-weighed properties of the stars in NGC 3521. Along the photometric major axis, we find: i) no age gradient for the stars in the bulge, and a negative age gradient for the stars in the disc; ii) negative metallicity gradients and sub-solar $\alpha$-enhancement for both the bulge and the disc. We propose the following picture for the formation of NGC 3521: initial formation a long time ago ($\geq 7$ Gyr), followed by a second burst of star formation or a merger ($\approx$ 3 Gyrs ago), which contributed predominantly to the mass-build up of the bulge. Recently ($\leq 1$ Gyr), the disc of NGC 3521 experienced an additional episode of star formation that started in the innermost regions.


title: Directing Chemotaxis-Based Spatial Self-Organization via Biased, Random Initial Conditions

abstract: Inspired by the chemotaxis interaction of living cells, we have developed an agent-based approach for self-organizing shape formation. Since all our simulations begin with a different uniform random configuration and our agents move stochastically, it has been observed that the self-organization process may form two or more stable final configurations. These differing configurations may be characterized via statistical moments of the agents' locations. In order to direct the agents to robustly form one specific configuration, we generate biased initial conditions whose statistical moments are related to moments of the desired configuration. With this approach, we are able to successfully direct the aggregating swarms to produced a desired macroscopic shape, starting from randomized initial conditions with controlled statistical properties.


title: Accurately predicting the escape fraction of ionizing photons using restframe ultraviolet absorption lines

abstract: The fraction of ionizing photons that escape high-redshift galaxies sensitively determines whether galaxies reionized the early universe. However, this escape fraction cannot be measured from high-redshift galaxies because the opacity of the intergalactic medium is large at high redshifts. Without methods to indirectly measure the escape fraction of high-redshift galaxies, it is unlikely that we will know what reionized the universe. Here, we analyze the far-ultraviolet (UV) H I (Lyman series) and low-ionization metal absorption lines of nine low-redshift, confirmed Lyman continuum emitting galaxies. We use the H I covering fractions, column densities, and dust attenuations measured in a companion paper to predict the escape fraction of ionizing photons. We find good agreement between the predicted and observed Lyman continuum escape fractions (within $1.4\sigma$) using both the H I and ISM absorption lines. The ionizing photons escape through holes in the H I, but we show that dust attenuation reduces the fraction of photons that escape galaxies. This means that the average high-redshift galaxy likely emits more ionizing photons than low-redshift galaxies. Two other indirect methods accurately predict the escape fractions: the Ly$\alpha$ escape fraction and the optical [O III]/[O II] flux ratio. We use these indirect methods to predict the escape fraction of a sample of 21 galaxies with rest-frame UV spectra but without Lyman continuum observations. Many of these galaxies have low escape fractions ($f_{\rm esc} \le 1$\%), but 11 have escape fractions $>1$\%. The methods presented here will measure the escape fractions of high-redshift galaxies, enabling future telescopes to determine whether star-forming galaxies reionized the early universe.


title: Nucleon Quark Distribution Functions from the Dyson-Schwinger Equations

abstract: We present results for the nucleon's leading-twist spin-independent valence parton distribution functions obtained from a theoretical framework based on the Dyson-Schwinger equations (DSEs) of QCD that previously gave an excellent description of nucleon electromagnetic form factors. We employ the rainbow-ladder truncation of the DSEs and utilize nucleon bound state amplitudes from the Poincar\'e-covariant Faddeev equation, where the dominant scalar and axial-vector quark-quark correlations are included. This DSE framework is used to numerically evaluate the first 20 moments of the valence $u$ and $d$ quark distribution functions, from which the $x$-dependence of the distributions is found to be well constrained. We find good agreement with empirical parameterizations of experimental data and make the prediction that the $d/u$ ratio in the $x\to 1$ limit, invariant under scale evolution, takes the value $d/u \to 0.087 \pm 0.010$. We find that this ratio is rather sensitive to the strength of axial-vector diquark correlations. However, contrary to a naive expectation, our result for the $d/u$ ratio in the $x\to 1$ limit does not vanish when only scalar diquark correlations are present, although it is an order of magnitude smaller than our $d/u$ result that also includes axial-vector diquarks. The valence quark distribution results are set in a broader context via a simple pion cloud model estimate of sea-quark light-cone momenta and gluon light-cone momentum.


title: Quantum simulation of partially distinguishable boson sampling

abstract: Boson Sampling is the problem of sampling from the same output probability distribution as a collection of indistinguishable single photons input into a linear interferometer. It has been shown that, subject to certain computational complexity conjectures, in general the problem is difficult to solve classically, motivating optical experiments aimed at demonstrating quantum computational "supremacy". There are a number of challenges faced by such experiments, including the generation of indistinguishable single photons. We provide a quantum circuit that simulates bosonic sampling with arbitrarily distinguishable particles. This makes clear how distinguishabililty leads to decoherence in the standard quantum circuit model, allowing insight to be gained. At the heart of the circuit is the quantum Schur transform, which follows from a representation theoretic approach to the physics of distinguishable particles in first quantisation. The techniques are quite general and have application beyond boson sampling.


title: NNLO QCD corrections to the polarized top quark decay $t(\uparrow) \to X_b+W^+$

abstract: We compute the next-to-next-to-leading order (NNLO) QCD corrections to the decay $t(\uparrow) \to X_b +W^+$ of a polarized top quark. The spin-momentum correlation in this quasi two-body decay is described by the polar angle distribution $\mathrm{d}\Gamma/\mathrm{d}\cos\theta_P=\frac{\Gamma}{2}(1+P_t\, \alpha_P\, \cos\theta_P)$ where $P_t$ is the polarization of the top quark and $\alpha_P$ denotes the asymmetry parameter of the decay. For the latter we find $\alpha^{\mathrm{NNLO}}_P=0.3792\pm 0.0037$.


title: Listing Maximal Subgraphs in Strongly Accessible Set Systems

abstract: Algorithms for listing the subgraphs satisfying a given property (e.g.,being a clique, a cut, a cycle, etc.) fall within the general framework of set systems. A set system (U, F) uses a ground set U (e.g., the network nodes) and an indicator F, subset of 2^U, of which subsets of U have the required property. For the problem of listing all sets in F maximal under inclusion, the ambitious goal is to cover a large class of set systems, preserving at the same time the efficiency of the enumeration. Among the existing algorithms, the best-known ones list the maximal subsets in time proportional to their number but may require exponential space. In this paper we improve the state of the art in two directions by introducing an algorithmic framework that, under standard suitable conditions, simultaneously (i) extends the class of problems that can be solved efficiently to strongly accessible set systems, and (ii) reduces the additional space usage from exponential in |U| to stateless, thus accounting for just O(q) space, where q <= |U| is the largest size of a maximal set in F


title: A note on Perelman's no shrinking breather theorem

abstract: As an application of his entropy formula, Perelman proved that every compact shrinking breather is a shrinking gradient Ricci soliton. We give a proof for the complete noncompact case by using Perelman's $\mathcal{L}$-geometry. Our proof follows the argument in Lu and Zheng of constructing an ancient solution, and removes a technical assumption made by them.


title: Competition between delta isobars and hyperons and properties of compact stars

abstract: The $\Delta$-isobar degrees of freedom are included in the covariant density functional (CDF) theory to study the equation of state (EoS) and composition of dense matter in compact stars. In addition to $\Delta$'s we include the full octet of baryons, which allows us to study the interplay between the onset of delta isobars and hyperonic degrees of freedom. Using both the Hartree and Hartree-Fock approximation we find that $\Delta$'s appear already at densities slightly above the saturation density of nuclear matter for a wide range of the meson-$\Delta$ coupling constants. This delays the appearance of hyperons and significantly affects the gross properties of compact stars. Specifically, $\Delta$'s soften the EoS at low densities but stiffen it at high densities. This softening reduces the radius of a canonical $1.4 M_\odot$ star by up to 2~km for a reasonably attractive $\Delta$ potential in matter, while the stiffening results in larger maximum masses of compact stars. We conclude that the hypernuclear CDF parametrizations that satisfy the 2$M_\odot$ maximum mass constraint remain valid when $\Delta$ isobars are included, with the important consequence that the resulting stellar radii are shifted toward lower values, which is in agreement with the analysis of neutron star radii.


title: Hate Speech Detection: A Solved Problem? The Challenging Case of Long Tail on Twitter

abstract: In recent years, the increasing propagation of hate speech on social media and the urgent need for effective counter-measures have drawn significant investment from governments, companies, and researchers. A large number of methods have been developed for automated hate speech detection online. This aims to classify textual content into non-hate or hate speech, in which case the method may also identify the targeting characteristics (i.e., types of hate, such as race, and religion) in the hate speech. However, we notice significant difference between the performance of the two (i.e., non-hate v.s. hate). In this work, we argue for a focus on the latter problem for practical reasons. We show that it is a much more challenging task, as our analysis of the language in the typical datasets shows that hate speech lacks unique, discriminative features and therefore is found in the 'long tail' in a dataset that is difficult to discover. We then propose Deep Neural Network structures serving as feature extractors that are particularly effective for capturing the semantics of hate speech. Our methods are evaluated on the largest collection of hate speech datasets based on Twitter, and are shown to be able to outperform the best performing method by up to 5 percentage points in macro-average F1, or 8 percentage points in the more challenging case of identifying hateful content.


title: Disconnected Cuts in Claw-free Graphs

abstract: A disconnected cut of a connected graph is a vertex cut that itself also induces a disconnected subgraph. The decision problem whether a graph has a disconnected cut is called Disconnected Cut. This problem is closely related to several homomorphism and contraction problems, and fits in an extensive line of research on vertex cuts with additional properties. It is known that Disconnected Cut is NP-hard on general graphs, while polynomial-time algorithms are known for several graph classes. However, the complexity of the problem on claw-free graphs remained an open question. Its connection to the complexity of the problem to contract a claw-free graph to the 4-vertex cycle $C_4$ led Ito et al. (TCS 2011) to explicitly ask to resolve this open question. We prove that Disconnected Cut is polynomial-time solvable on claw-free graphs, answering the question of Ito et al. The centerpiece of our result is a novel decomposition theorem for claw-free graphs of diameter 2, which we believe is of independent interest and expands the research line initiated by Chudnovsky and Seymour (JCTB 2007-2012) and Hermelin et al. (ICALP 2011). On our way to exploit this decomposition theorem, we characterize how disconnected cuts interact with certain cobipartite subgraphs, and prove two further novel algorithmic results, namely Disconnected Cut is polynomial-time solvable on circular-arc graphs and line graphs.


title: Automating Reading Comprehension by Generating Question and Answer Pairs

abstract: Neural network-based methods represent the state-of-the-art in question generation from text. Existing work focuses on generating only questions from text without concerning itself with answer generation. Moreover, our analysis shows that handling rare words and generating the most appropriate question given a candidate answer are still challenges facing existing approaches. We present a novel two-stage process to generate question-answer pairs from the text. For the first stage, we present alternatives for encoding the span of the pivotal answer in the sentence using Pointer Networks. In our second stage, we employ sequence to sequence models for question generation, enhanced with rich linguistic features. Finally, global attention and answer encoding are used for generating the question most relevant to the answer. We motivate and linguistically analyze the role of each component in our framework and consider compositions of these. This analysis is supported by extensive experimental evaluations. Using standard evaluation metrics as well as human evaluations, our experimental results validate the significant improvement in the quality of questions generated by our framework over the state-of-the-art. The technique presented here represents another step towards more automated reading comprehension assessment. We also present a live system \footnote{Demo of the system is available at \url{https://www.cse.iitb.ac.in/~vishwajeet/autoqg.html}.} to demonstrate the effectiveness of our approach.


title: Syntax-Aware Language Modeling with Recurrent Neural Networks

abstract: Neural language models (LMs) are typically trained using only lexical features, such as surface forms of words. In this paper, we argue this deprives the LM of crucial syntactic signals that can be detected at high confidence using existing parsers. We present a simple but highly effective approach for training neural LMs using both lexical and syntactic information, and a novel approach for applying such LMs to unparsed text using sequential Monte Carlo sampling. In experiments on a range of corpora and corpus sizes, we show our approach consistently outperforms standard lexical LMs in character-level language modeling; on the other hand, for word-level models the models are on a par with standard language models. These results indicate potential for expanding LMs beyond lexical surface features to higher-level NLP features for character-level models.


title: Standing Wave Decomposition Gaussian Process

abstract: We propose a Standing Wave Decomposition (SWD) approximation to Gaussian Process regression (GP). GP involves a costly matrix inversion operation, which limits applicability to large data analysis. For an input space that can be approximated by a grid and when correlations among data are short-ranged, the kernel matrix inversion can be replaced by analytic diagonalization using the SWD. We show that this approach applies to uni- and multi-dimensional input data, extends to include longer-range correlations, and the grid can be in a latent space and used as inducing points. Through simulations, we show that our approximate method applied to the squared exponential kernel outperforms existing methods in predictive accuracy per unit time in the regime where data are plentiful. Our SWD-GP is recommended for regression analyses where there is a relatively large amount of data and/or there are constraints on computation time.


title: Co-occurrence of the Benford-like and Zipf Laws Arising from the Texts Representing Human and Artificial Languages

abstract: We demonstrate that large texts, representing human (English, Russian, Ukrainian) and artificial (C++, Java) languages, display quantitative patterns characterized by the Benford-like and Zipf laws. The frequency of a word following the Zipf law is inversely proportional to its rank, whereas the total numbers of a certain word appearing in the text generate the uneven Benford-like distribution of leading numbers. Excluding the most popular words essentially improves the correlation of actual textual data with the Zipfian distribution, whereas the Benford distribution of leading numbers (arising from the overall amount of a certain word) is insensitive to the same elimination procedure. The calculated values of the moduli of slopes of double logarithmical plots for artificial languages (C++, Java) are markedly larger than those for human ones.


title: Coupling tunable D-band directional coupler for millimeter-wave applications

abstract: A coupling tunable D-band directional coupler is designed based on a novel coupling grid structure proposed in this letter. The designed directional coupler has excellent performance with ultra-wideband. The coupling can be tuned from -28.2 dB to -33.2 dB at 140 GHz by changing the angle of the coupling grid, and the dynamic range of the coupling is about 5 dB. The return loss is smaller than -15 dB in the whole D-band from 110 GHz to 170 GHz. A 3-dB coupler use the similar coupling structure is also designed. The coupling is 3.3144 dB at the center frequency of 140 GHz.


title: Provably robust estimation of modulo 1 samples of a smooth function with applications to phase unwrapping

abstract: Consider an unknown smooth function $f: [0,1]^d \rightarrow \mathbb{R}$, and say we are given $n$ noisy mod 1 samples of $f$, i.e., $y_i = (f(x_i) + \eta_i)\mod 1$, for $x_i \in [0,1]^d$, where $\eta_i$ denotes the noise. Given the samples $(x_i,y_i)_{i=1}^{n}$, our goal is to recover smooth, robust estimates of the clean samples $f(x_i) \bmod 1$. We formulate a natural approach for solving this problem, which works with angular embeddings of the noisy mod 1 samples over the unit circle, inspired by the angular synchronization framework. This amounts to solving a smoothness regularized least-squares problem -- a quadratically constrained quadratic program (QCQP) -- where the variables are constrained to lie on the unit circle. Our approach is based on solving its relaxation, which is a trust-region sub-problem and hence solvable efficiently. We provide theoretical guarantees demonstrating its robustness to noise for adversarial, and random Gaussian and Bernoulli noise models. To the best of our knowledge, these are the first such theoretical results for this problem. We demonstrate the robustness and efficiency of our approach via extensive numerical simulations on synthetic data, along with a simple least-squares solution for the unwrapping stage, that recovers the original samples of $f$ (up to a global shift). It is shown to perform well at high levels of noise, when taking as input the denoised modulo $1$ samples. Finally, we also consider two other approaches for denoising the modulo 1 samples that leverage tools from Riemannian optimization on manifolds, including a Burer-Monteiro approach for a semidefinite programming relaxation of our formulation. For the two-dimensional version of the problem, which has applications in radar interferometry, we are able to solve instances of real-world data with a million sample points in under 10 seconds, on a personal laptop.


title: IcoRating: A Deep-Learning System for Scam ICO Identification

abstract: Cryptocurrencies (or digital tokens, digital currencies, e.g., BTC, ETH, XRP, NEO) have been rapidly gaining ground in use, value, and understanding among the public, bringing astonishing profits to investors. Unlike other money and banking systems, most digital tokens do not require central authorities. Being decentralized poses significant challenges for credit rating. Most ICOs are currently not subject to government regulations, which makes a reliable credit rating system for ICO projects necessary and urgent. In this paper, we introduce IcoRating, the first learning--based cryptocurrency rating system. We exploit natural-language processing techniques to analyze various aspects of 2,251 digital currencies to date, such as white paper content, founding teams, Github repositories, websites, etc. Supervised learning models are used to correlate the life span and the price change of cryptocurrencies with these features. For the best setting, the proposed system is able to identify scam ICO projects with 0.83 precision. We hope this work will help investors identify scam ICOs and attract more efforts in automatically evaluating and analyzing ICO projects.


title: Unique dynamic crossover in supercooled x,3-dihydroxypropyl acrylate (x = 1, 2) isomers mixture

abstract: The previtreous dynamics in glass forming monomer, glycerol monoacrylate (GMA), using broadband dielectric spectroscopy (BDS) was tested. Measurements revealed the clear dynamic crossover at temperature $T_B = 254$ K and the time scale $\tau(T_B) = 5.4$ ns for the primary (structural) relaxation time and no hallmarks for the crossover for the DC electric conductivity $\sigma_{DC}$. This result was revealed via the derivative-based and distortions-sensitive analysis $dln{H_{a}}/d(1/T)$ vs. $1/T$, where $H_a$ is for the apparent activation energy. Subsequent tests of the fractional Debye-Stokes-Einsten relation $\sigma_{DC}(\tau_{\alpha})^S = const$ showed that the crossover is associated with $S = 1$ (for $T>T_B$)->$S = 0.84$ (for $T<T_B$). The crossover is associated with the emergence of the secondary beta relaxation which smoothly develops deeply into the solid amorphous phase below the glass temperature $T_g$.


title: Competitive Machine Learning: Best Theoretical Prediction vs Optimization

abstract: Machine learning is often used in competitive scenarios: Participants learn and fit static models, and those models compete in a shared platform. The common assumption is that in order to win a competition one has to have the best predictive model, i.e., the model with the smallest out-sample error. Is that necessarily true? Does the best theoretical predictive model for a target always yield the best reward in a competition? If not, can one take the best model and purposefully change it into a theoretically inferior model which in practice results in a higher competitive edge? How does that modification look like? And finally, if all participants modify their prediction models towards the best practical performance, who benefits the most? players with inferior models, or those with theoretical superiority? The main theme of this paper is to raise these important questions and propose a theoretical model to answer them. We consider a study case where two linear predictive models compete over a shared target. The model with the closest estimate gets the whole reward, which is equal to the absolute value of the target. We characterize the reward function of each model, and using a basic game theoretic approach, demonstrate that the inferior competitor can significantly improve his performance by choosing optimal model coefficients that are different from the best theoretical prediction. This is a preliminary study that emphasizes the fact that in many applications where predictive machine learning is at the service of competition, much can be gained from practical (back-testing) optimization of the model compared to static prediction improvement.


title: Fault Localization Models in Debugging

abstract: Debugging is considered as a rigorous but important feature of software engineering process. Since more than a decade, the software engineering research community is exploring different techniques for removal of faults from programs but it is quite difficult to overcome all the faults of software programs. Thus, it is still remains as a real challenge for software debugging and maintenance community. In this paper, we briefly introduced software anomalies and faults classification and then explained different fault localization models using theory of diagnosis. Furthermore, we compared and contrasted between value based and dependencies based models in accordance with different real misbehaviours and presented some insight information for the debugging process. Moreover, we discussed the results of both models and manifested the shortcomings as well as advantages of these models in terms of debugging and maintenance.


title: Sequential Outlier Detection based on Incremental Decision Trees

abstract: We introduce an online outlier detection algorithm to detect outliers in a sequentially observed data stream. For this purpose, we use a two-stage filtering and hedging approach. In the first stage, we construct a multi-modal probability density function to model the normal samples. In the second stage, given a new observation, we label it as an anomaly if the value of aforementioned density function is below a specified threshold at the newly observed point. In order to construct our multi-modal density function, we use an incremental decision tree to construct a set of subspaces of the observation space. We train a single component density function of the exponential family using the observations, which fall inside each subspace represented on the tree. These single component density functions are then adaptively combined to produce our multi-modal density function, which is shown to achieve the performance of the best convex combination of the density functions defined on the subspaces. As we observe more samples, our tree grows and produces more subspaces. As a result, our modeling power increases in time, while mitigating overfitting issues. In order to choose our threshold level to label the observations, we use an adaptive thresholding scheme. We show that our adaptive threshold level achieves the performance of the optimal pre-fixed threshold level, which knows the observation labels in hindsight. Our algorithm provides significant performance improvements over the state of the art in our wide set of experiments involving both synthetic as well as real data.


title: Active Galactic Nuclei Feedback in an Elliptical Galaxy with the Most Updated AGN Physics (II): High-Angular Momentum Case

abstract: This is the second paper of our series of works of studying the effects of active galactic nuclei (AGN) feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. In these simulations, the inner boundary is chosen so that the Bondi radius is resolved. Physical processes like star formation, SNe Ia and II are taken into account. Compared to previous works, the main improvements is that we adopt the most updated AGN physics, which is described in detail in the first paper of this series (Yuan et al. 2018, Paper I). These improvements include the discrimination of the two accretion modes of the central AGN and the most updated descriptions of the wind and radiation in the two modes. In Paper I, we consider the case that the specific angular momentum of the gas in the galaxy is very low. In this paper, we consider the case that the specific angular momentum of the gas is high. In the galactic scale, we adopt the gravitational torques raised due to non-axisymmetric structure in the galaxy as the mechanism of the transfer of angular momentum of gas, as proposed in some recent works. Since our simulations are axisymmetric, we make use of a parameterized prescription to mimic this mechanism. Same as Paper I, we investigate the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. Special attention is paid to the effects of specific angular momentum of the galaxy on these properties. We find that some results are qualitatively similar to those shown in Paper I, while some results such as star formation and black hole growth do show a significant difference due to the mass concentration in the galactic disk as a consequence of galactic rotation.


title: Existence and Construction of Galilean invariant $z\neq2$ Theories

abstract: We prove a no-go theorem for the construction of a Galilean boost invariant and $z\neq2$ anisotropic scale invariant field theory with a finite dimensional basis of fields. Two point correlators in such theories, we show, grow unboundedly with spatial separation. Correlators of theories with an infinite dimensional basis of fields, for example, labeled by a continuous parameter, do not necessarily exhibit this bad behavior. Hence, such theories behave effectively as if in one extra dimension. Embedding the symmetry algebra into the conformal algebra of one higher dimension also reveals the existence of an internal continuous parameter. Consideration of isometries shows that the non-relativistic holographic picture assumes a canonical form, where the bulk gravitational theory lives in a space-time with one extra dimension. This can be contrasted with the original proposal by Balasubramanian and McGreevy, and by Son, where the metric of a $d+2$ dimensional space-time is proposed to be dual of a $d$ dimensional field theory. We provide explicit examples of theories living at fixed point with anisotropic scaling exponent $z=\frac{2\ell}{\ell+1}\,,\ell\in \mathbb{Z}$

categories: hep-th cond-mat.quant-gas cond-mat.str-el math-ph math.MP


title: Nonparametric Risk Assessment and Density Estimation for Persistence Landscapes

abstract: This paper presents approximate confidence intervals for each function of parameters in a Banach space based on a bootstrap algorithm. We apply kernel density approach to estimate the persistence landscape. In addition, we evaluate the quality distribution function estimator of random variables using integrated mean square error (IMSE). The results of simulation studies show a significant improvement achieved by our approach compared to the standard version of confidence intervals algorithm. In the next step, we provide several algorithms to solve our model. Finally, real data analysis shows that the accuracy of our method compared to that of previous works for computing the confidence interval.


title: Early-type galaxies in the Antlia Cluster: Catalogue and isophotal analysis

abstract: We present a statistical isophotal analysis of 138 early-type galaxies in the Antlia cluster, located at a distance of ${\sim} 35$ Mpc. The observational material consists of CCD images of four $36$ arcmin ${\times} 36$ arcmin fields obtained with the MOSAIC II camera at the Blanco 4-m telescope at CTIO. Our present work supersedes previous Antlia studies in the sense that the covered area is four times larger, the limiting magnitude is $M_{B} {\sim} -9.6$ mag, and the surface photometry parameters of each galaxy are derived from S\'ersic model fits extrapolated to infinity. In a companion previous study we focused on the scaling relations obtained by means of surface photometry, and now we present the data, on which the previous paper is based, the parameters of the isophotal fits as well as an isophotal analysis. For each galaxy, we derive isophotal shape parameters along the semi-major axis and search for correlations within different radial bins. Through extensive statistical tests, we also analyse the behaviour of these values against photometric and global parameters of the galaxies themselves. While some galaxies do display radial gradients in their ellipticity (${\epsilon}$) and/or their Fourier coefficients, differences in mean values between adjacent regions are not statistically significant. Regarding Fourier coefficients, dwarf galaxies usually display gradients between all adjacent regions, while non-dwarfs tend to show this behaviour just between the two outermost regions. Globally, there is no obvious correlation between Fourier coefficients and luminosity for the whole magnitude range ($-12 {\gtrsim} M_{V} {\gtrsim} -22$); however, dwarfs display much higher dispersions at all radii.


title: Twist Angle-Dependent Bands and Valley Inversion in 2D Materials/hBN Heterostructures

abstract: The use of relative twist angle between adjacent atomic layers in a van der Waals heterostructure, has emerged as a new degree of freedom to tune electronic and optoelectronic properties of devices based on 2D materials. Using ABA-stacked trilayer (TLG) graphene as the model system, we show that, contrary to conventional wisdom, the band structures of 2D materials are systematically tunable depending on their relative alignment angle between hexagonal BN (hBN), even at very large twist angles. Moreover, addition or removal of the hBN substrate results in an inversion of the K and K' valley in TLG's lowest Landau level (LL). Our work illustrates the critical role played by substrates in van der Waals heterostructures and opens the door towards band structure modification and valley control via substrate and twist angle engineering.


title: Modulus metrics on networks

abstract: The concept of $p$-modulus gives a way to measure the richness of a family of objects on a graph. In this paper, we investigate the families of connecting walks between two fixed nodes and show how to use $p$-modulus to form a parametrized family of graph metrics that generalize several well-known and widely-used metrics. We also investigate a characteristic of metrics called the "antisnowflaking exponent" and present some numerical findings supporting a conjecture about the new metrics. We end with explicit computations of the new metrics on some selected graphs.


title: Geometric and LP-based heuristics for the quadratic travelling salesman problem

abstract: A generalization of the classical TSP is the so-called quadratic travelling salesman problem (QTSP), in which a cost coefficient is associated with the transition in every vertex, i.e. with every pair of edges traversed in succession. In this paper we consider two geometrically motivated special cases of the QTSP known from the literature, namely the angular-metric TSP, where transition costs correspond to turning angles in every vertex, and the angular-distance-metric TSP, where a linear combination of turning angles and Euclidean distances is considered. At first we introduce a wide range of heuristic approaches, motivated by the typical geometric structure of optimal solutions. In particular, we exploit lens-shaped neighborhoods of edges and a decomposition of the graph into layers of convex hulls, which are then merged into a tour by a greedy-type procedure or by utilizing an ILP model. Secondly, we consider an ILP model for a standard linearization of QTSP and compute fractional solutions of a relaxation. By rounding we obtain a collection of subtours, paths and isolated points, which are combined into a tour by various strategies, all of them involving auxiliary ILP models. Finally, different improvement heuristics are proposed, most notably a matheuristic which locally reoptimizes the solution for rectangular sectors of the given point set by an ILP approach. Extensive computational experiments for benchmark instances from the literature and extensions thereof illustrate the Pareto-efficient frontier of algorithms in a (running time, objective value)-space. It turns out that our new methods clearly dominate the previously published heuristics.


title: Network Traffic Driven Storage Repair

abstract: Recently we constructed an explicit family of locally repairable and locally regenerating codes. Their existence was proven by Kamath et al. but no explicit construction was given. Our design is based on HashTag codes that can have different sub-packetization levels. In this work we emphasize the importance of having two ways to repair a node: repair only with local parity nodes or repair with both local and global parity nodes. We say that the repair strategy is network traffic driven since it is in connection with the concrete system and code parameters: the repair bandwidth of the code, the number of I/O operations, the access time for the contacted parts and the size of the stored file. We show the benefits of having repair duality in one practical example implemented in Hadoop. We also give algorithms for efficient repair of the global parity nodes.


title: Nonperturbative Renormalization Group for the Landau-de Gennes Model

abstract: We studied the nematic isotropic phase transition by applying the functional renormalization group to the Landau-de Gennes model. We derived the flow equations for the effective potential as well as the cubic and quartic "couplings" and the anomalous dimension. We then solved the coupled flow equations on a grid using Newton Raphson method. A first order phase transition is observed. We also investigated the nematic isotropic puzzle (the NI puzzle) in this paper. We obtained the NI transition temperature difference ${T_c-T^*}=5.85K$ with sizable improvement over previous results.


title: Scoring Formulation for Multi-Condition Joint PLDA

abstract: The joint PLDA model, is a generalization of PLDA where the nuisance variable is no longer considered independent across samples, but potentially shared (tied) across samples that correspond to the same nuisance condition. The original work considered a single nuisance condition, deriving the EM and scoring formulas for this scenario. In this document, we show how to obtain likelihood ratios for scoring when multiple nuisance conditions are allowed in the model.


title: Lattice Diagrams of Knots and Diagrams of Lattice Stick Knots

abstract: We give a simple example showing that a knot or link diagram that lies in the ${\mathbb{Z}}^2$ lattice is not necessarily the projection of a lattice stick knot or link in the ${\mathbb{Z}}^3$ lattice, and we give a necessary and sufficient condition for when a knot or link diagram that lies in the ${\mathbb{Z}}^2$ lattice is in fact the projection of a lattice stick knot or link.