title: Dynamic Precursors of Flares in Active Region NOAA 10486

abstract: Four different methods are applied here to study the precursors of flare activity in the Active Region NOAA 10486. Two approaches track the temporal behaviour of suitably chosen features (one, the weighted horizontal gradient WGM, is generalised form the horizontal gradient of the magnetic field, GM; another is the sum of the horizontal gradient of the magnetic field, GS, for all sunspot pairs). WGM is a photospheric indicator that is a proxy measure of magnetic non-potentiality of a specific area of the active region, i.e. it captures the temporal variation of the weighted horizontal gradient of magnetic flux summed up for the region where opposite magnetic polarities are highly mixed. The third one, referred to as the separateness parameter, S(lf), considers the overall morphology. Further, GS and S(lf) are photospheric newly defined quick-look indicators of the polarity mix of the entire active region. The fourth method is tracking the temporal variation of small x-ray flares, their times of succession and their energies observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager instrument. All approaches yield specific pre-cursory signatures for the imminence of flares.


title: The divisible sandpile at critical density

abstract: The divisible sandpile starts with i.i.d. random variables ("masses") at the vertices of an infinite, vertex-transitive graph, and redistributes mass by a local toppling rule in an attempt to make all masses at most 1. The process stabilizes almost surely if m<1 and it almost surely does not stabilize if m>1, where $m$ is the mean mass per vertex. The main result of this paper is that in the critical case m=1, if the initial masses have finite variance, then the process almost surely does not stabilize. To give quantitative estimates on a finite graph, we relate the number of topplings to a discrete biLaplacian Gaussian field.


title: When maximum-principle functions cease to exist

abstract: We consider geometric flow equations for contracting and expanding normal velocities, including powers of the Gauss curvature, of the mean curvature, and of the norm of the second fundamental form, and ask whether - after appropriate rescaling - closed strictly convex surfaces converge to spheres. To prove this, many authors use certain functions of the principal curvatures, which we call maximum-principle functions. We show when such functions cease to exist and exist, while presenting newly discovered maximum-principle functions.


title: Pre-relaxation in weakly interacting models

abstract: We consider time evolution in models close to integrable points with hidden symmetries that generate infinitely many local conservation laws that do not commute with one another. The system is expected to (locally) relax to a thermal ensemble if integrability is broken, or to a so-called generalised Gibbs ensemble if unbroken. In some circumstances expectation values exhibit quasi-stationary behaviour long before their typical relaxation time. For integrability-breaking perturbations, these are also called pre-thermalisation plateaux, and emerge e.g. in the strong coupling limit of the Bose-Hubbard model. As a result of the hidden symmetries, quasi-stationarity appears also in integrable models, for example in the Ising limit of the XXZ model. We investigate a weak coupling limit, identify a time window in which the effects of the perturbations become significant and solve the time evolution through a mean-field mapping. As an explicit example we study the XYZ spin-$\frac{1}{2}$ chain with additional perturbations that break integrability. One of the most intriguing results of the analysis is the appearance of persistent oscillatory behaviour. To unravel its origin, we study in detail a toy model: the transverse-field Ising chain with an additional nonlocal interaction proportional to the square of the transverse spin per unit length [Phys. Rev. Lett. 111, 197203 (2013)]. Despite being nonlocal, this belongs to a class of models that emerge as intermediate steps of the mean-field mapping and shares many dynamical properties with the weakly interacting models under consideration.


title: Fermion mass and mixing pattern in a minimal T7 flavor 331 model

abstract: We present a model based on the $SU(3)_{C}\otimes SU(3)_{L}\otimes U(1)_{X}$ gauge symmetry having an extra $T_{7}\otimes Z_{2}\otimes Z_{3}\otimes Z_{14} $ flavor group, which successfully describes the observed SM fermion mass and mixing pattern. In this framework, the light active neutrino masses arise via double seesaw mechanism and the observed charged fermion mass and quark mixing hierarchy is a consequence of the $Z_{2}\otimes Z_{3}\otimes Z_{14}$ symmetry breaking at very high energy. In our minimal $T_{7}$ flavor 331 model, the spectrum of neutrinos includes very light active neutrinos as well as heavy and very heavy sterile neutrinos. The model has in total 16 effective free parameters, which are fitted to reproduce the experimental values of the 18 physical observables in the quark and lepton sectors. The obtained physical observables for both quark and lepton sectors are compatible with their experimental values. The model predicts the effective Majorana neutrino mass parameter of neutrinoless double beta decay to be $% m_{\beta \beta }=$ 3 and 40 meV for the normal and the inverted neutrino spectrum, respectively. Furthermore, our model features a vanishing leptonic Dirac CP violating phase.


title: Suppression of large-scale perturbations by stiff solid

abstract: Evolution of large-scale scalar perturbations in the presence of stiff solid (solid with pressure to energy density ratio > 1/3) is studied. If the solid dominated the dynamics of the universe long enough, the perturbations could end up suppressed by as much as several orders of magnitude. To avoid too steep large-angle power spectrum of CMB, radiation must have prevailed over the solid long enough before recombination.


title: Broad frequency and amplitude control of vibration in freestanding graphene via scanning tunneling microscopy with calculated dynamic pseudo-magnetic fields

abstract: A technique to locally generate mechanical vibrations in freestanding graphene using scanning tunneling microscopy (STM) is presented. The frequency of the mechanical vibrations is tuned over nearly four decades and is centered around 10 Hz. The amplitude of the vibrations also changes over nearly three decades centered on 1 nm. The oscillating motion is generated in two ways: first, by scanning the STM tip on the surface and second, by scanning the bias voltage on the STM tip. The frequency and amplitude of the displaced freestanding graphene is quantitatively transformed to the frequency and strength of the locally generated pseudo-magnetic field for our specific geometry.


title: The Relative Lyapunov Indicators: Theory and Application to Dynamical Astronomy

abstract: A recently introduced chaos detection method, the Relative Lyapunov Indicator (RLI) is investigated in the cases of symplectic mappings and continuous Hamiltonian systems. It is shown that the RLI is an efficient numerical tool in determining the true nature of individual orbits, and in separating ordered and chaotic regions of the phase space of dynamical systems. A comparison between the RLI and some other variational indicators are presented, as well as the recent applications of the RLI to various problems of dynamical astronomy.


title: Molecfit: A general tool for telluric absorption correction II. Quantitative evaluation on ESO-VLT X-Shooter spectra

abstract: Context: Absorption by molecules in the Earth's atmosphere strongly affects ground-based astronomical observations. The resulting absorption line strength and shape depend on the highly variable physical state of the atmosphere, i.e. pressure, temperature, and mixing ratio of the different molecules involved. Usually, supplementary observations of so-called telluric standard stars (TSS) are needed to correct for this effect, which is expensive in terms of telescope time. We have developed the software package molecfit to provide synthetic transmission spectra based on parameters obtained by fitting narrow ranges of the observed spectra of scientific objects. These spectra are calculated by means of the radiative transfer code LBLRTM and an atmospheric model. In this way, the telluric absorption correction for suitable objects can be performed without any additional calibration observations of TSS. Aims: We evaluate the quality of the telluric absorption correction using molecfit with a set of archival ESO-VLT X-Shooter visible and near-infrared spectra. Methods: Thanks to the wavelength coverage from the U to the K band, X-Shooter is well suited to investigate the quality of the telluric absorption correction with respect to the observing conditions, the instrumental set-up, input parameters of the code, the signal-to-noise of the input spectrum, and the atmospheric profiles. These investigations are based on two figures of merit, I_off and I_res, that describe the systematic offsets and the remaining small-scale residuals of the corrections. We also compare the quality of the telluric absorption correction achieved with moelcfit to the classical method based on a telluric standard star. (Abridged)


title: Tuned, driven, and active soft matter

abstract: One characteristic feature of soft matter systems is their strong response to external stimuli. As a consequence they are comparatively easily driven out of their ground state and out of equilibrium, which leads to many of their fascinating properties. Here, we review illustrative examples. This review is structured by an increasing distance from the equilibrium ground state. On each level, examples of increasing degree of complexity are considered. In detail, we first consider systems that are quasi-statically tuned or switched to a new state by applying external fields. These are common liquid crystals, liquid crystalline elastomers, or ferrogels and magnetic elastomers. Next, we concentrate on systems steadily driven from outside e.g. by an imposed flow field. In our case, we review the reaction of nematic liquid crystals, of bulk-filling periodically modulated structures such as block copolymers, and of localized vesicular objects to an imposed shear flow. Finally, we focus on systems that are "active" and "self-driven". Here our range spans from idealized self-propelled point particles, via sterically interacting particles like granular hoppers, via microswimmers such as self-phoretically driven artificial Janus particles or biological microorganisms, via deformable self-propelled particles like droplets, up to the collective behavior of insects, fish, and birds. As we emphasize, similarities emerge in the features and behavior of systems that at first glance may not necessarily appear related. We thus hope that our overview will further stimulate the search for basic unifying principles underlying the physics of these soft materials out of their equilibrium ground state.


title: Probabilistic approach to the distribution of primes and to the proof of Legendre and Elliott-Halberstam conjectures

abstract: Probabilistic models for the distribution of primes in the natural numbers are constructed in the article. The author found and proved the probabilistic estimates of the deviation $R(x)=|\pi(x)- Li(x)|$. The author has analyzed the probabilistic models of the distribution of primes in the natural numbers and affirmed the validity of the probabilistic estimates of proved deviations $R(x)$ stronger than the estimates made under the assumption of Riemann conjecture. Legendre's conjecture was proved in this paper with probability arbitrarily close to 1 based on the probability estimates. Probabilistic models for the distribution of primes in the arithmetic progression $ki+l, (k,l)=1$ are also built in this paper. The author has proved the probability estimates for the deviation $R(x,k,l)=|\pi(x,k, l)-Li(x)/\varphi(k)|$. He has analyzed the probability models of the distribution of primes in the arithmetic progression and affirmed the validity of probabilistic estimates of proved deviations $R(x,k,l)$ stronger than the estimates made under the assumption of the extended Riemann conjecture. Elliott-Halberstam conjecture $\sum_{1 \leq k \leq x^a} {\max_{(k,l)=1}[R(x,k,l)]} \leq C x/\ln^A(x)$ was proved in this paper with probability arbitrarily close to 1 for all $0<a<1$ and $A>0$, based on the probability estimates.


title: The Flavour Portal to Dark Matter

abstract: We present a class of models in which dark matter (DM) is a fermionic singlet under the Standard Model (SM) gauge group but is charged under a symmetry of flavour that acts as well on the SM fermions. Interactions between DM and SM particles are mediated by the scalar fields that spontaneously break the flavour symmetry, the so-called flavons. In the case of gauged flavour symmetries, the interactions are also mediated by the flavour gauge bosons. We first discuss the construction and the generic features of this class of models. Then a concrete example with an abelian flavour symmetry is considered. We compute the complementary constraints from the relic abundance, direct detection experiments and flavour observables, showing that wide portions of the parameter space are still viable. Other possibilities like non-abelian flavour symmetries can be analysed within the same framework.


title: SDSS J013127.34-032100.1: a candidate blazar with a 11 billion solar mass black hole at $z$=5.18

abstract: The radio-loud quasar SDSS J013127.34-032100.1at a redshift z=5.18 is one of the most distant radio-loud objects. The radio to optical flux ratio (i.e. the radio-loudness) of the source is large, making it a promising blazar candidate. Its overall spectral energy distribution, completed by the X-ray flux and spectral slope derived through Target of Opportunity Swift/XRT observations, is interpreted by a non-thermal jet plus an accretion disc and molecular torus model. We estimate that its black hole mass is (1.1+-0.2)1e10 Msun. for an accretion efficiency eta=0.08, scaling roughly linearly with eta. Although there is a factor ~2 of systematic uncertainty, this black hole mass is the largest found at these redshifts in a radio loud object. We derive a viewing angle between 3 and 5 degrees. This implies that there must be other (hundreds) sources with the same black hole mass of SDSS J013127.34-032100.1, but whose jets are pointing away from Earth. We discuss the problems posed by the existence of such large black hole masses at such redshifts, especially in jetted quasars. In fact, if they are associated to rapidly spinning black holes, the accretion efficiency is high, implying a slower pace of black hole growth with respect to radio-quiet quasars.


title: New Pre-main-Sequence Stars in the Upper Scorpius Subgroup of Sco-Cen

abstract: We present 237 new spectroscopically confirmed pre-main-sequence K and M-type stars in the young Upper Scorpius subgroup of the Sco-Cen association, the near- est region of recent massive star formation. Using the Wide-Field Spectrograph at the Australian National University 2.3 m telescope at Siding Spring, we observed 397 kinematically and photometrically selected candidate members of Upper Scorpius, and identified new members by the presence of Lithium absorption. The HR-diagram of the new members shows a spread of ages, ranging from ~3-20 Myr, which broadly agrees with the current age estimates of ~5-10 Myr. We find a significant range of Li 6708 equivalent widths among the members, and a minor dependence of HR-diagram position on the measured equivalent width of the Li 6708A line, with members that appear younger having more Lithium. This could indicate the presence of either popu- lations of different age, or a spread of ages in Upper Scorpius. We also use Wide-Field Infrared Survey Explorer data to infer circumstellar disk presence in 25 of the members on the basis of infrared excesses, including two candidate transition disks. We find that 11.2+-3.4% of the M0-M2 spectral type (0.4-0.8 M_sun) Upper Sco stars display an excess that indicates the presence of a gaseous disk.


title: Gas Density Fluctuations in the Perseus Cluster: Clumping Factor and Velocity Power Spectrum

abstract: X-ray surface brightness fluctuations in the core of the Perseus Cluster are analyzed, using deep observations with the Chandra observatory. The amplitude of gas density fluctuations on different scales is measured in a set of radial annuli. It varies from 8 to 12 per cent on scales of ~10-30 kpc within radii of 30-160 kpc from the cluster center and from 9 to 7 per cent on scales of ~20-30 kpc in an outer, 60-220 kpc annulus. Using a statistical linear relation between the observed amplitude of density fluctuations and predicted velocity, the characteristic velocity of gas motions on each scale is calculated. The typical amplitudes of the velocity outside the central 30 kpc region are 90-140 km/s on ~20-30 kpc scales and 70-100 km/s on smaller scales ~7-10 kpc. The velocity power spectrum is consistent with cascade of turbulence and its slope is in a broad agreement with the slope for canonical Kolmogorov turbulence. The gas clumping factor estimated from the power spectrum of the density fluctuations is lower than 7-8 per cent for radii ~30-220 kpc from the center, leading to a density bias of less than 3-4 per cent in the cluster core. Uncertainties of the analysis are examined and discussed. Future measurements of the gas velocities with the Astro-H, Athena and Smart-X observatories will directly measure the gas density-velocity perturbation relation and further reduce systematic uncertainties in these quantities.


title: Massive holographic QCD in the Veneziano limit

abstract: QCD at finite, flavor independent quark mass is analyzed by using bottom-up holography in the Veneziano limit, where the backreaction of quarks to the gluon dynamics is fully included. The dependence on the quark mass of observables such as the bound state masses, the chiral condensate, the S-parameter, and the critical temperatures is studied. Many of the results are argued to be universal, i.e., independent of the details of the holographic model, and compared to explicit computations in the V-QCD models. The effect of adding four-fermion operators in QCD is also discussed.


title: Universal spin-triplet superconducting correlations of Majorana fermions

abstract: We establish that Majorana fermions on the boundary of topological superconductors have only spin-triplet superconducting correlations independent of whether the bulk superconducting gap is spin singlet or triplet. This is universal for time-reversal broken (respected) topological superconductors (TSCs) with an odd number of (pairs of) Majorana fermions on the boundary. Consequently, resonant Andreev reflection induced by Majorana fermions only occurs in spin-triplet channels and always injects spin-triplet Cooper pairs into the leads. This spin-triplet condensate results in a spin-orbit coupling (SOC) controlled oscillatory critical current without a $0$-$\pi$ transition in the TSC/SOC-semiconductor/TSC Josephson junction. The observation of this unique current-phase relation can serve as the definitive signal for Majorana fermions. Our study shows a technique for manipulating Majorana fermions based on their spin-triplet superconducting correlations.


title: Testing General Relativity with Present and Future Astrophysical Observations

abstract: One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.


title: Extending LHC Coverage to Light Pseudoscalar Mediators and Coy Dark Sectors

abstract: Many dark matter models involving weakly interacting massive particles (WIMPs) feature new, relatively light pseudoscalars that mediate dark matter pair annihilation into Standard Model fermions. In particular, simple models of this type can explain the gamma ray excess originating in the Galactic Center as observed by the Fermi Large Area Telescope. In many cases the pseudoscalar's branching ratio into WIMPs is suppressed, making these states challenging to detect at colliders through standard dark matter searches. Here, we study the prospects for observing these light mediator states at the LHC without exploiting missing energy techniques. While existing searches effectively probe pseudoscalars with masses between 5 - 14 GeV and above 90 GeV, the LHC reach can be extended to cover much of the interesting parameter space in the intermediate 20 - 80 GeV mass range in which the mediator can have appreciable Yukawa-like couplings to Standard Model fermions but would have escaped detection by LEP and other experiments. Models explaining the Galactic Center excess via a light pseudoscalar mediator can give rise to a promising signal in this regime through the associated production of the mediator with bottom quarks while satisfying all other existing constraints. We perform an analysis of the backgrounds and trigger efficiencies, detailing the cuts that can be used to extract the signal. A significant portion of the otherwise unconstrained parameter space of these models can be conclusively tested at the 13 TeV LHC with 100 fb$^{-1}$, and we encourage the ATLAS and CMS collaborations to extend their existing searches to this mass range.


title: IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes

abstract: In the extreme violence of merger and mass accretion, compact objects like black holes and neutron stars are thought to launch some of the most luminous outbursts of electromagnetic and gravitational wave energy in the Universe. Modeling these systems realistically is a central problem in theoretical astrophysics, but has proven extremely challenging, requiring the development of numerical relativity codes that solve Einstein's equations for the spacetime, coupled to the equations of general relativistic (ideal) magnetohydrodynamics (GRMHD) for the magnetized fluids. Over the past decade, the Illinois Numerical Relativity (ILNR) Group's dynamical spacetime GRMHD code has proven itself as a robust and reliable tool for theoretical modeling of such GRMHD phenomena. However, the code was written "by experts and for experts" of the code, with a steep learning curve that would severely hinder community adoption if it were open-sourced. Here we present IllinoisGRMHD, which is an open-source, highly-extensible rewrite of the original closed-source GRMHD code of the ILNR Group. Reducing the learning curve was the primary focus of this rewrite, with the goal of facilitating community involvement in the code's use and development, as well as the minimization of human effort in generating new science. IllinoisGRMHD also saves computer time, generating roundoff-precision identical output to the original code on adaptive-mesh grids, but nearly twice as fast at scales of hundreds to thousands of cores.


title: Resonances in retrograde circumbinary discs

abstract: We analyse the interaction of an eccentric binary with a circular coplanar circumbinary disc that rotates in a retrograde sense with respect to the binary. In the circular binary case, no Lindblad resonances lie within the disc and no Lindblad resonant torques are produced, as was previously known. By analytic means, we show that when the binary orbit is eccentric, there exist components of the gravitational potential of the binary which rotate in a retrograde sense to the binary orbit and so rotate progradely with respect to this disc, allowing a resonant interaction to occur between the binary and the disc. The resulting resonant torques distinctly alter the disc response from the circular binary case. We describe results of three-dimensional hydrodynamic simulations to explore this effect and categorise the response of the disc in terms of modes whose strengths vary as a function of binary mass ratio and eccentricity. These mode strengths are weak compared to the largest mode strengths expected in the prograde case where the binary and disc rotate in the same sense. However, for sufficiently high binary eccentricity, resonant torques open a gap in a retrograde circumbinary disc, while permitting gas inflow on to the binary via gas streams. The inflow results in a time varying accretion rate on to the binary that is modulated over the binary orbital period, as was previously found to occur in the prograde case.


title: Fluctuation effects in rotating Bose-Einstein condensates with broken $\mathrm{SU}(2)$ and $\mathrm{U}(1)\times \mathrm{U}(1)$ symmetries in the presence of intercomponent density-density interactions

abstract: Thermal fluctuations and melting transitions for rotating single-component superfluids have been intensively studied and are well understood. In contrast, the thermal effects on vortex states for two-component superfluids with density-density interaction, which have a much richer variety of vortex ground states, have been much less studied. Here, we investigate the thermal effects on vortex matter in superfluids with $\mathrm{U(1)}\times \mathrm{U(1)}$ broken symmetries and intercomponent density-density interactions, as well as the case with a larger $\mathrm{SU(2)}$ broken symmetry obtainable from the $\mathrm{U(1)}\times \mathrm{U(1)}$-symmetric case by tuning scattering lengths. In the former case we find that, in addition to first-order melting transitions, the system exhibits thermally driven phase transitions between square and hexagonal lattices. Our main result, however, concerns the case where the condensate exhibits $\mathrm{SU(2)}$-symmetry, and where vortices are not topological. At finite temperature, the system exhibits effects which do not have a counter-part in single component systems. Namely, it has a state where thermally averaged quantities show no regular vortex lattice, yet the system retains superfluid coherence along the axis of rotation. In such a state, the thermal fluctuations result in transitions between different (nearly)-degenerate vortex states without undergoing a melting transition. Our results apply to multi-component Bose-Einstein condensates, and we suggest how to experimentally detect some of these unusual effects in such systems.


title: Heating and cooling of the neutral ISM in the NGC4736 circumnuclear ring

abstract: The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736.


title: Spin Current in Junctions composed of Multi-band Superconductors with a Spin-density Wave

abstract: We calculate a nondissipative spin current and show that it can flow with or without a charge current. We consider a two-band model which can be applied to the description of Fe-based pnictides in coexistence regime of superconductivity and spin-density wave. Using quasiclassical Green's functions approach and tunneling Hamiltonian method we show that there exists a possibility to switch off the Josephson current while leaving the spin current finite. Moreover, it is possible to have the critical Josephson current and the critical spin current being proportional to each other, thus giving a possibility to measure the spin current via the Josephson current. The underlying mechanism is the interfering hopping of electron--hole pairs between different bands of the superconductors composing the junction. This is an intrinsic property of the system and provides a unique and natural way to utilize junctions made solely of pnictides in promising applications in spintronics devices.


title: On the possibility of observable signatures of leptonium from astrophysical sources

abstract: The formation of Ps in our Galaxy is well measured, and has led to important and unanswered questions on the origin of the positrons. In principle it should be possible to form analogous systems from mu and tau leptons, viz. true muonium and true tauonium. However the probability of formation for these systems is greatly reduced due to the intrinsically short lifetimes of the mu and tau leptons. Likewise, the decay of the atoms is hastened by the high probability of the constituent particles decaying. Nevertheless, if sufficient numbers of mu and tau pairs are produced in high energy astrophysical environments there may be significant production of true muonium and true tauonium, despite the small probabilities. This paper addresses this possibility. We have calculated the pair production spectra of mu and tau leptons from photon-photon annihilation and electron-positron annihilation in astrophysical environments. We have computed the cross-sections for radiative recombination and direct annihilation of the pairs, and the decay constants for the various allowable decays, and the wavelengths and energies of the recombination and annihilation signatures. In this way we have calculated the probabilities for the formation of true muonium and true tauonium, and the branching ratios for the various observable signatures. We have estimated the expected fluxes from accretion discs around microquasars and active galactic nuclei, and from interactions of jets with clouds and stars. We find that accretion discs around stellar mass black holes in our own Galaxy should have observable signatures at X-ray and gamma-ray energies that are in principle observable with current observatories.


title: A Semi-Analytical Line Transfer (SALT) model to interpret the spectra of galaxy outflows

abstract: We present a Semi-Analytical Line Transfer model, SALT, to study the absorption and re-emission line profiles from expanding galactic envelopes. The envelopes are described as a superposition of shells with density and velocity varying with the distance from the center. We adopt the Sobolev approximation to describe the interaction between the photons escaping from each shell and the remaining of the envelope. We include the effect of multiple scatterings within each shell, properly accounting for the atomic structure of the scattering ions. We also account for the effect of a finite circular aperture on actual observations. For equal geometries and density distributions, our models reproduce the main features of the profiles generated with more complicated transfer codes. Also, our SALT line profiles nicely reproduce the typical asymmetric resonant absorption line profiles observed in star-forming/starburst galaxies whereas these absorption profiles cannot be reproduced with thin shells moving at a fixed outflow velocity. We show that scattered resonant emission fills in the resonant absorption profiles, with a strength that is different for each transition. Observationally, the effect of resonant filling depends on both the outflow geometry and the size of the outflow relative to the spectroscopic aperture. Neglecting these effects will lead to incorrect values of gas covering fraction and column density. When a fluorescent channel is available, the resonant profiles alone cannot be used to infer the presence of scattered re-emission. Conversely, the presence of emission lines of fluorescent transitions reveals that emission filling cannot be neglected.


title: The Bubble-like Interior of the Core-Collapse Supernova Remnant Cassiopeia A

abstract: The death of massive stars is believed to involve aspheric explosions initiated by the collapse of an iron core. The specifics of how these catastrophic explosions proceed remain uncertain due, in part, to limited observational constraints on various processes that can introduce asymmetries deep inside the star. Here we present near-infrared observations of the young Milky Way supernova remnant Cassiopeia A, descendant of a type IIb core-collapse explosion, and a three-dimensional map of its interior, unshocked ejecta. The remnant's interior has a bubble-like morphology that smoothly connects to and helps explain the multi-ringed structures seen in the remnant's bright reverse shocked main shell of expanding debris. This internal structure may have originated from turbulent mixing processes that encouraged the development of outwardly expanding plumes of radioactive 56Ni-rich ejecta. If this is true, substantial amounts of its decay product, 56Fe, may still reside in these interior cavities.


title: Variational analysis of driven-dissipative Rydberg gases

abstract: We study the non-equilibrium steady state arising from the interplay between coherent and dissipative dynamics in strongly interacting Rydberg gases using a recently introduced variational method [H. Weimer, Phys. Rev. Lett. 114, 040402 (2015)]. We give a detailed discussion of the properties of this novel approach, and we provide a comparison with methods related to the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy. We find that the variational approach offers some intrinsic advantages, and we also show that it is able to explain the experimental results obtained in an ultracold Rydberg gas on an unprecedented quantitative level.


title: Topological Phases Protected By Reflection Symmetry and Cross-cap States

abstract: Twisting symmetries provides an efficient method to diagnose symmetry-protected topological (SPT) phases. In this paper, edge theories of (2+1)-dimensional topological phases protected by reflection as well as other symmetries are studied by twisting reflection symmetry, which effectively puts the edge theories on an unoriented spacetime, such as the Klein bottle. A key technical step taken in this paper is the use of the so-called cross-cap states, which encode entirely the unoriented nature of spacetime, and can be obtained by rearranging the spacetime geometry and exchanging the role of space and time coordinates. When the system is in a non-trivial SPT phase, we find that the corresponding cross-cap state is non-invariant under the action of the symmetries of the SPT phase, but acquires an anomalous phase. This anomalous phase, with a proper definition of a reference state, on which symmetry acts trivially, reproduces the known classification of (2+1)-dimensional bosonic and fermionic SPT phases protected by reflection symmetry, including in particular the Z_8 classification of topological crystalline superconductors protected by reflection and time-reversal symmetries.


title: Planetary Candidates Observed by Kepler V: Planet Sample from Q1-Q12 (36 Months)

abstract: The Kepler mission discovered 2842 exoplanet candidates with 2 years of data. We provide updates to the Kepler planet candidate sample based upon 3 years (Q1-Q12) of data. Through a series of tests to exclude false-positives, primarily caused by eclipsing binary stars and instrumental systematics, 855 additional planetary candidates have been discovered, bringing the total number known to 3697. We provide revised transit parameters and accompanying posterior distributions based on a Markov Chain Monte Carlo algorithm for the cumulative catalogue of Kepler Objects of Interest. There are now 130 candidates in the cumulative catalogue that receive less than twice the flux the Earth receives and more than 1100 have a radius less than 1.5 Rearth. There are now a dozen candidates meeting both criteria, roughly doubling the number of candidate Earth analogs. A majority of planetary candidates have a high probability of being bonafide planets, however, there are populations of likely false-positives. We discuss and suggest additional cuts that can be easily applied to the catalogue to produce a set of planetary candidates with good fidelity. The full catalogue is publicly available at the NASA Exoplanet Archive.


title: Coexistence of charge-density-wave and pair-density-wave orders in underdoped cuprates

abstract: We analyze incommensurate charge-density-wave (CDW) and pair-density-wave (PDW) orders with transferred momenta $(\pm Q,0)$/$(0,\pm Q)$ in underdoped cuprates within the spin-fermion model. Both orders appear due to exchange of spin fluctuations before magnetic order develops.We argue that the ordered state with the lowest energy has non-zero CDW and PDW components with the same momentum. Such a state breaks $C_4$ lattice rotational symmetry, time-reversal symmetry, and mirror symmetries. We argue that the feedback from CDW/PDW order on fermionic dispersion is consistent with ARPES data. We discuss the interplay between the CDW/PDW order and $d_{x^2-y^2}$ superconductivity and make specific predictions for experiments.


title: A Search for Fast X-ray Variability from Active Galactic Nuclei using Swift

abstract: Blazars are a class of active galactic nuclei (AGNs) known for their very rapid variabilty in the high energy regions of the electromagnetic spectrum. Despite this known fast variability, X-ray observations have generally not revealed variability in blazars with rate doubling or halving timescales less than approximately 15 min. Since its launch, the Swift X-ray Telescope has obtained 0.2-10 keV X-ray data on 143 AGNs, including blazars, through intense target of opportunity observations that can be analyzed in a multiwavelength context and used to model jet parameters, particularly during flare states. We have analyzed this broad Swift data set in a search for short timescale variability in blazars that could limit the size of the emission region in the blazar jet. While we do find several low-significance possible flares with potential indications of rapid variability, we find no strong evidence for rapid ($<$15 minutes) doubling or halving times in flares in the soft X-ray energy band for the AGNs analyzed.


title: Shocked Gas in IRAS F17207-0014: ISM Collisions and Outflows

abstract: We combine optical and near-infrared AO-assisted integral field observations of the merging ULIRG IRAS F17207-0014 from the Wide-Field Spectrograph (WiFeS) and Keck/OSIRIS. The optical emission line ratios [N II]/H$\alpha$, [S II]/H$\alpha$, and [O I]/H$\alpha$ reveal a mixing sequence of shocks present throughout the galaxy, with the strongest contributions coming from large radii (up to 100% at $\sim$5 kpc in some directions), suggesting galactic-scale winds. The near-infrared observations, which have approximately 30 times higher spatial resolution, show that two sorts of shocks are present in the vicinity of the merging nuclei: low-level shocks distributed throughout our field-of-view evidenced by an H$_{2}$/Br$\gamma$ line ratio of $\sim$0.6-4, and strong collimated shocks with a high H$_{2}$/Br$\gamma$ line ratio of $\sim$4-8, extending south from the two nuclear disks approximately 400 pc ($\sim$0.5 arcsec). Our data suggest that the diffuse shocks are caused by the collision of the interstellar media associated with the two progenitor galaxies and the strong shocks trace the base of a collimated outflow coming from the nucleus of one of the two disks.


title: Yukawa Unification and Sparticle Spectroscopy in Gauge Mediation Models

abstract: We explore the implications of t-b-tau (and b-tau) Yukawa coupling unification condition on the fundamental parameter space and sparticle spectroscopy in the minimal gauge mediated supersymmetry breaking (mGMSB) model. We find that this scenario prefers values of the CP-odd Higgs mass m_A > 1 TeV, with all colored sparticle masses above 3 TeV. These predictions will be hard to test at LHC13 but they may be testable at HE-LHC 33 TeV or a 100 TeV collider. Both t-b-tau and b-tau Yukawa coupling unifications prefer a relatively light gravitino with mass < 30 eV, which makes it a candidate hot dark matter particle. However, it cannot account for more than 15 % of the observed dark matter density.


title: Polarization effects in double open-charm production at LHCb

abstract: Double open-charm production is one of the most promising channels to disentangle single from double parton scattering (DPS) and study different properties of DPS. Several studies of the DPS contributions have been made. A missing ingredient so far has been the study of polarization effects, arising from spin correlations between the two partons inside an unpolarized proton. We investigate the impact polarization has on the double open-charm cross section. We show that the longitudinally polarized gluons can give significant contributions to the cross section, but for most of the considered kinematic region only have a moderate effect on the shape. We compare our findings to the LHCb data in the D0D0 final state, identify observables where polarization does have an impact on the distribution of the final state particles, and suggest measurements which could lead to first experimental indications of, or limits on, polarization in DPS.


title: Atlas of solar hidden photon emission

abstract: Hidden photons, gauge bosons of a U(1) symmetry of a hidden sector, can constitute the dark matter of the universe and a smoking gun for large volume compactifications of string theory. In the sub-eV mass range, a possible discovery experiment consists on searching the copious flux of these particles emitted from the Sun in a helioscope setup \`a la Sikivie. In this paper, we compute the flux of transversely polarised HPs from the Sun, a necessary ingredient for interpreting such experiments. We provide a detailed exposition of photon-HP oscillations in inhomogenous media, with special focus on resonance oscillations, which play a leading role in many cases. The region of the Sun emitting HPs resonantly is a thin spherical shell for which we justify an averaged-emission formula and which implies a distinctive morphology of the angular distribution of HPs on Earth in many cases. Low mass HPs with energies in the visible and IR have resonances very close to the photosphere where the solar plasma is not fully ionised and requires building a detailed model of solar refraction and absorption. We present results for a broad range of HP masses (from 0-1 keV) and energies (from the IR to the X-ray range), the most complete atlas of solar HP emission to date.


title: Accelerate micromagnetic simulations with GPU programming in MATLAB

abstract: A finite-difference Micromagnetic simulation code written in MATLAB is presented with Graphics Processing Unit (GPU) acceleration. The high performance of Graphics Processing Unit (GPU) is demonstrated compared to a typical Central Processing Unit (CPU) based code. The speed-up of GPU to CPU is shown to be greater than 30 for problems with larger sizes on a mid-end GPU in single precision. The code is less than 200 lines and suitable for new algorithm developing.


title: Role of the coin in the spectrum of quantum walks

abstract: The most elementary quantum walk is characterized by a 2-dimensional unitary coin flip matrix, which can be parameterized by 4 real variables. The influence of the choice of the coin flip matrix on the time evolution operator is analysed in a systematic way. By changing the coin parameters, the dispersion and asymmetry of eigenvalues of the time evolution operator can be tuned in a controlled way. The reduced eigenvectors in coin space are distributed along trajectories on the surface or inside the Bloch sphere, depending on the degeneracy of the spectrum. At certain values of the coin parameters the spectrum of the time evolution operator becomes 2-fold degenerate, but there might exist unique eigenvalues at the top and bottom of each quasi-energy band. The eigenstates corresponding to such eigenvalues are robust against arbitrary temporal variations in the bias parameter of the coin, as long as rest of the parameters remain unchanged.


title: Very large magnetoresistance in Fe$_{0.28}$TaS$_{2}$ single crystals

abstract: Magnetic moments intercalated into layered transition metal dichalcogenides are an excellent system for investigating the rich physics associated with magnetic ordering in a strongly anisotropic, strong spin-orbit coupling environment. We examine electronic transport and magnetization in Fe$_{0.28}$TaS$_{2}$, a highly anisotropic ferromagnet with a Curie temperature $T_{\mathrm{C}} \sim 68.8~$K. We find anomalous Hall data confirming a dominance of spin-orbit coupling in the magnetotransport properties of this material, and a remarkably large field-perpendicular-to-plane MR exceeding 60% at 2 K, much larger than the typical MR for bulk metals, and comparable to state-of-the-art GMR in thin film heterostructures, and smaller only than CMR in Mn perovskites or high mobility semiconductors. Even within the Fe$_x$TaS$_2$ series, for the current $x$ = 0.28 single crystals the MR is nearly $100\times$ higher than that found previously in the commensurate compound Fe$_{0.25}$TaS$_{2}$. After considering alternatives, we argue that the large MR arises from spin disorder scattering in the strong spin-orbit coupling environment, and suggest that this can be a design principle for materials with large MR.


title: Discrete solitons of spin-orbit coupled Bose-Einstein condensates in optical lattices

abstract: We study localized nonlinear excitations of a dilute Bose-Einstein condensate (BEC) with spin-orbit coupling in a deep optical lattice (OL). We use Wannier functions to derive a tight-binding model that includes the spin-orbit coupling (SOC) at the discrete level in the form of a generalized discrete nonlinear Sch\"odinger equation. Spectral properties are investigated and the existence and stability of discrete solitons and breathers with different symmetry properties with respect to the OL is demonstrated. We show that the symmetry of the modes can be changed from on-site to inter-site and to asymmetric modes simply by changing the interspecies interaction. Asymmetric modes appear to be novel modes intrinsic of the SOC.


title: Multivariate Stop loss Mixed Erlang Reinsurance risk: Aggregation, Capital allocation and Default risk

abstract: In this paper, we address the aggregation of dependent stop loss reinsurance risks where the dependence among the ceding insurer(s) risks is governed by the Sarmanov distribution and each individual risk belongs to the class of Erlang mixtures. We investigate the effects of the ceding insurer(s) risk dependencies on the reinsurer risk profile by deriving a closed formula for the distribution function of the aggregated stop loss reinsurance risk. Furthermore, diversification effects from aggregating reinsurance risks are examined by deriving a closed expression for the risk capital needed for the whole portfolio of the reinsurer and also the allocated risk capital for each business unit under the TVaR capital allocation principle. Moreover, given the risk capital that the reinsurer holds, we express the default probability of the reinsurer analytically. In case the reinsurer is in default, we determine analytical expressions for the amount of the aggregate reinsured unpaid losses and the unpaid losses of each reinsured line of business of the ceding insurer(s). These results are illustrated by numerical examples.


title: A hybrid method without extrapolation step for solving variational inequality problems

abstract: In this paper, we introduce a new method for solving variational inequality problems with monotone and Lipschitz-continuous mapping in Hilbert space. The iterative process is based on two well-known projection method and the hybrid (or outer approximation) method. However we do not use an extrapolation step in the projection method. The absence of one projection in our method is explained by slightly different choice of sets in hybrid method. We prove a strong convergence of the sequences generated by our method.


title: Effective Matsusaka's Theorem for surfaces in characteristic p

abstract: We obtain an effective version of Matsusaka's theorem for arbitrary smooth algebraic surfaces in positive characteristic, which provides an effective bound on the multiple which makes an ample line bundle D very ample. The proof for pathological surfaces is based on a Reider-type theorem. As a consequence, a Kawamata-Viehweg-type vanishing theorem is proved for arbitrary smooth algebraic surfaces in positive characteristic.


title: CR regular embeddings and immersions of compact orientable 4-manifolds into C3

abstract: We show that a compact orientable 4-manifold M has a CR regular immersion into C3 if and only if both its first Pontryagin class and its Euler characteristic vanish, and has a CR regular embedding into C3 if and only if in addition the second Stiefel-Whitney class of M vanishes.


title: Patterns of ideals of numerical semigroups

abstract: This article introduces patterns of ideals of numerical semigroups, thereby unifying previous definitions of patterns of numerical semigroups. Several results of general interest are proved. More precisely, this article presents results on the structure of the image of patterns of ideals, and also on the structure of the sets of patterns admitted by a given ideal.


title: g-Factor Modification in a Bulk InGaAs Epilayer by an In-plane Electric Field

abstract: We report on the modification of the g-factor by an in-plane electric field in an In$_{0.031}$Ga$_{0.969}$As epilayer. We performed external magnetic field scans of the Kerr rotation of the InGaAs film in order to independently determine the g-factor and the spin-orbit fields. The g-factor increases from $-0.4473\pm0.0001$ at 0 V/cm to $-0.4419\pm0.0001$ at 25 V/cm applied along the [1$\overline{1}$0] crystal axis. In addition, spatially-resolved spin measurements show a g-factor dependence on diffusive velocity. The change in g-factor with electric field can have a large effect on the determination of the internal spin-orbit and nuclear fields from Larmor precession frequency measurements.


title: Inverse problem with transmission eigenvalues for the discrete Schr\"odinger equation

abstract: The discrete Schr\"odinger equation with the Dirichlet boundary condition is considered on a half-line lattice when the potential is real valued and compactly supported. The inverse problem of recovery of the potential from the so-called transmission eigenvalues is analyzed. The Marchenko method and the Gel'fand-Levitan method are used to solve the inverse problem uniquely, except in one "unusual" case where the sum of the transmission eigenvalues is equal to a certain integer related to the support of the potential. It is shown that in the unusual case there may be a unique solution corresponding to certain sets of transmission eigenvalues, there may be a finite number of distinct potentials for some sets of transmission eigenvalues, or there may be infinitely many potentials for some sets of transmission eigenvalues. The theory presented is illustrated with several explicit examples.


title: The Beauty of Capturing Faces: Rating the Quality of Digital Portraits

abstract: Digital portrait photographs are everywhere, and while the number of face pictures keeps growing, not much work has been done to on automatic portrait beauty assessment. In this paper, we design a specific framework to automatically evaluate the beauty of digital portraits. To this end, we procure a large dataset of face images annotated not only with aesthetic scores but also with information about the traits of the subject portrayed. We design a set of visual features based on portrait photography literature, and extensively analyze their relation with portrait beauty, exposing interesting findings about what makes a portrait beautiful. We find that the beauty of a portrait is linked to its artistic value, and independent from age, race and gender of the subject. We also show that a classifier trained with our features to separate beautiful portraits from non-beautiful portraits outperforms generic aesthetic classifiers.


title: Gravitational Correction to Vacuum Polarization

abstract: We consider the gravitational correction to (electronic) vacuum polarization in the presence of a gravitational background field. The Dirac propagators for the virtual fermions are modified to include the leading gravitational correction (potential term) which corresponds to a coordinate-dependent fermion mass. The mass term is assumed to be uniform over a length scale commensurate with the virtual electron-positron pair. The on-mass shell renormalization condition ensures that the gravitational correction vanishes on the mass shell of the photon, i.e., the speed of light is unaffected by the quantum field theoretical loop correction, in full agreement with the equivalence principle. Nontrivial corrections are obtained for off-shell, virtual photons. We compare our findings to other works on generalized Lorentz transformations and combined quantum-electrodynamic gravitational corrections to the speed of light which have recently appeared in the literature.


title: Ten scenarios from early radiation to late time acceleration with a minimally coupled dark energy

abstract: We consider General Relativity with matter, radiation and a minimally coupled dark energy defined by an equation of state w. Using dynamical system method, we find the equilibrium points of such a theory assuming an expanding Universe and a positive dark energy density. Two of these points correspond to classical radiation and matter dominated epochs for the Universe. For the other points, dark energy mimics matter, radiation or accelerates Universe expansion. We then look for possible sequences of epochs describing a Universe starting with some radiation dominated epoch(s) (mimicked or not by dark energy), then matter dominated epoch(s) (mimicked or not by dark energy) and ending with an accelerated expansion. We find ten sequences able to follow this Universe history without singular behaviour of w at some saddle points. Most of them are new in dark energy literature. To get more than these ten sequences, w has to be singular at some specific saddle equilibrium points. This is an unusual mathematical property of the equation of state in dark energy literature, whose physical consequences tend to be discarded by observations. This thus distinguishes the ten above sequences from an infinity of ways to describe Universe expansion.


title: Measurement-device-independent quantum key distribution: from idea towards application

abstract: We assess the overall performance of our quantum key distribution (QKD) system implementing the measurement-device-independent (MDI) protocol using components with varying capabilities such as different single photon detectors and qubit preparation hardware. We experimentally show that superconducting nanowire single photon detectors allow QKD over a channel featuring 60 dB loss, and QKD with more than 600 bits of secret key per second (not considering finite key effects) over a 16 dB loss channel. This corresponds to 300 km and 80 km of standard telecommunication fiber, respectively. We also demonstrate that the integration of our QKD system into FPGA-based hardware (instead of state-of-the-art arbitrary waveform generators) does not impact on its performance. Our investigation allows us to acquire an improved understanding of the trade-offs between complexity, cost and system performance, which is required for future customization of MDI-QKD. Given that our system can be operated outside the laboratory over deployed fiber, we conclude that MDI-QKD is a promising approach to information-theoretic secure key distribution.


title: Generation of large-scale winds in horizontally anisotropic convection

abstract: We simulate three-dimensional, horizontally periodic Rayleigh-B\'enard convection between free-slip horizontal plates, rotating about a distant horizontal axis. When both the temperature difference between the plates and the rotation rate are sufficiently large, a strong horizontal wind is generated that is perpendicular to both the rotation vector and the gravity vector. The wind is turbulent, large-scale, and vertically sheared. Horizontal anisotropy, engendered here by rotation, appears necessary for such wind generation. Most of the kinetic energy of the flow resides in the wind, and the vertical turbulent heat flux is much lower on average than when there is no wind.


title: Arbitrary Transform Telescopes: The Generalization of Interferometry

abstract: The basic principle of astronomical interferometry is to derive the angular distribution of radiation in the sky from the Fourier transform of the electric field on the ground. What is so special about the Fourier transform? Nothing, it turns out. I consider the possibility of performing other transforms on the electric field with digital technology. The Fractional Fourier Transform (FrFT) is useful for interpreting observations of sources that are close to the interferometer (in the atmosphere for radio interferometers). Essentially, applying the FrFT focuses the array somewhere nearer than infinity. Combined with the other Linear Canonical Transforms, any homogeneous linear optical system with thin elements can be instantiated. The time variation of the electric field can also be decomposed into other bases besides the Fourier modes, which is especially useful for dispersed transients or quick pulses. I discuss why the Fourier basis is so commonly used, and suggest it is partly because most astrophysical sources vary slowly in time.


title: Iterative reconstruction of the detector response for medical gamma cameras

abstract: Statistical event reconstruction techniques can give better results for gamma cameras than the traditional centroid method. However, implementation of such techniques requires detailed knowledge of the PMT light response functions. Here we describe an iterative technique which allows to obtain the response functions from flood irradiation data without imposing strict requirements on the spatial uniformity of the event distribution. A successful application of the technique for medical gamma cameras is demonstrated using both simulated and experimental data. We show that this technique can be used for monitoring of the photomultiplier gain variations. An implementation of the iterative reconstruction technique capable of operating in real-time is also presented.


title: Two surface plasmon decay of plasma oscillations

abstract: The interaction of ultra-intense lasers with solid foils can be used to accelerate ions to high energies well exceeding 60 MeV. The non-linear relativistic motion of electrons in the intense laser radiation leads to their acceleration and later to the acceleration of ions. Ions can be accelerated from the front surface, the foil interior region, and the foil rear surface (TNSA, most widely used), or the foil may be accelerated as a whole if sufficiently thin (RPA). Here, we focus on the most widely used mechanism for laser ion-acceleration of TNSA. Starting from perfectly flat foils we show by simulations how electron filamentation at or inside the solid leads to a spatial modulations in the ions. The exact dynamics depend very sensitively on the chosen initial parameters which has a tremendous effect on electron dynamics. In the case of step-like density gradients we find evidence that suggests a two-surface-plasmon decay of plasma oscillations triggering a Raileigh-Taylor-like instability.


title: Experimental Characterization of Magnetic Materials for the Magnetic Shielding of Cryomodules in Particle Accelerators

abstract: The magnetic properties of two important passive magnetic shielding materials (A4K and Amumetal) for accelerator applications, subjected to various processing and heat treatment conditions are studied comprehensively over a wide range of temperatures: from Cryogenic to room temperature. We analyze the effect of processing on the extent of degradation of the magnetic properties of both materials and investigate the possibility of restoring these properties by re-annealing.


title: Characterizing the Star Formation of the Low-Mass SHIELD Galaxies from Hubble Space Telescope Imaging

abstract: The Survey of HI in Extremely Low-mass Dwarfs (SHIELD) is an on-going multi-wavelength program to characterize the gas, star formation, and evolution in gas-rich, very low-mass galaxies that populate the faint end of the galaxy luminosity function. The galaxies were selected from the first ~10% of the HI ALFALFA survey based on their low HI mass and low baryonic mass. Here, we measure the star-formation properties from optically resolved stellar populations for 12 galaxies using a color-magnitude diagram fitting technique. We derive lifetime average star-formation rates (SFRs), recent SFRs, stellar masses, and gas fractions. Overall, the recent SFRs are comparable to the lifetime SFRs with mean birthrate parameter of 1.4, with a surprisingly narrow standard deviation of 0.7. Two galaxies are classified as dwarf transition galaxies (dTrans). These dTrans systems have star-formation and gas properties consistent with the rest of the sample, in agreement with previous results that some dTrans galaxies may simply be low-luminosity dIrrs. We do not find a correlation between the recent star-formation activity and the distance to the nearest neighboring galaxy, suggesting that the star-formation process is not driven by gravitational interactions, but regulated internally. Further, we find a broadening in the star-formation and gas properties (i.e., specific SFRs, stellar masses, and gas fractions) compared to the generally tight correlation found in more massive galaxies. Overall, the star-formation and gas properties indicate these very low-mass galaxies host a fluctuating, non-deterministic, and inefficient star-formation process.


title: On a relation between the self-linking number and the braid index of closed braids in open books

abstract: We prove a generalization of the Jones-Kawamuro conjecture that relates the self-linking number and the braid index of closed braids, for planar open books with certain additional conditions and modifications. We show that our result is optimal in some sense by giving several counter examples for naive generalizations of the Jones-Kawamuro conjecture.


title: Per-Block-Convex Data Modeling by Accelerated Stochastic Approximation

abstract: Applications involving dictionary learning, non-negative matrix factorization, subspace clustering, and parallel factor tensor decomposition tasks motivate well algorithms for per-block-convex and non-smooth optimization problems. By leveraging the stochastic approximation paradigm and first-order acceleration schemes, this paper develops an online and modular learning algorithm for a large class of non-convex data models, where convexity is manifested only per-block of variables whenever the rest of them are held fixed. The advocated algorithm incurs computational complexity that scales linearly with the number of unknowns. Under minimal assumptions on the cost functions of the composite optimization task, without bounding constraints on the optimization variables, or any explicit information on bounds of Lipschitz coefficients, the expected cost evaluated online at the resultant iterates is provably convergent with quadratic rate to an accumulation point of the (per-block) minima, while subgradients of the expected cost asymptotically vanish in the mean-squared sense. The merits of the general approach are demonstrated in two online learning setups: (i) Robust linear regression using a sparsity-cognizant total least-squares criterion; and (ii) semi-supervised dictionary learning for network-wide link load tracking and imputation with missing entries. Numerical tests on synthetic and real data highlight the potential of the proposed framework for streaming data analytics by demonstrating superior performance over block coordinate descent, and reduced complexity relative to the popular alternating-direction method of multipliers.


title: The Joint Embedding Property and Maximal Models

abstract: We introduce the notion of a `pure` Abstract Elementary Class to block trivial counterexamples. We study classes of models of bipartite graphs and show: Main Theorem (cf. Theorem 3.5.2 and Corollary 3.5.6): If $(\lambda_i : i \le \alpha<\aleph_1)$ is a strictly increasing sequence of characterizable cardinals (Definition 2.1) whose models satisfy JEP$(<\lambda_0)$, there is an $L_{\omega_1,\omega}$ -sentence $\psi$ whose models form a pure AEC and (1) The models of $\psi$ satisfy JEP$(<\lambda_0)$, while JEP fails for all larger cardinals and AP fails in all infinite cardinals. (2) There exist $2^{\lambda_i^+}$ non-isomorphic maximal models of $\psi$ in $\lambda_i^+$, for all $i \le \alpha$, but no maximal models in any other cardinality; and (3) $\psi$ has arbitrarily large models. In particular this shows the Hanf number for JEP and the Hanf number for maximality for pure AEC with Lowenheim number $\aleph_0$ are at least $\beth_{\omega_1}$. We show that although AP$(\kappa)$ for each $\kappa$ implies the full amalgamation property, JEP$(\kappa)$ for each \kappa does not imply the full joint embedding property. We show the main combinatorial device of this paper cannot be used to extend the main theorem to a complete sentence.


title: Efficient Scheme of Experimental Quantifying non-Markovianity in High-Dimension Systems

abstract: The non-Markovianity is a prominent concept of the dynamics of the open quantum systems, which is of fundamental importance in quantum mechanics and quantum information. Despite of lots of efforts, the experimentally measuring of non-Markovianity of an open system is still limited to very small systems. Presently, it is still impossible to experimentally quantify the non-Markovianity of high dimension systems with the widely used Breuer-Laine-Piilo (BLP) trace distance measure. In this paper, we propose a method, combining experimental measurements and numerical calculations, that allow quantifying the non-Markovianity of a $N$ dimension system only scaled as $N^2$, successfully avoid the exponential scaling with the dimension of the open system in the current method. After the benchmark with a two-dimension open system, we demonstrate the method in quantifying the non-Markovanity of a high dimension open quantum random walk system.


title: Relaxed singular vectors, Jack symmetric functions and fractional level $\widehat{\mathfrak{sl}}(2)$ models

abstract: The fractional level models are (logarithmic) conformal field theories associated with affine Kac-Moody (super)algebras at certain levels $k \in \mathbb{Q}$. They are particularly noteworthy because of several longstanding difficulties that have only recently been resolved. Here, Wakimoto's free field realisation is combined with the theory of Jack symmetric functions to analyse the fractional level $\widehat{\mathfrak{sl}}(2)$ models. The first main results are explicit formulae for the singular vectors of minimal grade in relaxed Wakimoto modules. These are closely related to the minimal grade singular vectors in relaxed (parabolic) Verma modules. Further results include an explicit presentation of Zhu's algebra and an elegant new proof of the classification of simple relaxed highest weight modules over the corresponding vertex operator algebra. These results suggest that generalisations to higher rank fractional level models are now within reach.


title: Virtual Full-Duplex Buffer-Aided Relaying in the Presence of Inter-Relay Interference

abstract: In this paper, we study virtual full-duplex (FD) buffer-aided relaying to recover the loss of multiplexing gain caused by half-duplex (HD) relaying in a multiple relay network, where each relay is equipped with a buffer and multiple antennas, through joint opportunistic relay selection (RS) and beamforming (BF) design. The main idea of virtual FD buffer-aided relaying is that the source and one of the relays simultaneously transmit their own information to another relay and the destination, respectively. In such networks, inter-relay interference (IRI) is a crucial problem which has to be resolved like self-interference in the FD relaying. In contrast to previous work that neglected IRI, we propose joint RS and BF schemes taking IRI into consideration by using multiple antennas at the relays. In order to maximize average end-to-end rate, we propose a weighted sum-rate maximization strategy assuming that adaptive rate transmission is employed in both the source to relay and relay to destination links. Then, we propose several BF schemes cancelling or suppressing IRI in order to maximize the weighted sum-rate. Numerical results show that our proposed optimal, zero forcing, and minimum mean square error BF-based RS schemes asymptotically approach the ideal FD relaying upper bound when increasing the number of antennas and/or the number of relays.


title: Tensor Factorization via Matrix Factorization

abstract: Tensor factorization arises in many machine learning applications, such knowledge base modeling and parameter estimation in latent variable models. However, numerical methods for tensor factorization have not reached the level of maturity of matrix factorization methods. In this paper, we propose a new method for CP tensor factorization that uses random projections to reduce the problem to simultaneous matrix diagonalization. Our method is conceptually simple and also applies to non-orthogonal and asymmetric tensors of arbitrary order. We prove that a small number random projections essentially preserves the spectral information in the tensor, allowing us to remove the dependence on the eigengap that plagued earlier tensor-to-matrix reductions. Experimentally, our method outperforms existing tensor factorization methods on both simulated data and two real datasets.


title: Scalable Probabilistic Frames

abstract: We consider the problem of rescaling the lengths of a finite frame thereby transforming it into a tight one. Such frames are called scalable and have received a lot of attention in recent years. In this note we investigate the question in terms of probabilistic frames and give conditions under which a (discrete) probabilistic frame is scalable.


title: Self-force: Computational Strategies

abstract: Building on substantial foundational progress in understanding the effect of a small body's self-field on its own motion, the past 15 years has seen the emergence of several strategies for explicitly computing self-field corrections to the equations of motion of a small, point-like charge. These approaches broadly fall into three categories: (i) mode-sum regularization, (ii) effective source approaches and (iii) worldline convolution methods. This paper reviews the various approaches and gives details of how each one is implemented in practice, highlighting some of the key features in each case.


title: Performance Analysis of Raptor Codes under Maximum-Likelihood (ML) Decoding

abstract: Raptor codes have been widely used in many multimedia broadcast/multicast applications. However, our understanding of Raptor codes is still incomplete due to the insufficient amount of theoretical work on the performance analysis of Raptor codes, particularly under maximum-likelihood (ML) decoding, which provides an optimal benchmark on the system performance for the other decoding schemes to compare against. For the first time, this paper provides an upper bound and a lower bound, on the packet error performance of Raptor codes under ML decoding, which is measured by the probability that all source packets can be successfully decoded by a receiver with a given number of successfully received coded packets. Simulations are conducted to validate the accuracy of the analysis. More specifically, Raptor codes with different degree distribution and pre-coders, are evaluated using the derived bounds with high accuracy.


title: The Asaeda-Haagerup fusion categories

abstract: The classification of subfactors of small index revealed several new subfactors. The first subfactor above index 4, the Haagerup subfactor, is increasingly well understood and appears to lie in a (discrete) infinite family of subfactors where the Z/3Z symmetry is replaced by other finite Abelian groups. The goal of this paper is to give a similarly good description of the Asaeda-Haagerup subfactor which emerged from our study of its Brauer-Picard groupoid. More specifically, we construct a new subfactor S which is a Z/4Z x Z/2Z analogue of the Haagerup subfactor and we show that the even parts of the Asaeda-Haagerup subfactor are higher Morita equivalent to an orbifold quotient of S. This gives a new construction of the Asaeda-Haagerup subfactor which is much more symmetric and easier to work with than the original construction. As a consequence, we can settle many open questions about the Asaeda-Haagerup subfactor: calculating its Drinfel'd center, classifying all extensions of the Asaeda-Haagerup fusion categories, finding the full higher Morita equivalence class of the Asaeda-Haagerup fusion categories, and finding intermediate subfactor lattices for subfactors coming from the Asaeda-Haagerup categories. The details of the applications will be given in subsequent papers.


title: Self-deflecting plasmonic lattice solitons and sur-face modes in chirped plasmonic arrays

abstract: We show that chirped metal-dielectric waveguide arrays with focusing cubic nonlinearity can support plasmonic lattice solitons that undergo self-deflection in the transverse plane. Such lattice solitons are deeply-subwavelength self-sustained excitations, although they cover several periods of the array. Upon propagation,the excitations accelerate in the transverse plane and follow trajectories curved in the direction in which the separation between neighboring metallic layers decreases, a phenomenon that yields considerable deflection angles. The deflection angle can be controlled by varying the array chirp. We also reveal the existence of surface modes at the boundary of truncated plasmonic chirped arraythat form even in the absence of nonlinearity.


title: A solution to the Bernstein problem in the three-dimensional Heisenberg group via loop groups

abstract: In this note we present a short alternative proof for the Bernstein problem in the three-dimensional Heisenberg group ${\rm Nil}_3$ by using the loop group technique.


title: Spectra of barium, radium, and element 120; application of the combined correlation potential, singles-doubles, and configuration interaction ab initio method

abstract: We apply a version of the recently developed approach combining the correlation potential, linearized singles-doubles coupled-cluster, and the configuration interaction methods to the spectra of the heavy alkaline earths barium, radium, and element 120. Quantum electrodynamics radiative corrections are included. We have found unprecedented agreement between ab initio theory and experiment for the spectra of barium and radium, and we make accurate predictions for missing and unreliable data for all three atoms.


title: On the Convergence of Massive MIMO Systems

abstract: In this paper we examine convergence properties of massive MIMO systems with the aim of determining the number of antennas required for massive MIMO gains. We consider three characteristics of a channel matrix and study their asymptotic behaviour. Furthermore, we derive ZF SNR and MF SINR for a scenario of unequal receive powers. In our results we include the effects of spatial correlation. We show that the rate of convergence of channel metrics is much slower than that of the ZF/MF precoder properties.


title: A Big Data Architecture Design for Smart Grids Based on Random Matrix Theory

abstract: Model-based analysis tools, built on assumptions and simplifications, are difficult to handle smart grids with data characterized by 4Vs data. This paper, using random matrix theory (RMT), motivates data-driven tools to perceive the complex grids in highdimension; meanwhile, an architecture with detailed procedures is proposed. In algorithm perspective, the architecture performs a high-dimensional analysis, and compares the findings with RMT predictions to conduct anomaly detections. Mean Spectral Radius (MSR), as a statistical indicator, is defined to reflect the correlations of system data in different dimensions. In management mode perspective, a group-work mode is discussed for smart grids operation. This mode breaks through regional limitations for energy flows and data flows, and makes advanced big data analyses possible. For a specific large-scale zone-dividing system with multiple connected utilities, each site, operating under the group-work mode, is able to work out the regional MSR only with its own measured/simulated data. The large-scale interconnected system, in this way, is naturally decoupled from statistical parameters perspective, rather than from engineering models perspective. Furthermore, a comparative analysis of these distributed MSRs, even with imperceptible different raw data, will produce a contour line to detect the event and locate the source. It demonstrates that the architecture is compatible with the block calculation only using the regional small database; beyond that, this architecture, as a data-driven solution, is sensitive to system situation awareness, and practical for real large-scale interconnected systems. Five case studies and their visualizations validate the designed architecture in various fields of power systems. To our best knowledge, this study is the first attempt to apply big data technology into smart grids.


title: Overspinning a Kerr black hole: the effect of self-force

abstract: We study the scenario in which a massive particle is thrown into a rapidly rotating Kerr black hole in an attempt to spin it up beyond its extremal limit, challenging weak cosmic censorship. We work in black-hole perturbation theory, and focus on non-spinning, uncharged particles sent in on equatorial orbits. We first identify the complete parameter-space region in which overspinning occurs when back-reaction effects from the particle's self-gravity are ignored. We find, in particular, that overspinning can be achieved only with particles sent in from infinity. Gravitational self-force effects may prevent overspinning by radiating away a sufficient amount of the particle's angular momentum ("dissipative effect"), and/or by increasing the effective centrifugal repulsion, so that particles with suitable parameters never get captured ("conservative effect"). We analyze the full effect of the self-force, thereby completing previous studies by Jacobson and Sotiriou (who neglected the self-force) and by Barausse, Cardoso and Khanna (who considered the dissipative effect on a subset of orbits). Our main result is an inequality, involving certain self-force quantities, which describes a necessary and sufficient condition for the overspinning scenario to be overruled. This "censorship" condition is formulated on a certain one-parameter family of geodesics in an extremal Kerr geometry. We find that the censorship condition is insensitive to the dissipative effect (within the first-order self-force approximation used here), except for a subset of perfectly fine-tuned orbits, for which a separate censorship condition is derived. We do not obtain here the self-force input needed to evaluate either of our two conditions, but discuss the prospects for producing the necessary data using state-of-the-art numerical codes.


title: Homogenisation of a Row of Dislocation Dipoles from Discrete Dislocation Dynamics

abstract: Conventional discrete-to-continuum approaches have seen their limitation in describing the collective behaviour of the multi-polar configurations of dislocations, which are widely observed in crystalline materials. The reason is that dislocation dipoles, which play an important role in determining the mechanical properties of crystals, often get smeared out when traditional homogenisation methods are applied. To address such difficulties, the collective behaviour of a row of dislocation dipoles is studied by using matched asymptotic techniques. The discrete-to-continuum transition is facilitated by introducing two field variables respectively describing the dislocation pair density potential and the dislocation pair width. It is found that the dislocation pair width evolves much faster than the pair density. Such hierarchy in evolution time scales enables us to describe the dislocation dynamics at the coarse-grained level by an evolution equation for the slowly varying variable (the pair density) coupled with an equilibrium equation for the fast varying variable (the pair width). The time-scale separation method adopted here paves a way for properly incorporating dipole-like (zero net Burgers vector but non-vanishing) dislocation structures, known as the statistically stored dislocations (SSDs) into macroscopic models of crystal plasticity in three dimensions. Moreover, the natural transition between different equilibrium patterns found here may also shed light on understanding the emergence of the persistent slip bands (PSBs) in fatigue metals induced by cyclic loads.


title: Pointwise estimates and regularity in geometric optics and other Generated Jacobian Equations

abstract: The study of reflector surfaces in geometric optics necessitates the analysis of certain nonlinear equations of Monge-Amp\`ere type known as generated Jacobian equations. These equations, whose general existence theory has been recently developed by Trudinger go beyond the framework of optimal transport. We obtain pointwise estimates for weak solutions of such equations under a condition analogous to the A3w condition of Ma, Trudinger and Wang. Estimates of this type have played an important role in the regularity theory for optimal transport maps and were previously unknown in this context, including the important case of the near field reflector problem.


title: 1.60 m Telescope P-E

abstract: A comparative the telescope twin P-E, one is in Observat\'orio Pico dos Dias, Brazil, other in Observatoire Mont-M\'egantic, Qu\'ebec, Canada. The same project, the same beginning configuration but different stories, different solutions for the problems by the years.


title: Finite anticanonical transformations in field-antifield formalism

abstract: We study the role of arbitrary (finite) anticanonical transformations in the field-antifield formalism, and the gauge-fixing procedure based on the use of these transformations. The properties of the generating functionals of the Green functions subjected to finite anticanonical transformations are considered.


title: $L(1,1)-$ Labeling of Direct Product of Cycles

abstract: An $L(1,1)$-labeling of a graph $G$ is an assignment of labels from $\{0,1 \cdots, k \}$ to the vertices of $G$ such that two vertices that are adjacent or have a common neighbor receive distinct labels. The $\lambda_1^1-$ number, $\lambda_1^1(G)$ of $G$ is the minimum value $k$ such that $G$ admits an $L(1,1)$ labeling. We establish the $\lambda_1^1-$ numbers for direct product of cycles $C_m \times C_n$ for all positive $m, n \geq 3$, where both $m,n$ are even or when one of them is even and the other odd.


title: Generalized Simplified Variable-Scaled Min Sum LDPC decoder for irregular LDPC Codes

abstract: In this paper, we propose a novel low complexity scaling strategy of min-sum decoding algorithm for irregular LDPC codes. In the proposed method, we generalize our previously proposed simplified Variable Scaled Min-Sum (SVS-min-sum) by replacing the sub-optimal starting value and heuristic update for the scaling factor sequence by optimized values. Density evolution and Nelder-Mead optimization are used offline, prior to the decoding, to obtain the optimal starting point and per iteration updating step size for the scaling factor sequence of the proposed scaling strategy. The optimization of these parameters proves to be of noticeable positive impact on the decoding performance. We used different DVB-T2 LDPC codes in our simulation. Simulation results show the superior performance (in both WER and latency) of the proposed algorithm to other Min-Sum based algorithms. In addition to that, generalized SVS-min-sum algorithm has very close performance to LLR-SPA with much lower complexity.


title: Notes on Schubert, Grothendieck and Key Polynomials

abstract: We introduce common generalization of (double) Schubert, Grothendieck, Demazure, dual and stable Grothendieck polynomials, and Di Francesco-Zinn-Justin polynomials. Our approach is based on the study of algebraic and combinatorial properties of the reduced rectangular plactic algebra and associated Cauchy kernels.


title: On Vectorization of Deep Convolutional Neural Networks for Vision Tasks

abstract: We recently have witnessed many ground-breaking results in machine learning and computer vision, generated by using deep convolutional neural networks (CNN). While the success mainly stems from the large volume of training data and the deep network architectures, the vector processing hardware (e.g. GPU) undisputedly plays a vital role in modern CNN implementations to support massive computation. Though much attention was paid in the extent literature to understand the algorithmic side of deep CNN, little research was dedicated to the vectorization for scaling up CNNs. In this paper, we studied the vectorization process of key building blocks in deep CNNs, in order to better understand and facilitate parallel implementation. Key steps in training and testing deep CNNs are abstracted as matrix and vector operators, upon which parallelism can be easily achieved. We developed and compared six implementations with various degrees of vectorization with which we illustrated the impact of vectorization on the speed of model training and testing. Besides, a unified CNN framework for both high-level and low-level vision tasks is provided, along with a vectorized Matlab implementation with state-of-the-art speed performance.


title: Dimension reduction for the Landau-de Gennes model in planar nematic thin films

abstract: We use the method of $\Gamma$-convergence to study the behavior of the Landau-de Gennes model for a nematic liquid crystalline film in the limit of vanishing thickness. In this asymptotic regime, surface energy plays a greater role and we take particular care in understanding its influence on the structure of the minimizers of the derived two-dimensional energy. We assume general weak anchoring conditions on the top and the bottom surfaces of the film and the strong Dirichlet boundary conditions on the lateral boundary of the film. The constants in the weak anchoring conditions are chosen so as to enforce that a surface-energy-minimizing nematic $Q$-tensor has the normal to the film as one of its eigenvectors. We establish a general convergence result and then discuss the limiting problem in several parameter regimes.


title: Sequential Probability Assignment with Binary Alphabets and Large Classes of Experts

abstract: We analyze the problem of sequential probability assignment for binary outcomes with side information and logarithmic loss, where regret---or, redundancy---is measured with respect to a (possibly infinite) class of experts. We provide upper and lower bounds for minimax regret in terms of sequential complexities of the class. These complexities were recently shown to give matching (up to logarithmic factors) upper and lower bounds for sequential prediction with general convex Lipschitz loss functions (Rakhlin and Sridharan, 2015). To deal with unbounded gradients of the logarithmic loss, we present a new analysis that employs a sequential chaining technique with a Bernstein-type bound. The introduced complexities are intrinsic to the problem of sequential probability assignment, as illustrated by our lower bound. We also consider an example of a large class of experts parametrized by vectors in a high-dimensional Euclidean ball (or a Hilbert ball). The typical discretization approach fails, while our techniques give a non-trivial bound. For this problem we also present an algorithm based on regularization with a self-concordant barrier. This algorithm is of an independent interest, as it requires a bound on the function values rather than gradients.


title: Affine connections and singularities of tangent surfaces to space curves

abstract: We give the complete solution to the local diffeomorphism classification problem of generic singularities which appear in tangent surfaces, in as wider situations as possible. We interpret tangent geodesics as tangent lines whenever a (semi-)Riemannian metric, or, more generally, an affine connection is given in an ambient space of arbitrary dimension. Then, given an immersed curve, or, more generally a directed curve or a frontal curve which has well-defined tangent directions along the curve, we define the tangent surface as the ruled surface by tangent geodesics to the curve. We apply the characterization of frontal singularities found by Kokubu, Rossman, Saji, Umehara, Yamada, and Fujimori, Saji, Umehara, Yamada, and found by the first author related to the procedure of openings of singularities.


title: Only accessible information is useful: insights from gradient-mediated patterning

abstract: Information theory is gaining popularity as a tool to characterize performance of biological systems. However, information is commonly quantified without reference to whether or how a system could extract and use it; as a result, information-theoretic quantities are easily misinterpreted. Here we take the example of pattern-forming developmental systems which are commonly structured as cascades of sequential gene expression steps. Such a multi-tiered structure appears to constitute sub-optimal use of the positional information provided by the input morphogen because noise is added at each tier. However, the conventional theory fails to distinguish between the total information in a morphogen and information that can be usefully extracted and interpreted by downstream elements. We demonstrate that quantifying the information that is _accessible_ to the system naturally explains the prevalence of multi-tiered network architectures as a consequence of the noise inherent to the control of gene expression. We support our argument with empirical observations from patterning along the major body axis of the fruit fly embryo. Our results exhibit the limitations of the standard information-theoretic characterization of biological signaling and illustrate how they can be resolved.


title: Matching densities for Galois representations

abstract: Given a pair of n-dimensional complex Galois representations over Q, we define their matching density to be the density, if it exists, of the set of places at which the traces of Frobenius of the two Galois representations are equal. We will show that the set of matching densities of such pairs of irreducible Galois representations (for all n) is dense in the interval [0, 1]. Under the strong Artin conjecture, this also implies the corresponding statement for cuspidal automorphic representations.


title: Interatomic potentials for ionic systems with density functional accuracy based on charge densities obtained by a neural network

abstract: Based on an analysis of the short range chemical environment of each atom in a system, standard machine learning based approaches to the construction of interatomic potentials aim at determining directly the central quantity which is the total energy. This prevents for instance an accurate description of the energetics of systems where long range charge transfer is important as well as of ionized systems. We propose therefore not to target directly with machine learning methods the total energy but an intermediate physical quantity namely the charge density, which then in turn allows to determine the total energy. By allowing the electronic charge to distribute itself in an optimal way over the system, we can describe not only neutral but also ionized systems with unprecedented accuracy. We demonstrate the power of our approach for both neutral and ionized NaCl clusters where charge redistribution plays a decisive role for the energetics. We are able to obtain chemical accuracy, i.e. errors of less than a milli Hartree per atom compared to the reference density functional results. The introduction of physically motivated quantities which are determined by the short range atomic environment via a neural network leads also to an increased stability of the machine learning process and transferability of the potential.


title: Maximal $\bf L^p$ analysis of finite element solutions for parabolic equations with nonsmooth coefficients in convex polyhedra

abstract: The paper is concerned with Galerkin finite element solutions for parabolic equations in a convex polygon or polyhehron with a diffusion coefficient in $W^{1,N+\epsilon}$ for some $\epsilon>0$, where $N$ denotes the dimension of the domain. We prove the analyticity of the semigroup generated by the discrete elliptic operator, the discrete maximal $L^p$ regularity and the optimal $L^p$ error estimate of the finite element solution for the parabolic equation.


title: Superconducting properties of sulfur-doped iron selenide

abstract: The recent discovery of high-temperature superconductivity in single-layer iron selenide has generated significant experimental interest for optimizing the superconducting properties of iron-based superconductors through the lattice modification. For simulating the similar effect by changing the chemical composition due to S doping, we investigate the superconducting properties of high-quality single crystals of FeSe$_{1-x}$S$_{x}$ ($x$=0, 0.04, 0.09, and 0.11) using magnetization, resistivity, the London penetration depth, and low temperature specific heat measurements. We show that the introduction of S to FeSe enhances the superconducting transition temperature $T_{c}$, anisotropy, upper critical field $H_{c2}$, and critical current density $J_{c}$. The upper critical field $H_{c2}(T)$ and its anisotropy are strongly temperature dependent, indicating a multiband superconductivity in this system. Through the measurements and analysis of the London penetration depth $\lambda _{ab}(T)$ and specific heat, we show clear evidence for strong coupling two-gap $s$-wave superconductivity. The temperature-dependence of $\lambda _{ab}(T)$ calculated from the lower critical field and electronic specific heat can be well described by using a two-band model with $s$-wave-like gaps. We find that a $d$-wave and single-gap BCS theory under the weak-coupling approach can not describe our experiments. The change of specific heat induced by the magnetic field can be understood only in terms of multiband superconductivity.


title: On A Relative Mumford-Newstead Theorem

abstract: In this paper we prove a relative version of the classical Mumford-Newstead theorem for a family of smooth curves degenerating to a reducible curve with a simple node. We also prove a Torelli-type theorem by showing that certain moduli spaces of torsion-free sheaves on a reducible curve allows us to recover the curve from the moduli space.


title: Finding Connected Dense $k$-Subgraphs

abstract: Given a connected graph $G$ on $n$ vertices and a positive integer $k\le n$, a subgraph of $G$ on $k$ vertices is called a $k$-subgraph in $G$. We design combinatorial approximation algorithms for finding a connected $k$-subgraph in $G$ such that its density is at least a factor $\Omega(\max\{n^{-2/5},k^2/n^2\})$ of the density of the densest $k$-subgraph in $G$ (which is not necessarily connected). These particularly provide the first non-trivial approximations for the densest connected $k$-subgraph problem on general graphs.


title: Structure-Based Self-Triggered Consensus in Networks of Multiagents with Switching Topologies

abstract: In this paper, we propose a new self-triggered consensus algorithm in networks of multi-agents. Different from existing works, which are based on the observation of states, here, each agent determines its next update time based on its coupling structure. Both centralized and distributed approaches of the algorithms have been discussed. By transforming the algorithm to a proper discrete-time systems without self delays, we established a new analysis framework to prove the convergence of the algorithm. Then we extended the algorithm to networks with switching topologies, especially stochastically switching topologies. Compared to existing works, our algorithm is easier to understand and implement. It explicitly provides positive lower and upper bounds for the update time interval of each agent based on its coupling structure, which can also be independently adjusted by each agent according to its own situation. Our work reveals that the event/self triggered algorithms are essentially discrete and more suitable to a discrete analysis framework. Numerical simulations are also provided to illustrate the theoretical results.