title: First Data Release of the All-sky NOAO Source Catalog

abstract: Most of the sky has been imaged with NOAO's telescopes from both hemispheres. While the large majority of these data were obtained for PI-led projects and almost all of the images are publicly available, only a small fraction have been released to the community via well-calibrated and easily accessible catalogs. We are remedying this by creating a catalog of sources from most of the public data taken on the CTIO-4m+DECam and the KPNO-4m+Mosaic3. This catalog, called the NOAO Source Catalog (NSC), contains over 2.9 billion unique objects, 34 billion individual source measurements, covers ~30,000 square degrees of the sky, has depths of ~23rd magnitude in most broadband filters with ~1-2% photometric precision, and astrometric accuracy of ~7 mas. In addition, ~2 billion objects and ~21,000 square degrees of sky have photometry in three or more bands. The NSC will be useful for exploring stellar streams, dwarf satellite galaxies, QSOs, high-proper motion stars, variable stars and other transients. The NSC catalog is publicly available via the NOAO Data Lab service.


title: Weak Lensing Peaks in Simulated Light-Cones: Investigating the Coupling between Dark Matter and Dark Energy

abstract: In this paper, we study the statistical properties of weak lensing peaks in light-cones generated from cosmological simulations. In order to assess the prospects of such observable as a cosmological probe, we consider simulations that include interacting Dark Energy (hereafter DE) models with coupling term between DE and Dark Matter. Cosmological models that produce a larger population of massive clusters have more numerous high signal-to-noise peaks; among models with comparable numbers of clusters those with more concentrated haloes produce more peaks. The most extreme model under investigation shows a difference in peak counts of about $20\%$ with respect to the reference $\mathrm{\Lambda}$CDM model. We find that peak statistics can be used to distinguish a coupling DE model from a reference one with the same power spectrum normalisation. The differences in the expansion history and the growth rate of structure formation are reflected in their halo counts, non-linear scale features and, through them, in the properties of the lensing peaks. For a source redshift distribution consistent with the expectations of future space-based wide field surveys, we find that typically seventy percent of the cluster population contributes to weak-lensing peaks with signal-to-noise ratios larger than two, and that the fraction of clusters in peaks approaches one-hundred percent for haloes with redshift z$\leq$0.5. Our analysis demonstrates that peak statistics are an important tool for disentangling DE models by accurately tracing the structure formation processes as a function of the cosmic time.


title: Dark matter contraction and stellar-mass-to-light ratio gradients in massive early-type galaxies

abstract: We present models for the dark and luminous mass structure of 12 strong lensing early-type galaxies (ETGs). We combine pixel-based modelling of multiband HST/ACS imaging with Jeans modelling of kinematics obtained from Keck/ESI spectra to disentangle the dark and luminous contributions to the mass. Assuming a gNFW profile for the dark matter halo and a spatially constant stellar-mass-to-light ratio $\Upsilon_{\star}$ for the baryonic mass, we infer distributions for $\Upsilon_{\star}$ consistent with IMFs that are heavier than the Milky Way's (with a global mean mismatch parameter relative to a Chabrier IMF $\mu_{\alpha c} = 1.80 \pm 0.14$) and halo inner density slopes which span a large range but are generally cuspier than the dark-matter-only prediction ($\mu_{\gamma'} = 2.01_{-0.22}^{+0.19}$). We investigate possible reasons for overestimating the halo slope, including the neglect of spatially varying stellar-mas-to-light ratios and/or stellar orbital anisotropy, and find that a quarter of the systems prefer radially declining stellar-mass-to-light ratio gradients, but that the overall effect on our inference on the halo slope is small. We suggest a coherent explanation of these results in the context of inside-out galaxy growth, and that the relative importance of different baryonic processes in shaping the dark halo may depend on halo environment.


title: Polarization of K-shell dielectronic recombination satellite lines of Fe XIX-XXV and its application for diagnostics of anisotropies of hot plasmas

abstract: We present a systematic measurement of the X-ray emission asymmetries in the K-shell dielectronic, trielectronic, and quadruelectronic recombination of free electrons into highly charged ions. Iron ions in He-like through O-like charge states were produced in an electron beam ion trap, and the electron-ion collision energy was scanned over the recombination resonances. Two identical X-ray detectors mounted head-on and side-on with respect to the electron beam propagation recorded X-rays emitted in the decay of resonantly populated states. The degrees of linear polarization of X-rays inferred from observed emission asymmetries benchmark distorted-wave predictions of the Flexible Atomic Code (FAC) for several dielectronic recombination satellite lines. The present method also demonstrates its applicability for diagnostics of energy and direction of electron beams inside hot anisotropic plasmas. Both experimental and theoretical data can be used for modeling of hot astrophysical and fusion plasmas.

categories: astro-ph.HE astro-ph.IM astro-ph.SR physics.atom-ph physics.plasm-ph


title: Multiband counterparts of two eclipsing ultraluminous X-ray sources in M 51

abstract: We present the discovery and interpretation of ionized nebulae around two ultraluminous X-ray sources in M 51; both sources share the rare property of showing X-ray eclipses by their companion stars, and are therefore prime targets for follow-up studies. Using archival Hubble Space Telescope images, we found an elongated, 100-pc-long emission-line structure associated with one X-ray source (CXOM51 J132940.0+471237; ULX-1 for simplicity), and a more circular, ionized nebula at the location of the second source (CXOM51 J132939.5+471244; ULX-2 for simplicity). We observed both nebulae with the Large Binocular Telescope's Multi-Object Double Spectrograph. From our analysis of the optical spectra, we argue that the gas in the ULX-1 bubble is shock-ionized, consistent with the effect of a jet with a kinetic power of $\approx$2 $\times 10^{39}$ erg s$^{-1}$. Additional X-ray photo-ionization may also be present, to explain the strength of high-ionization lines such as He II $\lambda 4686$ and [Ne V] $\lambda 3426$. On the other hand, the emission lines from the ULX-2 bubble are typical for photoionization by normal O stars suggesting that the nebula is actually an H II region not physically related to the ULX but is simply a chance alignment. From archival Very Large Array data, we also detect spatially extended, steep-spectrum radio emission at the location of the ULX-1 bubble (consistent with its jet origin), but no radio counterpart for ULX-2 (consistent with the lack of shock-ionized gas around that source).


title: The KBC Void: Consistency with Supernovae Type Ia and the Kinematic SZ Effect in a $\Lambda$LTB Model

abstract: There is substantial and growing observational evidence from the normalized luminosity density in the near-infrared that the local universe is under-dense on scales of several hundred Megaparsecs. We test whether our parameterization of the observational data of such a "void" is compatible with the latest supernovae type Ia data and with constraints from line-of-sight peculiar velocity motions of galaxy clusters with respect to the cosmic microwave background rest frame, known as the linear kinematic Sunyaev-Zel'dovich (kSZ) effect. Our study is based on the large local void (LLV) radial profile observed by Keenan, Barger, and Cowie (KBC) and a theoretical void description based on the Lema\^itre-Tolman-Bondi model with a nonzero cosmological constant ($\Lambda$LTB). We find consistency with the measured luminosity distance-redshift relation on radial scales relevant to the KBC LLV through a comparison with 217 low-redshift supernovae type Ia over the redshift range $0.0233 < z < 0.15$. We assess the implications of the KBC LLV in light of the tension between "local" and "cosmic" measurements of the Hubble constant, $H_{0}$. We find that when the existence of the KBC LLV is fully accounted for, this tension is reduced from $3.4\sigma$ to $2.75\sigma$. We find that previous linear kSZ constraints, as well as new ones from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT), are fully compatible with the existence of the KBC LLV.


title: Beacons into the Cosmic Dark Ages: Boosted transmission of Ly$\alpha$ from UV bright galaxies at $z \gtrsim 7$

abstract: Recent detections of Lyman alpha (Ly$\alpha$) emission from $z>7.5$ galaxies were somewhat unexpected given a dearth of previous non-detections in this era when the intergalactic medium (IGM) is still highly neutral. But these detections were from UV bright galaxies, which preferentially live in overdensities which reionize early, and have significantly Doppler-shifted Ly$\alpha$ line profiles emerging from their interstellar media (ISM), making them less affected by the global IGM state. Using a combination of reionization simulations and empirical ISM models we show, as a result of these two effects, UV bright galaxies in overdensities have $>2\times$ higher transmission through the $z\sim7$ IGM than typical field galaxies, and this boosted transmission is enhanced as the neutral fraction increases. The boosted transmission is not sufficient to explain the observed high Ly$\alpha$ fraction of $M_\mathrm{UV} \lesssim -22$ galaxies (Stark et al. 2017), suggesting Ly$\alpha$ emitted by these galaxies must be stronger than expected due to enhanced production and/or selection effects. Despite the bias of UV bright galaxies to reside in overdensities we show Ly$\alpha$ observations of such galaxies can accurately measure the global neutral hydrogen fraction, particularly when Ly$\alpha$ from UV faint galaxies is extinguished, making them ideal candidates for spectroscopic follow-up into the cosmic Dark Ages.


title: CHANG-ES X: Spatially-resolved Separation of Thermal Contribution from Radio Continuum Emission in Edge-on Galaxies

abstract: We analyze the application of star formation rate (SFR) calibrations using H$\alpha$ and 22 micron infrared imaging data in predicting the thermal radio component for a test sample of 3 edge-on galaxies (NGC 891, NGC 3044, and NGC 4631) in the Continuum Halos in Nearby Galaxies -- an EVLA Survey (CHANG-ES). We use a mixture of H$\alpha$ and 24 micron calibration from Calzetti et al. (2007), and a linear 22 micron only calibration from Jarrett et al. (2013) on the test sample. We apply these relations on a pixel-to-pixel basis to create thermal prediction maps in the two CHANG-ES bands: L- and C-band (1.5 GHz and 6.0 GHz, respectively). We analyze the resulting non-thermal spectral index maps, and find a characteristic steepening of the non-thermal spectral index with vertical distance from the disk after application of all methods. We find possible evidence of extinction in the 22 micron data as compared to 70 micron Spitzer Multband Imaging Photometer (MIPS) imaging in NGC 891. We analyze a larger sample of edge-on and face-on galaxy 25 micron to 100 micron flux ratios, and find that the ratios for edge-ons are systematically lower by a factor of 1.36, a result we attribute to excess extinction in the mid-IR in edge-ons. We introduce a new calibration for correcting the H$\alpha$ luminosity for dust when galaxies are edge-on or very dusty.


title: The high energy fate of the minimal Goldstone Higgs

abstract: We consider a minimal model where the Higgs boson arises as an elementary pseudo-Nambu-Goldstone boson. The model is based on an extended scalar sector with global SO(5)/SO(4) symmetry. To achieve the correct electroweak symmetry breaking pattern, the model is augmented either with an explicit symmetry breaking term or an extra singlet scalar field. We consider separately both of these possibili- ties. We fit the model with the known particle spectrum at the electroweak scale and extrapolate to high energies using renormalization group. We find that the model can remain stable and perturbative up to the Planck scale provided that the heavy beyond Standard Model scalar states have masses in a narrow interval around 3 TeV.


title: Active processes in one dimension

abstract: We consider the thermal and athermal overdamped motion of particles in 1D geometries where discrete internal degrees of freedom (spin) are coupled with the translational motion. Adding a driving velocity that depends on the time-dependent spin constitutes the simplest model of active particles (run-and-tumble processes) where the violation of the equipartition principle and of the Sutherland-Einstein relation can be studied in detail even when there is generalized reversibility. We give an example (with four spin values) where the irreversibility of the translational motion manifests itself only in higher-order (than two) time correlations. We derive a generalized telegraph equation as the Smoluchowski equation for the spatial density for an arbitrary number of spin values. We also investigate the Arrhenius exponential law for run-and-tumble particles, due to their activity the slope of the potential becomes important in contrast to the passive diffusion case and activity enhances the escape from a potential well (if that slope is high enough). Finally, in the absence of a driving velocity, the presence of internal currents such as in the chemistry of molecular motors may be transmitted to the translational motion and the internal activity is crucial for the direction of the emerging spatial current.


title: The Stagger-grid: A grid of 3D stellar atmosphere models V. Synthetic stellar spectra and broad-band photometry

abstract: Context. The surface structures and dynamics of cool stars are characterized by the presence of convective motions and turbulent flows which shape the emergent spectrum. Aims. We used realistic three-dimensional radiative hydrodynamical simulations from the Stagger-grid to calculate synthetic spectra with the radiative transfer code Optim3D for stars with different stellar parameters to predict photometric colors and convective velocity shifts. Methods. We calculated spectra from 1000 to 200 000 A with a constant resolving power of 20 000 and from 8470 and 8710 A (Gaia Radial Velocity Spectrometer - RVS - spectral range), with a constant resolving power of 300 000. Results. We used synthetic spectra to compute theoretical colors in the Johnson-Cousins UBV(RI), SDSS, 2MASS, Gaia, SkyMapper, Stromgren systems, and HST-WFC3. We showed that 1D versus 3D differences are limited to a small percent except for the narrow filters that span the optical and UV region of the spectrum. In addition, we derived the effect of the convective velocity fields on selected Fe I lines. We found the overall convective shift for 3D simulations with respect to the reference 1D hydrostatic models, revealing line shifts of between -0.235 and +0.361 km/s. We showed a net correlation of the convective shifts with the effective temperature: lower effective temperatures denote redshifts and higher effective temperatures denote blueshifts. We conclude that the extraction of accurate radial velocities from RVS spectra need an appropriate wavelength correction from convection shifts. Conclusions. The use of realistic 3D hydrodynamical stellar atmosphere simulations has a small but significant impact on the predicted photometry compared with classical 1D hydrostatic models for late-type stars. We make all the spectra publicly available for the community through the POLLUX database.


title: Verifying and Synthesizing Constant-Resource Implementations with Types

abstract: We propose a novel type system for verifying that programs correctly implement constant-resource behavior. Our type system extends recent work on automatic amortized resource analysis (AARA), a set of techniques that automatically derive provable upper bounds on the resource consumption of programs. We devise new techniques that build on the potential method to achieve compositionality, precision, and automation. A strict global requirement that a program always maintains constant resource usage is too restrictive for most practical applications. It is sufficient to require that the program's resource behavior remain constant with respect to an attacker who is only allowed to observe part of the program's state and behavior. To account for this, our type system incorporates information flow tracking into its resource analysis. This allows our system to certify programs that need to violate the constant-time requirement in certain cases, as long as doing so does not leak confidential information to attackers. We formalize this guarantee by defining a new notion of resource-aware noninterference, and prove that our system enforces it. Finally, we show how our type inference algorithm can be used to synthesize a constant-time implementation from one that cannot be verified as secure, effectively repairing insecure programs automatically. We also show how a second novel AARA system that computes lower bounds on resource usage can be used to derive quantitative bounds on the amount of information that a program leaks through its resource use. We implemented each of these systems in Resource Aware ML, and show that it can be applied to verify constant-time behavior in a number of applications including encryption and decryption routines, database queries, and other resource-aware functionality.


title: Power Kripke-Platek set theory and the axiom of choice

abstract: Whilst Power Kripke-Platek set theory, KPP, shares many properties with ordinary Kripke-Platek set theory, KP, in several ways it behaves quite differently from KP. This is perhaps most strikingly demonstrated by a result, due to Mathias, to the effect that adding the axiom of constructibility to KPP gives rise to a much stronger theory, whereas in the case of KP the constructible hierarchy provides an inner model, so that KP and KP+V=L have the same strength. This paper will be concerned with the relationship between KPP and KPP plus the axiom of choice or even the global axiom of choice, GAC. Since L is the standard vehicle to furnish a model in which this axiom holds, the usual argument for demonstrating that the addition of AC or GAC to KPP does not increase proof-theoretic strength does not apply in any obvious way. Among other tools, the paper uses techniques from ordinal analysis to show that KPP+GAC has the same strength as KPP, thereby answering a question of Mathias. Moreover, it is shown that KPP+GAC is conservative over KPP for Pi-1-4 statements of analysis. The method of ordinal analysis for theories with power set was developed in an earlier paper. The technique allows one to compute witnessing information from infinitary proofs, providing bounds for the transfinite iterations of the power set operation that are provable in a theory. As the theory KPP+GAC provides a very useful tool for defining models and realizability models of other theories that are hard to construct without access to a uniform selection mechanism, it is desirable to determine its exact proof-theoretic strength. This knowledge can for instance be used to determine the strength of Feferman's operational set theory with power set operation as well as constructive Zermelo-Fraenkel set theory with the axiom of choice.


title: Natural proof of the characterization of relatively compact families in $L^p-$spaces on locally compact groups

abstract: In the paper we look for an elegant proof of the characterization of relatively compact families in $L^p-$spaces. At first glance, the suggested approach may seem convoluted and lengthy, but we spare no effort to argue that our proof is in fact more natural than the ones existent in the literature. The key idea is that the three properties which characterize relative compactness in $L^p-$spaces ($L^p-$boundedness, $L^p-$equicontinuity and $L^p-$equivanishing) are "preserved" (or rather "inherited") when the family $\Ffamily\subset L^p$ is convolved with a continuous and compactly supported function. This new family turns out to be relatively compact in $C_0-$space and it remains to be demonstrated that relative compactness in $C_0-$space implies relative compactness of the original family $\Ffamily.$


title: Clustering with Outlier Removal

abstract: Cluster analysis and outlier detection are strongly coupled tasks in data mining area. Cluster structure can be easily destroyed by few outliers; on the contrary, outliers are defined by the concept of cluster, which are recognized as the points belonging to none of the clusters. Unfortunately, most existing studies do not notice the coupled relationship between these two task and handle them separately. In light of this, we consider the joint cluster analysis and outlier detection problem, and propose the Clustering with Outlier Removal (COR) algorithm. Generally speaking, the original space is transformed into the binary space via generating basic partitions in order to define clusters. Then an objective function based Holoentropy is designed to enhance the compactness of each cluster with a few outliers removed. With further analyses on the objective function, only partial of the problem can be handled by K-means optimization. To provide an integrated solution, an auxiliary binary matrix is nontrivally introduced so that COR completely and efficiently solves the challenging problem via a unified K-means-- with theoretical supports. Extensive experimental results on numerous data sets in various domains demonstrate the effectiveness and efficiency of COR significantly over state-of-the-art methods in terms of cluster validity and outlier detection. Some key factors in COR are further analyzed for practical use. Finally, an application on flight trajectory is provided to demonstrate the effectiveness of COR in the real-world scenario.


title: Knowledge-based Word Sense Disambiguation using Topic Models

abstract: Word Sense Disambiguation is an open problem in Natural Language Processing which is particularly challenging and useful in the unsupervised setting where all the words in any given text need to be disambiguated without using any labeled data. Typically WSD systems use the sentence or a small window of words around the target word as the context for disambiguation because their computational complexity scales exponentially with the size of the context. In this paper, we leverage the formalism of topic model to design a WSD system that scales linearly with the number of words in the context. As a result, our system is able to utilize the whole document as the context for a word to be disambiguated. The proposed method is a variant of Latent Dirichlet Allocation in which the topic proportions for a document are replaced by synset proportions. We further utilize the information in the WordNet by assigning a non-uniform prior to synset distribution over words and a logistic-normal prior for document distribution over synsets. We evaluate the proposed method on Senseval-2, Senseval-3, SemEval-2007, SemEval-2013 and SemEval-2015 English All-Word WSD datasets and show that it outperforms the state-of-the-art unsupervised knowledge-based WSD system by a significant margin.


title: On randomized counting versus randomised decision

abstract: We study the question of which counting problems admit f.p.r.a.s., under a structural complexity perspective. Since problems in #P with NP-complete decision version do not admit f.p.r.a.s. (unless NP = RP), we study subclasses of #P, having decision version either in P or in RP. We explore inclusions between these subclasses and we present all possible worlds with respect to NP v.s. RP and RP v.s. P.


title: Automated data processing architecture for the Gemini Planet Imager Exoplanet Survey

abstract: The Gemini Planet Imager Exoplanet Survey (GPIES) is a multi-year direct imaging survey of 600 stars to discover and characterize young Jovian exoplanets and their environments. We have developed an automated data architecture to process and index all data related to the survey uniformly. An automated and flexible data processing framework, which we term the Data Cruncher, combines multiple data reduction pipelines together to process all spectroscopic, polarimetric, and calibration data taken with GPIES. With no human intervention, fully reduced and calibrated data products are available less than an hour after the data are taken to expedite follow-up on potential objects of interest. The Data Cruncher can run on a supercomputer to reprocess all GPIES data in a single day as improvements are made to our data reduction pipelines. A backend MySQL database indexes all files, which are synced to the cloud, and a front-end web server allows for easy browsing of all files associated with GPIES. To help observers, quicklook displays show reduced data as they are processed in real-time, and chatbots on Slack post observing information as well as reduced data products. Together, the GPIES automated data processing architecture reduces our workload, provides real-time data reduction, optimizes our observing strategy, and maintains a homogeneously reduced dataset to study planet occurrence and instrument performance.


title: Local Mixing Time: Distributed Computation and Applications

abstract: The mixing time of a graph is an important metric, which is not only useful in analyzing connectivity and expansion properties of the network, but also serves as a key parameter in designing efficient algorithms. We introduce a new notion of mixing of a random walk on a (undirected) graph, called local mixing. Informally, the local mixing with respect to a given node $s$, is the mixing of a random walk probability distribution restricted to a large enough subset of nodes --- say, a subset of size at least $n/\beta$ for a given parameter $\beta$ --- containing $s$. The time to mix over such a subset by a random walk starting from a source node $s$ is called the local mixing time with respect to $s$. The local mixing time captures the local connectivity and expansion properties around a given source node and is a useful parameter that determines the running time of algorithms for partial information spreading, gossip etc. Our first contribution is formally defining the notion of local mixing time in an undirected graph. We then present an efficient distributed algorithm which computes a constant factor approximation to the local mixing time with respect to a source node $s$ in $\tilde{O}(\tau_s)$ rounds, where $\tau_s$ is the local mixing time w.r.t $s$ in an $n$-node regular graph. This bound holds when $\tau_s$ is significantly smaller than the conductance of the local mixing set (i.e., the set where the walk mixes locally); this is typically the interesting case where the local mixing time is significantly smaller than the mixing time (with respect to $s$). We also present a distributed algorithm that computes the exact local mixing time in $\tilde{O}(\tau_s \mathcal{D})$ rounds, where $\mathcal{D} =\min\{\tau_s, D\}$ and $D$ is the diameter of the graph. We further show that local mixing time tightly characterizes the complexity of partial information spreading.


title: Phonon networks with SiV centers in diamond waveguides

abstract: We propose and analyze a novel realization of a solid-state quantum network, where separated silicon-vacancy centers are coupled via the phonon modes of a quasi-1D diamond waveguide. In our approach, quantum states encoded in long-lived electronic spin states can be converted into propagating phonon wavepackets and be reabsorbed efficiently by a distant defect center. Our analysis shows that under realistic conditions, this approach enables the implementation of high-fidelity, scalable quantum communication protocols within chip-scale spin-qubit networks. Apart from quantum information processing, this setup constitutes a novel waveguide QED platform, where strong-coupling effects between solid-state defects and individual propagating phonons can be explored at the quantum level.


title: Gamma-ray Emitting Narrow Line Seyfert 1 Galaxies in The Sloan Digital Sky Survey

abstract: The detection of significant $\gamma$-ray emission from radio-loud narrow line Seyfert 1 (NLSy1s) galaxies enables us to study jets in environments different than those in blazars. However, due to the small number of known $\gamma$-ray emitting NLSy1 ($\gamma$-NLSy1) galaxies, a comprehensive study could not be performed. Here we report the first detection of significant $\gamma$-ray emission from four active galactic nuclei (AGN), recently classified as NLSy1 from their Sloan Digital Sky Survey (SDSS) optical spectrum. Three flat spectrum radio quasars (FSRQs) present in the third Large Area Telescope AGN catalog (3LAC) are also found as $\gamma$-NLSy1 galaxies. Comparing the $\gamma$-ray properties of these objects with 3LAC blazars reveals their spectral shapes to be similar to FSRQs, however, with low $\gamma$-ray luminosity ($\lesssim10^{46-47}$ erg s$^{-1}$). In the Wide-field Infrared Survey Explorer color-color diagram, these objects occupy a region mainly populated by FSRQs. Using the H$_{\beta}$ emission line parameters, we find that on average $\gamma$-NLSy1 have smaller black hole masses than FSRQs at similar redshifts. In the low-resolution SDSS image of one of the $\gamma$-NLSy1 source, we find the evidence of an extended structure. We conclude by noting that overall many observational properties of $\gamma$-NLSy1 sources are similar to FSRQs and therefore, these objects could be their low black hole mass counterparts, as predicted in the literature.


title: Rankin-Cohen brackets and Serre derivatives as Poincar\'e series

abstract: We give expressions for the Serre derivatives of Eisenstein and Poincar\'e series as well as their Rankin-Cohen brackets with arbitrary modular forms in terms of the Poincar\'e averaging construction, and derive several identities for the Ramanujan tau function as applications.


title: Pressurized rf cavities in ionizing beams

abstract: A muon collider or Higgs factory requires significant reduction of the six dimensional emittance of the beam prior to acceleration. One method to accomplish this involves building a cooling channel using high pressure gas filled radio frequency cavities. The performance of such a cavity when subjected to an intense particle beam must be investigated before this technology can be validated. To this end, a high pressure gas filled radio frequency (rf) test cell was built and placed in a 400 MeV beam line from the Fermilab linac to study the plasma evolution and its effect on the cavity. Hydrogen, deuterium, helium and nitrogen gases were studied. Additionally, sulfur hexafluoride and dry air were used as dopants to aid in the removal of plasma electrons. Measurements were made using a variety of beam intensities, gas pressures, dopant concentrations, and cavity rf electric fields, both with and without a 3 T external solenoidal magnetic field. Energy dissipation per electron-ion pair, electron-ion recombination rates, ion-ion recombination rates, and electron attachment times to $SF_6$ and $O_2$ were measured.


title: Structural Logic and Abstract Elementary Classes with Intersection

abstract: We give a syntactic characterization of abstract elementary classes (AECs) closed under intersections using a new logic with a quantifier for isomorphism types that we call structural logic: we prove that AECs with intersections correspond to classes of models of a universal theory in structural logic. This generalizes Tarski's syntactic characterization of universal classes. As a corollary, we obtain that any AEC with countable L\"owenheim-Skolem number is axiomatizable in $\mathbb{L}_{\infty, \omega} (Q)$, where $Q$ is the quantifier "there exists uncountably many".


title: Particles co-orbital to Janus and Epimetheus: a firefly planetary ring

abstract: The Cassini spacecraft found a new and unique ring that shares the trajectory of Janus and Epimetheus, co-orbital satellites of Saturn. Performing image analysis, we found this to be a continuous ring. Its width is between 30% and 50% larger than previously announced. We also verified that the ring behaves like a firefly. It can only be seen from time to time, when Cassini, the ring and the Sun are arranged in a particular geometric configuration, in very high phase angles. Otherwise, it remains "in the dark", not visible to Cassini's cameras. Through numerical simulations, we found a very short lifetime for the ring particles, less than a couple of decades. Consequently, the ring needs to be constantly replenished. Using a model of particles production due to micrometeorites impacts on the surfaces of Janus and Epimetheus, we reproduce the ring, explaining its existence and the "firefly" behavior.


title: Quantum to Classical One Way Function and Its Applications in Quantum Money Authentication

abstract: In 2013, Farid and Vasiliev [arXiv:quant-ph/1310.4922] for the first time proposed a way to construct a protocol for the realisation of "{\em Classical to Quantum}" one-way hash function, a derivative of the Quantum one-way function as defined by Gottesman and Chuang [Technical Report arXiv:quant-ph/0105032] and used it for constructing quantum digital signatures. We, on the other hand, for the first time, propose the idea of a different kind of one-way function, which is "{\em quantum-classical}" in nature, that is, it takes an $n$-qubit quantum state of a definite kind as its input and produces a classical output. We formally define such a one-way function and propose a way to construct and realise it. The proposed one-way function turns out to be very useful in authenticating a quantum state in any quantum money scheme and so we can construct many different quantum money schemes based on such a one-way function. Later in the paper, we also give explicit constructions of some interesting quantum money schemes like quantum bitcoins and quantum currency schemes, solely based on the proposed one-way function. The security of such schemes can be explained on the basis of the security of the underlying one-way functions.


title: Descriptive Cellular Homology

abstract: This article introduces descriptive cellular homology on cell complexes, which is an extension of J.H.C. Whitehead's CW topology. A main result is that a descriptive cellular complex is a topology on fibres in a fibre bundle. An application of two forms of cellular homology is given in terms of the persistence of shapes in CW spaces.


title: Curvature and topology dependency of the cosmological spectra

abstract: In this article we investigate dependency of the adiabatic and entropy spectral indices of the cosmological perturbations on the geometry and topology of the background universe. Our discussion includes the post-inflationary universe i.e. radiation-dust mixture era. For this purpose, we first extract an explicit equation describing evolution of the comoving curvature perturbation in the FLRW universe with arbitrary spatial sectional curvature. We may percept when $K\neq 0$, curvature scale would be as significant as the perturbations scales to recognize the behavior of the spectral indices. We also focus on the entropy perturbation in order to extract behavior of the isocurvature spectral index in terms of the curvature index and time. Our analysis shows that spectra of curvature and entropy perturbations in sub-horizon scales could be function of topology. Moreover, an accurate analysis makes clear that time-average of isocurvature index in case $K=0$ is about zero,so that imprint of entropy perturbation in time duration may be negligible. We also consider evolution of the cosmological indices for super-curvature modes in the case $K=-1$. In the most results dependency to curvature, initial conditions and scale modes are thoroughly vivid.


title: Multiplex core-periphery organization of the human connectome

abstract: The behavior of many complex systems is determined by a core of densely interconnected units. While many methods are available to identify the core of a network when connections between nodes are all of the same type, a principled approach to define the core when multiple types of connectivity are allowed is still lacking. Here we introduce a general framework to define and extract the core-periphery structure of multi-layer networks by explicitly taking into account the connectivity of the nodes at each layer. We show how our method works on synthetic networks with different size, density, and overlap between the cores at the different layers. We then apply the method to multiplex brain networks whose layers encode information both on the anatomical and the functional connectivity among regions of the human cortex. Results confirm the presence of the main known hubs, but also suggest the existence of novel brain core regions that have been discarded by previous analysis which focused exclusively on the structural layer. Our work is a step forward in the identification of the core of the human connectome, and contributes to shed light to a fundamental question in modern neuroscience.


title: The gravitational wave signal from core-collapse supernovae

abstract: We study gravitational waves (GWs) from a set of two-dimensional multi-group neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Our goal is to systematize the current knowledge about the post-bounce CCSN GW signal and recognize the templatable features that could be used by the ground-based laser interferometers. We demonstrate that starting from ~400ms after core bounce the dominant GW signal represents the fundamental quadrupole (l=2) oscillation mode (f-mode) of the proto-neutron star (PNS), which can be accurately reproduced by a linear perturbation analysis of the angle-averaged PNS profile. Before that, in the time interval between ~200 and ~400ms after bounce, the dominant mode has two radial nodes and represents a g-mode. We associate the high-frequency noise in the GW spectrograms above the main signal with p-modes, while below the dominant frequency there is a region with very little power. The collection of models presented here summarizes the dependence of the CCSN GW signal on the progenitor mass, equation of state, many-body corrections to the neutrino opacity, and rotation. Weak dependence of the dominant GW frequency on the progenitor mass motivates us to provide a simple fit for it as a function of time, which can be used as a prior when looking for CCSN candidates in the LIGO data.


title: An application of the partial r-Bell polynomials on some family of bivariate polynomials

abstract: The aim of this paper is to give some combinatorial relations linked polynomials generalizing those of Appell type to the partial r-Bell polynomials. We give an inverse relation, recurrence relations involving some family of polynomials and their exact expressions at rational values in terms of the partial r-Bell polynomials. We illustrate the obtained results by various comprehensive examples.


title: Fitting a round peg into a round hole: asympotically correcting the generalized gradient approximation for correlation

abstract: We revisit the two derivations of the PBE correlation functional: The real-space cut-off of the exchange-correlation hole and the imposition of exact conditions. These differ in the Lieb-Simon limit, exemplified by the scaling of neutral atoms to large $N$ and $Z$, in which LDA becomes relatively exact. We use the leading correction to this limit for neutral atoms to design an asymptotically corrected correlation GGA as a compromise between these two constructions which becomes relatively more accurate for atoms with increasing atomic number. When paired with a similar correction for exchange, this acGGA satisfies more exact conditions than PBE. Combined with the known $r_s$-dependence of the gradient expansion for correlation, this correction accurately reproduces correlation energies of closed shell atoms down to Be. We test this acGGA for atoms and molecules, finding substantial improvements over PBE, but also showing that optimal global hybrids of acGGA do not improve upon PBE0, and are similar to meta-GGA values. We discuss the relevance of these results to Jacob's ladder of non-empirical density functional construction.


title: Multiplicative operators in the spectral problem of integrable systems

abstract: We consider the spectral problem of the Lax pair associated to periodic integrable partial differential equations. We assume this spectral problem to be a polynomial of degree $d$ in the spectral parameter $\lambda$. From this assumption, we find the conservation laws as well as the hyperelliptic curve required to solve the periodic inverse problem. A recursion formula is developed, as well as $d$ additional conditions which give additional information to integrate the equations under consideration. We also include two examples to show how the techniques developed work. For the Korteweg-deVries (KdV) equation, the degree of the multiplicative equation is $d=1$. Hence, we only have one condition and one recursion formula. The condition gives in each degree of the recursion the conserved densities for KdV equation, recovering the Lax hierarchy. For the Nonlinear Schr\"odinger (NLS) equation, the degree of the multiplicative operator is $d=2$. Hence, we have a couple of conditions that we use to deduce conserved density constants and the Lax hierarchy for such equation. Additionally, we explicitely write down the hyper-elliptic curve associted to the NLS equation. Our approach can be use for other completely integrable differential equations, as long as we have a polynomial multiplicative operator associated to them.


title: Homogenization of some degenerate pseudoparabolic variational inequalities

abstract: Multiscale analysis of a degenerate pseudoparabolic variational inequality, modelling the two-phase flow with dynamical capillary pressure in a perforated domain, is the main topic of this work. Regularisation and penalty operator methods are applied to show the existence of a solution of the nonlinear degenerate pseudoparabolic variational inequality defined in a domain with microscopic perforations, as well as to derive a priori estimates for solutions of the microscopic problem. The main challenge is the derivation of a priori estimates for solutions of the variational inequality, uniformly with respect to the regularisation parameter and to the small parameter defining the scale of the microstructure. The method of two-scale convergence is used to derive the corresponding macroscopic obstacle problem.


title: The global and radial stellar mass assembly of Milky Way-sized galaxies

abstract: We study the global and radial stellar mass assembly of eight zoomed-in MW-sized galaxies produced in Hydrodynamics cosmological simulations. The disk-dominated galaxies (4) show a fast initial stellar mass growth in the innermost parts, driven mostly by in-situ SF, but since $z\sim2-1$ the SF enters in a long-term quenching phase. The outer regions follow this trend but more gentle as more external they are. As the result, the radial stellar mass growth is highly inside-out due to both the inside-out structural growth and inside-out SF quenching. The half-mass radius evolves fast; for instance, $R_{0.5}$($z=1$)$<0.5$$R_{0.5}$($z=0$). Two other runs resemble lenticular galaxies. One shows also a pronounced inside-out growth and the other one presents a nearly uniform radial mass assembly. The other two galaxies suffered late major mergers. Their normalized radial mass growth histories (MGHs) are nearly close among them but with periods of outside-in assembly during or after the mergers. For all the simulations, the archaeological radial MGHs calculated from the $z = 0$ stellar-particles age distribution are similar to the current MGHs, which evidences that the mass assembly by ex-situ stars and the radial mass transport do not change significantly their radial mass distributions. Our results agree qualitatively with observational inferences from the fossil record method applied to a survey of local galaxies and from look-back observations of progenitors of MW-sized galaxies. However, the inside-out growth mode is more pronounced and the $R_{0.5}$ growth is faster in simulations than in observational inferences.


title: Explosions of Thorne-Zytkow objects

abstract: We propose that massive Thorne-Zytkow objects can explode. A Thorne-Zytkow object is a theoretically predicted star that has a neutron core. When nuclear reactions supporting a massive Thorne-Zytkow object terminate, a strong accretion occurs towards the central neutron core. The accretion rate is large enough to sustain a super-Eddington accretion towards the neutron core. The neutron core may collapse to a black hole after a while. A strong large-scale outflow or a jet can be launched from the super-Eddington accretion disk and the collapsing Thorne-Zytkow object can be turned into an explosion. The ejecta have about 10 Msun but the explosion energy depends on when the accretion is suppressed. We presume that the explosion energy could be as low as ~ 1e47 erg and such a low-energy explosion could be observed like a failed supernova. The maximum possible explosion energy is ~ 1e52 erg and such a high-energy explosion could be observed as an energetic Type II supernova or a superluminous supernova. Explosions of Thorne-Zytkow objects may provide a new path to spread lithium and other heavy elements produced through the irp process such as molybdenum in the Universe.


title: Urban Explorations: Analysis of Public Park Usage using Mobile GPS Data

abstract: This study analyzes mobile phone data derived from 10 million daily active users across the United States to better understand the spatio-temporal activity patterns of users in Central Park, New York. The aim of this initial investigation is to create quantifiable measures for understanding public space usage in regions of the city that have no natural data source for measuring activity. We analyze the trip behaviors of users across time and different regions in the park to find patterns of co-location and shared time and, thus, potential social interaction. We find that regions with established amenities and points of interest exhibit a higher percentage of shared experiences, indicating that institutional amenities act as 'beacons' for users' experiences in the park.


title: Vectorization of Line Drawings via PolyVector Fields

abstract: Image tracing is a foundational component of the workflow in graphic design, engineering, and computer animation, linking hand-drawn concept images to collections of smooth curves needed for geometry processing and editing. Even for clean line drawings, modern algorithms often fail to faithfully vectorize junctions, or points at which curves meet; this produces vector drawings with incorrect connectivity. This subtle issue undermines the practical application of vectorization tools and accounts for hesitance among artists and engineers to use automatic vectorization software. To address this issue, we propose a novel image vectorization method based on state-of-the-art mathematical algorithms for frame field processing. Our algorithm is tailored specifically to disambiguate junctions without sacrificing quality.


title: Constraining the speed of sound inside neutron stars with chiral effective field theory interactions and observations

abstract: The dense matter equation of state (EOS) determines neutron star (NS) structure but can be calculated reliably only up to one to two times the nuclear saturation density, using accurate many-body methods that employ nuclear interactions from chiral effective field theory constrained by scattering data. In this work, we use physically motivated ansatzes for the speed of sound $c_S$ at high density to extend microscopic calculations of neutron-rich matter to the highest densities encountered in stable NS cores. We show how existing and expected astrophysical constraints on NS masses and radii from X-ray observations can constrain the speed of sound in the NS core. We confirm earlier expectations that $c_S$ is likely to violate the conformal limit of $c_S^2\leq c^2/3 $, possibly reaching values closer to the speed of light $c$ at a few times the nuclear saturation density, independent of the nuclear Hamiltonian. If QCD obeys the conformal limit, we conclude that the rapid increase of $c_S$ required to accommodate a $2 $ M$_\odot$ NS suggests a form of strongly interacting matter where a description in terms of nucleons will be unwieldy, even between one and two times the nuclear saturation density. For typical NSs with masses in the range $1.2-1.4~$ M$_\odot$, we find radii between $10$ and $14$ km, and the smallest possible radius of a $1.4$ M$_{\odot}$ NS consistent with constraints from nuclear physics and observations is $8.4$ km. We also discuss how future observations could constrain the EOS and guide theoretical developments in nuclear physics.


title: Green matrix estimates of block Jacobi matrices I: Unbounded gap in the essential spectrum

abstract: This work deals with decay bounds for Green matrices and generalized eigenvectors of block Jacobi matrices when the real part of the spectral parameter lies in an infinite gap of the operator's essential spectrum. We consider the cases of commutative and noncommutative matrix entries separately. An example of a block Jacobi operator with noncommutative entries and nonnegative essential spectrum is given to illustrate the results.


title: The reciprocal sum of primitive nondeficient numbers

abstract: We investigate the reciprocal sum of primitive nondeficient numbers, or pnds. In 1934, Erdos showed that the reciprocal sum of pnds converges, which he used to prove that the abundant numbers have a natural density. We show the reciprocal sum of pnds is between 0.348 and 0.380.


title: Logarithmic Corrections to Black Hole Entropy: the Non-BPS Branch

abstract: We compute the leading logarithmic correction to black hole entropy on the non-BPS branch of 4D ${\cal N}\geq 2$ supergravity theories. This branch corresponds to finite temperature black holes whose extremal limit does not preserve supersymmetry, such as the $D0-D6$ system in string theory. Starting from a black hole in minimal Kaluza-Klein theory, we discuss in detail its embedding into ${\cal N}=8, 6, 4, 2$ supergravity, its spectrum of quadratic fluctuations in all these environments, and the resulting quantum corrections. We find that the $c$-anomaly vanishes only when ${\cal N}\geq 6$, in contrast to the BPS branch where $c$ vanishes for all ${\cal N}\geq 2$. We briefly discuss potential repercussions this feature could have in a microscopic description of these black holes.


title: Thermal stitching: Extending the reach of quantum fermion solvers

abstract: For quantum fermion problems, many accurate solvers are limited by the temperature regime in which they can be usefully applied. The Mermin theorem implies the uniqueness of an effective potential from which both the exact density and free energy at a target temperature can be found, via a calculation at a different, reference temperature. We derive exact expressions for both the potential and the free energy in such a calculation, and introduce three controllable approximations that reduce the cost of such calculations. We illustrate the effective potential and its free energy, and test the approximations, on the asymmetric two-site Hubbard model at finite temperature.


title: Tensor Train decomposition on TensorFlow (T3F)

abstract: Tensor Train decomposition is used across many branches of machine learning. We present T3F -- a library for Tensor Train decomposition based on TensorFlow. T3F supports GPU execution, batch processing, automatic differentiation, and versatile functionality for the Riemannian optimization framework, which takes into account the underlying manifold structure to construct efficient optimization methods. The library makes it easier to implement machine learning papers that rely on the Tensor Train decomposition. T3F includes documentation, examples and 94% test coverage.


title: Role of core-collapse supernovae in explaining Solar System abundances of p nuclides

abstract: The production of the heavy stable proton-rich isotopes between 74Se and 196Hg -- the p nuclides -- is due to the contribution from different nucleosynthesis processes, activated in different types of stars. Whereas these processes have been subject to various studies, their relative contributions to Galactic Chemical Evolution are still a matter of debate. We investigate for the first time the nucleosynthesis of p nuclides in GCE by including metallicity and progenitor mass-dependent yields of core-collapse supernovae (ccSNe) into a chemical evolution model. We used a grid of metallicities and progenitor masses from two different sets of stellar yields and followed the contribution of ccSNe to the Galactic abundances as function of time. In combination with previous studies on p-nucleus production in thermonuclear supernovae (SNIa), and using the same GCE description, this allows us to compare the respective roles of SNeIa and ccSNe in the production of p-nuclei in the Galaxy. The gamma process in ccSN is very efficient for a wide range of progenitor masses (13M_sun_- 25M_sun) at solar metallicity. Since it is a secondary process with its efficiency depending on the initial abundance of heavy elements, its contribution is strongly reduced below solar metallicity. This makes it challenging to explain the inventory of the p nuclides in the Solar System by the contribution from ccSNe alone. In particular, we find that ccSNe contribute less than 10% of the solar p nuclide abundances, with only a few exceptions. Due to the uncertain contribution from other nucleosynthesis sites in ccSNe, such as neutrino winds or $\alpha$-rich freeze out, we conclude that the light p-nuclides 74$Se,78$Kr, 84Sr, and 92Mo may either still be completely or only partially produced in ccSNe.


title: Universal avalanche statistics and triggering close to failure in a mean field model of rheological fracture

abstract: The hypothesis of critical failure relates the presence of an ultimate stability point in the structural constitutive equation of materials to a divergence of characteristic scales in the microscopic dynamics responsible for deformation. Avalanche models involving critical failure have determined common universality classes for stick-slip processes and fracture. However, not all empirical failure processes exhibit the trademarks of criticality. The rheological properties of materials introduce dissipation, usually reproduced in conceptual models as a hardening of the coarse grained elements of the system. Here, we investigate the effects of transient hardening on (i) the activity rate and (ii) the statistical properties of avalanches. We find the explicit representation of transient hardening in the presence of generalized viscoelasticity and solve the corresponding mean field model of fracture. In the quasistatic limit, the accelerated energy release is invariant with respect to rheology and the avalanche propagation can be reinterpreted in terms of a stochastic counting process. A single universality class can be defined from such analogy, and all statistical properties depend only on the distance to criticality. We also prove that inter-event correlations emerge due to the hardening --- even in the quasistatic limit --- that can be interpreted as "aftershocks" and "foreshocks".


title: Magnetic oscillations in a model of graphene

abstract: We consider a quantum graph as a model of graphene in constant magnetic field and describe the density of states in terms of relativistic Landau levels satisfying a Bohr--Sommerfeld quantization condition. That provides semiclassical corrections (with the magnetic flux as the semiclassical parameter) in the study of magnetic oscillations.


title: Tempest: Temporal Dynamics in Anonymity Systems

abstract: Many recent proposals for anonymous communication omit from their security analyses a consideration of the effects of time on important system components. In practice, many components of anonymity systems, such as the client location and network structure, exhibit changes and patterns over time. In this paper, we focus on the effect of such temporal dynamics on the security of anonymity networks. We present Tempest, a suite of novel attacks based on (1) client mobility, (2) usage patterns, and (3) changes in the underlying network routing. Using experimental analysis on real-world datasets, we demonstrate that these temporal attacks degrade user privacy across a wide range of anonymity networks, including deployed systems such as Tor; path-selection protocols for Tor such as DeNASA, TAPS, and Counter-RAPTOR; and network-layer anonymity protocols for Internet routing such as Dovetail and HORNET. The degradation is in some cases surprisingly severe. For example, a single host failure or network route change could quickly and with high certainty identify the client's ISP to a malicious host or ISP. The adversary behind each attack is relatively weak - generally passive and in control of one network location or a small number of hosts. Our findings suggest that designers of anonymity systems should rigorously consider the impact of temporal dynamics when analyzing anonymity.


title: Improved Style Transfer by Respecting Inter-layer Correlations

abstract: A popular series of style transfer methods apply a style to a content image by controlling mean and covariance of values in early layers of a feature stack. This is insufficient for transferring styles that have strong structure across spatial scales like, e.g., textures where dots lie on long curves. This paper demonstrates that controlling inter-layer correlations yields visible improvements in style transfer methods. We achieve this control by computing cross-layer, rather than within-layer, gram matrices. We find that (a) cross-layer gram matrices are sufficient to control within-layer statistics. Inter-layer correlations improves style transfer and texture synthesis. The paper shows numerous examples on "hard" real style transfer problems (e.g. long scale and hierarchical patterns); (b) a fast approximate style transfer method can control cross-layer gram matrices; (c) we demonstrate that multiplicative, rather than additive style and content loss, results in very good style transfer. Multiplicative loss produces a visible emphasis on boundaries, and means that one hyper-parameter can be eliminated.


title: Universality results for kinetically constrained spin models in two dimensions

abstract: Kinetically constrained models (KCM) are reversible interacting particle systems on $\mathbb Z^d$ with continuous time Markov dynamics of Glauber type, which represent a natural stochastic (and non-monotone) counterpart of the family of cellular automata known as $\mathcal U$-bootstrap percolation. KCM also display some of the peculiar features of the so-called "glassy dynamics", and as such they are extensively used in the physics literature to model the liquid-glass transition, a major and longstanding open problem in condensed matter physics. We consider two-dimensional KCM with update rule $\mathcal U$, and focus on proving universality results for the mean infection time of the origin, in the same spirit as those recently established in the setting of $\mathcal U$-bootstrap percolation. We first identify what we believe are the correct universality classes, which turn out to be different from those of $\mathcal U$-bootstrap percolation. We then prove universal upper bounds on the mean infection time within each class, which we conjecture to be sharp up to logarithmic corrections. In certain cases, including all supercritical models, and the well-known Duarte model, our conjecture has recently been confirmed in [MMT]. In fact, in these cases our upper bound is sharp up to a constant factor in the exponent. For certain classes of update rules, it turns out that the infection time of the KCM diverges much faster than for the corresponding $\mathcal U$-bootstrap process when the equilibrium density of infected sites goes to zero. This is due to the occurrence of energy barriers which determine the dominant behaviour for KCM, but which do not matter for the monotone bootstrap dynamics.


title: Davydov-Chaban Hamiltonian within the formalism of deformation-dependent effective mass for Davidson potential

abstract: In this work, we modify the Davydov-Chaban Hamiltonian describing the collective motion of a $\gamma$-rigid atomic nucleus by allowing the mass to depend on nuclear deformation. Exact analytical expressions are derived for energy spectra as well as normalized wave functions for Davidson potential. The model, called Z(4)-DDMD (Deformation Dependent Mass with Davidson potential), is achieved by using the Asymptotic Iteration Method (AIM). The numerical calculations for energy spectra and B(E2) transition probabilities are compared to the experimental data of $^{192-196}$Pt isotopes. The obtained results show an overall agreement with the experiment and an important improvement in respect to other models.


title: Limits on Ionized Gas in M81's Globular Clusters

abstract: We use NSF's Karl G. Jansky Very Large Array to constrain the mass of ionized gas in 206 globular star clusters (GCs) in M81, a nearby spiral galaxy. We detect none of the GCs and impose a typical gas-mass upper limit of 550 solar masses (3-sigma). These findings bear on GC evolution in M81.


title: A Comprehensive Survey of Ontology Summarization: Measures and Methods

abstract: The Semantic Web is becoming a large scale framework that enables data to be published, shared, and reused in the form of ontologies. The ontology which is considered as basic building block of semantic web consists of two layers including data and schema layer. With the current exponential development of ontologies in both data size and complexity of schemas, ontology understanding which is playing an important role in different tasks such as ontology engineering, ontology learning, etc., is becoming more difficult. Ontology summarization as a way to distill knowledge from an ontology and generate an abridge version to facilitate a better understanding is getting more attention recently. There are various approaches available for ontology summarization which are focusing on different measures in order to produce a proper summary for a given ontology. In this paper, we mainly focus on the common metrics which are using for ontology summarization and meet the state-of-the-art in ontology summarization.


title: A note on Arithmetic Diophantine series

abstract: We consider some asymptotic analysis for series related to the work of Hardy and Littlewood on Diophantine approximation, as well as Davenport. In particular, we expand on ideas from some previous work on arithmetic series and the RH.


title: Canonical Schottky barrier heights of the transition metal dichalcogenide monolayers in contact with a metal

abstract: The transition metal dichalcogenide ($MX_{2}$, where $M$=Mo, W and $X$=S, Se, Te) monolayers are of high interest for semiconducting applications at the nanoscale level; this interest is due to both their direct band gaps and high charge mobilities. In this regard, an in-depth understating of the related Schottky barrier heights, associated with the incorporation of $MX_{2}$ sheets into novel low-dimensional metal-semiconductor junctions, is of crucial importance. Herein, we generate and provide analysis of the Schottky barrier heights behavior to account for the metal-induced gap states concept as its explanation. In particular, the present investigations concentrate on the estimation of the charge neutrality levels directly by employing the primary theoretical model, $i.e.$ the cell-averaged Green's function formalism combined the complex band structure technique. The results presented herein place charge neutrality levels in the vicinity of the mid-gap; this is in agreement with previous reports and analogous to the behavior of three-dimensional semiconductors. The calculated canonical Schottky barrier heights are also found to be in agreement with other computational and experimental values in cases where the difference between electronegativities of the semiconductor and metal contact is small. Moreover, the influence of the spin-orbit effects is herein considered and supports that Schottky barrier heights have metal-induced gap state-derived character, regardless whether spin-orbit coupling interactions are considered. The results presented within this report constitute a direct and vital verification of the importance of metal-induced gap states in explaining the behavior of observed Schottky barrier heights at $MX_{2}$-metal junctions.


title: A spiral interface with positive Alt-Caffarelli-Friedman limit at the origin

abstract: We give an example of a pair of nonnegative subharmonic functions with disjoint support for which the Alt-Caffarelli-Friedman monotonicity formula has strictly positive limit at the origin, and yet the interface between their supports lacks a (unique) tangent there. This clarifies a remark appearing in the literature (see \cite{cs05}) that the positivity of the limit of the ACF formula implies unique tangents; this is true under some additional assumptions, but false in general. In our example, blow-ups converge to the expected piecewise linear two-plane function along subsequences, but the limiting function depends on the subsequence due to the spiraling nature of the interface.


title: Dirac electrons in quantum rings

abstract: We consider quantum rings realized in materials where the dynamics of charge carriers mimics that of two-dimensional (2D) Dirac electrons. A general theoretical description of the ring-subband structure is developed that applies to a range of currently available 2D systems, including graphene, transition-metal dichalcogenides, and narrow-gap semiconductor quantum wells. We employ the scattering-matrix approach to calculate the electronic two-terminal conductance through the ring and investigate how it is affected by Dirac-electron interference. The interplay of pseudo-spin chirality and hard-wall confinement is found to distinctly affect the geometric phase that is experimentally accessible in mesoscopic-conductance measurements. We derive an effective Hamiltonian for the azimuthal motion of charge carriers in the ring that yields deeper insight into the physical origin of the observed transport effects, including the unique behavior exhibited by the lowest ring subband in the normal and topological (i.e., band-inverted) regimes. Our work provides a unified approach to characterizing confined Dirac electrons, which can be used to explore the design of valley- and spintronic devices based on quantum interference and the confinement-tunable geometric phase.


title: Vanishing theorems for representation homology and the derived cotangent complex

abstract: Let $G$ be a reductive affine algebraic group defined over a field $k$ of characteristic zero. In this paper, we study the cotangent complex of the derived $G$-representation scheme $ {\rm DRep}_G(X)$ of a pointed connected topological space $X$. We use an (algebraic version of) unstable Adams spectral sequence relating the cotangent homology of $ {\rm DRep}_G(X) $ to the representation homology $ {\rm HR}_*(X,G) := \pi_*{\mathcal O}[{\rm DRep}_G(X)] $ to prove some vanishing theorems for groups and geometrically interesting spaces. Our examples include virtually free groups, Riemann surfaces, link complements in $ {\mathbb R}^3 $ and generalized lens spaces. In particular, for any f.g. virtually free group $ \Gamma $, we show that $\, {\rm HR}_i({\rm B}\Gamma, G) = 0 \,$ for all $ i > 0 $. For a closed Riemann surface $\Sigma_g $ of genus $ g \ge 1 $, we have $\, {\rm HR}_i(\Sigma_g, G) = 0 \,$ for all $ i > \dim G $. The sharp vanishing bounds for $ \Sigma_g $ depend actually on the genus: we conjecture that if $ g = 1 $, then $\, {\rm HR}_i(\Sigma_g, G) = 0 \,$ for $ i > {\rm rank}\,G $, and if $ g \ge 2 $, then $\, {\rm HR}_i(\Sigma_g, G) = 0 \,$ for $ i > \dim\,{\mathcal Z}(G) \,$, where $ {\mathcal Z}(G) $ is the center of $G$. We prove these bounds locally on the smooth locus of the representation scheme $ {\rm Rep}_G[\pi_1(\Sigma_g)]\,$ in the case of complex connected reductive groups. One important consequence of our results is the existence of a well-defined $K$-theoretic virtual fundamental class for $ {\rm DRep}_G(X)$ in the sense of Ciocan-Fontanine and Kapranov. We give a new `Tor formula' for this class in terms of functor homology.


title: Models of bright nickel-free supernovae from stripped massive stars with circumstellar shells

abstract: The nature of an emerging class of rapidly fading supernovae (RFSNe)--characterized by their short-lived light curve duration, but varying widely in peak brightness--remains puzzling. Whether the RFSNe arise from low-mass thermonuclear eruptions on white dwarfs or from the core collapse of massive stars is still a matter of dispute. We explore the possibility that the explosion of hydrogen-free massive stars could produce bright but rapidly fading transients if the effective pre-supernova radii are large and if little or no radioactive nickel is ejected. The source of radiation is then purely due to shock cooling. We study this model of RFSNe using spherically symmetric hydrodynamics and radiation transport calculations of the explosion of stripped stars embedded in helium-dominated winds or shells of various masses and extent. We present a parameter study showing how the properties of the circumstellar envelopes affect the dynamics of the explosion and can lead to a diversity of light curves. We also explore the dynamics of the fallback of the innermost stellar layers, which might be able to remove radioactive nickel from the ejecta, making the rapid decline in the late time light curve possible. We provide scaling relations that describe how the duration and luminosity of these events depend on the supernova kinetic energy and the mass and radius of the circumstellar material.


title: Audio Adversarial Examples: Targeted Attacks on Speech-to-Text

abstract: We construct targeted audio adversarial examples on automatic speech recognition. Given any audio waveform, we can produce another that is over 99.9% similar, but transcribes as any phrase we choose (recognizing up to 50 characters per second of audio). We apply our white-box iterative optimization-based attack to Mozilla's implementation DeepSpeech end-to-end, and show it has a 100% success rate. The feasibility of this attack introduce a new domain to study adversarial examples.


title: Quadrilaterals inscribed in convex curves

abstract: We classify the set of quadrilaterals that can be inscribed in convex Jordan curves, in the continuous as well as in the smooth case. This answers a question of Makeev in the special case of convex curves. The difficulty of this problem comes from the fact that standard topological arguments to prove the existence of solutions do not apply here due to the lack of sufficient symmetry. Instead, the proof makes use of an area argument of Karasev and Tao, which we furthermore simplify and elaborate on. The continuous case requires an additional analysis of the singular points, and a small miracle, which then extends to show that the problems of inscribing isosceles trapezoids in smooth curves and in piecewise $C^1$ curves are equivalent.


title: Ground-state properties of light kaonic nuclei signaling symmetry energy at high densities

abstract: A sensitive correlation between the ground-state properties of light kaonic nuclei and the symmetry energy at high densities is constructed under the framework of relativistic mean-field theory. Taking oxygen isotopes as an example, we see that a high-density core is produced in kaonic oxygen nuclei, due to the strongly attractive antikaon-nucleon interaction. It is found that the $1S_{1/2}$ state energy in the high-density core of kaonic nuclei can directly probe the variation of the symmetry energy at supranormal nuclear density, and a sensitive correlation between the neutron skin thickness and the symmetry energy at supranormal density is established directly. Meanwhile, the sensitivity of the neutron skin thickness to the low-density slope of the symmetry energy is greatly increased in the corresponding kaonic nuclei. These sensitive relationships are established upon the fact that the isovector potential in the central region of kaonic nuclei becomes very sensitive to the variation of the symmetry energy. These findings might provide another perspective to constrain high-density symmetry energy, and await experimental verification in the future.


title: RadVel: The Radial Velocity Modeling Toolkit

abstract: RadVel is an open source Python package for modeling Keplerian orbits in radial velocity (RV) time series. RadVel provides a convenient framework to fit RVs using maximum a posteriori optimization and to compute robust confidence intervals by sampling the posterior probability density via Markov Chain Monte Carlo (MCMC). RadVel allows users to float or fix parameters, impose priors, and perform Bayesian model comparison. We have implemented realtime MCMC convergence tests to ensure adequate sampling of the posterior. RadVel can output a number of publication-quality plots and tables. Users may interface with RadVel through a convenient command-line interface or directly from Python. The code is object-oriented and thus naturally extensible. We encourage contributions from the community. Documentation is available at http://radvel.readthedocs.io.


title: Why Markets are Inefficient: A Gambling "Theory" of Financial Markets For Practitioners and Theorists

abstract: The purpose of this article is to propose a new "theory," the Strategic Analysis of Financial Markets (SAFM) theory, that explains the operation of financial markets using the analytical perspective of an enlightened gambler. The gambler understands that all opportunities for superior performance arise from suboptimal decisions by humans, but understands also that knowledge of human decision making alone is not enough to understand market behavior --- one must still model how those decisions lead to market prices. Thus are there three parts to the model: gambling theory, human decision making, and strategic problem solving. A new theory is necessary because at this writing in 2017, there is no theory of financial markets acceptable to both practitioners and theorists. Theorists' efficient market theory, for example, cannot explain bubbles and crashes nor the exceptional returns of famous investors and speculators such as Warren Buffett and George Soros. At the same time, a new theory must be sufficiently quantitative, explain market "anomalies" and provide predictions in order to satisfy theorists. It is hoped that the SAFM framework will meet these requirements.


title: Face Flashing: a Secure Liveness Detection Protocol based on Light Reflections

abstract: Face authentication systems are becoming increasingly prevalent, especially with the rapid development of Deep Learning technologies. However, human facial information is easy to be captured and reproduced, which makes face authentication systems vulnerable to various attacks. Liveness detection is an important defense technique to prevent such attacks, but existing solutions did not provide clear and strong security guarantees, especially in terms of time. To overcome these limitations, we propose a new liveness detection protocol called Face Flashing that significantly increases the bar for launching successful attacks on face authentication systems. By randomly flashing well-designed pictures on a screen and analyzing the reflected light, our protocol has leveraged physical characteristics of human faces: reflection processing at the speed of light, unique textual features, and uneven 3D shapes. Cooperating with working mechanism of the screen and digital cameras, our protocol is able to detect subtle traces left by an attacking process. To demonstrate the effectiveness of Face Flashing, we implemented a prototype and performed thorough evaluations with large data set collected from real-world scenarios. The results show that our Timing Verification can effectively detect the time gap between legitimate authentications and malicious cases. Our Face Verification can also differentiate 2D plane from 3D objects accurately. The overall accuracy of our liveness detection system is 98.8\%, and its robustness was evaluated in different scenarios. In the worst case, our system's accuracy decreased to a still-high 97.3\%.


title: High Dimensional Elliptical Sliced Inverse Regression in non-Gaussian Distributions

abstract: Sliced inverse regression (SIR) is the most widely-used sufficient dimension reduction method due to its simplicity, generality and computational efficiency. However, when the distribution of the covariates deviates from the multivariate normal distribution, the estimation efficiency of SIR is rather low. In this paper, we propose a robust alternative to SIR - called elliptical sliced inverse regression (ESIR) for analysing high dimensional, elliptically distributed data. There are wide range of applications of the elliptically distributed data, especially in finance and economics where the distribution of the data is often heavy-tailed. To tackle the heavy-tailed elliptically distributed covariates, we novelly utilize the multivariate Kendall's tau matrix in a framework of so-called generalized eigenvector problem for sufficient dimension reduction. Methodologically, we present a practical algorithm for our method. Theoretically, we investigate the asymptotic behavior of the ESIR estimator and obtain the corresponding convergence rate under high dimensional setting. Quantities of simulation results show that ESIR significantly improves the estimation efficiency in heavy-tailed scenarios. A stock exchange data analysis also demonstrates the effectiveness of our method. Moreover, ESIR can be easily extended to most other sufficient dimension reduction methods.


title: Initial Beam Dynamics Simulations of a High-Average-Current Field-Emission Electron Source in a Superconducting RadioFrequency Gun

abstract: High-power electron beams are sought-after tools in support to a wide array of societal applications. This paper investigates the production of high-power electron beams by combining a high-current field-emission electron source to a superconducting radio-frequency (SRF) cavity. We especially carry out beam-dynamics simulations that demonstrate the viability of the scheme to form $\sim$ 300 kW average-power electron beam using a 1+1/2-cell SRF gun.


title: Generating Neural Networks with Neural Networks

abstract: Hypernetworks are neural networks that generate weights for another neural network. We formulate the hypernetwork training objective as a compromise between accuracy and diversity, where the diversity takes into account trivial symmetry transformations of the target network. We explain how this simple formulation generalizes variational inference. We use multi-layered perceptrons to form the mapping from the low dimensional input random vector to the high dimensional weight space, and demonstrate how to reduce the number of parameters in this mapping by parameter sharing. We perform experiments and show that the generated weights are diverse and lie on a non-trivial manifold.


title: Adversarial Perturbation Intensity Achieving Chosen Intra-Technique Transferability Level for Logistic Regression

abstract: Machine Learning models have been shown to be vulnerable to adversarial examples, ie. the manipulation of data by a attacker to defeat a defender's classifier at test time. We present a novel probabilistic definition of adversarial examples in perfect or limited knowledge setting using prior probability distributions on the defender's classifier. Using the asymptotic properties of the logistic regression, we derive a closed-form expression of the intensity of any adversarial perturbation, in order to achieve a given expected misclassification rate. This technique is relevant in a threat model of known model specifications and unknown training data. To our knowledge, this is the first method that allows an attacker to directly choose the probability of attack success. We evaluate our approach on two real-world datasets.


title: All-dielectric metasurface for wavefront control at terahertz frequencies

abstract: Recently, metasurfaces have gained popularity due to their ability to offer a spatially varying phase response, low intrinsic losses and high transmittance. Here, we demonstrate numerically and experimentally a silicon metasurface at THz frequencies that converts a Gaussian beam into a Vortex beam independent of the polarization of the incident beam. The metasurface consists of an array of sub-wavelength silicon cross resonators made of a high refractive index material on substrates such as sapphire and CaF$_2$ that are transparent at IR-THz spectral range. With these substrates, it is possible to create phase elements for a specific spectral range including at the molecular finger printing around 10$~\mu$m as well as at longer THz wavelengths where secondary molecular structures can be revealed. This device offers high transmittance and a phase coverage of 0 to $2\pi$. The transmittance phase is tuned by varying the dimensions of the meta-atoms. To demonstrate wavefront engineering, we used a discretized spiraling phase profile to convert the incident Gaussian beam to vortex beam. To realize this, we divided the metasurface surface into eight angular sectors and chose eight different dimensions for the crosses providing successive phase shifts spaced by $\pi/4$ radians for each of these sectors. Photolithography and reactive ion etching (RIE) were used to fabricate these silicon crosses as the dimensions of these cylinders range up to few hundreds of micrometers. Large 1-cm-diameter optical elements were successfully fabricated and characterised by optical profilometry.


title: Nodeless superconductivity and time-reversal symmetry breaking in the noncentrosymmetric superconductor Re24Ti5

abstract: The noncentrosymmetric superconductor Re$_{24}$Ti$_{5}$, a time-reversal symmetry (TRS) breaking candidate with $T_c = 6$\,K, was studied by means of muon-spin rotation/relaxation ($\mu$SR) and tunnel-diode oscillator (TDO) techniques. At a macroscopic level, its bulk superconductivity was investigated via electrical resistivity, magnetic susceptibility, and heat capacity measurements. The low-temperature penetration depth, superfluid density and electronic heat capacity all evidence an $s$-wave coupling with an enhanced superconducting gap. The spontaneous magnetic fields revealed by zero-field $\mu$SR below $T_c$ indicate a time-reversal symmetry breaking and thus the unconventional nature of superconductivity in Re$_{24}$Ti$_{5}$. The concomitant occurrence of TRS breaking also in the isostructural Re$_6$(Zr,Hf) compounds, hints at its common origin in this superconducting family and that an enhanced spin-orbital coupling does not affect pairing symmetry.


title: Evidence for GeV Cosmic Rays from White Dwarfs in the Local Cosmic Ray Spectra and in the Gamma-ray Emissivity of the Inner Galaxy

abstract: Recent observations found that electrons are accelerated to $\sim$10 GeV and emit synchrotron hard X-rays in two magnetic white dwarfs (WDs), also known as cataclysmic variables (CVs). In nova outbursts of WDs, multi-GeV gamma-rays were detected inferring that protons are accelerated to 100 GeV or higher. In recent optical surveys, the WD density is found to be higher near the Sun than in the Galactic disk by a factor $\sim$2.5. The cosmic rays (CR) produced by local CVs and novae will accumulate in the local bubble for $10^6$ - $10^7$ yrs. On these findings, we search for CRs from historic CVs and novae in the observed CR spectra. We model the CR spectra at the heliopause as sums of Galactic and local components based on observational data as much as possible. The initial Galactic CR electron and proton spectra are deduced from the gamma-ray emissivity, the local electron spectrum from the hard X-ray spectra at the CVs, and the local proton spectrum inferred by gamma-ray spectrum at novae. These spectral shapes are then expressed in a simple set of polynomial functions of CR energy and regressively fitted until the high-energy ($>$100 GeV) CR spectra near Earth and the Voyager-1 spectra at the heliopause are reproduced. We then extend the modeling to nuclear CR spectra and find that one spectral shape fits all local nuclear CRs and the apparent hardening of the nuclear CR spectra is caused by the roll-down of local nuclear spectra around 100 - 200 GeV. All local CR spectra populate in a limited energy band below 100 - 200 GeV and enhance gamma-ray emissivity below $\sim$10 GeV. Such an enhancement is observed in the inner Galaxy, suggesting the CR fluxes from CVs and novae are substantially higher there.


title: From Eliza to XiaoIce: Challenges and Opportunities with Social Chatbots

abstract: Conversational systems have come a long way since their inception in the 1960s. After decades of research and development, we've seen progress from Eliza and Parry in the 60's and 70's, to task-completion systems as in the DARPA Communicator program in the 2000s, to intelligent personal assistants such as Siri in the 2010s, to today's social chatbots like XiaoIce. Social chatbots' appeal lies not only in their ability to respond to users' diverse requests, but also in being able to establish an emotional connection with users. The latter is done by satisfying users' need for communication, affection, as well as social belonging. To further the advancement and adoption of social chatbots, their design must focus on user engagement and take both intellectual quotient (IQ) and emotional quotient (EQ) into account. Users should want to engage with a social chatbot; as such, we define the success metric for social chatbots as conversation-turns per session (CPS). Using XiaoIce as an illustrative example, we discuss key technologies in building social chatbots from core chat to visual awareness to skills. We also show how XiaoIce can dynamically recognize emotion and engage the user throughout long conversations with appropriate interpersonal responses. As we become the first generation of humans ever living with AI, we have a responsibility to design social chatbots to be both useful and empathetic, so they will become ubiquitous and help society as a whole.


title: Measurement of $\Gamma_{ee}(J/\psi)$ with KEDR detector

abstract: The product of the electronic width of the $J/\psi$ meson and the branching fractions of its decay to hadrons and electrons has been measured using the KEDR detector at the VEPP-4M $e^+e^-$ collider. The obtained values are: $\Gamma_{ee}(J/\psi) = 5.550 \pm 0.056 \pm 0.089 \, \text{keV}, $ $\Gamma_{ee}(J/\psi) \cdot \mathcal{B}_\text{hadrons}(J/\psi) = 4.884 \pm 0.048 \pm 0.078 \, \text{keV}, $ $\Gamma_{ee}(J/\psi) \cdot \mathcal{B}_{ee}(J/\psi) = 0.3331 \pm 0.0066 \pm 0.0040 \, \text{keV}.$ The uncertainties shown are statistical and systematic, respectively. Using the result presented and the world-average value of the electronic branching fraction, one obtains the total width of the $J/\psi$ meson: $\Gamma = 92.94 \pm 1.83 \, \text{keV}.$ These results are consistent with the previous experiments.


title: Frame-based Sparse Analysis and Synthesis Signal Representations and Parseval K-SVD

abstract: Frames are the foundation of the linear operators used in the decomposition and reconstruction of signals, such as the discrete Fourier transform, Gabor, wavelets, and curvelet transforms. The emergence of sparse representation models has shifted of the emphasis in frame theory toward sparse l1-minimization problems. In this paper, we apply frame theory to the sparse representation of signals in which a synthesis dictionary is used for a frame and an analysis dictionary is used for a dual frame. We sought to formulate a novel dual frame design in which the sparse vector obtained through the decomposition of any signal is also the sparse solution representing signals based on a reconstruction frame. Our findings demonstrate that this type of dual frame cannot be constructed for over-complete frames, thereby precluding the use of any linear analysis operator in driving the sparse synthesis coefficient for signal representation. Nonetheless, the best approximation to the sparse synthesis solution can be derived from the analysis coefficient using the canonical dual frame. In this study, we developed a novel dictionary learning algorithm (called Parseval K-SVD) to learn a tight-frame dictionary. We then leveraged the analysis and synthesis perspectives of signal representation with frames to derive optimization formulations for problems pertaining to image recovery. Our preliminary, results demonstrate that the images recovered using this approach are correlated to the frame bounds of dictionaries, thereby demonstrating the importance of using different dictionaries for different applications.


title: Edge optical scattering of two-dimensional materials

abstract: Rayleigh scattering has shown powerful abilities to study electron resonances of nanomaterials regardless of the specific shapes. In analogy to Rayleigh scattering, here we demonstrate that edge optical scattering from two-dimensional (2D) materials also has the similar advantage. Our result shows that, in visible spectral range, as long as the lateral size of a 2D sample is larger than 2 {\mu}m, its edge scattering is essentially a knife-edge diffraction, and the intensity distribution of the high-angle scattering in $\mathbf{k}$ space is nearly independent of the lateral size and the shape of the 2D samples. The high-angle edge scattering spectra are purely determined by the intrinsic dielectric properties of the 2D materials. As an example, we experimentally verify this feature in single-layer $MoS_{2}$, in which A and B excitons are clearly detected in the edge scattering spectra, and the scattering images in $\mathbf{k}$ space and real space are consistent with our theoretical model. This study shows that the edge scattering is a highly practical and efficient method for optical studies of various 2D materials as well as thin films with clear edges.


title: Compressive sensing adaptation for polynomial chaos expansions

abstract: Basis adaptation in Homogeneous Chaos spaces rely on a suitable rotation of the underlying Gaussian germ. Several rotations have been proposed in the literature resulting in adaptations with different convergence properties. In this paper we present a new adaptation mechanism that builds on compressive sensing algorithms, resulting in a reduced polynomial chaos approximation with optimal sparsity. The developed adaptation algorithm consists of a two-step optimization procedure that computes the optimal coefficients and the input projection matrix of a low dimensional chaos expansion with respect to an optimally rotated basis. We demonstrate the attractive features of our algorithm through several numerical examples including the application on Large-Eddy Simulation (LES) calculations of turbulent combustion in a HIFiRE scramjet engine.


title: Expansion of Iterated Stratonovich Stochastic Integrals of Multiplicity 2 Based on Double Fourier-Legendre Series Summarized by Pringsheim Method

abstract: The article is devoted to the expansion of iterated Stratonovich stochastic integrals of second multiplicity into the double series of products of standard Gaussian random variables. The proof of expansion is based on the application of double Fourier-Legendre series and double trigonometric Fourier series summarized by Pringsheim method. The results of the article are generalized to the case of an arbitrary complete orthonormal system of functions in the space $L_2([t, T])$ and $\psi_1(\tau), \psi_2(\tau)\in L_2([t, T]).$ The considered expansion can be applied to the numerical integration of Ito stochastic differential equations. Some recent results on the expansion of iterated Stratonovich stochastic integrals of multiplicities 3 to 6 are given.


title: Cluster algebra structures on Poisson nilpotent algebras

abstract: Various coordinate rings of varieties appearing in the theory of Poisson Lie groups and Poisson homogeneous spaces belong to the large, axiomatically defined class of symmetric Poisson nilpotent algebras, e.g. coordinate rings of Schubert cells for symmetrizable Kac--Moody groups. We prove that every Poisson nilpotent algebra satisfying a mild condition on certain scalars is canonically isomorphic to a cluster algebra which coincides with the corresponding upper cluster algebra, without additional localizations by frozen variables. The constructed cluster structure is compatible with the Poisson structure in the sense of Gekhtman, Shapiro and Vainshtein. All Poisson nilpotent algebras are proved to be equivariant Poisson Unique Factorization Domains. Their seeds are constructed from sequences of Poisson-prime elements for chains of Poisson UFDs; mutation matrices are effectively determined from linear systems in terms of the underlying Poisson structure. Uniqueness, existence, mutation, and other properties are established for these sequences of Poisson-prime elements.


title: Pairing-energy coefficients of neutron-rich fragments in spallation reactions

abstract: The ratio of pairing-energy coefficient to temperature ($a_{p}/T$) of neutron-rich fragments produced in spallation reactions has been investigated by adopting an isobaric yield ratio method deduced in the framework of a modified Fisher model. A series of spallation reactions, 0.5$A$ and 1$A$ GeV $^{208}$Pb + $p$, 1$A$ GeV $^{238}$U + $p$, 0.5$A$ GeV $^{136}$Xe + $d$, 0.2$A$, 0.5$A$ and 1$A$ GeV $^{136}$Xe + $p$, and $^{56}$Fe + $p$ with incident energy ranging from 0.3$A$ to 1.5$A$ GeV, has been analysed. An obvious odd-even staggering is shown in the fragments with small neutron excess ($I\equiv N - Z$), and in the relatively small-$A$ fragments which have large $I$. The values of $a_{p}/T$ for the fragments, with $I$ from 0 to 36, have been found to be in a range from -4 to 4, and most values of $a_{p}/T$ fall in the range from -1 to 1. It is suggested that a small pairing-energy coefficient should be considered in predicting the cross sections of fragments in spallation reactions. It is also concluded that the method proposed in this article is not good for fragments with $A/A_{s} >$ 85\% (where $A_{s}$ is the mass number of the spallation system).


title: Restrictions on Potential Automatic Structures on Thompson's Group F

abstract: We show that a large class of languages in the standard finite generating set X = {$x_0, x_1, x_0^{-1}, x_1^{-1}$} cannot be part of an automatic structure for Thompson's Group F. These languages are ones that accept at least one representative of each element of F of word length that is within a fixed constant of a geodesic representative of the element. To accomplish this, we look at a specific element of F and trace two different paths through the Cayley graph to that element. We show that staying within the length restrictions along these two paths would force that element to have contradictory properties.


title: Complete conformal classification of the Friedmann-Lemaitre-Robertson-Walker solutions with a linear equation of state

abstract: We completely classify Friedmann-Lema\^{i}tre-Robertson-Walker solutions with spatial curvature $K=0,\pm 1$ and equation of state $p=w\rho$, according to their conformal structure, singularities and trapping horizons. We do not assume any energy conditions and allow $\rho < 0$, thereby going beyond the usual well-known solutions. For each spatial curvature, there is an initial spacelike big-bang singularity for $w>-1/3$ and $\rho>0$, while no big-bang singularity for $w<-1$ and $\rho>0$. For $K=0$ or $-1$, $-1<w<-1/3$ and $\rho>0$, there is an initial null big-bang singularity. For each spatial curvature, there is a final spacelike future big-rip singularity for $w<-1$ and $\rho>0$, with null geodesics being future complete for $-5/3\le w<-1$ but incomplete for $w<-5/3$. For $w=-1/3$, the expansion speed is constant. For $-1<w<-1/3$ and $K=1$, the universe contracts from infinity, then bounces and expands back to infinity. For $K=0$, the past boundary consists of timelike infinity and a regular null hypersurface for $-5/3<w<-1$, while it consists of past timelike and past null infinities for $w\le -5/3$. For $w<-1$ and $K=1$, the spacetime contracts from an initial spacelike past big-rip singularity, then bounces and blows up at a final spacelike future big-rip singularity. For $w<-1$ and $K=-1$, the past boundary consists of a regular null hypersurface. The trapping horizons are timelike, null and spacelike for $w\in (-1,1/3)$, $w\in \{1/3, -1\}$ and $w\in (-\infty,-1)\cup (1/3,\infty)$, respectively. A negative energy density ($\rho <0$) is possible only for $K=-1$. In this case, for $w>-1/3$, the universe contracts from infinity, then bounces and expands to infinity; for $-1<w<-1/3$, it starts from a big-bang singularity and contracts to a big-crunch singularity; for $w<-1$, it expands from a regular null hypersurface and contracts to another regular null hypersurface.