title: On semilinear representations of the infinite symmetric group

abstract: In this note the smooth (i.e. with open stabilizers) linear and {\sl semilinear} representations of certain permutation groups (such as infinite symmetric group or automorphism group of an infinite-dimensional vector space over a finite field) are studied. Many results here are well-known to the experts, at least in the case of {\sl linear representations} of symmetric group. The presented results suggest, in particular, that an analogue of Hilbert's Theorem 90 should hold: in the case of faithful action of the group on the base field the irreducible smooth semilinear representations are one-dimensional (and trivial in appropriate sense).


title: Towards a Lagrange-Newton approach for PDE constrained shape optimization

abstract: The novel Riemannian view on shape optimization developed in [Schulz, FoCM, 2014] is extended to a Lagrange-Newton approach for PDE constrained shape optimization problems. The extension is based on optimization on Riemannian vector space bundles and exemplified for a simple numerical example.


title: Exact Recovery in the Stochastic Block Model

abstract: The stochastic block model (SBM) with two communities, or equivalently the planted bisection model, is a popular model of random graph exhibiting a cluster behaviour. In the symmetric case, the graph has two equally sized clusters and vertices connect with probability $p$ within clusters and $q$ across clusters. In the past two decades, a large body of literature in statistics and computer science has focused on providing lower-bounds on the scaling of $|p-q|$ to ensure exact recovery. In this paper, we identify a sharp threshold phenomenon for exact recovery: if $\alpha=pn/\log(n)$ and $\beta=qn/\log(n)$ are constant (with $\alpha>\beta$), recovering the communities with high probability is possible if $\frac{\alpha+\beta}{2} - \sqrt{\alpha \beta}>1$ and impossible if $\frac{\alpha+\beta}{2} - \sqrt{\alpha \beta}<1$. In particular, this improves the existing bounds. This also sets a new line of sight for efficient clustering algorithms. While maximum likelihood (ML) achieves the optimal threshold (by definition), it is in the worst-case NP-hard. This paper proposes an efficient algorithm based on a semidefinite programming relaxation of ML, which is proved to succeed in recovering the communities close to the threshold, while numerical experiments suggest it may achieve the threshold. An efficient algorithm which succeeds all the way down to the threshold is also obtained using a partial recovery algorithm combined with a local improvement procedure.


title: Multilayer formation and evaporation of deuterated ices in prestellar and protostellar cores

abstract: Extremely large deuteration of several molecules has been observed towards prestellar cores and low-mass protostars for a decade. New observations performed towards low-mass protostars suggest that water presents a lower deuteration in the warm inner gas than in the cold external envelope. We coupled a gas-grain astrochemical model with a one-dimension model of collapsing core to properly follow the formation and the deuteration of interstellar ices as well as their subsequent evaporation in the low-mass protostellar envelopes with the aim of interpreting the spatial and temporal evolutions of their deuteration. The astrochemical model follows the formation and the evaporation of ices with a multilayer approach and also includes a state-of-the-art deuterated chemical network by taking the spin states of H$_2$ and light ions into account. Because of their slow formation, interstellar ices are chemically heterogeneous and show an increase of their deuterium fractionation towards the surface. The differentiation of the deuteration in ices induces an evolution of the deuteration within protostellar envelopes. The warm inner region is poorly deuterated because it includes the whole molecular content of ices while the deuteration predicted in the cold external envelope scales with the highly deuterated surface of ices. We are able to reproduce the observed evolution of water deuteration within protostellar envelopes but we are still unable to predict the super-high deuteration observed for formaldehyde and methanol. Finally, the extension of this study to the deuteration of complex organics (COMs), important for the prebiotic chemistry, shows a good agreement with the observations, suggesting that we can use the deuteration to retrace their mechanisms and their moments of formation.


title: Limiting mechanism for critical current in topologically frustrated Josephson junctions

abstract: Eutectic Sr$_2$RuO$_4$-Ru samples with $\mu$m-sized Ru-metal inclusions support inhomogeneous superconductivity above the bulk transition of Sr$_2$RuO$_4$ in the so-called 3-Kelvin phase. In Pb/Ru/Sr$_2$RuO$_4$ Josephson junctions as realized by Maeno et al., a Pb film is indirectly coupled to the superconductor Sr$_2$RuO$_4$ mediated by the proximity-induced superconducting Ru-inclusions, yielding an extended Josephson contact through the interface between Ru and Sr$_2$RuO$_4$. Motivated by this experimental setup, we formulate a sine-Gordon model for the Josephson phase of the interface, assuming a simple cylindrical shape for the Ru-inclusion hosting the proximity-induced $s$-wave superconducting phase. Considering the Sr$_2$RuO$_4$ as a chiral $p$-wave superconductor, we discuss two types of Josephson junctions, a frustrated one due to the nature of the order parameter in Sr$_2$RuO$_4$, and an unfrustrated one for the topologically trivial 3-Kelvin phase. While the latter situation displays standard junction behavior, the former yields an unusual limiting mechanism for the critical current, based on a pinning-depinning transition of a spontaneously induced magnetic flux driven by an externally applied current. We analyze different coupling limits and show that different critical currents can arise for the two topologies. This concept fits well to recent experimental data obtained for the above setup showing an anomalous temperature dependence of the critical current at the transition temperature $T_c $ of bulk Sr$_2$RuO$_4$.


title: Spatial statistics and attentional dynamics in scene viewing

abstract: In humans and in foveated animals visual acuity is highly concentrated at the center of gaze, so that choosing where to look next is an important example of online, rapid decision making. Computational neuroscientists have developed biologically-inspired models of visual attention, termed saliency maps, which successfully predict where people fixate on average. Using point process theory for spatial statistics, we show that scanpaths contain, however, important statistical structure, such as spatial clustering on top of distributions of gaze positions. Here we develop a dynamical model of saccadic selection that accurately predicts the distribution of gaze positions as well as spatial clustering along individual scanpaths. Our model relies on, first, activation dynamics via spatially- limited (foveated) access to saliency information, and, second, a leaky memory process controlling the re-inspection of target regions. This theoretical framework models a form of context-dependent decision-making, linking neural dynamics of attention to behavioral gaze data.


title: Matchings in Benjamini-Schramm convergent graph sequences

abstract: We introduce the matching measure of a finite graph as the uniform distribution on the roots of the matching polynomial of the graph. We analyze the asymptotic behavior of the matching measure for graph sequences with bounded degree. A graph parameter is said to be estimable if it converges along every Benjamini-Schramm convergent sparse graph sequence. We prove that the normalized logarithm of the number of matchings is estimable. We also show that the analogous statement for perfect matchings already fails for d-regular bipartite graphs for any fixed d at least 3. The latter result relies on analyzing the probability that a randomly chosen perfect matching contains a particular edge. However, for any sequence of d-regular bipartite graphs converging to the d-regular tree, we prove that the normalized logarithm of the number of perfect matchings converges. This applies to random d-regular bipartite graphs. We show that the limit equals to the exponent in Schrijver's lower bound on the number of perfect matchings. Our analytic approach also yields a short proof for the Nguyen-Onak (also Elek--Lippner) theorem saying that the matching ratio is estimable. In fact, we prove the slightly stronger result that the independence ratio is estimable for claw-free graphs.


title: Fast and Fuzzy Private Set Intersection

abstract: Private Set Intersection (PSI) is usually implemented as a sequence of encryption rounds between pairs of users, whereas the present work implements PSI in a simpler fashion: each set only needs to be encrypted once, after which each pair of users need only one ordinary set comparison. This is typically orders of magnitude faster than ordinary PSI at the cost of some ``fuzziness" in the matching, which may nonetheless be tolerable or even desirable. This is demonstrated in the case where the sets consist of English words processed with WordNet.


title: The description of strong correlation within self-consistent Green's function second-order perturbation theory

abstract: We report an implementation of self-consistent Green's function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Green's function and self-energy are built on the imaginary frequency and imaginary time domain respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H$_{32}$ finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.


title: Correlations in the Monte Carlo Glauber model

abstract: Event-by-event fluctuations of observables are often modeled using the Monte Carlo Glauber model, in which the energy is initially deposited in sources associated with wounded nucleons. In this paper, we analyze in detail the correlations between these sources in proton-nucleus and nucleus-nucleus collisions. There are correlations arising from nucleon-nucleon correlations within each nucleus, and correlations due to the collision mechanism, which we dub twin correlations. We investigate this new phenomenon in detail. At the RHIC and LHC energies, correlations are found to have modest effects on size and eccentricity fluctuations, such that the Glauber model produces to a good approximation a collection of independent sources.


title: Invariants of Legendrian and transverse knots in monopole knot homology

abstract: We use the contact invariant defined in [2] to construct a new invariant of Legendrian knots in Kronheimer and Mrowka's monopole knot homology theory (KHM), following a prescription of Stipsicz and V\'ertesi. Our Legendrian invariant improves upon an earlier Legendrian invariant in KHM defined by the second author in several important respects. Most notably, ours is preserved by negative stabilization. This fact enables us to define a transverse knot invariant in KHM via Legendrian approximation. It also makes our invariant a more likely candidate for the monopole Floer analogue of the "LOSS" invariant in knot Floer homology. Like its predecessor, our Legendrian invariant behaves functorially with respect to Lagrangian concordance. We show how this fact can be used to compute our invariant in several examples. As a byproduct of our investigations, we provide the first infinite family of nonreversible Lagrangian concordances between prime knots.


title: Study of the C IV ($\lambda\lambda$1548.187 - 1550.772), Si IV ($\lambda\lambda$1393.755 - 1402.770) and O IV ($\lambda$1401.156) Regions of the QSO J021327.25-001446.93

abstract: Broad Absorption Line Regions - BALR are composed of a number of successive independent absorbing density layers. Using the GR model, we analyze the UV Si IV ({\lambda}{\lambda}1393.755 - 1402.770), O IV ({\lambda}1401.156) and C IV ({\lambda}{\lambda}1548.187 - 1550.772) resonance lines in the spectra of a certain QSO and discuss the results concerning its kinematic properties (rotational, radial and random velocities).


title: Free fall misconceptions: results of a graph based pre test of sophomore civil engineering students

abstract: A partially unusual behaviour was found among 14 sophomore students of civil engineering who took a pre test for a free fall laboratory session, in the context of a general mechanics course. An analysis contemplating mathematics models and physics models consistency was made. In all cases, the students presented evidence favoring a correct free fall acceleration model, whilst their position component versus time, and velocity component versus time graphs revealed complex misconceptions both on the physical phenomenon and it's implicit mathematics consistency. The last suggests an inability to make satisfactory connections through definitions between graphed variables. In other words, evidence strongly suggests that students are perfectly able to memorize the free fall acceleration model, whilst not understanding it's significance at any level. This small study originated the develope and validation of a tutorial on free fall graphs for position, velocity and acceleration models, as part of a following cross universities major project.


title: Instanton Floer homology and contact structures

abstract: We define an invariant of contact 3-manifolds with convex boundary using Kronheimer and Mrowka's sutured instanton Floer homology theory. To the best of our knowledge, this is the first invariant of contact manifolds -- with or without boundary -- defined in the instanton Floer setting. We prove that our invariant vanishes for overtwisted contact structures and is nonzero for contact manifolds with boundary which embed into Stein fillable contact manifolds. Moreover, we propose a strategy by which our contact invariant might be used to relate the fundamental group of a closed contact 3-manifold to properties of its Stein fillings. Our construction is inspired by a reformulation of a similar invariant in the monopole Floer setting defined by the authors in [1].


title: Scenario for delocalization in translation invariant systems

abstract: We investigate the possibility of Many-Body Localization in translation invariant Hamiltonian systems, which was recently brought up by several authors. A key feature of Many-Body Localized disordered systems is recovered, namely the fact that resonant spots are rare and far-between. However, we point out that resonant spots are mobile, unlike in models with strong quenched disorder, and that these mobile spots constitute a possible mechanism for delocalization, albeit possibly only on very long timescales. In some models, this argument for delocalization can be made very explicit in first order of perturbation theory in the hopping. For models where this does not work, we present instead a non-perturbative argument that relies solely on ergodicity inside the resonant spots.


title: The Logic of Identity: Distinguishability and Indistinguishability in Classical and Quantum Physics

abstract: The suggestion that particles of the same kind may be indistinguishable in a fundamental sense, even so that challenges to traditional notions of individuality and identity may arise, has first come up in the context of classical statistical mechanics. In particular, the Gibbs paradox has sometimes been interpreted as a sign of the untenability of the classical concept of a particle and as a premonition that quantum theory is needed. This idea of a quantum connection stubbornly persists in the literature, even though it has also been criticized frequently. Here we shall argue that although this criticism is justified, the proposed alternative solutions have often been wrong and have not put the paradox in its right perspective. In fact, the Gibbs paradox is unrelated to fundamental issues of particle identity; only distinguishability in a pragmatic sense plays a role (in this we develop ideas of van Kampen [10]), and in principle the paradox always is there as long as the concept of a particle applies at all. In line with this we show that the paradox survives even in quantum mechanics, in spite of the quantum mechanical (anti-)symmetrization postulates.


title: The Star Formation Histories of Local Group Dwarf Galaxies II. Searching For Signatures of Reionization

abstract: We search for signatures of reionization in the star formation histories (SFHs) of 38 Local Group dwarf galaxies (10$^4$ $<$ M$_{\star}$ $<$ 10$^9$ M$_{\odot}$). The SFHs are derived from color-magnitude diagrams using archival Hubble Space Telescope/Wide Field Planetary Camera 2 imaging. Only five quenched galaxies (And V, And VI, And XIII, Leo IV, Hercules) are consistent with forming the bulk of their stars before reionization, when full uncertainties are considered. Observations of 13 of the predicted `true fossils' identified by Bovill & Ricotti show that only two (Hercules and Leo IV) indicate star formation quenched by reionization. However, both are within the virial radius of the Milky Way and evidence of tidal disturbance complicates this interpretation. We argue that the late-time gas capture scenario posited by Ricotti for the low mass, gas-rich, and star-forming fossil candidate Leo T is observationally indistinguishable from simple gas retention. Given the ambiguity between environmental effects and reionization, the best reionization fossil candidates are quenched low mass field galaxies (e.g., KKR 25).


title: How to Ask for a Favor: A Case Study on the Success of Altruistic Requests

abstract: Requests are at the core of many social media systems such as question & answer sites and online philanthropy communities. While the success of such requests is critical to the success of the community, the factors that lead community members to satisfy a request are largely unknown. Success of a request depends on factors like who is asking, how they are asking, when are they asking, and most critically what is being requested, ranging from small favors to substantial monetary donations. We present a case study of altruistic requests in an online community where all requests ask for the very same contribution and do not offer anything tangible in return, allowing us to disentangle what is requested from textual and social factors. Drawing from social psychology literature, we extract high-level social features from text that operationalize social relations between recipient and donor and demonstrate that these extracted relations are predictive of success. More specifically, we find that clearly communicating need through the narrative is essential and that that linguistic indications of gratitude, evidentiality, and generalized reciprocity, as well as high status of the asker further increase the likelihood of success. Building on this understanding, we develop a model that can predict the success of unseen requests, significantly improving over several baselines. We link these findings to research in psychology on helping behavior, providing a basis for further analysis of success in social media systems.


title: Accreting planets as dust dams in `transition' discs

abstract: We investigate under what circumstances an embedded planet in a protoplanetary disc may sculpt the dust distribution such that it observationally presents as a `transition' disc. We concern ourselves with `transition' discs that have large holes ($\gtrsim 10$ AU) and high accretion rates ($\sim 10^{-9}-10^{-8}$ M$_\odot$ yr$^{-1}$). Particularly, those discs which photoevaporative models struggle to explain. Assuming the standard picture for how massive planets sculpt their parent discs, along with the observed accretion rates in `transition' discs, we find that the accretion luminosity from the forming planet is significant, and can dominate over the stellar luminosity at the gap edge. This planetary accretion luminosity can apply a significant radiation pressure to small ($s\lesssim 1\mu$m) dust particles provided they are suitably decoupled from the gas. Secular evolution calculations that account for the evolution of the gas and dust components in a disc with an embedded, accreting planet, show that only with the addition of the radiation pressure can we explain the full observed characteristics of a `transition' disc (NIR dip in the SED, mm cavity and high accretion rate). At suitably high planet masses ($\gtrsim 3-4$ M$_J$), radiation pressure from the accreting planet is able to hold back the small dust particles, producing a heavily dust-depleted inner disc that is optically thin (vertically and radially) to Infra-Red radiation. We use our models to calculate synthetic observations and present a observational evolutionary scenario for a forming planet, sculpting its parent disc. The planet-disc system will present as a `transition' disc with a dip in the SED, only when the planet mass and planetary accretion rate is high enough. At other times it will present as a disc with a primordial SED, but with a cavity in the mm, as observed in a handful of protoplanetary discs.


title: PynPoint Code for Exoplanet Imaging

abstract: We announce the public release of PynPoint, a Python package that we have developed for analysing exoplanet data taken with the angular differential imaging observing technique. In particular, PynPoint is designed to model the point spread function of the central star and to subtract its flux contribution to reveal nearby faint companion planets. The current version of the package does this correction by using a principal component analysis method to build a basis set for modelling the point spread function of the observations. We demonstrate the performance of the package by reanalysing publicly available data on the exoplanet beta Pictoris b, which consists of close to 24,000 individual image frames. We show that PynPoint is able to analyse this typical data in roughly 1.5 minutes on a Mac Pro, when the number of images is reduced by co-adding in sets of 5. The main computational work parallelises well as a result of a reliance on SciPy and NumPy functions. For this calculation the peak memory load is 6Gb, which can be run comfortably on most workstations. A simpler calculation, by co-adding over 50, takes 3 seconds with a peak memory usage of 600 Mb. This can be performed easily on a laptop. In developing the package we have modularised the code so that we will be able to extend functionality in future releases, through the inclusion of more modules, without it affecting the users application programming interface. We distribute the PynPoint package through the central PyPi sever, and the documentation is available online (http://pynpoint.ethz.ch).


title: The Effect of Orbital Eccentricity on the Dynamical Evolution of Star Clusters

abstract: We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker mean tidal field. However, we find that perigalactic passes and tidal heating due to an eccentric orbit can partially compensate for the decreased mean tidal field by energizing stars to higher velocities and stripping additional stars from the cluster, accelerating the relaxation process. We find that the corresponding circular orbit which best describes the evolution of a cluster on an eccentric orbit is much less than its semi-major axis or time averaged galactocentric distance. Since clusters spend the majority of their lifetimes near apogalacticon, the properties of clusters which appear very dynamically evolved for a given galactocentric distance can be explained by an eccentric orbit. Additionally we find that the evolution of the slope of the mass function within the core radius is roughly orbit-independent, so it could place additional constraints on the initial mass and initial size of globular clusters with solved orbits. We use our results to demonstrate how the orbit of Milky Way globular clusters can be constrained given standard observable parameters like galactocentric distance and the slope of the mass function. We then place constraints on the unsolved orbits of NGC 1261,NGC 6352, NGC 6496, and NGC 6304 based on their positions and mass functions.


title: The Onset of Massive Star Formation: The Evolution of Temperature and Density Structure in an Infrared Dark Cloud

abstract: We present new NH3 (1,1), (2,2), and (4,4) observations from the Karl G. Jansky Very Large Array (VLA) compiled with work in the literature to explore the range of conditions observed in young, massive star-forming regions. To sample the effects of evolution independent from those of distance/resolution, abundance, and large-scale environment, we compare clumps in different evolutionary stages within a single Infrared Dark Cloud (IRDC), G32.02+0.06. We find that the early stages of clustered star formation are characterized by dense, parsec-scale filamentary structures interspersed with complexes of dense cores (<0.1 pc cores clustered in complexes separated by ~1 pc) with masses from about 10 to 100 Msun. The most quiescent core is the most extended while the star-forming cores are denser and more compact, showing very similar column density structure before and shortly after the onset of massive star formation, with peak surface densities >~ 1 g cm^-2. Quiescent cores and filaments show smoothly varying temperatures from 10-20 K, rising to over 40 K in star-forming cores. We calculate the virial parameters for 16 cores and find that the level of support provided by turbulence is generally insufficient to support them against gravitational collapse (alpha_vir ~ 0.6). The star-forming filaments show smooth velocity fields, punctuated by discontinuities at the sites of active star formation. We discuss the Massive Molecular Filament (MMF; M > 10^5 Msun, l > 60 pc) hosting the IRDC, hypothesizing that it may have been shaped by previous generations of massive stars.


title: Testing the universality of the IMF with Bayesian statistics: young clusters

abstract: The universality of the stellar initial mass function (IMF) is tested using Bayesian statistics with a sample of eight young Galactic stellar clusters (IC 348, ONC, NGC 2024, NGC 6611, NGC 2264, $\rho$ Ophiuchi, Chameleon I, and Taurus). We infer the posterior probability distribution function (pPDF) of the IMF parameters when the likelihood function is described by a tapered power law function, a lognormal distribution at low masses coupled to a power law at higher masses, and a multi-component power law function. The inter-cluster comparison of the pPDFs of the IMF parameters for each likelihood function shows that these distributions do not overlap within the $1\sigma$ uncertainty level. Furthermore, the most probable values of the IMF parameters for most of the clusters deviate substantially from their values for the Galactic field stellar IMF. We also quantify the effects of taking into account the completeness correction as well as the uncertainties on the measured masses. The inclusion of the former affects the inferred pPDFs of the slope of the IMF at the low mass end while considering the latter affects the pPDFs of the slope of the IMF in the intermediate- to high mass regime. As variations are observed in all of the IMF parameters at once and for each of the considered likelihood functions, even for completeness corrected samples, we argue that the observed variations are real and significant, at least for the sample of eight clusters considered in this work. The results presented here clearly show that the IMF is not universal.


title: TASTE IV. Refining ephemeris and orbital parameters for HAT-P-20b and WASP-1b

abstract: We present four new light curves of transiting exoplanets WASP-1b and HAT-P-20b, observed within the TASTE (The Asiago Search for Transit timing variations of Exoplanets) project. We re-analyzed light curves from the literature in a homogeneous way, calculating a refined ephemeris and orbital-physical parameters for both objects. WASP-1b does not show any significant Transit Timing Variation signal at the 120 s-level. As for HAT-P-20b, we detected a deviation from our re-estimated linear ephemeris that could be ascribed to the presence of a perturber or, more probably, to a previously unnoticed high level of stellar activity. The rotational period of HAT-P-20 A we obtained from archival data (P_rot ~ 14.5 days), combined with its optical variability and strong emission of CaII H&K lines, is consistent with a young stellar age (< 1 Gyr) and support the hypothesis that stellar activity may be responsible of the measured deviations of the transit times.


title: Implementation of the Dicke lattice model in hybrid quantum system arrays

abstract: Generalized Dicke models can be implemented in hybrid quantum systems built from ensembles of nitrogen-vacancy (NV) centers in diamond coupled to superconducting microwave cavities. By engineering cavity assisted Raman transitions between two spin states of the NV defect, a fully tunable model for collective light-matter interactions in the ultra-strong coupling limit can be obtained. Our analysis of the resulting non-equilibrium phases for a single cavity and for coupled cavity arrays shows that different superradiant phase transitions can be observed using existing experimental technologies, even in the presence of large inhomogeneous broadening of the spin ensemble. The phase diagram of the Dicke lattice model displays distinct features induced by dissipation, which can serve as a genuine experimental signature for phase transitions in driven open quantum systems.


title: Permuted Random Walk Exits Typically in Linear Time

abstract: Given a permutation sigma of the integers {-n,-n+1,...,n} we consider the Markov chain X_{sigma}, which jumps from k to sigma (k\pm 1) equally likely if k\neq -n,n. We prove that the expected hitting time of {-n,n} starting from any point is Theta(n) with high probability when sigma is a uniformly chosen permutation. We prove this by showing that with high probability, the digraph of allowed transitions is an Eulerian expander; we then utilize general estimates of hitting times in directed Eulerian expanders.


title: Gravito-Turbulent Disks in 3D: Turbulent Velocities vs. Depth

abstract: Characterizing turbulence in protoplanetary disks is crucial for understanding how they accrete and spawn planets. Recent measurements of spectral line broadening promise to diagnose turbulence, with different lines probing different depths. We use 3D local hydrodynamic simulations of cooling, self-gravitating disks to resolve how motions driven by "gravito-turbulence" vary with height. We find that gravito-turbulence is practically as vigorous at altitude as at depth: even though gas at altitude is much too rarefied to be itself self-gravitating, it is strongly forced by self-gravitating overdensities at the midplane. The long-range nature of gravity means that turbulent velocities are nearly uniform vertically, increasing by just a factor of 2 from midplane to surface, even as the density ranges over nearly three orders of magnitude. The insensitivity of gravito-turbulence to height contrasts with the behavior of disks afflicted by the magnetorotational instability (MRI); in the latter case, non-circular velocities increase by at least a factor of 15 from midplane to surface, with various non-ideal effects only magnifying this factor. The distinct vertical profiles of gravito-turbulence vs. MRI turbulence may be used in conjunction with measurements of non-thermal linewidths at various depths to identify the source of transport in protoplanetary disks.


title: Learning with many experts: model selection and sparsity

abstract: Experts classifying data are often imprecise. Recently, several models have been proposed to train classifiers using the noisy labels generated by these experts. How to choose between these models? In such situations, the true labels are unavailable. Thus, one cannot perform model selection using the standard versions of methods such as empirical risk minimization and cross validation. In order to allow model selection, we present a surrogate loss and provide theoretical guarantees that assure its consistency. Next, we discuss how this loss can be used to tune a penalization which introduces sparsity in the parameters of a traditional class of models. Sparsity provides more parsimonious models and can avoid overfitting. Nevertheless, it has seldom been discussed in the context of noisy labels due to the difficulty in model selection and, therefore, in choosing tuning parameters. We apply these techniques to several sets of simulated and real data.


title: Elastic free energy drives the shape of prevascular solid tumors

abstract: It is well established that the mechanical environment influences cell functions in health and disease. Here, we address how the mechanical environment influences tumor growth, in particular, the shape of solid tumors. In an in vitro tumor model, which isolates mechanical interactions between tumor cells and a hydrogel, we find that tumors grow as ellipsoids, resembling the same, oft-reported observation of in vivo tumors. Specifically, an oblate ellipsoidal tumor shape robustly occurs when the tumors grow in hydrogels that are stiffer than the tumors, but when they grow in more compliant hydrogels they remain closer to spherical in shape. Using large scale, nonlinear elasticity computations we show that the oblate ellipsoidal shape minimizes the elastic free energy of the tumor-hydrogel system. Having eliminated a number of other candidate explanations, we hypothesize that minimization of the elastic free energy is the reason for predominance of the experimentally observed ellipsoidal shape. This result may hold significance for explaining the shape progression of early solid tumors in vivo and is an important step in understanding the processes underlying solid tumor growth.


title: Neutrino diffusive transport in hot quark matter: a detailed analysis

abstract: We perform an extensive analysis of neutrino diffusion in quark matter within the MIT bag model at arbitrary temperature and degeneracy. We examine in detail the contribution of each of the relevant weak interaction processes to the total neutrino opacity and evaluate the effect of the strange quark mass, the bag constant, and the QCD perturbative corrections to the MIT bag model. We also investigate the anisotropic contribution to the neutrino distribution function in scatterings, the mean energy transfer and the mean scattering angle. The density and temperature dependence of the diffusion coefficients $D_n$ that govern the cooling and deleptonization of a compact star is shown in detail. Finally, our numerical results for the neutrino mean free paths are compared against known analytic approximations. We conclude that neutrino scattering constitutes a significant portion of the total neutrino opacity in leptonized quark matter and neutrino-quark scattering is, in general, very similar to neutrino-electron scattering with respect to both mean energy transfer per scattering and mean scattering angle.


title: Effects of Sampling Methods on Prediction Quality. The Case of Classifying Land Cover Using Decision Trees

abstract: Clever sampling methods can be used to improve the handling of big data and increase its usefulness. The subject of this study is remote sensing, specifically airborne laser scanning point clouds representing different classes of ground cover. The aim is to derive a supervised learning model for the classification using CARTs. In order to measure the effect of different sampling methods on the classification accuracy, various experiments with varying types of sampling methods, sample sizes, and accuracy metrics have been designed. Numerical results for a subset of a large surveying project covering the lower Rhine area in Germany are shown. General conclusions regarding sampling design are drawn and presented.


title: Algorithm Instance Games

abstract: This paper introduces algorithm instance games (AIGs) as a conceptual classification applying to games in which outcomes are resolved from joint strategies algorithmically. For such games, a fundamental question asks: How do the details of the algorithm's description influence agents' strategic behavior? We analyze two versions of an AIG based on the set-cover optimization problem. In these games, joint strategies correspond to instances of the set-cover problem, with each subset (of a given universe of elements) representing the strategy of a single agent. Outcomes are covers computed from the joint strategies by a set-cover algorithm. In one variant of this game, outcomes are computed by a deterministic greedy algorithm, and the other variant utilizes a non-deterministic form of the greedy algorithm. We characterize Nash equilibrium strategies for both versions of the game, finding that agents' strategies can vary considerably between the two settings. In particular, we find that the version of the game based on the deterministic algorithm only admits Nash equilibrium in which agents choose strategies (i.e., subsets) containing at most one element, with no two agents picking the same element. On the other hand, in the version of the game based on the non-deterministic algorithm, Nash equilibrium strategies can include agents with zero, one, or every element, and the same element can appear in the strategies of multiple agents.


title: Emergent Spacetime in Stochastically Evolving Dimensions

abstract: Changing the dimensionality of the space-time at the smallest and largest distances has manifold theoretical advantages. If the space is lower dimensional in the high energy regime, then there are no ultraviolet divergencies in field theories, it is possible to quantize gravity, and the theory of matter plus gravity is free of divergencies or renormalizable. If the space is higher dimensional at cosmological scales, then some cosmological problems (including the cosmological constant problem) can be attacked from a completely new perspective. In this paper, we construct an explicit model of "evolving dimensions" in which the dimensions open up as the temperature of the universe drops. We adopt the string theory framework in which the dimensions are fields that live on the string worldsheet, and add temperature dependent mass terms for them. At the Big Bang, all the dimensions are very heavy and are not excited. As the universe cools down, dimensions open up one by one. Thus, the dimensionality of the space we live in depends on the energy or temperature that we are probing. In particular, we provide a kinematic Brandenberger-Vafa argument for how a discrete {\it causal set}, and eventually a continuum $(3+1)$-dim spacetime along with Einstein gravity emerge in the Infrared from the worldsheet action. The $(3+1)$-dim Planck mass and the string scale become directly related, {\it without any} compactification. Amongst other predictions, we argue that LHC might be blind to new physics even if it comes at the TeV scale. In contrast, cosmic ray experiments, especially those that can register the very beginning of the shower, and collisions with high multiplicity and density of particles, might be sensitive to the dimensional cross-over.


title: Fatou and relative Fatou theorem for subordinate Brownian motions with Gaussian components on smooth domains

abstract: We prove relative Fatou's theorem for nonnegative harmonic functions with respect to a large class of killed subordinate Brownian motions with Gaussian components in bounded $C^{1,1}$ open sets in $\mathbb{R}^{d}$, $d\geq 2$, which asserts the existence of nontangential limit of the ratio of two harmonic functions with respect to the killed processes. When $D=B(x_{0},r)$ is a ball we prove Fatou theorem. That is, we establish the existence of nontangential limit of a single nonnegative harmonic function. We also prove this is the best result possible by showing that there is a nonnegative harmonic function which does not have a tangential limit a.e. when $d=2$ and $D=B(0,1)$.


title: A $U_A(1)$ symmetry restoration scenario supported by the generalized Witten-Veneziano relation and its analytic solution

abstract: The Witten-Veneziano relation, or, alternatively, its generalization proposed by Shore, facilitates understanding and describing the complex of eta and eta' mesons. We present an analytic, closed-form solution to Shore's equations which gives results on the eta-eta' complex in full agreement with results previously obtained numerically. Although the Witten-Veneziano relation and Shore's equations are related, the ways they were previously used in the context of dynamical models to calculate eta and eta' properties, were rather different. However, with the analytic solution, the calculation can be formulated similarly to the approach through the Witten-Veneziano relation, and with some conceptual improvements. In the process, one strengthens the arguments in favor of a possible relation between the U_A(1) and SU_A(3) chiral symmetry breaking and restoration. To test this scenario, the experiments such as those at RHIC, NICA and FAIR, which extend the RHIC (and LHC) high-temperature scans also to the finite-density parts of the QCD phase diagram, should pay particular attention to the signatures from the eta'-eta complex indicating the symmetry restoration.


title: Quantum impurity in a Tomonaga-Luttinger liquid: continuous-time quantum Monte Carlo approach

abstract: We develop a continuous-time quantum Monte Carlo (CTQMC) method for quantum impurities coupled to interacting quantum wires described by a Tomonaga-Luttinger liquid. The method is negative-sign free for any values of the Tomonaga-Luttinger parameter, which is rigorously proved, and thus, efficient low-temperature calculations are possible. Duality between electrons and bosons in one dimensional systems allows us to construct a simple formula for the CTQMC algorithm in these systems. We show that the CTQMC for Tomonaga-Luttinger liquids can be implemented with only minor modifications of previous CTQMC codes developed for impurities coupled to non-interacting fermions. We apply this method to the Kane-Fisher model of a potential scatterer in a spin-less quantum wire and to a single spin coupled with the edge state of a two-dimensional topological insulator assuming an anisotropic XXZ coupling. Various dynamical response functions such as the electron Green's function and spin-spin correlation functions are calculated numerically and their scaling properties are discussed.


title: Heavy-quark production and elliptic flow in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV

abstract: We present measurements of electrons and positrons from the semileptonic decays of heavy-flavor hadrons at midrapidity ($|y|<$ 0.35) in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV. The data were collected in 2010 by the PHENIX experiment that included the new hadron-blind detector. The invariant yield of electrons from heavy-flavor decays is measured as a function of transverse momentum in the range $1<p_T^e<5$ GeV/$c$. The invariant yield per binary collision is slightly enhanced above the $p$$+$$p$ reference in Au$+$Au 0%--20%, 20%--40% and 40%--60% centralities at a comparable level. This may be a result of the interplay between initial-state Cronin effects, final-state flow, and energy loss for heavy-quark production at this low beam energy. The $v_2$ of electrons from heavy-flavor decays is nonzero when averaged between $1.3<p_T^e<2.5$ GeV/$c$ from $0<{\rm centrality}<40$% collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV. For 20%--40% centrality collisions, the $v_2$ at $\sqrt{s_{_{NN}}}=62.4$ GeV is smaller than that for heavy flavor decays at $\sqrt{s_{_{NN}}}=200$ GeV. The $v_2$ of the electrons from heavy-flavor decay at the lower beam energy is also smaller than $v_2$ for pions. Both results indicate that the heavy-quarks interact with the medium formed in these collisions, but they may not be at the same level of thermalization with the medium as observed at $\sqrt{s_{_{NN}}}=200$ GeV.


title: Computational homogenization of fibrous piezoelectric materials

abstract: Flexible piezoelectric devices made of polymeric materials are widely used for micro- and nano-electro-mechanical systems. In particular, numerous recent applications concern energy harvesting. Due to the importance of computational modeling to understand the influence that microscale geometry and constitutive variables exert on the macroscopic behavior, a numerical approach is developed here for multiscale and multiphysics modeling of thin piezoelectric sheets made of aligned arrays of polymeric nanofibers, manufactured by electrospinning. At the microscale, the representative volume element consists in piezoelectric polymeric nanofibers, assumed to feature a piezoelastic behavior and subjected to electromechanical contact constraints. The latter are incorporated into the virtual work equations by formulating suitable electric, mechanical and coupling potentials and the constraints are enforced by using the penalty method. From the solution of the micro-scale boundary value problem, a suitable scale transition procedure leads to identifying the performance of a macroscopic thin piezoelectric shell element.


title: Mass and magnification maps for the Hubble Space Telescope Frontier Fields clusters: implications for high redshift studies

abstract: Extending over three Hubble Space Telescope (HST) cycles, the Hubble Frontier Fields (HFF) initiative constitutes the largest commitment ever of HST time to the exploration of the distant Universe via gravitational lensing by massive galaxy clusters. We here present models of the mass distribution in the six HFF cluster lenses, derived from a joint strong- and weak-lensing analysis anchored by a total of 88 multiple-image systems identified in existing HST data. The resulting maps of the projected mass distribution and of the gravitational magnification effectively calibrate the HFF clusters as gravitational telescopes. Allowing the computation of search areas in the source plane, these maps are provided to the community to facilitate the exploitation of forthcoming HFF data for quantitative studies of the gravitationally lensed population of background galaxies. Our models of the gravitational magnification afforded by the HFF clusters allow us to quantify the lensing-induced boost in sensitivity over blank-field observations and predict that galaxies at $z>10$ and as faint as m(AB)=32 will be detectable, up to 2 magnitudes fainter than the limit of the Hubble Ultra Deep Field.


title: The emergence of classical behavior in magnetic adatoms

abstract: A wide class of nanomagnets shows striking quantum behavior, known as quantum spin tunneling (QST): instead of two degenerate ground states with opposite magnetizations, a bonding-antibonding pair forms, resulting in a splitting of the ground state doublet with wave functions linear combination of two classically opposite magnetic states, leading to the quenching of their magnetic moment. Here we study how QST is destroyed and classical behavior emerges in the case of magnetic adatoms, as the strength of their coupling, either to the substrate or to each other, is increased. Both spin-substrate and spin-spin coupling renormalize the QST splitting to zero allowing the environmental decoherence to eliminate superpositions between classical states, leading to the emergence of spontaneous magnetization.


title: The obstacle-mass constraint problem for hyperbolic conservation laws. Solvability

abstract: In this work we introduce the obstacle-mass constraint problem for a multidimensional scalar hyperbolic conservation law. We prove existence of an entropy solution to this problem by a penalization/viscosity method. The mass constraint introduces a nonlocal Lagrange multiplier in the penalized equation, giving rise to a nonlocal parabolic problem. We determine conditions on the initial data and on the obstacle function which ensure global in time existence of solution. These are not smoothness conditions, but relate to the propagation of the support of the initial data.


title: Ehrenfest Scheme of Higher Dimensional Topological AdS Black Holes in The Third Order Lovelock-Born-Infeld Gravity

abstract: Interpreting the cosmological constant as a thermodynamic pressure and its conjugate quantity as a thermodynamic volume, we reconsider the investigation of P-V critical behaviors of (1+n)-dimensional topological AdS black holes in Lovelock-Born-Infeld gravity. In particular, we give an explicit expression of the universal number \chi=\frac{P_c v_c}{T_c} in terms of the space dimension $n$. Then, we examine the phase transitions at the critical points of such topological black holes for 6 \leq n \leq 11 as required by the physical condition of the thermodynamical quantities. More precisely, the Ehrenfest equations have been checked revealing that the black hole system undergoes a second phase transition at the critical points.


title: Interface Problems for Dispersive equations

abstract: The interface problem for the linear Schr\"odinger equation in one-dimensional piecewise homogeneous domains is examined by providing an explicit solution in each domain. The location of the interfaces is known and the continuity of the wave function and a jump in their derivative at the interface are the only conditions imposed. The problem of two semi-infinite domains and that of two finite-sized domains are examined in detail. The problem and the method considered here extend that of an earlier paper by Deconinck, Pelloni and Sheils (2014). The dispersive nature of the problem presents additional difficulties that are addressed here.


title: The Comparison of Physical Properties Derived from Gas and Dust in a Massive Star-Forming Region

abstract: We explore the relationship between gas and dust in massive star-forming regions by comparing physical properties derived from each. We compare the temperatures and column densities in a massive star-forming Infrared Dark Cloud (IRDC, G32.02+0.05), which shows a range of evolutionary states, from quiescent to active. The gas properties were derived using radiative transfer modeling of the (1,1), (2,2), and (4,4) transitions of NH3 on the Karl G. Jansky Very Large Array (VLA), while the dust temperatures and column densities were calculated using cirrus-subtracted, modified blackbody fits to Herschel data. We compare the derived column densities to calculate an NH3 abundance, 4.6 x 10^-8. In the coldest star-forming region, we find that the measured dust temperatures are lower than the measured gas temperatures (mean and standard deviations T_dust ~ 11.6 +/- 0.2 K vs. T_gas ~ 15.2 +/- 1.5 K), which may indicate that the gas and dust are not well-coupled in the youngest regions (~0.5 Myr) or that these observations probe a regime where the dust and/or gas temperature measurements are unreliable. Finally, we calculate millimeter fluxes based on the temperatures and column densities derived from NH3 which suggest that millimeter dust continuum observations of massive star-forming regions, such as the Bolocam Galactic Plane Survey or ATLASGAL, can probe hot cores, cold cores, and the dense gas lanes from which they form, and are generally not dominated by the hottest core.


title: The effect of quenched bond disorder on first-order phase transitions

abstract: We investigate the effect of quenched bond disorder on the two-dimensional three-color Ashkin-Teller model, which undergoes a first-order phase transition in the absence of impurities. This is one of the simplest and striking models in which quantitative numerical simulations can be carried out to investigate emergent criticality due to disorder rounding of first-order transition. Utilizing extensive cluster Monte Carlo simulations on large lattice sizes of up to $128 \times 128$ spins, each of which is represented by three colors taking values $\pm 1$, we show that the rounding of the first-order phase transition is an emergent criticality. We further calculate the correlation length critical exponent, $\nu$, and the magnetization critical exponent, $\beta$, from finite size scaling analysis. We find that the critical exponents, $\nu$ and $\beta$, change as the strength of disorder or the four-spin coupling varies, and we show that the critical exponents appear not to be in the Ising universality class. We know of no analytical approaches that can explain our non-perturbative results. However our results should inspire further work on this important problem, either numerical or analytical.


title: Pathways for virus assembly around nucleic acids

abstract: Understanding the pathways by which viral capsid proteins assemble around their genomes could identify key intermediates as potential drug targets. In this work we use computer simulations to characterize assembly over a wide range of capsid protein-protein interaction strengths and solution ionic strengths. We find that assembly pathways can be categorized into two classes, in which intermediates are either predominantly ordered or disordered. Our results suggest that estimating the protein-protein and the protein-genome binding affinities may be sufficient to predict which pathway occurs. Furthermore, the calculated phase diagrams suggest that knowledge of the dominant assembly pathway and its relationship to control parameters could identify optimal strategies to thwart or redirect assembly to block infection. Finally, analysis of simulation trajectories suggests that the two classes of assembly pathways can be distinguished in single molecule fluorescence correlation spectroscopy or bulk time resolved small angle x-ray scattering experiments.


title: Beta Reduction is Invariant, Indeed (Long Version)

abstract: Slot and van Emde Boas' weak invariance thesis states that reasonable machines can simulate each other within a polynomially overhead in time. Is $\lambda$-calculus a reasonable machine? Is there a way to measure the computational complexity of a $\lambda$-term? This paper presents the first complete positive answer to this long-standing problem. Moreover, our answer is completely machine-independent and based over a standard notion in the theory of $\lambda$-calculus: the length of a leftmost-outermost derivation to normal form is an invariant cost model. Such a theorem cannot be proved by directly relating $\lambda$-calculus with Turing machines or random access machines, because of the size explosion problem: there are terms that in a linear number of steps produce an exponentially long output. The first step towards the solution is to shift to a notion of evaluation for which the length and the size of the output are linearly related. This is done by adopting the linear substitution calculus (LSC), a calculus of explicit substitutions modelled after linear logic and proof-nets and admitting a decomposition of leftmost-outermost derivations with the desired property. Thus, the LSC is invariant with respect to, say, random access machines. The second step is to show that LSC is invariant with respect to the $\lambda$-calculus. The size explosion problem seems to imply that this is not possible: having the same notions of normal form, evaluation in the LSC is exponentially longer than in the $\lambda$-calculus. We solve such an impasse by introducing a new form of shared normal form and shared reduction, deemed useful. Useful evaluation avoids those steps that only unshare the output without contributing to $\beta$-redexes, i.e., the steps that cause the blow-up in size.


title: Alexandrov spaces with large volume growth

abstract: Let $(X,d)$ be an $n$-dimensional Alexandrov space whose Hausdorff measure $\mathcal{H}^n$ satisfies a condition giving the metric measure space $(X,d,\mathcal{H}^n)$ a notion of having nonnegative Ricci curvature. We examine the influence of large volume growth on these spaces and generalize some classical arguments from Riemannian geometry showing that when the volume growth is sufficiently large, then $(X,d,\mathcal{H}^n)$ has finite topological type.


title: Intrinsic flat convergence with bounded Ricci curvature

abstract: In this paper we address the relationship between Gromov-Hausdorff limits and intrinsic flat limits of complete Riemannian manifolds. In \cite{SormaniWenger2010, SormaniWenger2011}, Sormani-Wenger show that for a sequence of Riemannian manifolds with nonnegative Ricci curvature, a uniform upper bound on diameter, and non-collapsed volume, the intrinsic flat limit exists and agrees with the Gromov-Hausdorff limit. This can be viewed as a non-cancellation theorem showing that for such sequences, points don't cancel each other out in the limit. Here we prove a similar no-cancellation theorem, replacing the assumption of nonnegative Ricci curvature with a two-sided bound on Ricci curvature. This version corrects a mistake in the previous version of this paper (where we assume only an arbitrary lower Ricci bound) which was due to a crucial error in one of our supporting theorems for that argument.


title: Linked symplectic forms and limit linear series in rank 2 with special determinant

abstract: We generalize the prior linked symplectic Grassmannian construction, applying it to to prove smoothing results for rank-2 limit linear series with fixed special determinant on chains of curves. We apply this general machinery to prove new results on nonemptiness and dimension of rank-2 Brill-Noether loci in a range of degrees.


title: Designing anti-cancer drugs and directing anti-cancer therapy

abstract: A prototype for a web application was designed and implemented as a guide to be used by clinicians when designing the best drug therapy for a specific cancer patient, given biological data derived from the patients tumor tissue biopsy. A representation of the patients metabolic pathways is displayed as a graph in the application, with nodes as substrates and products and edges as enzymes. The top metabolically active sub- paths in the pathway, ranked using an algorithm based on both the patients biological data and the graph topology, are also displayed and can be individually highlighted to examine potential enzymatic sites to be disrupted by a drug.


title: Optimal Exploration-Exploitation in a Multi-Armed-Bandit Problem with Non-stationary Rewards

abstract: In a multi-armed bandit (MAB) problem a gambler needs to choose at each round of play one of K arms, each characterized by an unknown reward distribution. Reward realizations are only observed when an arm is selected, and the gambler's objective is to maximize his cumulative expected earnings over some given horizon of play T. To do this, the gambler needs to acquire information about arms (exploration) while simultaneously optimizing immediate rewards (exploitation); the price paid due to this trade off is often referred to as the regret, and the main question is how small can this price be as a function of the horizon length T. This problem has been studied extensively when the reward distributions do not change over time; an assumption that supports a sharp characterization of the regret, yet is often violated in practical settings. In this paper, we focus on a MAB formulation which allows for a broad range of temporal uncertainties in the rewards, while still maintaining mathematical tractability. We fully characterize the (regret) complexity of this class of MAB problems by establishing a direct link between the extent of allowable reward "variation" and the minimal achievable regret. Our analysis draws some connections between two rather disparate strands of literature: the adversarial and the stochastic MAB frameworks.


title: The effect of radial migration on galactic disks

abstract: We study the radial migration of stars driven by recurring multi-arm spiral features in an exponential disk embedded in a dark matter halo. The spiral perturbations redistribute angular momentum within the disk and lead to substantial radial displacements of individual stars, in a manner that largely preserves the circularity of their orbits and that results, after 5 Gyr (~40 full rotations at the disk scalelength), in little radial heating and no appreciable changes to the vertical or radial structure of the disk. Our results clarify a number of issues related to the spatial distribution and kinematics of migrators. In particular, we find that migrators are a heavily biased subset of stars with preferentially low vertical velocity dispersions. This "provenance bias" for migrators is not surprising in hindsight, for stars with small vertical excursions spend more time near the disk plane and thus respond more readily to non-axisymmetric perturbations. We also find that the vertical velocity dispersion of outward migrators always decreases, whereas the opposite holds for inward migrators. To first order, newly arrived migrators simply replace stars that have migrated off to other radii, thus inheriting the vertical bias of the latter. Extreme migrators might therefore be recognized, if present, by the unexpectedly small amplitude of their vertical excursions. Our results show that migration, understood as changes in angular momentum that preserve circularity, can affect strongly the thin disk, but cast doubts on models that envision the Galactic thick disk as a relic of radial migration.


title: Adaptive Monte Carlo via Bandit Allocation

abstract: We consider the problem of sequentially choosing between a set of unbiased Monte Carlo estimators to minimize the mean-squared-error (MSE) of a final combined estimate. By reducing this task to a stochastic multi-armed bandit problem, we show that well developed allocation strategies can be used to achieve an MSE that approaches that of the best estimator chosen in retrospect. We then extend these developments to a scenario where alternative estimators have different, possibly stochastic costs. The outcome is a new set of adaptive Monte Carlo strategies that provide stronger guarantees than previous approaches while offering practical advantages.


title: Fully Bayesian Logistic Regression with Hyper-Lasso Priors for High-dimensional Feature Selection

abstract: High-dimensional feature selection arises in many areas of modern science. For example, in genomic research we want to find the genes that can be used to separate tissues of different classes (e.g. cancer and normal) from tens of thousands of genes that are active (expressed) in certain tissue cells. To this end, we wish to fit regression and classification models with a large number of features (also called variables, predictors). In the past decade, penalized likelihood methods for fitting regression models based on hyper-LASSO penalization have received increasing attention in the literature. However, fully Bayesian methods that use Markov chain Monte Carlo (MCMC) are still in lack of development in the literature. In this paper we introduce an MCMC (fully Bayesian) method for learning severely multi-modal posteriors of logistic regression models based on hyper-LASSO priors (non-convex penalties). Our MCMC algorithm uses Hamiltonian Monte Carlo in a restricted Gibbs sampling framework; we call our method Bayesian logistic regression with hyper-LASSO (BLRHL) priors. We have used simulation studies and real data analysis to demonstrate the superior performance of hyper-LASSO priors, and to investigate the issues of choosing heaviness and scale of hyper-LASSO priors.


title: Self-biased Reconfigurable Graphene Stacks for Terahertz Plasmonics

abstract: The gate-controllable complex conductivity of graphene offers unprecedented opportunities for reconfigurable plasmonics at THz and mid-IR frequencies. However, the requirement of a gating electrode close to graphene and the single `control knob' that this approach offers for graphene conductivity limits the practical implementation and performance of graphene-controllable plasmonic devices. Herein, we report on graphene stacks composed of two or more graphene monolayers separated by electrically thin dielectrics and present a simple and rigorous theoretical framework for their characterization. In a first implementation, two graphene layers gate each other, thereby behaving as a controllable single equivalent layer but without any additional gating structure. Second, we show that adding an additional gate --a third graphene layer or an external gate-- allows independent control of the complex conductivity of each layer within the stack and hence provides enhanced control on the stack equivalent complex conductivity. The proposed concepts are first theoretically studied and then demonstrated experimentally via a detailed procedure allowing extraction of the parameters of each layer independently and for arbitrary pre-doping. These results are believed to be instrumental to the development of THz and mid-IR plasmonic devices with enhanced performance and reconfiguration capabilities.


title: Natural scalars in the NMSSM

abstract: In the motivated hypothesis that the scalar bosons of the Next-to-Minimal Supersymmetric Standard Model (NMSSM) be the lightest new particles around, a possible strategy to search for signs of the extra CP-even states is outlined. It is shown how the measurements of the couplings of the 126 GeV Higgs boson constrain the region of the physical parameters in a generic NMSSM which minimises the fine-tuning of the electroweak scale. We also determine the cross section for the production of a heavier CP-even scalar, together with its most relevant branching ratios.


title: Inapproximability of Nash Equilibrium

abstract: We prove that finding an $\epsilon$-approximate Nash equilibrium is PPAD-complete for constant $\epsilon$ and a particularly simple class of games: polymatrix, degree 3 graphical games, in which each player has only two actions. As corollaries, we also prove similar inapproximability results for Bayesian Nash equilibrium in a two-player incomplete information game with a constant number of actions, for relative $\epsilon$-Well Supported Nash Equilibrium in a two-player game, for market equilibrium in a non-monotone market, for the generalized circuit problem defined by Chen, Deng, and Teng [CDT'09], and for approximate competitive equilibrium from equal incomes with indivisible goods.


title: Large Code Base Change Ripple Management in C++: My thoughts on how a new Boost C++ Library could help

abstract: C++ 98/03 already has a reputation for overwhelming complexity compared to other programming languages. The raft of new features in C++ 11/14 suggests that the complexity in the next generation of C++ code bases will overwhelm still further. The planned C++ 17 will probably worsen matters in ways difficult to presently imagine. Countervailing against this rise in software complexity is the hard de-exponentialisation of computer hardware capacity growth expected no later than 2020, and which will have even harder to imagine consequences on all computer software. WG21 C++ 17 study groups SG2 (Modules), SG7 (Reflection), SG8 (Concepts), and to a lesser extent SG10 (Feature Test) and SG12 (Undefined Behaviour), are all fundamentally about significantly improving complexity management in C++ 17, yet WG21's significant work on improving C++ complexity management is rarely mentioned explicitly. This presentation pitches a novel implementation solution for some of these complexity scaling problems, tying together SG2 and SG7 with parts of SG3 (Filesystem): a standardised but very lightweight transactional graph database based on Boost.ASIO, Boost.AFIO and Boost.Graph at the very core of the C++ runtime, making future C++ codebases considerably more tractable and affordable to all users of C++.


title: Representations of the alternating group which are irreducible over subgroups. II

abstract: We prove that non-trivial representations of the alternating group $A_n$ are reducible over a primitive proper subgroup which is isomorphic to some alternating group $A_m$.


title: Plasmon-assisted quantum control of distant emitters

abstract: We show how to generate and control the correlations in a set of two distant quantum emitters coupled to a one-dimensional dissipative plasmonic waveguide. An external laser field enhances the dimer's steady-state correlations and allows an active control (switching on/off) of nonclassical correlations. The plasmon-assisted dipolar-interacting qubits exhibit persistent correlations, which in turn can be decoupled and made to evolve independently from each other. The setup enables long-distance ($\sim 1\, \mu$m) qubit control that works for both resonant and detuned emitters. For suitable emitter initialization, we also show that the quantum correlation is always greater than the classical one.


title: Modular Representation Theory of Symmetric Groups

abstract: We review some recent advances in modular representation theory of symmetric groups and related Hecke algebras. We discuss connections with Khovanov-Lauda-Rouquier algebras and gradings on the blocks of the group algebras $F\Sigma_n$, which these connections reveal; graded categorification and connections with quantum groups and crystal bases; modular branching rules and the Mullineaux map; graded cellular structure and graded Specht modules; cuspidal systems for affine KLR algebras and imaginary Schur-Weyl duality, which connects representation theory of these algebras to the usual Schur algebras of smaller rank.


title: Hydrodynamic limit for conservative spin systems with super-quadratic, partially inhomogeneous single-site potential

abstract: We consider an interacting unbounded spin system, with conservation of the mean spin. We derive quantitative rates of convergence to the hydrodynamic limit provided the single-site potential is a bounded perturbation of a strictly convex function with polynomial growth, and with an additional random inhomogeneous linear term. This additional linear term models the impact of a random chemical potential. The argument adapts the two-scale approach of Grunewald, Otto, Villani and Westdickenberg from the quadratic to the general case. The main ingredient is the derivation of a covariance estimate that is uniform in the system size. We also show that this covariance estimate can be used to change the iterative argument of [MO] for deducing the optimal scaling LSI for the canonical ensemble into a two-scale argument in the sense of [GOVW]. We also prove the LSI for canonical ensembles with an inhomogeneous linear term.


title: Affine highest weight categories and affine quasihereditary algebras

abstract: Koenig and Xi introduced {\em affine cellular algebras}. Kleshchev and Loubert showed that an important class of {\em infinite dimensional} algebras, the KLR algebras $R(\Gamma)$ of finite Lie type $\Gamma$, are (graded) affine cellular; in fact, the corresponding affine cell ideals are idempotent. This additional property is reminiscent of the properties of {\em quasihereditary algebras} of Cline-Parshall-Scott in a {\em finite dimensional} situation. A fundamental result of Cline-Parshall-Scott says that a finite dimensional algebra $A$ is quasihereditary if and only if the category of finite dimensional $A$-modules is a {\em highest weight category}. On the other hand, S. Kato and Brundan-Kleshchev-McNamara proved that the category of {\em finitely generated graded} $R(\Gamma)$-modules has many features reminiscent of those of a highest weight category. The goal of this paper is to axiomatize and study the notions of an {\em affine quasihereditary algebra} and an {\em affine highest weight category}. In particular, we prove an affine analogue of the Cline-Parshall-Scott Theorem. We also develop {\em stratified} versions of these notions.


title: The Dirichlet problem for elliptic systems with data in K\"othe function spaces

abstract: We show that the boundedness of the Hardy-Littlewood maximal operator on a K\"othe function space ${\mathbb{X}}$ and on its K\"othe dual ${\mathbb{X}}'$ is equivalent to the well-posedness of the $\mathbb{X}$-Dirichlet and $\mathbb{X}'$-Dirichlet problems in $\mathbb{R}^{n}_{+}$ in the class of all second-order, homogeneous, elliptic systems, with constant complex coefficients. As a consequence, we obtain that the Dirichlet problem for such systems is well-posed for boundary data in Lebesgue spaces, variable exponent Lebesgue spaces, Lorentz spaces, Zygmund spaces, as well as their weighted versions. We also discuss a version of the aforementioned result which contains, as a particular case, the Dirichlet problem for elliptic systems with data in the classical Hardy space $H^1$, and the Beurling-Hardy space ${\rm HA}^p$ for $p\in(1,\infty)$. Based on the well-posedness of the $L^p$-Dirichlet problem we then prove the uniqueness of the Poisson kernel associated with such systems, as well as the fact that they generate a strongly continuous semigroup in natural settings. Finally, we establish a general Fatou type theorem guaranteeing the existence of the pointwise nontangential boundary trace for null-solutions of such systems.


title: Containment: A Variation of Cops and Robbers

abstract: We consider "Containment": a variation of the graph pursuit game of Cops and Robber in which cops move from edge to adjacent edge, the robber moves from vertex to adjacent vertex (but cannot move along an edge occupied by a cop), and the cops win by "containing" the robber---that is, by occupying all $\deg(v)$ of the edges incident with a vertex $v$ while the robber is at $v$. We develop bounds that relate the minimal number of cops, $\xi(G)$, required to contain a robber to the well-known "cop-number" $c(G)$ in the original game: in particular, $c(G) {\le} \xi(G) {\le} \gamma(G) \Delta(G)$. We note that $\xi(G) {\geq} \delta(G)$ for all graphs $G$, and analyze several families of graphs in which equality holds, as well as several in which the inequality is strict. We also give examples of graphs which require an unbounded number of cops in order to contain a robber, and note that there exist cubic graphs with $\xi(G) \geq \Omega(n^{1/6})$.


title: Graphic Processing Unit Simulation of Axon Growth and Guidance through Cue Diffusion on Massively Parallel Processors

abstract: Neural development represents not only an exciting and complex field of study, with ongoing progress, but it also became the epicentre of neuroscience and developmental biology, as it strives to describe the underlying cellular and molecular mechanisms by which the central nervous system emerges during the various levels of embryonic development phases. The nervous system is a dynamic entity, where the genetic information plays an important role in shaping the intra- and extracellular environments, which in turn offer a reliable foundation for the stem cell precursors to divide and form neurons. Throughout the embryonic development stages, the neurons undergo different processes: migration at an immature level from the initial place in the embryo to a predefined final position, axonal differentiation and guidance of the motile growth cone towards a postsynaptic target, synaptic formation between axons and target, and lastly long-term synaptic changes which underlie learning and memory. In order to gain a better understanding of how the nervous system develops, mathematical and computational models have been created and expanded in order to bridge the gap between system-level dynamics and lower level cellular and molecular processes. This research paper aims to illustrate the potential of theoretical mathematical and computational models for analysing one important stage of neural development - axonal growth and guidance mechanisms in the presence of diffusion cues, through a visual simulation which is optimized via the graphic processing unit and parallel programming techniques.


title: Cobalt-56 gamma-ray emission lines from the type Ia supernova 2014J

abstract: A type Ia supernova is thought to be a thermonuclear explosion of either a single carbon-oxygen white dwarf or of a pair of merging white dwarfs. The explosion fuses a large amount of radioactive 56Ni. After the explosion, the decay chain from 56Ni to 56Co to 56Fe generates gamma-ray photons, which are reprocessed in the expanding ejecta and give rise to powerful optical emission. Here we report the detection of 56Co lines at energies of 847 and 1238 keV and a gamma-ray continuum in the 200-400 keV band from the type Ia supernova 2014J in the nearby galaxy M82. The line fluxes suggest that about 0.6 +/- 0.1 solar masses of radioactive 56Ni were synthesized during the explosion. The line broadening gives a characteristic mass-weighted ejecta expansion velocity of 10000 +/- 3000 km/s. The observed gamma-ray properties are in broad agreement with the canonical model of an explosion of a white dwarf just massive enough to be unstable to gravitational collapse, but do not immediately exclude more complicated merger scenarios, which fuse comparable amount of 56Ni.


title: Mapping neuron positions to curved cortical surfaces

abstract: Real cortical tissue curves and folds according to experimental data. However, our current simulations only use unfolded cortical layers. This project extends the cortical slice model in order to be able to specify arbitrary curvatures to the cortical layers and investigate the effect of curvature and folding on neuron density and layer morphology.


title: Common Connectome Constraints: From C. elegans and Drosophila to Homo sapiens

abstract: Neural systems show a modular and typically also a hierarchical organisation across different levels and across different species. Topology relates to function, but it is also influences dynamics as earlier studies showed its effect on synchrony, oscillation, and activity propagation. Understanding the link between the hierarchical organisation and processing (e.g. does consciousness structurally correlate with the top level of the hierarchy and where is the 'top' in a network?) remains one of the main challenges of the field. In addition, although neuron nodes are often treated as uniform entities, they can differ in terms of function (e.g. inhibitory vs. excitatory), morphology, or gene expression pattern.


title: A hierarchical network organization helps to retain comparable oscillation patterns in rats and human-sized brains

abstract: Activity in coupled systems is often oscillatory, for example, the firing pattern of neuronal populations. Whereas these oscillations have been studied predominantly in local circuits, here we show how the topology of large-scale networks, leading to large feedback loops, influences oscillations in the resting state. We find that the hierarchical modular organization of neuronal networks supports distinct spectral patterns of neural rhythms similar to those observed experimentally in different species such as rat and human. For individual neurons, multiple peak frequencies with non-integer ratios between subsequent peaks occurred. These ratios occurred both for models with the spatial size of the rat as well as for the human brain. We argue that the potential influence of longer connections, and thus longer delays, are balanced by a reduced number of long-distance connections in larger brain networks. In conclusion, we show that a hierarchical neuronal network provides a scalable backbone for multiple brain rhythms. This structural backbone could complement well-studied regional and cellular mechanisms for the generation or prevention of rhythms.


title: Auslander-Reiten quiver of type A and generalized quantum affine Schur-Weyl duality

abstract: The quiver Hecke algebra $R$ can be also understood as a generalization of the affine Hecke algebra of type $A$ in the context of the quantum affine Schur-Weyl duality by the results of Kang, Kashiwara and Kim. On the other hand, it is well-known that the Auslander-Reiten(AR) quivers $\Gamma_Q$ of finite simply-laced types have a deep relation with the positive roots systems of the corresponding types. In this paper, we present explicit combinatorial descriptions for the AR-quivers $\Gamma_Q$ of finite type $A$. Using the combinatorial descriptions, we can investigate relations between finite dimensional module categories over the quantum affine algebra $U'_q(A_n^{(i)})$ $(i=1,2)$ and finite dimensional graded module categories over the quiver Hecke algebra $R_{A_n}$ associated to $A_n$ through the generalized quantum affine Schur-Weyl duality functor.


title: Thermal effects on bipartite and multipartite correlations in fiber coupled cavity arrays

abstract: We investigate the thermal influence of fibers on the dynamics of bipartite and multipartite correlations in fiber coupled cavity arrays where each cavity is resonantly coupled to a two-level atom. The atom-cavity systems connected by fibers can be considered as polaritonic qubits. We first derive a master equation to describe the evolution of the atom-cavity systems. The bipartite (multipartite) correlations is measured by concurrence and discord (spin squeezing). Then, we solve the master equation numerically and study the thermal effects on the concurrence, discord, and spin squeezing of qubits. On the one hand, at zero temperature, there are steady-state bipartite and multipartite correlations. One the other hand, the thermal fluctuations of a fiber may blockade the generation of entanglement of two qubits connected directly by the fiber while the discord can be generated and stored for a long time. This thermal-induced blockade effects of bipartite correlations may be useful for quantum information processing. The bipartite correlations of a longer chain of qubits is more robust than a shorter one in the presence of thermal fluctuations.


title: Finite range and upper branch effects on itinerant ferromagnetism in repulsive Fermi gases: Bethe-Goldstone ladder resummation approach

abstract: We investigate the ferromagnetic transition in repulsive Fermi gases at zero temperature with upper branch and effective range effects. Based on a general effective Lagrangian that reproduces precisely the two-body $s$-wave scattering phase shift, we obtain a nonperturbative expression of the energy density as a function of the polarization by using the Bethe-Goldstone ladder resummation. For hard sphere potential, the predicted critical gas parameter $k_{\rm F}a=0.816$ and the spin susceptibility agree well with the results from fixed-node diffusion Monte Carlo calculations. In general, positive and negative effective ranges have opposite effects on the critical gas parameter $k_{\rm F}a$: While a positive effective range reduces the critical gas parameter, a negative effective range increases it. For attractive potential or Feshbach resonance model, the many-body upper branch exhibits an energy maximum at $k_{\rm F}a=\alpha$ with $\alpha=1.34$ from the Bethe-Goldstone ladder resummation, which is qualitatively consistent with experimental results. The many-body T-matrix has a positive-energy pole for $k_{\rm F}a>\alpha$ and it becomes impossible to distinguish the bound state and the scattering state. These positive-energy bound states become occupied and therefore the upper branch reaches an energy maximum at $k_{\rm F}a=\alpha$. In the zero range limit, there exists a narrow window ($0.86<k_{\rm F}a<1.56$) for the ferromagnetic phase. At sufficiently large negative effective range, the ferromagnetic phase disappears. On the other hand, the appearance of positive-energy bound state resonantly enhances the two-body decay rate around $k_{\rm F}a=\alpha$ and may prevent the study of equilibrium phases and ferromagnetism of the upper branch Fermi gas.


title: The Historic Set of Ergodic Averages in Some Nonuniformly Hyperbolic Systems

abstract: This article is devoted to the study of the historic set of ergodic averages in some nonuniformly hyperbolic systems. In particular, our results hold for the robust classes of multidimensional nonuniformly expanding local diffeomorphisms and Viana maps.


title: Post-selection point and interval estimation of signal sizes in Gaussian samples

abstract: We tackle the problem of the estimation of a vector of means from a single vector-valued observation $y$. Whereas previous work reduces the size of the estimates for the largest (absolute) sample elements via shrinkage (like James-Stein) or biases estimated via empirical Bayes methodology, we take a novel approach. We adapt recent developments by Lee et al (2013) in post selection inference for the Lasso to the orthogonal setting, where sample elements have different underlying signal sizes. This is exactly the setup encountered when estimating many means. It is shown that other selection procedures, like selecting the $K$ largest (absolute) sample elements and the Benjamini-Hochberg procedure, can be cast into their framework, allowing us to leverage their results. Point and interval estimates for signal sizes are proposed. These seem to perform quite well against competitors, both recent and more tenured. Furthermore, we prove an upper bound to the worst case risk of our estimator, when combined with the Benjamini-Hochberg procedure, and show that it is within a constant multiple of the minimax risk over a rich set of parameter spaces meant to evoke sparsity.


title: Existence of a regular unimodular triangulation of the edge polytopes of finite graphs

abstract: In this paper we give several criteria for the edge polytope of a fundamental FHM-graph to possess a regular unimodular triangulation in terms of some simple data of the the graph. We further apply our criteria to several examples of graphs and show that their edge polytopes possess a regular unimodular triangulation.


title: An Agent-based Modeling Framework for Sociotechnical Simulation of Water Distribution Contamination Events

abstract: In the event that a bacteriological or chemical toxin is intro- duced to a water distribution network, a large population of consumers may become exposed to the contaminant. A contamination event may be poorly predictable dynamic process due to the interactions of consumers and utility managers during an event. Consumers that become aware of a threat may select protective actions that change their water demands from typical demand patterns, and new hydraulic conditions can arise that differ from conditions that are predicted when demands are considered as exogenous inputs. Consequently, the movement of the contaminant plume in the pipe network may shift from its expected trajectory. A sociotechnical model is developed here to integrate agent-based models of consumers with an engineering water distribution system model and capture the dynamics between consumer behaviors and the water distribution system for predicting contaminant transport and public exposure. Consumers are simulated as agents with behaviors defined for water use activities, mobility, word-of-mouth communication, and demand reduction, based on a set of rules representing an agents autonomy and reaction to health impacts, the environment, and the actions of other agents. As consumers decrease their water use, the demand exerted on the water distribution system is updated; as the flow directions and volumes shift in response, the location of the contaminant plume is updated and the amount of contaminant consumed by each agent is calculated. The framework is tested through simulating realistic contamination scenarios for a virtual city and water distribution system.


title: Integrated silicon optomechanical transducers and their application in atomic force microscopy

abstract: This chapter describes the basic design, simulation, and fabrication for fully Si integrated, waveguide coupled, optomechanical force and displacement sensors. The approach of full Si integration of all stationary nanophotonic components with mechanically separated movable components creates the opportunity to independently engineer these two parts for a variety of MEMS and NEMS sensing applications that require high precision, high bandwidth, and small footprint. Further integration of actuators for static and dynamic actuation is also possible.


title: Regularisation Paths for Conditional Logistic Regression: the clogitL1 package

abstract: We apply the cyclic coordinate descent algorithm of Friedman, Hastie and Tibshirani (2010) to the fitting of a conditional logistic regression model with lasso ($\ell_1$) and elastic net penalties. The sequential strong rules of Tibshirani et al (2012) are also used in the algorithm and it is shown that these offer a considerable speed up over the standard coordinate descent algorithm with warm starts. Once implemented, the algorithm is used in simulation studies to compare the variable selection and prediction performance of the conditional logistic regression model against that of its unconditional (standard) counterpart. We find that the conditional model performs admirably on datasets drawn from a suitable conditional distribution, outperforming its unconditional counterpart at variable selection. The conditional model is also fit to a small real world dataset, demonstrating how we obtain regularisation paths for the parameters of the model and how we apply cross validation for this method where natural unconditional prediction rules are hard to come by.


title: Double series representations for Schur's partition function and related identities

abstract: We prove new double summation hypergeometric $q$-series representations for several families of partitions, including those that appear in the famous product identities of G\"ollnitz, Gordon, and Schur. We give several different proofs for our results, using bijective partitions mappings and modular diagrams, the theory of $q$-difference equations and recurrences, and the theories of summation and transformation for $q$-series. We also consider a general family of similar double series and highlight a number of other interesting special cases.


title: On the Cluster Physics of Sunyaev-Zel'dovich and X-ray Surveys IV: Characterizing Density and Pressure Clumping due to Infalling Substructures

abstract: Understanding the outskirts of galaxy clusters at the virial radius (R200) and beyond is critical for an accurate determination of cluster masses and to ensure unbiased cosmological parameter estimates from cluster surveys. This problem has drawn renewed interest due to recent determinations of gas mass fractions beyond R200, which appear to be considerably larger than the cosmic mean, and because the clusters' total Sunyaev-Zel'dovich flux receives a significant contribution from these regions. Here, we use a large suite of cosmological hydrodynamical simulations to study the clumpiness of density and pressure and employ different variants of simulated physics, including radiative gas physics and thermal feedback by active galactic nuclei. We find that density and pressure clumping closely trace each other as a function of radius, but the bias on density remains on average < 20% within the virial radius R200. At larger radius, clumping increases steeply due to the continuous infall of coherent structures that have not yet passed the accretion shock. Density and pressure clumping increase with cluster mass and redshift, which probes on average dynamically younger objects that are still in the process of assembling. The angular power spectra of gas density and pressure show that the clumping signal is dominated by comparably large substructures with scales >R200/5, signaling the presence of gravitationally-driven "super-clumping". In contrast, the angular power spectrum of the dark matter (DM) shows an almost uniform size distribution due to unimpeded subhalos. The quadrupolar anisotropy dominates the signal and correlates well across different radii as a result of the prolateness of the DM potential. We provide a synopsis of the radial dependence of the clusters' non-equilibrium measures (kinetic pressure support, ellipticity, and clumping) that all increase sharply beyond R200.


title: Self-dual codes and quadratic residue codes over the ring $\mathbb{Z}_9+u\mathbb{Z}_9$

abstract: In this paper, we introduce a new definitions of the Gray weight and the Gray map for linear codes over $\mathbb{Z}_9+u\mathbb{Z}_9$ with $u^2=u$. Some results on self-dual codes over this ring are investigated. Further, the structural properties of quadratic residue codes are also considered. Two self-dual codes with parameters $[22,11,5]$ and $[24,12,9]$ over $\mathbb{Z}_9$ are obtained.


title: Phosphorene nanoribbon as a promising candidate for thermoelectric applications

abstract: In this work, the electronic properties of phosphorene nanoribbons with different width and edge configurations are studied by using density functional theory. It is found that the armchair phosphorene nanoribbons are semiconducting while the zigzag nanoribbons are metallic. The band gaps of armchair nanoribbons decrease monotonically with increasing ribbon width. By passivating the edge phosphorus atoms with hydrogen, the zigzag series also become semiconducting, while the armchair series exhibit a larger band gap than their pristine counterpart. The electronic transport properties of these phosphorene nanoribbons are then investigated using Boltzmann theory and relaxation time approximation. We find that all the semiconducting nanoribbons exhibit very large values of Seebeck coefficient and can be further enhanced by hydrogen passivation at the edge. Taking armchair nanoribbon with width N=7 as an example, we calculate the lattice thermal conductivity with the help of phonon Boltzmann transport equation. Due to significantly enhanced Seebeck coefficient and decreased thermal conductivity, the phosphorene nanoribbon exhibit a very high figure of merit (ZT value) of 4.0 at room temperature, which suggests its appealing thermoelectric applications.


title: Ferromagnetism and quantum anomalous Hall effect in one-side-saturated buckled honeycomb lattices

abstract: The recently synthesized silicene as well as theoretically discussed germanene are examples of buckled honeycomb structures. The buckled structures allow one to manipulate asymmetry between two underlying sublattices of honeycomb structures. Here by taking germanene as a prototype of buckled honeycomb lattices, we explore magnetism induced by breaking sublattice symmetry through saturating chemical bonds on one-side of the buckled honeycomb lattice. It is shown that when fractions of chemical bonds on one-side are saturated, two narrow bands always exist at half filling. Furthermore, the narrow bands generally support flat band ferromagnetism in the presence of the Hubbard $U$ interaction. The induced magnetization is directly related to the saturation fraction and is thus controllable in magnitude through the saturation fraction. Most importantly, we find that depending on the saturation fraction, the ground state of an one-side saturated germanene may become a quantum anomalous Hall (QAH) insulator characterized by a Chern number that vanishes for larger magnetization. The non-vanishing Chern number for smaller magnetization implies that the associated quantum Hall effect tends to survive at high temperatures. Our findings provide a potential method to engineer buckled honeycomb structures into high-temperature QAH insulators.


title: Heat invariants of the perturbed polyharmonic Steklov problem

abstract: For a given bounded domain $\Omega$ with smooth boundary in a smooth Riemannian manifold $(\mathcal{M},g)$, we establish a procedure to get all the coefficients of the asymptotic expansion of the trace of the heat kernel associated with the perturbed polyharmonic Dirichlet-to-Neumann operator $\Lambda_m$ ($m\ge 1$) as $t\to 0^+$. We also explicitly calculate the first four coefficients of this asymptotic expansion. These coefficients (i.e., heat invariants) provide precise information for the area and curvatures of the boundary $\partial \Omega$ in terms of the spectrum of the perturbed polyharmonic Steklov problem. In particular, when $m=1$ and $q\equiv 0$ our work recovers the previous corresponding results in \cite{PS} and \cite{Liu3}.


title: Group-based Sparse Representation for Image Restoration

abstract: Traditional patch-based sparse representation modeling of natural images usually suffer from two problems. First, it has to solve a large-scale optimization problem with high computational complexity in dictionary learning. Second, each patch is considered independently in dictionary learning and sparse coding, which ignores the relationship among patches, resulting in inaccurate sparse coding coefficients. In this paper, instead of using patch as the basic unit of sparse representation, we exploit the concept of group as the basic unit of sparse representation, which is composed of nonlocal patches with similar structures, and establish a novel sparse representation modeling of natural images, called group-based sparse representation (GSR). The proposed GSR is able to sparsely represent natural images in the domain of group, which enforces the intrinsic local sparsity and nonlocal self-similarity of images simultaneously in a unified framework. Moreover, an effective self-adaptive dictionary learning method for each group with low complexity is designed, rather than dictionary learning from natural images. To make GSR tractable and robust, a split Bregman based technique is developed to solve the proposed GSR-driven minimization problem for image restoration efficiently. Extensive experiments on image inpainting, image deblurring and image compressive sensing recovery manifest that the proposed GSR modeling outperforms many current state-of-the-art schemes in both PSNR and visual perception.