title: Comparative study of afterpulsing behavior and models in single photon counting avalanche photo diode detectors

abstract: Single-photon detectors, like Avalanche Photo Diodes (APDs), have a great importance in many fields like quantum key distribution, laser ranging, florescence microscopy, etc. Afterpulsing is a typical non ideal behavior of APDs, operated in Geiger mode, which adversely affects any application which measures the number or timing of detection events. Several conflicting studies have tried to model afterpulsing behavior and link it to fundamental semiconductor physics. Here we describe experiments we performed on at least 2 different detectors from 3 manufacturers, all commonly used in quantum optics setups. We report on the timing structure of the counts. By fitting different models to these data, we found statistically significant differences between the detectors used. Furthermore, we report on the presence of high order afterpulses which are not accounted for in any of the standard models of afterpulse behavior. This might have significant implication e.g. in quantum key distribution systems, as the modeling of the devices in use will have to be part of the security analysis for some systems. Significant differences between manufacturers and even between seemingly identical detectors of the same make and batch suggest that every individual detector needs to be calibrated. We describe a reliable procedure to perform such a calibration and present our findings on three different commercial products.


title: Axisymmetric Charge-Conservative Electromagnetic Particle Simulation Algorithm on Unstructured Grids: Application to Microwave Vacuum Electronic Devices

abstract: We present a 2.5-dimensional charge-conservative electromagnetic particle-in-cell (EM-PIC) algorithm optimized for the analysis of vacuum electronic devices (VED) with cylindrical symmetry (axisymmetry). We explore the axisymmetry present in the device geometry, fields, and sources to reduce the dimensionality of the problem from 3D to 2D. Further, we explore `transformation optics' principles to map the original problem in polar coordinates to an equivalent problem on Cartesian coordinates with an effective (artificial) inhomogeneous medium introduced. The resulting problem in the meridian plane is discretized using an unstructured 2D mesh considering TE-polarized fields and properly scaled charges. EM field and source variables (node-based charges and edge-based currents) are expressed as differential forms of various degrees, and discretized using Whitney forms. Using leapfrog time integration, we obtain a mixed finite-element time-domain scheme for the full-discrete Maxwell's equations. We achieve a local and explicit time-update for the field equations by employing the sparse approximate inverse (SPAI) algorithm. Interpolating field values to particles' positions for solving Newton-Lorentz equations of motion is also done via Whitney forms. Particles are advanced using the Boris algorithm with a relativistic correction. In the scatter step, we apply a radial scaling factor on top of a charge-conserving scatter scheme tailored for 2-dimensional unstructured grids. As validation examples, we demonstrate simulations that investigate the physical performance of VEDs designed to harness particle bunching effects arising from the coherent (resonance) Cerenkov electron beam interactions within micromachined slow-wave structures.


title: Vapor Compression Cycle Control for Automotive Air Conditioning Systems with a Linear Parameter Varying Approach

abstract: This paper investigates an output tracking problem for the vapor compression cycle in automotive Air Conditioning (A/C) systems using Linear Parameter Varying (LPV) techniques. Stemming from a recently developed first-principle A/C model, Jacobian linearization is first exploited to develop an LPV-based model that is nonlinearly dependent on time-varying system parameters such as evaporator pressure and superheat temperature. To facilitate the control implementation, a Tensor Product (TP) model transformation is applied to transform the LPV-based model to a TP-type convex polytopic model. LPV controllers are then designed to guarantee system stability, robustness and H-infinity performance. Simulations are presented to demonstrate the efficacy of the developed framework.


title: Rotating NSs/QSs and recent astrophysical observations

abstract: Both fast and slow configurations of rotating neutron stars (NSs) are studied with the recently-constructed unified NS EoSs. The calculations for pure quark stars (QSs) and hybrid stars (HSs) are also done, using several updated quark matter EoSs and Gibbs construction for obtaining hadron-quark mixed phase. All three types of EoSs fulfill the recent 2-solar-mass constrain. By confronting the glitch observations with the theoretical calculations for the crustal moment of inertia (MoI), we find that the glitch crisis is still present in Vela-like pulsars. An upcoming accurate MoI measurement (eg., a possible 10\% accuracy for pulsar PSR J0737-3039A) could distinguish QSs from NSs, since the MoIs of QSs are generally $> \sim 1.5 $ times larger than NSs and HSs, no matter the compactness and the mass of the stars. Using tabulated EoSs, we compute stationary and equilibrium sequences of rapidly rotating, relativistic stars in general relativity from the well-tested $rns$ code, assuming the matter comprising the star to be a perfect fluid. All three observed properties of the short gamma-ray bursts (SGRBs) internal plateaus sample are simulated using the rotating configurations of NSs/QSs as basic input. We finally argue that for some characteristic SGRBs, the post-merger products of NS-NS mergers are probably supramassive QSs rather than NSs, and NS-NS mergers are a plausible location for quark deconfinement and the formation of QSs.


title: A spectral-Galerkin turbulent channel flow solver for large-scale simulations

abstract: A fully (pseudo-)spectral solver for direct numerical simulations of large-scale turbulent channel flows is described. The solver utilizes the Chebyshev base functions suggested by J. Shen [SIAM J. Sci. Comput., 16, 1, 1995], that lead to stable and robust numerical schemes, even at very large scale. New and fast algorithms for the direct solution of the linear systems are devised, and algorithms and matrices for all required scalar products and transforms are provided. We validate the solver for very high Reynolds numbers. Specifically, the solver is shown to reproduce the first order statistics of Hoyas and Jim\'{e}nez [Phys. Fluids, 18(1), 2006], for a channel flow at $Re_{\tau}=2000$. The solver is available through the open source project spectralDNS [https://github.com/spectralDNS].


title: Raman spectral shift versus strain and composition in GeSn layers with: 6 to 15% Sn contents

abstract: GeSn alloys are the subject of intense research activities as these group IV semiconductors present direct bandgap behaviors for high Sn contents. Today, the control of strain becomes an important challenge to improve GeSn devices. Strain micro-measurements are usually performed by Raman spectroscopy. However, different relationships linking the Raman spectral shifts to the built-in strain can be found in the literature. They were deduced from studies on low Sn content GeSn layers (i.e. xSn<8%) or on GeSiSn layers. In this work, we have calibrated the GeSn Raman relationship for really high Sn content GeSn binaries (6<xSn<15%). We have used fully strained GeSn layers and fully relaxed GeSn under-etched microstructures to clearly differentiate the contributions of strain and chemical composition on the Ge-Ge Raman spectral shift. We have shown that the GeSn Raman-strain coefficient for high Sn contents is higher compared to that for pure Ge.


title: CO Spectral Line Energy Distributions in Galactic Sources: Empirical Interpretation of Extragalactic Observations

abstract: The relative populations in rotational transitions of CO can be useful for inferring gas conditions and excitation mechanisms at work in the interstellar medium. We present CO emission lines from rotational transitions observed with Herschel/HIFI in the star-forming cores Orion S, Orion KL, Sgr B2(M), and W49N. Integrated line fluxes from these observations are combined with those from Herschel/PACS observations of the same sources to construct CO spectral line energy distributions (SLEDs) from $5\leq J_u\leq48$. These CO SLEDs are compared to those reported in other galaxies, with the intention of empirically determining which mechanisms dominate excitation in such systems. We find that CO SLEDs in Galactic star-forming cores cannot be used to reproduce those observed in other galaxies, although the discrepancies arise primarily as a result of beam filling factors. The much larger regions sampled by the Herschel beams at distances of several Mpc contain significant amounts of cooler gas which dominate the extragalactic CO SLEDs, in contrast to observations of Galactic star-forming regions which are focused specifically on cores containing primarily hot molecular gas.


title: Dynamics and Morphology of the Milky Way Spiral Arms from the Metallicity Distribution and Radial Mixing

abstract: Albeit radial migration must be a ubiquitous process in disc galaxies, its significance in the evolution of stellar discs is not always reflected through global trends. However, there are other key observables, such as the metallicity distribution function (MDF), that may shed some light in this matter. We argue that the shape of the MDF not only tells us whether the stellar disc experienced radial migration, but it also contains important clues on the structure that triggered it. Specifically, the MDF contains information about the dynamics and morphology of the spiral pattern. To constrain the spiral parameters, we have included a detailed chemical tagging in our simulations; this allows us to produce a restriction of the structural parameters of the spiral arms in the Milky Way as well as a method to constrain chemical evolution models towards the center of the Galactic disc, where no chemical model provides information. We also found that it is unlikely that the Sun was formed near its current galactocentric position, therefore it might be inaccurate to consider the Sun as representative of the chemical abundances in the solar neighborhood. We also show that a stellar disc of the Milky Way, after evolving dynamically and chemically for 5 Gyr, preserves 80% of its original global metallicity gradient despite having suffered important heating and radial migration; this means that the presence of a metallicity gradient in a given galaxy, does not guarantee that radial mixing has not played a role in its evolution.


title: Deep 230-470 MHz VLA Observations of the mini-halo in the Perseus Cluster

abstract: We present a low-frequency view of the Perseus cluster with new observations from the Karl G. Jansky Very Large Array (JVLA) at 230-470 MHz. The data reveal a multitude of new structures associated with the mini-halo. The mini-halo seems to be influenced both by the AGN activity as well as by the sloshing motion of the cool core cluster's gas. In addition, it has a filamentary structure similar to that seen in radio relics found in merging clusters. We present a detailed description of the data reduction and imaging process of the dataset. The depth and resolution of the observations allow us to conduct for the first time a detailed comparison of the mini-halo structure with the X-ray structure as seen in the Chandra X-ray images. The resulting image shows very clearly that the mini-halo emission is mostly contained behind the cold fronts, similar to that predicted by simulations of gas sloshing in galaxy clusters. However, due to the proximity of the Perseus cluster, as well as the quality of the data at low radio frequencies and at X-ray wavelengths, we also find evidence of fine structure. This structure includes several radial radio filaments extending in different directions, a concave radio structure associated with the southern X-ray bay and sharp edges that correlate with X-ray edges. Mini-halos are therefore not simply diffuse, uniform radio sources, but are rather filled with a rich variety of complex structures. These results illustrate the high-quality images that can be obtained with the new JVLA at low radio-frequencies, as well as the necessity to obtain deeper, higher-fidelity radio images of mini-halos and halos in clusters to further understand their origin.


title: Not so lumpy after all: modeling the depletion of dark matter subhalos by Milky Way-like galaxies

abstract: Among the most important goals in cosmology is detecting and quantifying small ($M_{\rm halo}\simeq10^{6-9}~\mathrm{M}_\odot$) dark matter (DM) subhalos. Current probes around the Milky Way (MW) are most sensitive to such substructure within $\sim20$ kpc of the halo center, where the galaxy contributes significantly to the potential. We explore the effects of baryons on subhalo populations in $\Lambda$CDM using cosmological zoom-in baryonic simulations of MW-mass halos from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two halos run using (1) DM-only (DMO), (2) full baryonic physics, and (3) DM with an embedded disk potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain $\sim2\times$ as many subhalos within 100 kpc of the halo center; this excess is $\gtrsim5\times$ within 25 kpc. At $z=0$, the baryonic simulations are completely devoid of subhalos down to $3\times10^6~\mathrm{M}_\odot$ within $15$ kpc of the MW-mass galaxy, and fewer than 20 surviving subhalos have orbital pericenters <20 kpc. Despite the complexities of baryonic physics, the simple addition of an embedded central disk potential to DMO simulations reproduces this subhalo depletion, including trends with radius, remarkably well. Thus, the additional tidal field from the central galaxy is the primary cause of subhalo depletion. Subhalos on radial orbits that pass close to the central galaxy are preferentially destroyed, causing the surviving subhalo population to have tangentially biased orbits compared to DMO predictions. Our method of embedding a disk potential in DMO simulations provides a fast and accurate alternative to full baryonic simulations, thus enabling suites of cosmological simulations that can provide accurate and statistical predictions of substructure populations.


title: Interaction effects in a microscopic quantum wire model with strong spin-orbit interaction

abstract: We investigate the effect of strong interactions on the spectral properties of quantum wires with strong Rashba spin-orbit interaction in a magnetic field, using a combination of Matrix Product State and bosonization techniques. Quantum wires with strong Rashba spin-orbit interaction and magnetic field exhibit a partial gap in one-half of the conducting modes. Such systems have attracted wide-spread experimental and theoretical attention due to their unusual physical properties, among which are spin-dependent transport, or a topological superconducting phase when under the proximity effect of an s-wave superconductor. As a microscopic model for the quantum wire we study an extended Hubbard model with spin-orbit interaction and Zeeman field. We obtain spin resolved spectral densities from the real-time evolution of excitations, and calculate the phase diagram. We find that interactions increase the pseudo gap at $k = 0$ and thus also enhance the Majorana-supporting phase and stabilize the helical spin order. Furthermore, we calculate the optical conductivity and compare it with the low energy spiral Luttinger Liquid result, obtained from field theoretical calculations. With interactions, the optical conductivity is dominated by an excotic excitation of a bound soliton-antisoliton pair known as a breather state. We visualize the oscillating motion of the breather state, which could provide the route to their experimental detection in e.g. cold atom experiments.


title: Studying distant galaxies: A Handbook of Methods and Analyses

abstract: Distant galaxies encapsulate the various stages of galaxy evolution and formation from over 95% of the development of the universe. As early as twenty-five years ago, little was known about them, however since the first systematic survey was completed in the 1990s, increasing amounts of resources have been devoted to their discovery and research. This book summarises for the first time the numerous techniques used for observing, analysing, and understanding the evolution and formation of these distant galaxies. In this rapidly expanding research field, this text is an every-day companion handbook for graduate students and active researchers. It provides guidelines in sample selection, imaging, integrated spectroscopy and 3D spectroscopy, which help to avoid the numerous pitfalls of observational and analysis techniques in use in extragalactic astronomy. It also paves the way for establishing relations between fundamental properties of distant galaxies. At each step, the reader is assisted with numerous practical examples and ready-to-use methodology to help understand and analyse research.


title: The Eccentric Kozai-Lidov mechanism for Outer Test Particle

abstract: The secular approximation of the hierarchical three body systems has been proven to be very useful in addressing many astrophysical systems, from planets, stars to black holes. In such a system two objects are on a tight orbit, and the tertiary is on a much wider orbit. Here we study the dynamics of a system by taking the tertiary mass to zero and solve the hierarchical three body system up to the octupole level of approximation. We find a rich dynamics that the outer orbit undergoes due to gravitational perturbations from the inner binary. The nominal result of the precession of the nodes is mostly limited for the lowest order of approximation, however, when the octupole-level of approximation is introduced the system becomes chaotic, as expected, and the tertiary oscillates below and above 90deg, similarly to the non-test particle flip behavior (e.g., Naoz 2016). We provide the Hamiltonian of the system and investigate the dynamics of the system from the quadrupole to the octupole level of approximations. We also analyze the chaotic and quasi-periodic orbital evolution by studying the surfaces of sections. Furthermore, including general relativity, we show case the long term evolution of individual debris disk particles under the influence of a far away interior eccentric planet. We show that this dynamics can naturally result in retrograde objects and a puffy disk after a long timescale evolution (few Gyr) for initially aligned configuration.


title: Randomized benchmarking of barrier versus tilt control of a singlet-triplet qubit

abstract: Decoherence due to charge noise is one of the central challenges in using spin qubits in semiconductor quantum dots as a platform for quantum information processing. Recently, it has been experimentally demonstrated in both Si and GaAs singlet-triplet qubits that the effects of charge noise can be suppressed if qubit operations are implemented using symmetric barrier control instead of the standard tilt control. Here, we investigate the key issue of whether the benefits of barrier control persist over the entire set of single-qubit gates by performing randomized benchmarking simulations. We find the surprising result that the improvement afforded by barrier control depends sensitively on the amount of spin noise: for the minimal nuclear spin noise levels present in Si, the coherence time improves by more than two orders of magnitude whereas in GaAs, by contrast the coherence time is essentially the same for barrier and tilt control. However, we establish that barrier control becomes beneficial if qubit operations are performed using a new family of composite pulses that reduce gate times by up to 90%. With these optimized pulses, barrier control is the best way to achieve high-fidelity quantum gates in singlet-triplet qubits.


title: Upper limit set by causality on the tidal deformability of a neutron star

abstract: A principal goal of gravitational-wave astronomy is to constrain the neutron star equation of state (EOS) by measuring the tidal deformability of neutron stars. The tidally induced departure of the waveform from that of point-particle (or spinless binary black hole (BBH)) increases with the stiffness of the EOS. We show that causality (the requirement that the speed of sound is less than the speed of light for a perfect fluid satisfying a one-parameter equation of state) places an upper bound on tidal deformability as a function of mass. Like the upper mass limit, the limit on deformabity is obtained by using an EOS with $v_{sound} = c$ for high densities and matching to a low density (candidate) EOS at a matching density of order nuclear saturation density. We use these results and those of [B.D. Lackey et al., Phys. Rev. D 89, 043009 (2014)] to estimate the resulting upper limit on the gravitational-wave phase shift of a black hole-neutron star (BHNS) binary relative to a BBH. Even for assumptions weak enough to allow an maximum mass of $4\ M_\odot$ (a match at nuclear saturation density to an unusually stiff low-density candidate EOS), the upper limit on dimensionless tidal deformability is stringent. It leads to a still more stringent estimated upper limit on the maximum tidally induced phase shift prior to merger. We comment in an appendix on the relation between causality, the condition $v_{sound} < c$, and the condition $dp/d\epsilon < 1$ for the effective EOS governing the equilibrium star.


title: Constraints on the optical polarization source in the luminous non-blazar quasar 3C 323.1 (PG 1545+210) from the photometric and polarimetric variability

abstract: We examine the optical photometric and polarimetric variability of the luminous type 1 non-blazar quasar 3C 323.1 (PG 1545+210). Two optical spectro-polarimetric measurements taken during the periods 1996$-$98 and 2003 combined with a $V$-band imaging polarimetric measurement taken in 2002 reveal that (1) as noted in the literature, the polarization of 3C 323.1 is confined only to the continuum emission, that is, the emission from the broad line region is unpolarized; (2) the polarized flux spectra show evidence of a time-variable broad absorption feature in the wavelength range of the Balmer continuum and other recombination lines; (3) weak variability in the polarization position angle ($PA$) of $\sim$ 4 deg over a time-scale of 4$-$6 years is observed; and (4) the V-band total flux and the polarized flux show highly correlated variability over a time-scale of one year. Taking the above-mentioned photometric and polarimetric variability properties and the results from previous studies into consideration, we propose a geometrical model for the polarization source in 3C 323.1, in which an equatorial absorbing region and an axi-asymmetric equatorial electron-scattering region are assumed to be located between the accretion disc and the broad line region. The scattering/absorbing regions can perhaps be attributed to the accretion disc wind or flared disc surface, but further polarimetric monitoring observations for 3C~323.1 and other quasars with continuum-confined polarization are needed to probe the true physical origins of these regions.


title: Central elements of the Jennings basis and certain Morita invariants

abstract: From Morita theoretic viewpoint, computing Morita invariants is important. We prove that the intersection of the center and the $n$th (right) socle $ZS^n(A) := Z(A) \cap \operatorname{Soc}^n(A)$ of a finite-dimensional algebra $A$ is a Morita invariant; This is a generalization of important Morita invariants --- the center $Z(A)$ and the Reynolds ideal $ZS^1(A)$. As an example, we also studied $ZS^n(FG)$ for the group algebra $FG$ of a finite $p$-group $G$ over a field $F$ of positive characteristic $p$. Such an algebra has a basis along the socle filtration, known as the Jennings basis. We prove certain elements of the Jennings basis are central and hence form a linearly independent set of $ZS^n(FG)$. In fact, such elements form a basis of $ZS^n(FG)$ for every integer $1 \le n \le p$ if $G$ is powerful. As a corollary we have $\operatorname{Soc}^p(FG) \subseteq Z(FG)$ if $G$ is powerful.


title: Primordial black holes as seeds of magnetic fields in the universe

abstract: Although it is assumed that magnetic fields in accretion disks are dragged from the interstellar medium, the idea is likely not applicable to primordial black holes (PBHs) formed in the early universe. Here we show that magnetic fields can be generated in initially unmagnetized accretion disks around PBHs through the Biermann battery mechanism, and therefore provide the small scale seeds of magnetic field in the universe. The radial temperature and vertical density profiles of these disks provide the necessary conditions for the battery to operate naturally. The generated seed fields have a toroidal structure with opposite sign in the upper and lower half of the disk. In the case of a thin accretion disk around a rotating PBH, the field generation rate increases with increasing PBH spin. At a fixed $r/r_{\rm isco}$, where $r$ is the radial distance from the PBH and $r_{\rm isco}$ is the radius of the innermost stable circular orbit, the battery scales as $M^{-9/4}$, where $M$ is the PBH's mass. The very weak dependency of the battery on accretion rate, makes this mechanism a viable candidate to provide seed fields in an initially unmagnetized accretion disk, following which the magnetorotational instability could take over.


title: From Blackbirds to Black Holes: Investigating Capture-Recapture Methods for Time Domain Astronomy

abstract: In time domain astronomy, recurrent transients present a special problem: how to infer total populations from limited observations. Monitoring observations may give a biassed view of the underlying population due to limitations on observing time, visibility and instrumental sensitivity. A similar problem exists in the life sciences, where animal populations (such as migratory birds) or disease prevalence, must be estimated from sparse and incomplete data. The class of methods termed Capture-Recapture is used to reconstruct population estimates from time-series records of encounters with the study population. This paper investigates the performance of Capture-Recapture methods in astronomy via a series of numerical simulations. The Blackbirds code simulates monitoring of populations of transients, in this case accreting binary stars (neutron star or black hole accreting from a stellar companion) under a range of observing strategies. We first generate realistic light-curves for populations of binaries with contrasting orbital period distributions. These models are then randomly sampled at observing cadences typical of existing and planned monitoring surveys. The classical capture-recapture methods, Lincoln-Peterson, Schnabel estimators, related techniques, and newer methods implemented in the Rcapture package are compared. A general exponential model based on the radioactive decay law is introduced, and demonstrated to recover (at 95% confidence) the underlying population abundance and duty cycle, in a fraction of the observing visits (10-50%) required to discover all the sources in the simulation. Capture-Recapture is a promising addition to the toolbox of time domain astronomy, and methods implemented in R by the biostats community can be readily called from within Python.


title: Chemical Abundances of Two Stars in the Large Magellanic Cloud Globular Cluster NGC~1718

abstract: Detailed chemical abundances of two stars in the intermediate-age Large Magellanic Cloud (LMC) globular cluster NGC~1718 are presented, based on high resolution spectroscopic observations with the MIKE spectrograph. The detailed abundances confirm NGC~1718 to be a fairly metal-rich cluster, with an average [Fe/H] ~ -0.55+/-0.01. The two red giants appear to have primordial O, Na, Mg, and Al abundances, with no convincing signs of a composition difference between the two stars---hence, based on these two stars, NGC~1718 shows no evidence for hosting multiple populations. The Mg abundance is lower than Milky Way field stars, but is similar to LMC field stars at the same metallicity. The previous claims of very low [Mg/Fe] in NGC~1718 are therefore not supported in this study. Other abundances (Si, Ca, Ti, V, Mn, Ni, Cu, Rb, Y, Zr, La, and Eu) all follow the LMC field star trend, demonstrating yet again that (for most elements) globular clusters trace the abundances of their host galaxy's field stars. Similar to the field stars, NGC~1718 is found to be mildly deficient in explosive $\alpha$-elements, but moderately to strongly deficient in O, Na, Mg, Al, and Cu, elements which form during hydrostatic burning in massive stars. NGC~1718 is also enhanced in La, suggesting that it was enriched in ejecta from metal-poor AGB stars.


title: Blue Supergiant X-ray Binaries in the Nearby Dwarf Galaxy IC 10

abstract: In young starburst galaxies, the X-ray population is expected to be dominated by the relics of the most massive and short-lived stars, black-hole and neutron-star high mass X-ray binaries (XRBs). In the closest such galaxy, IC 10, we have made a multi-wavelength census of these objects. Employing a novel statistical correlation technique, we have matched our list of 110 X-ray point sources, derived from a decade of Chandra observations, against published photometric data. We report an 8 sigma correlation between the celestial coordinates of the two catalogs, with 42 X-ray sources having an optical counterpart. Applying an optical color-magnitude selection to isolate blue supergiant (SG) stars in IC 10, we find 16 matches. Both cases show a statistically significant overabundance versus the expectation value for chance alignments. The blue objects also exhibit systematically higher fx/fv ratios than other stars in the same magnitude range. Blue SG-XRBs include a major class of progenitors of double-degenerate binaries, hence their numbers are an important factor in modeling the rate of gravitational wave sources. We suggest that the anomalous features of the IC 10 stellar population are explained if the age of the IC 10 starburst is close to the time of the peak of interaction for massive binaries.


title: The nearest ultra diffuse galaxy: UGC2162

abstract: We describe the structural, stellar population and gas properties of the nearest Ultra Diffuse Galaxy (UDG) discovered so far: UGC2162 (z=0.00392; R$_{e,g}$=1.7$(\pm$0.2) kpc; $\mu_g(0)$=24.4$\pm$0.1 mag/arcsec$^2$; g-i=0.33$\pm$0.02). This galaxy, located at a distance of 12.3($\pm$1.7) Mpc, is a member of the M77 group. UGC2162 has a stellar mass of $\sim$2($^{+2}_{-1}$)$\times$10$^7$ M$_\odot$ and is embedded within a cloud of HI gas $\sim$10 times more massive: $\sim$1.9($\pm$0.6)$\times$10$^8$ M$_\odot$. Using the width of its HI line as a dynamical proxy, the enclosed mass within the inner R$\sim$5 kpc is $\sim$4.6($\pm$0.8)$\times$10$^9$ M$_\odot$ (i.e. M/L$\sim$200). The estimated virial mass from the cumulative mass curve is $\sim$8($\pm$2)$\times$10$^{10}$ M$_\odot$. Ultra-deep imaging from the IAC Stripe82 Legacy Project show that the galaxy is irregular and has many star forming knots, with a gas-phase metallicity around one-third of the solar value. Its estimated Star Formation Rate (SFR) is $\sim$0.01 M$_\odot$/yr. This SFR would double the stellar mass of the object in $\sim$2 Gyr. If the object were to stop forming stars at this moment, after a passive evolution, its surface brightness would become extremely faint: $\mu_g(0)$$\sim$27 mag/arcsec$^2$ and its size would remain large R$_{e,g}$$\sim$ 1.8 kpc. Such faintness would make it almost undetectable to most present-day surveys. This suggests that there could be an important population of M$_{\star}$$\sim$10$^7$ M$_\odot$ "dark galaxies" in rich environments (depleted of HI gas) waiting to be discovered by current and future ultra-deep surveys.


title: Engineering negative stress-energy densities with quantum energy teleportation

abstract: We use quantum energy teleportation in the light-matter interaction as an operational means to create quantum field states that violate energy conditions and have negative local stress-energy densities. We show that the protocol is optimal in the sense that it scales in a way that saturates the quantum interest conjecture.


title: Non-naturally reductive Einstein metrics on $\mathrm{SO}(n)$

abstract: In this article, we prove that every compact simple Lie group $SO(n)$ for $n\geq 10$ admits at least $2\left([\frac{n-1}{3}]-2\right)$ non-naturally reductive left-invariant Einstein metrics.


title: A Multi-Planet System Transiting the $V$ = 9 Rapidly Rotating F-Star HD 106315

abstract: We report the discovery of a multi-planet system orbiting HD 106315, a rapidly rotating mid F-type star, using data from the K2 mission. HD 106315 hosts a $2.51\pm0.12\,R_\oplus$ sub-Neptune in a 9.5 day orbit, and a $4.31_{-0.27}^{+0.24}\,R_\oplus$ super-Neptune in a 21 day orbit. The projected rotational velocity of HD 106315 (12.9 km s$^{-1}$) likely precludes precise measurements of the planets' masses, but could enable a measurement of the sky-projected spin-orbit obliquity for the outer planet via Doppler tomography. The eccentricities of both planets were constrained to be consistent with 0, following a global modeling of the system that includes a Gaia distance and dynamical arguments. The HD 106315 system is one of few multi-planet systems hosting a Neptune-sized planet for which orbital obliquity measurements are possible, making it an excellent test-case for formation mechanisms of warm-Neptunian systems. The brightness of the host star also makes HD 106315 c a candidate for future transmission spectroscopic follow-up studies.


title: Limitations of spin-fermion models in studying underdoped cuprates

abstract: A microscopic basis is provided for the spin-fermion model used to describe the physics of the underdoped cuprates. In this way, a spin-fermion coupling is shown to take care of the local no double occupancy constraint, which is ignored in the weakly coupled regime. This ingredient is proved to be essential to capture the physics of strong correlations, however. We elaborate further on our recent proposal for the strong-coupling version of the spin-fermion model that may prove to be the correct starting point to deal with the strong correlation physics displayed by the underdoped cuprates.


title: Multi-band polarimetry of post-asymptotic giant branch stars stars. Optical measurements

abstract: We present new optical broad-band (UBVRI) aperture polarimetric observations of 53 post-asymptotic giant branch (AGB) stars selected to exhibit a large near-infrared excess. 24 out of the 53 stars (45% of our sample) are presented for the first time. A statistical analysis shows four distinctive groups of polarized post-AGB stars: unpolarized or very lowly polarized (degree of polarization or DoP < 1%), lowly polarized (1% < DoP < 4%), moderately polarized (4% < DoP < 8%) and highly polarized (DoP > 8%). 23 out of the 53 (66%) belong to the first group, 10 (19%) to the second, five (9%) to the third and only three (6%) to the last group. Approximately, 34% of our sample was found to be unpolarized objects, which is close to the percentage of round planetary nebulae. On average, the low and moderate groups show a wavelength-dependent polarization that increases towards shorter wavelength, implying an intrinsic origin of the polarization, which signifies a Rayleigh-like scattering spectrum typical for non-symmetrical envelopes composed principally of small dust grains. The moderately polarized stars exhibit higher K-W3 and W1-W3 colour indices compared with the group of lowly polarized stars suggesting a possible relation between DoP and mass-loss rate. Moreover, they are found to be systematically colder (redder in B-V), which may be associated with the condensation process close to these stars that results in higher degree of polarization. We also provide evidence that multiple scattering in optically thin polar outflows is the mechanism that gives high DoP in post-AGB stars with a bipolar or multi-polar envelopes.


title: A simplified version of Persson's multiscale theory for rubber friction due to viscoelastic losses

abstract: We show the full multiscale Persson's theory for rubber friction due to viscoelastic losses can be approximated extremely closely to simpler models, like that suggested by Persson in 1998 and similarly by Popov in his 2010 book (but we don't make any use of the so-called "Method of Dimensionality Reduction"), so it is essentially a single scale model at the so called large wavevector cutoff. The dependence on the entire spectrum of roughness is therefore only confusing, and we confirm this with actual exact calculations and reference to recent Lorenz et al (2015) data. The multiscale aspect is irrelevant with respect to the real critical issue: the choice of the "free parameter" best fit truncation cutoff, which shows the models are mainly fitting equations. In this sense, we provide at least a very simple one.


title: Two Small Transiting Planets and a Possible Third Body Orbiting HD 106315

abstract: The masses, atmospheric makeups, spin-orbit alignments, and system architectures of extrasolar planets can be best studied when the planets orbit bright stars. We report the discovery of three bodies orbiting HD 106315, a bright (V = 8.97 mag) F5 dwarf targeted by our K2 survey for transiting exoplanets. Two small, transiting planets have radii of 2.23 (+0.30/-0.25) R_Earth and 3.95 (+0.42/-0.39) R_Earth and orbital periods of 9.55 d and 21.06 d, respectively. A radial velocity (RV) trend of 0.3 +/- 0.1 m/s/d indicates the likely presence of a third body orbiting HD 106315 with period >160 d and mass >45 M_Earth. Transits of this object would have depths of >0.1% and are definitively ruled out. Though the star has v sin i = 13.2 km/s, it exhibits short-timescale RV variability of just 6.4 m/s, and so is a good target for RV measurements of the mass and density of the inner two planets and the outer object's orbit and mass. Furthermore, the combination of RV noise and moderate v sin i makes HD 106315 a valuable laboratory for studying the spin-orbit alignment of small planets through the Rossiter-McLaughlin effect. Space-based atmospheric characterization of the two transiting planets via transit and eclipse spectroscopy should also be feasible. This discovery demonstrates again the power of K2 to find compelling exoplanets worthy of future study.


title: The Cauchy-Leray integral: counter-examples to the $L^p$-theory

abstract: We prove the optimality of the hypotheses guaranteeing the $L^p$-boundedness for the Cauchy-Leray integral in $\mathbb C^n$, $n\geq 2$, obtained in [LS-4]. Two domains, both elementary in nature, show that the geometric requirement of strong $\mathbb C$-linear convexity, together with regularity of order 2, are both necessary.


title: Quantum and Super-quantum enhancements to two-sender, two-receiver channels

abstract: We study the consequences of 'super-quantum non-local correlations' as represented by the PR-box model of Popescu and Rohrlich, and show PR-boxes can enhance the capacity of noisy interference channels between two senders and two receivers. PR-box correlations violate Bell/CHSH inequalities and are thus stronger -- more non-local -- than quantum mechanics; yet weak enough to respect special relativity in prohibiting faster-than-light communication. Understanding their power will yield insight into the non-locality of quantum mechanics. We exhibit two proof-of-concept channels: first, we show a channel between two sender-receiver pairs where the senders are not allowed to communicate, for which a shared super-quantum bit (a PR-box) allows perfect communication. This feat is not achievable with the best classical (senders share no resources) or quantum entanglement-assisted (senders share entanglement) strategies. Second, we demonstrate a class of channels for which a tunable parameter achieves a double separation of capacities; for some range of \epsilon, the super-quantum assisted strategy does better than the entanglement-assisted strategy, which in turn does better than the classical one.


title: Inflation due to non-minimal coupling of $f(R)$ gravity to a scalar field

abstract: In this work we investigate a inflationary scenario generated by a large scalar field $\phi$ that non-minimally couples to a $f(R)$ modified gravity model. For a Starobinsky's like model, it is found that along a particular flat direction, the scalar potential takes a simple form $V = (M_p^4/4) [V(\phi)/\alpha(\phi)^2]$ where $\alpha(\phi)$ is a non-minimal coupling to Ricci scalar $R$ in the model. The inflation, therefore, is effectively represented as a single field inflaton scenario. For a specific example, such as a scalar potential $V(\phi) = \mu_1\phi^2 + \mu_2\phi^4$, we found that the predictions match nicely in the $1\sigma$ confidence level of \textit{Plank} TT, TE, EE+lowP combination data of Planck 2015 CMB data for $ 0<\mu_3 \leq 100$, where $\mu_3 := |\mu_1|/(\mu_2 M_p^2)$. For example, taking $\mu_3 = 0.01$ the scalar-to-tensor ratio $r=0.0004$ and and spectral index $n_s = 0.96985$ for $N_\ast = 50$ while taking $\mu_3 = 100.0$ produces $r= 0.03$ and $n_s = 0.96359$ for $N_\ast = 60$.


title: Twin lead ballistic conductor based on nanoribbon edge transport

abstract: If a device like a graphene nanoribbon (GNR) has all its four corners attached to electric current leads, the device becomes a quantum junction through which two electrical circuits can interact. We study such system theoretically for stationary currents. The 4-point energy-dependent conductance matrix of the nanostructure and the classical resistors in the circuits are parameters of the model. The two bias voltages in the circuits are the control variables of the studied system while the electrochemical potentials at the device's terminals are non-trivially dependent on the voltages. For the special case of the linear-response regime analytical formulae for the operation of the coupled quantum-classical device are derived and applied. For higher bias voltages numerical solutions are obtained. The effects of non-equilibrium Fermi levels are captured using a recursive algorithm in which self-consistency between the electrochemical potentials and the currents is reached within few iterations. The developed approach allows to study scenarios ranging from independent circuits to strongly coupled ones. For the chosen model of the GNR with highly conductive zigzag edges we determine the regime in which the single device carries two almost independent currents.


title: The distribution of old stars around the Milky Way's central black hole I: Star counts

abstract: (abridged) In this paper we revisit the problem of inferring the innermost structure of the Milky Way's nuclear star cluster via star counts, to clarify whether it displays a core or a cusp around the central black hole. Through image stacking and improved PSF fitting we push the completeness limit about one magnitude deeper than in previous, comparable work. Contrary to previous work, we analyse the stellar density in well-defined magnitude ranges in order to be able to constrain stellar masses and ages. The RC and brighter giant stars display a core-like surface density profile within a projected radius R<0.3 pc of the central black hole, in agreement with previous studies, but show a cusp-like surface density distribution at larger R. The surface density of the fainter stars can be described well by a single power-law at R<2 pc. The cusp-like profile of the faint stars persists even if we take into account the possible contamination of stars in this brightness range by young pre-main sequence stars. The data are inconsistent with a core-profile for the faint stars.Finally, we show that a 3D Nuker law provides a very good description of the cluster structure. We conclude that the observed stellar density at the Galactic Centre, as it can be inferred with current instruments, is consistent with the existence of a stellar cusp around the Milky Way's central black hole, Sgr A*. This cusp is well developed inside the influence radius of about 3 pc of Sgr A* and can be described by a single three-dimensional power-law with an exponent gamma=1.23+-0.05. The apparent lack of RC stars and brighter giants at projected distances of R < 0.3 pc (R<8") of the massive black hole may indicate that some mechanism has altered their distribution or intrinsic luminosity.


title: The distribution of stars around the Milky Way's central black hole II: Diffuse light from sub-giants and dwarfs

abstract: This is the second of three papers that search for the predicted stellar cusp around the Milky Way's central black hole, Sagittarius A*, with new data and methods. We aim to infer the distribution of the faintest stellar population currently accessible through observations around Sagittarius A*. We use adaptive optics assisted high angular resolution images obtained with the NACO instrument at the ESO VLT. Through optimised PSF fitting we remove the light from all detected stars above a given magnitude limit. Subsequently we analyse the remaining, diffuse light density. The analysed diffuse light arises from sub-giant and main-sequence stars with KS ~ 19 - 20 with masses of 1 - 2 Msol . These stars can be old enough to be dynamically relaxed. The observed power-law profile and its slope are consistent with the existence of a relaxed stellar cusp around the Milky Way's central black hole. We find that a Nuker law provides an adequate description of the nuclear cluster's intrinsic shape (assuming spherical symmetry). The 3D power-law slope near Sgr A* is \gamma = 1.23 +- 0.05. At a distance of 0.01 pc from the black hole, we estimate a stellar mass density of 2.3 +- 0.3 x 10^7 Msol pc^-3 and a total enclosed stellar mass of 180 +- 20 Msol. These estimates assume a constant mass-to-light ratio and do not take stellar remnants into account. The fact that no cusp is observed for bright (Ks 16) giant stars at projected distances of roughly 0.1-0.3 pc implies that some mechanism has altered their appearance or distribution.


title: The distribution of stars around the Milky Way's black hole III: Comparison with simulations

abstract: The distribution of stars around a massive black hole (MBH) has been addressed in stellar dynamics for the last four decades by a number of authors. Because of its proximity, the centre of the Milky Way is the only observational test case where the stellar distribution can be accurately tested. Past observational work indicated that the brightest giants in the Galactic Centre (GC) may show a density deficit around the central black hole, not a cusp-like distribution, while we theoretically expect the presence of a stellar cusp. We here present a solution to this long-standing problem. We performed direct-summation $N-$body simulations of star clusters around massive black holes and compared the results of our simulations with new observational data of the GC's nuclear cluster. We find that after a Hubble time, the distribution of bright stars as well as the diffuse light follow power-law distributions in projection with slopes of $\Gamma \approx 0.3$ in our simulations. This is in excellent agreement with what is seen in star counts and in the distribution of the diffuse stellar light extracted from adaptive-optics (AO) assisted near-infrared observations of the GC. Our simulations also confirm that there exists a missing giant star population within a projected radius of a few arcsec around Sgr A*. Such a depletion of giant stars in the innermost 0.1 pc could be explained by a previously present gaseous disc and collisions, which means that a stellar cusp would also be present at the innermost radii, but in the form of degenerate compact cores.


title: HI vs. H$\alpha$ - Comparing the Kinematic Tracers in Modeling the Initial Conditions of the Mice

abstract: We explore the effect of using different kinematic tracers (HI and H$\alpha$) on reconstructing the encounter parameters of the Mice major galaxy merger (NGC 4676A/B). We observed the Mice using the SparsePak Integral Field Unit (IFU) on the WIYN telescope, and compared the H$\alpha$ velocity map with VLA HI observations. The relatively high spectral resolution of our data (R $\approx$ 5000) allows us to resolve more than one kinematic component in the emission lines of some fibers. We separate the H$\alpha$-[N II] emission of the star-forming regions from shocks using their [N II]/H$\alpha$ line ratio and velocity dispersion. We show that the velocity of star-forming regions agree with that of the cold gas (HI), particularly, in the tidal tails of the system. We reconstruct the morphology and kinematics of these tidal tails utilizing an automated modeling method based on the Identikit software package. We quantify the goodness of fit and the uncertainties of the derived encounter parameters. Most of the initial conditions reconstructed using H$\alpha$ and HI are consistent with each other, and qualitatively agree with the results of previous works. For example, we find 210$\pm^{50}_{40}$ Myrs, and 180$\pm^{50}_{40}$ Myrs for the time since pericenter, when modeling H$\alpha$ and HI kinematics, respectively. This confirms that in some cases, H$\alpha$ kinematics can be used instead of HI kinematics for reconstructing the initial conditions of galaxy mergers, and our automated modeling method is applicable to some merging systems.


title: DARHT Axis-I Diode Simulations II: Geometric Scaling

abstract: The dual-axis radiography for hydrodynamic testing (DARHT) facility at Los Alamos National Laboratory uses two electron linear-induction accelerators (LIA) to produce the source spots for perpendicular flash radiographs of a dynamic experiment.Manipulating the beam current is a means for adjusting the dose, and one way to do this is to change the size of the cathode. This article describes simulations undertaken to develop scaling laws for use as design tools in changing the Axis-1 beam current by changing the cathode size.


title: Gradient-free two-points optimal method for non smooth stochastic convex optimization problem with additional small noise

abstract: Using double-smoothing technique and stochastic mirror descent with inexact oracle we built an optimal algorithm (up to a multiplicative factor) for two-points gradient-free non-smooth stochastic convex programming. We investigate how much can be the level of noise (the nature of this noise isn't necessary stochastic) for the rate of convergence to be maintained (up to a multiplicative factor).


title: Efficient estimation in the Topp-Leone distribution

abstract: In the current paper, the estimation of the probability density function and the cumulative distribution function of the Topp-Leone distribution is considered. We derive the following estimators: maximum likelihood estimator, uniformly minimum variance unbiased estimator, percentile estimator, least squares estimator and weighted least squares estimator. A simulation study shows that the maximum likelihood estimator is more efficient than the others estimators.


title: Strong pseudoconvexity in Banach spaces

abstract: Having been unclear how to define that a domain is strictly pseudoconvex in the infinite-dimensional setting, we develop a general theory having Banach spaces in mind. We first focus on finite dimension and eliminate the need of two degrees of differentiability of the boundary of a domain, since differentiable functions are difficult to find in infinite dimension. We introduce $\ell$-strict pseudoconvexity for $\ell\geq 1$, $1$-strict pseudoconvexity at the boundary, $\ell$-uniform pseudoconvexity for $\ell\geq 0$ and finally strong pseudoconvexity. Defining $\ell$-strict pseudoconvexity and $\ell$-uniform pseudoconvexity for $\ell<2$ depends on extending a notion of strict plurisubharmonicity to cases lacking $C^2$-smoothness, first studying it in the sense of distribution and then considering it in infinite dimension. Examples of strictly plurisubharmonic functions as well as strongly pseudoconvex domains are presented, which end up related to important classical Banach spaces. Finally, some solutions to the inhomogeneous Cauchy-Riemann equations for $\overline{\partial}$-closed $(0,1)$-forms in infinite-dimensional domains are shown, giving new information about domains affinely isomorphic to the ball of $\ell_1$ which appeared in the study of strong pseudoconvexity and some more domains biholomorphically related to open and convex domains of $\ell_1$ such as its ball.


title: Emittance Growth in Linear Induction Accelerators

abstract: The Dual-Axis Radiographic Hydrotest (DARHT) facility uses bremsstrahlung radiation source spots produced by the focused electron beams from two linear induction accelerators (LIAs) to radiograph large hydrodynamic experiments driven by high explosives. Radiographic resolution is determined by the size of the source spot, and beam emittance is the ultimate limitation to spot size. Some of the possible causes for the emittance growth in the DARHT LIA have been investigated using particle-in-cell (PIC) codes, and are discussed in this article. The results suggest that the most likely source of emittance growth is a mismatch of the beam to the magnetic transport, which can cause beam halo.


title: Coulomb Branch for A-type Balanced Quivers in 3d $\mathcal{N}=4$ gauge theories

abstract: We study Coulomb branch moduli spaces of a class of three dimensional $\mathcal{N}=4$ gauge theories whose quiver satisfies the balance condition. The Coulomb branch is described by dressed monopole operators which can be counted using the Monopole formula. We mainly focus on A-type quivers in this paper, using Hilbert Series to study their moduli spaces, and present the interesting pattern which emerges. All of these balanced A-type quiver gauge theories can be realized on brane intervals in Type IIB string theory, where mirror symmetry acts by exchanging the five branes and induces an equivalence between Coulomb branch and Higgs branch of mirror pairs. For each theory, we explicitly discuss the gauge invariant generators on the Higgs branch and the relations they satisfy. Finally, some analysis on $D_4$ balanced quivers also presents an interesting structure of their moduli spaces.


title: Streaming k-Means Clustering with Fast Queries

abstract: We present methods for k-means clustering on a stream with a focus on providing fast responses to clustering queries. Compared to the current state-of-the-art, our methods provide substantial improvement in the query time for cluster centers while retaining the desirable properties of provably small approximation error and low space usage. Our algorithms rely on a novel idea of "coreset caching" that systematically reuses coresets (summaries of data) computed for recent queries in answering the current clustering query. We present both theoretical analysis and detailed experiments demonstrating their correctness and efficiency


title: The g-Good-Neighbor Conditional Diagnosability of Locally Twisted Cubes

abstract: In the work of Peng et al. in 2012, a new measure was proposed for fault diagnosis of systems: namely, g-good-neighbor conditional diagnosability, which requires that any fault-free vertex has at least g fault-free neighbors in the system. In this paper, we establish the g-good-neighbor conditional diagnosability of locally twisted cubes under the PMC model and the MM^* model.


title: Surface Plasmon Assisted Gentle Ablation of Nanostructures by Femtosecond Oscillator

abstract: We experimentally demonstrate the use of subwavelength optical nanoantennae to assist the gentle ablation of nanostructures directly using ultralow fluence from a Ti: sapphire oscillator through the excitation of surface plasmon waves. We show that this ablation mechanism is the same for metal and dielectric. The analytical solutions of ablation threshold are in excellent agreement with the experiment estimations. Surface plasmon assisted locally enhanced ablation at nanoscale provides a method for nanomachining, manipulation and modification the nanostructures without collateral thermal damage to the materials. It is also shown that this ablation can deposit low-density high quality thin nano film.


title: Emittance growth in the DARHT Axis-II Downstream Transport

abstract: Using a particle-in-cell (PIC) code, we investigated the possibilities for emittance growth through the quadrupole magnets of the system used to transport the high-current electron beam from an induction accelerator to the bremsstrahlung converter target used for flash radiography. We found that even highly mismatched beams exhibited little emittance growth (< 6%), which we attribute to softening of their initial hard edge current distributions. We also used this PIC code to evaluate the accuracy of emittance measurements using a solenoid focal scan following the quadrupole magnets. If the beam is round after the solenoids, the simulations indicate that the measurement is highly accurate, but it is substantially inaccurate for elliptical beams.


title: Revealing long-range multi-particle collectivity in small collision systems via subevent cumulants

abstract: Multi-particle azimuthal cumulants, often used to study collective flow in high-energy heavy-ion collisions, have recently been applied in small collision systems such as $pp$ and $p$+A to extract the second-order azimuthal harmonic flow $v_2$. Recent observation of four-, six- and eight-particle cumulants with "correct sign" $c_2\{4\}<0, c_2\{6\}>0, c_2\{8\}<0$ and approximate equality of the inferred single-particle harmonic flow, $v_2\{4\}\approx v_2\{6\}\approx v_2\{8\}$, have been used as strong evidence for a collective emission of all soft particles produced in the collisions. We show that these relations in principle could be violated due to the non-Gaussianity in the event-by-event fluctuation of flow and/or non-flow. Furthermore, we show, using $pp$ events generated with the PYTHIA model, that $c_2\{2k\}$ obtained with standard cumulant method are dominated by non-flow from dijets. An alternative cumulant method based on two or more $\eta$-separated subevents is proposed to suppress the dijet contribution. The new method is shown to be able to recover a flow signal as low as 4\% imposed on the PYTHIA events, independently of how the event activity class is defined. Therefore the subevent cumulant method offers a more robust way of studying collectivity based on the existence of long-range azimuthal correlations between multiple distinct $\eta$ ranges. The prospect of using the subevent cumulants to study collective flow in A+A collisions, in particular its longitudinal dynamics, is discussed.


title: Delay-Optimal Scheduling for Queueing Systems with Switching Overhead

abstract: We study the scheduling polices for asymptotically optimal delay in queueing systems with switching overhead. Such systems consist of a single server that serves multiple queues, and some capacity is lost whenever the server switches to serve a different set of queues. The capacity loss due to this switching overhead can be significant in many emerging applications, and needs to be explicitly addressed in the design of scheduling policies. For example, in 60GHz wireless networks with directional antennas, base stations need to train and reconfigure their beam patterns whenever they switch from one client to another. Considerable switching overhead can also be observed in many other queueing systems such as transportation networks and manufacturing systems. While the celebrated Max-Weight policy achieves asymptotically optimal average delay for systems without switching overhead, it fails to preserve throughput-optimality, let alone delay-optimality, when switching overhead is taken into account. We propose a class of Biased Max-Weight scheduling policies that explicitly takes switching overhead into account. The Biased Max-Weight policy can use either queue length or head-of-line waiting time as an indicator of the system status. We prove that our policies not only are throughput-optimal, but also can be made arbitrarily close to the asymptotic lower bound on average delay. To validate the performance of the proposed policies, we provide extensive simulation with various system topologies and different traffic patterns. We show that the proposed policies indeed achieve much better delay performance than that of the state-of-the-art policy.


title: The Plateau-Rayleigh instability in solids is a simple phase separation

abstract: A long elastic cylinder, radius $a$ and shear-modulus $\mu$, becomes unstable given sufficient surface tension $\gamma$. We show this instability can be simply understood by considering the energy, $E(\lambda)$, of such a cylinder subject to a homogenous longitudinal stretch $\lambda$. Although $E(\lambda)$ has a unique minimum, if surface tension is sufficient ($\Gamma\equiv\gamma/(a\mu)>\sqrt{32}$) it looses convexity in a finite region. We use a Maxwell construction to show that, if stretched into this region, the cylinder will phase separate into two segments with different stretches $\lambda_1$ and $\lambda_2$. Our model thus explains why the instability has infinite wavelength, and allows us to calculate the instability's sub-critical hysteresis loop (as a function of imposed stretch), showing that instability proceeds with constant amplitude and at constant (positive) tension as the cylinder is stretched between $\lambda_1$ and $\lambda_2$. We use full nonlinear finite-element calculations to verify these predictions, and to characterize the interface between the two phases. Near $\Gamma=\sqrt{32}$ the length of such an interface diverges introducing a new length-scale and allowing us to construct a 1-D effective theory. This treatment yields an analytic expression for the interface itself, revealing its characteristic length grows as $l_{wall}\sim a/\sqrt{\Gamma-\sqrt{32}}$.


title: An automated system to measure the quantum efficiency of CCDs for astronomy

abstract: We describe a system to measure the Quantum Efficiency in the wavelength range of 300 nm to 1100 nm of 40x40 mm n-channel CCD sensors for the construction of the 3.2 gigapixel LSST focal plane. The technique uses a series of instrument to create a very uniform flux of photons of controllable intensity in the wavelength range of interest across the face the sensor. This allows the absolute Quantum Efficiency to be measured with an accuracy in the 7% range. This system will be part of a production facility at Brookhaven National Lab for the basic component of the LSST camera.


title: On Higher Order Positive Differential Energy Operator

abstract: The higher order differential energy operator (DEO), denoted via $\Upsilon_k(x)$, is an extension to the second order famous Teager-Kaiser operator. The DEO helps measuring the higher order gauge of energy of a signal which is useful for AM-FM demodulation. However, the energy criterion defined by the DEO is not compliant with the presumption of positivity of energy. In this paper we introduce a higher order operator called Positive Differential Energy Operator (PDEO). This operator which can be obtained using alternative recursive relations, resolves the energy sign problem. The simulations demonstrate that the proposed operator can outperform DEOs in terms of Average Signal to Error Ratio (ASER) in AM/FM demodulation.


title: Additional cases of positive twisted torus knots

abstract: A twisted torus knot is a knot obtained from a torus knot by twisting adjacent strands by full twists. The twisted torus knots lie in $F$, the genus 2 Heegaard surface for $S^3$. Primitive/primitive and primitive/Seifert knots lie in $F$ in a particular way. Dean gives sufficient conditions for the parameters of the twisted torus knots to ensure they are primitive/primitive or primitive/Seifert. Using Dean's conditions, Doleshal shows that there are infinitely many twisted torus knots that are fibered and that there are twisted torus knots with distinct primitive/Seifert representatives with the same slope in $F$. In this paper, we extend Doleshal's results to show there is a four parameter family of positive twisted torus knots. Additionally, we provide new examples of twisted torus knots with distinct representatives with the same surface slope in $F$.


title: Formal Analysis of SEU Mitigation for Early Dependability and Performability Analysis of FPGA-based Space Applications

abstract: SRAM-based FPGAs are increasingly popular in the aerospace industry due to their field programmability and low cost. However, they suffer from cosmic radiation induced Single Event Upsets (SEUs). In safety-critical applications, the dependability of the design is a prime concern since failures may have catastrophic consequences. An early analysis of the relationship between dependability metrics, performability-area trade-off, and different mitigation techniques for such applications can reduce the design effort while increasing the design confidence. This paper introduces a novel methodology based on probabilistic model checking, for the analysis of the reliability, availability, safety and performance-area tradeoffs of safety-critical systems for early design decisions. Starting from the high-level description of a system, a Markov reward model is constructed from the Control Data Flow Graph (CDFG) and a component characterization library targeting FPGAs. The proposed model and exhaustive analysis capture all the failure states (based on the fault detection coverage) and repairs possible in the system. We present quantitative results based on an FIR filter circuit to illustrate the applicability of the proposed approach and to demonstrate that a wide range of useful dependability and performability properties can be analyzed using the proposed methodology. The modeling results show the relationship between different mitigation techniques and fault detection coverage, exposing their direct impact on the design for early decisions.


title: AC-Biased Shift Registers as Fabrication Process Benchmark Circuits and Flux Trapping Diagnostic Tool

abstract: We develop an ac-biased shift register introduced in our previous work (V.K. Semenov et al., IEEE Trans. Appl. Supercond., vol. 25, no. 3, 1301507, June 2015) into a benchmark circuit for evaluation of superconductor electronics fabrication technology. The developed testing technique allows for extracting margins of all individual cells in the shift register, which in turn makes it possible to estimate statistical distribution of Josephson junctions in the circuit. We applied this approach to successfully test registers having 8, 16, 36, and 202 thousand cells and, respectively, about 33000, 65000, 144000, and 809000 Josephson junctions. The circuits were fabricated at MIT Lincoln Laboratory, using a fully planarized process, 0.4 {\mu}m inductor linewidth, and 1.33x10^6 cm^-2 junction density. They are presently the largest operational superconducting SFQ circuits ever made. The developed technique distinguishes between hard defects (fabrication-related) and soft defects (measurement-related) and locates them in the circuit. The soft defects are specific to superconducting circuits and caused by magnetic flux trapping either inside the active cells or in the dedicated flux-trapping moats near the cells. The number and distribution of soft defects depend on the ambient magnetic field and vary with thermal cycling even if done in the same magnetic environment.


title: Single-Sided Field-Free Line Generator Magnet for Multidimensional Magnetic Particle Imaging

abstract: Magnetic Particle Imaging (MPI) is an emerging medical imaging modality that is not yet adopted by clinical practice. Most of the working MPI prototypes including commercial-grade research MPI scanners utilize cylindrical bores that limit the access to the scanner and the imaging volume. Recently a single-sided or an asymmetric device that is based on a field-free point (FFP) coplanar coil topology has been introduced that shows promise in alleviating access constraint issues. In this paper, we present a simulation study of selection coils for a novel design of a single-sided MPI device that has an advantage of a more sensitive field-free line (FFL) topology.


title: Exploring Model Predictive Control to Generate Optimal Control Policies for HRI Dynamical Systems

abstract: We model Human-Robot-Interaction (HRI) scenarios as linear dynamical systems and use Model Predictive Control (MPC) with mixed integer constraints to generate human-aware control policies. We motivate the approach by presenting two scenarios. The first involves an assistive robot that aims to maximize productivity while minimizing the human's workload, and the second involves a listening humanoid robot that manages its eye contact behavior to maximize "connection" and minimize social "awkwardness" with the human during the interaction. Our simulation results show that the robot generates useful behaviors as it finds control policies to minimize the specified cost function. Further, we implement the second scenario on a humanoid robot and test the eye contact scenario with 48 human participants to demonstrate and evaluate the desired controller behavior. The humanoid generated 25% more eye contact when it was told to maximize connection over when it was told to maximize awkwardness. However, despite showing the desired behavior, there was no statistical difference between the participant's perceived connection with the humanoid.


title: Performance Analysis of Low-Density Parity-Check Codes over 2D Interference Channels via Density Evolution

abstract: The theoretical analysis of detection and decoding of low-density parity-check (LDPC) codes transmitted over channels with two-dimensional (2D) interference and additive white Gaussian noise (AWGN) is provided in this paper. The detection and decoding system adopts the joint iterative detection and decoding scheme (JIDDS) in which the log-domain sum-product algorithm is adopted to decode the LDPC codes. The graph representations of the JIDDS are explained. Using the graph representations, we prove that the message-flow neighborhood of the detection and decoding system will be tree-like for a sufficiently long code length. We further confirm that the performance of the JIDDS will concentrate around the performance in which message-flow neighborhood is tree-like. Based on the tree-like message-flow neighborhood, we employ a modified density evolution algorithm to track the message densities during the iterations. A threshold is calculated using the density evolution algorithm which can be considered as the theoretical performance limit of the system. Simulation results demonstrate that the modified density evolution is effective in analyzing the performance of 2D interference systems.


title: Treatment of the singularities of the Stermheimer function for a layered electron gas

abstract: Isolating the singularity in the Green's function solution of the inhomogeneous, differential equation for the Stermheimer function of a layered electron gas, permits the construction of an approximate solution of the Stermheimer function which includes the exact contributions from its singularities.


title: First-principles-based method for electron localization: Application to monolayer hexagonal boron nitride

abstract: We present a first-principles-based many-body typical medium dynamical cluster approximation method for characterizing electron localization in disordered structures. This method applied to monolayer hexagonal boron nitride shows that the presence of a boron vacancies could turn this wide-gap insulator into a correlated metal. Depending on the strength of the electron interactions, these calculations suggest that conduction could be obtained at a boron vacancy concentration as low as $1.0\%$. We also explore the distribution of the local density of states, a fingerprint of spatial variations, which allows localized and delocalized states to be distinguished. The presented method enables the study of disorder-driven insulator-metal transitions not only in $h$-BN but also in other physical materials.


title: Relations between Schramm spaces and generalized Wiener classes

abstract: We give necessary and sufficient conditions for the embeddings $\Lambda\text{BV}^{(p)}\subseteq \Gamma\text{BV}^{(q_n\uparrow q)}$ and $\Phi\text{BV}\subseteq\text{BV}^{(q_n\uparrow q)}$. As a consequence, a number of results in the literature, including a fundamental theorem of Perlman and Waterman, are simultaneously extended.


title: On Gauge Invariant Cosmological Perturbations in UV-modified Horava Gravity

abstract: We consider gauge invariant cosmological perturbations in UV-modified, z=3 Horava gravity with one scalar matter field, which has been proposed as a renormalizable gravity theory without the ghost problem in four dimensions. In order to exhibit its dynamical degrees of freedom, we consider the Hamiltonian reduction method and find that, by solving "all" the constraint equations, the degrees of freedom are the same as those of Einstein gravity: One scalar and two tensor (graviton) modes when a scalar matter field presents. However, we confirm that there is no extra graviton modes and general relativity is recovered in IR, which achieves the consistency of the model. From the UV-modification terms which break the detailed balance condition in UV, we obtain scale-invariant power spectrums for "non"-inflationary backgrounds, like the power-law expansions, without knowing the details of early expansion history of Universe. This could provide a new framework for the Big Bang cosmology. Moreover, we find that "tensor and scalar fluctuations travel differently in UV, generally". We present also some clarifying remarks about confusing points in the literatures.


title: $\beta$-decay half-lives and $\beta$-delayed neutron emission probabilities for several isotopes of Au, Hg, Tl, Pb and Bi, beyond N=126

abstract: $Background:$ Previous measurements of $\beta$-delayed neutron emitters comprise around 230 nuclei, spanning from the $^{8}$He up to $^{150}$La. Apart from $^{210}$Tl, with a minuscule branching ratio of 0.007\%, no other neutron emitter is measured yet beyond $A=150$. Therefore new data are needed, particularly in the heavy mass region around N=126, in order to guide theoretical models and to understand the formation of the third r-process peak at $A\sim195$. $Purpose:$ To measure both, $\beta$-decay half-lives and neutron branching ratios of several neutron-rich Au, Hg, Tl, Pb and Bi isotopes beyond $N=126$. $Method:$ Ions of interest are produced by fragmentation of a $^{238}$U beam, selected and identified via the GSI-FRS fragment separator. A stack of segmented silicon detectors (SIMBA) is used to measure ion-implants and $\beta$-decays. An array of 30 $^3$He tubes embedded in a polyethylene matrix (BELEN) is used to detect neutrons with high efficiency and selectivity. A self-triggered digital system is employed to acquire data and to enable time-correlations. The latter are analyzed with an analytical model and results for the half-lives and neutron-branching ratios are derived using the binned Maximum-Likelihood method. $Results:$ Twenty new $\beta$-decay half-lives are reported for $^{204-206}$Au, $^{208-211}$Hg,$^{211-216}$Tl,$^{215-218}$Pb and $^{218-220}$Bi, nine of them for the first time. Neutron emission probabilities are reported for $^{210,211}$Hg and $^{211-216}$Tl. $Conclusions:$ The new $\beta$-decay half-lives are in good agreement with previous measurements in this region. The measured neutron emission probabilities are comparable or smaller than values predicted by global models like RHB+RQRPA.


title: Finitely forcible graph limits are universal

abstract: The theory of graph limits represents large graphs by analytic objects called graphons. Graph limits determined by finitely many graph densities, which are represented by finitely forcible graphons, arise in various scenarios, particularly within extremal combinatorics. Lovasz and Szegedy conjectured that all such graphons possess a simple structure, e.g., the space of their typical vertices is always finite dimensional; this was disproved by several ad hoc constructions of complex finitely forcible graphons. We prove that any graphon is a subgraphon of a finitely forcible graphon. This dismisses any hope for a result showing that finitely forcible graphons possess a simple structure, and is surprising when contrasted with the fact that finitely forcible graphons form a meager set in the space of all graphons. In addition, since any finitely forcible graphon represents the unique minimizer of some linear combination of densities of subgraphs, our result also shows that such minimization problems, which conceptually are among the simplest kind within extremal graph theory, may in fact have unique optimal solutions with arbitrarily complex structure.


title: The index of singular zeros of harmonic mappings of anti-analytic degree one

abstract: We study harmonic mappings of the form $f(z) = h(z) - \overline{z}$, where $h$ is an analytic function. In particular we are interested in the index (a generalized multiplicity) of the zeros of such functions. Outside the critical set of $f$, where the Jacobian of $f$ is non-vanishing, it is known that this index has similar properties as the classical multiplicity of zeros of analytic functions. Little is known about the index of zeros on the critical set, where the Jacobian vanishes; such zeros are called singular zeros. Our main result is a characterization of the index of singular zeros, which enables one to determine the index directly from the power series of $h$.


title: Multiplication and Presence of Shielding Material from Time-Correlated Pulse-Height Measurements of Subcritical Plutonium Assemblies

abstract: We present the results from the first measurements of the Time-Correlated Pulse-Height (TCPH) distributions from 4.5 kg sphere of $\alpha$-phase weapons-grade plutonium metal in five configurations: bare, reflected by 1.27 cm and 2.54 cm of tungsten, and 2.54 cm and 7.62 cm of polyethylene. A new method for characterizing source multiplication and shielding configuration is also demonstrated. The method relies on solving for the underlying fission chain timing distribution that drives the spreading of the measured TCPH distribution. We found that a gamma distribution fits the fission chain timing distribution well and that the fit parameters correlate with both multiplication (rate parameter) and shielding material types (shape parameter). The source-to-detector distance was another free parameter that we were able to optimize, and proved to be the most well constrained parameter. MCNPX-PoliMi simulations were used to complement the measurements and help illustrate trends in these parameters and their relation to multiplication and the amount and type of material coupled to the subcritical assembly.


title: Deep Neural Networks for Czech Multi-label Document Classification

abstract: This paper is focused on automatic multi-label document classification of Czech text documents. The current approaches usually use some pre-processing which can have negative impact (loss of information, additional implementation work, etc). Therefore, we would like to omit it and use deep neural networks that learn from simple features. This choice was motivated by their successful usage in many other machine learning fields. Two different networks are compared: the first one is a standard multi-layer perceptron, while the second one is a popular convolutional network. The experiments on a Czech newspaper corpus show that both networks significantly outperform baseline method which uses a rich set of features with maximum entropy classifier. We have also shown that convolutional network gives the best results.


title: Recent Developments in CMS Fast Simulation

abstract: CMS has developed a fast detector simulation package, which serves as a fast and reliable alternative to the detailed GEANT4-based (full) simulation, and enables efficient simulation of large numbers of standard model and new physics events. Fast simulation becomes particularly important with the current increase in the LHC luminosity. Here, I will discuss the basic principles behind the CMS fast simulation framework, and how they are implemented in the different detector components in order to simulate and reconstruct sufficiently accurate physics objects for analysis. I will focus on recent developments in tracking and geometry interface, which improve the flexibility and emulation performance of the framework, and allow a better synchronization with the full simulation. I will then show how these developments have led to an improved agreement of basic analysis objects and event variables between fast and full simulation.


title: Potential relationship of epistemic games to group dynamics and learning orientations towards physics problem solving

abstract: Current investigations into pedagogical goals of introductory algebra-based physics students at the University of Central Arkansas, by learning orientation towards an in-class metacognitive group problem solving task, seek to determine possible relationships with attitudinal shifts and course performance. Students thus far have been untreated with known group-based learning pedagogies, so as to establish trends of common group habits, and ultimately to properly inform implementation of group-based pedagogies in reaction to these trends. However, students' group dynamics and learning orientations prove difficult to map to group-based measurements; an estimate of group learning orientation and preferred working group dynamic is here explored as a potential means of interpreting students' use of problem solving strategies. A means of "sampling" audiovisual data in a live classroom of several simultaneous groups is also presented as a way to estimate the frequency of chosen strategies to this end.


title: The multidegree of the multi-image variety

abstract: The multi-image variety is a subvariety of Gr(1,P^3)^n that models taking pictures with n rational cameras. We compute its cohomology class in the cohomology of Gr(1,P^3)^n, and from there its multidegree as a subvariety of (P^5)^n under the Pl\"ucker embedding.


title: Crossovers due to anisotropy and disorder in quantum critical fluctuations

abstract: Scaling relations are used to study cross-overs, due to anisotropic spin interactions or single ion anisotropy, and due to disorder, in the thermodynamics and correlation functions near quantum-critical transitions. The principal results are simple with a wide range of applications. The region of attraction to the stable anisotropic fixed point in the quantum-critical region is exponentially enhanced by the dynamical critical exponent $z$ compared to the region of attraction to the fixed point in the quantum disordered region. The result implies that, even for small anisotropy, the region of attraction to the stable incommensurate Ising or planar metallic anti-ferromagnetic critical points, which belong to the universality class of the XY model with $z \to \infty$, covers the entire quantum-critical region. In crossovers due to disorder, the instability of the pure fixed point in the quantum disordered region is exponentially enhanced by $z$ compared to that in the quantum-critical region. This result suggests that for some classes of disorder and for large enough $z$, one may find singularities in the correlations as a function of frequency and temperature down to very low temperatures even though the correlation length in space remains short range.


title: PRSP: A Plugin-based Framework for RDF Stream Processing

abstract: In this paper, we propose a plugin-based framework for RDF stream processing named PRSP. Within this framework, we can employ SPARQL query engines to process C-SPARQL queries with maintaining the high performance of those engines in a simple way. Taking advantage of PRSP, we can process large-scale RDF streams in a distributed context via distributed SPARQL engines. Besides, we can evaluate the performance and correctness of existing SPARQL query engines in handling RDF streams in a united way, which amends the evaluation of them ranging from static RDF (i.e., RDF graph) to dynamic RDF (i.e., RDF stream). Finally, within PRSP, we experimently evaluate the correctness and the performance on YABench. The experiments show that PRSP can still maintain the high performance of those engines in RDF stream processing although there are some slight differences among them.


title: Location Inference from Tweets using Grid-based Classification

abstract: The impact of social media and its growing association with the sharing of ideas and propagation of messages remains vital in everyday communication. Twitter is one effective platform for the dissemination of news and stories about recent events happening around the world. It has a continually growing database currently adopted by over 300 million users. In this paper we propose a novel grid-based approach employing supervised Multinomial Naive Bayes while extracting geographic entities from relevant user descriptions metadata which gives a spatial indication of the user location. To the best of our knowledge our approach is the first to make location inference from tweets using geo-enriched grid-based classification. Our approach performs better than existing baselines achieving more than 57% accuracy at city-level granularity. In addition we present a novel framework for content-based estimation of user locations by specifying levels of granularity required in pre-defined location grids.


title: The flip Markov chain for connected regular graphs

abstract: Mahlmann and Schindelhauer (2005) defined a Markov chain which they called $k$-Flipper, and showed that it is irreducible on the set of all connected regular graphs of a given degree (at least 3). We study the 1-Flipper chain, which we call the flip chain, and prove that the flip chain converges rapidly to the uniform distribution over connected $2r$-regular graphs with $n$ vertices, where $n\geq 8$ and $r = r(n)\geq 2$. Formally, we prove that the distribution of the flip chain will be within $\varepsilon$ of uniform in total variation distance after $\text{poly}(n,r,\log(\varepsilon^{-1}))$ steps. This polynomial upper bound on the mixing time is given explicitly, and improves markedly on a previous bound given by Feder et al.(2006). We achieve this improvement by using a direct two-stage canonical path construction, which we define in a general setting. This work has applications to decentralised networks based on random regular connected graphs of even degree, as a self-stabilising protocol in which nodes spontaneously perform random flips in order to repair the network.


title: Multichannel Molecular State and Rectified Short-range Boundary Condition for Spin-orbit Coupled Ultracold Fermions Near p-wave Resonances

abstract: We study the interplay of spin-orbit coupling (SOC) and strong p-wave interaction to the scattering property of spin-1/2 ultracold Fermi gases. Based on a two-channel square-well potential generating p-wave resonance, we show that the presence of an isotropic SOC, even for its length much longer than the potential range, can greatly modify the p-wave short-range boundary condition(BC). As a result, the conventional p-wave BC cannot predict the induced molecules near p-wave resonance, which can be fully destroyed to vanish due to strong interference between s- and p-wave channels. By analyzing the intrinsic reasons for the breakdown of conventional BC, we propose a new p-wave BC that can excellently reproduce the exact molecule solutions and also equally apply for a wide class of single-particle potentials besides SOC. This work reveals the significant effect of SOC to both the short- and long-range properties of fermions near p-wave resonance, paving the way for future exploring interesting few- and many-body physics in such system.


title: Semialgebraic metric spcaes and resolution of singularities of definable sets

abstract: Consider the semialgebraic structure over the real field. More generally, let an ominimal structure be over a real closed field. We show that a definable metric space X with a definable metric d is embedded into a Euclidean space so that its closure is compact and the metric on the image induced by d is extended to a definable metric on the closure if and only if the limit of d(r(t);r(t)) is 0 as t converges to 0 for any definable continuous curve r from (0, 1] to X (Theorem 1). We also find two compact semialgebraic metric spaces over the real field which are isometric but not semialgebraically isometric (Theorem 2). A version of blow up is the key to the proof of Theorem 1. Using it in the same way, we prove a resolution of singularities of definable sets (Theorem 3). We prove the theorems by a constructive procedure.


title: Characterizing the Gaussian coherence-breaking channel and its property with assistant entanglement inputs

abstract: In the paper, we give a characterization of arbitrary $n$-mode Gaussian coherence breaking channels (GCBCs) and show the tensor product of a GCBC and arbitrary a Gaussian channel maps all input states into product states. The inclusion relation among GCBCs, Gaussian positive partial transpose channels (GPPTC), entanglement breaking channels (GEBC) and Gaussian classical-quantum channels (GCQC) is displayed. Finally, we prove that Gaussian classical ($\chi$-) capacity for the tensor products of GCBCs and arbitrary Gaussian channels is additive.


title: Infinite-dimensional Stochastic Differential Equations with Symmetry

abstract: We review recent progress in the study of infinite-dimensional stochastic differential equations with symmetry. This paper contains examples arising from random matrix theory.


title: Network Inference from a Link-Traced Sample using Approximate Bayesian Computation

abstract: We present a new inference method based on approximate Bayesian computation for estimating parameters governing an entire network based on link-traced samples of that network. To do this, we first take summary statistics from an observed link-traced network sample, such as a recruitment network of subjects in a hard-to-reach population. Then we assume prior distributions, such as multivariate uniform, for the distribution of some parameters governing the structure of the network and behaviour of its nodes. Then, we draw many independent and identically distributed values for these parameters. For each set of values, we simulate a population network, take a link-traced sample from that network, and find the summary statistics for that sample. The statistics from the sample, and the parameters that eventually led to that sample, are collectively treated as a single point. We take a Kernel Density estimate of the points from many simulations, and observe the density across the hyperplane coinciding with the statistic values of the originally observed sample. This density function is treat as a posterior estimate of the paramaters of the network that provided the observed sample. We also apply this method to a network of precedence citations between legal documents, centered around cases overseen by the Supreme Court of Canada, is observed. The features of certain cases that lead to their frequent citation are inferred, and their effects estimated by ABC. Future work and extensions are also briefly discussed.


title: Sign-changing solutions for the stationary Kirchhoff problems involving the fractional Laplacian in \mathcal{R}^N

abstract: In this paper, we study the existence of least energy sign-changing solutions for a Kirchhoff-type problem involving the fractional Laplacian operator. By using the constraint variational method and quantitative deformation lemma, we obtain a least energy nodal solution ub for the given problem. Moreover, we show that the energy of ub is strictly larger than twice the ground state energy. We also give a convergence property of u_b as b goes to 0, where b is regarded as a parameter.


title: Breaking Locality Accelerates Block Gauss-Seidel

abstract: Recent work by Nesterov and Stich showed that momentum can be used to accelerate the rate of convergence for block Gauss-Seidel in the setting where a fixed partitioning of the coordinates is chosen ahead of time. We show that this setting is too restrictive, constructing instances where breaking locality by running non-accelerated Gauss-Seidel with randomly sampled coordinates substantially outperforms accelerated Gauss-Seidel with any fixed partitioning. Motivated by this finding, we analyze the accelerated block Gauss-Seidel algorithm in the random coordinate sampling setting. Our analysis captures the benefit of acceleration with a new data-dependent parameter which is well behaved when the matrix sub-blocks are well-conditioned. Empirically, we show that accelerated Gauss-Seidel with random coordinate sampling provides speedups for large scale machine learning tasks when compared to non-accelerated Gauss-Seidel and the classical conjugate-gradient algorithm.


title: 3D B_2 Model for Radiative Transfer Equation Part I: Modelling

abstract: We extend to three-dimensional space the approximate M_2 model for the slab geometry studied in our previous paper. The B_2 model therein, as a special case of the second order extended quadrature method of moments (EQMOM), is proved to be globally hyperbolic. The model we proposed here extends EQMOM to multiple dimensions following the idea to approximate the maximum entropy closure for the slab geometry case. Like the M_2 closure, the ansatz of the new model has the capacity to capture both isotropic and beam-like solutions, while the new model has fluxes in closed-form, thus is applicable to practical numerical simulations. The rotational invariance, realizability, and hyperbolicity of the model are studied.


title: Effects of an eccentric inner Jupiter on the dynamical evolution of icy body reservoirs in a planetary scattering scenario

abstract: We analyze the process of planetary scattering around M0-type stars. To do this, we carry out N-body simulations with three Jupiter-mass planets close to their instability limit together with an outer planetesimal disk. This paper focuses on the analysis of systems in which a single Jupiter-mass planet survives after the dynamical instability event. The small body reservoirs show different dynamical behaviors. In fact, our simulations produce particles on prograde and retrograde orbits, as well as particles whose orbital plane flips from prograde to retrograde and back again along their evolution. Such particles are called "Type-F particles". We find strong correlations between the inclination $i$ and the ascending node longitude $\Omega$ of such particles. First, $\Omega$ librates around 90$^{\circ}$ or/and 270$^{\circ}$. This property is very important since it represents a necessary and sufficient condition for the flipping of an orbit. Moreover, the libration periods of $i$ and $\Omega$ are equal and they are out to phase by a quarter period. We also remark that the larger the libration amplitude of $i$, the larger the libration amplitude of $\Omega$. Finally, we analyze the initial conditions of Type-F particles of all our simulations immediately after the dynamical instability event, when a single Jupiter-mass planet survives in the system. We carry out this study with the goal to determine the parameter space that lead to the flipping of an orbit. Our results suggest that the orbit of a test particle can flip for any value of its initial eccentricity, although we found only two Type-F particles with initial inclinations $i <$ 17$^{\circ}$. Moreover, our study indicates that the minimum value of the inclination of the Type-F particles in a given system decreases with an increase in the eccentricity of the giant planet.


title: Long Timescale Credit Assignment in NeuralNetworks with External Memory

abstract: Credit assignment in traditional recurrent neural networks usually involves back-propagating through a long chain of tied weight matrices. The length of this chain scales linearly with the number of time-steps as the same network is run at each time-step. This creates many problems, such as vanishing gradients, that have been well studied. In contrast, a NNEM's architecture recurrent activity doesn't involve a long chain of activity (though some architectures such as the NTM do utilize a traditional recurrent architecture as a controller). Rather, the externally stored embedding vectors are used at each time-step, but no messages are passed from previous time-steps. This means that vanishing gradients aren't a problem, as all of the necessary gradient paths are short. However, these paths are extremely numerous (one per embedding vector in memory) and reused for a very long time (until it leaves the memory). Thus, the forward-pass information of each memory must be stored for the entire duration of the memory. This is problematic as this additional storage far surpasses that of the actual memories, to the extent that large memories on infeasible to back-propagate through in high dimensional settings. One way to get around the need to hold onto forward-pass information is to recalculate the forward-pass whenever gradient information is available. However, if the observations are too large to store in the domain of interest, direct reinstatement of a forward pass cannot occur. Instead, we rely on a learned autoencoder to reinstate the observation, and then use the embedding network to recalculate the forward-pass. Since the recalculated embedding vector is unlikely to perfectly match the one stored in memory, we try out 2 approximations to utilize error gradient w.r.t. the vector in memory.