title: Tidal debris morphology and the orbits of satellite galaxies

abstract: How do galaxies move relative to one another? While we can examine the motion of dark matter subhalos around their hosts in simulations of structure formation, determining the orbits of satellites around their parent galaxies from observations is impossible except for a small number of nearby cases. In this work we outline a novel approach to probing the orbital distributions of infalling satellite galaxies using the morphology of tidal debris structures. It has long been understood that the destruction of satellites on near-radial orbits tends to lead to the formation of shells of debris, while those on less eccentric orbits produce tidal streams. We combine an understanding of the scaling relations governing the orbital properties of debris with a simple model of how these orbits phase-mix over time to produce a `morphology metric' that more rigorously quantifies the conditions required for shells to be apparent in debris structures as a function of the satellite's mass and orbit and the interaction time. Using this metric we demonstrate how differences in orbit distributions can alter the relative frequency of shells and stream structures observed around galaxies. These experiments suggest that more detailed modeling and careful comparisons with current and future surveys of low surface brightness features around nearby galaxies should be capable of actually constraining orbital distributions and provide new insights into our understanding of structure formation.


title: Super-Eddington Stellar Winds Driven by Near-Surface Energy Deposition

abstract: We develop analytic and numerical models of the properties of super-Eddington stellar winds, motivated by phases in stellar evolution when super-Eddington energy deposition (via, e.g., unstable fusion, wave heating, or a binary companion) heats a region near the stellar surface. This appears to occur in luminous blue variables (LBVs), Type IIn supernovae progenitors, classical novae, and X-ray bursts. We show that when the wind kinetic power exceeds Eddington, the photons are trapped and behave like a fluid. Convection does not play a significant role in the wind energy transport. The wind properties depend on the ratio of a characteristic speed in the problem vc ~ (Edot G)^{1/5} (where Edot is the heating rate) to the stellar escape speed near the heating region vesc(r_h). For vc > vesc(r_h) the wind kinetic power at large radii Edot_w ~ Edot. For vc < vesc(r_h), most of the energy is used to unbind the wind material and thus Edot_w < Edot. Multidimensional hydrodynamic simulations without radiation diffusion using FLASH and one-dimensional hydrodynamic simulations with radiation diffusion using MESA are in good agreement with the analytic predictions. The photon luminosity from the wind is itself super-Eddington but in many cases the photon luminosity is likely dominated by `internal shocks' in the wind. We discuss the application of our models to eruptive mass loss from massive stars and argue that the wind models described here can account for the broad properties of LBV outflows and the enhanced mass loss in the years prior to Type IIn core-collapse supernovae.


title: The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the correlation function of LOWZ and CMASS galaxies in Data Release 12

abstract: We present distance scale measurements from the baryon acoustic oscillation signal in the CMASS and LOWZ samples from the Data Release 12 of the Baryon Oscillation Spectroscopic Survey (BOSS). The total volume probed is 14.5 Gpc$^3$, a 10 per cent increment from Data Release 11. From an analysis of the spherically averaged correlation function, we infer a distance to $z=0.57$ of $D_V(z)r^{\rm fid}_{\rm d}/r_ {\rm d}=2028\pm21$ Mpc and a distance to $z=0.32$ of $D_V(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1264\pm22$ Mpc assuming a cosmology in which $r^{\rm fid}_{\rm d}=147.10$ Mpc. From the anisotropic analysis, we find an angular diameter distance to $z=0.57$ of $D_{\rm A}(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1401\pm21$ Mpc and a distance to $z=0.32$ of $981\pm20$ Mpc, a 1.5 per cent and 2.0 per cent measurement respectively. The Hubble parameter at $z=0.57$ is $H(z)r_{\rm d}/r^{\rm fid}_{\rm d}=100.3\pm3.7$ km s$^{-1}$ Mpc$^{-1}$ and its value at $z=0.32$ is $79.2\pm5.6$ km s$^{-1}$ Mpc$^{-1}$, a 3.7 per cent and 7.1 per cent measurement respectively. These cosmic distance scale constraints are in excellent agreement with a $\Lambda$CDM model with cosmological parameters released by the recent Planck 2015 results.


title: The Baryon Cycle at High Redshifts: Effects of Galactic Winds on Galaxy Evolution in Overdense and Average Regions

abstract: We employ high-resolution cosmological zoom-in simulations focusing on a high-sigma peak and an average cosmological field at $z\sim 6-12$, in order to investigate the influence of environment and baryonic feedback on galaxy evolution in the reionization epoch. Strong feedback, e.g., galactic winds, caused by elevated star formation rates (SFRs) is expected to play an important role in this evolution. We compare different outflow prescriptions: (i) constant wind velocity (CW), (ii) variable wind scaling with galaxy properties (VW), and (iii) no outflows (NW). The overdensity leads to accelerated evolution of dark matter and baryonic structures, absent in the "normal" region, and to shallow galaxy stellar mass functions at the low-mass end. Although CW shows little dependence on both environments, the more physically motivated VW model does exhibit this effect. In addition, VW can reproduce the observed specific SFR (sSFR) and the sSFR-stellar mass relation, which CW and NW fail to satisfy simultaneously. Winds also differ substantially in affecting the state of the intergalactic medium (IGM). The difference lies in volume-filling factor of hot, high-metallicity gas which is near unity for CW, while it remains confined in massive filaments for VW, and locked up in galaxies for NW. Such gas is nearly absent in the normal region. Although all wind models suffer from deficiencies, the VW model seems to be promising in correlating the outflow properties to those of host galaxies. Further constraints on the state of the IGM at high-$z$ are needed to separate different wind models.


title: The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: BAO measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies

abstract: [abridged] We present an anisotropic analysis of the baryonic acoustic oscillation (BAO) scale in the twelfth and final data release of the Baryonic Oscillation Spectroscopic Survey (BOSS). We independently analyse the LOWZ and CMASS galaxy samples: the LOWZ sample contains contains 361 762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$; the CMASS sample consists of 777 202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. We extract the BAO peak position from the monopole power spectrum moment, $\alpha_0$, and from the $\mu^2$ moment, $\alpha_2$, where $\mu$ is the cosine of the angle to the line-of-sight. The $\mu^2$-moment provides equivalent information to that available in the quadrupole but is simpler to analyse. After applying a reconstruction algorithm to reduce the BAO suppression by bulk motions, we measure the BAO peak position in the monopole and $\mu^2$-moment, which are related to radial and angular shifts in scale. We report $H(z_{\rm LOWZ})r_s(z_d)=(11.60\pm0.60)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.66\pm0.16$ with a cross-correlation coefficient of $r_{HD_A}=0.41$, for the LOWZ sample; and $H(z_{\rm CMASS})r_s(z_d)=(14.56\pm0.37)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.42\pm0.13$ with a cross-correlation coefficient of $r_{HD_A}=0.47$, for the CMASS sample. We combine these results with the measurements of the BAO peak position in the monopole and quadrupole correlation function of the same dataset \citep[][companion paper]{Cuestaetal2015} and report the consensus values: $H(z_{\rm LOWZ})r_s(z_d)=(11.63\pm0.69)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.67\pm0.15$ with $r_{HD_A}=0.35$ for the LOWZ sample; $H(z_{\rm CMASS})r_s(z_d)=(14.67\pm0.42)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.47\pm0.12$ with $r_{HD_A}=0.52$ for the CMASS sample.


title: Sharpening the Weak Gravity Conjecture with Dimensional Reduction

abstract: We investigate the behavior of the Weak Gravity Conjecture (WGC) under toroidal compactification and RG flows, finding evidence that WGC bounds for single photons become weaker in the infrared. By contrast, we find that a photon satisfying the WGC will not necessarily satisfy it after toroidal compactification when black holes charged under the Kaluza-Klein photons are considered. Doing so either requires an infinite number of states of different charges to satisfy the WGC in the original theory or a restriction on allowed compactification radii. These subtleties suggest that if the Weak Gravity Conjecture is true, we must seek a stronger form of the conjecture that is robust under compactification. We propose a "Lattice Weak Gravity Conjecture" that meets this requirement: a superextremal particle should exist for every charge in the charge lattice. The perturbative heterotic string satisfies this conjecture. We also use compactification to explore the extent to which the WGC applies to axions. We argue that gravitational instanton solutions in theories of axions coupled to dilaton-like fields are analogous to extremal black holes, motivating a WGC for axions. This is further supported by a match between the instanton action and that of wrapped black branes in a higher-dimensional UV completion.


title: Hadron Collider Tests of Neutrino Mass-Generating Mechanisms

abstract: The Standard Model of particle physics (SM) is presently the best description of nature at small distances and high energies. However, with tiny but nonzero neutrino masses, a Higgs boson mass unstable under radiative corrections, and little guidance on understanding the hierarchy of fermion masses, the SM remains an unsatisfactory description of nature. Well-motivated scenarios that resolve these issues exist but also predict extended gauge (e.g., Left-Right Symmetric Models), scalar (e.g., Supersymmetry), and/or fermion sectors (e.g., Seesaw Models). Hence, discovering such new states would have far-reaching implications. After reviewing basic tenets of the SM and collider physics, several beyond the SM (BSM) scenarios that alleviate these shortcomings are investigated. Emphasis is placed on the production of a heavy Majorana neutrinos at hadron colliders in the context of low-energy, effective theories that simultaneously explain the origin of neutrino masses and their smallness compared to other elementary fermions, the so-called Seesaw Mechanisms. As probes of new physics, rare top quark decays to Higgs bosons in the context of the SM, the Types I and II Two Higgs Doublet Model (2HDM), and the semi-model independent framework of Effective Field Theory (EFT) have also been investigated. Observation prospects and discovery potentials of these models at current and future collider experiments are quantified.


title: Detecting Effects of Filaments on Galaxy Properties in the Sloan Digital Sky Survey III

abstract: We study the effects of filaments on galaxy properties in the Sloan Digital Sky Survey (SDSS) Data Release 12 using filaments from the `Cosmic Web Reconstruction' catalogue (Chen et al. 2016), a publicly available filament catalogue for SDSS. Since filaments are tracers of medium-to-high density regions, we expect that galaxy properties associated with the environment are dependent on the distance to the nearest filament. Our analysis demonstrates that a red galaxy or a high-mass galaxy tend to reside closer to filaments than a blue or low-mass galaxy. After adjusting the effect from stellar mass, on average, early-forming galaxies or large galaxies have a shorter distance to filaments than late-forming galaxies or small galaxies. For the Main galaxy sample (MGS), all signals are very significant ($>6\sigma$). For the LOWZ and CMASS sample, the stellar mass and size are significant ($>2 \sigma$). The filament effects we observe persist until $z = 0.7$ (the edge of the CMASS sample). Comparing our results to those using the galaxy distances from redMaPPer galaxy clusters as a reference, we find a similar result between filaments and clusters. Moreover, we find that the effect of clusters on the stellar mass of nearby galaxies depends on the galaxy's filamentary environment. Our findings illustrate the strong correlation of galaxy properties with proximity to density ridges, strongly supporting the claim that density ridges are good tracers of filaments.


title: The Formation of Submillimetre-Bright Galaxies from Gas Infall over a Billion Years

abstract: Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence).


title: Vector boson scattering at the LHC. A study of the WW -> WW channels with the Warsaw cut

abstract: We study W boson scattering in the same- and opposite-sign channels under the assumption that no resonances are present in the collider processes $pp\rightarrow l^{\pm}\nu_{l}l^{\pm}\nu_{l}jj$ and $pp\rightarrow l^{\pm}\nu_{l}l^{\mp}\nu_{l}jj$, respectively. Basic selection cuts together with a restriction on the combination of the final lepton and jet momenta (the Warsaw cut) makes it possible to argue that at the LHC a luminosity of 100 fb$^{-1}$ and a center-of-mass energy of $\sqrt{s}= 13$ TeV will allow to constrain the leading effective lagrangian coefficients at the permil level. We also discuss limits on the other coefficients of the effective lagrangian as well as stronger constraints provided by higher energy and luminosity. We show that the same-sign $WW \rightarrow WW$ channel suffices in providing the most stringent constraints.


title: Early Dust Formation and a Massive Progenitor for SN 2011ja?

abstract: SN 2011ja was a bright (I = -18.3) Type II supernova occurring in the nearby edge on spiral galaxy NGC 4945. Flat-topped and multi-peaked H-alpha and H-beta spectral emission lines appear between 64 - 84 days post-explosion, indicating interaction with a disc-like circumstellar medium inclined 30-45 degrees from edge-on. After day 84 an increase in the H- and K-band flux along with heavy attenuation of the red wing of the emission lines are strong indications of early dust formation, likely located in the cool dense shell created between the forward shock of the SN ejecta and the reverse shock created as the ejecta plows into the existing CSM. Radiative transfer modeling reveals both ~1.5 x 10^-4 Msun of pre-existing dust located ~ 10^16.7 cm away and ~ 5 x 10^-5 Msun of newly formed dust. Spectral observations after 1.5 years reveal the possibility that the fading SN is located within a young (3-6 Myr) massive stellar cluster, which when combined with tentative 56Ni mass estimates of 0.2 Msun may indicate a massive (> 25 Msun) progenitor for SN 2011ja.


title: Exchange-interaction of two spin qubits mediated by a superconductor

abstract: Entangling two quantum bits by letting them interact is the crucial requirements for building a quantum processor. For qubits based on the spin of the electron, these two-qubit gates are typically performed by exchange interaction of the electrons captured in two nearby quantum dots. Since the exchange interaction relies on tunneling of the electrons, the range of interaction for conventional approaches is severely limited as the tunneling amplitude decays exponentially with the length of the tunneling barrier. Here, we present a novel approach to couple two spin qubits via a superconducting coupler. In essence, the superconducting coupler provides a tunneling barrier for the electrons which can be tuned with exquisite precision. We show that as a result exchange couplings over a distance of several microns become realistic, thus enabling flexible designs of multi-qubit systems.


title: Probing the Physical Properties of z=4.5 Lyman Alpha Emitters with Spitzer

abstract: We present the results from a stellar population modeling analysis of a sample of 162 z=4.5, and 14 z=5.7 Lyman alpha emitting galaxies (LAEs) in the Bootes field, using deep Spitzer/IRAC data at 3.6 and 4.5 um from the Spitzer Lyman Alpha Survey, along with Hubble Space Telescope NICMOS and WFC3 imaging at 1.1 and 1.6 um for a subset of the LAEs. This represents one of the largest samples of high-redshift LAEs imaged with Spitzer IRAC. We find that 30/162 (19%) of the z=4.5 LAEs and 9/14 (64%) of the z=5.7 LAEs are detected at >3-sigma in at least one IRAC band. Individual z=4.5 IRAC-detected LAEs have a large range of stellar mass, from 5x10^8 to 10^11 Msol. One-third of the IRAC-detected LAEs have older stellar population ages of 100 Myr - 1 Gyr, while the remainder have ages < 100 Myr. A stacking analysis of IRAC-undetected LAEs shows this population to be primarily low mass (8 -- 20 x 10^8 Msol) and young (64 - 570 Myr). We find a correlation between stellar mass and the dust-corrected ultraviolet-based star-formation rate (SFR) similar to that at lower redshifts, in that higher mass galaxies exhibit higher SFRs. However, the z=4.5 LAE correlation is elevated 4-5 times in SFR compared to continuum-selected galaxies at similar redshifts. The exception is the most massive LAEs which have SFRs similar to galaxies at lower redshifts suggesting that they may represent a different population of galaxies than the traditional lower-mass LAEs, perhaps with a different mechanism promoting Lyman alpha photon escape.


title: Physical processes in protoplanetary disks

abstract: This review introduces physical processes in protoplanetary disks relevant to accretion and the initial stages of planet formation. After a brief overview of the observational context, I introduce the elementary theory of disk structure and evolution, review the gas-phase physics of angular momentum transport through turbulence and disk winds, and discuss possible origins for the episodic accretion observed in Young Stellar Objects. Turning to solids, I review the evolution of single particles under aerodynamic forces, and describe the conditions necessary for the development of collective gas-particle instabilities. Observations show that disks can exhibit pronounced large-scale structure, and I discuss the types of structures that may form from gas and particle interactions at ice lines, vortices and zonal flows, prior to the formation of large planetary bodies. I conclude with disk dispersal.


title: Completeness on the worm domain and the M\"untz-Sz\'asz problem for the Bergman space

abstract: In this paper we are concerned with the problem of completeness in the Bergman space of the worm domain $\mathcal{W}_\mu$ and its truncated version $\mathcal{W}'_\mu$. We determine some orthogonal systems and show that they are not complete, while showing that the union of two particular of such systems is complete. In order to prove our completeness result we introduce the Muentz-Szasz problem for the 1-dimensional Bergman space of the disk $\{\zeta : |\zeta-1|<1\}$ and find a sufficient condition for its solution.


title: The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Effect of smoothing of density field on reconstruction and anisotropic BAO analysis

abstract: The reconstruction algorithm introduced by \cite{Eis07}, which is widely used in clustering analysis, is based on the inference of the first order Lagrangian displacement field from the Gaussian smoothed galaxy density field in redshift space. The \modif2{smoothing scale} applied to the density field affects the inferred displacement field that is used to move {the galaxies}, and partially \modif2{erases} the nonlinear evolution {of the density field}. In this article, we explore this crucial step \modif2{in} the reconstruction algorithm. We study the performance of the reconstruction technique using two metrics: first, we study the performance using the anisotropic clustering, extending previous studies focused on isotropic clustering; second, we study its effect on the displacement field. We find that smoothing has a strong effect in the quadrupole of the correlation function and affects the accuracy and precision \modif2{with} which we can measure $D_A (z)$ and $H(z)$. We find that the optimal smoothing scale to use in the reconstruction algorithm applied to BOSS-CMASS is between 5-10 $h^{-1}$Mpc. Varying from the "usual" 15$h^{-1}$Mpc to $5 h^{-1}$Mpc \modif2{shows} $\sim$ 0.3\% variations in $D_A(z)$ and $\sim$ 0.4\% $H(z)$ and uncertainties are also reduced by 40\% and 30\% respectively. We also find that the accuracy of velocity field reconstruction depends strongly on the smoothing scale used for the density field. We measure the bias and uncertainties associated with different choices of smoothing length.


title: The JCMT Gould Belt Survey: A Quantitative Comparison Between SCUBA-2 Data Reduction Methods

abstract: Performing ground-based submillimetre observations is a difficult task as the measurements are subject to absorption and emission from water vapour in the Earth's atmosphere and time variation in weather and instrument stability. Removing these features and other artifacts from the data is a vital process which affects the characteristics of the recovered astronomical structure we seek to study. In this paper, we explore two data reduction methods for data taken with the Submillimetre Common-User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). The JCMT Legacy Reduction 1 (JCMT LR1) and The Gould Belt Legacy Survey Legacy Release 1 (GBS LR1) reduction both use the same software, Starlink, but differ in their choice of data reduction parameters. We find that the JCMT LR1 reduction is suitable for determining whether or not compact emission is present in a given region and the GBS LR1 reduction is tuned in a robust way to uncover more extended emission, which better serves more in-depth physical analyses of star-forming regions. Using the GBS LR1 method, we find that compact sources are recovered well, even at a peak brightness of only 3 times the noise, whereas the reconstruction of larger objects requires much care when drawing boundaries around the expected astronomical signal in the data reduction process. Incorrect boundaries can lead to false structure identification or it can cause structure to be missed. In the JCMT LR1 reduction, the extent of the true structure of objects larger than a point source is never fully recovered.


title: The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies

abstract: We measure and analyse the clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) relative to the line-of-sight (LOS), for LOWZ and CMASS galaxy samples drawn from the final Data Release 12 (DR12). The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm lowz}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm cmass}=0.57$. From the power spectrum monopole and quadrupole moments around the LOS, we measure the growth of structure parameter $f$ times the amplitude of dark matter density fluctuations $\sigma_8$ by modeling the Redshift-Space Distortion signal. When the geometrical Alcock-Paczynski effect is also constrained from the same data, we find joint constraints on $f\sigma_8$, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon $D_A(z)/r_s(z_d)$. We find $f(z_{\rm lowz})\sigma_8(z_{\rm lowz})=0.394\pm0.062$, $D_A(z_{\rm lowz})/r_s(z_d)=6.35\pm0.19$, $H(z_{\rm lowz})r_s(z_d)=(11.41\pm 0.56)\,{10^3\rm km}s^{-1}$ for the LOWZ sample, and $f(z_{\rm cmass})\sigma_8(z_{\rm cmass})=0.444\pm0.038$, $D_A(z_{\rm cmass})/r_s(z_d)=9.42\pm0.15$, $H(z_{\rm cmass})r_s(z_d)=(13.92 \pm 0.44)\, {10^3\rm km}s^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 measurements. Assuming the Hubble parameter and angular distance parameter are fixed at fiducial $\Lambda$CDM values, we find $f(z_{\rm lowz})\sigma_8(z_{\rm lowz})=0.485\pm0.044$ and $f(z_{\rm cmass})\sigma_8(z_{\rm cmass})=0.436\pm0.022$ for the LOWZ and CMASS samples, respectively.


title: Smooth gauge and Wannier functions for topological band structures in arbitrary dimensions

abstract: The construction of exponentially localized Wannier functions for a set of bands requires a choice of Bloch-like functions that span the space of the bands in question, and are smooth and periodic functions of k in the entire Brillouin zone. For bands with nontrivial topology, such smooth Bloch functions can only be chosen such that they do not respect the symmetries that protect the topology. This symmetry breaking is a necessary, but not sufficient condition for smoothness, and, in general, finding smooth Bloch functions for topological bands is a complicated task. We present a generic technique for finding smooth Bloch functions and constructing exponentially localized Wannier functions in the presence of nontrivial topology, given that the net Chern number of the bands in question vanishes. The technique is verified against known results in the Kane-Mele model. It is then applied to the topological insulator Bi2Se3, where the topological state is protected by two symmetries: time reversal and inversion. The resultant exponentially localized Wannier functions break both these symmetries. Finally, we illustrate how the calculation of the Chern-Simons orbital magnetoelectric response is facilitated by the proposed smooth gauge construction.


title: Effects of Pb doping on structural and electronics properties of Bi$_2$Sr$_2$Ca$_2$Cu$_3$O$_{10}$

abstract: Pb doping effect in the Bi$_2$Sr$_2$Ca$_2$Cu$_3$O$_{10}$ compound (Bi2223) on the structural and electronic properties were investigated, using the Local Density (LDA) and Virtual Crystal (VCA) approximations within the framework of the Density Functional Theory (DFT), taking as reference the procedure implemented by H.Lin {\it et al.} in the Bi2212 compound [{\it Phys. Rev. Lett.} {\bf 96} (2006) 097001]. Results show that, the incorporation of Pb-dopant in Bi2223 lead a rigid displacement of the Bi/Pb-O bands toward higher energies, with a null contribution at the Fermi level, around the high symmetry point $\overline{\text{M}}$ in the irreducible Brillouin zone, for Pb doping concentration equal to or more than 26\%, avoiding the presence of the so-called Bi-O {\it pockets} in the Fermi surface, in good agreement with angle-resolved photoemission spectroscopy (ARPES) and nuclear magnetic resonance (NMR) experiments, although a slight metallic character of the Bi-O bonds is still observed which would disagree with some experimental reports. The calculations show that the changes on the structural properties are associated to the presence or absence of the Bi-O {\it pockets} in the Fermi surface


title: Approximation of small-amplitude weakly coupled oscillators with discrete nonlinear Schrodinger equations

abstract: Small-amplitude weakly coupled oscillators of the Klein-Gordon lattices are approximated by equations of the discrete nonlinear Schrodinger type. We show how to justify this approximation by two methods, which have been very popular in the recent literature. The first method relies on a priori energy estimates and multi-scale decompositions. The second method is based on a resonant normal form theorem. We show that although the two methods are different in the implementation, they produce equivalent results as the end product. We also discuss applications of the discrete nonlinear Schrodinger equation in the context of existence and stability of breathers of the Klein--Gordon lattice.


title: Exchange-Repairs: Managing Inconsistency in Data Exchange

abstract: In a data exchange setting with target constraints, it is often the case that a given source instance has no solutions. In such cases, the semantics of target queries trivialize. The aim of this paper is to introduce and explore a new framework that gives meaningful semantics in such cases by using the notion of exchange-repairs. Informally, an exchange-repair of a source instance is another source instance that differs minimally from the first, but has a solution. Exchange-repairs give rise to a natural notion of exchange-repair certain answers (XR-certain answers) for target queries. We show that for schema mappings specified by source-to-target GAV dependencies and target equality-generating dependencies (egds), the XR-certain answers of a target conjunctive query can be rewritten as the consistent answers (in the sense of standard database repairs) of a union of conjunctive queries over the source schema with respect to a set of egds over the source schema, making it possible to use a consistent query-answering system to compute XR-certain answers in data exchange. We then examine the general case of schema mappings specified by source-to-target GLAV constraints, a weakly acyclic set of target tgds and a set of target egds. The main result asserts that, for such settings, the XR-certain answers of conjunctive queries can be rewritten as the certain answers of a union of conjunctive queries with respect to the stable models of a disjunctive logic program over a suitable expansion of the source schema.


title: Evidence That Hydra I is a Tidally Disrupting Milky Way Dwarf Galaxy

abstract: The Eastern Banded Structure (EBS) and Hydra~I halo overdensity are very nearby (d $\sim$ 10 kpc) objects discovered in SDSS data. Previous studies of the region have shown that EBS and Hydra I are spatially coincident, cold structures at the same distance, suggesting that Hydra I may be the EBS's progenitor. We combine new wide-field DECam imaging and MMT/Hectochelle spectroscopic observations of Hydra I with SDSS archival spectroscopic observations to quantify Hydra I's present-day chemodynamical properties, and to infer whether it originated as a star cluster or dwarf galaxy. While previous work using shallow SDSS imaging assumed a standard old, metal-poor stellar population, our deeper DECam imaging reveals that Hydra~I has a thin, well-defined main sequence turnoff of intermediate age ($\sim 5-6$ Gyr) and metallicity ([Fe/H] = $-0.9$ dex). We measure statistically significant spreads in both the iron and alpha-element abundances of $\sigma_{[Fe/H]} = 0.13 \pm 0.02$ dex and $\sigma_{[\alpha/{\rm Fe}]} = 0.09 \pm 0.03$ dex, respectively, and place upper limits on both the rotation and its proper motion. Hydra~I's intermediate age and [Fe/H] -- as well as its low [$\alpha$/Fe], apparent [Fe/H] spread, and present-day low luminosity -- suggest that its progenitor was a dwarf galaxy, which subsequently lost more than $99.99\%$ of its stellar mass.


title: The Transverse Momentum Dependent Fragmentation Function at NNLO

abstract: We calculate the unpolarized transverse momentum dependent fragmentation function (TMDFF) at next-to-next-to-leading order (NNLO), evaluating separately TMD soft factor and TMD collinear correlator. For the first time the cancellation of spurious rapidity divergences in a properly defined individual TMD beyond the first non-trivial order is shown. This represents a strong check of the given TMD definition. We extract the matching coefficient necessary to perform the transverse momentum resummation at next-to-next-to-next-to-leading-logarithmic accuracy. The universal character of the soft function, which enters the definition of all (un)polarized TMD distribution/fragmentation functions, facilitates the future calculation of all the other TMDs and their coefficients at NNLO, pushing forward the accuracy of theoretical predictions for the current and next generation of high energy colliders.


title: Decomposing highly edge-connected graphs into paths of any given length

abstract: In 2006, Bar\'at and Thomassen posed the following conjecture: for each tree $T$, there exists a natural number $k_T$ such that, if $G$ is a $k_T$-edge-connected graph and $|E(G)|$ is divisible by $|E(T)|$, then $G$ admits a decomposition into copies of $T$. This conjecture was verified for stars, some bistars, paths of length $3$, $5$, and $2^r$ for every positive integer $r$. We prove that this conjecture holds for paths of any fixed length.


title: Semidefinite programming relaxations for linear semi-infinite polynomial programming

abstract: This paper studies a class of so-called linear semi-infinite polynomial programming (LSIPP) problems. It is a subclass of linear semi-infinite programming problems whose constraint functions are polynomials in parameters and index sets are basic semialgebraic sets. We present a hierarchy of semidefinite programming (SDP) relaxations for LSIPP problems. Convergence rate analysis of the SDP relaxations is established based on some existing results. We show how to verify the compactness of feasible sets of LSIPP problems. In the end, we extend the SDP relaxation method to more general semi-infinite programming problems.


title: A hybrid recursive multilevel incomplete factorization preconditioner for solving general linear systems

abstract: In this paper we introduce an algebraic recursive multilevel incomplete factorization preconditioner, based on a distributed Schur complement formulation, for solving general linear systems. The novelty of the proposed method is to combine factorization techniques of both implicit and explicit type, recursive combinatorial algorithms, multilevel mechanisms and overlapping strategies to maximize sparsity in the inverse factors and consequently reduce the factorization costs. Numerical experiments demonstrate the good potential of the proposed solver to precondition effectively general linear systems, also against other state-of-the-art iterative solvers of both implicit and explicit form.


title: Wide-Band, Low-Frequency Pulse Profiles of 100 Radio Pulsars with LOFAR

abstract: LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

categories: astro-ph.HE astro-ph.GA astro-ph.IM astro-ph.SR


title: SnapVX: A Network-Based Convex Optimization Solver

abstract: SnapVX is a high-performance Python solver for convex optimization problems defined on networks. For these problems, it provides a fast and scalable solution with guaranteed global convergence. SnapVX combines the capabilities of two open source software packages: Snap.py and CVXPY. Snap.py is a large scale graph processing library, and CVXPY provides a general modeling framework for small-scale subproblems. SnapVX offers a customizable yet easy-to-use interface with out-of-the-box functionality. Based on the Alternating Direction Method of Multipliers (ADMM), it is able to efficiently store, analyze, and solve large optimization problems from a variety of different applications. Documentation, examples, and more can be found on the SnapVX website at http://snap.stanford.edu/snapvx.


title: Isometric Embeddability of Snowflakes

abstract: We show that a snowflake of a metric space with positive Hausdorff dimension does not admit an isometric embedding into euclidean space.


title: A Note on a Multilinear Local $Tb$ Theorem for Calder\'on-Zygmund Operators

abstract: In this short note, we extend a local $Tb$ theorem that was proved in \cite{GHO} to a full multilinear local $Tb$ theorem.


title: The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: mock galaxy catalogues for the BOSS Final Data Release

abstract: We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $\sigma$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.


title: Moving point charge as a soliton in nonlinear electrodynamics

abstract: The field of a moving pointlike charge is determined in nonlinear local electrodynamics. As a model Lagrangian for the latter we take the one whose nonlinearity is the Euler-Heisenberg Lagrangian of quantum electrodynamics truncated at the leading term of its expansion in powers of the first field invariant. The total energy of the field produced by a point charge is finite in that model; thereby its field configuration is a soliton. We define a finite energy-momentum vector of this field configuration to demonstrate that its components satisfy the standard mechanical relation characteristic of a free moving massive particle.


title: Creature forcing and topological Ramsey spaces

abstract: This article introduces a line of investigation into connections between creature forcings and topological Ramsey spaces. Three examples of sets of pure candidates for creature forcings are shown to contain dense subsets which are actually topological Ramsey spaces. A new variant of the product tree Ramsey theorem is proved in order to obtain the pigeonhole principles for two of these examples.


title: Consistency bands for the mean excess function and application to graphical goodness of fit test for financial data

abstract: In this paper, we use the modern setting of functional empirical processes and recent techniques on uniform estimation for non parametric objects to derive consistency bands for the mean excess function in the i.i.d. case. We apply our results for modelling financial data, in particular Dow Jones data basis to see how good the Generalized hyperbolic distribution models fit monthly data.


title: The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Modeling the clustering and halo occupation distribution of BOSS-CMASS galaxies in the Final Data Release

abstract: We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $\Lambda$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$\sigma$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$\sigma$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.


title: Submaximally symmetric CR-structures

abstract: Hypersurface type CR-structures with non-degenerate Levi form on a manifold of dimension $(2n+1)$ have maximal symmetry dimension $n^2+4n+3$. We prove that the next (submaximal) possible dimension for a (local) symmetry algebra is $n^2+4$ for Levi-indefinite structures and $n^2+3$ for Levi-definite structures when $n>1$. In the exceptional case of CR-dimension $n=1$, the submaximal symmetry dimension 3 was computed by E.\,Cartan.


title: Contraction of Hamiltonian $K$-spaces

abstract: In the spirit of recent work of Harada-Kaveh and Nishinou-Nohara-Ueda, we study the symplectic geometry of Popov's horospherical degenerations of complex algebraic varieties with the action of a complex linearly reductive group. We formulate an intrinsic symplectic contraction of a Hamiltonian space, which is a surjective, continuous map onto a new Hamiltonian space that is a symplectomorphism on an explicitly defined dense open subspace. This map is given by a precise formula, using techniques from the theory of symplectic reduction and symplectic implosion. We then show, using the Vinberg monoid, that the gradient-Hamiltonian flow for a horospherical degeneration of an algebraic variety gives rise to this contraction from a general fiber to the special fiber. We apply this construction to branching problems in representation theory, and finally we show how the Gel'fand-Tsetlin integrable system can be understood to arise this way.


title: What Powers the 3000-Day Light Curve of SN 2006gy?

abstract: SN 2006gy was the most luminous SN ever observed at the time of its discovery and the first of the newly defined class of superluminous supernovae (SLSNe). The extraordinary energetics of SN 2006gy and all SLSNe (>10^51 erg) require either atypically large explosion energies (e.g., pair-instability explosion) or the efficient conversion of kinetic into radiative energy (e.g., shock interaction). The mass-loss characteristics can therefore offer important clues regarding the progenitor system. For the case of SN 2006gy, both a scattered and thermal light echo from circumstellar material (CSM) have been reported at later epochs (day ~800), ruling out the likelihood of a pair-instability event and leading to constraints on the characteristics of the CSM. Owing to the proximity of the SN to the bright host-galaxy nucleus, continued monitoring of the light echo has not been trivial, requiring the high resolution offered by the Hubble Space Telescope (HST) or ground-based adaptive optics (AO). Here we report detections of SN 2006gy using HST and Keck AO at ~3000 days post-explosion and consider the emission mechanism for the very late-time light curve. While the optical light curve and optical spectral energy distribution are consistent with a continued scattered-light echo, a thermal echo is insufficient to power the K'-band emission by day 3000. Instead, we present evidence for late-time infrared emission from dust that is radiatively heated by CSM interaction within an extremely dense dust shell, and we consider the implications on the CSM characteristics and progenitor system.


title: Sections of the regular simplex - Volume formulas and estimates

abstract: We state a general formula to compute the volume of the intersection of the regular $n$-simplex with some $k$-dimensional subspace. It is known that for central hyperplanes the one through the centroid containing $n-1$ vertices gives the maximal volume. We show that, for fixed small distances of a hyperplane to the centroid, the hyperplane containing $n-1$ vertices is still volume maximizing. The proof also yields a new and short argument for the result on central sections. With the same technique we give a partial result for the minimal central hyperplane section. Finally, we obtain a bound for $k$-dimensional sections.


title: Isotropic and Anisotropic Regimes of the Field-Dependent Spin Dynamics in Sr2IrO4: Raman Scattering Studies

abstract: A major focus of experimental interest in Sr2IrO4 has been to clarify how the magnetic excitations of this strongly spin-orbit coupled system differ from the predictions of anisotropic 2D spin-1/2 Heisenberg model and to explore the extent to which strong spin-orbit coupling affects the magnetic properties of iridates. Here, we present a high-resolution inelastic light (Raman) scattering study of the low energy magnetic excitation spectrum of Sr2IrO4 and doped Eu-doped Sr2IrO4 as functions of both temperature and applied magnetic field. We show that the high-field (H>1.5 T) in-plane spin dynamics of Sr2IrO4 are isotropic and governed by the interplay between the applied field and the small in-plane ferromagnetic spin components induced by the Dzyaloshinskii-Moriya interaction. However, the spin dynamics of Sr2IrO4 at lower fields (H<1.5 T) exhibit important effects associated with interlayer coupling and in-plane anisotropy, including a spin-flop transition at Hc in Sr2IrO4 that occurs either discontinuously or via a continuous rotation of the spins, depending upon the in-plane orientation of the applied field. These results show that in-plane anisotropy and interlayer coupling effects play important roles in the low-field magnetic and dynamical properties of Sr2IrO4.


title: Sharper periodicity and stabilization maps for configuration spaces of closed manifolds

abstract: In this note we study the homology of configuration spaces of closed manifolds. We sharpen the eventual periodicity results of Nagpal and Cantero-Palmer, prove integral homological stability for configuration spaces of odd-dimensional manifolds and introduce a stabilization map on the homology with $Z[1/2]$-coefficients of configuration spaces of odd-dimensional manifolds.


title: From Quantum Chaos and Eigenstate Thermalization to Statistical Mechanics and Thermodynamics

abstract: This review gives a pedagogical introduction to the eigenstate thermalization hypothesis (ETH), its basis, and its implications to statistical mechanics and thermodynamics. In the first part, ETH is introduced as a natural extension of ideas from quantum chaos and random matrix theory (RMT). To this end, we present a brief overview of classical and quantum chaos, as well as RMT and some of its most important predictions. The latter include the statistics of energy levels, eigenstate components, and matrix elements of observables. Building on these, we introduce the ETH and show that it allows one to describe thermalization in isolated chaotic systems without invoking the notion of an external bath. We examine numerical evidence of eigenstate thermalization from studies of many-body lattice systems. We also introduce the concept of a quench as a means of taking isolated systems out of equilibrium, and discuss results of numerical experiments on quantum quenches. The second part of the review explores the implications of quantum chaos and ETH to thermodynamics. Basic thermodynamic relations are derived, including the second law of thermodynamics, the fundamental thermodynamic relation, fluctuation theorems, and the Einstein and Onsager relations. In particular, it is shown that quantum chaos allows one to prove these relations for individual Hamiltonian eigenstates and thus extend them to arbitrary stationary statistical ensembles. We then show how one can use these relations to obtain nontrivial universal energy distributions in continuously driven systems. At the end of the review, we briefly discuss the relaxation dynamics and description after relaxation of integrable quantum systems, for which ETH is violated. We introduce the concept of the generalized Gibbs ensemble, and discuss its connection with ideas of prethermalization in weakly interacting systems.

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


title: Generating maps on surfaces

abstract: We describe procedures for generating all 2-cell embedded simple graphs with up to a fixed number of vertices on a given surface. We also modify these procedures to generate closed 2-cell embeddings and polyhedral embeddings. We give results of computer implementations of these procedures for seven surfaces: the sphere, the torus, the double torus, the projective plane, the Klein bottle, the triple cross surface, and the quadruple cross surface.


title: Constraining cosmology with pairwise velocity estimator

abstract: In this paper, we develop a full statistical method for the pairwise velocity estimator previously proposed, and apply Cosmicflows-2 catalogue to this method to constrain cosmology. We first calculate the covariance matrix for line-of-sight velocities for a given catalogue, and then simulate the mock full-sky surveys from it, and then calculate the variance for the pairwise velocity field. By applying the $8315$ independent galaxy samples and compressed $5224$ group samples from Cosmicflows-2 catalogue to this statistical method, we find that the joint constraint on $\Omega^{0.6}_{\rm m}h$ and $\sigma_{8}$ is completely consistent with the WMAP 9-year and Planck 2015 best-fitting cosmology. Currently, there is no evidence for the modified gravity models or any dynamic dark energy models from this practice, and the error-bars need to be reduced in order to provide any concrete evidence against/to support $\Lambda$CDM cosmology.


title: Stability and periodicity in the modular representation theory of symmetric groups

abstract: We study asymptotic properties of the modular representation theory of symmetric groups and investigate modular analogs of stabilization phenomena in characteristic zero. The main results are equivalences of categories between certain abelian subcategories of representations of $S_n$ and $S_m$ for different $n$ and $m$. We apply these results to obtain a structural result for $FI$-modules, and to prove a result conjectured by Deligne in a recent letter to Ostrik.


title: The AU Mic Debris Disk: far-infrared and submillimeter resolved imaging

abstract: We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500 and 850 micron. The disk is resolved at 70, 160 and 450 micron. In addition to the planetesimal belt, we detect thermal emission from AU Mic's halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is $3.9 \times 10^{-4}$ and its mm-grain dust mass is 0.01 MEarth (+/- 20%). We create a simple spatial model that reconciles the disk SED as a blackbody of 53 +/- 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best fit model is consistent with the "birth ring" model explored in earlier works, i.e., an edge-on dust belt extending from 8.8-40 AU, but with an additional halo component with an $r^{-1.5}$ surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass loss rates of 10-100x solar, compact (zero porosity) grains can only be removed if they are very small, consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 micron can be removed via corpuscular forces (i.e., the stellar wind).


title: Cosmic-ray acceleration in young protostars

abstract: The main signature of the interaction between cosmic rays and molecular clouds is the high ionisation degree. This decreases towards the densest parts of a cloud, where star formation is expected, because of energy losses and magnetic effects. However recent observations hint to high levels of ionisation in protostellar systems, therefore leading to an apparent contradiction that could be explained by the presence of energetic particles accelerated within young protostars. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient particle acceleration through the diffusive shock acceleration mechanism. We find that jet shocks can be strong accelerators of protons which can be boosted up to relativistic energies. Another possibly efficient acceleration site is located at protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate protons. Our results demonstrate the possibility of accelerating particles during the early phase of a proto-Solar-like system and can be used as an argument to support available observations. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the star and planet formation process as well as on the formation of pre-biotic molecules.


title: The Scale-Dependence of Halo Assembly Bias

abstract: The two-point clustering of dark matter halos is influenced by halo properties besides mass, a phenomenon referred to as halo assembly bias. Using the depth of the gravitational potential well, $V_{\rm max}$, as our secondary halo property, in this paper we present the first study of the scale-dependence assembly bias. In the large-scale linear regime, $r\geq10h^{-1}{\rm Mpc},$ our findings are in keeping with previous results. In particular, at the low-mass end ($M_{\rm vir}<M_{\rm coll}\approx10^{12.5}{\rm M}_{\odot}$), halos with high-$V_{\rm max}$ show stronger large-scale clustering relative to halos with low-$V_{\rm max}$ of the same mass, this trend weakens and reverses for $M_{\rm vir}\geq M_{\rm coll}.$ In the nonlinear regime, assembly bias in low-mass halos exhibits a pronounced scale-dependent "bump" at $500h^{-1}{\rm kpc}-5h^{-1}{\rm Mpc},$ a new result. This feature weakens and eventually vanishes for halos of higher mass. We show that this scale-dependent signature can primarily be attributed to a special subpopulation of ejected halos, defined as present-day host halos that were previously members of a higher-mass halo at some point in their past history. A corollary of our results is that galaxy clustering on scales of $r\sim1-2h^{-1}{\rm Mpc}$ can be impacted by up to $\sim15\%$ by the choice of the halo property used in the halo model, even for stellar mass-limited samples.


title: Information-theoretic bounds for exact recovery in weighted stochastic block models using the Renyi divergence

abstract: We derive sharp thresholds for exact recovery of communities in a weighted stochastic block model, where observations are collected in the form of a weighted adjacency matrix, and the weight of each edge is generated independently from a distribution determined by the community membership of its endpoints. Our main result, characterizing the precise boundary between success and failure of maximum likelihood estimation when edge weights are drawn from discrete distributions, involves the Renyi divergence of order $\frac{1}{2}$ between the distributions of within-community and between-community edges. When the Renyi divergence is above a certain threshold, meaning the edge distributions are sufficiently separated, maximum likelihood succeeds with probability tending to 1; when the Renyi divergence is below the threshold, maximum likelihood fails with probability bounded away from 0. In the language of graphical channels, the Renyi divergence pinpoints the information-theoretic capacity of discrete graphical channels with binary inputs. Our results generalize previously established thresholds derived specifically for unweighted block models, and support an important natural intuition relating the intrinsic hardness of community estimation to the problem of edge classification. Along the way, we establish a general relationship between the Renyi divergence and the probability of success of the maximum likelihood estimator for arbitrary edge weight distributions. Finally, we discuss consequences of our bounds for the related problems of censored block models and submatrix localization, which may be seen as special cases of the framework developed in our paper.


title: Florigens and antiflorigens: a molecular genetic understanding

abstract: Florigens, the leaf-derived signals that initiate flowering, have been described as mysterious, elusive and the Holy Grail of plant biology.


title: A Multi-Agent System Approach to Load-Balancing and Resource Allocation for Distributed Computing

abstract: In this research we use a decentralized computing approach to allocate and schedule tasks on a massively distributed grid. Using emergent properties of multi-agent systems, the algorithm dynamically creates and dissociates clusters to serve the changing resource demands of a global task queue. The algorithm is compared to a standard First-in First-out (FIFO) scheduling algorithm. Experiments done on a simulator show that the distributed resource allocation protocol (dRAP) algorithm outperforms the FIFO scheduling algorithm on time to empty queue, average waiting time and CPU utilization. Such a decentralized computing approach holds promise for massively distributed processing scenarios like SETI@home and Google MapReduce.


title: Another dual of MacMahon's theorem on plane partitions

abstract: In this paper we introduce a counterpart structure to the shamrocks studied in the paper "A dual of Macmahon's theorem on plane partitions" by M. Ciucu and C. Krattenthaler (Proc. Natl. Acad. Sci. USA, vol. 110 (2013), 4518-4523), which, just like the latter, can be included at the center of a lattice hexagon on the triangular lattice so that the region obtained from the hexagon by removing it has its number of lozenge tilings given by a simple product formula. The new structure, called a fern, consists of an arbitrary number of equilateral triangles of alternating orientations lined up along a lattice line. The shamrock and the fern seem to be the only structures with this property. It would be interesting to understand why these are the only two such structures.


title: Quasi-MLE for quadratic ARCH model with long memory

abstract: We discuss parametric quasi-maximum likelihood estimation for quadratic ARCH process with long memory introduced in Doukhan et al. (2015) and Grublyt\.e and \v{S}karnulis (2015) with conditional variance given by a strictly positive quadratic form of observable stationary sequence. We prove consistency and asymptotic normality of the corresponding QMLE estimates, including the estimate of long memory parameter $0< d < 1/2$. A simulation study of empirical MSE is included.


title: Leveraging waveform complexity for confident detection of gravitational waves

abstract: The recent completion of Advanced LIGO suggests that gravitational waves (GWs) may soon be directly observed. Past searches for gravitational-wave transients have been impacted by transient noise artifacts, known as glitches, introduced into LIGO data due to instrumental and environmental effects. In this work, we explore how waveform complexity, instead of signal-to-noise ratio, can be used to rank event candidates and distinguish short duration astrophysical signals from glitches. We test this framework using a new hierarchical pipeline that directly compares the Bayesian evidence of explicit signal and glitch models. The hierarchical pipeline is shown to have strong performance, and in particular, allows high-confidence detections of a range of waveforms at realistic signal-to-noise ratio with a two detector network.


title: Twisted Simplicial Groups and Twisted Homology of Categories

abstract: Let $A$ be either a simplicial complex $K$ or a small category $\mathcal C$ with $V(A)$ as its set of vertices or objects. We define a twisted structure on $A$ with coefficients in a simplicial group $G$ as a function $$ \delta\colon V(A)\longrightarrow \operatorname{End}(G), \quad v\mapsto \delta_v $$ such that $\delta_v\circ \delta_w=\delta_w\circ \delta_v$ if there exists an edge in $A$ joining $v$ with $w$ or an arrow either from $v$ to $w$ or from $w$ to $v$. We give a canonical construction of twisted simplicial group as well as twisted homology for $A$ with a given twisted structure. Also we determine the homotopy type of of this simplicial group as the loop space over certain twisted smash product.


title: Triviality properties of principal bundles on singular curves

abstract: We show that principal bundles for a semisimple group on an arbitrary affine curve over an algebraically closed field are trivial, provided the order of $\pi_1$ of the group is invertible in the ground field, or if the curve has semi-normal singularities. Several consequences and extensions of this result (and method) are given. As an application, we realize conformal blocks bundles on moduli stacks of stable curves as push forwards of line bundles on (relative) moduli stacks of principal bundles on the universal curve.


title: An integral formulation for wave propagation on weakly non-uniform potential flows

abstract: An integral formulation for acoustic radiation in moving flows is presented. It is based on a potential formulation for acoustic radiation on weakly non-uniform subsonic mean flows. This work is motivated by the absence of suitable kernels for wave propagation on non-uniform flow. The integral solution is formulated using a Green's function obtained by combining the Taylor and Lorentz transformations. Although most conventional approaches based on either transform solve the Helmholtz problem in a transformed domain, the current Green's function and associated integral equation are derived in the physical space. A dimensional error analysis is developed to identify the limitations of the current formulation. Numerical applications are performed to assess the accuracy of the integral solution. It is tested as a means of extrapolating a numerical solution available on the outer boundary of a domain to the far field, and as a means of solving scattering problems by rigid surfaces in non-uniform flows. The results show that the error associated with the physical model deteriorates with increasing frequency and mean flow Mach number. However, the error is generated only in the domain where mean flow non-uniformities are significant and is constant in regions where the flow is uniform.


title: Proximity induced quantum anomalous Hall effect in graphene/EuO hetero-structures

abstract: In a heterostructure of graphene and the ferromagnetic insulator EuO, the Eu atoms induce proximity exchange and inter-valley interactions in the graphene layer. Constrained by the lattice symmetries, and guided by ab initio calculations, a model Hamiltonian is constructed that describes the low-energy bands. Band parameters such as proximity induced exchange splitting, spin orbit coupling, and inter-valley interaction are determined. Calculations of the Chern number identify the conditions under which the hetero-structures exhibit topologically non-trivial bands that give rise to the quantum anomalous Hall effect with a Hall conductivity of $\sigma_{xy} = 2 e^2/h$.


title: Large-Scale Mode Identification and Data-Driven Sciences

abstract: Bump-hunting or mode identification is a fundamental problem that arises in almost every scientific field of data-driven discovery. Surprisingly, very few data modeling tools are available for automatic (not requiring manual case-by-base investigation), objective (not subjective), and nonparametric (not based on restrictive parametric model assumptions) mode discovery, which can scale to large data sets. This article introduces LPMode--an algorithm based on a new theory for detecting multimodality of a probability density. We apply LPMode to answer important research questions arising in various fields from environmental science, ecology, econometrics, analytical chemistry to astronomy and cancer genomics.


title: On Computational Paths and the Fundamental Groupoid of a Type

abstract: The main objective of this work is to study mathematical properties of computational paths. Originally proposed by de Queiroz \& Gabbay (1994) as `sequences of rewrites', computational paths can be seen as the grounds on which the propositional equality between two computational objects stand. Using computational paths and categorical semantics, we take any type $A$ of type theory and construct a groupoid for this type. We call this groupoid the fundamental groupoid of a type $A$, since it is similar to the one obtained using the homotopical interpretation of the identity type. The main difference is that instead of being just a semantical interpretation, computational paths are entities of the syntax of type theory. We also expand our results, using computational paths to construct fundamental groupoids of higher levels.


title: Parallel algorithms and concentration bounds for the Lovasz Local Lemma via witness DAGs

abstract: The Lov\'{a}sz Local Lemma (LLL) is a cornerstone principle in the probabilistic method of combinatorics, and a seminal algorithm of Moser & Tardos (2010) provides an efficient randomized algorithm to implement it. This can be parallelized to give an algorithm that uses polynomially many processors and runs in $O(\log^3 n)$ time on an EREW PRAM, stemming from $O(\log n)$ adaptive computations of a maximal independent set (MIS). Chung et al. (2014) developed faster local and parallel algorithms, potentially running in time $O(\log^2 n)$, but these algorithms require more stringent conditions than the LLL. We give a new parallel algorithm that works under essentially the same conditions as the original algorithm of Moser & Tardos but uses only a single MIS computation, thus running in $O(\log^2 n)$ time on an EREW PRAM. This can be derandomized to give an NC algorithm running in time $O(\log^2 n)$ as well, speeding up a previous NC LLL algorithm of Chandrasekaran et al. (2013). We also provide improved and tighter bounds on the run-times of the sequential and parallel resampling-based algorithms originally developed by Moser & Tardos. These apply to any problem instance in which the tighter Shearer LLL criterion is satisfied.


title: Limits on the anomalous $Wtq$ couplings

abstract: Within the model-independent framework of $SU(3)\times SU(2)\times U(1)$ gauge-invariant dimension-six operators, we study flavor-changing $Wtq$ couplings ($q=d,s$) and related four-quark contact interactions involving the top. We obtain bounds on those couplings from Tevatron and LHC data for single-top production and branching fractions in top decays, as well as other experimental results on flavor-changing neutral-current processes including $B\to X_{q}\gamma$ and $Z\to b\bar{q}$ decays ($q=d,s$). We also update the bounds on flavor-diagonal $Wtb$ couplings using the most recent measurement of $W$-helicity fractions in top decays from top-pair production.


title: Magnetic Catalysis in Graphene

abstract: One of the most important developments in condensed matter physics in recent years has been the discovery and characterization of graphene. A two-dimensional layer of Carbon arranged in a hexagonal lattice, graphene exhibits many interesting electronic properties, most notably that the low energy excitations behave as massless Dirac fermions. These excitations interact strongly via the Coulomb interaction and thus non-perturbative methods are necessary. Using methods borrowed from lattice QCD, we study the graphene effective theory in the presence of an external magnetic field. Graphene, along with other $(2+1)$-dimensional field theories, has been predicted to undergo spontaneous breaking of flavor symmetry including the formation of a gap as a result of the external magnetic field. This phenomenon is known as magnetic catalysis. Our study investigates magnetic catalysis using a fully non-perturbative approach.


title: Topological Susceptibility from Slabs

abstract: In quantum field theories with topological sectors, a non-perturbative quantity of interest is the topological susceptibility chi_t. In principle it seems straightforward to measure chi_t by means of Monte Carlo simulations. However, for local update algorithms and fine lattice spacings, this tends to be difficult, since the Monte Carlo history rarely changes the topological sector. Here we test a method to measure chi_t even if data from only one sector are available. It is based on the topological charges in sub-volumes, which we denote as slabs. Assuming a Gaussian distribution of these charges, this method enables the evaluation of chi_t, as we demonstrate with numerical results for non-linear sigma-models.


title: Efficient Reassembling of Graphs, Part 1: The Linear Case

abstract: The reassembling of a simple connected graph G = (V,E) is an abstraction of a problem arising in earlier studies of network analysis. Its simplest formulation is in two steps: (1) We cut every edge of G into two halves, thus obtaining a collection of n=|V| one-vertex components. (2) We splice the two halves of every edge together, not of all the edges at once, but in some ordering \Theta of the edges that minimizes two measures that depend on the edge-boundary degrees of assembled components. The edge-boundary degree of a component A (subset of V) is the number of edges in G with one endpoint in A and one endpoint in V-A. We call the maximum edge-boundary degree encountered during the reassembling process the alpha-measure of the reassembling, and the sum of all edge-boundary degrees is its beta-measure. The alpha-optimization (resp. beta-optimization) of the reassembling of G is to determine an order \Theta for splicing the edges that minimizes its alpha-measure (resp. beta-measure). There are different forms of reassembling. We consider only cases satisfying the condition that if the an edge between disjoint components A and B is spliced, then all the edges between A and B are spliced at the same time. In this report, we examine the particular case of linear reassembling, which requires that the next edge to be spliced must be adjacent to an already spliced edge. We delay other forms of reassembling to follow-up reports. We prove that alpha-optimization of linear reassembling and minimum-cutwidth linear arrangment (CutWidth) are polynomially reducible to each other, and that beta-optimization of linear reassembling and minimum-cost linear arrangement (MinArr) are polynomially reducible to each other.


title: Spectral analysis of stable processes on the positive half-line

abstract: We study the spectral expansion of the semigroup of a general stable process killed on the first exit from the positive half-line. Starting with the Wiener-Hopf factorization we obtain the q-resolvent density for the killed process, from which we derive the spectral expansion of the semigroup via the inverse Laplace transform. The eigenfunctions and co-eigenfunctions are given rather explicitly in terms of the double sine function and they give rise to a pair of integral transforms which generalize the classical Fourier sine transform. Our results provide the first explicit example of a spectral expansion of the semigroup of a non-symmetric Levy process killed on the first exit form the positive half-line.


title: Analysis of the ${\frac{1}{2}}^{\pm}$ pentaquark states in the diquark-diquark-antiquark model with QCD sum rules

abstract: In this article, we construct both the axialvector-diquark-axialvector-diquark-antiquark type and axialvector-diquark-scalar-diquark-antiquark type interpolating currents, then calculate the contributions of the vacuum condensates up to dimension-10 in the operator product expansion, and study the masses and pole residues of the $J^P={\frac{1}{2}}^\pm$ hidden-charm pentaquark states with the QCD sum rules in a systematic way. In calculations, we use the formula $\mu=\sqrt{M^2_{P}-(2{\mathbb{M}}_c)^2}$ to determine the energy scales of the QCD spectral densities. We take into account the $SU(3)$ breaking effects of the light quarks, and obtain the masses of the hidden charm pentaquark states with the strangeness $S=0,\,-1,\,-2,\,-3$, which can be confronted with the experimental data in the future.


title: Hyperbolic Dimension and Decomposition Complexity

abstract: The aim of this paper is to provide some new tools to aid the study of decomposition complexity, a notion introduced by Guentner, Tessera and Yu. In this paper, three equivalent definitions for decomposition complexity are established. We prove that metric spaces with finite hyperbolic dimension have finite (weak) decomposition complexity, and we prove that the collection of metric families that are coarsely embeddable into Hilbert space is closed under decomposition. A method for showing that certain metric spaces do not have finite decomposition complexity is also discussed.


title: Detecting nanohertz gravitational waves with pulsar timing arrays

abstract: Complementary to ground-based laser interferometers, pulsar timing array experiments are being carried out to search for nanohertz gravitational waves. Using the world's most powerful radio telescopes, three major international collaborations have collected $\sim$10-year high precision timing data for tens of millisecond pulsars. In this paper we give an overview on pulsar timing experiments, gravitational wave detection in the nanohertz regime, and recent results obtained by various timing array projects.


title: Exciton band structure in layered MoSe2: from a monolayer to the bulk limit

abstract: We present the micro-photoluminescence ($\mu$PL) and micro-reflectance contrast spectroscopy studies on thin films of MoSe2 with layer thicknesses ranging from a monolayer (1L) up to 5L. The thickness dependent evolution of the ground and excited state excitonic transitions taking place at various points of the Brillouin zone is determined. Temperature activated energy shifts and linewidth broadenings of the excitonic resonances in 1L, 2L and 3L flakes are accounted for by using standard formalisms previously developed for semiconductors. A peculiar shape of the optical response of the ground state (A) exciton in monolayer MoSe2 is tentatively attributed to the appearance of Fano-type resonance. Rather trivial and clearly decaying PL spectra of monolayer MoSe2 with temperature confirm that the ground state exciton in this material is optically bright in contrast to a dark exciton ground state in monolayer WSe2.


title: Measurement of jpsi production in p+p collisions at \sqrts = 500 GeV at STAR experiment

abstract: Quarkonium measurements in heavy-ion collisions play an essential role in understanding the hot, dense medium created in such collisions. As a reference, their production mechanism in p+p collisions needs to be thoroughly understood. In this paper, we report the measurement of inclusive cross section of jpsi with transverse momentum (pT) above 4 GeV/c at mid-rapidity in p+p collisions at \sqrts = 500 GeV by the STAR experiment. The ratio of the yield of psi(2S) to jpsi integrated over 4 < pT < 12 GeV/c is also presented. Furthermore, the jpsi yields are studied in different event multiplicity bins in different jpsi pT regions, where the low pT measurement is enabled by the newly installed Muon Telescope Detector. A strong increase of the relative jpsi yield with the event multiplicity is observed for all pT with significant pT dependence.


title: On Stickelberger Elements for $\mathbb{Q}(\zeta_{p^{n+1}})^+$ and $p$-adic $L$-functions

abstract: We give a survey of a couple known constructions of $p$-adic $L$-functions including Iwasawa's construction from classical Stickelberger elements. We then construct "real" Stickelberger elements, i.e., explicit elements in the Galois group ring with $\mathbb{Z}_p$ coefficients that annihilate the Sylow $p$-subgroup of the ideal class group of $\mathbb{Q}(\zeta_{p^{n+1}})^+$. In analogy with Iwasawa's work, we show that these elements are coherent in $\mathbb{Z}_p$-towers and give rise to twisted $p$-adic $L$-functions.


title: A new approach for analyzing panel AR(1) series with application to the unit root test

abstract: This paper derives several novel tests to improve on the t-test for testing AR(1) coefficients of panel time series, i.e., of multiple time series, when each has a small number of observations. These tests can determine the acceptance or the rejection of each hypothesis individually while controlling the average type one error. Strikingly, the testing statistics derived by the empirical Bayes approach can be approximated by a simple form similar to the t-statistic; the only difference is that the means and the variances are estimated by shrinkage estimators. Simulations demonstrate that the proposed tests have higher average power than the t-test in all settings we examine including those when the priors are miss-specified and the cross section series are dependent.


title: Cosmic Web Reconstruction through Density Ridges: Catalogue

abstract: We construct a catalogue for filaments using a novel approach called SCMS (subspace constrained mean shift; Ozertem & Erdogmus 2011; Chen et al. 2015). SCMS is a gradient-based method that detects filaments through density ridges (smooth curves tracing high-density regions). A great advantage of SCMS is its uncertainty measure, which allows an evaluation of the errors for the detected filaments. To detect filaments, we use data from the Sloan Digital Sky Survey, which consist of three galaxy samples: the NYU main galaxy sample (MGS), the LOWZ sample and the CMASS sample. Each of the three dataset covers different redshift regions so that the combined sample allows detection of filaments up to z = 0.7. Our filament catalogue consists of a sequence of two-dimensional filament maps at different redshifts that provide several useful statistics on the evolution cosmic web. To construct the maps, we select spectroscopically confirmed galaxies within 0.050 < z < 0.700 and partition them into 130 bins. For each bin, we ignore the redshift, treating the galaxy observations as a 2-D data and detect filaments using SCMS. The filament catalogue consists of 130 individual 2-D filament maps, and each map comprises points on the detected filaments that describe the filamentary structures at a particular redshift. We also apply our filament catalogue to investigate galaxy luminosity and its relation with distance to filament. Using a volume-limited sample, we find strong evidence (6.1$\sigma$ - 12.3$\sigma$) that galaxies close to filaments are generally brighter than those at significant distance from filaments.


title: Borell's generalized Pr\'ekopa-Leindler inequality: A simple proof

abstract: We present a simple proof of Christer Borell's general inequality in the Brunn-Minkowski theory. We then discuss applications of Borell's inequality to the log-Brunn-Minkowski inequality of B\"or\"oczky, Lutwak, Yang and Zhang.


title: A note on closedness of the sum of two closed subspaces in a Banach space

abstract: Let $X$ be a Banach space, and $M,N$ be two closed subspaces of $X$. We present several necessary and sufficient conditions for the closedness of $M+N$ ($M+N$ is not necessarily direct sum).


title: Toward superconducting critical current by design

abstract: We present the new paradigm of critical current by design. Analogous to materials by design, it aims at predicting the optimal defect landscape in a superconductor for targeted applications by elucidating the vortex dynamics responsible for the bulk critical current. To highlight this approach, we demonstrate the synergistic combination of critical current measurements on commercial high-temperature superconductors containing self-assembled and irradiation tailored correlated defects by using large-scale time-dependent Ginzburg-Landau simulations for vortex dynamics.


title: Complementary magnetic localized surface plasmons

abstract: Magnetic localized surface plasmons (LSPs) supported on metallic structures corrugated by very long and curved grooves have been recently proposed and demonstrated on an extremely thin metallic spiral structure (MSS) in the microwave regime. However, the mode profile for the magnetic LSPs was demonstrated by measuring only the electric field, not the magnetic field.


title: Steady state thermodynamics in population dynamics

abstract: We report that population dynamics in fluctuating environment accompanies mathematically equivalent structure to steady state thermodynamics. By employing the structure, population growth in fluctuating environment is decomposed into housekeeping and excess parts. The housekeeping part represents the integral of stationary growth rate for each condition during a history of the environmental change. The excess part accounts for the excess growth generated when environment is switched. Focusing on the excess growth, we obtain Clausius inequality, which gives the upper bound of the excess growth. The equality is shown to be achieved in quasistatic environmental changes. We also clarify that this bound can be evaluated by "lineage fitness" that is an experimentally observable quantity.


title: Efficient Neighborhood Selection for Gaussian Graphical Models

abstract: This paper addresses the problem of neighborhood selection for Gaussian graphical models. We present two heuristic algorithms: a forward-backward greedy algorithm for general Gaussian graphical models based on mutual information test, and a threshold-based algorithm for walk summable Gaussian graphical models. Both algorithms are shown to be structurally consistent, and efficient. Numerical results show that both algorithms work very well.


title: Flexibility of nucleic acids: from DNA to RNA

abstract: The structural flexibility of nucleic acids plays a key role in many fundamental life processes, such as gene replication and expression, DNA-protein recognition, and gene regulation. To obtain a thorough understanding of nucleic acid flexibility, extensive studies have been performed using various experimental methods and theoretical models. In this review, we will introduce the progress that has been made in understanding the flexibility of nucleic acids including DNAs and RNAs, and will emphasize the experimental findings and the effects of salt, temperature, and sequence. Finally, we will discuss the major unanswered questions in understanding the flexibility of nucleic acids.

categories: physics.bio-ph cond-mat.soft cond-mat.stat-mech physics.chem-ph


title: From Facial Parts Responses to Face Detection: A Deep Learning Approach

abstract: In this paper, we propose a novel deep convolutional network (DCN) that achieves outstanding performance on FDDB, PASCAL Face, and AFW. Specifically, our method achieves a high recall rate of 90.99% on the challenging FDDB benchmark, outperforming the state-of-the-art method by a large margin of 2.91%. Importantly, we consider finding faces from a new perspective through scoring facial parts responses by their spatial structure and arrangement. The scoring mechanism is carefully formulated considering challenging cases where faces are only partially visible. This consideration allows our network to detect faces under severe occlusion and unconstrained pose variation, which are the main difficulty and bottleneck of most existing face detection approaches. We show that despite the use of DCN, our network can achieve practical runtime speed.


title: Non-Hermitian neutrino oscillations in matter with PT symmetric Hamiltonians

abstract: We introduce and develop a novel approach to extend the ordinary two-flavor neutrino oscillation formalism in matter using a non-Hermitian PT symmetric effective Hamiltonian. The condition of PT symmetry is weaker and less mathematical than that of hermicity, but more physical, and such an extension of the formalism can give rise to sub-leading effects in neutrino flavor transitions similar to the effects by so-called non-standard neutrino interactions. We derive the necessary conditions for the spectrum of the effective Hamiltonian to be real as well as the mappings between the fundamental and effective parameters. We find that the real spectrum of the effective Hamiltonian will depend on all new fundamental parameters introduced in the non-Hermitian PT symmetric extension of the usual neutrino oscillation formalism and that either i) the spectrum is exact and the effective leptonic mixing must always be maximal or ii) the spectrum is approximate and all new fundamental parameters must be small.


title: Complex ordering in spin networks: Critical role of adaptation rate for dynamically evolving interactions

abstract: Many complex systems can be represented as networks of dynamical elements whose states evolve in response to interactions with neighboring elements, noise and external stimuli. The collective behavior of such systems can exhibit remarkable ordering phenomena such as chimera order corresponding to coexistence of ordered and disordered regions. Often, the interactions in such systems can also evolve over time responding to changes in the dynamical states of the elements. Link adaptation inspired by Hebbian learning, the dominant paradigm for neuronal plasticity, has been earlier shown to result in structural balance by removing any initial frustration in a system that arises through conflicting interactions. Here we show that the rate of the adaptive dynamics for the interactions is crucial in deciding the emergence of different ordering behavior (including chimera) and frustration in networks of Ising spins. In particular, we observe that small changes in the link adaptation rate about a critical value result in the system exhibiting radically different energy landscapes, viz., smooth landscape corresponding to balanced systems seen for fast learning, and rugged landscapes corresponding to frustrated systems seen for slow learning.


title: Stabilization of helical macromolecular phases by confined bending

abstract: By means of extensive replica-exchange simulations of generic coarse-grained models for helical polymers, we systematically investigate the structural transitions into all possible helical phases for flexible and semiflexible elastic polymers with self-interaction under the influence of torsion barriers. The competing interactions lead to a variety of conformational phases including disordered helical arrangements, single helices, and ordered, tertiary helix bundles. Most remarkably, we find that a bending restraint entails a clear separation and stabilization of the helical phases. This aids in understanding why semiflexible polymers such as double-stranded DNA tend to form pronounced helical structures and proteins often exhibit an abundance of helical structures, such as helix bundles, within their tertiary structure.


title: Two-gap superconducting properties of alkaline-earth intercalated $A_{x}(NH_{3})Fe_{2}Se_{2}$ (A = Ba or Sr)

abstract: Superconducting properties were studied on high quality superconductors $Ba_{x}(NH_{3})Fe_{2}Se_{2}$ ($T_{c}$ = 39 K) and $Sr_{x}(NH_{3})Fe_{2}Se_{2}$ ($T_{c}$ = 44 K) prepared by intercalating Ba/Sr atoms into tetragonal $\beta$-FeSe by liquid ammonia. The elongated c-axis and almost unchanged a-axis of $Ba_{x}(NH_{3})Fe_{2}Se_{2}$, comparing with $\beta$-FeSe, suggested an unchanged intra-$Fe_{2}Se_{2}$-layer structure and the $T_{c}$ enhancement is due to a 3D to 2D-like Fermi surface transformation. The superconducting coherent lengths $\xi$(0), Ginzburg-Landau parameters $\kappa$ and penetration depths $\lambda$(0) obtained from the extrapolated lower and upper critical fields $B_{c1}$(0) and $B_{c2}$(0) indicates that both compounds are typical type-II superconductors. The temperature dependence of 1/$\lambda^{2}$(T) of $Ba_{x}(NH_{3})Fe_{2}Se_{2}$ deduced from the low field magnetic susceptibility shows a two-gap s-wave behaviour with superconducting gaps of $\Delta_{1}$ = 6.47 meV and $\Delta_{2}$ = 1.06 meV.


title: Production of Dynamic Frozen Waves: Controlling shape, location (and speed) of diffraction-resistant beams

abstract: In recent times, we experimentally realized a quite efficient modeling of the shape of diffraction-resistant optical beams; thus generating for the first time the so-called Frozen Waves (FW), whose longitudinal intensity pattern can be arbitrarily chosen, within a prefixed space interval of the propagation axis. Such waves possess a host of potential applications: in medicine, biomedical optics, optical tweezers, atom guiding, remote sensing, tractor beams, optical communications or metrology, and other topics in photonic areas. In this work, we extend our theory of FWs -- which led to beams endowed with a static envelope -- through a dynamic modeling of the FWs, whose shape is now allowed to evolve in time in a predetermined way. And we experimentally create such dynamic FWs in Optics, via a computational holographic technique and a spatial light modulator. Experimental results are here presented for two cases of dynamic FWs, one of the zeroth and the other of higher order, the last one being the most interesting, consisting in a cylindrical surface of light whose geometry changes in space and time.


title: Identifying collusion groups using spectral clustering

abstract: In an illiquid stock, traders can collude and place orders on a predetermined price and quantity at a fixed schedule. This is usually done to manipulate the price of the stock or to create artificial liquidity in the stock, which may mislead genuine investors. Here, the problem is to identify such group of colluding traders. We modeled the problem instance as a graph, where each trader corresponds to a vertex of the graph and trade corresponds to edges of the graph. Further, we assign weights on edges depending on total volume, total number of trades, maximum change in the price and commonality between two vertices. Spectral clustering algorithms are used on the constructed graph to identify colluding group(s). We have compared our results with simulated data to show the effectiveness of spectral clustering to detecting colluding groups. Moreover, we also have used parameters of real data to test the effectiveness of our algorithm.


title: Harmonic Extension

abstract: In this paper, we consider the harmonic extension problem, which is widely used in many applications of machine learning. We find that the transitional method of graph Laplacian fails to produce a good approximation of the classical harmonic function. To tackle this problem, we propose a new method called the point integral method (PIM). We consider the harmonic extension problem from the point of view of solving PDEs on manifolds. The basic idea of the PIM method is to approximate the harmonicity using an integral equation, which is easy to be discretized from points. Based on the integral equation, we explain the reason why the transitional graph Laplacian may fail to approximate the harmonicity in the classical sense and propose a different approach which we call the volume constraint method (VCM). Theoretically, both the PIM and the VCM computes a harmonic function with convergence guarantees, and practically, they are both simple, which amount to solve a linear system. One important application of the harmonic extension in machine learning is semi-supervised learning. We run a popular semi-supervised learning algorithm by Zhu et al. over a couple of well-known datasets and compare the performance of the aforementioned approaches. Our experiments show the PIM performs the best.


title: Stochastic gradient descent methods for estimation with large data sets

abstract: We develop methods for parameter estimation in settings with large-scale data sets, where traditional methods are no longer tenable. Our methods rely on stochastic approximations, which are computationally efficient as they maintain one iterate as a parameter estimate, and successively update that iterate based on a single data point. When the update is based on a noisy gradient, the stochastic approximation is known as standard stochastic gradient descent, which has been fundamental in modern applications with large data sets. Additionally, our methods are numerically stable because they employ implicit updates of the iterates. Intuitively, an implicit update is a shrinked version of a standard one, where the shrinkage factor depends on the observed Fisher information at the corresponding data point. This shrinkage prevents numerical divergence of the iterates, which can be caused either by excess noise or outliers. Our sgd package in R offers the most extensive and robust implementation of stochastic gradient descent methods. We demonstrate that sgd dominates alternative software in runtime for several estimation problems with massive data sets. Our applications include the wide class of generalized linear models as well as M-estimation for robust regression.


title: Cash-Flow Based Dynamic Inventory Management

abstract: Small-to-medium size enterprises (SMEs), including many startup firms, need to manage interrelated flows of cash and inventories of goods. In this paper, we model a firm that can finance its inventory (ordered or manufactured) with loans in order to meet random demand which in general may not be time stationary. The firm earns interest on its cash on hand and pays interest on its debt. The objective is to maximize the expected value of the firm's %working capital at the end of a finite planning horizon. Our study shows that the optimal ordering policy is characterized by a pair of threshold variables for each period as function of the initial state of the period. Further, upper and lower bounds for the threshold values are developed using two simple-to-compute ordering policies. Based on these bounds, we provide an efficient algorithm to compute the two threshold values. Since the underlying state space is two-dimensional which leads to high computational complexity of the optimization algorithm, we also derive upper bounds for the optimal value function by reducing the optimization problem to one dimension. Subsequently, it is shown that policies of similar structure are optimal when the loan and deposit interest rates are piecewise linear functions, when there is a maximal loan limit and when unsatisfied demand is backordered. Finally, further managerial insights are provided with numerical studies.