title: The mass-metallicity relation of zCOSMOS galaxies at z ~ 0.7, its dependence on SFR, and the existence of massive low-metallicity galaxies

abstract: (Abridged) The knowledge of the number and of the physical nature of low-metallicity massive galaxies is crucial for the determination and interpretation of the mass-metallicity relation (MZR). Using VLT-ISAAC near-infrared (NIR) spectroscopy of 39 zCOSMOS z~0.7 galaxies, we have measured Halpha and [NII] emission line fluxes for galaxies with [OII], Hbeta and [OIII] available from VIMOS optical spectroscopy. The NIR spectroscopy enables us to break the degeneracy of the R23 method to derive unambiguously O/H gas metallicities, and also SFRs from extinction corrected Halpha. Using, as a benchmark, the position in the D4000 vs. [OIII]/Hbeta diagram of galaxies with reliable O/Hs from NIR spectroscopy, we were able to break the lower/upper branch R23 degeneracy of additional 900 zCOSMOS z~0.7 galaxies. Additionally, the Halpha-based SFR measurements were used to find the best SFR calibration based on [OII] for the zCOSMOS z~0.7 galaxies without Halpha measurements. We find a fraction of 19% of lower mass 9.5<logM/Msun<10.3 zCOSMOS galaxies which shows a larger evolution of the MZR relation, compared to higher mass galaxies, being more metal poor at a given mass by a factor of 2-3 compared to SDSS. This indicates that the low-mass MZR slope is getting steeper at z~0.7 compared to local galaxies. The existence of these metal-poor galaxies at z~0.7 can be interpreted as the chemical version of galaxy downsizing. Moreover, the sample of zCOSMOS galaxies shows direct evidence that SFR influences the MZR at these redshifts. The comparison of the measured metallicities for the zCOSMOS sample with the values expected for a non-evolving fundamental metallicity relation (FMR) shows broadly agreement, and reveals that also galaxies with lower metallicities and typically higher (specific) SFRs, as found in our zCOSMOS sample at z~0.7, are in agreement with the predictions of a non-evolving Z(M,SFR).


title: Magnetic Fields in High-Mass Infrared Dark Clouds

abstract: High-mass Stars are cosmic engines known to dominate the energetics in the Milky Way and other galaxies. However, their formation is still not well understood. Massive, cold, dense clouds, often appearing as Infrared Dark Clouds (IRDCs), are the nurseries of massive stars. No measurements of magnetic fields in IRDCs in a state prior to the onset of high-mass star formation (HMSF) have previously been available, and prevailing HMSF theories do not consider strong magnetic fields. Here, we report observations of magnetic fields in two of the most massive IRDCs in the Milky Way. We show that IRDCs G11.11-0.12 and G0.253+0.016 are strongly magnetized and that the strong magnetic field is as important as turbulence and gravity for HMSF. The main dense filament in G11.11-0.12 is perpendicular to the magnetic field, while the lower density filament merging onto the main filament is parallel to the magnetic field. The implied magnetic field is strong enough to suppress fragmentation sufficiently to allow HMSF. Other mechanisms reducing fragmentation, such as the entrapment of heating from young stars via high mass surface densities, are not required to facilitate HMSF.


title: The metallicity of galactic winds

abstract: The abundance evolution of galaxies depends critically on the balance between the mixing of metals in their interstellar medium, the inflow of new gas and the outflow of enriched gas. We study these processes in gas columns perpendicular to a galactic disk using sub-parsec resolution simulations that track stellar ejecta with the Flash code. We model a simplified interstellar medium stirred and enriched by supernovae and their progenitors. We vary the density distribution of the gas column and integrate our results over an exponential disk to predict wind and ISM enrichment properties for disk galaxies. We find that winds from more massive galaxies are hotter and more highly enriched, in stark contrast to that which is often assumed in galaxy formation models. We use these findings in a simple model of galactic enrichment evolution, in which the metallicity of forming galaxies is the result of accretion of nearly pristine gas and outflow of enriched gas along an equilibrium sequence. We compare these predictions to the observed mass-metallicity relation, and demonstrate how the galaxy's gas fraction is a key controlling parameter. This explains the observed flattening of the mass-metallicity relation at higher stellar masses.


title: Relativistic simulations of black hole-neutron star coalescence: the jet emerges

abstract: We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of a binary black hole-neutron star on a quasicircular orbit that undergoes merger. The binary mass ratio is 3:1, the black hole initial spin parameter $a/m=0.75$ ($m$ is the black hole Christodoulou mass) aligned with the orbital angular momentum, and the neutron star is an irrotational $\Gamma=2$ polytrope. About two orbits prior to merger (at time $t=t_B$), we seed the neutron star with a dynamically weak interior dipole magnetic field that extends into the stellar exterior. At $t=t_B$ the exterior has a low-density atmosphere with constant plasma parameter $\beta\equiv P_{\rm gas}/P_{\rm mag}$. Varying $\beta$ at $t_B$ in the exterior from $0.1$ to $0.01$, we find that at a time $\sim 4000{\rm M} \sim 100(M_{\rm NS}/1.4M_\odot)$ms [M is the total (ADM) mass] following the onset of accretion of tidally disrupted debris, magnetic winding above the remnant black hole poles builds up the magnetic field sufficiently to launch a mildly relativistic, collimated outflow - an incipient jet. The duration of the accretion and the lifetime of the jet is $\Delta t\sim 0.5(M_{\rm NS}/1.4M_\odot)$s. Our simulations furnish the first explicit examples in GRMHD which show that a jet can emerge following a black hole - neutron star merger.


title: Cross-correlating 21cm intensity maps with Lyman Break Galaxies in the post-reionization era

abstract: We investigate the cross-correlation between the spatial distribution of Lyman Break Galaxies (LBGs) and the 21cm intensity mapping signal at $z\sim[3-5]$. At these redshifts, galactic feedback is supposed to only marginally affect the matter power spectrum, and the neutral hydrogen distribution is independently constrained by quasar spectra. Using a high resolution N-body simulation, populated with neutral hydrogen a posteriori, we forecast for the expected LBG-21cm cross-spectrum and its error for a 21cm field observed by the Square Kilometre Array (SKA1-LOW and SKA1-MID), combined with a spectroscopic LBG survey with the same volume. The cross power can be detected with a signal-to-noise ratio (SNR) up to ~10 times higher (and down to ~4 times smaller scales) than the 21cm auto-spectrum for this set-up, with the SNR depending only very weakly on redshift and the LBG population. We also show that while both the 21cm auto- and LBG-21cm cross-spectra can be reliably recovered after the cleaning of smooth-spectrum foreground contamination, only the cross-power is robust to problematic non-smooth foregrounds like polarized synchrotron emission.


title: The Evolution of Post-Starburst Galaxies from $z\sim1$ to the Present

abstract: Post-starburst galaxies are in the transitional stage between blue, star-forming galaxies and red, quiescent galaxies, and therefore hold important clues for our understanding of galaxy evolution. In this paper, we systematically searched for and identified a large sample of post-starburst galaxies from the spectroscopic dataset of the Sloan Digital Sky Survey (SDSS) Data Release 9. In total, we found more than 6000 objects with redshifts between $z\sim0.05$ and $z\sim1.3$, making this the largest sample of post-starburst galaxies in the literature. We calculated the luminosity function of the post-starburst galaxies using two uniformly selected subsamples: the SDSS Main Galaxy Sample and the Baryon Oscillation Spectroscopic Survey CMASS Sample. The luminosity functions are reasonably fit by half-Gaussian functions. The peak magnitudes shift as a function of redshift from $M\sim-23.5$ at $z\sim0.8$ to $M\sim-20.3$ at $z\sim0.1$. This is consistent with the downsizing trend, whereby more massive galaxies form earlier than low-mass galaxies. We compared the mass of the post-starburst stellar population found in our sample to the decline of the global star-formation rate and found that only a small amount ($\sim1\%$) of all star-formation quenching in the redshift range $z=0.2-0.7$ results in post-starburst galaxies in the luminosity range our sample is sensitive to. Therefore, luminous post-starburst galaxies are not the place where most of the decline in star-formation rate of the universe is happening.


title: WFIRST Ultra-Precise Astrometry II: Asteroseismology

abstract: WFIRST microlensing observations will return high-precision parallaxes, sigma(pi) < 0.3 microarcsec, for the roughly 1 million stars with H<14 in its 2.8 deg^2 field toward the Galactic bulge. Combined with its 40,000 epochs of high precision photometry (~0.7 mmag at H_vega=14 and ~0.1 mmag at H=8), this will yield a wealth of asteroseismic data of giant stars, primarily in the Galactic bulge but including a substantial fraction of disk stars at all Galactocentric radii interior to the Sun. For brighter stars, the astrometric data will yield an external check on the radii derived from the two asteroseismic parameters, <Delta nu> and nu_max, while for the fainter ones, it will enable a mass measurement from the single measurable asteroseismic parameter nu_max. Simulations based on Kepler data indicate that WFIRST will be capable of detecting oscillations in stars from slightly less luminous than the red clump to the tip of the red giant branch, yielding roughly 1 million detections.


title: Linearized iterative least-squares (LIL): A parameter fitting algorithm for component separation in multifrequency CMB experiments such as Planck

abstract: We present an efficient algorithm for the least squares parameter fitting optimized for component separation in multi-frequency CMB experiments. We sidestep some of the problems associated with non-linear optimization by taking advantage of the quasi-linear nature of the foreground model. We demonstrate our algorithm, linearized iterative least-squares (LIL), on the publicly available Planck sky model FFP6 simulations and compare our result with the other algorithms. We work at full Planck resolution and show that degrading the resolution of all channels to that of the lowest frequency channel is not necessary. Finally we present results for the publicly available Planck data. Our algorithm is extremely fast, fitting 6 parameters to 7 lowest Planck channels at full resolution (50 million pixels) in less than 160 CPU-minutes (or few minutes running in parallel on few tens of cores). LIL is therefore easily scalable to future experiments which may have even higher resolution and more frequency channels. We also naturally propagate the uncertainties in different parameters due to noise in the maps as well as degeneracies between the parameters to the final errors on the parameters using Fisher matrix. One indirect application of LIL could be a front-end for Bayesian parameter fitting to find the maximum of the likelihood to be used as the starting point for the Gibbs sampling. We show for rare components, such as the carbon-monoxide emission, present in small fraction of sky, the optimal approach should combine parameter fitting with model selection. LIL may also be useful in other astrophysical applications which satisfy the quasi-linearity criteria.


title: WISE photometry for 400 million SDSS sources

abstract: We present photometry of images from the Wide-Field Infrared Survey Explorer (WISE; Wright et al. 2010) of over 400 million sources detected by the Sloan Digital Sky Survey (SDSS; York et al. 2000). We use a "forced photometry" technique, using measured SDSS source positions, star-galaxy separation and galaxy profiles to define the sources whose fluxes are to be measured in the WISE images. We perform photometry with The Tractor image modeling code, working on our "unWISE" coaddds and taking account of the WISE point-spread function and a noise model. The result is a measurement of the flux of each SDSS source in each WISE band. Many sources have little flux in the WISE bands, so often the measurements we report are consistent with zero. However, for many sources we get three- or four-sigma measurements; these sources would not be reported by the WISE pipeline and will not appear in the WISE catalog, yet they can be highly informative for some scientific questions. In addition, these small-signal measurements can be used in stacking analyses at catalog level. The forced photometry approach has the advantage that we measure a consistent set of sources between SDSS and WISE, taking advantage of the resolution and depth of the SDSS images to interpret the WISE images; objects that are resolved in SDSS but blended together in WISE still have accurate measurements in our photometry. Our results, and the code used to produce them, are publicly available at http://unwise.me.


title: Clumpy Galaxies in CANDELS. I. The Definition of UV Clumps and the Fraction of Clumpy Galaxies at 0.5<z<3

abstract: Although giant clumps of stars are crucial to galaxy formation and evolution, the most basic demographics of clumps are still uncertain, mainly because the definition of clumps has not been thoroughly discussed. In this paper, we study the basic demographics of clumps in star-forming galaxies (SFGs) at 0.5<z<3, using our proposed physical definition that UV-bright clumps are discrete star-forming regions that individually contribute more than 8% of the rest-frame UV light of their galaxies. Clumps defined this way are significantly brighter than the HII regions of nearby large spiral galaxies, either individually or blended, when physical spatial resolution and cosmological dimming are considered. Under this definition, we measure the fraction of SFGs that contain at least one off-center clump (Fclumpy) and the contributions of clumps to the rest-frame UV light and star formation rate of SFGs in the CANDELS/GOODS-S and UDS fields, where our mass-complete sample consists of 3239 galaxies with axial ratio q>0.5. The redshift evolution of Fclumpy changes with the stellar mass (M*) of the galaxies. Low-mass (log(M*/Msun)<9.8) galaxies keep an almost constant Fclumpy of about 60% from z~3.0 to z~0.5. Intermediate-mass and massive galaxies drop their Fclumpy from 55% at z~3.0 to 40% and 15%, respectively, at z~0.5. We find that (1) the trend of disk stabilization predicted by violent disk instability matches the Fclumpy trend of massive galaxies; (2) minor mergers are a viable explanation of the Fclumpy trend of intermediate-mass galaxies at z<1.5, given a realistic observability timescale; and (3) major mergers are unlikely responsible for the Fclumpy trend in all masses at z<1.5. The clump contribution to the rest-frame UV light of SFGs shows a broad peak around galaxies with log(M*/Msun)~10.5 at all redshifts, possibly linked to the molecular gas fraction of the galaxies. (Abridged)


title: Constraining duty cycles through a Bayesian technique

abstract: The duty cycle (DC) of astrophysical sources is generally defined as the fraction of time during which the sources are active. However, DCs are generally not provided with statistical uncertainties, since the standard approach is to perform Monte Carlo bootstrap simulations to evaluate them, which can be quite time consuming for a large sample of sources. As an alternative, considerably less time-consuming approach, we derived the theoretical expectation value for the DC and its error for sources whose state is one of two possible, mutually exclusive states, inactive (off) or flaring (on), as based on a finite set of independent observational data points. Following a Bayesian approach, we derived the analytical expression for the posterior, the conjugated distribution adopted as prior, and the expectation value and variance. We applied our method to the specific case of the inactivity duty cycle (IDC) for supergiant fast X-ray transients. We also studied IDC as a function of the number of observations in the sample. Finally, we compare the results with the theoretical expectations. We found excellent agreement with our findings based on the standard bootstrap method. Our Bayesian treatment can be applied to all sets of independent observations of two-state sources, such as active galactic nuclei, X-ray binaries, etc. In addition to being far less time consuming than bootstrap methods, the additional strength of this approach becomes obvious when considering a well-populated class of sources ($N_{\rm src} \geq 50$) for which the prior can be fully characterized by fitting the distribution of the observed DCs for all sources in the class, so that, through the prior, one can further constrain the DC of a new source by exploiting the information acquired on the DC distribution derived from the other sources. [Abridged]


title: The colors of satellite galaxies in the Illustris Simulation

abstract: Observationally, the fraction of blue satellite galaxies decreases steeply with host halo mass, and their radial distribution around central galaxies is significantly shallower in massive (M_* >10e11M_sun) than in Milky Way like systems. Theoretical models, based primarily on semi-analytical techniques, have had a long-standing problem with reproducing these trends, instead predicting too few blue satellites in general but also estimating a radial distribution that is too shallow, regardless of primary mass. In this Letter, we use the Illustris cosmological simulation to study the properties of satellite galaxies around isolated primaries. For the first time, we find good agreement between theory and observations. We identify the main source of this success relative to earlier work to be a consequence of the large gas contents of satellites at infall, a factor ~5-10 times larger than in semi-analytical models. Because of their relatively large gas reservoirs, satellites can continue to form stars long after infall, with a typical timescale for star-formation to be quenched ~2 Gyr in groups but more than ~5 Gyr for satellites around Milky Way like primaries. The gas contents we infer are consistent with z=0 observations of HI gas in galaxies, although we find large discrepancies among reported values in the literature. A testable prediction of our model is that the gas-to-stellar mass ratio of satellite progenitors should vary only weakly with cosmic time.


title: The Subaru high-z quasar survey: discovery of faint z~6 quasars

abstract: We present the discovery of one or two extremely faint z~6 quasars in 6.5 deg^2 utilizing a unique capability of the wide-field imaging of the Subaru/Suprime-Cam. The quasar selection was made in (i'-z_B) and (z_B-z_R) colors, where z_B and z_R are bandpasses with central wavelengths of 8842A and 9841A, respectively. The color selection can effectively isolate quasars at z~6 from M/L/T dwarfs without the J-band photometry down to z_R<24.0, which is 3.5 mag. deeper than SDSS. We have selected 17 promising quasar candidates. The follow-up spectroscopy for seven targets identified one apparent quasar at z=6.156 with M_1450=-23.10. We also identified one possible quasar at z=6.041 with a faint continuum of M_1450=-22.58 and a narrow Lyman-alpha emission with HWHM=427 km/s, which cannot be distinguished from Lyman-alpha emitters. We derive the quasar luminosity function at z~6 by combining our faint quasar sample with the bright quasar samples by SDSS and CFHQS. Including our data points invokes a higher number density in the faintest bin of the quasar luminosity function than the previous estimate employed. This suggests a steeper faint-end slope than lower-z, though it is yet uncertain based on a small number of spectroscopically identified faint quasars and several quasar candidates are still remain to be diagnosed. The steepening of the quasar luminosity function at the faint-end does increase the expected emission rate of the ionizing photon, however, it only changes by a factor of ~2-6. This was found to be still insufficient for the required photon budget of reionization at z~6.


title: Understanding the Unique Assembly History of Central Group Galaxies

abstract: Central Galaxies (CGs) in massive halos live in unique environments with formation histories closely linked to that of the host halo. In local clusters they have larger sizes ($R_e$) and lower velocity dispersions (sigma) at fixed stellar mass M_star, and much larger R_e at a fixed $\sigma$ than field and satellite galaxies (non-CGs). Using spectroscopic observations of group galaxies selected from the COSMOS survey, we compare the dynamical scaling relations of early-type CGs and non-CGs at z~0.6, to distinguish possible mechanisms that produce the required evolution. CGs are systematically offset towards larger R_e at fixed $\sigma$ compared to non-CGs with similar M_star. The CG R_e-M_star relation also shows differences, primarily driven by a sub-population (~15%) of galaxies with large $R_e$, while the M_star-sigma relations are indistinguishable. These results are accentuated when double Sersic profiles, which better fit light in the outer regions of galaxies, are adopted. They suggest that even group-scale CGs can develop extended components by these redshifts that can increase total $R_e$ and M_star estimates by factors of ~2. To probe the evolutionary link between our sample and cluster CGs, we also analyze two cluster samples at z~0.6 and z~0. We find similar results for the more massive halos at comparable z, but much more distinct CG scaling relations at low-z. Thus, the rapid, late-time accretion of outer components, perhaps via the stripping and accretion of satellites, would appear to be a key feature that distinguishes the evolutionary history of CGs.


title: Eight billion asteroids in the Oort cloud

abstract: The Oort cloud is usually thought of as a collection of icy comets inhabiting the outer reaches of the Solar system, but this picture is incomplete. We use simulations of the formation of the Oort cloud to show that ~4% of the small bodies in the Oort cloud should have formed within 2.5 au of the Sun, and hence be ice-free rock-iron bodies. If we assume these Oort cloud asteroids have the same size distribution as their cometary counterparts, the Large Synoptic Survey Telescope should find roughly a dozen Oort cloud asteroids during ten years of operations. Measurement of the asteroid fraction within the Oort cloud can serve as an excellent test of the Solar system's formation and dynamical history. Oort cloud asteroids could be of particular concern as impact hazards as their high mass density, high impact velocity, and low visibility make them both hard to detect and hard to divert or destroy. However, they should be a rare class of object, and we estimate globally catastrophic collisions should only occur about once per billion years.


title: Maximally Informative Hierarchical Representations of High-Dimensional Data

abstract: We consider a set of probabilistic functions of some input variables as a representation of the inputs. We present bounds on how informative a representation is about input data. We extend these bounds to hierarchical representations so that we can quantify the contribution of each layer towards capturing the information in the original data. The special form of these bounds leads to a simple, bottom-up optimization procedure to construct hierarchical representations that are also maximally informative about the data. This optimization has linear computational complexity and constant sample complexity in the number of variables. These results establish a new approach to unsupervised learning of deep representations that is both principled and practical. We demonstrate the usefulness of the approach on both synthetic and real-world data.


title: Overmassive black holes in the $M_{\rm BH} - \sigma$ diagram do not belong to over (dry) merged galaxies

abstract: Semi-analytical models in a $\Lambda$CDM cosmology have predicted the presence of outlying, "overmassive" black holes at the high-mass end of the (black hole mass -- galaxy velocity dispersion) $M_{\rm BH} - \sigma$ diagram (which we update here with a sample of 89 galaxies). They are a consequence of having experienced more dry mergers -- thought not to increase a galaxy's velocity dispersion -- than the "main-sequence" population. Wet mergers and gas-rich processes, on the other hand, preserve the main correlation. Due to the scouring action of binary supermassive black holes, the extent of these dry mergers (since the last significant wet merger) can be traced by the ratio between the central stellar mass deficit and the black hole mass ($M_{\rm def,*} / M_{\rm BH}$). However, in a sample of 23 galaxies with partially depleted cores, including central cluster galaxies, we show that the "overmassive" black holes are actually hosted by galaxies that appear to have undergone the lowest degree of such merging. In addition, the rotational kinematics of 37 galaxies in the $M_{\rm BH} - \sigma$ diagram reveals that fast and slow rotators are not significantly offset from each other, also contrary to what is expected if these two populations were the product of wet and dry mergers respectively. The observations are thus not in accordance with model predictions and further investigation is required.


title: The instability spectrum of weakly-magnetized SU(2) Reissner-Nordstr\"om black holes

abstract: It is well known that the U(1) Reissner-Nordstr\"om black hole is stable within the framework of the Einstein-Maxwell theory. However, the SU(2) Reissner-Nordstr\"om black-hole solution of the coupled Einstein-Yang-Mills equations is known to be unstable. In fact, this magnetically charged black hole is characterized by an infinite set of unstable (growing in time) perturbation modes. In the present paper we study {\it analytically} the instability resonance spectrum of weakly-magnetized SU(2) Reissner-Nordstr\"om black holes. In particular, we obtain explicit analytical expressions for the infinite set $\{\omega_n\}_{n=0}^{n=\infty}$ of imaginary eigenvalues that characterize the instability growth rates of the perturbation modes. We discuss the role played by these unstable eigenvalues as critical exponents in the gravitational collapse of the Yang-Mills field. Finally, it is shown that our analytical formulas for the characteristic black-hole instability spectrum agree with new numerical data that recently appeared in the literature.


title: Quasi Many-body Localization in Translation Invariant Systems

abstract: It is typically assumed that disorder is essential to realize Anderson localization. Recently, a number of proposals have suggested that an interacting, translation invariant system can also exhibit localization. We examine these claims in the context of a one-dimensional spin ladder. At intermediate time scales, we find slow growth of entanglement entropy consistent with the phenomenology of many-body localization. However, at longer times, all finite wavelength spin polarizations decay in a finite time, independent of system size. We identify a single length scale which parametrically controls both the eventual spin transport times and the divergence of the susceptibility to spin glass ordering. We dub this long pre-thermal dynamical behavior, intermediate between full localization and diffusion, quasi-many body localization.

categories: cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech cond-mat.str-el quant-ph


title: Complexified boost invariance and holographic heavy ion collisions

abstract: At strong coupling holographic studies have shown that heavy ion collisions do not obey normal boost invariance. Here we study a modified boost invariance through a complex shift in time, and show that this leads to surprisingly good agreement with numerical holographic computations. When including perturbations the agreement becomes even better, both in the hydrodynamic and the far-from-equilibrium regime. One of the main advantages is an analytic formulation of the stress-energy tensor of the longitudinal dynamics of holographic heavy ion collisions.


title: Making the Most of the Relic Density for Dark Matter Searches at the LHC 14 TeV Run

abstract: As the LHC continues to search for new weakly interacting particles, it is important to remember that the search is strongly motivated by the existence of dark matter. In view of a possible positive signal, it is essential to ask whether the newly discovered weakly interacting particle can be be assigned the label "dark matter". Within a given set of simplified models and modest working assumptions, we reinterpret the relic abundance bound as a relic abundance range, and compare the parameter space yielding the correct relic abundance with projections of the Run II exclusion regions. Assuming that dark matter is within the reach of the LHC, we also make the comparison with the potential 5$\sigma$ discovery regions. Reversing the logic, relic density calculations can be used to optimize dark matter searches by motivating choices of parameters where the LHC can probe most deeply into the dark matter parameter space. In the event that DM is seen outside of the region giving the correct relic abundance, we will learn that either thermal relic DM is ruled out in that model, or the DM-quark coupling is suppressed relative to the DM coupling strength to other SM particles.


title: Convergence rate, location and $\partial_z^2$ condition for fully bubbling solutions to SU(n+1) Toda Systems

abstract: It is well known that the study of $SU(n+1)$ Toda systems is important not only to Chern-Simons models in Physics, but also to the understanding of holomorphic curves, harmonic sequences or harmonic maps from Riemann surfaces to $\mathbb C\mathbb P^n$. One major goal in the study of $SU(n+1)$ Toda system on Riemann surfaces is to completely understand the asymptotic behavior of fully bubbling solutions. In this article we use a unified approach to study fully bubbling solutions to general $SU(n+1)$ Toda systems and we prove three major sharp estimates important for constructing bubbling solutions: the closeness of blowup solutions to entire solutions, the location of blowup points and a $\partial_z^2$ condition.


title: Ground state entanglement constrains low-energy excitations

abstract: For a general quantum many-body system, we show that its ground-state entanglement imposes a fundamental constraint on the low-energy excitations. For two-dimensional systems, our result implies that any system that supports anyons must have a nonvanishing topological entanglement entropy. We demonstrate the generality of this argument by applying it to three-dimensional quantum many-body systems, and showing that there is a pair of ground state topological invariants that are associated to their physical boundaries. From the pair, one can determine whether the given boundary can or cannot absorb point-like or line-like excitations.


title: The far-infrared/radio correlation and radio spectral index of galaxies in the SFR-M* plane up to z 2

abstract: [Abridged] We study the evolution of the radio spectral index and far-infrared/radio correlation (FRC) across the star-formation rate-stellar masse (i.e. SFR-M*) plane up to z 2. We start from a M*-selected sample of galaxies with reliable SFR and redshift estimates. We then grid the SFR-M* plane in several redshift ranges and measure the infrared luminosity, radio luminosity, radio spectral index, and ultimately the FRC index (i.e. qFIR) of each SFR-M*-z bin. The infrared luminosities of our SFR-M*-z bins are estimated using their stacked far-infrared flux densities inferred from observations obtained with Herschel. Their radio luminosities and radio spectral indices (i.e. alpha, where Snu nu^-alpha) are estimated using their stacked 1.4GHz and 610MHz flux densities from the VLA and GMRT, respectively. Our far-infrared and radio observations include the most widely studied blank extragalactic fields -GOODS-N/S, ECDFS, and COSMOS- covering a sky area of 2deg^2. Using this methodology, we constrain the radio spectral index and FRC index of star-forming galaxies with M*>10^10Msun and 0<z<2.3. We find that alpha^1.4GHz_610MHz does not evolve significantly with redshift or with the distance of a galaxy with respect to the main sequence (MS) of the SFR-M* plane (i.e. Delta_log(SSFR)_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). Instead, star-forming galaxies have a radio spectral index consistent with a canonical value of 0.8, which suggests that their radio spectra are dominated by non-thermal optically thin synchrotron emission. We find that qFIR displays a moderate but statistically significant redshift evolution as qFIR(z)=(2.35+/-0.08)*(1+z)^(-0.12+/-0.04), consistent with some previous literature. Finally, we find no significant correlation between qFIR and Delta_log(SSFR)_MS, though a weak positive trend, as observed in one of our redshift bins, cannot be firmly ruled out using our dataset.


title: The definition of environment and its relation to the quenching of galaxies at z=1-2 in a hierarchical Universe

abstract: A well calibrated method to describe the environment of galaxies at all redshifts is essential for the study of structure formation. Such a calibration should include well understood correlations with halo mass, and the possibility to identify galaxies which dominate their potential well (centrals), and their satellites. Focusing on z = 1 and 2 we propose a method of environmental calibration which can be applied to the next generation of low to medium resolution spectroscopic surveys. Using an up-to-date semi-analytic model of galaxy formation, we measure the local density of galaxies in fixed apertures on different scales. There is a clear correlation of density with halo mass for satellite galaxies, while a significant population of low mass centrals is found at high densities in the neighbourhood of massive haloes. In this case the density simply traces the mass of the most massive halo within the aperture. To identify central and satellite galaxies, we apply an observationally motivated stellar mass rank method which is both highly pure and complete, especially in the more massive haloes where such a division is most meaningful. Finally we examine a test case for the recovery of environmental trends: the passive fraction of galaxies and its dependence on stellar and halo mass for centrals and satellites. With careful calibration, observationally defined quantities do a good job of recovering known trends in the model. This result stands even with reduced redshift accuracy, provided the sample is deep enough to preserve a wide dynamic range of density.


title: Fast Function to Function Regression

abstract: We analyze the problem of regression when both input covariates and output responses are functions from a nonparametric function class. Function to function regression (FFR) covers a large range of interesting applications including time-series prediction problems, and also more general tasks like studying a mapping between two separate types of distributions. However, previous nonparametric estimators for FFR type problems scale badly computationally with the number of input/output pairs in a data-set. Given the complexity of a mapping between general functions it may be necessary to consider large data-sets in order to achieve a low estimation risk. To address this issue, we develop a novel scalable nonparametric estimator, the Triple-Basis Estimator (3BE), which is capable of operating over datasets with many instances. To the best of our knowledge, the 3BE is the first nonparametric FFR estimator that can scale to massive datasets. We analyze the 3BE's risk and derive an upperbound rate. Furthermore, we show an improvement of several orders of magnitude in terms of prediction speed and a reduction in error over previous estimators in various real-world data-sets.


title: Weak-value amplification as an optimal metrological protocol

abstract: The implementation of weak-value amplification requires the pre- and post-selection of states of a quantum system, followed by the observation of the response of the meter, which interacts weakly with the system. Data acquisition from the meter is conditioned to successful post-selection events. Here we derive an optimal post-selection procedure for estimating the coupling constant between system and meter, and show that it leads both to weak-value amplification and to the saturation of the quantum Fisher information, under conditions fulfilled by all previously reported experiments on the amplification of weak signals. For most of the pre-selected states, full information on the coupling constant can be extracted from the meter data set alone, while for a small fraction of the space of pre-selected states, it must be obtained from the post-selection statistics.


title: The level four braid group

abstract: By evaluating the Burau representation at t=-1, we obtain a symplectic representation of the braid group. We define the congruence subgroups of the braid group to be the preimages of the principal congruence subgroups of the symplectic group. Our main result is that the level four congruence subgroup of the braid group is equal to the group generated by squares of Dehn twists and is also equal to the group generated by squares of pure braids. We also compute the image of the point pushing subgroup under the symplectic representation.


title: Framed correspondences and the Milnor-Witt K-theory

abstract: The article is to construct a graded ring isomorphism between $H_0(ZF(\Delta^{\bullet}_k,\mathbb{G}_m^{\wedge *}))$ and the Milnor-Witt K-theory ring $K^{MW}_{*\geqslant 0}(k)$, where $k$ is a field of characteristic zero and $ZF_*(k)$ is the category of linear framed correspondences of algebraic k-varieties, introduced by Garkusha and Panin. As it was shown by Garkusha and Panin, this partially recovers the computation of the motivic cohomotopy groups ${\pi}^{n,n}(\Sigma^{\infty}_{S^1}\Sigma^{\infty}_{\mathbb{G}_m}S^0)(k)$, which is originally due to Morel.


title: Detection of discretized single-shell penetration in mesoscopic vortex matter

abstract: We investigated configurational changes in mesoscopic vortex matter with less than thousand vortices during flux penetration in freestanding 50 $\mu$m diameter disks of Bi$_{2}$Sr$_{2}$CaCu$_{2}$O$_{8 +\delta}$. High-resolution AC and DC local magnetometry data reveal oscillations in the transmittivity echoed in peaks in the third-harmonics magnetic signal fainting on increasing vortex density. By means of extra experimental evidence and a simple geometrical analysis we show that these features fingerprint the discretized entrance of single-shells of vortices having a shape that mimics the sample edge.


title: Cohomology classes of interval positroid varieties and a conjecture of Liu

abstract: To each finite subset of $\mathbb{Z}^2$ (a diagram), one can associate a subvariety of a complex Grassmannian (a diagram variety), and a representation of a symmetric group (a Specht module). Liu has conjectured that the cohomology class of a diagram variety is represented by the Frobenius characteristic of the corresponding Specht module. We give a counterexample to this conjecture. However, we show that for the diagram variety of a permutation diagram, Liu's conjectured cohomology class $\sigma$ is at least an upper bound on the actual class $\tau$, in the sense that $\sigma - \tau$ is a nonnegative linear combination of Schubert classes. To do this, we exhibit the appropriate diagram variety as a component in a degeneration of one of Knutson's interval positroid varieties (up to Grassmann duality). A priori, the cohomology classes of these interval positroid varieties are represented by affine Stanley symmetric functions. We give a different formula for these classes as ordinary Stanley symmetric functions, one with the advantage of being Schur-positive and compatible with inclusions between Grassmannians.


title: Quantum Limit in a Magnetic Field for Triplet Superconductivity in a Quasi-One-Dimensional Conductor

abstract: We theoretically consider the upper critical magnetic field, perpendicular to a conducting axis in a triplet quasi-one-dimensional superconductor. In particular, we demonstrate that, at high magnetic fields, the orbital effects against superconductivity in a magnetic field are reversible and, therefore, superconductivity can restore. It is important that the above mentioned quantum limit can be achieved in presumably triplet quasi-one-dimensional superconductor Li$_{0.9}$Mo$_6$O$_{17}$ [J.-F. Mercure et al., Phys. Rev. Lett. {\bf 108}, 187003 (2012)] at laboratory available pulsed magnetic fields of the order of $H = 500-700 \ T$.


title: Taming the pion condensation in QCD at finite baryon density

abstract: In the Monte Carlo study of QCD at finite baryon density based upon the phase reweighting method, the pion condensation in the phase-quenched theory and associated zero-mode prevent us to go to the low-temperature high-density region. We propose a method to circumvent them by a simple modification of the density of state method. We first argue that the standard version of the density of state method, which is invented to solve the overlapping problem, is effective only for a certain `good' class of observables. We then modify it so as to solve the overlap problem for `bad' observables as well. While, in the standard version of the density of state method, we usually constrain an observable we are interested in, we fix a different observable in our new method which has a sharp peak at some particular value characterizing the correct vacuum of the target theory. In the finite-density QCD, such an observable is the pion condensate. The average phase becomes vanishingly small as the value of the pion condensate becomes large, hence it is enough to consider configurations with small values of pion condensate, where the zero mode does not appear. We demonstrate an effectiveness of our method by using a toy model (the chiral random matrix theory) which captures the properties of finite-density QCD qualitatively. We also argue how to apply our method to other theories including finite-density QCD. Although the example we study numerically is based on the phase reweighting method, the same idea can be applied to more general reweighting methods and we show how this idea can be applied to find a possible QCD critical point.


title: Pulsar Observations Using the First Station of the Long Wavelength Array and the LWA Pulsar Data Archive

abstract: We present initial pulsar results from the first station of the Long Wavelength Array (LWA1) obtained during the commissioning period of LWA1 and early science results. We present detections of periodic emission from 44 previously known pulsars, including 3 millisecond pulsars (MSPs). The effects of the interstellar medium on pulsar emission are significantly enhanced at the low frequencies of the LWA1 band (10--88 MHz), making LWA1 a very sensitive instrument for characterizing changes in dispersion measures (DM) and other effects from the interstellar medium. Pulsars also often have significant evolution in their pulse profile at low frequency and a break in their spectral index. We report DM measurements for 44 pulsars, mean flux density measurements for 36 pulsars, and multi-frequency component spacing and widths for 15 pulsars with more than one profile component. For 27 pulsars, we report spectral index measurements within our frequency range. We also introduce the LWA1 Pulsar Data Archive, which stores reduced data products from LWA1 pulsar observations. Reduced data products for the observations presented here can be found on the archive. Reduced data products from future LWA1 pulsar observations will also be made available through the archive.


title: Packing odd $T$-joins with at most two terminals

abstract: Take a graph $G$, an edge subset $\Sigma\subseteq E(G)$, and a set of terminals $T\subseteq V(G)$ where $|T|$ is even. The triple $(G,\Sigma,T)$ is called a signed graft. A $T$-join is odd if it contains an odd number of edges from $\Sigma$. Let $\nu$ be the maximum number of edge-disjoint odd $T$-joins. A signature is a set of the form $\Sigma\triangle \delta(U)$ where $U\subseteq V(G)$ and $|U\cap T)$ is even. Let $\tau$ be the minimum cardinality a $T$-cut or a signature can achieve. Then $\nu\leq \tau$ and we say that $(G,\Sigma,T)$ packs if equality holds here. We prove that $(G,\Sigma,T)$ packs if the signed graft is Eulerian and it excludes two special non-packing minors. Our result confirms the Cycling Conjecture for the class of clutters of odd $T$-joins with at most two terminals. Corollaries of this result include, the characterizations of weakly and evenly bipartite graphs, packing two-commodity paths, packing $T$-joins with at most four terminals, and a new result on covering edges with cuts.


title: Breakdown of statistical inference from some random experiments

abstract: Many experiments can be interpreted in terms of random processes operating according to some internal protocols. When experiments are costly or cannot be repeated only one or a few finite samples are available. In this paper we study data generated by pseudo-random computer experiments operating according to particular internal protocols. We show that the standard statistical analysis performed on a sample, containing 100000 data points or more, may sometimes be highly misleading and statistical errors largely underestimated. Our results confirm in a dramatic way the dangers of standard asymptotic statistical inference if a sample is not homogenous. We demonstrate that analyzing various subdivisions of samples by multiple chi-square tests and chi-square frequency graphs is very effective in detecting sample inhomogeneity. Therefore to assure correctness of the statistical inference the above mentioned chi-square tests and other non-parametric sample homogeneity tests should be incorporated in any statistical analysis of experimental data. If such tests are not performed the reported conclusions and estimates of the errors cannot be trusted.


title: Limit behaviour of $\mu-$equicontinuous cellular automata

abstract: The concept of $\mu-$equicontinuity was introduced by Gilman to classify cellular automata. We show that under some conditions the sequence of Cesaro averages of a measure $\mu,$ converge under the actions of a $\mu -$equicontinuous CA. We address questions raised by Blanchard-Tisseur on whether the limit measure is either shift-ergodic, a uniform Bernoulli measure or ergodic with respect to the CA. Many of our results hold for CA on multidimensional subshifts.


title: Superfluid (Amplitude) Fluctuations Above $T_c$ in a Unitary Fermi Gas

abstract: We study the transport properties of a Fermi gas with strong attractive interactions close to the unitary limit. In particular, we compute the spin diffusion lifetime of the Fermi gas due to superfluid fluctuations above the BCS transition temperature $T_c$. To calculate the spin diffusion lifetime we need the scattering amplitudes. The scattering amplitudes are dominated by the superfluid fluctuations at temperatures just above $T_c$. The normal scattering amplitudes are calculated from the Landau parameters. These Landau parameters are obtained from the local version of the induced interaction model for computing Landau parameters. We also calculate the leading order finite temperature correction to the diffusion lifetime. A calculation of the spin diffusion coefficient is presented in the end. Upon choosing a proper value of $F_0^a$, we are able to present a good match between the theoretical result and the experimental measurement which indicates the presence of the superfluid fluctuations near $T_c$.


title: Multiple solutions for Schrodinger equations under the partially sublinear case

abstract: In this paper, we prove the infinitely many solutions to a class of sublinear Schr\"{o}dinger-Poisson equations by using an extension of Clark's theorem established by Zhaoli Liu and Zhi-Qiang Wang.


title: Tidal evolution in multiple planet systems: application to Kepler-62 and Kepler-186

abstract: A large number of observed exoplanets are part of multiple planet systems. Most of these systems are sufficiently close-in to be tidally evolving. In such systems, there is a competition between the excitation caused by planet-planet interactions and tidal damping. Using as an example two multiple planet systems, which host planets in the surface liquid water habitable zone (HZ): Kepler-62 and Kepler-186, we show the importance and effect of both planetary and stellar tides on the dynamical evolution of planets and on the climate of the HZ planets.


title: Higher-order MRFs based image super resolution: why not MAP?

abstract: A trainable filter-based higher-order Markov Random Fields (MRFs) model - the so called Fields of Experts (FoE), has proved a highly effective image prior model for many classic image restoration problems. Generally, two options are available to incorporate the learned FoE prior in the inference procedure: (1) sampling-based minimum mean square error (MMSE) estimate, and (2) energy minimization-based maximum a posteriori (MAP) estimate. This letter is devoted to the FoE prior based single image super resolution (SR) problem, and we suggest to make use of the MAP estimate for inference based on two facts: (I) It is well-known that the MAP inference has a remarkable advantage of high computational efficiency, while the sampling-based MMSE estimate is very time consuming. (II) Practical SR experiment results demonstrate that the MAP estimate works equally well compared to the MMSE estimate with exactly the same FoE prior model. Moreover, it can lead to even further improvements by incorporating our discriminatively trained FoE prior model. In summary, we hold that for higher-order natural image prior based SR problem, it is better to employ the MAP estimate for inference.


title: Hysteretic Magnetotransport in SmB6 at Low Magnetic Fields

abstract: Utilizing Corbino disc structures, we have examined the magnetic field response of resistivity for the surface states of SmB6 on different crystalline surfaces at low temperatures. Our results reveal a hysteretic behavior whose magnitude depends on the magnetic field sweep rate and temperature. Although this feature becomes smaller when the field sweep is slower, a complete elimination or saturation is not observed in our slowest sweep-rate measurements, which is much slower than a typical magnetotransport trace. These observations cannot be explained by quantum interference corrections such as weak anti-localization. Instead, they are consistent with behaviors of glassy surface magnetic ordering, whose magnetic origin is most likely from samarium oxide (Sm2O3) forming on the surface during exposure to ambient conditions.


title: Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates

abstract: We extend the positivity-preserving method of Zhang & Shu (2010, JCP, 229, 3091-3120) to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for non-equilibrium transport problems in general relativity and also in science and engineering applications with specific geometries. The method achieves high-order accuracy using Discontinuous Galerkin (DG) discretization of phase space and strong stability-preserving, Runge-Kutta (SSP-RK) time integration. Special care in taken to ensure that the method preserves strict bounds for the phase space distribution function $f$; i.e., $f\in[0,1]$. The combination of suitable CFL conditions and the use of the high-order limiter proposed in Zhang & Shu (2010) is sufficient to ensure positivity of the distribution function. However, to ensure that the distribution function satisfies the upper bound, the discretization must, in addition, preserve the divergence-free property of the phase space flow. Proofs that highlight the necessary conditions are presented for general curvilinear coordinates, and the details of these conditions are worked out for some commonly used coordinate systems (i.e., spherical polar spatial coordinates in spherical symmetry and cylindrical spatial coordinates in axial symmetry, both with spherical momentum coordinates). Results from numerical experiments --- including one example in spherical symmetry adopting the Schwarzschild metric --- demonstrate that the method achieves high-order accuracy and that the distribution function satisfies the maximum principle.


title: Optical Multi-Channel Intensity Interferometry - or: How To Resolve O-Stars in the Magellanic Clouds

abstract: Intensity interferometry, based on the Hanbury Brown-Twiss effect, is a simple and inexpensive method for optical interferometry at microarcsecond angular resolutions; its use in astronomy was abandoned in the 1970s because of low sensitivity. Motivated by recent technical developments, we argue that the sensitivity of large modern intensity interferometers can be improved by factors up to approximately 25,000, corresponding to 11 photometric magnitudes, compared to the pioneering Narrabri Stellar Interferometer. This is made possible by (i) using avalanche photodiodes (APD) as light detectors, (ii) distributing the light received from the source over multiple independent spectral channels, and (iii) use of arrays composed of multiple large light collectors. Our approach permits the construction of large (with baselines ranging from few kilometers to intercontinental distances) optical interferometers at the cost of (very) long-baseline radio interferometers. Realistic intensity interferometer designs are able to achieve limiting R-band magnitudes as good as ~14, sufficient for spatially resolved observations of main-sequence O-type stars in the Magellanic Clouds. Multi-channel intensity interferometers can address a wide variety of science cases: (i) linear radii, effective temperatures, and luminosities of stars; (ii) mass-radius relationships of compact stellar remnants; (iii) stellar rotation; (iv) stellar convection and the interaction of stellar photospheres and magnetic fields; (v) the structure and evolution of multiple stars; (vi) direct measurements of interstellar distances; (vii) the physics of gas accretion onto supermassive black holes; and (viii) calibration of amplitude interferometers by providing a sample of calibrator stars.


title: On the rigid string contribution to the interquark potential

abstract: Owing to Lorentz invariance, the leading terms in the effective action describing the low-energy dynamics of flux tubes in a confining gauge theory are universal. Besides the terms corresponding to the Nambu-Goto action, Lorentz invariance also allows terms describing a "rigid string", proposed long ago by Kleinert and Polyakov as a model for confinement in a dual superconductor scenario. In this contribution we discuss some of the non-trivial properties of this type of string and evaluate analytically the corrections to the interquark potential induced by extrinsic curvature terms appearing in the rigid string action.


title: Alpha clusters and collective flow in ultra-relativistic carbon - heavy nucleus collisions

abstract: We investigate ultrarelativistic collisions of the 12C nucleus with heavy targets and show that the harmonic flow measures based on ratios of cumulant moments are particularly suited to study the intrinsic deformation of the 12C wave function. That way one can probe the expected alpha clusterization in the ground state, which leads to large initial triangularity in the shape of the fireball in the transverse plane. We show that the clusterization effect results in very characteristic behavior of the ratios of the cumulant moments as functions of the number of participant nucleons, both for the elliptic and triangular deformations. Thus the experimental event-by-event studies of harmonic flow in ultrarelativistic light-heavy collisions may offer a new window to look at the ground-state structure of light nuclei.


title: On the Expansion Rate, Age, and Distance of the Supernova Remnant G266.2-1.2 (Vela Jr.)

abstract: An analysis of Chandra ACIS data for two relatively bright and narrow portions of the northwestern rim of G266.2-1.2 (a.k.a. RX J0852.0-4622 or Vela Jr.) reveal evidence of a radial displacement of 2.40 +/- 0.56 arcsec between 2003 and 2008. The corresponding expansion rate (0.42 +/- 0.10 arcsec/yr or 13.6 +/- 4.2%/kyr) is about half the rate reported for an analysis of XMM-Newton data from a similar, but not identical, portion of the rim over a similar, but not identical, time interval (0.84 +/- 0.23 arcsec/yr, Katsuda et al. 2008a). If the Chandra rate is representative of the remnant as a whole, then the results of a hydrodynamic analysis suggest that G266.2-1.2 is between 2.4 and 5.1 kyr old if it is expanding into a uniform ambient medium (whether or not it was produced by a Type Ia or Type II event). If the remnant is expanding into the material shed by a steady stellar wind, then the age could be as much as 50% higher. The Chandra expansion rate and a requirement that the shock speed be greater than or equal to 1000 km/s yields a lower limit on the distance of 0.5 kpc. An analysis of previously-published distance estimates and constraints suggests G266.2-1.2 is no further than 1.0 kpc. This range of distances is consistent with the distance to the nearer of two groups of material in the Vela Molecular Ridge (0.7 +/- 0.2 kpc, Liseau et al. 1992) and to the Vel OB1 association (0.8 kpc, Eggen 1982).


title: Effective string description of the interquark potential in the 3D U(1) lattice gauge theory

abstract: The U(1) lattice gauge theory in three dimensions is a perfect laboratory to study the properties of the confining string. On the one hand, thanks to the mapping to a Coulomb gas of monopoles, the confining properties of the model can be studied semi-classically. On the other hand, high-precision numerical estimates of Polyakov loop correlators can be obtained via a duality map to a spin model. This allowed us to perform high-precision tests of the universal behavior of the effective string and to find macroscopic deviations with respect to the expected Nambu-Goto predictions. These corrections could be fitted with very good precision including a contribution (which is consistent with Lorentz symmetry) proportional to the square of the extrinsic curvature in the effective string action, as originally suggested by Polyakov. Performing our analysis at different values of $\beta$ we were able to show that this term scales as expected by Polyakov's solution and dominates in the continuum. We also discuss the interplay between the extrinsic curvature contribution and the boundary correction induced by the Polyakov loops.


title: Pairwise balanced designs covered by bounded flats

abstract: We prove that for any $K$ and $d$, there exist, for all sufficiently large admissible $v$, a pairwise balanced design PBD$(v,K)$ of dimension $d$ for which all $d$-point-generated flats are bounded by a constant independent of $v$. We also tighten a prior upper bound for $K = \{3,4,5\}$, in which case there are no divisibility restrictions on the number of points. One consequence of this latter result is the construction of latin squares `covered' by small subsquares.


title: The rise and fall of a challenger: the Bullet Cluster in $\Lambda$ Cold Dark Matter simulations

abstract: The Bullet Cluster has provided some of the best evidence for the $\Lambda$ cold dark matter ($\Lambda\rm{CDM}$) model via direct empirical proof of the existence of collisionless dark matter, while posing a serious challenge owing to the unusually high inferred pairwise velocities of its progenitor clusters. Here we investigate the probability of finding such a high-velocity pair in large-volume N-body simulations, particularly focusing on differences between halo finding algorithms. We find that algorithms that do not account for the kinematics of infalling groups yield vastly different statistics and probabilities. When employing the ROCKSTAR halo finder that considers particle velocities, we find numerous Bullet-like pair candidates that closely match not only the high pairwise velocity, but also the mass, mass ratio, separation distance, and collision angle of the initial conditions that have been shown to produce the Bullet Cluster in non-cosmological hydrodynamic simulations. The probability of finding a high pairwise velocity pair among haloes with $M_{\rm halo}\geq10^{14} M_{\odot}$ is $4.6\times 10^{-4}$ using ROCKSTAR, while it is $\approx 34\times$ lower using a friends-of-friends (FOF) based approach as in previous studies. This is because the typical spatial extent of Bullet progenitors is such that FOF tends to group them into a single halo despite clearly distinct kinematics. Further requiring an appropriately high average mass among the two progenitors, we find the comoving number density of potential Bullet-like candidates to be on the order of $\approx10^{-10} \,{\rm Mpc}^{-3}$. Our findings suggest that $\Lambda\rm{CDM}$ straightforwardly produces massive, high relative velocity halo pairs analogous to Bullet Cluster progenitors, and hence the Bullet Cluster does not present a challenge to the $\Lambda\rm{CDM}$ model.


title: The connection between the UV colour of early type galaxies and the stellar initial mass function revisited

abstract: We extend our initial study of the connection between the UV colour of galaxies and both the inferred stellar mass-to-light ratio, $\Upsilon_*$, and a mass-to-light ratio referenced to Salpeter initial mass function (IMF) models of the same age and metallicity, $\Upsilon_*/\Upsilon_{Sal}$, using new UV magnitude measurements for a much larger sample of early-type galaxies, ETGs, with dynamically determined mass-to-light ratios. We confirm the principal empirical finding of our first study, a strong correlation between the GALEX FUV-NUV colour and $\Upsilon_*$. We show that this finding is not the result of spectral distortions limited to a single passband (eg. metallicity-dependent line-blanketing in the NUV band), or of the analysis methodology used to measure $\Upsilon_*$, or of the inclusion or exclusion of the correction for stellar population effects as accounted for using $\Upsilon_*/\Upsilon_{Sal}$. The sense of the correlation is that galaxies with larger $\Upsilon_*$, or larger $\Upsilon_*/\Upsilon_{Sal}$, are bluer in the UV. We conjecture that differences in the low mass end of the stellar initial mass function, IMF, are related to the nature of the extreme horizontal branch stars generally responsible for the UV flux in ETGs. If so, then UV color can be used to identify ETGs with particular IMF properties and to estimate $\Upsilon_*$. We also demonstrate that UV colour can be used to decrease the scatter about the Fundamental Plane and Manifold, and to select peculiar galaxies for follow-up with which to further explore the cause of variations in $\Upsilon_*$ and UV colour.


title: Constant terms of Eisenstein series over a totally real field

abstract: In this paper, we compute constant terms of Eisenstein series defined over a totally real field, at various cusps. In his paper published in 2003, M. Ohta computed the constant terms of Eisenstein series of weight two over the field of rational numbers, at all equivalence classes of cusps. As for Eisenstein series defined over a totally real field, S. Dasgupta, H. Darmon and R. Pollack calculated the constant terms at particular (not all) equivalence classes of cusps in 2011. We compute constant terms of Eisenstein series defined over a general totally real field at all equivalence classes of cusps, and describe them explicitly in terms of Hecke $L$-functions. This investigation is motivated by M. Ohta's work on congruence modules related to Eisenstein series defined over the field of rational numbers.


title: Diagonality and idempotents with applications to problems in operator theory and frame theory

abstract: We prove that a nonzero idempotent is zero-diagonal if and only if it is not a Hilbert-Schmidt perturbation of a projection, along with other useful equivalences. Zero-diagonal operators are those whose diagonal entries are identically zero in some basis. We also prove that any bounded sequence appears as the diagonal of some idempotent operator, thereby providing a characterization of inner products of dual frame pairs in infinite dimensions. Furthermore, we show that any absolutely summable sequence whose sum is a positive integer appears as the diagonal of a finite rank idempotent.


title: Hamiltonian anomalies from extended field theories

abstract: We develop a proposal by Freed to see anomalous field theories as relative field theories, namely field theories taking value in a field theory in one dimension higher, the anomaly field theory. We show that when the anomaly field theory is extended down to codimension 2, familiar facts about Hamiltonian anomalies can be naturally recovered, such as the fact that the anomalous symmetry group admits only a projective representation on the state space, or that the latter is really an abelian gerbe rather than an ordinary Hilbert space. We include in the discussion the case of non-invertible anomaly field theories, which is relevant to six-dimensional (2,0) superconformal theories. In this case, we show that the Hamiltonian anomaly is characterized by a degree 2 non-abelian group cohomology class, associated to the non-abelian gerbe playing the role of the state space of the anomalous theory. We construct Dai-Freed theories, governing the anomalies of chiral fermionic theories, and Wess-Zumino theories, governing the anomalies of Wess-Zumino terms, as extended field theories down to codimension 2.


title: Spectral functions of charmonium with 2+1 flavours of dynamical quarks

abstract: Finite temperature charmonium spectral functions in the pseudoscalar(PS) and vector(V) channels are studied in lattice QCD with 2+1 flavours of dynamical Wilson quarks, on fine isotropic lattices (with a lattice spacing of 0.057fm), with a non-physical pion mass of 545MeV. The highest temperature studied is approximately 1.4Tc. Up to this temperature no significant variation of the spectral function is seen in the PS channel. The V channel shows some temperature dependence, which seems to be consistent with a temperature dependent low frequency peak related to heavy quark transport, plus a temperature independent term at omega > 0. These results are in accord with previous calculations using the quenched approximation.


title: Dynamical Evolution of the Earth-Moon Progenitors - Whence Theia?

abstract: We present integrations of a model Solar System with five terrestrial planets (beginning ~30-50 Myr after the formation of primitive Solar System bodies) in order to determine the preferred regions of parameter space leading to a giant impact that resulted in the formation of the Moon. Our results indicate which choices of semimajor axes and eccentricities for Theia (the proto-Moon) at this epoch can produce a late Giant Impact, assuming that Mercury, Venus, and Mars are near the current orbits. We find that the likely semimajor axis of Theia, at the epoch when our simulations begin, depends on the assumed mass ratio of Earth-Moon progenitors (8/1, 4/1, or 1/1). The low eccentricities of the terrestrial planets are most commonly produced when the progenitors have similar semimajor axes at the epoch when our integrations commence. Additionally, we show that mean motion resonances among the terrestrial planets and perturbations from the giant planets can affect the dynamical evolution of the system leading to a late Giant Impact.


title: Physical Galaxy Pairs and Their Effects on Star Formation

abstract: We present 776 truly physical galaxy pairs, 569 of them are close pairs and 208 false pairs from Karachentsev (1972) and Reduzzi & Rampazzo (1995) catalogues, which contains 1012 galaxy pairs. Also we carried out star formation activity through the far-infrared emission (FIR) in physical (truly) interacting galaxies in some galaxy pairs and compared them with projection (optical) interacting galaxy pairs. We focused on the triggering of star formation by interactions and analyzed the enhancement of star formation activity in terms of truly physical galaxy pairs. The large fraction of star formation activity is probably due to the activity in the exchange of matter between the truly companions. The star formation rate (SFR) of galaxies in truly galaxy pairs is found to be more enhanced than the apparent pairs.


title: On Weak Hamiltonicity of a Random Hypergraph

abstract: A {\it weak (Berge) cycle} is an alternating sequence of vertices and (hyper)edges $C=(v_0, e_1, v_1, ..., v_{\ell-1}, e_\ell, v_{\ell}=v_0)$ such that the vertices $v_0, ..., v_{\ell-1}$ are distinct with $v_k, v_{k+1} \in e_{k}$ for each $k$, but the edges $e_1, ..., e_\ell$ are not necessarily distinct. We prove that the main barrier to the random $d$-uniform hypergraph $H_d(n,p),$ where each of the potential edges of cardinality $d$ is present with probability $p$, developing a weak Hamilton cycle is the presence of isolated vertices. In particular, for $d \geq 3$ fixed and $p=(d-1)! \frac{\ln n + c}{n^{d-1}}$, the probability that $H_d(n, p)$ has a weak Hamilton cycle tends to $e^{-e^{-c}}$, which is also the limiting probability that $H_d(n,p)$ has no isolated vertices. As a consequence, the probability that the random hypergraph $H_d(n, m=\frac{n(\ln n + c)}{d}),$ where $m$ potential edges are chosen uniformly at random to be present, is weak Hamiltonian also tends to $e^{-e^{-c}}$.


title: Phase Diagram of Fractional Quantum Hall Effect of Composite Fermions in Multi-Component Systems

abstract: While the integer quantum Hall effect of composite fermions manifests as the prominent fractional quantum Hall effect (FQHE) of electrons, the FQHE of composite fermions produces further, more delicate states, arising from a weak residual interaction between composite fermions. We study the spin phase diagram of these states, motivated by the recent experimental observation by Liu {\em et al.} \cite{Liu14a,Liu14b} of several spin-polarization transitions at 4/5, 5/7, 6/5, 9/7, 7/9, 8/11 and 10/13 in GaAs systems. We show that the FQHE of composite fermions is much more prevalent in multicomponent systems, and consider the feasibility of such states for systems with ${\cal N}$ components for an SU(${\cal N}$) symmetric interaction. Our results apply to GaAs quantum wells, wherein electrons have two components, to AlAs quantum wells and graphene, wherein electrons have four components (two spins and two valleys), and to an H-terminated Si(111) surface, which can have six components. The aim of this article is to provide a fairly comprehensive list of possible incompressible fractional quantum Hall states of composite fermions, their SU(${\cal N}$) spin content, their energies, and their phase diagram as a function of the generalized "Zeeman" energy. We obtain results at three levels of approximation: from ground state wave functions of the composite fermion theory, from composite fermion diagonalization, and, whenever possible, from exact diagonalization. Effects of finite quantum well thickness and Landau level mixing are neglected in this study. We compare our theoretical results with the experiments of Liu {\em et al.} \cite{Liu14a,Liu14b} as well as of Yeh {\em et al.} \cite{Yeh99} for a two component system.


title: Synchronization of Heterogeneous Kuramoto Oscillators with Arbitrary Topology

abstract: We study synchronization of coupled Kuramoto oscillators with heterogeneous inherent frequencies and general underlying connectivity. We provide conditions on the coupling strength and the initial phases which guarantee the existence of a Positively Invariant Set (PIS) and lead to synchronization. Unlike previous works that focus only on analytical bounds, here we introduce an optimization approach to provide a computational-analytical bound that can further exploit the particular features of each individual system such as topology and frequency distribution. Examples are provided to illustrate our results as well as the improvement over previous existing bounds.


title: Square-free values of multivariate polynomials over function fields in linear sparse sets

abstract: Let f be a square-free polynomial in Fq[t][x] where Fq is a field of q elements. We view f as a polynomial in the variable x with coefficients in the ring Fq[t]. We study squarefree values of f in sparse subsets of Fq[t] which are given by a linear condition. The motivation for our study is an analogue problem of representing square-free integers by integer polynomials, where it is conjectured that setting aside some simple exceptional cases, a square-free polynomial f in Z[x] takes infinitely many square-free values. Let c(t) be a polynomial in Fq[t] of degree less than m, and let k < m be coprime to q. A consequence of the main result we show, is that if q is sufficiently large with respect to m and the degrees of f in t and x, then there exist $\beta_1,\beta_2$ in Fq such that $f(t,c(t)+\beta_1t^k+\beta_2)$ is square-free. Moreover, as q tends to infinity, the last is true for almost all $\beta_1$ and $\beta_2$ in Fq. The main result shows that a similar result holds also for other cases. We then generalize the results to multivariate polynomials.


title: Frenetic aspects of second order response

abstract: Starting from second order around thermal equilibrium, the response of a statistical mechanical system to an external stimulus is not only governed by dissipation and depends explicitly on dynamical details of the system. The so called frenetic contribution in second order around equilibrium is illustrated in different physical examples, such as for non-thermodynamic aspects in the coupling between system and reservoir, for the dependence on disorder in dielectric response and for the nonlinear correction to the Sutherland--Einstein relation. More generally, the way in which a system's dynamical activity changes by the pertubation is visible (only) from nonlinear response.


title: Emergence of chaos in interacting communities

abstract: We introduce a simple dynamical model of two interacting communities whose elements are subject to stochastic discrete-time updates governed by only bilinear interactions. When the intra- and inter-couplings are cooperative, the two communities reach asymptotically an equilibrium state. However, when the intra- or inter-couplings are anti-cooperative, the system may remain in perpetual oscillations and, when the coupling values belong to certain intervals, two possible scenarios arise, characterized either by erratic aperiodic trajectories and high sensitiveness to small changes of the couplings, or by chaotic trajectories and bifurcation cascades. Quite interestingly, we find out that even a moderate consensus in one single community can remove the chaos. Connections of the model with interacting stock markets are discussed.


title: Consensus Message Passing for Layered Graphical Models

abstract: Generative models provide a powerful framework for probabilistic reasoning. However, in many domains their use has been hampered by the practical difficulties of inference. This is particularly the case in computer vision, where models of the imaging process tend to be large, loopy and layered. For this reason bottom-up conditional models have traditionally dominated in such domains. We find that widely-used, general-purpose message passing inference algorithms such as Expectation Propagation (EP) and Variational Message Passing (VMP) fail on the simplest of vision models. With these models in mind, we introduce a modification to message passing that learns to exploit their layered structure by passing 'consensus' messages that guide inference towards good solutions. Experiments on a variety of problems show that the proposed technique leads to significantly more accurate inference results, not only when compared to standard EP and VMP, but also when compared to competitive bottom-up conditional models.


title: GMWB Riders in a Binomial Framework - Pricing, Hedging, and Diversification of Mortality Risk

abstract: We construct a binomial model for a guaranteed minimum withdrawal benefit (GMWB) rider to a variable annuity (VA) under optimal policyholder behaviour. The binomial model results in explicitly formulated perfect hedging strategies funded using only periodic fee income. We consider the separate perspectives of the insurer and policyholder and introduce a unifying relationship. Decompositions of the VA and GMWB contract into term-certain payments and options representing the guarantee and early surrender features are extended to the binomial framework. We incorporate an approximation algorithm for Asian options that significantly improves efficiency of the binomial model while retaining accuracy. Several numerical examples are provided which illustrate both the accuracy and the tractability of the binomial model. We extend the binomial model to include policy holder mortality and death benefits. Pricing, hedging, and the decompositions of the contract are extended to incorporate mortality risk. We prove limiting results for the hedging strategies and demonstrate mortality risk diversification. Numerical examples are provided which illustrate the effectiveness of hedging and the diversification of mortality risk under capacity constraints with finite pools.


title: Deep Structured learning for mass segmentation from Mammograms

abstract: In this paper, we present a novel method for the segmentation of breast masses from mammograms exploring structured and deep learning. Specifically, using structured support vector machine (SSVM), we formulate a model that combines different types of potential functions, including one that classifies image regions using deep learning. Our main goal with this work is to show the accuracy and efficiency improvements that these relatively new techniques can provide for the segmentation of breast masses from mammograms. We also propose an easily reproducible quantitative analysis to as- sess the performance of breast mass segmentation methodologies based on widely accepted accuracy and running time measurements on public datasets, which will facilitate further comparisons for this segmentation problem. In particular, we use two publicly available datasets (DDSM-BCRP and INbreast) and propose the computa- tion of the running time taken for the methodology to produce a mass segmentation given an input image and the use of the Dice index to quantitatively measure the segmentation accuracy. For both databases, we show that our proposed methodology produces competitive results in terms of accuracy and running time.


title: Parallel training of DNNs with Natural Gradient and Parameter Averaging

abstract: We describe the neural-network training framework used in the Kaldi speech recognition toolkit, which is geared towards training DNNs with large amounts of training data using multiple GPU-equipped or multi-core machines. In order to be as hardware-agnostic as possible, we needed a way to use multiple machines without generating excessive network traffic. Our method is to average the neural network parameters periodically (typically every minute or two), and redistribute the averaged parameters to the machines for further training. Each machine sees different data. By itself, this method does not work very well. However, we have another method, an approximate and efficient implementation of Natural Gradient for Stochastic Gradient Descent (NG-SGD), which seems to allow our periodic-averaging method to work well, as well as substantially improving the convergence of SGD on a single machine.


title: Legacy ExtraGalactic UV Survey (LEGUS) with The Hubble Space Telescope. I. Survey Description

abstract: The Legacy ExtraGalactic UV Survey (LEGUS) is a Cycle 21 Treasury program on the Hubble Space Telescope, aimed at the investigation of star formation and its relation with galactic environment in nearby galaxies, from the scales of individual stars to those of ~kpc-size clustered structures. Five-band imaging, from the near-ultraviolet to the I-band, with the Wide Field Camera 3, plus parallel optical imaging with the Advanced Camera for Surveys, is being collected for selected pointings of 50 galaxies within the local 12 Mpc. The filters used for the observations with the Wide Field Camera 3 are: F275W(2,704 A), F336W(3,355 A), F438W(4,325 A), F555W(5,308 A), and F814W(8,024 A); the parallel observations with the Advanced Camera for Surveys use the filters: F435W(4,328 A), F606W(5,921 A), and F814W(8,057 A). The multi-band images are yielding accurate recent (<~50 Myr) star formation histories from resolved massive stars and the extinction-corrected ages and masses of star clusters and associations. The extensive inventories of massive stars and clustered systems will be used to investigate the spatial and temporal evolution of star formation within galaxies. This will, in turn, inform theories of galaxy evolution and improve the understanding of the physical underpinning of the gas-star formation relation and the nature of star formation at high redshift. This paper describes the survey, its goals and observational strategy, and the initial science results. Because LEGUS will provide a reference survey and a foundation for future observations with JWST and with ALMA, a large number of data products are planned for delivery to the community.


title: Evolution and nucleosynthesis of asymptotic giant branch stellar models of low metallicity

abstract: We present stellar evolutionary tracks and nucleosynthetic predictions for a grid of stellar models of low- and intermediate-mass asymptotic giant branch (AGB) stars at $Z=0.001$ ([Fe/H]$=-1.2$). The models cover an initial mass range from 1 M$_{\odot}$ to 7 M$_{\odot}$. Final surface abundances and stellar yields are calculated for all elements from hydrogen to bismuth as well as isotopes up to the iron group. We present the first study of neutron-capture nucleosynthesis in intermediate-mass AGB models, including a super-AGB model, of [Fe/H] = $-1.2$. We examine in detail a low-mass AGB model of 2 M$_{\odot}$ where the $^{13}$C($\alpha$,$n$)$^{16}$O reaction is the main source of neutrons. We also examine an intermediate-mass AGB model of 5 M$_{\odot}$ where intershell temperatures are high enough to activate the $^{22}$Ne neutron source, which produces high neutron densities up to $\sim 10^{14}$ n cm$^{-3}$. Hot bottom burning is activated in models with $M \geq 3$ M$_{\odot}$. With the 3 M$_{\odot}$ model we investigate the effect of varying the extent in mass of the region where protons are mixed from the envelope into the intershell at the deepest extent of each third dredge-up. We compare the results of the low-mass models to three post-AGB stars with a metallicity of [Fe/H] $\simeq -1.2$. The composition is a good match to the predicted neutron-capture abundances except for Pb and we confirm that the observed Pb abundances are lower than what is calculated by AGB models.


title: Invariant four-variable automorphic kernel functions

abstract: Let $F$ be a number field, let $\mathbb{A}_F$ be its ring of adeles, and let $g_1,g_2,h_1,h_2 \in \mathrm{GL}_2(\mathbb{A}_F)$. Previously the author provided an absolutely convergent geometric expression for the four variable kernel function $$ \sum_{\pi} K_{\pi}(g_1,g_2)K_{\pi^{\vee}}(h_1,h_2)L(s,(\pi \times \pi^{\vee})^S), $$ where the sum is over isomorphism classes of cuspidal automorphic representations $\pi$ of $\mathrm{GL}_2(\mathbb{A}_F)$. Here $K_{\pi}$ is the typical kernel function representing the action of a test function on the space of the cuspidal automorphic representation $\pi$. In this paper we show how to use ideas from the circle method to provide an alternate expansion for the four variable kernel function that is visibly invariant under the natural action of $\mathrm{GL}_2(F) \times \mathrm{GL}_2(F)$.


title: On the Sum of a Prime and a Square-free Number

abstract: We prove that every integer greater than two may be written as the sum of a prime and a square-free number.


title: Throughput Optimization in Multi-Channel Cognitive Radios with Hard Deadline Constraints

abstract: In a cognitive radio scenario we consider a single secondary user (SU) accessing a multi-channel system. The SU senses the channels sequentially to detect if a primary user (PU) is occupying the channels, and stops its search to access a channel if it offers a significantly high throughput. The optimal stopping rule and power control problem is considered. The problem is formulated as a SU's throughput-maximization problem under a power, interference and packet delay constraints. We first show the effect of the optimal stopping rule on the packet delay, then solve this optimization problem for both the overlay system where the SU transmits only at the spectrum holes as well as the underlay system where tolerable interference (or tolerable collision probability) is allowed. We provide closed-form expressions for the optimal stopping rule, and show that the optimal power control strategy for this multi-channel problem is a modified water-filling approach. We extend the work to multiple SU scenario and show that when the number of SUs is large the complexity of the solution becomes smaller than that of the single SU case. We discuss the application of this problem in typical networks where packets arrive simultaneously and have the same departure deadline. We further propose an online adaptation policy to the optimal stopping rule that meets the packets' hard-deadline constraint and, at the same time, gives higher throughput than the offline policy.


title: Attosecond Timing in Optical-to-Electrical Conversion

abstract: The most frequency-stable sources of electromagnetic radiation are produced optically, and optical frequency combs provide the means for high fidelity frequency transfer across hundreds of terahertz and into the microwave domain. A critical step in this photonic-based synthesis of microwave signals is the optical-to-electrical conversion process. Here we show that attosecond (as) timing stability can be preserved across the opto-electronic interface of a photodiode, despite an intrinsic temporal response that is more than six orders of magnitude slower. The excess timing noise in the photodetection of a periodic train of ultrashort optical pulses behaves as flicker noise (1/f) with amplitude of 4 as/Sqrt(Hz) at 1 Hz offset. The corresponding fractional frequency fluctuations are 1.4x10-17 at 1 second and 5.5x10-20 at 1000 seconds. These results demonstrate that direct photodetection, as part of frequency-comb-based microwave synthesis, can support the timing performance of the best optical frequency standards, and thereby opens the possibility for generating microwave signals with significantly better stability than any existing source.


title: Derived categories of Grassmannians over integers and modular representation theory

abstract: In this paper we study the derived categories of coherent sheaves on Grassmannians $\operatorname{Gr}(k,n),$ defined over the ring of integers. We prove that the category $D^b(\operatorname{Gr}(k,n))$ has a semi-orthogonal decomposition, with components being full subcategories of the derived category of representations of $GL_k.$ This in particular implies existence of a full exceptional collection, which is a refinement of Kapranov's collection \cite{Kap}, which was constructed over a field of characteristic zero. We also describe the right dual semi-orthogonal decomposition which has a similar form, and its components are full subcategories of the derived category of representations of $GL_{n-k}.$ The resulting equivalences between the components of the two decompositions are given by a version of Koszul duality for strict polynomial functors. We also construct a tilting vector bundle on $\operatorname{Gr}(k,n).$ We show that its endomorphism algebra has two natural structures of a split quasi-hereditary algebra over $\mathbb{Z},$ and we identify the objects of $D^b(\operatorname{Gr}(k,n)),$ which correspond to the standard and costandard modules in both structures. All the results automatically extend to the case of arbitrary commutative base ring and the category of perfect complexes on the Grassmannian, by extension of scalars (base change). Similar results over fields of arbitrary characteristic were obtained independently in \cite{BLVdB}, by different methods.


title: Some remarks on stability of cones for the one-phase free boundary problem

abstract: We show that stable cones for the one-phase free boundary problem are hyperplanes in dimension $4$. As a corollary, both one and two-phase energy minimizing hypersurfaces are smooth in dimension $4$.


title: HIGHz: A Survey of the Most HI-Massive Galaxies at z~0.2

abstract: We present the results of the HIGHz Arecibo survey, which measured the HI content of 39 galaxies at redshift $z>0.16$ selected from the Sloan Digital Sky Survey. These are all actively star-forming, disk-dominated systems in relatively isolated environments, with stellar and HI masses larger than $10^{10}$ M$_\odot$ and redshifts $0.17\leq z\leq 0.25$. Our sample includes not only the highest-redshift detections of HI emission from individual galaxies to date, but also some of the most HI-massive systems known. Despite being exceptionally large, the HI reservoirs of these galaxies are consistent with what is expected from their ultraviolet and optical properties. This, and the fact that the galaxies lie on the baryonic Tully-Fisher relation, suggests that HIGHz systems are rare, scaled-up versions of local disk galaxies. We show that the most HI-massive galaxies discovered in the Arecibo Legacy Fast ALFA survey are the local analogues of HIGHz, and discuss the possible connection between our sample and the turbulent, gas-rich disks identified at $z\sim 1$. The HIGHz sample provides a first glimpse into the properties of the massive, HI-rich galaxies that will be detected at higher redshifts by the next-generation HI surveys with the Square Kilometer Array and its pathfinders.


title: Zero modes in vortex-fermion system with compact extra space

abstract: The existence of fermionic zero modes is shown in the presence of vortex configuration of pure $SU(2)$ gauge field on the manifold $M_4 \times S^2$. From the perspective of four-dimensional effective theory, these zero modes are almost the same as the Jackiw-Rossi type zero modes of the vortex-fermion system.


title: Concurrency Models with Causality and Events as Psi-calculi

abstract: Psi-calculi are a parametric framework for nominal calculi, where standard calculi are found as instances, like the pi-calculus, or the cryptographic spi-calculus and applied-pi. Psi-calculi have an interleaving operational semantics, with a strong foundation on the theory of nominal sets and process algebras. Much of the expressive power of psi-calculi comes from their logical part, i.e., assertions, conditions, and entailment, which are left quite open thus accommodating a wide range of logics. We are interested in how this expressiveness can deal with event-based models of concurrency. We thus take the popular prime event structures model and give an encoding into an instance of psi-calculi. We also take the recent and expressive model of Dynamic Condition Response Graphs (in which event structures are strictly included) and give an encoding into another corresponding instance of psi-calculi. The encodings that we achieve look rather natural and intuitive. Additional results about these encodings give us more confidence in their correctness.


title: Toward Sequentializing Overparallelized Protocol Code

abstract: In our ongoing work, we use constraint automata to compile protocol specifications expressed as Reo connectors into efficient executable code, e.g., in C. We have by now studied this automata based compilation approach rather well, and have devised effective solutions to some of its problems. Because our approach is based on constraint automata, the approach, its problems, and our solutions are in fact useful and relevant well beyond the specific case of compiling Reo. In this short paper, we identify and analyze two such rather unexpected problems.


title: Decay and persistence of spatial coherence during phonon-assisted relaxation in double quantum dots

abstract: We present a theoretical study of the evolution of spatial coherence during relaxation between delocalized exciton states in a pair of vertically stacked semiconductor quantum dots coupled to acoustic phonons. We show that spatial coherence can be transferred to the ground state even in a system of uncoupled non-identical quantum dots if a particular kind of degeneracy between the inter-level energy splittings is present. The phonon-assisted mechanism of coherence transfer leads to a dependence of the amount of the resulting coherence on the inter-dot distance and temperature. We analyze also the impact of carrier-phonon dynamics on a coupled system, where spatial coherence is present in the delocalized ground state.


title: On-the-fly Probabilistic Model Checking

abstract: Model checking approaches can be divided into two broad categories: global approaches that determine the set of all states in a model M that satisfy a temporal logic formula f, and local approaches in which, given a state s in M, the procedure determines whether s satisfies f. When s is a term of a process language, the model checking procedure can be executed "on-the-fly", driven by the syntactical structure of s. For certain classes of systems, e.g. those composed of many parallel components, the local approach is preferable because, depending on the specific property, it may be sufficient to generate and inspect only a relatively small part of the state space. We propose an efficient, on-the-fly, PCTL model checking procedure that is parametric with respect to the semantic interpretation of the language. The procedure comprises both bounded and unbounded until modalities. The correctness of the procedure is shown and its efficiency is compared with a global PCTL model checker on representative applications.


title: The Boolean Algebra of Cubical Areas as a Tensor Product in the Category of Semilattices with Zero

abstract: In this paper we describe a model of concurrency together with an algebraic structure reflecting the parallel composition. For the sake of simplicity we restrict to linear concurrent programs i.e. the ones with no loops nor branching. Such programs are given a semantics using cubical areas. Such a semantics is said to be geometric. The collection of all these cubical areas enjoys a structure of tensor product in the category of semi-lattice with zero. These results naturally extend to fully fledged concurrent programs up to some technical tricks.


title: From Orchestration to Choreography through Contract Automata

abstract: We study the relations between a contract automata and an interaction model. In the former model, distributed services are abstracted away as automata - oblivious of their partners - that coordinate with each other through an orchestrator. The interaction model relies on channel-based asynchronous communication and choreography to coordinate distributed services. We define a notion of strong agreement on the contract model, exhibit a natural mapping from the contract model to the interaction model, and give conditions to ensure that strong agreement corresponds to well-formed choreography.


title: A note on two notions of compliance

abstract: We establish a relation between two models of contracts: binary session types, and a model based on event structures and game-theoretic notions. In particular, we show that compliance in session types corresponds to the existence of certain winning strategies in game-based contracts.


title: Current driven asymmetric magnetization switching in perpendicularly magnetized CoFeB/MgO heterostructures

abstract: The flow of in-plane current through ultrathin magnetic heterostructures can cause magnetization switching or domain wall nucleation owing to bulk and interfacial effects. Within the magnetic layer, the current can create magnetic instabilities via spin transfer torques (STT). At interface(s), spin current generated from the spin Hall effect in a neighboring layer can exert torques, referred to as the spin Hall torques, on the magnetic moments. Here, we study current induced magnetization switching in perpendicularly magnetized CoFeB/MgO heterostructures with a heavy metal (HM) underlayer. Depending on the thickness of the HM underlayer, we find distinct differences in the inplane field dependence of the threshold switching current. The STT is likely responsible for the magnetization reversal for the thinner underlayer films whereas the spin Hall torques cause the switching for thicker underlayer films. For the latter, we find differences in the switching current for positive and negative currents and initial magnetization directions. We find that the growth process during the film deposition introduces an anisotropy that breaks the symmetry of the system and causes the asymmetric switching. The presence of such symmetry breaking anisotropy enables deterministic magnetization switching at zero external fields.


title: Loss-induced suppression and revival of lasing

abstract: Controlling and reversing the effects of loss are major challenges in optical systems. For lasers losses need to be overcome by a sufficient amount of gain to reach the lasing threshold. We show how to turn losses into gain by steering the parameters of a system to the vicinity of an exceptional point (EP), which occurs when the eigenvalues and the corresponding eigenstates of a system coalesce. In our system of coupled microresonators, EPs are manifested as the loss-induced suppression and revival of lasing. Below a critical value, adding loss annihilates an existing Raman laser. Beyond this critical threshold, lasing recovers despite the increasing loss, in stark contrast to what would be expected from conventional laser theory. Our results exemplify the counterintuitive features of EPs and present an innovative method for reversing the effect of loss.