title: Robust and sparse k-means clustering for high-dimensional data

abstract: In real-world application scenarios, the identification of groups poses a significant challenge due to possibly occurring outliers and existing noise variables. Therefore, there is a need for a clustering method which is capable of revealing the group structure in data containing both outliers and noise variables without any pre-knowledge. In this paper, we propose a $k$-means-based algorithm incorporating a weighting function which leads to an automatic weight assignment for each observation. In order to cope with noise variables, a lasso-type penalty is used in an objective function adjusted by observation weights. We finally introduce a framework for selecting both the number of clusters and variables based on a modified gap statistic. The conducted experiments on simulated and real-world data demonstrate the advantage of the method to identify groups, outliers, and informative variables simultaneously.


title: Identifiability of linear compartmental models: the singular locus

abstract: This work addresses the problem of identifiability, that is, the question of whether parameters can be recovered from data, for linear compartmental models. Using standard differential algebra techniques, the question of whether a given model is generically locally identifiable is equivalent to asking whether the Jacobian matrix of a certain coefficient map, arising from input-output equations, is generically full rank. A natural next step is to study the set of parameter values where the Jacobian matrix drops in rank, which we refer to as the locus of non-identifiable parameter values, or, for short, the singular locus. In this work, we give a formula for coefficient maps in terms of acyclic subgraphs of the model's underlying directed graph and, then, study the case when the singular locus is defined by a single equation, the singular-locus equation. We prove that the singular-locus equation can be used to determine when submodels are generically locally identifiable. We also determine the singular-locus equation for two families of linear compartmental models, cycle and mammillary (star) models with input and output in a single compartment. We also state a conjecture for the corresponding equation for a third family: catenary (path) models. Finally, we introduce the identifiability degree, which is the number of parameter values that map to a generic input-output data vector. This degree was previously computed for mammillary and catenary models, and here we determine this degree for cycle models.


title: From (e,e'p) to Neutrino Scattering

abstract: Since a long time electron scattering has been envisaged as a powerful and preferential tool to investigate nuclear properties. In particular, the (e,e'p) knockout reaction has provided a wealth of information on the single particle (s.p.) aspects of nuclear structure, on the validity and the limit of the independent particle shell model. The work done for electron scattering is extremely useful also for the analysis and the interpretation of neutrino oscillation experiments, where nuclei are used as neutrino detectors and it is crucial that nuclear effects in neutrino-nucleus interactions are well under control. In this contribution it is discussed if and how the work done for (e,e'p) can be exploited for the analysis of neutrino-nucleus scattering data.


title: Analysis and mitigation of interface losses in trenched superconducting coplanar waveguide resonators

abstract: Improving the performance of superconducting qubits and resonators generally results from a combination of materials and fabrication process improvements and design modifications that reduce device sensitivity to residual losses. One instance of this approach is to use trenching into the device substrate in combination with superconductors and dielectrics with low intrinsic losses to improve quality factors and coherence times. Here we demonstrate titanium nitride coplanar waveguide resonators with mean quality factors exceeding two million and controlled trenching reaching 2.2 $\mu$m into the silicon substrate. Additionally, we measure sets of resonators with a range of sizes and trench depths and compare these results with finite-element simulations to demonstrate quantitative agreement with a model of interface dielectric loss. We then apply this analysis to determine the extent to which trenching can improve resonator performance.


title: Parallel lives: A local-realistic interpretation of "nonlocal" boxes

abstract: We carry out a thought experiment in an imaginary world. Our world is both local and realistic, yet it violates a Bell inequality more than does quantum theory. This serves to debunk the myth that equates local realism with local hidden variables in the simplest possible manner. Along the way, we reinterpret the celebrated 1935 argument of Einstein, Podolsky and Rosen, and come to the conclusion that they were right in their questioning the completeness of the Copenhagen version of quantum theory, provided one believes in a local-realistic universe. Throughout our journey, we strive to explain our views from first principles, without expecting mathematical sophistication nor specialized prior knowledge from the reader.


title: Critical behavior of SU(3) lattice gauge theory with 12 light flavors

abstract: It is expected that when the number of light flavors of gauge theories is increased near or beyond some critical value, new and interesting behavior occurs. We discuss the qualitative properties of the RG flows for a local $SU(3)$ theory with $N_f$ light fundamental flavors for $N_f$ near 12. We discuss the realization of the chiral symmetry and remind that the Wigner mode seems an unlikely alternative. We propose to use a linear sigma model to describe the light $\sigma$ masses found in recent lattice calculations for $N_f$=8 and 12. Progress made after the conference are reported in arXiv:1709.09264, where it is claimed that the breaking of the axial $U(1)_A$ is a key ingredient to get a small $\sigma$ mass. For SU(3) lattice gauge theory with 12 flavors of unimproved staggered fermions, the scaling of the imaginary part of the zeros of the partition function in the complex coupling plane is consistent with a first order phase transition for small values of the mass. Zech Gelzer and Diego Floor are investigating the scaling near the endpoint of the line of first order phase transition in the mass-coupling plane where a light and weakly interacting scalar is expected. We briefly discuss recent calculations of the second-order R\'enyi entanglement entropy and estimations of the central charge using the Tensor Renormalization Group method for the two-dimensional O(2) model and the possibility to extend these calculations to higher dimensions.


title: Non-equilibrium transport in the pseudospin-1 Dirac-Weyl system

abstract: Solid state materials hosting pseudospin-1 quasiparticles have attracted a great deal of recent attention. In these materials, the energy band contains of a pair of Dirac cones and a flat band through the connecting point of the cones. As the "caging" of carriers with a zero group velocity, the flat band itself has zero conductivity. However, in a non-equilibrium situation where a constant electric field is suddenly switched on, the flat band can enhance the resulting current in both the linear and nonlinear response regimes through distinct physical mechanisms. Using the ($2+1$) dimensional pseudospin-$1$ Dirac-Weyl system as a concrete setting, we demonstrate that, in the weak field regime, the interband current is about twice larger than that for pseudospin-1/2 system due to the interplay between the flat band and the negative band, with the scaling behavior determined by the Kubo formula. In the strong field regime, the intraband current is $\sqrt{2}$ times larger than that in the pseudospin-1/2 system, due to the additional contribution from particles residing in the flat band. In this case, the current and field follows the scaling law associated with Landau-Zener tunneling. These results provide a better understanding of the role of the flat band in non-equilibrium transport and are experimentally testable using electronic or photonic systems.


title: Topological $\mathbb{Z}_2$ Resonating-Valence-Bond Spin Liquid on the Square Lattice

abstract: A one-parameter family of long-range resonating valence bond (RVB) state on the square lattice was previously proposed to describe a critical spin liquid (SL) phase of the spin-$1/2$ frustrated Heisenberg model. We provide evidence that this RVB state in fact also realises a topological (long-range entangled) $\mathbb{Z}_2$ SL, limited by two transitions to critical SL phases. The topological phase is naturally connected to the $\mathbb{Z}_2$ gauge symmetry of the local tensor. This work shows that, on one hand, spin-$1/2$ topological SL with $C_{4v}$ point group symmetry and $SU(2)$ spin rotation symmetry exists on the square lattice and, on the other hand, criticality and nonbipartiteness are compatible. We also point out that, strong similarities between our phase diagram and the ones of classical interacting dimer models suggest both can be described by similar Kosterlitz-Thouless transitions. This scenario is further supported by the analysis of the one-dimensional boundary state.


title: Landau Damping with Electron Lenses in Space-Charge Dominated Beams

abstract: Progress on the Intensity Frontier of high energy physics critically depends on record high intensity charged particles accelerators. Beams in such machines become operationally limited by coherent beam instabilities, particularly enhanced in the regime of strong space charge (SC). Usual methods to control the instabilities, such as octupole magnets, beam feedback dampers and employment of chromatic effects, become less effective and insufficient. In [1] it was proposed to employ electron lenses for introduction of sufficient spread in particle oscillation frequencies needed for beam stabilization and in [2] it was shown that electron lenses are uniquely effective for Landau damping of transverse beam instabilities in high energy particle accelerators and their employment does not compromise incoherent (single particle) stability, dynamic aperture and the beam lifetime. Here we consider an important issue of effectiveness of the Landau damping with electron lenses in space-charge dominated beams and demonstrate that the desired stability can be assured with proper choice of the electron beam parameters and current distributions.


title: The distinguishing chromatic number of bipartite graphs of girth at least six

abstract: The distinguishing number $D(G)$ of a graph $G$ is the least integer $d$ such that $G$ has a vertex labeling with $d$ labels that is preserved only by a trivial automorphism. The distinguishing chromatic number $\chi_{D}(G)$ of $G$ is defined similarly, where, in addition, $f$ is assumed to be a proper labeling. Motivated by a conjecture in \cite{colins}, we prove that if $G$ is a bipartite graph of girth at least six with the maximum degree $\Delta (G)$, then $\chi_{D}(G)\leq \Delta (G)+1$. We also obtain an upper bound for $\chi_{D}(G)$ where $G$ is a graph with at most one cycle. Finally, we state a relationship between the distinguishing chromatic number of a graph and its spanning subgraphs.


title: Using a high-stability quartz-crystal microbalance to measure and model the chemical kinetics for gases in and on metals: oxygen in gold

abstract: This paper describes the use of a high-stability quartz-crystal microbalance (QCM) to measure the mass of a gas absorbed on and in the metal electrode on the quartz oscillator, when the gas pressure is low and the gas can be considered as rigidly attached to the metal, so viscosity effects are negligible. This provides an absolute measure of the total mass of gas uptake as a function of time, which can be used to model the kinetic processes involved. The technique can measure diffusion parameters of gases in metals close to room temperature at gas pressures much below one atmosphere, as relevant to surface processes such as atomic layer deposition and model studies of heterogeneous catalysis, whereas traditional diffusion measurements require temperatures over 400oC at gas pressures of at least a few Torr. A strong aspect of the method is the ability to combine the bulk measurement of absorbed mass by a QCM with a surface-sensitive technique such as Auger electron spectroscopy in the same vacuum chamber. The method is illustrated using atomic oxygen, formed under O2 gas at 6x10-5 Torr in the presence of a hot tungsten filament, interacting with the gold electrode on a QCM crystal held at 52 to 120oC. Some of the incident oxygen forms a surface oxide which eventually blocks more uptake, and the rest, about 80%, indiffuses. Surprisingly, the rate of oxygen uptake initially increases with the amount of oxygen previously absorbed; therefore, the measured oxygen uptake with time is reproducible only if pre-adsorption of oxygen conditions the sample. Temperatures above 130oC are necessary for measurable thermal desorption, but all the oxygen can be removed by CO scavenging at all temperatures of these experiments. Simple kinetic models are developed for fitting the experimental data to extract relevant parameters.


title: Zagier duality for level $p$ weakly holomorphic modular forms

abstract: We prove Zagier duality between the Fourier coefficients of canonical bases for spaces of weakly holomorphic modular forms of prime level $p$ with $11 \leq p \leq 37$ with poles only at the cusp at $\infty$, and special cases of duality for an infinite class of prime levels. We derive generating functions for the bases for genus 1 levels.


title: Estimation of Peer Effects in Endogenous Social Networks: Control Function Approach

abstract: We propose a method of estimating the linear-in-means model of peer effects in which the peer group, defined by a social network, is endogenous in the outcome equation for peer effects. Endogeneity is due to unobservable individual characteristics that influence both link formation in the network and the outcome of interest. We propose two estimators of the peer effect equation that control for the endogeneity of the social connections using a control function approach. We leave the functional form of the control function unspecified and treat it as unknown. To estimate the model, we use a sieve semiparametric approach, and we establish asymptotics of the semiparametric estimator.


title: Peccei-Quinn Relaxion

abstract: The relaxation mechanism, which solves the electroweak hierarchy problem without relying on TeV scale new physics, crucially depends on how a Higgs-dependent back-reaction potential is generated. In this paper, we suggest a new scenario in which the scalar potential induced by the QCD anomaly is responsible both for the relaxation mechanism and the Peccei-Quinn mechanism to solve the strong CP problem. The key idea is to introduce the relaxion and the QCD axion whose cosmic evolutions become quite different depending on an inflaton-dependent scalar potential. Our scheme raises the cutoff scale of the Higgs mass up to 10^7 GeV, and allows reheating temperature higher than the electroweak scale as would be required for viable cosmology. In addition, the QCD axion can account for the observed dark matter of the universe as produced by the conventional misalignment mechanism. We also consider the possibility that the couplings of the Standard Model depend on the inflaton and become stronger during inflation. In this case, the relaxation can be implemented with a sub-Planckian field excursion of the relaxion for a cutoff scale below 10 TeV.


title: Optimization and Performance of Bifacial Solar Modules: A Global Perspective

abstract: With the rapidly growing interest in bifacial photovoltaics (PV), a worldwide map of their potential performance can help assess and accelerate the global deployment of this emerging technology. However, the existing literature only highlights optimized bifacial PV for a few geographic locations or develops worldwide performance maps for very specific configurations, such as the vertical installation. It is still difficult to translate these location- and configuration-specific conclusions to a general optimized performance of this technology. In this paper, we present a global study and optimization of bifacial solar modules using a rigorous and comprehensive modeling framework. Our results demonstrate that with a low albedo of 0.25, the bifacial gain of ground-mounted bifacial modules is less than 10% worldwide. However, increasing the albedo to 0.5 and elevating modules 1 m above the ground can boost the bifacial gain to 30%. Moreover, we derive a set of empirical design rules, which optimize bifacial solar modules across the world, that provide the groundwork for rapid assessment of the location-specific performance. We find that ground-mounted, vertical, east-west-facing bifacial modules will outperform their south-north-facing, optimally tilted counterparts by up to 15% below the latitude of 30 degrees, for an albedo of 0.5. The relative energy output is the reverse of this in latitudes above 30 degrees. A detailed and systematic comparison with experimental data from Asia, Europe, and North America validates the model presented in this paper. An online simulation tool (https://nanohub.org/tools/pub) based on the model developed in this paper is also available for a user to predict and optimize bifacial modules in any arbitrary location across the globe.


title: A Rigorous Construction of the Supersymmetric Path Integral Associated to a Compact Spin Manifold

abstract: We give a rigorous construction of the path integral in N=1/2 supersymmetry as an integral map for differential forms on the loop space of a compact spin manifold. It is defined on the space of differential forms which can be represented by extended iterated integrals in the sense of Chen and Getzler-Jones-Petrack. Via the iterated integral map, we compare our path integral to the non-commutative loop space Chern character of G\"uneysu and the second author. Our theory provides a rigorous background to various formal proofs of the Atiyah-Singer index theorem using supersymmetric path integrals, as investigated by Alvarez-Gaum\'e, Atiyah, Bismut and Witten.


title: The Fermionic integral on loop space and the Pfaffian line bundle

abstract: As the loop space of a Riemannian manifold is infinite-dimensional, it is a non-trivial problem to make sense of the "top degree component" of a differential form on it. In this paper, we show that a formula from finite dimensions generalizes to assign a sensible "top degree component" to certain composite forms, obtained by wedging with the exponential (in the exterior algebra) of the canonical 2-form on the loop space. The result is a section on the Pfaffian line bundle on the loop space. We then identify this with a section of the line bundle obtained by transgression of the spin lifting gerbe. These results are a crucial ingredient for defining the fermionic part of the supersymmetric path integral on the loop space.


title: Sparse High-Dimensional Regression: Exact Scalable Algorithms and Phase Transitions

abstract: We present a novel binary convex reformulation of the sparse regression problem that constitutes a new duality perspective. We devise a new cutting plane method and provide evidence that it can solve to provable optimality the sparse regression problem for sample sizes n and number of regressors p in the 100,000s, that is two orders of magnitude better than the current state of the art, in seconds. The ability to solve the problem for very high dimensions allows us to observe new phase transition phenomena. Contrary to traditional complexity theory which suggests that the difficulty of a problem increases with problem size, the sparse regression problem has the property that as the number of samples $n$ increases the problem becomes easier in that the solution recovers 100% of the true signal, and our approach solves the problem extremely fast (in fact faster than Lasso), while for small number of samples n, our approach takes a larger amount of time to solve the problem, but importantly the optimal solution provides a statistically more relevant regressor. We argue that our exact sparse regression approach presents a superior alternative over heuristic methods available at present.


title: Sparse Hierarchical Regression with Polynomials

abstract: We present a novel method for exact hierarchical sparse polynomial regression. Our regressor is that degree $r$ polynomial which depends on at most $k$ inputs, counting at most $\ell$ monomial terms, which minimizes the sum of the squares of its prediction errors. The previous hierarchical sparse specification aligns well with modern big data settings where many inputs are not relevant for prediction purposes and the functional complexity of the regressor needs to be controlled as to avoid overfitting. We present a two-step approach to this hierarchical sparse regression problem. First, we discard irrelevant inputs using an extremely fast input ranking heuristic. Secondly, we take advantage of modern cutting plane methods for integer optimization to solve our resulting reduced hierarchical $(k, \ell)$-sparse problem exactly. The ability of our method to identify all $k$ relevant inputs and all $\ell$ monomial terms is shown empirically to experience a phase transition. Crucially, the same transition also presents itself in our ability to reject all irrelevant features and monomials as well. In the regime where our method is statistically powerful, its computational complexity is interestingly on par with Lasso based heuristics. The presented work fills a void in terms of a lack of powerful disciplined nonlinear sparse regression methods in high-dimensional settings. Our method is shown empirically to scale to regression problems with $n\approx 10,000$ observations for input dimension $p\approx 1,000$.


title: Relations between permutation representations in positive characteristic

abstract: Given a finite group G and a field F, a G-set X gives rise to an F[G]-permutation module F[X]. This defines a map from the Burnside ring of G to its representation ring over F. It is an old problem in representation theory, with wide-ranging applications in algebra, number theory, and geometry, to give explicit generators of the kernel K_F(G) of this map, i.e. to classify pairs of G-sets X, Y such that F[X] is isomorphic to F[Y]. When F has characteristic 0, a complete description of K_F(G) is now known. In this paper, we give a similar description of K_F(G) when F is a field of characteristic p>0 in all but the most complicated case, which is when G has a subquotient that is a non-p-hypo-elementary (p,p)-Dress group.


title: High-energy resummation in heavy-quark pair photoproduction

abstract: We present our predictions for the inclusive production of two heavy quark-antiquark pairs, separated by a large rapidity interval, in the collision of (quasi-)real photons at the energies of LEP2 and of some future electron-positron colliders. We include in our calculation the full resummation of leading logarithms in the center-of-mass energy and a partial resummation of the next-to-leading logarithms, within the Balitsky-Fadin-Kuraev-Lipatov (BFKL) approach.


title: Nonuniqueness of weak solutions to the Navier-Stokes equation

abstract: For initial datum of finite kinetic energy, Leray has proven in 1934 that there exists at least one global in time finite energy weak solution of the 3D Navier-Stokes equations. In this paper we prove that weak solutions of the 3D Navier-Stokes equations are not unique in the class of weak solutions with finite kinetic energy. Moreover, we prove that Holder continuous dissipative weak solutions of the 3D Euler equations may be obtained as a strong vanishing viscosity limit of a sequence of finite energy weak solutions of the 3D Navier-Stokes equations.


title: Topological critical slowing down: variations on a toy model

abstract: Numerical simulations of lattice quantum field theories whose continuum counterparts possess classical solutions with non-trivial topology face a severe critical slowing down as the continuum limit is approached. Standard Monte-Carlo algorithms develop a loss of ergodicity, with the system remaining frozen in configurations with fixed topology. We analyze the problem in a simple toy model, consisting of the path integral formulation of a quantum mechanical particle constrained to move on a circumference. More specifically, we implement for this toy model various techniques which have been proposed to solve or alleviate the problem for more complex systems, like non-abelian gauge theories, and compare them both in the regime of low temperature and in that of very high temperature. Among the various techniques, we consider also a new algorithm which completely solves the freezing problem, but unfortunately is specifically tailored for this particular model and not easily exportable to more complex systems.


title: Prominence eruption initiated by helical kink-instability of an embedded flux rope

abstract: We study the triggering mechanism of a limb-prominence eruption and the associated coronal mass ejection near AR 12342 using SDO and LASCO/SOHO observations. The prominence is seen with an embedded flux thread (FT) at one end and bifurcates from the middle to a different footpoint location. The morphological evolution of the FT is similar to an unstable flux rope (FR), which we regard as prominence embedded FR. The FR twist exceeds the critical value. In addition, the morphology of the prominence plasma in 304\AA~images marks the helical nature of the magnetic skeleton with a total of 2.96 turns along arc length. The potential field extrapolation model indicates that the critical height of the background magnetic field gradient falls within the inner corona (105Mm) consistent with the extent of coronal plasma loops. These results suggest that the helical kink instability in the embedded FR caused the slow rise of the prominence to a height of the torus instability domain. Moreover, the differential emission measure analysis unveils heating of the prominence plasma to coronal temperatures during eruption, suggesting a reconnection-related heating underneath the upward rising embedded FR. The prominence starts with a slow rise motion of 10km/s, followed by fast and slow acceleration phases having an average acceleration of $28.9m/s^2$, $2.4m/s^2$ in C2, C3 field of view respectively. As predicted by previous numerical simulations, the observed synchronous kinematic profiles of the CME leading edge and the core supports the involved FR instability in the prominence initiation.


title: Distributed Robust Set-Invariance for Interconnected Linear Systems

abstract: We introduce a class of distributed control policies for networks of discrete-time linear systems with polytopic additive disturbances. The objective is to restrict the network-level state and controls to user-specified polyhedral sets for all times. This problem arises in many safety-critical applications. We consider two problems. First, given a communication graph characterizing the structure of the information flow in the network, we find the optimal distributed control policy by solving a single linear program. Second, we find the sparsest communication graph required for the existence of a distributed invariance-inducing control policy. Illustrative examples, including one on platooning, are presented.


title: Integer Forcing: Effective SNR Distribution and Practical Block-Based Schemes

abstract: Integer Forcing (IF) is a novel linear receiver architecture, where instead of separating the codewords sent by each transmitter, and decoding them individually, forces integer-valued linear combinations at each receive antenna, and decodes the linear combinations. The original codewords are obtained by inverting the integer-valued matrix. While demonstrating superior performance, IF requires complex optimization in order to find the optimal linear combinations, and demands either using multi-level nested lattice codes, or reducing the rates of all transmitters to equal the weakest one. Finally, the distribution of the resulting effective SNR is hard to evaluate. In this paper, we first give simple upper and lower bounds on the effective SNR of a single linear combination in IF. These expressions allow us to easily bound the distribution of the effective SNR for any given linear combination used. We then suggest two simple block-based IF schemes. These schemes, while sub-optimal, significantly reduce the complexity of the optimization process, and, more importantly, do not require reducing the rates of all transmitters, as decoding is done block-wise. Finally, we bound the distribution of the effective SNR of the decoding schemes, and show via simulations the superiority of block-wise schemes at low SNR.


title: Estimation of Graphical Models using the $L_{1,2}$ Norm

abstract: Gaussian graphical models are recently used in economics to obtain networks of dependence among agents. A widely-used estimator is the Graphical Lasso (GLASSO), which amounts to a maximum likelihood estimation regularized using the $L_{1,1}$ matrix norm on the precision matrix $\Omega$. The $L_{1,1}$ norm is a lasso penalty that controls for sparsity, or the number of zeros in $\Omega$. We propose a new estimator called Structured Graphical Lasso (SGLASSO) that uses the $L_{1,2}$ mixed norm. The use of the $L_{1,2}$ penalty controls for the structure of the sparsity in $\Omega$. We show that when the network size is fixed, SGLASSO is asymptotically equivalent to an infeasible GLASSO problem which prioritizes the sparsity-recovery of high-degree nodes. Monte Carlo simulation shows that SGLASSO outperforms GLASSO in terms of estimating the overall precision matrix and in terms of estimating the structure of the graphical model. In an empirical illustration using a classic firms' investment dataset, we obtain a network of firms' dependence that exhibits the core-periphery structure, with General Motors, General Electric and U.S. Steel forming the core group of firms.


title: Answering UCQs under updates and in the presence of integrity constraints

abstract: We investigate the query evaluation problem for fixed queries over fully dynamic databases where tuples can be inserted or deleted. The task is to design a dynamic data structure that can immediately report the new result of a fixed query after every database update. We consider unions of conjunctive queries (UCQs) and focus on the query evaluation tasks testing (decide whether an input tuple belongs to the query result), enumeration (enumerate, without repetition, all tuples in the query result), and counting (output the number of tuples in the query result). We identify three increasingly restrictive classes of UCQs which we call t-hierarchical, q-hierarchical, and exhaustively q-hierarchical UCQs. Our main results provide the following dichotomies: If the query's homomorphic core is t-hierarchical (q-hierarchical, exhaustively q-hierarchical), then the testing (enumeration, counting) problem can be solved with constant update time and constant testing time (delay, counting time). Otherwise, it cannot be solved with sublinear update time and sublinear testing time (delay, counting time), unless the OV-conjecture and/or the OMv-conjecture fails. We also study the complexity of query evaluation in the dynamic setting in the presence of integrity constraints, and we obtain according dichotomy results for the special case of small domain constraints (i.e., constraints which state that all values in a particular column of a relation belong to a fixed domain of constant size).


title: Persistence, coexistence and extinction in two species chemotaxis models on bounded heterogeneous environments

abstract: In this paper, we consider a two species chemotaxis system of parabolic-parabolic-elliptic type with Lotka-Volterra type competition terms in heterogeneous media. We first find various conditions on the parameters which guarantee the global existence and boundedness of classical solutions with nonnegative initial functions. Next, we find further conditions on the parameters which establish the persistence of the two species. Then, under the same set of conditions for the persistence of two species, we prove the existence of coexistence states. Finally we prove the extinction phenomena in the sense that one of the species dies out asymptotically and the other reaches its carrying capacity as time goes to infinity. The persistence in general two species chemotaxis systems is studied for the first time. Several important techniques are developed to study the persistence and coexistence of the two species chemotaxis systems. Many existing results on the persistence, coexistence, and extinction on two species competition systems without chemotaxis are recovered.


title: Bayesian Multi Plate High Throughput Screening of Compounds

abstract: High throughput screening of compounds (chemicals) is an essential part of drug discovery [7], involving thousands to millions of compounds, with the purpose of identifying candidate hits. Most statistical tools, including the industry standard B-score method, work on individual compound plates and do not exploit cross-plate correlation or statistical strength among plates. We present a new statistical framework for high throughput screening of compounds based on Bayesian nonparametric modeling. The proposed approach is able to identify candidate hits from multiple plates simultaneously, sharing statistical strength among plates and providing more robust estimates of compound activity. It can flexibly accommodate arbitrary distributions of compound activities and is applicable to any plate geometry. The algorithm provides a principled statistical approach for hit identification and false discovery rate control. Experiments demonstrate significant improvements in hit identification sensitivity and specificity over the B-score method, which is highly sensitive to threshold choice. The framework is implemented as an efficient R extension package BHTSpack and is suitable for large scale data sets.


title: A counterexample regarding labelled well-quasi-ordering

abstract: Korpelainen, Lozin, and Razgon conjectured that a hereditary property of graphs which is well-quasi-ordered by the induced subgraph order and defined by only finitely many minimal forbidden induced subgraphs is labelled well-quasi-ordered, a notion stronger than that of $n$-well-quasi-order introduced by Pouzet in the 1970s. We present a counterexample to this conjecture. In fact, we exhibit a hereditary property of graphs which is well-quasi-ordered by the induced subgraph order and defined by finitely many minimal forbidden induced subgraphs yet is not $2$-well-quasi-ordered. This counterexample is based on the widdershins spiral, which has received some study in the area of permutation patterns.


title: Revealing new dynamical patterns in a reaction-diffusion model with cyclic competition via a novel computational framework

abstract: Understanding how patterns and travelling waves form in chemical and biological reaction-diffusion models is an area which has been widely researched, yet is still experiencing fast development. Surprisingly enough, we still do not have a clear understanding about all possible types of dynamical regimes in classical reaction-diffusion models such as Lotka-Volterra competition models with spatial dependence. In this work, we demonstrate some new types of wave propagation and pattern formation in a classical three species cyclic competition model with spatial diffusion, which have been so far missed in the literature. These new patterns are characterised by a high regularity in space, but are different from patterns previously known to exist in reaction-diffusion models, and may have important applications in improving our understanding of biological pattern formation and invasion theory. Finding these new patterns is made technically possible by using an automatic adaptive finite element method driven by a novel a posteriori error estimate which is proven to provide a reliable bound for the error of the numerical method. We demonstrate how this numerical framework allows us to easily explore the dynamical patterns both in two and three spatial dimensions.


title: The force-free twisted magnetosphere of a neutron star - II: Degeneracies of the Grad-Shafranov equation

abstract: We extend our previous study of equilibrium solutions of non-rotating force-free magnetospheres of neutron stars. We show that multiple solutions exist for the same sets of parameters, implying that the solutions are degenerate. We are able to obtain configurations with disconnected field lines, however, in nearly all cases these correspond to degenerate higher energy solutions. We carry out a wide parametric search in order to understand the properties of the solutions. We confirm our previous results that the lower energy solutions have up to $\sim 25\%$ more energy than the vacuum case, helicity of the order of $\sim 5$ (in some defined units), maximum twist of $\sim 1.5$ rad, and a dipole strength that is up to $\sim 40\%$ larger than the vacuum dipole. Including the degenerate higher energy solutions allows for larger theoretical limits of up to $\sim 80\%$ more energy with respect to the vacuum case, helicity of the order of $\sim 8$, and a dipole strength that can now be up to four times that of the vacuum dipole, while the twist can be significantly larger and even diverge for configurations with disconnected domains. The higher energy solutions are probably unstable, therefore, it is unlikely that such magnetospheres exist under normal conditions in magnetars and high magnetic field pulsars.


title: Effective learning is accompanied by high dimensional and efficient representations of neural activity

abstract: A fundamental cognitive process is the ability to map value and identity onto objects as we learn about them. Exactly how such mental constructs emerge and what kind of space best embeds this mapping remains incompletely understood. Here we develop tools to quantify the space and organization of such a mapping, thereby providing a framework for studying the geometric representations of neural responses as reflected in functional MRI. Considering how human subjects learn the values of novel objects, we show that quick learners have a higher dimensional geometric representation than slow learners, and hence more easily distinguishable whole-brain responses to objects of different value. Furthermore, we find that quick learners display a more compact embedding of their neural responses and hence have a higher ratio of their task-based dimension to their embedding dimension -- consistent with a greater efficiency of cognitive coding. Lastly, we investigate the neurophysiological drivers of high dimensional patterns at both regional and voxel levels, and we complete our study with a complementary test of the distinguishability of associated whole-brain responses. Our results demonstrate a spatial organization of neural responses characteristic of learning, and offer a suite of geometric measures applicable to the study of efficient coding in higher-order cognitive processes more broadly.


title: Integrating quantum key distribution with classical communications in backbone fiber network

abstract: Quantum key distribution (QKD) provides information-theoretic security based on the laws of quantum mechanics. The desire to reduce costs and increase robustness in real-world applications has motivated the study of coexistence between QKD and intense classical data traffic in a single fiber. Previous works on coexistence in metropolitan areas have used wavelength-division multiplexing, however, coexistence in backbone fiber networks remains a great experimental challenge, as Tbps data of up to 20 dBm optical power is transferred, and much more noise is generated for QKD. Here we present for the first time, to the best of our knowledge, the integration of QKD with a commercial backbone network of 3.6 Tbps classical data at 21 dBm launch power over 66 km fiber. With 20 GHz pass-band filtering and large effective core area fibers, real-time secure key rates can reach 4.5 kbps and 5.1 kbps for co-propagation and counter-propagation at the maximum launch power, respectively. This demonstrates feasibility and represents an important step towards building a quantum network that coexists with the current backbone fiber infrastructure of classical communications.


title: A degenerate three-level laser coupled to a squeezed vacuum reservoir

abstract: Employing the master equation for a three-level laser driven by coherent light and coupled to a squeezed vacuum reservoir, we obtain stochastic differential equations associated with the normal ordering. Using the solutions of the stochastic differential equations, we calculate the quadrature variance, the squeezing spectrum, the mean photon number, and the variance of the photon number. It turns out that the degree of squeezing increases with the linear gain coefficient or the squeeze parameter. It is also found that the driving coherent light decreases the mean photon number.


title: Notes on rate equations in nonlinear continuum mechanics

abstract: The paper gives an introduction to rate equations in nonlinear continuum mechanics which should obey specific transformation rules. Emphasis is placed on the geometrical nature of the operations involved in order to clarify the different concepts. The paper is particularly concerned with common classes of constitutive equations based on corotational stress rates and their proper implementation in time for solving initial boundary value problems. Hypoelastic simple shear is considered as an example application for the derived theory and algorithms.


title: Macroscopic Schoen conjecture for manifolds with non-zero simplicial volume

abstract: We prove that given a hyperbolic manifold endowed with an auxiliary Riemannian metric whose sectional curvature is negative and whose volume is sufficiently small in comparison to the hyperbolic one, we can always find for any radius at least $1$ a ball in its universal cover whose volume is bigger than the hyperbolic one. This result is deduced from a non-sharp macroscopic version of a conjecture by R. Schoen about scalar curvature, whose proof is a variation of an argument due to M. Gromov and based on a smoothing technique. We take the opportunity of this work to present a full account of this technique which involves simplicial volume and deserves to be better known.


title: Reducing resources for verification of quantum computations

abstract: We present two verification protocols where the correctness of a "target" computation is checked by means of "trap" computations that can be efficiently simulated on a classical computer. Our protocols rely on a minimal set of noise-free operations (preparation of eight single-qubit states or measurement of four observables, both on a single plane of the Bloch sphere) and achieve linear overhead. To the best of our knowledge, our protocols are the least demanding techniques able to achieve linear overhead. They represent a step towards further reducing the quantum requirements for verification.


title: Prospects for Improved Understanding of Isotopic Reactor Antineutrino Fluxes

abstract: Predictions of antineutrino fluxes produced by fission isotopes in a nuclear reactor have recently received increased scrutiny due to observed differences in predicted and measured inverse beta decay (IBD) yields, referred to as the 'reactor antineutrino flux anomaly.' In this paper, global fits are applied to existing IBD yield measurements to produce constraints on antineutrino production by individual plutonium and uranium fission isotopes. We find that fits including measurements from highly 235U-enriched cores and fits including Daya Bay's new fuel evolution result produce discrepant best-fit IBD yields for 235U and 239Pu. This discrepancy can be alleviated in a global analysis of all datasets through simultaneous fitting of 239Pu, 235U, and 238U yields. The measured IBD yield of U238 in this analysis is (7.02 +/- 1.65) x 10^-43 cm2/fission, nearly two standard deviations below existing predictions. Future hypothetical IBD yield measurements by short-baseline reactor experiments are examined to determine their possible impact on global understanding of isotopic IBD yields. It is found that future improved short-baseline IBD yield measurements at both high-enriched and low-enriched cores can significantly improve constraints for 235U, 238U, and 239Pu, providing comparable or superior precision to existing conversion- and summation-based antineutrino flux predictions. Systematic and experimental requirements for these future measurements are also investigated.


title: Adding successor: A transfer theorem for separation and covering

abstract: Given a class C of word languages, the C-separation problem asks for an algorithm that, given as input two regular languages, decides whether there exists a third language in C containing the first language, while being disjoint from the second. Separation is usually investigated as a means to obtain a deep understanding of the class C. In the paper, we are mainly interested in classes defined by logical formalisms. Such classes are often built on top of each other: given some logic, one builds a stronger one by adding new predicates to its signature. A natural construction is to enrich a logic with the successor relation. In this paper, we present a transfer result applying to this construction: we show that for suitable logically defined classes, separation for the logic enriched with the successor relation reduces to separation for the original logic. Our theorem also applies to a problem that is stronger than separation: covering. Moreover, we actually present two reductions: one for languages of finite words and the other for languages of infinite words.


title: Graph Convolutional Networks for Named Entity Recognition

abstract: In this paper we investigate the role of the dependency tree in a named entity recognizer upon using a set of GCN. We perform a comparison among different NER architectures and show that the grammar of a sentence positively influences the results. Experiments on the ontonotes dataset demonstrate consistent performance improvements, without requiring heavy feature engineering nor additional language-specific knowledge.


title: Semi-Inclusive $\pi_0$ target and beam-target asymmetries from 6 GeV electron scattering with CLAS

abstract: We present precision measurements of the target and beam-target spin asymmetries from neutral pion electroproduction in deep-inelastic scattering (DIS) using the CEBAF Large Acceptance Spectrometer (CLAS) at Jefferson Lab. We scattered 6-GeV, longitudinally polarized electrons off longitudinally polarized protons in a cryogenic $^{14}$NH$_3$ target, and extracted double and single target spin asymmetries for $ep\rightarrow e^\prime\pi^0X$ in multidimensional bins in four-momentum transfer ($1.0<Q^2<3.2$ GeV$^2$), Bjorken-$x$ ($0.12<x<0.48$), hadron energy fraction ($0.4<z<0.7$), transverse pion momentum ($0<P_T<1.0$ GeV), and azimuthal angle $\phi_h$ between the lepton scattering and hadron production planes. We extracted asymmetries as a function of both $x$ and $P_T$, which provide access to transverse-momentum distributions of longitudinally polarized quarks. The double spin asymmetries depend weakly on $P_T$. The $\sin 2\phi_h$ moments are zero within uncertainties, which is consistent with the expected suppression of the Collins fragmentation function. The observed $\sin\phi_h$ moments suggest that quark gluon correlations are significant at large $x$.


title: Versatile engineering of multimode squeezed states by optimizing the pump spectral profile in spontaneous parametric down-conversion

abstract: We study the quantum correlations induced by spontaneous parametric down-conversion (SPDC) of a frequency comb. We derive a theoretical method to find the output state corresponding to a pump with an arbitrary spectral profile. After applying it to the relevant example of a spectrally chirped pump, we run an optimization algorithm to numerically find the pump profiles maximizing some target functions. These include the number of independently squeezed modes and the variances of nullifiers defining cluster states used in many continuous-variable quantum information protocols. To assess the advantages of pump-shaping in real experiments we take into account the physical limitations of the pulse shaper.


title: Introducing DeepBalance: Random Deep Belief Network Ensembles to Address Class Imbalance

abstract: Class imbalance problems manifest in domains such as financial fraud detection or network intrusion analysis, where the prevalence of one class is much higher than another. Typically, practitioners are more interested in predicting the minority class than the majority class as the minority class may carry a higher misclassification cost. However, classifier performance deteriorates in the face of class imbalance as oftentimes classifiers may predict every point as the majority class. Methods for dealing with class imbalance include cost-sensitive learning or resampling techniques. In this paper, we introduce DeepBalance, an ensemble of deep belief networks trained with balanced bootstraps and random feature selection. We demonstrate that our proposed method outperforms baseline resampling methods such as SMOTE and under- and over-sampling in metrics such as AUC and sensitivity when applied to a highly imbalanced financial transaction data. Additionally, we explore performance and training time implications of various model parameters. Furthermore, we show that our model is easily parallelizable, which can reduce training times. Finally, we present an implementation of DeepBalance in R.


title: Natural Time, Nowcasting and the Physics of Earthquakes: Estimation of Seismic Risk to Global Megacities

abstract: This paper describes the use of the idea of natural time to propose a new method for characterizing the seismic risk to the world's major cities at risk of earthquakes. Rather than focus on forecasting, which is the computation of probabilities of future events, we define the term seismic nowcasting, which is the computation of the current state of seismic hazard in a defined geographic region.


title: A DMRG study: Pair-density wave in spin-valley locked systems

abstract: Interest in modulated paired states, long sought since the first proposals by Fulde and Ferrell and by Larkin and Ovchinnikov, has grown recently in the context of strongly coupled superconductors under the name of pair density wave. However, there is little theoretical understanding of how such a state might arise out of strong coupling physics in simple models. Although density matrix renormalization group has been a powerful tool for exploring strong coupling modulation phenomena of spin and charge stripe in the Hubbard model and the t-J model, there has been no numerical evidence of PDW within these models using DMRG. Here we note that a system with inversion breaking, C3v point group symmetry may host a PDW-like state. Motivated by the fact that spin-valley locked band structure of hole-doped group VI transition metal dichalcogenides materializes such a setting, we use DMRG to study the superconducting tendencies in spin-valley locked systems with strong short-ranged repulsion. Remarkably we find robust evidence for a PDW and the first of such evidence within DMRG studies of a simple fermionic model.


title: Doppler shift oscillations from a hot line observed by IRIS

abstract: We present a detailed investigation of the Doppler shift oscillations in a hot loop during an M7.1 flare on 2014 October 27 observed by the Interface Region Imaging Spectrograph. The periodic oscillations are observed in the Doppler shift of Fe XXI 1354.09 A (logT~7.05), and the dominant period is about 3.1 minutes. However, such 3.1-min oscillations are not found in the line-integrated intensity of Fe XXI 1354.09 A, AIA EUV fluxes, or microwave emissions. SDO/AIA and Hinode/XRT imaging observations indicate that the Doppler shift oscillations locate at the hot loop-top region (> 11MK). Moreover, the differential emission measure (DEM) results show that the temperature is increasing rapidly when the Doppler shift oscillates, but the number density does not exhibit the corresponding increases nor oscillations, implying that the flare loop is likely to oscillate in an incompressible mode. All these facts suggest that the Doppler shift oscillations at the shorter period are most likely the standing kink oscillations in a flare loop. Meanwhile, a longer period of about 10 minutes is identified in the time series of Doppler shift and line-integrated intensity, GOES SXR fluxes and AIA EUV light curves, indicating the periodic energy release in this flare, which may be caused by a slow mode wave.


title: Measuring the Eccentricity of Items

abstract: The long-tail phenomenon tells us that there are many items in the tail. However, not all tail items are the same. Each item acquires different kinds of users. Some items are loved by the general public, while some items are consumed by eccentric fans. In this paper, we propose a novel metric, item eccentricity, to incorporate this difference between consumers of the items. Eccentric items are defined as items that are consumed by eccentric users. We used this metric to analyze two real-world datasets of music and movies and observed the characteristics of items in terms of eccentricity. The results showed that our defined eccentricity of an item does not change much over time, and classified eccentric and noneccentric items present significantly distinct characteristics. The proposed metric effectively separates the eccentric and noneccentric items mixed in the tail, which could not be done with the previous measures, which only consider the popularity of items.


title: Active Information Acquisition for Linear Optimization

abstract: We consider partially-specified optimization problems where the goal is to actively, but efficiently, acquire missing information about the problem in order to solve it. An algorithm designer wishes to solve a linear program (LP), $\max \mathbf{c}^T \mathbf{x}$ s.t. $\mathbf{A}\mathbf{x} \leq \mathbf{b}, \mathbf{x} \ge \mathbf{0}$, but does not initially know some of the parameters. The algorithm can iteratively choose an unknown parameter and gather information in the form of a noisy sample centered at the parameter's (unknown) value. The goal is to find an approximately feasible and optimal solution to the underlying LP with high probability while drawing a small number of samples. We focus on two cases. (1) When the parameters $\mathbf{c}$ of the objective are initially unknown, we take an information-theoretic approach and give roughly matching upper and lower sample complexity bounds, with an (inefficient) successive-elimination algorithm. (2) When the parameters $\mathbf{b}$ of the constraints are initially unknown, we propose an efficient algorithm combining techniques from the ellipsoid method for LP and confidence-bound approaches from bandit algorithms. The algorithm adaptively gathers information about constraints only as needed in order to make progress. We give sample complexity bounds for the algorithm and demonstrate its improvement over a naive approach via simulation.


title: Status of the TREX-DM experiment at the Canfranc Underground Laboratory

abstract: The TREX-DM experiment is conceived to look for low mass WIMPs by means of a gas time projection chamber equipped with novel micromegas readout planes at the Canfranc Underground Laboratory. The detector can hold 20 l of pressurized gas up to 10 bar, which corresponds to 0.30 kg of Ar, or alternatively, 0.16 kg of Ne. The micromegas will be read with a self-triggered acquisition, allowing for effective thresholds below 0.4 keV (electron equivalent). The preliminary background model, following a complete material screening program, points to levels of the order of 1-10 counts keV-1 kg-1 d-1 in the region of interest, making TREX-DM competitive. The status of the commissioning, description of the background model and the corresponding WIMP sensitivity will be presented here.


title: Finding Small Weight Isomorphisms with Additional Constraints is Fixed-Parameter Tractable

abstract: Lubiw showed that several variants of Graph Isomorphism are NP-complete, where the solutions are required to satisfy certain additional constraints [SICOMP 10, 1981]. One of these, called Isomorphism With Restrictions, is to decide for two given graphs $X_1=(V,E_1)$ and $X_2=(V,E_2)$ and a subset $R\subseteq V\times V$ of forbidden pairs whether there is an isomorphism $\pi$ from $X_1$ to $X_2$ such that $\pi(i)\neq j$ for all $(i,j)\in R$. We prove that this problem and several of its generalizations are in fact in FPT: - The problem of deciding whether there is an isomorphism between two graphs that moves k vertices and satisfies Lubiw-style constraints is in FPT, with k and the size of $R$ as parameters. The problem remains in FPT if a CNF of such constraints is allowed. It follows that the problem to decide whether there is an isomorphism that moves exactly k vertices is in FPT. This solves a question left open in our article on exact weight automorphisms [STACS 2017]. - When the weight and complexity are unrestricted, finding isomorphisms that satisfy a CNF of Lubiw-style constraints can be solved in FPT with access to a GI oracle. - Checking if there is an isomorphism $\pi$ between two graphs with complexity t is also in FPT with t as parameter, where the complexity of a permutation is the Cayley measure defined as the minimum number t such that $\pi$ can be expressed as a product of t transpositions. - We consider a more general problem in which the vertex set of a graph X is partitioned into Red and Blue, and we are interested in an automorphism that stabilizes Red and Blue and moves exactly k vertices in Blue, where k is the parameter. This problem was introduced by [Downey and Fellows 1999], and we showed [STACS 2017] that it is W[1]-hard even with color classes of size 4 inside Red. Now, for color classes of size at most 3 inside Red, we show the problem is in FPT.


title: Universal entanglement timescale for R\'enyi entropies

abstract: Recently it was shown that the growth of entanglement in an initially separable state, as measured by the purity of subsystems, can be characterized by a timescale that takes a universal form for any Hamiltonian. We show that the same timescale governs the growth of entanglement for all R\'enyi entropies. Since the family of R\'enyi entropies completely characterizes the entanglement of a pure bipartite state, our timescale is a universal feature of bipartite entanglement. The timescale depends only on the interaction Hamiltonian and the initial state.


title: An Axiomatic Study of Scoring Rule Markets

abstract: Prediction markets are well-studied in the case where predictions are probabilities or expectations of future random variables. In 2008, Lambert, et al. proposed a generalization, which we call "scoring rule markets" (SRMs), in which traders predict the value of arbitrary statistics of the random variables, provided these statistics can be elicited by a scoring rule. Surprisingly, despite active recent work on prediction markets, there has not yet been any investigation into the properties of more general SRMs. To initiate such a study, we ask the following question: in what sense are SRMs "markets"? We classify SRMs according to several axioms that capture potentially desirable qualities of a market, such as the ability to freely exchange goods (contracts) for money. Not all SRMs satisfy our axioms: once a contract is purchased in any market for prediction the median of some variable, there will not necessarily be any way to sell that contract back, even in a very weak sense. Our main result is a characterization showing that slight generalizations of cost-function-based markets are the only markets to satisfy all of our axioms for finite-outcome random variables. Nonetheless, we find that several SRMs satisfy weaker versions of our axioms, including a novel share-based market mechanism for ratios of expected values.


title: Empirical Bayes Shrinkage and False Discovery Rate Estimation, Allowing For Unwanted Variation

abstract: We combine two important ideas in the analysis of large-scale genomics experiments (e.g. experiments that aim to identify genes that are differentially expressed between two conditions). The first is use of Empirical Bayes (EB) methods to handle the large number of potentially-sparse effects, and estimate false discovery rates and related quantities. The second is use of factor analysis methods to deal with sources of unwanted variation such as batch effects and unmeasured confounders. We describe a simple modular fitting procedure that combines key ideas from both these lines of research. This yields new, powerful EB methods for analyzing genomics experiments that account for both sparse effects and unwanted variation. In realistic simulations, these new methods provide significant gains in power and calibration over competing methods. In real data analysis we find that different methods, while often conceptually similar, can vary widely in their assessments of statistical significance. This highlights the need for care in both choice of methods and interpretation of results. All methods introduced in this paper are implemented in the R package vicar available at https://github.com/dcgerard/vicar .


title: Universal and approximate relations for the gravitational-wave damping timescale of $f$-modes in neutron stars

abstract: Existing estimates of the gravitational-wave damping timescale of the dominant quadrupole oscillation mode in the case of rapidly rotating stars are based on using a Newtonian estimate for the energy of the mode, in combination with the lowest-order post-Newtonian quadrupole formula for estimating the gravitational-wave luminosity. We investigate a number of other choices for estimating the gravitational-wave damping timescale in the nonrotating limit and construct a highly accurate, empirically corrected formula that has a maximum relative error of only 3% with respect to the perturbative result in full general relativity. The expressions involved are sufficiently general to be extended to the case of rapidly rotating stars. We also present a new higher-order empirical relation for the gravitational-wave damping timescale of quadrupole oscillations that is accurate in the whole range of expected values for the compactness of neutron stars, without the need for involving the moment of inertia.


title: Hyperfine Structure and Abundances of Heavy Elements in 68 Tauri (HD 27962)

abstract: HD 27962, also known as 68 Tauri, is a Chemically Peculiar Am star member of the Hyades Open Cluster in the local arm of the Galaxy. We have modeled the high resolution SOPHIE (R=75000) spectrum of 68 Tauri using updated model atmosphere and spectrum synthesis to derive chemical abundances in its atmosphere. In particular, we have studied the effect of the inclusion of Hyperfine Structure of various Baryum isotopes on the determination of the Baryum abundance in 68 Tauri. We have also derived new abundances using updated accurate atomic parameters retrieved from the NIST database.


title: An Extended Analytic Model for the Heating of Bondwires

abstract: We present an extended analytic formula for the calculation of the temperature profile along a bondwire embedded in a package. The resulting closed formula is built by coupling the heat transfer equations of the bondwire and the surrounding moulding compound by means of auxiliary variables that stem from an \emph{ad-hoc} linearisation and mediate the wire-mould thermal interaction. The model, which corrects typical simplifications in previously introduced analytic models, is also optimised against carefully taken experimental samples representing fusing events of bondwires within real packages.


title: Resonating Valence Bond Theory of Superconductivity: Beyond Cuprates

abstract: Resonating valence bond (RVB) theory of high Tc superconductivity, an electron correlation based mechanism, began as an insightful response by Anderson, to Bednorz and Muller's discovery of high Tc superconductivity in cuprates in late 1986. Shortly a theoretical framework for quantum spin liquids and superconductivity was developed. This theory adresses a formidable strong coupling quantum manybody problem, in modern times. It is built on certain key experimental facts: i) survival of a dynamical Mott localization in a metallic state, ii) proliferation of bond singlets and iii) absence of fermi liquid quasi particles. After summarising RVB theory I will provide an aerial view of, mostly, new superconductors where I believe that, to a large degree RVB mechanism is at work and indicate prospects for even higher Tc's.


title: Correlation function diagnostics for type-I fracton phases

abstract: Fracton phases are recent entrants to the roster of topological phases in three dimensions. They are characterized by subextensively divergent topological degeneracy and excitations that are constrained to move along lower dimensional subspaces, including the eponymous fractons that are immobile in isolation. We develop correlation function diagnostics to characterize Type I fracton phases which build on their exhibiting {\it partial deconfinement}. These are inspired by similar diagnostics from standard gauge theories and utilize a generalized gauging procedure that links fracton phases to classical Ising models with subsystem symmetries. En route, we explicitly construct the spacetime partition function for the plaquette Ising model which, under such gauging, maps into the X-cube fracton topological phase. We numerically verify our results for this model via Monte Carlo calculations.

categories: cond-mat.str-el cond-mat.stat-mech hep-lat quant-ph


title: A Simple and Efficient MapReduce Algorithm for Data Cube Materialization

abstract: Data cube materialization is a classical database operator introduced in Gray et al.~(Data Mining and Knowledge Discovery, Vol.~1), which is critical for many analysis tasks. Nandi et al.~(Transactions on Knowledge and Data Engineering, Vol.~6) first studied cube materialization for large scale datasets using the MapReduce framework, and proposed a sophisticated modification of a simple broadcast algorithm to handle a dataset with a 216GB cube size within 25 minutes with 2k machines in 2012. We take a different approach, and propose a simple MapReduce algorithm which (1) minimizes the total number of copy-add operations, (2) leverages locality of computation, and (3) balances work evenly across machines. As a result, the algorithm shows excellent performance, and materialized a real dataset with a cube size of 35.0G tuples and 1.75T bytes in 54 minutes, with 0.4k machines in 2014.


title: Quantum limits on the time-bandwidth product of an optical resonator

abstract: A thought-provoking proposal by Tsakmakidis et al. [Science 356, 1260 (2017)] suggests that nonreciprocal optics can break a time-bandwidth limit to passive resonators. Here I quantize their resonator model and show that quantum mechanics does impose a limit, or requires extra noise to be added in the same fashion as amplified spontaneous emission in an active resonator. I also use thermodynamics to argue that extra dissipation or noise must be present in their proposed device.


title: Simulating realistically complex comparative effectiveness studies with time-varying covariates and right-censored outcomes

abstract: Simulation studies are useful for evaluating and developing statistical methods for the analyses of complex problems. Performance of methods may be affected by multiple complexities present in real scenarios. Generating sufficiently realistic data for this purpose, however, can be challenging. Our study of the comparative effectiveness of HIV protocols on the risk of cardiovascular disease -- involving the longitudinal assessment of HIV patients -- is such an example. The correlation structure across covariates and within subjects over time must be considered as well as right-censoring of the outcome of interest, time to myocardial infarction. A challenge in simulating the covariates is to incorporate a joint distribution for variables of mixed type -- continuous, binary or polytomous. An additional challenge is incorporating within-subject correlation where some variables may vary over time and others may remain static. To address these issues, we extend the work of Demirtas and Doganay (2012). Identifying a model from which to simulate the right-censored outcome as a function of these covariates builds on work developed by Sylvestre and Abrahamowicz (2007). In this paper, we describe a cohesive and user-friendly approach accompanied by R code to simulate comparative effectiveness studies with right-censored outcomes that are functions of time-varying covariates.


title: Tight Conditional Lower Bounds for Longest Common Increasing Subsequence

abstract: We consider the canonical generalization of the well-studied Longest Increasing Subsequence problem to multiple sequences, called $k$-LCIS: Given $k$ integer sequences $X_1,\dots,X_k$ of length at most $n$, the task is to determine the length of the longest common subsequence of $X_1,\dots,X_k$ that is also strictly increasing. Especially for the case of $k=2$ (called LCIS for short), several algorithms have been proposed that require quadratic time in the worst case. Assuming the Strong Exponential Time Hypothesis (SETH), we prove a tight lower bound, specifically, that no algorithm solves LCIS in (strongly) subquadratic time. Interestingly, the proof makes no use of normalization tricks common to hardness proofs for similar problems such as LCS. We further strengthen this lower bound (1) to rule out $O((nL)^{1-\varepsilon})$ time algorithms for LCIS, where $L$ denotes the solution size, (2) to rule out $O(n^{k-\varepsilon})$ time algorithms for $k$-LCIS, and (3) to follow already from weaker variants of SETH. We obtain the same conditional lower bounds for the related Longest Common Weakly Increasing Subsequence problem.


title: Manyfield Inflation in Random Potentials

abstract: We construct models of inflation with many randomly interacting fields and use these to study the generation of cosmological observables. We model the potentials as multi-dimensional Gaussian random fields (GRFs) and identify powerful algebraic simplifications that, for the first time, make it possible to access the manyfield limit of inflation in GRF potentials. Focussing on small-field, slow-roll, approximate saddle-point inflation in potentials with structure on sub-Planckian scales, we construct explicit examples involving up to 100 fields and generate statistical ensembles comprising of 164,000 models involving 5 to 50 fields. For the subset of these that support at least sixty e-folds of inflation, we use the 'transport method' and $\delta N$ formalism to determine the predictions for cosmological observables at the end of inflation, including the power spectrum and the local non-Gaussianities of the primordial perturbations. We find three key results: i) Planck compatibility is not rare, but future experiments may rule out this class of models; ii) In the manyfield limit, the predictions from these models agree well with, but are sharper than, previous results derived using potentials constructed through non-equilibrium Random Matrix Theory; iii) Despite substantial multifield effects, non-Gaussianities are typically very small: $f_{\rm nl}^{\rm loc} \ll 1$. We conclude that many of the 'generic predictions' of single-field inflation can be emergent features of complex inflation models.


title: Thread-Modular Static Analysis for Relaxed Memory Models

abstract: We propose a memory-model-aware static program analysis method for accurately analyzing the behavior of concurrent software running on processors with weak consistency models such as x86-TSO, SPARC-PSO, and SPARC-RMO. At the center of our method is a unified framework for deciding the feasibility of inter-thread interferences to avoid propagating spurious data flows during static analysis and thus boost the performance of the static analyzer. We formulate the checking of interference feasibility as a set of Datalog rules which are both efficiently solvable and general enough to capture a range of hardware-level memory models. Compared to existing techniques, our method can significantly reduce the number of bogus alarms as well as unsound proofs. We implemented the method and evaluated it on a large set of multithreaded C programs. Our experiments showthe method significantly outperforms state-of-the-art techniques in terms of accuracy with only moderate run-time overhead.


title: Modular Verification of Interrupt-Driven Software

abstract: Interrupts have been widely used in safety-critical computer systems to handle outside stimuli and interact with the hardware, but reasoning about interrupt-driven software remains a difficult task. Although a number of static verification techniques have been proposed for interrupt-driven software, they often rely on constructing a monolithic verification model. Furthermore, they do not precisely capture the complete execution semantics of interrupts such as nested invocations of interrupt handlers. To overcome these limitations, we propose an abstract interpretation framework for static verification of interrupt-driven software that first analyzes each interrupt handler in isolation as if it were a sequential program, and then propagates the result to other interrupt handlers. This iterative process continues until results from all interrupt handlers reach a fixed point. Since our method never constructs the global model, it avoids the up-front blowup in model construction that hampers existing, non-modular, verification techniques. We have evaluated our method on 35 interrupt-driven applications with a total of 22,541 lines of code. Our results show the method is able to quickly and more accurately analyze the behavior of interrupts.


title: A Comet Active Beyond the Crystallization Zone

abstract: We present observations showing in-bound long-period comet C/2017 K2 (PANSTARRS) to be active at record heliocentric distance. Nucleus temperatures are too low (60 K to 70 K) either for water ice to sublimate or for amorphous ice to crystallize, requiring another source for the observed activity. Using the Hubble Space Telescope we find a sharply-bounded, circularly symmetric dust coma 10$^5$ km in radius, with a total scattering cross section of $\sim$10$^5$ km$^2$. The coma has a logarithmic surface brightness gradient -1 over much of its surface, indicating sustained, steady-state dust production. A lack of clear evidence for the action of solar radiation pressure suggests that the dust particles are large, with a mean size $\gtrsim$ 0.1 mm. Using a coma convolution model, we find a limit to the apparent magnitude of the nucleus $V >$ 25.2 (absolute magnitude $H >$ 12.9). With assumed geometric albedo $p_V$ = 0.04, the limit to the nucleus circular equivalent radius is $<$ 9 km. Pre-discovery observations from 2013 show that the comet was also active at 23.7 AU heliocentric distance. While neither water ice sublimation nor exothermic crystallization can account for the observed distant activity, the measured properties are consistent with activity driven by sublimating supervolatile ices such as CO$_2$, CO, O$_2$ and N$_2$. Survival of supervolatiles at the nucleus surface is likely a result of the comet's recent arrival from the frigid Oort cloud.


title: Magnetic monopoles and nonassociative deformations of quantum theory

abstract: We examine certain nonassociative deformations of quantum mechanics and gravity in three dimensions related to the dynamics of electrons in uniform distributions of magnetic charge. We describe a quantitative framework for nonassociative quantum mechanics in this setting, which exhibits new effects compared to ordinary quantum mechanics with sourceless magnetic fields, and the extent to which these theoretical consequences may be experimentally testable. We relate this theory to noncommutative Jordanian quantum mechanics, and show that its underlying algebra can be obtained as a contraction of the alternative algebra of octonions. The uncontracted octonion algebra conjecturally describes a nonassociative deformation of three-dimensional quantum gravity induced by magnetic monopoles, which we propose is realised by a non-geometric Kaluza-Klein monopole background in M-theory.


title: The Stable Rank of Diagonally Constructed ASH Algebras

abstract: We introduce a class of recursive subhomogeneous algebras that we call diagonal subhomogeneous and we give a notion of diagonal maps between these algebras. We show that any simple limit of diagonal subhomogeneous algebras with diagonal maps has stable rank one. As an application we show that for any minimal homeomorphism of a compact Hausdorff space the associated crossed product has stable rank one.


title: Towards Optimally Decentralized Multi-Robot Collision Avoidance via Deep Reinforcement Learning

abstract: Developing a safe and efficient collision avoidance policy for multiple robots is challenging in the decentralized scenarios where each robot generate its paths without observing other robots' states and intents. While other distributed multi-robot collision avoidance systems exist, they often require extracting agent-level features to plan a local collision-free action, which can be computationally prohibitive and not robust. More importantly, in practice the performance of these methods are much lower than their centralized counterparts. We present a decentralized sensor-level collision avoidance policy for multi-robot systems, which directly maps raw sensor measurements to an agent's steering commands in terms of movement velocity. As a first step toward reducing the performance gap between decentralized and centralized methods, we present a multi-scenario multi-stage training framework to find an optimal policy which is trained over a large number of robots on rich, complex environments simultaneously using a policy gradient based reinforcement learning algorithm. We validate the learned sensor-level collision avoidance policy in a variety of simulated scenarios with thorough performance evaluations and show that the final learned policy is able to find time efficient, collision-free paths for a large-scale robot system. We also demonstrate that the learned policy can be well generalized to new scenarios that do not appear in the entire training period, including navigating a heterogeneous group of robots and a large-scale scenario with 100 robots. Videos are available at https://sites.google.com/view/drlmaca


title: Platooning in the Presence of a Speed Drop: A Generalized Control Model

abstract: The positive impacts of platooning on travel time reliability, congestion, emissions, and energy consumption have been shown for homogeneous roadway segments. However, speed limit changes frequently throughout the transportation network, due to either safety-related considerations (e.g., workzone operations) or congestion management schemes (e.g., speed harmonization systems). These abrupt changes in speed limit can result in shock- wave formation and cause travel time unreliability. Therefore, designing a platooning strategy for tracking a reference velocity profile is critical to enabling end-to-end platooning. Accordingly, this study introduces a generalized control model to track a desired velocity profile, while ensuring safety in the platoon of autonomous vehicles. We define appropriate natural error terms and the target curve in the state space of the control system, which is the set of points where all error terms vanish and corresponds to the case when all vehicles move with the desired velocities and in the minimum safe distance between them. In this way, we change the tracking velocity profile problem into a state- feedback stabilization problem with respect to the target curve. Under certain mild assumptions on the Lipschitz constant of the speed drop profile, we show that the stabilizing feedback can be obtained via introducing a natural dynamics for the maximum of the error terms for each vehicle. Moreover, we show that with this stabilizing feedback collisions will not occur if the initial state of the system of vehicles is sufficiently close to the target curve. We also show that the error terms remain bounded throughout the time and space. Two scenarios were simulated, with and without initial perturbations, and results confirmed the effectiveness of the proposed control model in tracking the speed drop while ensuring safety and string stability.


title: Spherical harmonics for multiparticle final states

abstract: We examine collider physics via the spectral power H_l, which utilizes much more detector information than sequential jet reconstruction. We use H_l to define a novel jet definition and use it to reconstruct both e+e- > qqg and e+e > tt final states. We find that a significant amount of pileup can be trivially subtracted from H_l jets, provided that the "shape" of pileup can be determined ahead of time. Finally, we find that future progress with H_l jets will require the addition of degrees of freedom accounting for jet "shape" (e.g. radial energy distribution). This should allow the measurement of jet shapes using much more information than traditional techniques.


title: Back-reaction on background fields: a coherent state approach

abstract: There are many situations in which a strong electromagnetic field may be approximated as a fixed background. Going beyond this approximation, i.e. accounting for the back-reaction of quantum process on the field, is however challenging. Here we develop an approach to this problem which is a straightforward extension of background field methods. The approach follows from the observation that scattering in an on-shell background is equivalent to scattering between coherent states; we show that by deforming these states one can model back-reaction. Focussing on intense laser-matter interactions, we provide examples which model beam depletion and, furthermore, introduce an extremisation principle with which to determine the level of depletion in a given scattering process.


title: The first 62 AGN observed with SDSS-IV MaNGA - I: their characterization and definition of a control sample

abstract: We report the characterization of the first $62$ MaNGA Active Galactic Nuclei (AGN) hosts in the Fifth Product Launch (MPL-5) and the definition of a control sample of non-active galaxies. This control sample - comprising two galaxies for each AGN - was selected in order to match the AGN hosts in terms of stellar mass, redshift, visual morphology and inclination. The stellar masses are in the range $9.4<\log(M/M_\odot)<11.5$, and most objects have redshifts $\leq 0.08$. The AGN sample is mostly comprised of low-luminosity AGN, with only 17 nuclei with $L([OIII]\lambda 5007)\geq 3.8\times 10^{40}$ erg s$^{-1}$ (that we call "strong AGN"). The stellar population of the control sample galaxies within the inner $1$-$3$ kpc is dominated by the old ($\sim$ $4$ - $13$ Gyr) age component, with a small contribution of intermediate age ($\sim 640$-$940$ Myr) and young stars ($\leq 40$ Myr) to the total light at $5700\AA$. While the weaker AGN show a similar age distribution to that of the control galaxies, the strong AGN show an increased contribution of younger stars and a decreased contribution of older stars. Examining the relationship between the AGN stellar population properties and $L([OIII])$, we find that with increasing $L([OIII])$, the AGN exhibit a decreasing contribution from the oldest ($>4$ Gyr) stellar population relative to control galaxies, but have an increasing contribution from the younger components with ages $\sim 40$ Myr. We also find a correlation of the mean age differences (AGN - control) with $L([OIII])$, in the sense that more luminous AGN are younger than the control objects, while the low-luminosity AGN seem to be older. These results support a connection between the growth of the galaxy bulge via formation of new stars and the growth of the Supermassive Black Hole via matter accretion in the AGN phase.


title: Learning Complex Dexterous Manipulation with Deep Reinforcement Learning and Demonstrations

abstract: Dexterous multi-fingered hands are extremely versatile and provide a generic way to perform a multitude of tasks in human-centric environments. However, effectively controlling them remains challenging due to their high dimensionality and large number of potential contacts. Deep reinforcement learning (DRL) provides a model-agnostic approach to control complex dynamical systems, but has not been shown to scale to high-dimensional dexterous manipulation. Furthermore, deployment of DRL on physical systems remains challenging due to sample inefficiency. Consequently, the success of DRL in robotics has thus far been limited to simpler manipulators and tasks. In this work, we show that model-free DRL can effectively scale up to complex manipulation tasks with a high-dimensional 24-DoF hand, and solve them from scratch in simulated experiments. Furthermore, with the use of a small number of human demonstrations, the sample complexity can be significantly reduced, which enables learning with sample sizes equivalent to a few hours of robot experience. The use of demonstrations result in policies that exhibit very natural movements and, surprisingly, are also substantially more robust.


title: Parametric Resonance in Quantum Electrodynamics Vacuum Birefringence

abstract: Vacuum magnetic birefringence is one of the most interesting non-linear phenomena in quantum electrodynamics because it is a pure photon-photon result of the theory and it directly signalizes the violation of the classical superposition principle of electromagnetic fields in the full quantum theory. We perform analytical and numerical calculations when an electromagnetic wave interacts with an oscillating external magnetic field. We find that in an ideal cavity, when the external field frequency is around the electromagnetic wave frequency, the normal and parallel components of the wave suffer parametric resonance at different rates, producing a vacuum birefringence effect growing in time. We also study the case where there is no cavity and the oscillating magnetic field is spatially localized in a region of length $L$. In both cases we find also a rotation of the elliptical axis.


title: Overcoming Exploration in Reinforcement Learning with Demonstrations

abstract: Exploration in environments with sparse rewards has been a persistent problem in reinforcement learning (RL). Many tasks are natural to specify with a sparse reward, and manually shaping a reward function can result in suboptimal performance. However, finding a non-zero reward is exponentially more difficult with increasing task horizon or action dimensionality. This puts many real-world tasks out of practical reach of RL methods. In this work, we use demonstrations to overcome the exploration problem and successfully learn to perform long-horizon, multi-step robotics tasks with continuous control such as stacking blocks with a robot arm. Our method, which builds on top of Deep Deterministic Policy Gradients and Hindsight Experience Replay, provides an order of magnitude of speedup over RL on simulated robotics tasks. It is simple to implement and makes only the additional assumption that we can collect a small set of demonstrations. Furthermore, our method is able to solve tasks not solvable by either RL or behavior cloning alone, and often ends up outperforming the demonstrator policy.


title: Blandford-Znajek process in vacuo and its holographic dual

abstract: Blandford and Znajek discovered a process by which a spinning black hole can transfer rotational energy to a plasma, offering a mechanism for energy and jet emissions from quasars. Here we describe a version of this mechanism that operates with only vacuum electromagnetic fields outside the black hole. The setting, which is not astrophysically realistic, involves either a cylindrical black hole or one that lives in 2+1 spacetime dimensions, and the field is given in simple, closed form for a wide class of metrics. For asymptotically Anti-de Sitter black holes in 2+1 dimensions the holographic dual of this mechanism is the transfer of angular momentum and energy, via a resistive coupling, from a rotating thermal state containing an electric field to an external charge density rotating more slowly than the thermal state. In particular, the entropy increase of the thermal state due to Joule heating matches the Bekenstein-Hawking entropy increase of the black hole.


title: Initiating and monitoring the evolution of single electrons within atom-defined structures

abstract: Using a non-contact atomic force microscope we track and manipulate the position of single electrons confined to atomic structures engineered from silicon dangling bonds (DBs) on the hydrogen terminated silicon surface. By varying the probe-sample separation we mechanically manipulate the equilibrium position of individual surface silicon atoms and use this to directly switch the charge state of individual DBs. Because this mechanism is based on short range interactions and can be performed without applied bias voltage, we maintain both site-specific selectivity and single-electron control. We extract the short range forces involved with this mechanism by subtracting the long range forces acquired on a dimer vacancy site. As a result of relaxation of the silicon lattice to accommodate negatively charged DBs we observe charge configurations of DB structures that remain stable for many seconds at 4.5 K. Subsequently we use charge manipulation to directly prepare the ground state and metastable charge configurations of DB structures composed of up to six atoms.


title: ABJM $\theta$-Bremsstrahlung at four loops and beyond: non-planar corrections

abstract: We consider the Bremsstrahlung function associated to a 1/6-BPS Wilson loop in ABJM theory, with a cusp in the couplings to scalar fields. We non-trivially extend its recent four-loop computation at weak coupling to include non-planar corrections. We have recently proposed a conjecture relating this object to supersymmetric circular Wilson loops with multiple windings, which can be computed via localization. We find agreement between this proposal and the perturbative computation of the Bremsstrahlung function, including color sub-leading corrections. This supports the conjecture and hints at its validity beyond the planar approximation.


title: A construction of quarter BPS coherent states and Brauer algebras

abstract: BPS coherent states closely resemble semiclassical states and they have gravity dual descriptions in terms of semiclassical geometries. The half BPS coherent states have been well studied, however less is known about quarter BPS coherent states. Here we provide a construction of quarter BPS coherent states. They are coherent states built with two matrix fields, generalizing the half BPS case. These states are both the eigenstates of annihilation operators and in the kernel of dilatation operator. Another useful labeling of quarter BPS states is by representations of Brauer algebras and their projection onto a subalgebra $\mathbb{C}[S_n\times S_m]$. Here, the Schur-Weyl duality for the Walled Brauer algebra plays an important role in organizing the operators. One interesting subclass of these Brauer states are labeled by representations involving two Young tableaux. We obtain the overlap between quarter BPS Brauer states and quarter BPS coherent states, where the Schur polynomials are used. We also derive superposition formulas transforming a truncated version of quarter BPS coherent states to quarter BPS Brauer states. The entanglement entropy of Brauer states as well as the overlap between Brauer states and squeezed states are also computed.

categories: hep-th gr-qc math-ph math.MP math.RT quant-ph


title: A simple anisotropic three-dimensional quantum spin liquid with fracton topological order

abstract: We present a three-dimensional cubic lattice spin model, anisotropic in the $\hat{z}$ direction, that exhibits fracton topological order. The latter is a novel type of topological order characterized by the presence of immobile pointlike excitations, named fractons, residing at the corners of an operator with two-dimensional support. As other recent fracton models, ours exhibits a subextensive ground state degeneracy: On an $L_x\times L_y\times L_z$ three-torus, it has a $2^{2L_z}$ topological degeneracy, and an additional non-topological degeneracy equal to $2^{L_xL_y-2}$. The fractons can be combined into composite excitations that move either in a straight line along the $\hat{z}$ direction, or freely in the $xy$ plane at a given height $z$. While our model draws inspiration from the toric code, we demonstrate that it cannot be adiabatically connected to a layered toric code construction. Additionally, we investigate the effects of imposing open boundary conditions on our system. We find zero energy modes on the surfaces perpendicular to either the $\hat{x}$ or $\hat{y}$ directions, and their absence on the surfaces normal to $\hat{z}$. This result can be explained using the properties of the two kinds of composite two-fracton mobile excitations.


title: Divergences in open quantum systems

abstract: We show that for cubic scalar field theories in five and more spacetime dimensions, and for the $T = 0$ limit of the Caldeira-Leggett model, the quantum master equation for long-wavelength modes initially unentangled from short-distance modes, and at second order in perturbation theory, contains divergences in the non-Hamiltonian terms. These divergences ensure that the equations of motion for expectation values of composite operators closes on expectation values of renormalized operators. Along the way we show that initial "jolt" singularities which occur in the equations of motion for operators linear in the fundamental variables persist for quadratic operators, and are removed if one chooses an initial state projected onto low energies, following the Born-Oppenheimer approximation.


title: Universal broadening of the light cone in low-temperature transport

abstract: We consider the low-temperature transport properties of critical one-dimensional systems which can be described, at equilibrium, by a Luttinger liquid. We focus on the prototypical setting where two semi-infinite chains are prepared in two thermal states at small but different temperatures and suddenly joined together. At large distances $x$ and times $t$, conformal field theory characterizes the energy transport in terms of a single light cone spreading at the sound velocity $v$. Energy density and current take different constant values inside the light cone, on its left, and on its right, resulting in a three-step form of the corresponding profiles as a function of $\zeta=x/t$. Here, using a non-linear Luttinger liquid description, we show that for generic observables this picture is spoiled as soon as a non-linearity in the spectrum is present. In correspondence of the transition points $x/t=\pm v$ a novel universal region emerges at infinite times, whose width is proportional to the temperatures on the two sides. In this region, expectation values have a different temperature dependence and show smooth peaks as a function of $\zeta$. We explicitly compute the universal function describing such peaks. In the specific case of interacting integrable models, our predictions are analytically recovered by the generalized hydrodynamic approach.


title: Modeling baryonic interactions with the Clausius-type equation of state

abstract: The quantum statistical Clausius-based equation of state is used to describe the system of interacting nucleons. The interaction parameters $a$, $b$, and $c$ of the model are fixed by the empirically known nuclear ground state properties and nuclear incompressibility modulus. The model is generalized to describe the baryon-baryon interactions in the hadron resonance gas (HRG). The predictions of such a Clausius-HRG model are confronted with the lattice QCD data at zero and at small chemical potentials, and are also contrasted with the standard van der Waals approach. It is found that the behavior of the lattice QCD observables in a high-temperature hadron gas is sensitive to the nuclear matter properties. An improved description of the nuclear incompressibility factor correlates with an improved description of the lattice QCD data in the crossover transition region.


title: Vafa-Witten theory and iterated integrals of modular forms

abstract: Vafa-Witten (VW) theory is a topologically twisted version of N=4 supersymmetric Yang-Mills theory. S-duality suggests that the partition function of VW theory with gauge group SU(N) transforms as a modular form under duality transformations. Interestingly, Vafa and Witten demonstrated the presence of a modular anomaly, when the theory has gauge group SU(2) and is considered on the complex projective plane P2. This modular anomaly could be expressed as an integral of a modular form, and also be traded for a holomorphic anomaly. We demonstrate that the modular anomaly for gauge group SU(3) involves an iterated integral of modular forms. Moreover, the modular anomaly for SU(3) can be traded for a holomorphic anomaly, which is shown to factor into a product of the partition functions for lower rank gauge groups. The SU(3) partition function is mathematically an example of a mock modular form of depth two.