title: Investigating the Magnetic Imprints of Major Solar Eruptions with SDO/HMI High-Cadence Vector Magnetograms

abstract: The solar active region photospheric magnetic field evolves rapidly during major eruptive events, suggesting appreciable feedback from the corona. Previous studies of these "magnetic imprints" are mostly based on line-of-sight only or lower-cadence vector observations; a temporally resolved depiction of the vector field evolution is hitherto lacking. Here, we introduce the high-cadence (90~s or 135~s) vector magnetogram dataset from the Helioseismic and Magnetic Imager (HMI) that is well suited for investigating the phenomenon. These observations allow quantitative characterization of the permanent, step-like changes that are most pronounced in the horizontal field component ($B_h$). A highly structured pattern emerges from analysis of an archetypical event, \texttt{SOL2011-02-15T01:56}, where $B_h$ near the main polarity inversion line increases significantly during the earlier phase of the associated flare with a time scale of several minutes, while $B_h$ in the periphery decreases at later times with smaller magnitudes and a slightly longer time scale. The dataset also allows effective identification of the "magnetic transient" artifact, where enhanced flare emission alters the Stokes profiles and the inferred magnetic field becomes unreliable. Our results provide insights on the momentum processes in solar eruptions. The dataset may also be useful to the study of sunquakes and data-driven modeling of the corona.


title: A Converse to Banach's Fixed Point Theorem and its CLS Completeness

abstract: Banach's fixed point theorem for contraction maps has been widely used to analyze the convergence of iterative methods in non-convex problems. It is a common experience, however, that iterative maps fail to be globally contracting under the natural metric in their domain, making the applicability of Banach's theorem limited. We explore how generally we can apply Banach's fixed point theorem to establish the convergence of iterative methods when pairing it with carefully designed metrics. Our first result is a strong converse of Banach's theorem, showing that it is a universal analysis tool for establishing global convergence of iterative methods to unique fixed points, and for bounding their convergence rate. In other words, we show that, whenever an iterative map globally converges to a unique fixed point, there exists a metric under which the iterative map is contracting and which can be used to bound the number of iterations until convergence. We illustrate our approach in the widely used power method, providing a new way of bounding its convergence rate through contraction arguments. We next consider the computational complexity of Banach's fixed point theorem. Making the proof of our converse theorem constructive, we show that computing a fixed point whose existence is guaranteed by Banach's fixed point theorem is CLS-complete. We thus provide the first natural complete problem for the class CLS, which was defined in [Daskalakis, Papadimitriou 2011] to capture the complexity of problems such as P-matrix LCP, computing KKT-points, and finding mixed Nash equilibria in congestion and network coordination games.


title: Chiral Vortical Effect for Bosons

abstract: The thermal contribution to the chiral vortical effect is believed to be related to the axial anomaly in external gravitational fields. We use the universality of the spin-gravity interaction to extend this idea to a wider set of phenomena. We consider the Kubo formula at weak coupling for the spin current of a vector field and derive a novel anomalous effect caused by the medium rotation: the chiral vortical effect for bosons. The effect consists in a spin current of vector bosons along the angular velocity of the medium. We argue that it has the same anomalous nature as in the fermionic case and show that this effect provides a mechanism for helicity transfer, from flow helicity to magnetic helicity.


title: FUSE Spectroscopy of the Accreting Hot Components in Symbiotic Variables

abstract: We have conducted a spectroscopic analysis of the far ultraviolet archival spectra of four symbiotic variables, EG And, AE Ara, CQ Dra and RW Hya. RW Hya and EG And have never had a recorded outburst while CQ Dra and AE Ara have outburst histories. We analyze these systems while they are in quiescence in order to help reveal the physical properties of their hot components via comparisons of the observations with optically thick accretion disk models and NLTE model white dwarf photospheres. We have extended the wavelength coverage down to the Lyman Limit with FUSE spectra. We find that the hot component in RW Hya is a low mass white dwarf with a surface temperature of 160,000K. We re-examine whether or not the symbiotic system CQ Dra is a triple system with a red giant transferring matter to a hot component made up of a cataclysmic variable in which the white dwarf has a surface temperature as low as $\sim$20,000K. The very small size of the hot component contributing to the shortest wavelengths of the FUSE spectrum of CQ Dra agrees with an optically thick and geometrically thin ($\sim$4\% of the WD surface) hot ($\sim 120,000$K) boundary layer. Our analysis of EG And reveals that its hot component is a hot, bare, low mass white dwarf with a surface temperature of 80-95,000K, with a surface gravity $\log(g)= 7.5$. For AE Ara, we also find that a low gravity ($\log(g) \sim 6$) hot ($T \sim 130,000$K) WD accounts for the hot component.


title: On the inducibility of cycles

abstract: In 1975 Pippenger and Golumbic proved that any graph on $n$ vertices admits at most $2e(n/k)^k$ induced $k$-cycles. This bound is larger by a multiplicative factor of $2e$ than the simple lower bound obtained by a blow-up construction. Pippenger and Golumbic conjectured that the latter lower bound is essentially tight. In the present paper we establish a better upper bound of $(128e/81) \cdot (n/k)^k$. This constitutes the first progress towards proving the aforementioned conjecture since it was posed.


title: Improving high-pass fusion method using wavelets

abstract: In an appropriate image fusion method, spatial information of the panchromatic image is injected into the multispectral images such that the spectral information is not distorted. The high-pass modulation method is a successful method in image fusion. However, the main drawback of this method is that this technique uses the boxcar filter to extract the high frequency information of the panchromatic image. Using the boxcar filter introduces the ringing effect into the fused image. To cope with this problem, we use the wavelet transform instead of boxcar filters. Then, the results of the proposed method and those of other methods such as, Brovey, IHS, and PCA ones are compared. Experiments show the superiority of the proposed method in terms of correlation coefficient and mutual information.


title: Variability of the symbiotic X-ray binary GX 1+4: Enhanced activity near periastron passage

abstract: Context. GX 1+4 belongs to a rare class of X-ray binaries with red giant donors, symbiotic X-ray binaries. The system has a history of complicated variability on multiple timescales in the optical light and X-rays. The nature of this variability remains poorly understood. Aims. We study variability of GX 1+4 on long time-scale in X-ray and optical bands. Methods. The presented X-ray observations are from INTEGRAL Soft Gamma-Ray Imager and RXTE All Sky Monitor. The optical observations are from INTEGRAL Optical Monitoring Camera. Results. The variability of GX 1+4 both in optical light and hard X-ray emission (>17 keV) is dominated by ~50-70d quasi-periodic changes. The amplitude of this variability is highest during the periastron passage, while during the potential neutron star eclipse the system is always at minimum, which confirms the 1161d orbital period that has had been proposed for the system based on radial velocity curve. Neither the quasi-periodic variability or the orbital period are detected in soft X-ray emission (1.3-12.2 keV), where the binary shows no apparent periodicity.


title: Warm Modified Chaplygin Gas Shaft Inflation

abstract: In this paper, we examine the possible realization of a new family of inflation called "shaft inflation" by assuming the modified Chaplygin gas model and tachyon scalar field. We also consider the special form of dissipative coefficient as $\Gamma={a_0}\frac{T^{3}}{\phi^{2 }}$ and calculate the various inflationary parameters in the scenario of strong and weak dissipative regimes. In order to examine the behavior of inflationary parameters, the planes of $n_s - \phi,~n_s - r$ and $n_s - \alpha_s$ (where $n_s,~\alpha_s,~r$ and $\phi$ represent spectral index, its running, tensor to scalar ratio and scalar field respectively) are being developed which lead to the constraints: $r< 0.11$, $n_s=0.96\pm0.025$ and $\alpha_s =-0.019\pm0.025$. It is quite interesting that these results of inflationary parameters are compatible with BICEP$2$, WMAP $(7+9)$ and recent Planck data.


title: Quantum Radiation Properties of General Nonstationary Black Hole

abstract: Using the generalized tortoise coordinate transformations the quantum radiation properties of Klein-Gordon scalar particles, Maxwell's electromagnetic field equations and Dirac equations are investigated in general non-stationary black hole. The locations of the event horizon and the Hawking temperature depend on both time and angles. A new extra coupling effect is observed in the thermal radiation spectrum of Maxwell's equations and Dirac equations which is absent in the thermal radiation spectrum of scalar particles. We also observe that the chemical potential derived from scalar particles is equal to the highest energy of the negative energy state of the scalar particle in the non-thermal radiation in general non-stationary black hole. Applying generalized tortoise coordinate transformation a constant term $\xi$ is produced in the expression of thermal radiation in general non-stationary black hole. It indicates that generalized tortoise coordinate transformation is more accurate and reliable in the study of thermal radiation of black hole.


title: Dispersion relation for hadronic light-by-light scattering: two-pion contributions

abstract: In this third paper of a series dedicated to a dispersive treatment of the hadronic light-by-light (HLbL) tensor, we derive a partial-wave formulation for two-pion intermediate states in the HLbL contribution to the anomalous magnetic moment of the muon $(g-2)_\mu$, including a detailed discussion of the unitarity relation for arbitrary partial waves. We show that obtaining a final expression free from unphysical helicity partial waves is a subtle issue, which we thoroughly clarify. As a by-product, we obtain a set of sum rules that could be used to constrain future calculations of $\gamma^*\gamma^*\to\pi\pi$. We validate the formalism extensively using the pion-box contribution, defined by two-pion intermediate states with a pion-pole left-hand cut, and demonstrate how the full known result is reproduced when resumming the partial waves. Using dispersive fits to high-statistics data for the pion vector form factor, we provide an evaluation of the full pion box, $a_\mu^{\pi\text{-box}}=-15.9(2)\times 10^{-11}$. As an application of the partial-wave formalism, we present a first calculation of $\pi\pi$-rescattering effects in HLbL scattering, with $\gamma^*\gamma^*\to\pi\pi$ helicity partial waves constructed dispersively using $\pi\pi$ phase shifts derived from the inverse-amplitude method. In this way, the isospin-$0$ part of our calculation can be interpreted as the contribution of the $f_0(500)$ to HLbL scattering in $(g-2)_\mu$. We argue that the contribution due to charged-pion rescattering implements corrections related to the corresponding pion polarizability and show that these are moderate. Our final result for the sum of pion-box contribution and its $S$-wave rescattering corrections reads $a_\mu^{\pi\text{-box}} + a_{\mu,J=0}^{\pi\pi,\pi\text{-pole LHC}}=-24(1)\times 10^{-11}$.


title: On the robustness of the H$\beta$ Lick index as a cosmic clock in passive early-type galaxies

abstract: We examine the H$\beta$ Lick index in a sample of $\sim 24000$ massive ($\rm log(M/M_{\odot})>10.75$) and passive early-type galaxies extracted from SDSS at z<0.3, in order to assess the reliability of this index to constrain the epoch of formation and age evolution of these systems. We further investigate the possibility of exploiting this index as "cosmic chronometer", i.e. to derive the Hubble parameter from its differential evolution with redshift, hence constraining cosmological models independently of other probes. We find that the H$\beta$ strength increases with redshift as expected in passive evolution models, and shows at each redshift weaker values in more massive galaxies. However, a detailed comparison of the observed index with the predictions of stellar population synthesis models highlights a significant tension, with the observed index being systematically lower than expected. By analyzing the stacked spectra, we find a weak [NII]$\lambda6584$ emission line (not detectable in the single spectra) which anti-correlates with the mass, that can be interpreted as a hint of the presence of ionized gas. We estimated the correction of the H$\beta$ index by the residual emission component exploiting different approaches, but find it very uncertain and model-dependent. We conclude that, while the qualitative trends of the observed H$\beta$-z relations are consistent with the expected passive and downsizing scenario, the possible presence of ionized gas even in the most massive and passive galaxies prevents to use this index for a quantitative estimate of the age evolution and for cosmological applications.


title: X-ray spectroscopy of the z=6.4 quasar J1148+5251

abstract: We present the 78-ks Chandra observations of the $z=6.4$ quasar SDSS J1148+5251. The source is clearly detected in the energy range 0.3-7 keV with 42 counts (with a significance $\gtrsim9\sigma$). The X-ray spectrum is best-fitted by a power-law with photon index $\Gamma=1.9$ absorbed by a gas column density of $\rm N_{\rm H}=2.0^{+2.0}_{-1.5}\times10^{23}\,\rm cm^{-2}$. We measure an intrinsic luminosity at 2-10 keV and 10-40 keV equal to $\sim 1.5\times 10^{45}~\rm erg~s^{-1}$, comparable with luminous local and intermediate-redshift quasar properties. Moreover, the X-ray to optical power-law slope value ($\alpha_{\rm OX}=-1.76\pm 0.14$) of J1148 is consistent with the one found in quasars with similar rest-frame 2500 \AA ~luminosity ($L_{\rm 2500}\sim 10^{32}~\rm erg~s^{-1}$\AA$^{-1}$). Then we use Chandra data to test a physically motivated model that computes the intrinsic X-ray flux emitted by a quasar starting from the properties of the powering black hole and assuming that X-ray emission is attenuated by intervening, metal-rich ($Z\geq \rm Z_{\odot}$) molecular clouds distributed on $\sim$kpc scales in the host galaxy. Our analysis favors a black hole mass $M_{\rm BH} \sim 3\times 10^9 \rm M_\odot$ and a molecular hydrogen mass $M_{\rm H_2}\sim 2\times 10^{10} \rm M_\odot$, in good agreement with estimates obtained from previous studies. We finally discuss strengths and limits of our analysis.


title: VALES: III. The calibration between the dust continuum and interstellar gas content of star-forming galaxies

abstract: We present the calibration between the dust continuum luminosity and interstellar gas content obtained from the Valpara\'{i}so ALMA Line Emission Survey (VALES) sample of 67 main-sequence star-forming galaxies at 0.02<$z$<0.35. We use CO(1-0) observations from the Atacama Large Millimetre/submillimetre Array (ALMA) to trace the molecular gas mass, $M_{\mathrm{H}_{2}}$, and estimate the rest-frame monochromatic luminosity at 850 $\mu$m, $L_{\nu_{850}}$, by extrapolating the dust continuum from MAGPHYS modelling of the far-ultraviolet to submillimetre spectral energy distribution sampled by the Galaxy And Mass Assembly (GAMA) survey. Adopting $\alpha_{\rm CO}$ = 6.5 (K km s$^{-1}$ pc$^{2}$)$^{-1}$, the average ratio of $L_{\nu_{850}}/M_{\mathrm{H}_{2}}$ = (6.4$\pm$1.4)$\times10^{19}$ erg s$^{-1}$ Hz$^{-1}$ $\mathrm{M}_{\odot}^{-1}$, in excellent agreement with literature values. We obtain a linear fit of $\log_{10}$ ($M_{\mathrm{H}_{2}}/\mathrm{M}_{\odot}$) = (0.92$\pm$0.02) $\log_{10}$ ($L_{\nu_{850}}$/erg s$^{-1}$ Hz$^{-1}$)-(17.31$\pm$0.59). We provide relations between $L_{\nu_{850}}$, $M_{\mathrm{H}_{2}}$ and $M_{\mathrm{ISM}}$ when combining the VALES and literature samples, and adopting a Galactic $\alpha_{\rm CO}$ value.


title: The Mass and Momentum Outflow Rates of Photoionized Galactic Outflows

abstract: Galactic outflows are believed to play an important role in regulating star formation in galaxies, but estimates of the outflowing mass and momentum have historically been based on uncertain assumptions. Here, we measure the mass, momentum, and energy outflow rates of seven nearby star-forming galaxies using ultraviolet absorption lines and observationally motivated estimates for the density, metallicity, and radius of the outflow. Low-mass galaxies generate outflows faster than their escape velocities with mass outflow rates up to twenty times larger than their star formation rates. These outflows from low-mass galaxies also have momenta larger than provided from supernovae alone, indicating that multiple momentum sources drive these outflows. Only 1-20\% of the supernovae energy is converted into kinetic energy, and this fraction decreases with increasing stellar mass such that low-mass galaxies drive more efficient outflows. We find scaling relations between the outflows and the stellar mass of their host galaxies (M$_\ast$) at the 2-3$\sigma$ significance level. The mass-loading factor, or the mass outflow rate divided by the star formation rate, scales as M$_\ast^{-0.4}$ and with the circular velocity as v$_\mathrm{circ}^{-1.6}$. The scaling of the mass-loading factor is similar to recent simulations, but the observations are a factor of five smaller, possibly indicating that there is a substantial amount of unprobed gas in a different ionization phase. The outflow momenta are consistent with a model where star formation drives the outflow while gravity counteracts this acceleration.


title: Temporary Capture of Asteroids by an Eccentric Planet

abstract: We have investigated the probability of temporary capture of asteroids in eccentric orbits by a planet in a circular or an eccentric orbit through analytical and numerical calculations. We found that in the limit of the circular orbit, the capture probability is $\sim 0.1\%$ of encounters to the planet's Hill sphere, independent of planetary mass and semimajor axis. In general, the temporary capture becomes more difficult as the planet's eccentricity ($e_{\rm p}$) increases. We found that the capture probability is almost independent of $e_{\rm p}$ until a critical value ($e_{\rm p}^{\rm c}$) that is given by $\simeq$ 5 times the Hill radius scaled by the planet's semimajor axis. For $e_{\rm p} > e_{\rm p}^{\rm c}$, the probability decreases approximately in proportion to $e_{\rm p}^{-1}$. The current orbital eccentricity of Mars is several times larger than $e_{\rm p}^{\rm c}$. However, since the range of secular change in Martian eccentricity overlaps $e_{\rm p}^{\rm c}$, the capture of minor bodies by the past Mars is not ruled out.


title: Detailed abundances from integrated-light spectroscopy: Milky Way globular clusters

abstract: We test the performance of our analysis technique for integrated-light spectra by applying it to seven well-studied Galactic GCs that span a wide range of metallicities. Integrated-light spectra were obtained by scanning the slit of the UVES spectrograph on the ESO Very Large Telescope across the half-light diameters of the clusters. We modelled the spectra using resolved HST colour-magnitude diagrams (CMDs), as well as theoretical isochrones, in combination with standard stellar atmosphere and spectral synthesis codes. The abundances of Fe, Na, Mg, Ca, Ti, Cr, and Ba were compared with literature data for individual stars in the clusters. The typical differences between iron abundances derived from our integrated-light spectra and those compiled from the literature are less than 0.1 dex. A larger difference is found for one cluster (NGC 6752), and is most likely caused primarily by stochastic fluctuations in the numbers of bright red giants within the scanned area. As expected, the alpha-elements (Ca, Ti) are enhanced by about 0.3 dex compared to the Solar-scaled composition, while the [Cr/Fe] ratios are close to Solar. When using up-to-date line lists, our [Mg/Fe] ratios also agree well with literature data. Our [Na/Fe] ratios are, on average, 0.08-0.14 dex lower than average values quoted in the literature, and our [Ba/Fe] ratios may be overestimated by 0.20-0.35 dex at the lowest metallicities. We find that analyses based on theoretical isochrones give very similar results to those based on resolved CMDs. Overall, the agreement between our integrated-light abundance measurements and the literature data is satisfactory. Refinements of the modelling procedure, such as corrections for stellar evolutionary and non-LTE effects, might further reduce some of the remaining offsets.


title: Polarization in Monte Carlo radiative transfer and dust scattering polarization signatures of spiral galaxies

abstract: Polarization is an important tool to further the understanding of interstellar dust and the sources behind it. In this paper we describe our implementation of polarization that is due to scattering of light by spherical grains and electrons in the dust Monte Carlo radiative transfer code SKIRT. In contrast to the implementations of other Monte Carlo radiative transfer codes, ours uses co-moving reference frames that rely solely on the scattering processes. It fully supports the peel-off mechanism that is crucial for the efficient calculation of images in 3D Monte Carlo codes. We develop reproducible test cases that push the limits of our code. The results of our program are validated by comparison with analytically calculated solutions. Additionally, we compare results of our code to previously published results. We apply our method to models of dusty spiral galaxies at near-infrared and optical wavelengths. We calculate polarization degree maps and show them to contain signatures that trace characteristics of the dust arms independent of the inclination or rotation of the galaxy.


title: Connecting the first galaxies with ultra faint dwarfs in the Local Group: chemical signatures of Population~III stars

abstract: We investigate the star formation history and chemical evolution of isolated analogues of Local Group (LG) ultra faint dwarf galaxies (UFDs; stellar mass range of 10^2 solar mass < M_star <10^5 solar mass) and gas rich, low mass dwarfs (Leo P analogs; stellar mass range of 10^5 solar mass < M_star <10^6 solar mass). We perform a suite of cosmological hydrodynamic zoom-in simulations to follow their evolution from the era of the first generation of stars down to z=0. We confirm that reionization, combined with supernova (SN) feedback, is primarily responsible for the truncated star formation in UFDs. Specifically, haloes with a virial mass of M_vir < 2 x 10^9 solar mass form> 90\% of stars prior to reionization. Our work further demonstrates the importance of Population~III (Pop~III) stars, with their intrinsically high $\rm [C/Fe]$ yields, and the associated external metal-enrichment, in producing low-metallicity stars ($\rm [Fe/H]\lesssim-4$) and carbon-enhanced metal-poor (CEMP) stars. We find that UFDs are composite systems, assembled from multiple progenitor haloes, some of which hosted only Population~II (Pop~II) stars formed in environments externally enriched by SNe in neighboring haloes, naturally producing, extremely low-metallicity Pop~II stars. We illustrate how the simulated chemical enrichment may be used to constrain the star formation histories (SFHs) of true observed UFDs. We find that Leo P analogs can form in haloes with M_vir ~ 4 x 10^9 solar mass (z=0). Such systems are less affected by reionization and continue to form stars until z=0, causing higher metallicity tails. Finally, we predict the existence of extremely low-metallicity stars in LG UFD galaxies that preserve the pure chemical signatures of Pop~III nucleosynthesis.


title: $\Upsilon$ and $\psi$ leptonic decays as probes of solutions to the $R_D^{(*)}$ puzzle

abstract: Experimental measurements of the ratios $R(D^{(*)})\equiv\frac{\Gamma(B\to D^{(*)}\tau\nu)}{\Gamma(B\to D^{(*)}\ell\nu)}$ ($\ell=e,\mu$) show a $3.9\sigma$ deviation from the Standard Model prediction. In the absence of light right-handed neutrinos, a new physics contribution to $b\to c\tau\nu$ decays necessarily modifies also $b\bar b\to\tau^+\tau^-$ and/or $c\bar c\to\tau^+\tau^-$ transitions. These contributions lead to violation of lepton flavor universality in, respectively, $\Upsilon$ and $\psi$ leptonic decays. We analyze the constraints resulting from measurements of the leptonic vector-meson decays on solutions to the $R(D^{(*)})$ puzzle. Available data from BaBar and Belle can already disfavor some of the new physics explanations of this anomaly. Further discrimination can be made by measuring $\Upsilon(1S,2S,3S)\to\tau\tau$ in the upcoming Belle II experiment.


title: AGN vs. host galaxy properties in the COSMOS field

abstract: The coeval AGN and galaxy evolution and the observed local relations between SMBHs and galaxy properties suggest some connection or feedback between SMBH growth and galaxy build-up. We looked for correlations between properties of X-ray detected AGN and their FIR detected host galaxies, to find quantitative evidences for this connection, highly debated in the latest years. We exploit the rich multi-wavelength data set available in the COSMOS field for a large sample (692 sources) of AGN and their hosts, in the redshift range $0.1<z<4$. We use X-ray data to select AGN and determine their properties (intrinsic luminosity and nuclear obscuration), and broad-band SED fitting to derive host galaxy properties (stellar mass $M_*$ and star formation rate SFR). We find that the AGN 2-10 keV luminosity ($L_{\rm X}$) and the host $8-1000~\mu m$ star formation luminosity ($L_{\rm IR}^{\rm SF}$) are significantly correlated. However, the average host $L_{\rm IR}^{\rm SF}$ has a flat distribution in bins of AGN $L_{\rm X}$, while the average AGN $L_{\rm X}$ increases in bins of host $L_{\rm IR}^{\rm SF}$, with logarithmic slope of $\sim0.7$, in the redshifts range $0.4<z<1.2$. We also discuss the comparison between the distribution of these two quantities and the predictions from hydro-dynamical simulations. Finally we find that the average column density ($N_H$) shows a positive correlation with the host $M_*$, at all redshifts, but not with the SFR (or $L_{\rm IR}^{\rm SF}$). This translates into a negative correlation with specific SFR. Our results are in agreement with the idea that BH accretion and SF rates are correlated, but occur with different variability time scales. The presence of a positive correlation between $N_H$ and host $M_*$ suggests that the X-ray $N_H$ is not entirely due to the circum-nuclear obscuring torus, but may also include a contribution from the host galaxy.


title: Optimal control of Rydberg lattice gases

abstract: We present optimal control protocols to prepare different many-body quantum states of Rydberg atoms in optical lattices. Specifically, we show how to prepare highly ordered many-body ground states, GHZ states as well as some superposition of symmetric excitation number Fock states, that inherit the translational symmetry from the Hamiltonian, within sufficiently short excitation times minimizing detrimental decoherence effects. For the GHZ states, we propose a two-step detection protocol to experimentally verify the optimal preparation of the target state based only on standard measurement techniques. Realistic experimental constraints and imperfections are taken into account by our optimization procedure making it applicable to ongoing experiments.


title: A light sneutrino rescues the light stop

abstract: Stop searches in supersymmetric frameworks with $R$-parity conservation usually assume the lightest neutralino to be the lightest supersymmetric particle. In this paper we consider an alternative scenario in which the left-handed tau sneutrino is lighter than neutralinos and stable at collider scales, but possibly unstable at cosmological scales. Moreover the (mostly right-handed) stop $\widetilde t$ is lighter than all electroweakinos, and heavier than the scalars of the third generation doublet, whose charged component, $\widetilde\tau$, is heavier than the neutral one, $\widetilde\nu$. The remaining supersymmetric particles are decoupled from the stop phenomenology. In most of the parameter space, the relevant stop decays are only into $t \widetilde\tau \tau$, $t\widetilde\nu\nu$ and $b \widetilde\nu \tau$ via off-shell electroweakinos. We constrain the branching ratios of these decays by recasting the most sensitive stop searches. Due to the "double invisible" kinematics of the $\widetilde t\to t\widetilde\nu\nu$ process, and the low efficiency in tagging the $t\widetilde\tau\tau$ decay products, light stops are generically allowed. In the minimal supersymmetric standard model with $\sim$ 100 GeV sneutrinos, stops with masses as small as $\sim$ 350 GeV turn out to be allowed at 95% CL.


title: Neural Decision Trees

abstract: In this paper we propose a synergistic melting of neural networks and decision trees (DT) we call neural decision trees (NDT). NDT is an architecture a la decision tree where each splitting node is an independent multilayer perceptron allowing oblique decision functions or arbritrary nonlinear decision function if more than one layer is used. This way, each MLP can be seen as a node of the tree. We then show that with the weight sharing asumption among those units, we end up with a Hashing Neural Network (HNN) which is a multilayer perceptron with sigmoid activation function for the last layer as opposed to the standard softmax. The output units then jointly represent the probability to be in a particular region. The proposed framework allows for global optimization as opposed to greedy in DT and differentiability w.r.t. all parameters and the input, allowing easy integration in any learnable pipeline, for example after CNNs for computer vision tasks. We also demonstrate the modeling power of HNN allowing to learn union of disjoint regions for final clustering or classification making it more general and powerful than standard softmax MLP requiring linear separability thus reducing the need on the inner layer to perform complex data transformations. We finally show experiments for supervised, semi-suppervised and unsupervised tasks and compare results with standard DTs and MLPs.


title: Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope using Sparse Modeling

abstract: We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the $\ell_1$-norm and total variation (TV) of the brightness distribution. In the proposed method, optimal regularization parameters, which represent the sparseness and effective spatial resolution of the image, are derived from data themselves using cross validation (CV). As an application of this technique, we present simulated observations of M87 with the EHT based on four physically motivated models. We confirm that $\ell_1$+TV regularization can achieve an optimal resolution of $\sim 20-30$% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer. With the proposed technique, the EHT can robustly and reasonably achieve super-resolution sufficient to clearly resolve the black hole shadow. These results make it promising for the EHT to provide an unprecedented view of the event-horizon-scale structure in the vicinity of the super-massive black hole in M87 and also the Galactic center Sgr A*.


title: In search of conformal theories

abstract: The conformal crossing equation puts very stringent constraints on the conformal data. We formulate it in way that makes the conformal symmetry more transparent. This allows for generalization of the crossing equation to arbitrary Lie group G. Using the crossing equation for SU(2) as a toy model, we find infinitely many solutions to the G-crossing equation. In particular, when G is specialized to the conformal group SO(d+1,1), we get infinitely many solutions to the conformal crossing equation.


title: On spectral sets of integers

abstract: Based on tiles and on the Coven-Meyerowitz property, we present some examples and some general constructions of spectral subsets of integers.


title: Nematic order on the surface of a three-dimensional topological insulator

abstract: We study the spontaneous breaking of rotational symmetry in the helical surface state of three-dimensional topological insulators due to strong electron-electron interactions, focusing on time-reversal invariant nematic order. Owing to the strongly spin-orbit coupled nature of the surface state, the nematic order parameter is linear in the electron momentum and necessarily involves the electron spin, in contrast with spin-degenerate nematic Fermi liquids. For a chemical potential at the Dirac point (zero doping), we find a first-order phase transition at zero temperature between isotropic and nematic Dirac semimetals. This extends to a thermal phase transition that changes from first to second order at a finite-temperature tricritical point. At finite doping, we find a transition between isotropic and nematic helical Fermi liquids that is second order even at zero temperature. Focusing on finite doping, we discuss various observable consequences of nematic order, such as anisotropies in transport and the spin susceptibility, the partial breakdown of spin-momentum locking, collective modes and induced spin fluctuations, and non-Fermi liquid behavior at the quantum critical point and in the nematic phase.


title: The primes are not metric Poissonian

abstract: It has been known since Vinogradov that, for irrational $\alpha$, the sequence of fractional parts $\{\alpha p\}$ is equidistributed in $\mathbb{R}/\mathbb{Z}$ as $p$ ranges over primes. There is a natural second-order equidistribution property, a pair correlation of such fractional parts, which has recently received renewed interest, in particular regarding its relation to additive combinatorics. In this paper we show that the primes do not enjoy this stronger equidistribution property.


title: A companion to "Knot invariants and M-theory I'' [arXiv:1608.05128]: proofs and derivations

abstract: We construct two distinct yet related M-theory models that provide suitable frameworks for the study of knot invariants. We then focus on the four-dimensional gauge theory that follows from appropriately compactifying one of these M-theory models. We show that this theory has indeed all required properties to host knots. Our analysis provides a unifying picture of the various recent works that attempt an understanding of knot invariants using techniques of four-dimensional physics. This is a companion paper to arXiv:1608.05128, covering all but section 3.3. It presents a detailed mathematical derivation of the main results there, as well as additional material. Among the new insights, those related to supersymmetry and the topological twist are highlighted. This paper offers an alternative, complementary formulation of the contents in~\cite{Dasgupta:2016rhc}, but is self-contained and can be read independently.


title: Stochastic Newton and Quasi-Newton Methods for Large Linear Least-squares Problems

abstract: We describe stochastic Newton and stochastic quasi-Newton approaches to efficiently solve large linear least-squares problems where the very large data sets present a significant computational burden (e.g., the size may exceed computer memory or data are collected in real-time). In our proposed framework, stochasticity is introduced in two different frameworks as a means to overcome these computational limitations, and probability distributions that can exploit structure and/or sparsity are considered. Theoretical results on consistency of the approximations for both the stochastic Newton and the stochastic quasi-Newton methods are provided. The results show, in particular, that stochastic Newton iterates, in contrast to stochastic quasi-Newton iterates, may not converge to the desired least-squares solution. Numerical examples, including an example from extreme learning machines, demonstrate the potential applications of these methods.


title: Reconfiguring motor circuits for a joint manual and BCI task

abstract: Designing brain-computer interfaces (BCIs) that can be used in conjunction with ongoing motor behavior requires an understanding of how neural activity co-opted for brain control interacts with existing neural circuits. For example, BCIs may be used to regain lost motor function after stroke. This requires that neural activity controlling unaffected limbs is dissociated from activity controlling the BCI. In this study we investigated how primary motor cortex accomplishes simultaneous BCI control and motor control in a task that explicitly required both activities to be driven from the same brain region (i.e. a dual-control task). Single-unit activity was recorded from intracortical, multi-electrode arrays while a non-human primate performed this dual-control task. Compared to activity observed during naturalistic motor control, we found that both units used to drive the BCI directly (control units) and units that did not directly control the BCI (non-control units) significantly changed their tuning to wrist torque. Using a measure of effective connectivity, we observed that control units decrease their connectivity. Through an analysis of variance we found that the intrinsic variability of the control units has a significant effect on task proficiency. When this variance is accounted for, motor cortical activity is flexible enough to perform novel BCI tasks that require active decoupling of natural associations to wrist motion. This study provides insight into the neural activity that enables a dual-control brain-computer interface.


title: Bach-flat isotropic gradient Ricci solitons

abstract: We construct examples of Bach-flat gradient Ricci solitons which are neither half conformally flat nor conformally Einstein.


title: sdA in SDSS DR12 are Overwhelmingly Not Extremely Low-Mass (ELM) White Dwarfs

abstract: In a search for new white dwarfs in DR12 of the Sloan Digital Sky Survey, Kepler et al. (2016) found atmospheric parameters for thousands of objects with effective temperatures below 20,000 K and surface gravities between 5.5 < log(g) < 6.5. They classified these objects as cool subdwarfs -- sdA -- and speculated that many may be extremely low-mass (ELM) white dwarfs (helium-core white dwarfs with masses below 0.3 Msun). We present evidence -- using radial velocities, photometric colors, and reduced proper motions -- that the vast majority (>99%) of these objects are unlikely to be ELM white dwarfs. Their true identity remains an interesting question.


title: Feasibility of Principal Component Analysis in hand gesture recognition system

abstract: Nowadays actions are increasingly being handled in electronic ways, instead of physical interaction. From earlier times biometrics is used in the authentication of a person. It recognizes a person by using a human trait associated with it like eyes (by calculating the distance between the eyes) and using hand gestures, fingerprint detection, face detection etc. Advantages of using these traits for identification are that they uniquely identify a person and cannot be forgotten or lost. These are unique features of a human being which are being used widely to make the human life simpler. Hand gesture recognition system is a powerful tool that supports efficient interaction between the user and the computer. The main moto of hand gesture recognition research is to create a system which can recognise specific hand gestures and use them to convey useful information for device control. This paper presents an experimental study over the feasibility of principal component analysis in hand gesture recognition system. PCA is a powerful tool for analyzing data. The primary goal of PCA is dimensionality reduction. Frames are extracted from the Sheffield KInect Gesture (SKIG) dataset. The implementation is done by creating a training set and then training the recognizer. It uses Eigen space by processing the eigenvalues and eigenvectors of the images in training set. Euclidean distance with the threshold value is used as similarity metric to recognize the gestures. The experimental results show that PCA is feasible to be used for hand gesture recognition system.


title: Surfing on protein waves: proteophoresis as a mechanism for bacterial genome partitioning

abstract: Efficient bacterial chromosome segregation typically requires the coordinated action of a three-components, fueled by adenosine triphosphate machinery called the partition complex. We present a phenomenological model accounting for the dynamic activity of this system. The model is obtained by coupling simple linear reaction-diffusion equations with a proteophoresis, or "volumetric" chemophoresis, force field. This minimal description, in the sense of Occam's Razor principle, captures most known experimental observations: dynamic oscillations of complex components, complex separation and subsequent symmetrical positioning. The predictions of our model are in phenomenological agreement with and provide substantial insights into recent experiments. From a non-linear physics view point, this system explores the active separation of matter at micrometric scales with a dynamical instability between static positioning and travelling wave regimes.


title: High order surface radiation conditions for time-harmonic waves in exterior domains

abstract: We formulate a new family of high order on-surface radiation conditions to approximate the outgoing solution to the Helmholtz equation in exterior domains. Motivated by the pseudo-differential expansion of the Dirichlet-to-Neumann operator developed by Antoine et al. (J. Math. Anal. Appl. 229:184-211, 1999), we design a systematic procedure to apply pseudo-differential symbols of arbitrarily high order. Numerical results are presented to illustrate the performance of the proposed method for solving both the Dirichlet and the Neumann boundary value problems. Possible improvements and extensions are also discussed.


title: Time series momentum and contrarian effects in the Chinese stock market

abstract: This paper concentrates on the time series momentum or contrarian effects in the Chinese stock market. We evaluate the performance of the time series momentum strategy applied to major stock indices in mainland China and explore the relation between the performance of time series momentum strategies and some firm-specific characteristics. Our findings indicate that there is a time series momentum effect in the short run and a contrarian effect in the long run in the Chinese stock market. The performances of the time series momentum and contrarian strategies are highly dependent on the look-back and holding periods and firm-specific characteristics.


title: Delta-net: Real-time Network Verification Using Atoms

abstract: Real-time network verification promises to automatically detect violations of network-wide reachability invariants on the data plane. To be useful in practice, these violations need to be detected in the order of milliseconds, without raising false alarms. To date, most real-time data plane checkers address this problem by exploiting at least one of the following two observations: (i) only small parts of the network tend to be affected by typical changes to the data plane, and (ii) many different packets tend to share the same forwarding behaviour in the entire network. This paper shows how to effectively exploit a third characteristic of the problem, namely: similarity among forwarding behaviour of packets through parts of the network, rather than its entirety. We propose the first provably amortized quasi-linear algorithm to do so. We implement our algorithm in a new real-time data plane checker, Delta-net. Our experiments with SDN-IP, a globally deployed ONOS software-defined networking application, and several hundred million IP prefix rules generated using topologies and BGP updates from real-world deployed networks, show that Delta-net checks a rule insertion or removal in approximately 40 microseconds on average, a more than 10X improvement over the state-of-the-art. We also show that Delta-net eliminates an inherent bottleneck in the state-of-the-art that restricts its use in answering Datalog-style "what if" queries.


title: Tides in a body librating about a spin-orbit resonance. Generalisation of the Darwin-Kaula theory

abstract: The Darwin-Kaula theory of bodily tides is intended for celestial bodies rotating without libration. We demonstrate that this theory, in its customary form, is inapplicable to a librating body. Specifically, in the presence of libration in longitude, the actual spectrum of Fourier tidal modes differs from the conventional spectrum rendered by the Darwin-Kaula theory for a non-librating celestial object. This necessitates derivation of formulae for the tidal torque and the tidal heating rate, that are applicable under libration. We derive the tidal spectrum for longitudinal forced libration with one and two main frequencies, generalisation to more main frequencies being straightforward. (By main frequencies we understand those emerging due to the triaxiality of the librating body.) Separately, we consider a case of free libration at one frequency (once again, generalisation to more frequencies being straightforward). We also calculate the tidal torque. This torque provides correction to the triaxiality-caused physical libration. Our theory is not self-consistent: we assume that the tidal torque is much smaller than the permanent-triaxiality-caused torque; so the additional libration due to tides is much weaker than the main libration due to the permanent triaxiality. Finally, we calculate the tidal dissipation rate in a body experiencing forced libration at the main mode, or free libration at one frequency, or superimposed forced and free librations.


title: Neutral gas heating by X-rays in primitive galaxies: Infrared observations of IZw18 with Herschel

abstract: (abridged) The dominant thermal mechanisms in the neutral interstellar medium, which acts as a star-forming gas reservoir, are uncertain in extremely metal-poor galaxies. Our objective is to identify the heating mechanisms in one such galaxy, IZw18, and assess the diagnostic value of fine-structure cooling lines. We also seek to constrain the mass of H$_2$, which, despite being an important catalyst and tracer of star formation, remains elusive in this object. Building on a previous photoionization model within a multi-sector topology, we provide additional constraints from the [CII] and [OI] lines and the dust mass recently measured with Herschel. The heating of the HI region appears to be mainly due to photoionization by radiation from a bright X-ray binary source, while photoelectric effect (PE) is negligible. The [CII] and [OI] lines imply an average X-ray luminosity of $4\times10^{40}$ erg s$^{-1}$, while the [NeV] upper limits bring strong constraints to the soft X-ray flux arising from the binary. A negligible amount of H$_2$ is predicted, but $\lesssim10^7$ M$_\odot$ of H$_2$ may be hidden in sufficiently dense clouds of order $\lesssim10$ pc in size. Regardless of the presence of significant amounts of H$_2$, [CII] and [OI] do not trace the so-called CO-dark gas, but the almost purely atomic medium. Although the [CII]+[OI]/TIR ratio is close to values found in more metal-rich sources, it cannot be safely used as a PE heating efficiency proxy. This ratio seems to be kept stable due to a correlation between the X-ray luminosity and the star-formation rate. We propose that X-ray heating is an important process in extremely metal-poor sources. The weak PE heating due to the low dust-to-gas ratio could be compensated for by the larger occurrence and power of X-ray binaries in low-metallicity galaxies. We speculate that X-ray heating may quench star formation.


title: Selection of the tagged photons by off axis heterodyne holography in ultrasound-modulated optical tomography

abstract: Ultrasound-modulated optical tomography (UOT) is a technique that images optical contrast deep inside scattering media. Heterodyne holography is a promising tool able to detect the UOT tagged photons with high efficiency. In this work, we describe theoretically the detection of the tagged photon in heterodyne holography based UOT, show how to filter the untagged photon, and discuss the effect of shot noise. The discussion considers also speckle decorrelation. We show that optimal detection sensitivity can be reached, if the frame exposure time of the camera used to perform the holographic detection is of the order of the decorrelation time.


title: A Complete Characterization of the 1-Dimensional Intrinsic Cech Persistence Diagrams for Metric Graphs

abstract: Metric graphs are special types of metric spaces used to model and represent simple, ubiquitous, geometric relations in data such as biological networks, social networks, and road networks. We are interested in giving a qualitative description of metric graphs using topological summaries. In particular, we provide a complete characterization of the 1-dimensional intrinsic Cech persistence diagrams for metric graphs using persistent homology. Together with complementary results by Adamaszek et. al, which imply results on intrinsic Cech persistence diagrams in all dimensions for a single cycle, our results constitute important steps toward characterizing intrinsic Cech persistence diagrams for arbitrary metric graphs across all dimensions.


title: Consequences of energy conservation violation: Late time solutions of $\Lambda(T) CDM$ subclass of $f(R,T)$ gravity using dynamical system approach

abstract: Very recently, the authors of [PRL {\bf 118} (2017) 021102] have shown that violation of energy-momentum tensor (EMT) could result in an accelerated expansion state via appearing an effective cosmological constant, in the context of unimodular gravity. Inspired by this outcome, in this paper we investigate cosmological consequences of violation of the EMT conservation in a particular class of $f(R,T)$ gravity when only the pressure-less fluid is present. In this respect, we focus on the late time solutions of models of the type $f(R,T)=R+\beta \Lambda(-T)$. As the first task, we study the solutions when the conservation of EMT is respected and then we proceed with those in which violation occurs. We have found, provided that the EMT conservation is violated, there generally exist two accelerated expansion solutions which their stability properties depend on the underlying model. More exactly, we obtain a dark energy solution for which the effective equation of state (EoS) depend on model parameters and a de Sitter solution. We present a method to parametrize $\Lambda(-T)$ function which is useful in dynamical system approach and has been employed in the herein model. Also, we discuss the cosmological solutions for models with $\Lambda(-T)=8\pi G(-T)^{\alpha}$ in the presence of the ultra relativistic matter.


title: Kondo-Ising and Tight-Binding Models for TmB4

abstract: In $TmB_4$, localized electrons with a large magnetic moment interact with metallic electrons in boron-derived bands. We examine the nature of $TmB_4$ using full-relativistic ab-initio density functional theory calculations, approximate tight-binding Hamiltonian results, and the development of an effective Kondo-Ising model for this system. Features of the Fermi surface relating to the anisotropic conduction of charge are discussed. The observed magnetic moment $\sim 6 \, \mu_B$ is argued to require a subtle crystal field effect in metallic systems, involving a flipped sign of the effective charges surrounding a Tm ion. The role of on-site quantum dynamics in the resulting Kondo-Ising type "impurity" model are highlighted. From this model, elimination of the conduction electrons will lead to spin-spin (RKKY-type) interaction of Ising character required to understand the observed fractional magnetization plateaus in $TmB_4$.


title: The Relation Between Fundamental Constants and Particle Physics Parameters

abstract: The observed constraints on the variability of the proton to electron mass ratio $\mu$ and the fine structure constant $\alpha$ are used to establish constraints on the variability of the Quantum Chromodynamic Scale and a combination of the Higgs Vacuum Expectation Value and the Yukawa couplings. Further model dependent assumptions provide constraints on the Higgs VEV and the Yukawa couplings separately. A primary conclusion is that limits on the variability of dimensionless fundamental constants such as $\mu$ and $\alpha$ provide important constraints on the parameter space of new physics and cosmologies.


title: Fulde-Ferrell States in Inverse Proximity Coupled Magnetically-Doped Topological Heterostructures

abstract: We study the superconducting properties of the thin film BCS superconductor proximity coupled to a magnetically doped topological insulator(TI). Using the mean field theory, we show that Fulde-Ferrell(FF) pairing can be induced in the conventional superconductor by having inverse proximity effect(IPE). This occurs when the IPE of the TI to the superconductor is large enough that the normal band of the superconductor possesses a proximity induced spin-orbit coupling and magnetization. We find that the energetics of the different pairings are dependent on the coupling strength between the TI and the BCS superconductor and the thickness of the superconductor film. As the thickness of the superconductor film is increased, we find a crossover from the FF pairing to the BCS pairing phase. This is a consequence of the increased number of the superconducting bands, which favor the BCS pairing, implying that the FF phase can only be observed in the thin-film limit. In addition, we also propose transport experiments that show distinct signatures of the FF phase.


title: Quantum information measures and their applications in quantum field theory

abstract: In the last decades, it has been understood that a wide variety of phenomena in quantum field theory (QFT) can be characterised using quantum information measures, such as the entanglement entropy of a state and the relative entropy between quantum states in the same Hilbert space. In this thesis, we use these and other tools from quantum information theory to study several interesting problems in quantum field theory. The topics analysed range from the study of the Aharonov-Bohm effect in QFT using entanglement entropy, to the consistence of the Ryu-Takayanagi formula (proposed in the context of the AdS/CFT duality) using properties of relative entropy. We show that relative entropy can also be used to obtain new interesting quantum energy inequalities, that constrain the spatial distribution of negative energy density.


title: From inflation to recent cosmic acceleration: The fermionic Elko field driving the evolution of the universe

abstract: In this paper we construct the complete evolution of the universe driven by the mass dimension one dark spinor called Elko, starting with inflation, passing by the matter dominated era and finishing with the recent accelerated expansion. The dynamic of the fermionic Elko field with a symmetry breaking type potential can reproduce all phases of the universe in a natural and elegant way. The dynamical equations in general case and slow roll conditions in the limit $H\ll m_{pl}$ are also presented for the Elko system. Numerical analysis for the number of e-foldings during inflation, energy density after inflation and for present time and also the actual size of the universe are in good agreement with the standard model of cosmology. An interpretation of the inflationary phase as a result of Pauli exclusion principle is also possible if the Elko field is treated as an average value of its quantum analogue.


title: Toward Streaming Synapse Detection with Compositional ConvNets

abstract: Connectomics is an emerging field in neuroscience that aims to reconstruct the 3-dimensional morphology of neurons from electron microscopy (EM) images. Recent studies have successfully demonstrated the use of convolutional neural networks (ConvNets) for segmenting cell membranes to individuate neurons. However, there has been comparatively little success in high-throughput identification of the intercellular synaptic connections required for deriving connectivity graphs. In this study, we take a compositional approach to segmenting synapses, modeling them explicitly as an intercellular cleft co-located with an asymmetric vesicle density along a cell membrane. Instead of requiring a deep network to learn all natural combinations of this compositionality, we train lighter networks to model the simpler marginal distributions of membranes, clefts and vesicles from just 100 electron microscopy samples. These feature maps are then combined with simple rules-based heuristics derived from prior biological knowledge. Our approach to synapse detection is both more accurate than previous state-of-the-art (7% higher recall and 5% higher F1-score) and yields a 20-fold speed-up compared to the previous fastest implementations. We demonstrate by reconstructing the first complete, directed connectome from the largest available anisotropic microscopy dataset (245 GB) of mouse somatosensory cortex (S1) in just 9.7 hours on a single shared-memory CPU system. We believe that this work marks an important step toward the goal of a microscope-pace streaming connectomics pipeline.


title: On the Entropy of Deformed Phase Space Black Hole and the Cosmological Constant

abstract: In this paper we study the effects of noncommutative phase space deformations on the Schwarzschild black hole. This idea has been previously studied in {Friedmann-Robertson-Walker (FRW)} cosmology, where this "noncommutativity" provides a simple mechanism that can explain the origin of the cosmological constant. In this paper we obtain the same relationship between the cosmological constant and the deformation parameter that appears in deformed phase space cosmology, but in the context of the deformed phase space black holes. This was achieved by comparing the entropy of the deformed Schwarzschild black hole with the entropy of the Schwarzschild-de Sitter black hole.


title: On Store Languages of Language Acceptors

abstract: It is well known that the "store language" of every pushdown automaton -- the set of store configurations (state and stack contents) that can appear as an intermediate step in accepting computations -- is a regular language. Here many models of language acceptors with various data structures are examined, along with a study of their store languages. For each model, an attempt is made to find the simplest model that accepts their store languages. Some connections between store languages of one-way and two-way machines generally are demonstrated, as with connections between nondeterministic and deterministic machines. A nice application of these store language results is also presented, showing a general technique for proving families accepted by many deterministic models are closed under right quotient with regular languages, resolving some open questions (and significantly simplifying proofs for others that are known) in the literature. Lower bounds on the space complexity for recognizing store languages for the languages to be non-regular are obtained.


title: Continuous-Time Visual-Inertial Odometry for Event Cameras

abstract: Event cameras are bio-inspired vision sensors that output pixel-level brightness changes instead of standard intensity frames. They offer significant advantages over standard cameras, namely a very high dynamic range, no motion blur, and a latency in the order of microseconds. However, due to the fundamentally different structure of the sensor's output, new algorithms that exploit the high temporal resolution and the asynchronous nature of the sensor are required. Recent work has shown that a continuous-time representation of the event camera pose can deal with the high temporal resolution and asynchronous nature of this sensor in a principled way. In this paper, we leverage such a continuous-time representation to perform visual-inertial odometry with an event camera. This representation allows direct integration of the asynchronous events with micro-second accuracy and the inertial measurements at high frequency. The event camera trajectory is approximated by a smooth curve in the space of rigid-body motions using cubic splines. This formulation significantly reduces the number of variables in trajectory estimation problems. We evaluate our method on real data from several scenes and compare the results against ground truth from a motion-capture system. We show that our method provides improved accuracy over the result of a state-of-the-art visual odometry method for event cameras. We also show that both the map orientation and scale can be recovered accurately by fusing events and inertial data. To the best of our knowledge, this is the first work on visual-inertial fusion with event cameras using a continuous-time framework.


title: Detecting Strong Ties Using Network Motifs

abstract: Detecting strong ties among users in social and information networks is a fundamental operation that can improve performance on a multitude of personalization and ranking tasks. Strong-tie edges are often readily obtained from the social network as users often participate in multiple overlapping networks via features such as following and messaging. These networks may vary greatly in size, density and the information they carry. This setting leads to a natural strong tie detection task: given a small set of labeled strong tie edges, how well can one detect unlabeled strong ties in the remainder of the network? This task becomes particularly daunting for the Twitter network due to scant availability of pairwise relationship attribute data, and sparsity of strong tie networks such as phone contacts. Given these challenges, a natural approach is to instead use structural network features for the task, produced by {\em combining} the strong and "weak" edges. In this work, we demonstrate via experiments on Twitter data that using only such structural network features is sufficient for detecting strong ties with high precision. These structural network features are obtained from the presence and frequency of small network motifs on combined strong and weak ties. We observe that using motifs larger than triads alleviate sparsity problems that arise for smaller motifs, both due to increased combinatorial possibilities as well as benefiting strongly from searching beyond the ego network. Empirically, we observe that not all motifs are equally useful, and need to be carefully constructed from the combined edges in order to be effective for strong tie detection. Finally, we reinforce our experimental findings with providing theoretical justification that suggests why incorporating these larger sized motifs as features could lead to increased performance in planted graph models.


title: Free-Space Entangled Quantum Carpets

abstract: The Talbot effect in quantum physics is known to produce intricate patterns in the probability distribution of a particle, known as "quantum carpets", corresponding to the revival and replication of the initial wave function. Recently, it was shown that one can encode a $D$-level qudit, in such a way that the Talbot effect can be used to process the $D$-dimensional quantum information [Far\'{\i}as et al, PRA (2015)]. Here we introduce a scheme to produce free-propagating "entangled quantum carpets" with pairs of photons produced by spontaneous parametric down-conversion. First we introduce an optical device that can be used to synthesize arbitrary superposition states of Talbot qudits. Sending spatially entangled photon pairs through a pair of these devices produces an entangled pair of qudits. As an application, we show how the Talbot effect can be used to test a $D$-dimensional Bell inequality. Numerical simulations show that violation of the Bell inequality depends strongly on the amount of spatial correlation in the initial two-photon state. We briefly discuss how our optical scheme might be adapted to matter wave experiments.


title: WaterGAN: Unsupervised Generative Network to Enable Real-time Color Correction of Monocular Underwater Images

abstract: This paper reports on WaterGAN, a generative adversarial network (GAN) for generating realistic underwater images from in-air image and depth pairings in an unsupervised pipeline used for color correction of monocular underwater images. Cameras onboard autonomous and remotely operated vehicles can capture high resolution images to map the seafloor, however, underwater image formation is subject to the complex process of light propagation through the water column. The raw images retrieved are characteristically different than images taken in air due to effects such as absorption and scattering, which cause attenuation of light at different rates for different wavelengths. While this physical process is well described theoretically, the model depends on many parameters intrinsic to the water column as well as the objects in the scene. These factors make recovery of these parameters difficult without simplifying assumptions or field calibration, hence, restoration of underwater images is a non-trivial problem. Deep learning has demonstrated great success in modeling complex nonlinear systems but requires a large amount of training data, which is difficult to compile in deep sea environments. Using WaterGAN, we generate a large training dataset of paired imagery, both raw underwater and true color in-air, as well as depth data. This data serves as input to a novel end-to-end network for color correction of monocular underwater images. Due to the depth-dependent water column effects inherent to underwater environments, we show that our end-to-end network implicitly learns a coarse depth estimate of the underwater scene from monocular underwater images. Our proposed pipeline is validated with testing on real data collected from both a pure water tank and from underwater surveys in field testing. Source code is made publicly available with sample datasets and pretrained models.


title: Controlling Parent Systems Through Swarms Using Abstraction

abstract: This study considers the control of parent-child systems where a parent system is acted on by a set of controllable child systems (i.e. a swarm). Examples of such systems include a swarm of robots pushing an object over a surface, a swarm of aerial vehicles carrying a large load, or a set of end effectors manipulating an object. In this paper, a general approach for decoupling the swarm from the parent system through a low-dimensional abstract state space is presented. The requirements of this approach are given along with how constraints on both systems propagate through the abstract state and impact the requirements of the controllers for both systems. To demonstrate, several controllers with hard state constraints are designed to track a given desired angle trajectory of a tilting plane with a swarm of robots driving on top. Both homogeneous and heterogeneous swarms of varying sizes and properties are considered to test the robustness of this architecture. The controllers are shown to be locally asymptotically stable and are demonstrated in simulation.


title: First passage properties in a crowded environment

abstract: We develop a model to compute the first-passage time of a random walker in a crowded environment. Hard-core particles with the same size and diffusion coefficient than the tracer diffuse, and the model allows to compute the first passage time of the tracer on euclidian lattices. The result is compared to classical Nakazato-Kitahara model, and extends previous results obtained for persistent random walker. The crowding in a confined media acts as a memory effect, and thus lead to a persistent-like behavior.


title: Boundary-Driven Anomalous Spirals in Oscillatory Media

abstract: We study a heretofore ignored class of spiral patterns for oscillatory media as characterized by the complex Landau-Ginzburg model. These spirals emerge from modulating the growth rate as a function of $r$, thereby turning off the instability. These spirals are uniquely determined by matching to those outer conditions, lifting a degeneracy in the set of steady-state solutions of the original equations. Unlike the well-studied spiral which acts a wave source, has a simple core structure and is insensitive to the details of the boundary on which no-flux conditions are imposed, these new spirals are wave sinks, have non-monotonic wavefront curvature near the core, and can be patterned by the form of the spatial boundary. We predict that these anomalous spirals could be produced in nonlinear optics experiments via spatially modulating the gain of the medium.


title: Hunt's Hypothesis (H) for the Sum of Two Independent Levy Processes

abstract: Which Levy processes satisfy Hunt's hypothesis (H) is a long-standing open problem in probabilistic potential theory. The study of this problem for one-dimensional Levy processes suggests us to consider (H) from the point of view of the sum of Levy processes. In this paper, we present theorems and examples on the validity of (H) for the sum of two independent Levy processes. We also give a novel condition on the Levy measure which implies (H) for a large class of one-dimensional Levy processes.


title: Singular FIOs in SAR Imaging, II: Transmitter and Receiver at Different Speeds

abstract: In this article, we consider two bistatic cases arising in synthetic aperture radar imaging: when the transmitter and receiver are both moving with different speeds along a single line parallel to the ground in the same direction or in the opposite directions. In both cases, we classify the forward operator $\Fc$ as a Fourier integral operator with fold/blowdown singularities. Next we analyze the normal operator $\Fc^*\Fc$ in both cases (where $\Fc^{*}$ is the $L^{2}$-adjoint of $\Fc$). When the transmitter and receiver move in the same direction, we prove that $\Fc^*\Fc$ belongs to a class of operators associated to two cleanly intersecting Lagrangians, $I^{p,l} (\Delta, C_1)$. When they move in opposite directions, $\Fc^*\Fc$ is a sum of such operators. In both cases artifacts appear and we show that they are, in general, as strong as the bona-fide part of the image. Moreover, we demonstrate that as soon as the source and receiver start to move in opposite directions, there is an interesting bifurcation in the type of artifact that appears in the image.


title: Deep Nonparametric Estimation of Discrete Conditional Distributions via Smoothed Dyadic Partitioning

abstract: We present an approach to deep estimation of discrete conditional probability distributions. Such models have several applications, including generative modeling of audio, image, and video data. Our approach combines two main techniques: dyadic partitioning and graph-based smoothing of the discrete space. By recursively decomposing each dimension into a series of binary splits and smoothing over the resulting distribution using graph-based trend filtering, we impose a strict structure to the model and achieve much higher sample efficiency. We demonstrate the advantages of our model through a series of benchmarks on both synthetic and real-world datasets, in some cases reducing the error by nearly half in comparison to other popular methods in the literature. All of our models are implemented in Tensorflow and publicly available at https://github.com/tansey/sdp .


title: An Optimal Algorithm for Computing the Spherical Depth of Points in the Plane

abstract: For a distribution function $F$ on $\mathbb{R}^d$ and a point $q\in \mathbb{R}^d$, the \emph{spherical depth} $\SphD(q;F)$ is defined to be the probability that a point $q$ is contained inside a random closed hyper-ball obtained from a pair of points from $F$. The spherical depth $\SphD(q;S)$ is also defined for an arbitrary data set $S\subseteq \mathbb{R}^d$ and $q\in \mathbb{R}^d$. This definition is based on counting all of the closed hyper-balls, obtained from pairs of points in $S$, that contain $q$. The significant advantage of using the spherical depth in multivariate data analysis is related to its complexity of computation. Unlike most other data depths, the time complexity of the spherical depth grows linearly rather than exponentially in the dimension $d$. The straightforward algorithm for computing the spherical depth in dimension $d$ takes $O(dn^2)$. The main result of this paper is an optimal algorithm that we present for computing the bivariate spherical depth. The algorithm takes $O(n \log n)$ time. By reducing the problem of \textit{Element Uniqueness}, we prove that computing the spherical depth requires $\Omega(n \log n)$ time. Some geometric properties of spherical depth are also investigated in this paper. These properties indicate that \emph{simplicial depth} ($\SD$) (Liu, 1990) is linearly bounded by spherical depth (in particular, $\SphD\geq \frac{2}{3}SD$). To illustrate this relationship between the spherical depth and the simplicial depth, some experimental results are provided. The obtained experimental bound ($\SphD\geq 2\SD$) indicates that, perhaps, a stronger theoretical bound can be achieved.


title: Horseshoe Regularization for Feature Subset Selection

abstract: Feature subset selection arises in many high-dimensional applications of statistics, such as compressed sensing and genomics. The $\ell_0$ penalty is ideal for this task, the caveat being it requires the NP-hard combinatorial evaluation of all models. A recent area of considerable interest is to develop efficient algorithms to fit models with a non-convex $\ell_\gamma$ penalty for $\gamma\in (0,1)$, which results in sparser models than the convex $\ell_1$ or lasso penalty, but is harder to fit. We propose an alternative, termed the horseshoe regularization penalty for feature subset selection, and demonstrate its theoretical and computational advantages. The distinguishing feature from existing non-convex optimization approaches is a full probabilistic representation of the penalty as the negative of the logarithm of a suitable prior, which in turn enables efficient expectation-maximization and local linear approximation algorithms for optimization and MCMC for uncertainty quantification. In synthetic and real data, the resulting algorithms provide better statistical performance, and the computation requires a fraction of time of state-of-the-art non-convex solvers.


title: An Approximation Method for the Scattering Data of One-Dimensional Soliton Equations under Arbitrary Rapidly-Decreasing Initial Pulses

abstract: We present a novel approximation method that can predict the number of solitons asymptotically appearing under arbitrary rapidly decreasing initial wave packets. The number of solitons can be estimated without integration of the original soliton equations. As an example, we take the one-dimensional nonlinear Schrodinger equation and estimate the behaviors of the scattering amplitude in detail. The results show good agreement compared with those obtained by direct numerical integration. The presented method is applicable to a wide class of one-dimensional soliton equations.


title: Weak Faddeev-Takhtajan-Volkov algebras; Lattice $W_n$ algebras

abstract: In this paper, we will start by looking through our project's historical general view and then we will try to construct a new Poisson bracket on our simplest example $sl_2$ and then we will try to give a universal construction based on our universal variables and then will try to construct lattice $W_2$ algebras which will play a key role in our other constructions on lattice $W_3$ algebras and finally we will try to find the only nontrivial dependent generator of our lattice $W_4$ algebras and so on for lattice $W_n$ algebras. And at the end of this paper, we will have appendix A, which will contain some parts of the Mathematica coding which we have used and have made for to find our algebra structures.


title: Making Asynchronous Distributed Computations Robust to Noise

abstract: We consider the problem of making distributed computations robust to noise, in particular to worst-case (adversarial) corruptions of messages. We give a general distributed interactive coding scheme which simulates any asynchronous distributed protocol while tolerating an optimal corruption of a $\Theta(1/n)$ fraction of all messages while incurring a moderate blowup of $O(n\log^2 n)$ in the communication complexity. Our result is the first fully distributed interactive coding scheme in which the topology of the communication network is not known in advance. Prior work required either a coordinating node to be connected to all other nodes in the network or assumed a synchronous network in which all nodes already know the complete topology of the network.


title: Renormalizability in $D$-dimensional higher-order gravity

abstract: A simple expression for calculating the classical potential concerning $D$-dimensional gravitational models is obtained through a method based on the generating functional. The prescription is then used as a mathematical tool to probe the conjecture that renormalizable higher-order gravity models --- which are, of course, nonunitary --- are endowed with a classical potential that is nonsingular at the origin. It is also shown that the converse of this statement is not true, which implies that the finiteness of the classical potential at the origin is a necessary but not a sufficient condition for the renormalizability of the model. The systems we have utilized to verify the conjecture were fourth- and sixth- order gravity models in $D$-dimensions. A discussion about the polemic question related to the renormalizability of new massive gravity, which Oda claimed to be renormalizable in 2009 and three years late was shown to be nonrenormalizable by Muneyuki and Ohta, is considered. We remark that the solution of this issue is straightforward if the aforementioned conjecture is employed. We point out that our analysis is restricted to local models in which the propagator has simple and real poles.


title: GapTV: Accurate and Interpretable Low-Dimensional Regression and Classification

abstract: We consider the problem of estimating a regression function in the common situation where the number of features is small, where interpretability of the model is a high priority, and where simple linear or additive models fail to provide adequate performance. To address this problem, we present GapTV, an approach that is conceptually related both to CART and to the more recent CRISP algorithm, a state-of-the-art alternative method for interpretable nonlinear regression. GapTV divides the feature space into blocks of constant value and fits the value of all blocks jointly via a convex optimization routine. Our method is fully data-adaptive, in that it incorporates highly robust routines for tuning all hyperparameters automatically. We compare our approach against CART and CRISP and demonstrate that GapTV finds a much better trade-off between accuracy and interpretability.


title: On permutations of order dividing a given integer

abstract: We give a detailed analysis of the proportion of elements in the symmetric group on $n$ points whose order divides $m$, for $n$ sufficiently large and $m \ge n$ with $m = O(n)$.


title: Binary quartic forms with bounded invariants and small Galois groups

abstract: In this paper, we consider integral and irreducible binary quartic forms whose Galois group is isomorphic to a subgroup of the dihedral group of order eight. We first show that the set of all such forms is a union of families indexed by integral binary quadratic forms $f(x,y)$ of non-zero discriminant. Then, we shall enumerate the $\operatorname{GL}_2(\mathbb{Z})$-equivalence classes of all such forms associated to a fixed $f(x,y)$.


title: Control methods for improved Fisher information with quantum sensing

abstract: Recently new approaches for sensing the frequency of time dependent Hamiltonians have been presented, and it was shown that the optimal Fisher information scales as $T^{4}.$ We present here our interpretation of this new scaling, where the relative phase is accumulated quadratically with time, and show that this can be produced by a variety of simple pulse sequences. Interestingly, this scaling has a limited duration, and we show that certain pulse sequences prolong the effect. The performance of these schemes is analyzed and we examine their relevance to state-of-the-art experiments. We analyze the $T^{3}$ scaling of the Fisher information which appears when multiple synchronized measurements are performed, and is the optimal scaling in the case of a finite coherence time.


title: Founsure 1.0: An Erasure Code Library with Efficient Repair and Update Features

abstract: Founsure is an open-source software library that implements a multi-dimensional graph-based erasure coding entirely based on fast exclusive OR (XOR) logic. Its implementation utilizes compiler optimizations and multi-threading to generate the right assembly code for the given multi-core CPU architecture with vector processing capabilities. Founsure possesses important features that shall find various applications in modern data storage, communication, and networked computer systems, in which the data needs protection against device, hardware, and node failures. As data size reached unprecedented levels, these systems have become hungry for network bandwidth, computational resources, and average consumed power. To address that, the proposed library provides a three-dimensional design space that trades off the computational complexity, coding overhead, and data/node repair bandwidth to meet different requirements of modern distributed data storage and processing systems. Founsure library enables efficient encoding, decoding, repairs/rebuilds, and updates while all the required data storage and computations are distributed across the network nodes.


title: Adiabatic corrections for velocity-gauge simulations of electron dynamics in periodic potentials

abstract: We show how to significantly reduce the number of energy bands required to model the interaction of light with crystalline solids in the velocity gauge. We achieve this by deriving analytical corrections to the electric current density. These corrections depend only on band energies, the matrix elements of the momentum operator, and the macroscopic vector potential. Thus, the corrections can be evaluated independently from modeling the interaction with light. In addition to improving the convergence of velocity-gauge calculations, our analytical approach overcomes the long-standing problem of divergences in expressions for linear and nonlinear susceptibilities.


title: Diameter Rigidity for K\"ahler manifolds with positive bisectional curvature

abstract: Let $M^n$ be a compact K\"ahler manifold with bisectional curvature bounded from below by $1$. If $diam(M) = \pi / \sqrt{2}$ and $vol(M)> vol(\mathbb{C}\mathbb{P}^n)/ 2^n$, we prove that $M$ is biholomorphically isometric to $\mathbb{C}\mathbb{P}^n$ with the standard Fubini-Study metric.


title: Continuation of homoclinic orbits in the suspension bridge equation: a computer-assisted proof

abstract: In this paper, we prove existence of symmetric homoclinic orbits for the suspension bridge equation $u""+\beta u" + e^u-1=0$ for all parameter values $\beta \in [0.5,1.9]$. For each $\beta$, a parameterization of the stable manifold is computed and the symmetric homoclinic orbits are obtained by solving a projected boundary value problem using Chebyshev series. The proof is computer-assisted and combines the uniform contraction theorem and the radii polynomial approach, which provides an efficient means of determining a set, centered at a numerical approximation of a solution, on which a Newton-like operator is a contraction.


title: Collective strong coupling of cold atoms to an all-fiber ring cavity

abstract: We experimentally demonstrate a ring geometry all-fiber cavity system for cavity quantum electrodynamics with an ensemble of cold atoms. The fiber cavity contains a nanofiber section which mediates atom-light interactions through an evanescent field. We observe well-resolved, vacuum Rabi splitting of the cavity transmission spectrum in the weak driving limit due to a collective enhancement of the coupling rate by the ensemble of atoms within the evanescent field, and we present a simple theoretical model to describe this. In addition, we demonstrate a method to control and stabilize the resonant frequency of the cavity by utilizing the thermal properties of the nanofiber.


title: The density of mid-sized Kuiper belt objects from ALMA thermal observations

abstract: The densities of mid-sized Kuiper belt objects are a key constraint into understanding the assembly of objects in the outer solar system. These objects are critical for understanding the currently unexplained transition from the smallest Kuiper belt objects with densities lower than that of water to the largest objects with significant rock content. Mapping this transition is made difficult by the uncertainties in the diameters of these objects, which maps into an even larger uncertainty in volume and thus density. The substantial collecting area of the Atacama Large Millimeter Array allows significantly more precise measurements of thermal emission from outer solar system objects and could potentially greatly improve the density measurements. Here we use new thermal observations of four objects with satellites to explore the improvements possible with millimeter data. We find that effects due to effective emissivity at millimeter wavelengths make it difficult to use the millimeter data directly to find diameters and thus volumes for these bodies. In addition, we find that when including the effects of model uncertainty, the true uncertainties on the sizes of outer solar system objects measured with radiometry are likely larger than those previously published. Substantial improvement in object sizes will likely require precise occultation measurements.


title: On Legendrian Embbeddings into Open Book Decompositions

abstract: We study Legendrian embeddings of a compact Legendrian submanifold $L$ sitting in a closed contact manifold $(M,\xi)$ whose contact structure is supported by a (contact) open book $\mathcal{OB}$ on $M$. We prove that if $\mathcal{OB}$ has Weinstein pages, then there exist a contact structure $\xi'$ on $M$, isotopic to $\xi$ and supported by $\mathcal{OB}$, and a contactomorphism $f:(M,\xi) \to (M,\xi')$ such that the image $f(L)$ of any such submanifold can be Legendrian isotoped so that it becomes disjoint from the closure of a page of $\mathcal{OB}$.


title: Comment on "Gamma-ray spectra from low-energy positron annihilation processes in molecules"

abstract: In the article by Ma~\emph{et al.}~[Phys.~Rev.~A {\bf 94}, 052709 (2016)], $\gamma$-ray spectra for positron annihilation on molecules were calculated in the independent-particle approximation with the positron wavefunction set to unity. Based on comparisons with experimental data they concluded that inner valence electrons play a dominant role in positron annihilation. These conclusions are incorrect and resulted from fallacious analysis that ignored the known effect of the positron wavefunction on the spectra.


title: Even Galois representations and the cohomology of GL(2,Z)

abstract: Let $\rho$ be a two-dimensional even Galois representation which is induced from a character $\chi$ of odd order of the absolute Galois group of a real quadratic field. After imposing some additional conditions on $\chi$, we attach $\rho$ to a Hecke eigenclass in the cohomology of ${\rm GL}(2,\mathbb Z)$ with coefficients in a certain infinite-dimensional vector space over a field of characteristic not equal to 2.


title: Hydrodynamical and radio evolution of young supernova remnant G1.9+0.3 based on the model of diffusive shock acceleration

abstract: The radio evolution of, so far the youngest known, Galactic supernova remnant (SNR) G1.9+0.3 is investigated by using three-dimensional (3D) hydrodynamic modelling and non-linear kinetic theory of cosmic ray (CR) acceleration in SNRs. We include consistent numerical treatment of magnetic field amplification (MFA) due to resonant streaming instability. Under the assumption that SNR G1.9+0.3 is the result of a Type Ia supernova explosion located near the Galactic Centre, using widely accepted values for explosion energy 10$^{51}$ erg and ejecta mass 1.4 $M_{\odot}$, the non-thermal continuum radio emission is calculated. The main purpose of this paper is to explain radio flux brightening measured over recent decades and also predict its future temporal evolution. We estimate that the SNR is now $\sim$ 120 yr old, expanding in an ambient density of 0.02 cm$^{-3}$, and explain its steep radio spectral index only by means of efficient non-linear diffusive shock acceleration (NLDSA). We also make comparison between simulations and observations of this young SNR, in order to test the models and assumptions suggested. Our model prediction of a radio flux density increase of $\sim$ 1.8 per cent yr$^{-1}$ during the past two decades agrees well with the measured values. We synthesize the synchrotron spectrum from radio to X-ray and it fits well the VLA, MOST, Effelsberg, Chandra and NuSTAR measurements. We also propose a simplified evolutionary model of the SNR in gamma rays and suggest it may be a promising target for gamma-ray observations at TeV energies with the future generation of instruments like Cherenkov Telescope Array. SNR G1.9+0.3 is the only known Galactic SNR with the increasing flux density and we present here the prediction that the flux density will start to decrease approximately 500 yr from now. We conclude that this is a general property of SNRs in free expansion phase.


title: On Uniqueness and Blowup Properties for a Class of Second Order SDEs

abstract: As the first step for approaching the uniqueness and blowup properties of the solutions of the stochastic wave equations with multiplicative noise, we analyze the conditions for the uniqueness and blowup properties of the solution $(X_t,Y_t)$ of the equations $dX_t= Y_tdt$, $dY_t = |X_t|^\alpha dB_t$, $(X_0,Y_0)=(x_0,y_0)$. In particular, we prove that solutions are nonunique if $0<\alpha<1$ and $(x_0,y_0)=(0,0)$ and unique if $1/2<\alpha<1$ and $(x_0,y_0)\neq(0,0)$. We also show that blowup in finite time holds if $\alpha>1$ and $(x_0,y_0)\neq(0,0)$.


title: Semiclassical second microlocalization at linear coisotropic submanifolds in the torus

abstract: We develop a semiclassical second microlocal calculus of pseudodifferential operators associated to linear coisotropic submanifolds $\mathcal{C}\subset T^* \mathbb{T}^n$, where $\mathbb{T}^n = \mathbb{R}^n / \mathbb{Z}^n$. First microlocalization is localization in phase space $T^* \mathbb{T}^n$; second microlocalization is finer localization near a submanifold of $T^* \mathbb{T}^n$. Our second microlocal operators test distributions on $\mathbb{T}^n$ (e.g., Laplace eigenfunctions) for a coisotropic wavefront set, a second microlocal measure of absence of \emph{coisotropic regularity}. This wavefront set tells us where, in the coisotropic, and in what directions, approaching the coisotropic, a distribution lacks coisotropic regularity. We prove propagation theorems for coisotropic wavefront that are analogous to H\"{o}rmander's theorem for pseudodifferential operators of real principal type. Furthermore, we study the propagation of coisotropic regularity for quasimodes of semiclassical pseudodifferential operators. We Taylor expand the relevant Hamiltonian vector field, partially in the characteristic variables, at the spherical normal bundle of the coisotropic. Provided the principal symbol is real valued and depends only on the fiber variables in the cotangent bundle, and the subprincipal symbol vanishes, we show that coisotropic wavefront is invariant under the first two terms of this expansion.


title: Contraction Analysis of Nonlinear DAE Systems

abstract: This paper studies the contraction properties of nonlinear differential-algebraic equation (DAE) systems. Specifically we develop scalable techniques for constructing the attraction regions associated with a particular stable equilibrium, by establishing the relation between the contraction rates of the original systems and the corresponding virtual extended systems. We show that for a contracting DAE system, the reduced system always contracts faster than the extended ones; furthermore, there always exists an extension with contraction rate arbitrarily close to that of the original system. The proposed construction technique is illustrated with a power system example in the context of transient stability assessment.


title: sourceR: Classification and Source Attribution of Infectious Agents among Heterogeneous Populations

abstract: Zoonotic diseases are a major cause of morbidity, and productivity losses in both humans and animal populations. Identifying the source of food-borne zoonoses (e.g. an animal reservoir or food product) is crucial for the identification and prioritisation of food safety interventions. For many zoonotic diseases it is difficult to attribute human cases to sources of infection because there is little epidemiological information on the cases. However, microbial strain typing allows zoonotic pathogens to be categorised, and the relative frequencies of the strain types among the sources and in human cases allows inference on the likely source of each infection. We introduce sourceR, an R package for quantitative source attribution, aimed at food-borne diseases. It implements a fully joint Bayesian model using strain-typed surveillance data from both human cases and source samples, capable of identifying important sources of infection. The model measures the force of infection from each source, allowing for varying survivability, pathogenicity and virulence of pathogen strains, and varying abilities of the sources to act as vehicles of infection. A Bayesian non-parametric (Dirichlet process) approach is used to cluster pathogen strain types by epidemiological behaviour, avoiding model overfitting and allowing detection of strain types associated with potentially high 'virulence'. sourceR is demonstrated using Campylobacter jejuni isolate data collected in New Zealand between 2005 and 2008. It enables straightforward attribution of cases of zoonotic infection to putative sources of infection by epidemiologists and public health decision makers. As sourceR develops, we intend it to become an important and flexible resource for food-borne disease attribution studies.


title: A generalized public goods game with coupling of individual ability and project benefit

abstract: Facing a heavy task, any single person can only make a limited contribution and team cooperation is needed. As one enjoys the benefit of the public goods, the potential benefits of the project are not always maximized and may be partly wasted. By incorporating individual ability and project benefit into the original public goods game, we study the coupling effect of the four parameters, the upper limit of individual contribution, the upper limit of individual benefit, the needed project cost and the upper limit of project benefit on the evolution of cooperation. Coevolving with the individual-level group size preferences, an increase in the upper limit of individual benefit promotes cooperation while an increase in the upper limit of individual contribution inhibits cooperation. The coupling of the upper limit of individual contribution and the needed project cost determines the critical point of the upper limit of project benefit, where the equilibrium frequency of cooperators reaches its highest level. Above the critical point, an increase in the upper limit of project benefit inhibits cooperation. The evolution of cooperation is closely related to the preferred group-size distribution. A functional relation between the frequency of cooperators and the dominant group size is found.


title: Building Usage Profiles Using Deep Neural Nets

abstract: To improve software quality, one needs to build test scenarios resembling the usage of a software product in the field. This task is rendered challenging when a product's customer base is large and diverse. In this scenario, existing profiling approaches, such as operational profiling, are difficult to apply. In this work, we consider publicly available video tutorials of a product to profile usage. Our goal is to construct an automatic approach to extract information about user actions from instructional videos. To achieve this goal, we use a Deep Convolutional Neural Network (DCNN) to recognize user actions. Our pilot study shows that a DCNN trained to recognize user actions in video can classify five different actions in a collection of 236 publicly available Microsoft Word tutorial videos (published on YouTube). In our empirical evaluation we report a mean average precision of 94.42% across all actions. This study demonstrates the efficacy of DCNN-based methods for extracting software usage information from videos. Moreover, this approach may aid in other software engineering activities that require information about customer usage of a product.