title: An integrable Lorentz-breaking deformation of two-dimensional CFTs

abstract: It has been recently shown that the deformation of an arbitrary two-dimensional conformal field theory by the composite irrelevant operator $T \bar T$, built from the components of the stress tensor, is solvable; in particular, the finite-size spectrum of the deformed theory can be obtained from that of the original CFT through a universal formula. We study a similarly universal, Lorentz-breaking deformation of two-dimensional CFTs that posess a conserved $U(1)$ current, $J$. The deformation takes the schematic form $J \bar T$ and is interesting because it preserves an $SL(2,\mathbb{R}) \times U(1)$ subgroup of the original global conformal symmetries. For the case of a purely (anti)chiral current, we find the finite-size spectrum of the deformed theory and study its thermodynamic properties. We test our predictions in a simple example involving deformed free fermions.


title: High Energy Resummation of Jet Observables

abstract: In this paper we investigate the extension of high energy resummation at LLx accuracy to jet observables. In particular, we present the high energy resummed expression of the transverse momentum distribution of the outgoing parton in the general partonic process $g(q) + g(q) \to g(q) + X$. In order to reach this result, several new ideas are introduced and exploited. First we prove that LLx resummation is achieved by dressing with hard radiation an off-shell gluon initiated LO process even if its on-shell limit is vanishing or trivial. Then we present a gauge-invariant framework where these calculations can be performed by using the modern helicity techniques. Finally, we show a possible way to restore gluon indistinguishability in the final state, which is otherwise lost in the resummation procedure, at all orders in $\alpha_s$ at LLx. All partonic channels are then resummed and cross-checked against fixed-order calculations up to $\mathcal{O}(\alpha_s^3)$


title: The dependence of convective core overshooting on stellar mass: semi-empirical determination using the diffusive approach with two different element mixtures

abstract: Convective core overshooting has a strong influence on the evolution of stars of moderate and high mass. Studies of double-lined eclipsing binaries and stellar oscillations have renewed interest in the possible dependence of overshooting on stellar mass, which has been poorly constrained by observations so far. Here we have used a sample of 29 well-studied double-lined eclipsing binaries in key locations of the H-R diagram to establish the explicit dependence of f(ov) on mass, where f(ov) is the free parameter in the diffusive approximation to overshooting. Measurements of the masses, radii, and temperatures of the binary components were compared against stellar evolution calculations based on the MESA code to infer semi-empirical values of f(ov) for each component. We find a clear mass dependence such that f(ov) rises sharply from zero in the range 1.2--2.0 solar masses, and levels off thereafter up to the 4.4 solar mass limit of our sample. Tests with two different element mixtures indicate the trend is the same, and we find it is also qualitatively similar to the one established in our previous study with the classical step-function implementation of overshooting characterized by the free parameter alpha(ov). Based on these measurements we infer an approximate relationship between the two overshooting parameters of alpha(ov)/f(ov) = 11.36 +/- 0.22, with a possible dependence on stellar properties.


title: Spheres, Charges, Instantons, and Bootstrap: A Five-Dimensional Odyssey

abstract: We combine supersymmetric localization and the conformal bootstrap to study five-dimensional superconformal field theories. To begin, we classify the admissible counter-terms and derive a general relation between the five-sphere partition function and the conformal and flavor central charges. Along the way, we discover a new superconformal anomaly in five dimensions. We then propose a precise triple factorization formula for the five-sphere partition function, that incorporates instantons and is consistent with flavor symmetry enhancement. We numerically evaluate the central charges for the rank-one Seiberg and Morrison-Seiberg theories, and find strong evidence for their saturation of bootstrap bounds, thereby determining the spectra of long multiplets in these theories. Lastly, our results provide new evidence for the $F$-theorem and possibly a $C$-theorem in five-dimensional superconformal theories.


title: Quantum Mechanics from Ergodic Average of Microstates

abstract: We formulate quantum mechanics as an effective theory of an underlying structure characterized by microstates $|{\mathcal M}^j(t)\rangle$, each one defined by the quantum state $|\Psi(t)\rangle$ and a complete set of commutative observables $O^j$. At any time $t$, $|{\mathcal M}^j(t)\rangle$ corresponds to a state $|O^j_k\rangle$, for some $k$ depending on $t$, and jumps after time intervals whose duration, of the order of the Compton time $\tau$, is proportional to the probability $|\langle O_k^j|\Psi(t)\rangle|^2$. This reproduces the Born rule and mimics the wave-particle duality. The theory is based on a partition of time whose flow is characterized by quantum probabilities. Ergodicity arises at ordinary quantum scales with the expectation values corresponding to time averaging over a period $\tau$. The measurement of $O^j$ provides a new partition of time and the outcome is the state $|O_k^j\rangle$ to which $|{\mathcal M}^j(t)\rangle$ corresponds at that time. The formulation, that shares some features with the path integral, can be tested by experiments involving time intervals of order $\tau$.

categories: quant-ph cond-mat.stat-mech hep-ph hep-th


title: A census of cool core galaxy clusters in IllustrisTNG

abstract: The thermodynamic structure of hot gas in galaxy clusters is sensitive to astrophysical processes and typically difficult to model with galaxy formation simulations. We explore the fraction of cool-core (CC) clusters in a large sample of $370$ clusters from IllustrisTNG, examining six common CC definitions. IllustrisTNG produces continuous CC criteria distributions, the extremes of which are classified as CC and non-cool-core (NCC), and the criteria are increasingly correlated for more massive clusters. At $z=0$, the CC fractions for $2$ criteria are in reasonable agreement with the observed fractions but the other $4$ CC fractions are lower than observed. This result is partly driven by systematic differences between the simulated and observed gas fraction profiles. The simulated CC fractions with redshift show tentative agreement with the observed fractions, but linear fits demonstrate that the simulated evolution is steeper than observed. The conversion of CCs to NCCs appears to begin later and act more rapidly in the simulations. Examining the fraction of CCs and NCCs defined as relaxed we find no evidence that CCs are more relaxed, suggesting that mergers are not solely responsible for disrupting CCs. A comparison of the median thermodynamic profiles defined by different CC criteria shows that the extent to which they evolve in the cluster core is dependent on the CC criteria. We conclude that the thermodynamic structure of galaxy clusters in IllustrisTNG shares many similarities with observations, but achieving better agreement most likely requires modifications of the underlying galaxy formation model.


title: The Mass Growth and Stellar Ages of Galaxies: Observations versus Simulations

abstract: Using observed stellar mass functions out to $z=5$, we measure the main progenitor stellar mass growth of descendant galaxies with masses of $\log{M_{*}/M_{\odot}}=11.5,11.0,10.5,10.0$ at $z\sim0.1$ using an evolving cumulative number density selection. From these mass growth histories, we are able to measure the time at which half the total stellar mass of the descendant galaxy was assembled, $t_{a}$, which, in order of decreasing mass corresponds to redshifts of $z_{a}=1.28, 0.92, 0.60$ and $0.51$. We compare this to the median light-weighted stellar age $t_{*}$ ($z_{*} = 2.08, 1.49, 0.82$ and $0.37$) of a sample of low redshift SDSS galaxies (from the literature) and find the timescales are consistent with more massive galaxies forming a higher fraction of their stars ex-situ compared to lower mass descendants. We find that both $t_{*}$ and $t_{a}$ strongly correlate with mass which is in contrast to what is found in the EAGLE hydrodynamical simulation which shows a flat relationship between $t_{a}$ and $M_{*}$. However, the semi-analytic model of \citet{henriques2015} is consistent with the observations in both $t_{a}$ and $t_{*}$ with $M_{*}$, showing the most recent semi-analytic models are better able to decouple the evolution of the baryons from the dark matter in lower-mass galaxies.


title: On the nature and physical conditions of the luminous Lya emitter CR7 and its rest-frame UV components

abstract: We present new HST/WFC3 observations and re-analyse VLT data to unveil the continuum, variability and rest-frame UV lines of the multiple UV clumps of the most luminous Ly$\alpha$ emitter at z=6.6, CR7. Our re-reduced, flux calibrated X-SHOOTER spectra of CR7 reveal a HeII emission line in observations obtained along the major axis of Lyman-alpha (Lya) emission with the best seeing conditions. HeII is spatially offset by +0.8'' from the peak of Lya emission, and it is found towards clump B. Our WFC3 grism spectra detects the UV continuum of CR7's clump A, yielding a power law with $\beta=-2.5^{+0.6}_{-0.7}$ and $M_{UV}=-21.87^{+0.25}_{-0.20}$. No significant variability is found for any of the UV clumps on their own, but there is tentative (~2.2$\sigma$) brightening of CR7 in F110W as a whole from 2012 to 2017. HST grism data fail to robustly detect rest-frame UV lines in any of the clumps, implying fluxes <2x10$^{-17}$ erg s$^{-1}$ cm$^{-2}$ (3 $\sigma$). We perform CLOUDY modelling to constrain the metallicity and the ionising nature of CR7. CR7 seems to be actively forming stars without any clear AGN activity in clump A, consistent with a metallicity of ~0.05-0.2 Z$_{\odot}$. Component C or an inter-clump component between B and C may host a high ionisation source. Our results highlight the need for spatially resolved information to study the formation and assembly of early galaxies.


title: Two-Face(s): ionized and neutral gas winds in the local Universe

abstract: We present a comprehensive study of the Na I $\lambda$5890, 5895 (Na I D) resonant lines in the Sloan Digital Sky Survey (SDSS, DR7) spectroscopic sample to look for neutral gas outflows in the local galaxies. Individual galaxy spectra are stacked in bins of M${\star}$ and SFR to investigate the dependence of galactic wind occurrence and velocity as a function of the galaxy position in the SFR-$M{\star}$ plane. In massive galaxies at the high SFR tail we find evidence of a significant blue-shifted Na I D absorption, which we interpret as evidence of neutral outflowing gas. The occurrence of the blue-shifted absorption is observed at the same significance for purely SF galaxies, AGN and composite systems at fixed SFR. In all classes of objects the blue-shift is the largest and the Na I D equivalent width the smallest for face-on galaxies while the absorption feature is at the systemic velocity for edge-on systems. This indicates that the neutral outflow is mostly perpendicular or biconical with respect to the galactic disk. We also compare the kinematics of the neutral gas with the ionized gas phase as traced by the [OIII]$\lambda$5007, H$\alpha$, [NII]$\lambda6548$ and [NII]$\lambda6584$ emission lines. Differently for the neutral gas phase, all the emission lines show evidence of perturbed kinematics only in galaxies with a significant level of nuclear activity and, they are independent from the disk inclination. In conclusion, we find that, in the local Universe, galactic winds show two faces which are related to two different ejection mechanisms, namely the neutral outflowing gas phase related to the SF activity along the galaxy disk and the ionized phase related to the AGN feedback. In both the neutral and ionized gas phases, the observed wind velocities suggest that the outflowing gas remains bound to the galaxy with no definitive effect on the gas reservoir.


title: Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions

abstract: We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact.


title: Hydrodynamics of electrons in graphene

abstract: Generic interacting many-body quantum systems are believed to behave as classical fluids on long time and length scales. Due to rapid progress in growing exceptionally pure crystals, we are now able to experimentally observe this collective motion of electrons in solid-state systems, including graphene. We present a review of recent progress in understanding the hydrodynamic limit of electronic motion in graphene, written for physicists from diverse communities. We begin by discussing the "phase diagram" of graphene, and the inevitable presence of impurities and phonons in experimental systems. We derive hydrodynamics, both from a phenomenological perspective and using kinetic theory. We then describe how hydrodynamic electron flow is visible in electronic transport measurements. Although we focus on graphene in this review, the broader framework naturally generalizes to other materials. We assume only basic knowledge of condensed matter physics, and no prior knowledge of hydrodynamics.

categories: cond-mat.str-el cond-mat.mes-hall hep-th physics.flu-dyn


title: Far-infrared to millimeter data of protoplanetary disks: dust growth in the Taurus, Ophiuchus, and Chamaeleon I star-forming regions

abstract: Far-infrared and (sub)millimeter fluxes can be used to study dust in protoplanetary disks, the building blocks of planets. Here, we combine observations from the Herschel Space Observatory with ancillary data of 284 protoplanetary disks in the Taurus, Chamaeleon I, and Ophiuchus star-forming regions, covering from the optical to mm/cm wavelengths. We analyze their spectral indices as a function of wavelength and determine their (sub)millimeter slopes when possible. Most disks display observational evidence of grain growth, in agreement with previous studies. No correlation is found between other tracers of disk evolution and the millimeter spectral indices. A simple disk model is used to fit these sources, and we derive posterior distributions for the optical depth at 1.3 mm and 10 au, the disk temperature at this same radius, and the dust opacity spectral index. We find the fluxes at 70 microns to correlate strongly with disk temperatures at 10 au, as derived from these simple models. We find tentative evidence for spectral indices in Chamaeleon I being steeper than those of disks in Taurus/Ophiuchus, although more millimeter observations are needed to confirm this trend and identify its possible origin. Additionally, we determine the median spectral energy distribution of each region and find them to be similar across the entire wavelength range studied, possibly due to the large scatter in disk properties and morphologies.


title: Clustering in the Stellar Abundances Space

abstract: We have studied the chemical enrichment history of the interstellar medium through an analysis of the n-dimensional stellar abundances space. This work is a non-parametric analysis of the stellar chemical abundance space. The main goal is to study the stars from their organization within this abundance space. Within this space, we seek to find clusters (in a statistical sense), that is, stars likely to share similar chemo-evolutionary history, using two methods: the hierarchical clustering and the principal component analysis. We analysed some selected abundance surveys available in the literature. For each sample, we labelled the group of stars according to its average abundance curve. In all samples, we identify the existence of a main enrichment pattern of the stars, which we call chemical enrichment flow. This flow is set by the structured and well defined mean rate at which the abundances of the interstellar medium increase, resulting from the mixture of the material ejected from the stars and stellar mass loss and interstellar medium gas. One of the main results of our analysis is the identification of subgroups of stars with peculiar chemistry. These stars are situated in regions outside of the enrichment flow in the abundance space. These peculiar stars show a mismatch in the enrichment rate of a few elements, such as Mg, Si, Sc and V, when compared to the mean enrichment rate of the other elements of the same stars. We believe the existence of these groups of stars with peculiar chemistry may be related to the accretion of planetary material onto stellar surfaces or may be due to production of the same chemical element by different nucleosynthetic sites.


title: Fate of superconductivity in three-dimensional disordered Luttinger semimetals

abstract: Superconducting instability can occur in three-dimensional quadratic band crossing semimetals only at a finite coupling strength, due to the vanishing of density of states at the quadratic band touching point. Since realistic materials are always disordered to some extent, we study the effect of short-ranged correlated disorder on this superconducting quantum critical point, using a controlled loop-expansion applying dimensional regularization. The renormalization group (RG) scheme allows us to determine the RG flows of the various interaction strengths, and shows that disorder destroys the superconducting quantum critical point. In fact, the system exhibits a runaway flow to strong disorder.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.supr-con hep-th


title: Probing the structure of the gas in the Milky Way through X-ray high resolution spectroscopy

abstract: We have developed a new X-ray absorption model, called {\tt IONeq}, which computes the optical depth $\tau(E)$ simultaneously for ions of all abundant elements, assuming ionization equilibrium and taking into account turbulent broadening. We use this model to analyze the interstellar medium (ISM) absorption features in the Milky Way for a sample of 18 galactic (LMXBs) and 42 extragalactic sources (mainly Blazars). The absorbing ISM was modeled as a combination of three components/phases - neutral ($T\lesssim1\times10^{4}$ K), warm ($T\sim5\times 10^{4}$ K) and hot ($T\sim2\times10^{6}$ K). We found that the spatial distribution of both, neutral and warm components, are difficult to describe using smooth profiles due to nonuniform distribution of the column densities over the sky. For the hot phase we used a combination of a flattened disk and a halo, finding comparable column densities for both spatial components, in the order of $\sim 6-7\times10^{18}\;{\rm cm^{-2}}$, although this conclusion depends on the adopted parametrization. If the halo component has sub-solar abundance $Z$, then the column density has be scaled up by a factor $\frac{Z_\odot}{Z}$. The vertically integrated column densities of the disk components suggests the following mass fractions for these three ISM phases in the Galactic disk: neutral $\sim~89\%$, warm $\sim 8\%$ and hot $\sim 3\%$ components, respectively. The constraints on the radial distribution of the halo component of the hot component are weak.


title: Light dark Higgs boson in minimal sub-GeV dark matter scenarios

abstract: Minimal scenarios with light (sub-GeV) dark matter whose relic density is obtained from thermal freeze-out must include new light mediators. In particular, a very well-motivated case is that of a new "dark" massive vector gauge boson mediator. The mass term for such mediator is most naturally obtained by a "dark Higgs mechanism" which leads to the presence of an often long-lived dark Higgs boson whose mass scale is the same as that of the mediator. We study the phenomenology and experimental constraints on two minimal, self-consistent dark sectors that include such a light dark Higgs boson. In one the dark matter is a pseudo-Dirac fermion, in the other a complex scalar. We find that the constraints from BBN and CMB are considerably relaxed in the framework of such minimal dark sectors. We present detection prospects for the dark Higgs boson in existing and projected proton beam-dump experiments. We show that future searches at experiments like Xenon1T or LDMX can probe all the relevant parameter space, complementing the various upcoming indirect constraints from astrophysical observations.


title: Probing new charged scalars with neutrino trident production

abstract: We investigate the possibility of using neutrino trident production to probe leptophilic charged scalars at future high intensity neutrino experiments. We show that under specific assumptions, this production process can provide competitive sensitivity for generic charged scalars as compared to common existing bounds. We also investigate how the recently proposed mixed-flavour production - where the two oppositely charged leptons in the final state need not be muon flavoured - can give a 20-50% increase in sensitivity for certain configurations of new physics couplings as compared to traditional trident modes. We then categorize all renormalizable leptophilic scalar extensions based on their representation under $SU(2)\times U(1)$, and discuss the Higgs triplet and Zee-Babu models as explicit UV realizations. We find that the inclusion of additional doubly charged scalars and the need to reproduce neutrino masses make trident production uncompetitive with current bounds for these specific UV completions. Our work represents the first application of neutrino trident production to study charged scalars, and of mixed-flavour final states to study physics beyond the Standard Model.


title: Limits on quantum gravity effects from Swift short gamma-ray bursts

abstract: The delay in the arrival times between high and low energy photons from cosmic sources can be used to test the violation of the Lorentz invariance (LIV), predicted by some quantum gravity theories, and to constrain its characteristic energy scale ${\rm E_{QG}}$ that is of the order of the Planck energy. Gamma-ray bursts (GRBs) and blazars are ideal for this purpose thanks to their broad spectral energy distribution and cosmological distances: at first order approximation, the constraints on ${\rm E_{QG}}$ are proportional to the photon energy separation and the distance of the source. However, the LIV tiny contribution to the total time delay can be dominated by intrinsic delays related to the physics of the sources: long GRBs typically show a delay between high and low energy photons related to their spectral evolution (spectral lag). Short GRBs have null intrinsic spectral lags and are therefore an ideal tool to measure any LIV effect. We considered a sample of $15$ short GRBs with known redshift observed by Swift and we estimate a limit on ${\rm E_{QG}}\gtrsim 1.5\times 10^{16}$ GeV. Our estimate represents an improvement with respect to the limit obtained with a larger (double) sample of long GRBs and is more robust than the estimates on single events because it accounts for the intrinsic delay in a statistical sense.


title: Cloud Atlas: Discovery of Rotational Spectral Modulations in a Low-mass, L-type Brown Dwarf Companion to a Star

abstract: Observations of rotational modulations of brown dwarfs and giant exoplanets allow the characterization of condensate cloud properties. As of now rotational spectral modulations have only been seen in three L-type brown dwarfs. We report here the discovery of rotational spectral modulations in LP261-75B, an L6-type intermediate surface gravity companion to an M4.5 star. As a part of the Cloud Atlas Treasury program we acquired time-resolved Wide Field Camera 3 grism spectroscopy (1.1--1.69~$\mu$m) of LP261-75B. We find gray spectral variations with the relative amplitude displaying only a weak wavelength dependence and no evidence for lower-amplitude modulations in the 1.4~$\mu$m water band than in the adjacent continuum. The likely rotational modulation period is 4.78$\pm$0.95 h, although the rotational phase is not well sampled. The minimum relative amplitude in the white light curve measured over the whole wavelength range is 2.41$\pm$0.14\%. We report an unusual light curve with seemingly three peaks approximately evenly distributed in rotational phase. The spectral modulations suggests that the upper atmosphere cloud properties in {LP261-75B} are similar to two other mid-L dwarfs of typical infrared colors, but differ from that of the extremely red L-dwarf WISE0047.


title: All the noncontextuality inequalities for arbitrary prepare-and-measure experiments with respect to any fixed sets of operational equivalences

abstract: Within the framework of generalized noncontextuality, we introduce a general technique for systematically deriving noncontextuality inequalities for any experiment involving finitely many preparations and finitely many measurements, each of which has a finite number of outcomes. Given any fixed sets of operational equivalences among the preparations and among the measurements as input, the algorithm returns a set of noncontextuality inequalities whose satisfaction is necessary and sufficient for a set of operational data to admit of a noncontextual model. Additionally, we show that the space of noncontextual data tables always defines a polytope. Finally, we provide a computationally efficient means for testing whether any set of numerical data admits of a noncontextual model, with respect to any fixed operational equivalences. Together, these techniques provide complete methods for characterizing arbitrary noncontextuality scenarios, both in theory and in practice. Because a quantum prepare-and-measure experiment admits of a noncontextual model if and only if it admits of a positive quasiprobability representation, our techniques also determine the necessary and sufficient conditions for the existence of such a representation.


title: The Parsec-Scale Relationship Between I_CO and A_V in Local Molecular Clouds

abstract: We measure the parsec-scale relationship between integrated CO intensity (I_CO) and visual extinction (A_V) in 24 local molecular clouds using maps of CO emission and dust optical depth from Planck. This relationship informs our understanding of CO emission across environments, but clean Milky Way measurements remain scarce. We find uniform I_CO for a given A_V, with the results bracketed by previous studies of the Pipe and Perseus clouds. Our measured I_CO-A_V relation broadly agrees with the standard Galactic CO-to-H2 conversion factor, the relation found for the Magellanic clouds at coarser resolution, and numerical simulations by Glover & Clark (2016). This supports the idea that CO emission primarily depends on shielding, which protects molecules from dissociating radiation. Evidence for CO saturation at high A_V and a threshold for CO emission at low A_V varies remains uncertain due to insufficient resolution and ambiguities in background subtraction. Resolution of order 0.1 pc may be required to measure these features. We use this I_CO-AV relation to predict how the CO-to-H2 conversion factor (X_CO) would change if the Solar Neighborhood clouds had different dust-to-gas ratio (metallicity). The calculations highlight the need for improved observations of the CO emission threshold and HI shielding layer depth. They are also sensitive to the shape of the column density distribution. Because local clouds collectively show a self-similar distribution, we predict a shallow metallicity dependence for X_CO down to a few tenths of solar metallicity. However, our calculations also imply dramatic variations in cloud-to-cloud X_CO at subsolar metallicity.


title: HyperMinHash: MinHash in LogLog space

abstract: In this extended abstract, we describe and analyze a lossy compression of MinHash from buckets of size $O(\log n)$ to buckets of size $O(\log\log n)$ by encoding using floating-point notation. This new compressed sketch, which we call HyperMinHash, as we build off a HyperLogLog scaffold, can be used as a drop-in replacement of MinHash. Unlike comparable Jaccard index fingerprinting algorithms in sub-logarithmic space (such as b-bit MinHash), HyperMinHash retains MinHash's features of streaming updates, unions, and cardinality estimation. For a multiplicative approximation error $1+ \epsilon$ on a Jaccard index $ t $, given a random oracle, HyperMinHash needs $O\left(\epsilon^{-2} \left( \log\log n + \log \frac{1}{ t \epsilon} \right)\right)$ space. HyperMinHash allows estimating Jaccard indices of 0.01 for set cardinalities on the order of $10^{19}$ with relative error of around 10\% using 64KiB of memory; MinHash can only estimate Jaccard indices for cardinalities of $10^{10}$ with the same memory consumption.


title: User-centric interdependent urban systems: using time-of-day electricity usage data to predict morning roadway congestion

abstract: Urban systems are interdependent as individuals' daily activities engage using those urban systems at certain time of day and locations. There may exist clear spatial and temporal correlations among usage patterns across all urban systems. This paper explores such a correlation among energy usage and roadway congestion. We propose a general framework to predict congestion starting time and congestion duration in the morning using the time-of-day electricity use data from anonymous households with no personally identifiable information. We show that using time-of-day electricity data from midnight to early morning from 322 households in the City of Austin, can make reliable prediction of congestion starting time of several highway segments, at the time as early as 2am. This predictor significantly outperforms a time-series predictor that uses only real-time travel time data up to 6am. We found that 8 out of the 10 typical electricity use patterns have statistically significant affects on morning congestion on highways in Austin. Some patterns have negative effects, represented by an early spike of electricity use followed by a drastic drop that could imply early departure from home. Others have positive effects, represented by a late night spike of electricity use possible implying late night activities that can lead to late morning departure from home.


title: Direct observations of different sunspot waves influenced by umbral flashes

abstract: We report the simultaneous presence of chromospheric umbral flashes and associated umbral waves, and propagating coronal disturbances, in a sunspot and related active region. We have analyzed time-distance maps obtained using the observations from Atmospheric Imaging Assembly (AIA) on-board Solar Dynamics Observatory (SDO). These maps show the simultaneous occurrence of different sunspot oscillation and waves such as umbral flashes, umbral waves, and coronal waves. Analysis of the original light curves, i.e., without implementing any Fourier filtering on them, show that the amplitudes of different sunspot waves observed at different atmospheric layers change in synchronization with the light curves obtained from the umbral flash region, thus demonstrating that these oscillations are modulated by umbral flashes. This study provides the first observational evidence of the influence of sunspot oscillations within the umbra on other sunspot waves extending up to the corona. The properties of these waves and oscillations can be utilized to study the inherent magnetic coupling among different layers of the solar atmosphere above sunspots.


title: Cooper Pair Induced Frustration and Nematicity of Two-Dimensional Magnetic Adatom Lattices

abstract: We propose utilizing the Cooper pair to induce magnetic frustration in systems of two-dimensional (2D) magnetic adatom lattices on s-wave superconducting surfaces. The competition between singlet electron correlations and the RKKY coupling is shown to lead to a variety of hidden order states that break the point-group symmetry of the 2D adatom lattice at finite temperature. The phase diagram is constructed using a newly developed effective bond theory [M. Schecter et al., Phys. Rev. Lett. 119, 157202 (2017)], and exhibits broad regions of long-range vestigial nematic order.


title: Topological phase transitions in small mesoscopic chiral p-wave superconductors

abstract: Spin-triplet chiral p-wave superconductivity is typically described by a two-component order parameter, and as such is prone to unique emergent effects when compared to the standard single-component superconductors. Here we present the equilibrium phase diagram for small mesoscopic chiral p-wave superconducting disks in the presence of magnetic field, obtained by solving the microscopic Bogoliubov-de Gennes equations self-consistently. In the ultra-small limit, the cylindrically-symmetric giant-vortex states are the ground state of the system. However, with increasing sample size, the cylindrical symmetry is broken as the two components of the order parameter segregate into domains, and the number of fragmented domain walls between them characterizes the resulting states. Such domain walls are topological defects unique for the p-wave order, and constitute a dominant phase in the mesoscopic regime. Moreover, we find two possible types of domain walls, identified by their chirality-dependent interaction with the edge states.


title: Observational Study on the Fine Structure and Dynamics of a Solar Jet. II. Energy Release Process Revealed by Spectral Analysis

abstract: We report a solar jet phenomenon associated with the C5.4 class flare on 2014 November 11. The data of jet was provided by Solar Dynamics Observatory (SDO), X-Ray Telescope (XRT) aboard Hinode, Interface Region Imaging Spectrograph (IRIS) and Domeless Solar Telescope (DST) at Hida Observatory, Kyoto University. These plentiful data enabled us to present this series of papers to discuss the entire processes of the observed phenomena including the energy storage, event trigger, and energy release. In this paper, we focus on the energy release process of the observed jet, and mainly describe our spectral analysis on the H-alpha data of DST to investigate the internal structure of the H-alpha jet and its temporal evolution. This analysis reveals that in the physical quantity distributions of the H-alpha jet, such as line-of-sight velocity and optical thickness, there is a significant gradient in the direction crossing the jet. We interpret this internal structure as the consequence of the migration of energy release site, based on the idea of ubiquitous reconnection. Moreover, by measuring the horizontal flow of the fine structures in jet, we succeeded in deriving the three-dimensional velocity field and the line-of-sight acceleration field of the H-alpha jet. The analysis result indicates a part of ejecta in the H-alpha jet experienced the additional acceleration after it had been ejected from the lower atmosphere. This secondary acceleration was specified to occur in the vicinity of the intersection between the trajectories of the H-alpha jet and X-ray jet observed by Hinode/XRT. We propose that a fundamental cause of this phenomenon is the magnetic reconnection involving the plasmoid in the observed jet.


title: Modular operads and Batalin-Vilkovisky geometry

abstract: This is a copy of the article published in IMRN (2007). I describe the noncommutative Batalin-Vilkovisky geometry associated naturally with arbitrary modular operad. The classical limit of this geometry is the noncommutative symplectic geometry of the corresponding tree-level cyclic operad. I show, in particular, that the algebras over the Feynman transform of a twisted modular operad P are in one-to-one correspondence with solutions to quantum master equation of Batalin-Vilkovisky geometry on the affine P-manifolds. As an application I give a construction of characteristic classes with values in the homology of the quotient of Deligne-Mumford moduli spaces. These classes are associated naturally with solutions to the quantum master equation on affine S[t]-manifolds, where S[t] is the twisted modular Det-operad constructed from symmetric groups, which generalizes the cyclic operad of associative algebras.


title: Cavity Waters Govern Insulin Association and Release: Inferences from Experimental Data and Molecular Dynamics Simulations

abstract: While a monomer of the ubiquitous hormone insulin is the biologically active form in the human body, its hexameric assembly acts as an efficient storage unit. However, the role of water molecules in the structure, stability and dynamics of the insulin hexamer is poorly understood. Here we combine experimental data with molecular dynamics simulations to investigate the shape, structure and stability of an insulin hexamer focusing on the role of water molecules. Both X-Ray analysis and computer simulations show that the core of the hexamer cavity is barrel-shaped, holding, on an average, sixteen water molecules. These encapsulated and constrained molecules impart structural stability to the hexamer. Apart from the electrostatic interactions with Zn2+ ions, an intricate hydrogen bond network amongst cavity water and neighboring protein residues stabilizes the hexameric association. These water molecules solvate six glutamate residues inside the cavity decreasing electrostatic repulsions amongst the negatively charged carboxylate groups. They also prevent association between glutamate residues and Zn2+ ions and maintain the integrity of the cavity. Simulations reveal that removal of these waters results in a collapse of the cavity. Subsequent analyses also show that the hydrogen bond network among these water molecules and protein residues that face the inner side of the cavity is more rigid with a slower relaxation as compared to that of the bulk solvent. Dynamics of cavity water reveal certain slow water molecules which form the back bone of the stable hydrogen bond network. The analysis presented here suggests a dominant role of structurally conserved water molecules in maintaining the integrity of the hexameric assembly and potentially modulating the dissociation of this assembly into the functional monomeric form.


title: Locally Nameless Permutation Types

abstract: We define "Locally Nameless Permutation Types", which fuse permutation types as used in Nominal Isabelle with the locally nameless representation. We show that this combination is particularly useful when formalizing programming languages where bound names may become free during execution ("extrusion"), common in process calculi. It inherits the generic definition of permutations and support, and associated lemmas, from the Nominal approach, and the ability to stay close to pencil-and-paper proofs from the locally nameless approach. We explain how to use cofinite quantification in this setting, show why reasoning about renaming is more important here than in languages without extrusion, and provide results about infinite support, necessary when reasoning about countable choice.


title: Laboratory experiments on the low temperature formation of carbonaceous grains in the ISM

abstract: The life-cycle of cosmic dust grains is far from being understood and the origin and evolution of interstellar medium (ISM) grains is still under debate. In the ISM, the cosmic dust destruction rate is faster than the production rate by stellar sources. However, observations of ISM refractory matter suggest that to maintain a steady amount of cosmic grains, some supplementary production mechanism takes place. In this context, we aimed to study possible re-formation mechanisms of cosmic grains taking place at low temperature directly in the ISM. The low temperature condensation of carbonaceous materials has been investigated in experiments mimicking the ISM conditions. Gas-phase carbonaceous precursors created by laser ablation of graphite were forced to accrete on cold substrates (T about 10 K) representing surviving dust grains. The growing and evolution of the condensing carbonaceous precursors have been monitored by MIR and UV spectroscopy under a number of experimental scenarios. It is demonstrated, for the first time, the possibility to form ISM carbonaceous grains in "situ". The condensation process is governed by carbon chains that first condense into small carbon clusters and finally into more stable carbonaceous materials, which structural characteristics are comparable to the material formed in gas-phase condensation experiments at very high temperature. We also show that the so-formed fullerene-like carbonaceous material is transformed into a more ordered material under VUV processing. The cold condensation mechanisms here discussed can give fundamental clues to fully understand the balance between the timescale for dust injection, destruction and re-formation in the ISM.


title: Many Paths to Equilibrium: GANs Do Not Need to Decrease a Divergence At Every Step

abstract: Generative adversarial networks (GANs) are a family of generative models that do not minimize a single training criterion. Unlike other generative models, the data distribution is learned via a game between a generator (the generative model) and a discriminator (a teacher providing training signal) that each minimize their own cost. GANs are designed to reach a Nash equilibrium at which each player cannot reduce their cost without changing the other players' parameters. One useful approach for the theory of GANs is to show that a divergence between the training distribution and the model distribution obtains its minimum value at equilibrium. Several recent research directions have been motivated by the idea that this divergence is the primary guide for the learning process and that every step of learning should decrease the divergence. We show that this view is overly restrictive. During GAN training, the discriminator provides learning signal in situations where the gradients of the divergences between distributions would not be useful. We provide empirical counterexamples to the view of GAN training as divergence minimization. Specifically, we demonstrate that GANs are able to learn distributions in situations where the divergence minimization point of view predicts they would fail. We also show that gradient penalties motivated from the divergence minimization perspective are equally helpful when applied in other contexts in which the divergence minimization perspective does not predict they would be helpful. This contributes to a growing body of evidence that GAN training may be more usefully viewed as approaching Nash equilibria via trajectories that do not necessarily minimize a specific divergence at each step.


title: Conformality Lost in Efimov Physics

abstract: A general mechanism for the loss of conformal invariance is the merger and disappearance of an infrared fixed point and an ultraviolet fixed point of a renormalization group flow. We show explicitly how this mechanism works in the case of identical bosons at unitarity as the spatial dimension $d$ is varied. For $d$ between the critical dimensions $d_{\rm 1}=2.30$ and $d_{\rm 2}=3.76$, there is loss of conformality as evidenced by the Efimov effect in the three-body sector. The beta function for an appropriate three-body coupling is a quadratic polynomial in that coupling. For $d<d_{\rm 1}$ and for $d>d_{\rm 2}$, the beta function has two real roots that correspond to infrared and ultraviolet fixed points. As $d$ approaches $d_{\rm 1}$ from below and as $d$ approaches $d_{\rm 2}$ from above, the fixed points merge and disappear into the complex plane. For $d_{\rm 1}<d<d_{\rm 2}$, the beta function has complex roots and the renormalization group flow for the three-body coupling is a limit cycle.


title: AdS$_3$ metric from UV/IR entanglement entropies of CFT$_2$

abstract: How to construct the $d+1$ dimensional geometry explicitly from the dual CFT$_d$ is a widely concerned problem. Specifically, given entanglement entropies of a CFT$_2$, which is purely expressed by two dimensional parameters, can we build the dual three dimensional geometry unambiguously? To do this, one must assume nothing is known about the three dimensional geometry and starts with the most general setup. In this paper, by identifying the UV and IR entanglement entropies of a perturbed CFT$_2$ with the geodesic lengths, we show that, the dual geometry is uniquely determined to be asymptotic AdS$_3$. The holographic dimension is generated by the energy cut-off of the CFT$_2$, according to the holographic principle. The pure AdS$_3$ is obtained by taking the massless limit. Our derivations apply to both static and covariant scenarios. Moreover, what deserves special attention is that the ratio of the numberical factors of the UV/IR entanglement entropies are crucial to have a dual geometry. We are led to conjecture a necessary condition of holographic CFT$_2$.

categories: hep-th cond-mat.stat-mech gr-qc quant-ph


title: BPS spectra from BPS graphs

abstract: I present a simple graphical method to find the BPS spectra of $A_1$ theories of class S. BPS graphs provide a bridge between spectral networks and BPS quivers, the two main frameworks for the study of BPS states. Here I show how to essentially read off from a BPS graph the quantum spectrum generator (or BPS monodromy), expressed as a product of quantum dilogarithms. Thanks to the framed wall-crossing phenomenon for line defects, the determination of the BPS spectrum reduces to the computation of quantum parallel transport across the edges of the BPS graph.


title: $\epsilon$-Monotone Fourier Methods for Optimal Stochastic Control in Finance

abstract: Stochastic control problems in finance often involve complex controls at discrete times. As a result numerically solving such problems, for example using methods based on partial differential or integro-differential equations, inevitably give rise to low order accuracy, usually at most second order. In many cases one can make use of Fourier methods to efficiently advance solutions between control monitoring dates and then apply numerical optimization methods across decision times. However Fourier methods are not monotone and as a result give rise to possible violations of arbitrage inequalities. This is problematic in the context of control problems, where the control is determined by comparing value functions. In this paper we give a preprocessing step for Fourier methods which involves projecting the Green's function onto the set of linear basis functions. The resulting algorithm is guaranteed to be monotone (to within a tolerance), $\ell_\infty$-stable and satisfies an $\epsilon$-discrete comparison principle. In addition the algorithm has the same complexity per step as a standard Fourier method while at the same time having second order accuracy for smooth problems.


title: Combining Lexical Features and a Supervised Learning Approach for Arabic Sentiment Analysis

abstract: The importance of building sentiment analysis tools for Arabic social media has been recognized during the past couple of years, especially with the rapid increase in the number of Arabic social media users. One of the main difficulties in tackling this problem is that text within social media is mostly colloquial, with many dialects being used within social media platforms. In this paper, we present a set of features that were integrated with a machine learning based sentiment analysis model and applied on Egyptian, Saudi, Levantine, and MSA Arabic social media datasets. Many of the proposed features were derived through the use of an Arabic Sentiment Lexicon. The model also presents emoticon based features, as well as input text related features such as the number of segments within the text, the length of the text, whether the text ends with a question mark or not, etc. We show that the presented features have resulted in an increased accuracy across six of the seven datasets we've experimented with and which are all benchmarked. Since the developed model out-performs all existing Arabic sentiment analysis systems that have publicly available datasets, we can state that this model presents state-of-the-art in Arabic sentiment analysis.


title: Mesoporous silica obtained with methyltriethoxysilane as co-precursor in alkaline medium

abstract: Mesoporous silica particles have been synthesized by sol-gel method from tetraethoxysilane (tetraethylorthosilicate, TEOS) and methyltriethoxysilane (MTES), in ethanol and water mixture, at different ratios of the of the silica precursors. Ammonia was used as catalyst at room temperature and hexadecyltrimethylammonium bromide (cetyltrimethylammonium bromide, CTAB) as the structure directing agent. Nitrogen sorption, X-ray diffraction and small-angle neutron scattering gave information on the evolution of the gel structure and pore morphologies in the function of MTES/TEOS molar ratio. Thermogravimetric and differential thermal analysis showed that with addition of MTES the exothermic peak indicating the oxidation of the low molecular weight organic fragments shift to higher temperature. A room-temperature, one-pot synthesis of MCM-41 type materials is presented, in which the variation of the MTES concentration allows to change the hydrophobicity, preserving the specific properties materials, like the ordered pore structure, large specific surface area and high porosity, making them suitable for selective uptake of guest species in drug loading applications. Specifically, the obtained materials had cylindrical pores, specific surface areas up to 1101 m2/g and total pore volumes up to 0.473 cm3/g. The obtained mesoporous materials are susceptible for further functionalization to improve their selective uptake of guest species in drug delivery applications.


abstract: This brief survey aims to set the stage and summarize some of the ideas under discussion at the Workshop on Singular Geometry and Higgs Bundles in String Theory, to be held at the American Institute of Mathematics from October 30th to November 3rd, 2017. One of the most interesting aspects of the duality revolution in string theory is the understanding that gauge fields and matter representations can be described by intersection of branes. Since gauge theory is at the heart of our description of physical interactions, it has opened the door to the geometric engineering of many physical systems, and in particular those involving Higgs bundles. This note presents a curated overview of some current advances and open problems in the area, with no intention of being a complete review of the whole subject.


title: 't Hooft mechanism, anomalous gauge U(1), and "invisible" axion from string

abstract: Among solutions of the strong CP problem, the "invisible" axion in the narrow axion window is argued to be the remaining possibility among natural solutions on the smallness of $\bar{\theta}$. Related to the gravity spoil of global symmetries, some prospective "invisible" axions from theory point of view are discussed. In all these discussions including the observational possibility and cosmological constraints, including a safe domain wall problem, must be included.


title: The Unbounded Integrality Gap of a Semidefinite Relaxation of the Traveling Salesman Problem

abstract: We study a semidefinite programming relaxation of the traveling salesman problem introduced by de Klerk, Pasechnik, and Sotirov [8] and show that their relaxation has an unbounded integrality gap. In particular, we give a family of instances such that the gap increases linearly with $n$. To obtain this result, we search for feasible solutions within a highly structured class of matrices; the problem of finding such solutions reduces to finding feasible solutions for a related linear program, which we do analytically. The solutions we find imply the unbounded integrality gap. Further, they imply several corollaries that help us better understand the semidefinite program and its relationship to other TSP relaxations. Using the same technique, we show that a more general semidefinite program introduced by de Klerk, de Oliveira Filho, and Pasechnik [7] for the $k$-cycle cover problem also has an unbounded integrality gap.


title: Magnetic Sensitivity of AlMn TESes and Shielding Considerations for Next-Generation CMB Surveys

abstract: In the next decade, new ground-based Cosmic Microwave Background (CMB) experiments such as Simons Observatory (SO), CCAT-prime, and CMB-S4 will increase the number of detectors observing the CMB by an order of magnitude or more, dramatically improving our understanding of cosmology and astrophysics. These projects will deploy receivers with as many as hundreds of thousands of transition edge sensor (TES) bolometers coupled to Superconducting Quantum Interference Device (SQUID)-based readout systems. It is well known that superconducting devices such as TESes and SQUIDs are sensitive to magnetic fields. However, the effects of magnetic fields on TESes are not easily predicted due to the complex behavior of the superconducting transition, which motivates direct measurements of the magnetic sensitivity of these devices. We present comparative four-lead measurements of the critical temperature versus applied magnetic field of AlMn TESes varying in geometry, doping, and leg length, including Advanced ACT (AdvACT) and POLARBEAR-2/Simons Array bolometers. Molybdenum-copper bilayer ACTPol TESes are also tested and are found to be more sensitive to magnetic fields than the AlMn devices. We present an observation of weak-link-like behavior in AlMn TESes at low critical currents. We also compare measurements of magnetic sensitivity for time division multiplexing SQUIDs and frequency division multiplexing microwave rf-SQUIDs. We discuss the implications of our measurements on the magnetic shielding required for future experiments that aim to map the CMB to near-fundamental limits.


title: Geometric Horizons

abstract: We discuss black hole spacetimes with a geometrically defined quasi-local horizon on which the curvature tensor is algebraically special relative to the alignment classification. Based on many examples and analytical results, we conjecture that a spacetime horizon is always more algebraically special (in all of the orders of specialization) than other regions of spacetime. Using recent results in invariant theory, such geometric black hole horizons can be identified by the alignment type $\textbf{II}$ or $\textbf{D}$ discriminant conditions in terms of scalar curvature invariants, which are not dependent on spacetime foliations. The above conjecture is, in fact, a suite of conjectures (isolated vs dynamical horizon; four vs higher dimensions; zeroth order invariants vs higher order differential invariants). However, we are particularly interested in applications in four dimensions and especially the location of a black hole in numerical computations.


title: NileTMRG at SemEval-2017 Task 4: Arabic Sentiment Analysis

abstract: This paper describes two systems that were used by the authors for addressing Arabic Sentiment Analysis as part of SemEval-2017, task 4. The authors participated in three Arabic related subtasks which are: Subtask A (Message Polarity Classification), Sub-task B (Topic-Based Message Polarity classification) and Subtask D (Tweet quantification) using the team name of NileTMRG. For subtask A, we made use of our previously developed sentiment analyzer which we augmented with a scored lexicon. For subtasks B and D, we used an ensemble of three different classifiers. The first classifier was a convolutional neural network for which we trained (word2vec) word embeddings. The second classifier consisted of a MultiLayer Perceptron, while the third classifier was a Logistic regression model that takes the same input as the second classifier. Voting between the three classifiers was used to determine the final outcome. The output from task B, was quantified to produce the results for task D. In all three Arabic related tasks in which NileTMRG participated, the team ranked at number one.


title: Computational Method for a Fractional Model of the Helium Burning Network

abstract: Stellar cores may be considered as a nuclear reactor that play important role in injecting new synthesized elements in the interstellar medium Helium burning is an important stage that contribute to the synthesis of key elements such as carbon, through the triple- {\alpha} process, and oxygen. In the present paper, we introduce a computational method for the fractional model of the nuclear helium burning in stellar cores. The system of fractional differential equations is solved simultaneously using series expansion method. The calculations are performed in the sense of modified Riemann-Liouville fractional derivative. Analytic expressions are obtained for the abundance of each element as a function of time. Comparing the abundances calculated at the fractional parameter \alpha=1 , which represents the integer solution, with the numerical solution revealed a good agreement with maximum error \epsilon =0.003. The product abundances are calculated at\alpha=0.25, 0.5, 0.75 to declare the effects of changing the fractional parameters on the calculations.


title: Serving deep learning models in a serverless platform

abstract: Serverless computing has emerged as a compelling paradigm for the development and deployment of a wide range of event based cloud applications. At the same time, cloud providers and enterprise companies are heavily adopting machine learning and Artificial Intelligence to either differentiate themselves, or provide their customers with value added services. In this work we evaluate the suitability of a serverless computing environment for the inferencing of large neural network models. Our experimental evaluations are executed on the AWS Lambda environment using the MxNet deep learning framework. Our experimental results show that while the inferencing latency can be within an acceptable range, longer delays due to cold starts can skew the latency distribution and hence risk violating more stringent SLAs.


title: The airglow layer emission altitude cannot be determined unambiguously from temperature comparison with lidars

abstract: I investigate the nightly mean emission height and width of the OH*(3-1) layer by comparing nightly mean temperatures measured by the ground-based spectrometer GRIPS 9 and the Na lidar at ALOMAR. The data set contains 42 coincident measurements between November 2010 and February 2014, when GRIPS 9 was in operation at the ALOMAR observatory (69.3$^\circ$N, 16.0$^\circ$E) in northern Norway. To closely resemble the mean temperature measured by GRIPS 9, I weight each nightly mean temperature profile measured by the lidar using Gaussian distributions with 40 different centre altitudes and 40 different full widths at half maximum. In principle, one can thus determine the altitude and width of an airglow layer by finding the minimum temperature difference between the two instruments. On most nights, several combinations of centre altitude and width yield a temperature difference of $\pm$2 K. The generally assumed altitude of 87 km and width of 8 km is never an unambiguous, good solution for any of the measurements. Even for a fixed width of $\sim$8.4 km, one can sometimes find several centre altitudes that yield equally good temperature agreement. Weighted temperatures measured by lidar are not suitable to determine unambiguously the emission height and width of an airglow layer. However, when actual altitude and width data are lacking, a comparison with lidars can provide an estimate of how representative a measured rotational temperature is of an assumed altitude and width. I found the rotational temperature to represent the temperature at the commonly assumed altitude of 87.4 km and width of 8.4 km to within $\pm$16 K, on average. This is not a measurement uncertainty.


title: What confines the rings of Saturn?

abstract: The viscous spreading of planetary rings is believed to be counteracted by satellite torques, either through an individual resonance or through overlapping resonances. For the A ring of Saturn, it has been commonly believed that the satellite Janus alone can prevent the ring from spreading via its 7:6 Lindblad resonance. We discuss this common misconception and show that, in reality, the A ring is confined by the contributions from the group of satellites Pan, Atlas, Prometheus, Pandora, Janus, Epimetheus, and Mimas, whose cumulative torques from various resonances gradually decrease the angular momentum flux transported outward through the ring via density and bending waves. We further argue that this decrease in angular momentum flux occurs through 'flux reversal'. Furthermore, we use the magnitude of the satellites' resonance torques to estimate the effective viscosity profile across the A ring, showing that it decreases with radius from ~50 cm2 s-1 to less than ~10 cm2 s-1. The gradual estimated decrease of the angular momentum flux and effective viscosity are roughly consistent with results obtained by balancing the shepherding torques from Pan and Daphnis with the viscous torque at the edges of the Encke and Keeler gaps, as well as the edge of the A ring. On the other hand, the Mimas 2:1 Lindblad resonance alone seems to be capable of confining the edge of the B ring, and contrary to the situation in the A ring, we show that the effective viscosity across the B ring is relatively constant at ~24-30 cm2 s-1.


title: How to encode a tree

abstract: We construct bijections giving three "codes" for trees. These codes follow naturally from the Matrix Tree Theorem of Tutte and have many advantages over the one produced by Prufer in 1918. One algorithm gives explicitly a bijection that is implicit in Orlin's manipulatorial proof of Cayley's formula (the formula was actually found first by Borchardt). Another is based on a proof of Knuth. The third is an implementation of Joyal's pseudo-bijective proof of the formula, and is equivalent to one previously found by Egecioglu and Remmel. In each case, we have at least two algorithms, one of which involves hands-on manipulations of the tree while the other involves a combinatorial and linear algebraic manipulation of a matrix.


title: Interpretable Machine Learning for Privacy-Preserving Pervasive Systems

abstract: Our everyday interactions with pervasive systems generate traces that capture various aspects of human behavior and enable machine learning algorithms to extract latent information about users. In this paper, we propose a machine learning interpretability framework that enables users to understand how these generated traces violate their privacy.


title: Complex Langevin Simulations of QCD at Finite Density -- Progress Report

abstract: We simulate lattice QCD at finite quark-number chemical potential to study nuclear matter, using the complex Langevin equation (CLE). The CLE is used because the fermion determinant is complex so that standard methods relying on importance sampling fail. Adaptive methods and gauge-cooling are used to prevent runaway solutions. Even then, the CLE is not guaranteed to give correct results. We are therefore performing extensive testing to determine under what, if any, conditions we can achieve reliable results. Our earlier simulations at $\beta=6/g^2=5.6$, $m=0.025$ on a $12^4$ lattice reproduced the expected phase structure but failed in the details. Our current simulations at $\beta=5.7$ on a $16^4$ lattice fail in similar ways while showing some improvement. We are therefore moving to even weaker couplings to see if the CLE might produce the correct results in the continuum (weak-coupling) limit, or, if it still fails, whether it might reproduce the results of the phase-quenched theory. We also discuss action (and other dynamics) modifications which might improve the performance of the CLE.


title: Optimal Control of Coefficients in Parabolic Free Boundary Problems Modeling Laser Ablation

abstract: Inverse Stefan problem arising in modeling of laser ablation of biomedical tissues is analyzed, where information on the coefficients, heat flux on the fixed boundary, and density of heat sources are missing and must be found along with the temperature and free boundary. Optimal control framework is employed, where the missing data and the free boundary are components of the control vector, and optimality criteria are based on the final moment measurement of the temperature and position of the free boundary. Discretization by finite differences is pursued, and convergence of the discrete optimal control problems to the original problem is proven.


title: Full-Duplex Heterogeneous Networks with Decoupled User Association: Rate Analysis and Traffic Scheduling

abstract: Full-duplex (FD) transmission in a point-to-point (P2P) link, wherein bidirectional traffic flows simultaneously share the same spectrum, has the capability of doubling the link rate by completely removing self-interferences. However, the rate gain of FD transmission in an interference-limited network is not as clear as a P2P link due to the interferences induced by complex FD and half-duplex (HD) transmission behaviors in the network. To thoroughly investigate the achievable link rates of users and base stations (BSs) in a heterogeneous network (HetNet) with decoupled user association, a general approach is proposed to model a HetNet in which all BSs and users can perform HD or FD transmission depending on their traffic patterns. We first characterize the decoupled rate-optimal user association scheme and use it to define the downlink and uplink rates in the HetNet. The tight lower bounds on the link rates of the FD users and BSs are found in a neat form that characterizes general channel fading, user and BS intensities and imperfect self-interference cancellation. These bounds outline the rate regions of the FD users that inspire proposing the opportunistic FD scheduling algorithms maximizing the sum rate of each bidirectional traffic and stabilize each of the queues in the HetNet.


title: Laws relating runs, long runs, and steps in gambler's ruin, with persistence in two strata

abstract: Define a certain gambler's ruin process $\mathbf{X}_{j}, \mbox{ \ }j\ge 0,$ such that the increments $\varepsilon_{j}:=\mathbf{X}_{j}-\mathbf{X}_{j-1}$ take values $\pm1$ and satisfy $P(\varepsilon_{j+1}=1|\varepsilon_{j}=1, |\mathbf{X}_{j}|=k)=P(\varepsilon_{j+1}=-1|\varepsilon_{j}=-1,|\mathbf{X}_{j}|=k)=a_k$, all $j\ge 1$, where $a_k=a$ if $ 0\le k\le f-1$, and $a_k=b$ if $f\le k<N$. Here $0<a, b <1$ denote persistence parameters and $ f ,N\in \mathbb{N} $ with $f<N$. The process starts at $\mathbf{X}_0=m\in (-N,N)$ and terminates when $|\mathbf{X}_j|=N$. Denote by ${\cal R}'_N$, ${\cal U}'_N$, and ${\cal L}'_N$, respectively, the numbers of runs, long runs, and steps in the meander portion of the gambler's ruin process. Define $X_N:=\left ({\cal L}'_N-\frac{1-a-b}{(1-a)(1-b)}{\cal R}'_N-\frac{1}{(1-a)(1-b)}{\cal U}'_N\right )/N$ and let $f\sim\eta N$ for some $0<\eta <1$. We show $\lim_{N\to\infty} E\{e^{itX_N}\}=\hat{\varphi}(t)$ exists in an explicit form. We obtain a companion theorem for the last visit portion of the gambler's ruin.


title: Charged kaon femtoscopic correlations in p$-$Pb collisions at $\mathbf{\sqrt{{\textit s}_{\rm NN}}}=5.02$ TeV with ALICE at the LHC

abstract: Particle correlations at small relative momenta can be used to measure the space$-$time characteristics of particle production on the femtoscopic (fm=10$^{-15}$ m) level in high-energy collisions. Such correlations arise due to quantum statistics and final state interactions. We report correlations of two charged identical kaons measured in p$-$Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV by the ALICE experiment at the LHC. The femtoscopic invariant radii and correlation strengths are extracted from the one-dimensional kaon correlation functions and are compared to those obtained in pp and Pb$-$Pb collisions at $\sqrt{s}=7$ TeV and $\sqrt{s_{\rm NN}}=2.76$ TeV, respectively. Kaon femtoscopy in p$-$Pb collisions is an important supplement to that in pp and Pb$-$Pb collisions because it allows one to understand the particle production mechanisms in different collision systems. It also complements the existing pion correlation results. The obtained radii increase with increasing multiplicity and decrease with increasing pair transverse momentum. At comparable multiplicity, the radii measured in p$-$Pb collisions are close to those observed in pp collisions, whereas it is difficult to compare them with the results from Pb$-$Pb collisions because of a big multiplicity gap.


title: Measurement of shower development and its $\textit{Moli\`ere radius}$ with a four-plane LumiCal test set-up

abstract: A prototype of a luminometer, designed for a future $e^{+}e^{-}$ collider detector, and consisting at present of a four-plane module, was tested in the CERN PS accelerator T9 beam. The objective of this beam test was to demonstrate a multi-plane operation, to study the development of the electromagnetic shower and to compare it with MC simulations. In addition, the effective $\textit{Moli\`ere radius}$ of this configuration is extracted.


title: Communication-avoiding Cholesky-QR2 for rectangular matrices

abstract: Scalable QR factorization algorithms for solving least squares and eigenvalue problems are critical given the increasing parallelism within modern machines. We introduce a more general parallelization of the CholeskyQR2 algorithm and show its effectiveness for a wide range of matrix sizes. Our algorithm executes over a 3D processor grid, the dimensions of which can be tuned to trade-off costs in synchronization, interprocessor communication, computational work, and memory footprint. We implement this algorithm, yielding a code that can achieve a factor of $\Theta(P^{1/6})$ less interprocessor communication on $P$ processors than any previous parallel QR implementation. Our performance study on Intel Knights-Landing and Cray XE supercomputers demonstrates the effectiveness of this CholeskyQR2 parallelization on a large number of nodes. Specifically, relative to ScaLAPACK's QR, on 1024 nodes of Stampede2, our CholeskyQR2 implementation is faster by 2.6x-3.3x in strong scaling tests and by 1.1x-1.9x in weak scaling tests.


title: The Hausdorff Mapping Is Nonexpanding

abstract: In the present paper we investigate the properties of the Hausdorff mapping $\mathcal{H}$, which takes each compact metric space to the space of its nonempty closed subspaces. It is shown that this mapping is nonexpanding (Lipschitz mapping with constant $1$). This paper gives several examples of classes of metric spaces, distances between which are preserved by the mapping $\mathcal{H}$. We also calculate distance between any connected metric space and any simplex with greater diameter than the former one. At the end of the paper we discuss some properties of the Hausdorff mapping which may help to prove that it is isometric


title: A Unified Framework for Long Range and Cold Start Forecasting of Seasonal Profiles in Time Series

abstract: Providing long-range forecasts is a fundamental challenge in time series modeling, which is only compounded by the challenge of having to form such forecasts when a time series has never previously been observed. The latter challenge is the time series version of the cold-start problem seen in recommender systems which, to our knowledge, has not been addressed in previous work. A similar problem occurs when a long range forecast is required after only observing a small number of time points --- a warm start forecast. With these aims in mind, we focus on forecasting seasonal profiles---or baseline demand---for periods on the order of a year in three cases: the long range case with multiple previously observed seasonal profiles, the cold start case with no previous observed seasonal profiles, and the warm start case with only a single partially observed profile. Classical time series approaches that perform iterated step-ahead forecasts based on previous observations struggle to provide accurate long range predictions; in settings with little to no observed data, such approaches are simply not applicable. Instead, we present a straightforward framework which combines ideas from high-dimensional regression and matrix factorization on a carefully constructed data matrix. Key to our formulation and resulting performance is leveraging (1) repeated patterns over fixed periods of time and across series, and (2) metadata associated with the individual series; without this additional data, the cold-start/warm-start problems are nearly impossible to solve. We demonstrate that our framework can accurately forecast an array of seasonal profiles on multiple large scale datasets.


title: About heavy neutrinos: Lepton-flavor violation in decays of charged leptons

abstract: The fundamental description of nature, beyond the Standard Model (SM), may include heavy neutrinos that mix and thus allow processes in which lepton flavor is not preserved. We investigate the impact of charged currents that couple heavy gauge bosons to heavy neutrinos and SM leptons on lepton-flavor-violating decays of SM leptons into three charged leptons, with no final-state neutrinos. We implement our expressions for the leading contributions to ${\rm Br}(l_\alpha\to l_\beta\,l_\sigma\,l_\sigma)$, which hold for either Dirac or Majorana neutrinos, to the trilepton decay $\mu\to3e$, of the muon, and so determine sets of masses of heavy neutrinos and the heavy gauge boson, within GeVs to few TeVs, that are consistent with the upper bounds provided by the SINDRUM Collaboration. We find, however, that constraints dictated by the upper bound on ${\rm Br}(\mu\to e\gamma)$, from the MEG Collaboration, are more stringent. We utilize such parameters to find that the contributions to tau decays are $\sim10^{-15}-10^{-13}$, well below bounds from $B$ factories. The mixing of heavy and SM charged bosons is also investigated. We find that current experimental data from MEG and SINDRUM would allow mixing angles as large as $\sim10^{-2}$, for a relatively light new charged boson, but the expected sensitivity of the Mu3e experiment would be capable of setting an upper bound on this angle as small as $\sim10^{-4}$ if the mass of this boson is within the range of few TeVs.


title: Composition of PPT Maps

abstract: M. Christandl conjectured that the composition of any trace preserving PPT map with itself is entanglement breaking. We prove that Christandl's conjecture holds asymptotically by showing that the distance between the iterates of any unital or trace preserving PPT map and the set of entanglement breaking maps tends to zero. Finally, for every graph we define a one-parameter family of maps on matrices and determine the least value of the parameter such that the map is variously, positive, completely positive, PPT and entanglement breaking in terms of properties of the graph. Our estimates are sharp enough to conclude that Christandl's conjecture holds for these families.


title: Platoon Stability and Safety Analysis of Cooperative Adaptive Cruise Control under Wireless Rician Fading Channels and Jamming Attacks

abstract: Cooperative Adaptive Cruise Control (CACC) is considered as a key enabling technology to automatically regulate the inter-vehicle distances in a vehicle platoon to improve traffic efficiency while maintaining safety. Although the wireless communication and physical processes in the existing CACC systems are integrated in one control framework, the coupling between wireless communication reliability and system states is not well modeled. Furthermore, the research on the impact of jamming attacks on the system stability and safety is largely open. In this paper, we conduct a comprehensive analysis on the stability and safety of the platoon under the wireless Rician fading channel model and jamming attacks. The effect of Rician fading and jamming on the communication reliability is incorporated in the modeling of string dynamics such that it captures its state dependency. Time-domain definition of string stability is utilized to delineate the impact of Rician fading and jamming on the CACC system's functionality and string stability. Attacker's possible locations at which it can destabilize the string is further studied based on the proposed model. From the safety perspective, reachable states (i.e., inter-vehicle distances) of the CACC system under unreliable wireless fading channels and jamming attacks is studied. Safety verification is investigated by examining the inter-vehicle distance trajectories. We propose a methodology to compute the upper and lower bounds of the trajectories of inter-vehicle distances between the lead vehicle and its follower. We conduct extensive simulations to evaluate the system stability and safety under jamming attacks in different scenarios. We identify that channel fading can degrade the performance of the CACC system, and the platoon's safety is highly sensitive to jamming attacks.


title: Detailed Abundances for the Old Population near the Galactic Center: I. Metallicity distribution of the Nuclear Star Cluster

abstract: We report the first high spectral resolution study of 17 M giants kinematically confirmed to lie within a few parsecs of the Galactic Center, using R=24,000 spectroscopy from Keck/NIRSPEC and a new linelist for the infrared K band. We consider their luminosities and kinematics, which classify these stars as members of the older stellar population and the central cluster. We find a median metallicity of <[Fe/H]>=-0.16 and a large spread from approximately -0.3 to +0.3 (quartiles). We find that the highest metallicities are [Fe/H]<+0.6, with most of the stars being at or below the Solar iron abundance. The abundances and the abundance distribution strongly resembles that of the Galactic bulge rather than disk or halo; in our small sample we find no statistical evidence for a dependence of velocity dispersion on metallicity.


title: Key Technologies and System Trade-Offs for Detection and Localization of Amateur Drones

abstract: The use of amateur drones (ADrs) is expected to significantly increase over the upcoming years. However, regulations do not allow such drones to fly over all areas, in addition to typical altitude limitations. As a result, there is an urgent need for ADrs surveillance solutions. These solutions should include means of accurate detection, classification, and localization of the unwanted drones in a no-fly zone. In this paper, we give an overview of promising techniques for modulation classification and signal strength based localization of ADrs by using surveillance drones (SDrs). By introducing a generic altitude dependent propagation model, we show how detection and localization performance depend on the altitude of SDrs. Particularly, our simulation results show a 25 dB reduction in the minimum detectable power or 10 times coverage enhancement of an SDr by flying at the optimum altitude. Moreover, for a target no-fly zone, the location estimation error of an ADr can be remarkably reduced by optimizing the positions of the SDrs. Finally, we conclude the paper with a general discussion about the future work and possible challenges of the aerial surveillance systems.


title: The velocity ellipsoid in the Galactic disc using Gaia DR1

abstract: The stellar velocity ellipsoid of the solar neighbour is re-examined using intermediate-old mono-abundances stellar groups with high quality chemistry data together with parallaxes and proper motions from Gaia DR1. We find the average velocity dispersion values for the three space velocity components for the thin and thick disc of (\sigma_{U},\sigma_{V},\sigma_{W})_{thin} = (33 \pm 4, 28 \pm 2, 23 \pm 2) and (\sigma_{U},\sigma_{V},\sigma_{W})_{thick} = (57 \pm 6, 38 \pm 5, 37 \pm 4) km s^{-1}, respectively. The mean values of the ratio between the semi-axes of the velocity ellipsoid for the thin disc are found to be, \sigma_{V}/\sigma_{U} = 0.70 \pm 0.13 and \sigma_{W}/\sigma_{U} is 0.64 \pm 0.08, while for the thick disc \sigma_{V}/\sigma_{U} = 0.67 \pm 0.11 and \sigma_{W}/\sigma_{U} is 0.66 \pm 0.11. Inputting these dispersions into the linear Str\"omberg relation for the thin disc groups, we find the Sun's velocity with respect to the LSR in Galactic rotation to be V_{\sun} = 13.9 \pm 3.4 km s^{-1}. A relation is found between the vertex deviation and the chemical abundances for the thin disc, ranging from -5 to +40^{\circ} as iron-abundance increases. For the thick disc we find a vertex deviation of l_{uv} \sim -15^{\circ}. The tilt angle (l_{uw}) in the U-W plane for the thin disc groups ranges from -10 to +15^\circ, but there is no evident relation between l_{uw} and the mean abundances. However we find a weak relation for l_{uw} as a function of iron abundances and \alpha-elements for most of the groups in the thick disc, where the tilt angle decreases from -5 to -20^\circ when [Fe/H] decreases and [\alpha/Fe] increases. The velocity anisotropy parameter is independent of the chemical group abundances and its value is nearly constant for both discs (\beta \sim 0.5), suggesting that the combined disc is dynamically relaxed.


title: The generalized Sierpi\'{n}ski Arrowhead Curve

abstract: We define special Hamiltonian-paths and special permutations of the up-facing dark tiles on a checked triangular grid related to the generalized Sierpi\'{n}ski Gasket. Our definitions and observations make possible the generalization of the Sierpi\'{n}ski Arrowhead Curve for all orders. We produce these symmetric recursive curves in many ways by two kinds of asymmetric paths which are in a bijective relation and unambiguously transformable into each other in any order. These node-rewriting and edge-rewriting recursive curves keep their self-avoiding and simple properties after the transformation and their cardinality specifies a new integer sequence. We show a transformation table to change the curves into each other and we give another table to change them into Lindenmayer-system strings both by the absolute direction codes of their edges.


title: Psyquandles, Singular Knots and Pseudoknots

abstract: We generalize the notion of biquandles to psyquandles and use these to define invariants of oriented singular links and pseudolinks. In addition to psyquandle counting invariants, we introduce Alexander psyquandles and corresponding invariants such as Alexander psyquandle polynomials and Alexander-Gr\"obner psyquandle invariants of oriented singular knots and links. We consider the relationship between Alexander psyquandle colorings of pseudolinks and p-colorings of pseudolinks. As a special case we define a generalization of the Alexander polynomial for oriented singular links and pseudolinks we call the Jablan polynomial and compute the invariant for all pseudoknots with up to five crossings and all 2-bouquet graphs with up to 6 classical crossings.


title: Search for superfluidity in supercooled liquid parahydrogen

abstract: The possible superfluid transition of supercooled liquid parahydrogen is investigated by quantum Monte Carlo simulations. The cooling protocol adopted here allows for the investigation of a fluid phase down to a temperature T=0.25 K. No evidence of superfluidity is found, as exchanges of identical particles are strongly suppressed even at the lowest temperature. Is shown that, contrary to a commonly held belief, it is not the well depth of the pair-wise interaction but rather its relatively large hard core diameter that physically hinders superfluidity in parahydrogen.


title: On the Relaxation of Hybrid Dynamical Systems

abstract: Hybrid dynamical systems have proven to be a powerful modeling abstraction, yet fundamental questions regarding the dynamical properties of these systems remain. In this paper, we develop a novel class of relaxations which we use to recover a number of classic systems theoretic properties for hybrid systems, such as existence and uniqueness of trajectories, even past the point of Zeno. Our relaxations also naturally give rise to a class of provably convergent numerical approximations, capable of simulating through Zeno. Using our methods, we are also able to perform sensitivity analysis about nominal trajectories undergoing a discrete transition -- a technique with many practical applications, such as assessing the stability of periodic orbits.


title: Homological stability of topological moduli spaces

abstract: Given a graded $E_1$-module over an $E_2$-algebra in spaces, we construct an augmented semi-simplicial space up to higher coherent homotopy over it, called its canonical resolution, whose graded connectivity yields homological stability for the graded pieces of the module with respect to constant and abelian coefficients. We furthermore introduce a notion of coefficient systems of finite degree in this context and show that, without further assumptions, the corresponding twisted homology groups stabilize as well. This generalizes a framework of Randal-Williams and Wahl for families of discrete groups. In many examples, the canonical resolution recovers geometric resolutions with known connectivity bounds. As a consequence, we derive new twisted homological stability results for e.g. moduli spaces of high-dimensional manifolds, unordered configuration spaces of manifolds with labels in a fibration, and moduli spaces of manifolds equipped with unordered embedded discs. This in turn implies representation stability for the ordered variants of the latter examples.


title: Making Faces: Universal Inverse Design of Surfaces with Thin Nematic Elastomer Sheets

abstract: Programmable shape-shifting materials can take different physical forms to achieve multifunctionality in a dynamic and controllable manner. Although morphing a shape from 2D to 3D via programmed inhomogeneous local deformations has been demonstrated in various ways, the inverse problem -- programming a sheet to take an arbitrary desired 3D shape -- is much harder yet critical to realize specific functions. Here, we address this inverse problem in thin liquid crystal elastomer (LCE) sheets, where the shape is preprogrammed by precise and local control of the molecular orientation of the liquid crystal monomers. We show how blueprints for arbitrary surface geometries as well as local extrinsic curvatures can be generated using approximate numerical methods. Backed by faithfully alignable and rapidly lockable LCE chemistry, we precisely embed our designs in LCE sheets using advanced top-down microfabrication techniques. We thus successfully produce flat sheets that, upon thermal activation, take an arbitrary desired shape, such as a face. The general design principles presented here for creating an arbitrary 3D shape will allow for exploration of unmet needs in flexible electronics, metamaterials, aerospace and medical devices, and more.


title: Decomposing graphs into edges and triangles

abstract: We prove the following 30-year old conjecture of Gy\H{o}ri and Tuza: the edges of every $n$-vertex graph $G$ can be decomposed into complete graphs $C_1,\ldots,C_\ell$ of orders two and three such that $|C_1|+\cdots+|C_\ell|\le (1/2+o(1))n^2$. This result implies the asymptotic version of the old result of Erd\H{o}s, Goodman and P\'osa that asserts the existence of such a decomposition with $\ell\le n^2/4$.


title: On the F-expanding of Homoclinic class

abstract: We establish a closing property for thin trapped homoclinic classes. Taking advantage of this property, we proved that if the homoclinic class $H(p)$ admits a dominated splitting $T_{H(p)}M=E\oplus_{<}F$, where $E$ is thin trapped (see Definition \ref{Def:TP}) and all periodic points homoclinically related to $p$ are uniformly $F$-expanding at the period (see Definition \ref{Def:expanding}), then $F$ is expanded (see Definition \ref{Def:TP}).


title: An FPT Algorithm Beating 2-Approximation for $k$-Cut

abstract: In the $k$-Cut problem, we are given an edge-weighted graph $G$ and an integer $k$, and have to remove a set of edges with minimum total weight so that $G$ has at least $k$ connected components. Prior work on this problem gives, for all $h \in [2,k]$, a $(2-h/k)$-approximation algorithm for $k$-cut that runs in time $n^{O(h)}$. Hence to get a $(2 - \varepsilon)$-approximation algorithm for some absolute constant $\varepsilon$, the best runtime using prior techniques is $n^{O(k\varepsilon)}$. Moreover, it was recently shown that getting a $(2 - \varepsilon)$-approximation for general $k$ is NP-hard, assuming the Small Set Expansion Hypothesis. If we use the size of the cut as the parameter, an FPT algorithm to find the exact $k$-Cut is known, but solving the $k$-Cut problem exactly is $W[1]$-hard if we parameterize only by the natural parameter of $k$. An immediate question is: \emph{can we approximate $k$-Cut better in FPT-time, using $k$ as the parameter?} We answer this question positively. We show that for some absolute constant $\varepsilon > 0$, there exists a $(2 - \varepsilon)$-approximation algorithm that runs in time $2^{O(k^6)} \cdot \widetilde{O} (n^4)$. This is the first FPT algorithm that is parameterized only by $k$ and strictly improves the $2$-approximation.


title: Scientific Synergy Between LSST and Euclid

abstract: Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating.


title: Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary

abstract: We solve the XXZ Gaudin model with generic boundary using the modified algebraic Bethe ansatz. The diagonal and triangular cases have been recovered in this general framework. We show that the model for odd or even lengths has two different behaviors. The corresponding Bethe equations are computed for all the cases. For the chain with even length, inhomogeneous Bethe equations are necessary. The higher spin Gaudin models with generic boundary is also treated.


title: Harnessing the Power of Many: Extensible Toolkit for Scalable Ensemble Applications

abstract: Many scientific problems require multiple distinct computational tasks to be executed in order to achieve a desired solution. We introduce the Ensemble Toolkit (EnTK) to address the challenges of scale, diversity and reliability they pose. We describe the design and implementation of EnTK, characterize its performance and integrate it with two distinct exemplar use cases: seismic inversion and adaptive analog ensembles. We perform nine experiments, characterizing EnTK overheads, strong and weak scalability, and the performance of two use case implementations, at scale and on production infrastructures. We show how EnTK meets the following general requirements: (i) implementing dedicated abstractions to support the description and execution of ensemble applications; (ii) support for execution on heterogeneous computing infrastructures; (iii) efficient scalability up to O(10^4) tasks; and (iv) fault tolerance. We discuss novel computational capabilities that EnTK enables and the scientific advantages arising thereof. We propose EnTK as an important addition to the suite of tools in support of production scientific computing.


title: High-Speed Demodulation of weak FBGs Based on Microwave Photonics and Chromatic Dispersion

abstract: A high speed quasi-distributed demodulation method based on the microwave photonics and the chromatic dispersion effect is designed and implemented for weak fiber Bragg gratings (FBGs). Due to the effect of dispersion compensation fiber (DCF), FBG wavelength shift leads to the change of the difference frequency signal at the mixer. With the way of crossing microwave sweep cycle, all wavelengths of cascade FBGs can be high speed obtained by measuring the frequencies change. Moreover, through the introduction of Chirp-Z and Hanning window algorithm, the analysis of difference frequency signal is achieved very well. By adopting the single-peak filter as a reference, the length disturbance of DCF caused by temperature can be also eliminated. Therefore, the accuracy of this novel method is greatly improved, and high speed demodulation of FBGs can easily realize. The feasibility and performance are experimentally demonstrated using 105 FBGs with 0.1% reflectivity, 1 m spatial interval. Results show that each grating can be distinguished well, and the demodulation rate is as high as 40 kHz, the accuracy is about 8 pm.


title: Predicting kinetics of polymorphic transformations from structure mapping and coordination analysis

abstract: To extend rational materials design and discovery into the space of metastable polymorphs, rapid and reliable assessment of their lifetimes is essential. Motivated by the early work of Buerger (1951), here we investigate the routes to predict kinetics of polymorphic transformations using solely crystallographic arguments. As part of this investigation we developed a general algorithm to map crystal structures onto each other and construct transformation pathways between them. The developed algorithm reproduces reliably known transformation pathways and reveals the critical role that the dissociation of chemical bonds along the pathway plays in choosing the best (low-energy barrier) mapping. By utilizing our structure-mapping algorithm we were able to quantitatively demonstrate the intuitive expectation that the minimal dissociation of chemical bonds along the pathway, or in ionic systems, the condition of coordination of atoms along the pathway not decreasing below the coordination in the end compounds, represents the requirement for fast transformation kinetics. We also demonstrate the effectiveness of the structure-mapping algorithm in combination with the coordination analysis in studying transformations between polymorphs for which a detailed atomic-level picture is presently elusive.


title: Estimating the radiative part of QED effects in superheavy nuclear quasimolecules

abstract: A method for calculating the electronic levels in the compact superheavy nuclear quasi-molecules, based on solving the two-center Dirac equation using the multipole expansion of two-center potential, is developed. For the internuclear distances up to $d\sim100$~fm such technique reveals a quite fast convergence and allows for computing the electronic levels in such systems with accuracy $\sim 10^{-6} $. The critical distances $R_{cr}$ between the nuclei for $1\sigma_g$ and $1\sigma_u$ electronic levels in the region $Z\simeq 87-100$ are calculated. By means of the same technique the shifts of electronic levels due to the effective interaction $\Delta U_{AMM}$ of the electron's magnetic anomaly with the Coulomb field of the closely spaced heavy nuclei are evaluated as a function of the internuclear distance and the charge of the nuclei, non-perturbatively both in $Z\alpha$ and (partially) in $\alpha/\pi$. It is shown, that the levels shifts near the lower continuum decrease with the enlarging size of the system of Coulomb sources both in the absolute units and in units of $Z^4 \alpha^5 / \pi n^3$. The last result is generalized to the whole self-energy contribution to the level shifts and so to the possible behavior of radiative part of QED-effects with virtual photon exchange near the lower continuum in the overcritical region.