title: The Old Host-Galaxy Environment of SSS17a, the First Electromagnetic Counterpart to a Gravitational Wave Source

abstract: We present an analysis of the host-galaxy environment of Swope Supernova Survey 2017a (SSS17a), the discovery of an electromagnetic counterpart to a gravitational wave source, GW170817. SSS17a occurred 1.9 kpc (in projection; 10.2") from the nucleus of NGC 4993, an S0 galaxy at a distance of 40 Mpc. We present a Hubble Space Telescope (HST) pre-trigger image of NGC 4993, Magellan optical spectroscopy of the nucleus of NGC 4993 and the location of SSS17a, and broad-band UV through IR photometry of NGC 4993. The spectrum and broad-band spectral-energy distribution indicate that NGC 4993 has a stellar mass of log (M/M_solar) = 10.49^{+0.08}_{-0.20} and star formation rate of 0.003 M_solar/yr, and the progenitor system of SSS17a likely had an age of >2.8 Gyr. There is no counterpart at the position of SSS17a in the HST pre-trigger image, indicating that the progenitor system had an absolute magnitude M_V > -5.8 mag. We detect dust lanes extending out to almost the position of SSS17a and >100 likely globular clusters associated with NGC 4993. The offset of SSS17a is similar to many short gamma-ray burst offsets, and its progenitor system was likely bound to NGC 4993. The environment of SSS17a is consistent with an old progenitor system such as a binary neutron star system.


title: The Unprecedented Properties of the First Electromagnetic Counterpart to a Gravitational Wave Source

abstract: We discovered Swope Supernova Survey 2017a (SSS17a) in the LIGO/Virgo Collaboration (LVC) localization volume of GW170817, the first detected binary neutron star (BNS) merger, only 10.9 hours after the trigger. No object was present at the location of SSS17a only a few days earlier, providing a qualitative spatial and temporal association with GW170817. Here we quantify this association, finding that SSS17a is almost certainly the counterpart of GW170817, with the chance of a coincidence being < 9 x 10^-6 (90% confidence). We arrive at this conclusion by comparing the optical properties of SSS17a to other known astrophysical transients, finding that SSS17a fades and cools faster than any other observed transient. For instance, SSS17a fades >5 mag in g within 7 days of our first data point while all other known transients of similar luminosity fade by <1 mag during the same time period. Its spectra are also unique, being mostly featureless, even as it cools. The rarity of "SSS17a-like" transients combined with the relatively small LVC localization volume and recent non-detection imply the extremely unlikely chance coincidence. We find that the volumetric rate of SSS17a-like transients is < 1.6 x 10^4 Gpc^-3 year^-1 and the Milky Way rate is <0.19 per century. A transient survey designed to discover similar events should be high cadence and observe in red filters. The LVC will likely detect substantially more BNS mergers than current optical surveys will independently discover SSS17a-like transients, however a 1-day cadence survey with LSST could discover an order of magnitude more events.


title: Maser, infrared and optical emission for late-type stars in the Galactic plane

abstract: Radio astrometric campaigns using VLBI have provided distances and proper motions for masers associated with young massive stars (BeSSeL survey). The ongoing BAaDE project plans to obtain astrometric information of SiO maser stars located in the inner Galaxy. These stars are associated with evolved, mass-losing stars. By overlapping optical (Gaia), infrared (2MASS, MSX and WISE) and radio (BAaDE) sources, we expect to obtain important clues on the intrinsic properties and population distribution of late-type stars. Moreover, a comparison of the Galactic parameters obtained with Gaia and VLBI can be done using radio observations on different targets: young massive stars (BeSSeL) and evolved stars (BAaDE).


title: The Diversity of Kilonova Emission in Short Gamma-Ray Bursts

abstract: The historic first joint detection of both gravitational wave and electromagnetic emission from a binary neutron star merger cemented the association between short gamma-ray bursts (SGRBs) and compact object mergers, as well as providing a well sampled multi-wavelength light curve of a radioactive kilonova (KN) for the first time. Here we compare the optical and near-infrared light curves of this KN, AT2017gfo, to the counterparts of a sample of nearby (z < 0.5) SGRBs to characterize their diversity in terms of their brightness distribution. Although at similar epochs AT2017gfo appears fainter than every SGRB-associated KN claimed so far, we find three bursts (GRBs 050509B, 061201 and 080905A) where, if the reported redshifts are correct, deep upper limits rule out the presence of a KN similar to AT2017gfo by several magnitudes. Combined with the properties of previously claimed KNe in SGRBs this suggests considerable diversity in the properties of KN drawn from compact object mergers, despite the similar physical conditions that are expected in many NS-NS mergers. We find that observer angle alone is not able to explain this diversity, which is likely a product of the merger type (NS-NS versus NS-BH) and the detailed properties of the binary (mass ratio, spins etc). Ultimately disentangling these properties should be possible through observations of SGRBs and gravitational wave sources, providing direct measurements of heavy element enrichment throughout the Universe.


title: Light Curves of the Neutron Star Merger GW170817/SSS17a: Implications for R-Process Nucleosynthesis

abstract: On 2017 August 17, gravitational waves were detected from a binary neutron star merger, GW170817, along with a coincident short gamma-ray burst, GRB170817A. An optical transient source, Swope Supernova Survey 17a (SSS17a), was subsequently identified as the counterpart of this event. We present ultraviolet, optical and infrared light curves of SSS17a extending from 10.9 hours to 18 days post-merger. We constrain the radioactively-powered transient resulting from the ejection of neutron-rich material. The fast rise of the light curves, subsequent decay, and rapid color evolution are consistent with multiple ejecta components of differing lanthanide abundance. The late-time light curve indicates that SSS17a produced at least ~0.05 solar masses of heavy elements, demonstrating that neutron star mergers play a role in r-process nucleosynthesis in the Universe.


title: The environment of the binary neutron star merger GW170817

abstract: We present Hubble Space Telescope and Chandra imaging, combined with Very Large Telescope MUSE integral field spectroscopy of the counterpart and host galaxy of the first binary neutron star merger detected via gravitational wave emission by LIGO & Virgo, GW170817. The host galaxy, NGC 4993, is an S0 galaxy at z=0.009783. There is evidence for large, face-on spiral shells in continuum imaging, and edge-on spiral features visible in nebular emission lines. This suggests that NGC 4993 has undergone a relatively recent (<1 Gyr) ``dry'' merger. This merger may provide the fuel for a weak active nucleus seen in Chandra imaging. At the location of the counterpart, HST imaging implies there is no globular or young stellar cluster, with a limit of a few thousand solar masses for any young system. The population in the vicinity is predominantly old with <1% of any light arising from a population with ages <500 Myr. Both the host galaxy properties and those of the transient location are consistent with the distributions seen for short-duration gamma-ray bursts, although the source position lies well within the effective radius (r_e ~ 3 kpc), providing an r_e-normalized offset that is closer than ~90% of short GRBs. For the long delay time implied by the stellar population, this suggests that the kick velocity was significantly less than the galaxy escape velocity. We do not see any narrow host galaxy interstellar medium features within the counterpart spectrum, implying low extinction, and that the binary may lie in front of the bulk of the host galaxy.


title: The first direct double neutron star merger detection: implications for cosmic nucleosynthesis

abstract: The astrophysical r-process site where about half of the elements heavier than iron are produced has been a puzzle for several decades. Here we discuss the role of neutron star mergers (NSMs) in the light of the first direct detection of such an event in both gravitational (GW) and electromagnetic (EM) waves. We analyse bolometric and NIR lightcurves of the first detected double neutron star merger and compare them to nuclear reaction network-based macronova models. The slope of the bolometric lightcurve is consistent with the radioactive decay of neutron star ejecta with $Y_e \lesssim 0.3$ (but not larger), which provides strong evidence for an r-process origin of the electromagnetic emission. This rules out in particular "nickel winds" as major source of the emission. We find that the NIR lightcurves can be well fitted either with or without lanthanide-rich ejecta. Our limits on the ejecta mass together with estimated rates directly confirm earlier purely theoretical or indirect observational conclusions that double neutron star mergers are indeed a major site of cosmic nucleosynthesis. If the ejecta mass was {\em typical}, NSMs can easily produce {\em all} of the estimated Galactic r-process matter, and --depending on the real rate-- potentially even more. This could be a hint that the event ejected a particularly large amount of mass, maybe due to a substantial difference between the component masses. This would be compatible with the mass limits obtained from the GW-observation. The recent observations suggests that NSMs are responsible for a broad range of r-process nuclei and that they are at least a major, but likely the dominant r-process site in the Universe.


title: An Ordinary Short Gamma-Ray Burst with Extraordinary Implications: Fermi-GBM Detection of GRB 170817A

abstract: On August 17, 2017 at 12:41:06 UTC the Fermi Gamma-ray Burst Monitor (GBM) detected and triggered on the short gamma-ray burst GRB 170817A. Approximately 1.7 s prior to this GRB, the Laser Interferometer Gravitational-Wave Observatory (LIGO) triggered on a binary compact merger candidate associated with the GRB. This is the first unambiguous coincident observation of gravitational waves and electromagnetic radiation from a single astrophysical source and marks the start of gravitational-wave multi-messenger astronomy. We report the GBM observations and analysis of this ordinary short GRB, which extraordinarily confirms that at least some short GRBs are produced by binary compact mergers.


title: The generating functional for the electromagnetic interaction in the strong gravitational field

abstract: This article is devoted problems of electromagnetic interaction in curved spacetime. Such problems exist, in particular, when we investigate electromagnetic quantum processes near black holes. The generalization of reduction formalism permits to find formulas for scattering matrices. For the free Dirac and electromagnetic fields corresponding generating functionals are calculated. Next we have found the generating functional for the interaction of these fields. This result holds the central position in our investigation. By means of generating functionals and reduction formulas we have obtained the scattering amplitudes for elementary electrodynamic processes: Compton scattering and the annihilation of electron-positron pair (for the tree-level approximation). On the base of these results we have formulated the generalized Feynman rules for the electromagnetic interaction in curved spacetime. Another electrodynamic processes can be studied by means of these Feynman rules and crossing-symmetry. The generating functional can be used to study problems of quantum statistics in curved spacetime.


title: GRB170817A associated with GW170817: multifrequency observations and modeling of prompt gamma-ray emission

abstract: We present our observations of electromagnetic transients associated with GW170817/GRB 170817A using optical telescopes of Chilescope observatory and Big Scanning Antenna (BSA) of Pushchino Radio Astronomy Observatory at 110MHz. The Chilescope observatory detected an optical transient of $\sim19^m$ on the third day in the outskirts of the galaxy NGC 4993; we continued observations following its rapid decrease. We put an upper limit of $1.5\times10^{4}$ Jy on any radio source with a duration of 10-60 s which may be associated with GW170817/GRB 170817A. The prompt gamma-ray emission consists of two distinctive components - a hard short pulse delayed by $\sim2$ seconds with respect to the LIGO signal and softer thermal pulse with $T\sim 10 $ keV lasting for another $\sim2$ seconds. The appearance of a thermal component at the end of the burst is unusual for short GRBs. Both the hard and the soft components do not satisfy the Amati relation, making GRB 170817A distinctively different from other short GRBs. Based on gamma-ray and optical observations, we develop a model for the prompt high-energy emission associated with GRB 170817A. The merger of two neutron stars creates an accretion torus of $\sim10^{-2} M_\odot$, which supplies the black hole with magnetic flux and confines the Blandford-Znajek-powered jet. We associate the hard prompt spike with the quasispherical breakout of the jet from the disk wind. As the jet plows through the wind with subrelativistic velocity, it creates a radiation-dominated shock that heats the wind material to tens of kiloelectron volts, producing the soft thermal component.


title: INTEGRAL Detection of the First Prompt Gamma-Ray Signal Coincident with the Gravitational Wave Event GW170817

abstract: We report the e INTernational Gamma-ray Astrophysics Laboratory (INTEGRAL) detection of the short gamma-ray burst GRB 170817A (discovered by Fermi-GBM) with a signal-to-noise ratio of 4.6, and, for the first time, its association with the gravitational waves (GWs) from binary neutron star (BNS) merging event GW170817 detected by the LIGO and Virgo observatories. The significance of association between the gamma-ray burst observed by INTEGRAL and GW170817 is 3.2 $\sigma$, while the association between the Fermi-GBM and INTEGRAL detections is 4.2 $\sigma$. GRB 170817A was detected by the SPI-ACS instrument about 2 s after the end of the gravitational wave event. We measure a fluence of $(1.4 \pm 0.4 \pm 0.6) \times$10$^{-7}$ erg cm$^{-2})$ (75--2000 keV), where, respectively, the statistical error is given at the 1 $\sigma$ confidence level, and the systematic error corresponds to the uncertainty in the spectral model and instrument response. We also report on the pointed follow-up observations carried out by INTEGRAL, starting 19.5 h after the event, and lasting for 5.4 days. We provide a stringent upper limit on any electromagnetic signal in a very broad energy range, from 3 keV to 8 MeV, constraining the soft gamma-ray afterglow flux to $<7.1\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ (80--300 keV). Exploiting the unique capabilities of INTEGRAL, we constrained the gamma-ray line emission from radioactive decays that are expected to be the principal source of the energy behind a kilonova event following a BNS coalescence. Finally, we put a stringent upper limit on any delayed bursting activity, for example from a newly formed magnetar.


title: Fermi-LAT observations of the LIGO/Virgo event GW170817

abstract: We present the Fermi Large Area Telescope (LAT) observations of the binary neutron star merger event GW170817 and the associated short gamma-ray burst (SGRB) GRB\,170817A detected by the Fermi Gamma-ray Burst Monitor. The LAT was entering the South Atlantic Anomaly at the time of the LIGO/Virgo trigger ($t_{\rm GW}$) and therefore cannot place constraints on the existence of high-energy (E $>$ 100 MeV) emission associated with the moment of binary coalescence. We focus instead on constraining high-energy emission on longer timescales. No candidate electromagnetic counterpart was detected by the LAT on timescales of minutes, hours, or days after the LIGO/Virgo detection. The resulting flux upper bound (at 95\% C.L.\/) from the LAT is $4.5\times$10$^{-10}$ erg cm$^{-2}$ s$^{-1}$ in the 0.1--1 GeV range covering a period from T0 + 1153 s to T0 + 2027 s. At the distance of GRB\,170817A, this flux upper bound corresponds to a luminosity upper bound of 9.7$\times10^{43}$ erg s$^{-1}$, which is 5 orders of magnitude less luminous than the only other LAT SGRB with known redshift, GRB\,090510. We also discuss the prospects for LAT detection of electromagnetic counterparts to future gravitational wave events from Advanced LIGO/Virgo in the context of GW170817/GRB\,170817A.


title: A generalization of moderated statistics to data adaptive semiparametric estimation in high-dimensional biology

abstract: The widespread availability of high-dimensional biological data has made the simultaneous screening of many biological characteristics a central problem in computational biology and allied sciences. While the dimensionality of such datasets continues to grow, so too does the complexity of biomarker identification from exposure patterns in health studies measuring baseline confounders; moreover, doing so while avoiding model misspecification remains an issue only partially addressed. Efficient estimators capable of incorporating flexible, data adaptive regression techniques in estimating relevant components of the data-generating distribution provide an avenue for avoiding model misspecification; however, in the context of high-dimensional problems that require the simultaneous estimation of numerous parameters, standard variance estimators have proven unstable, resulting in unreliable Type-I error control even under standard multiple testing corrections. We present a general approach for applying empirical Bayes shrinkage to variance estimators of a family of efficient, asymptotically linear estimators of population intervention causal effects. Our generalization of shrinkage-based variance estimators increases inferential stability in high-dimensional settings, facilitating the application of these estimators for deriving nonparametric variable importance measures in high-dimensional biological datasets with modest sample sizes. The result is a data adaptive approach for robustly uncovering stable causal associations in high-dimensional data in studies with limited samples. Our generalized variance estimator is evaluated against alternative variance estimators in numerical experiments. Identification of biomarkers with the proposed methodology is demonstrated in an analysis of high-dimensional DNA methylation data from an observational study on the epigenetic effects of tobacco smoking.


title: Swope Supernova Survey 2017a (SSS17a), the Optical Counterpart to a Gravitational Wave Source

abstract: On 2017 August 17, the Laser Interferometer Gravitational-wave Observatory (LIGO) and the Virgo interferometer detected gravitational waves emanating from a binary neutron star merger, GW170817. Nearly simultaneously, the Fermi and INTEGRAL telescopes detected a gamma-ray transient, GRB 170817A. 10.9 hours after the gravitational wave trigger, we discovered a transient and fading optical source, Swope Supernova Survey 2017a (SSS17a), coincident with GW170817. SSS17a is located in NGC 4993, an S0 galaxy at a distance of 40 megaparsecs. The precise location of GW170817 provides an opportunity to probe the nature of these cataclysmic events by combining electromagnetic and gravitational-wave observations.


title: A Neutron Star Binary Merger Model for GW170817/GRB170817a/SSS17a

abstract: The merging neutron star gravitational wave event GW170817 has been observed throughout the entire electromagnetic spectrum from radio waves to $\gamma$-rays. The resulting energetics, variability, and light curves are shown to be consistent with GW170817 originating from the merger of two neutron stars, in all likelihood followed by the prompt gravitational collapse of the massive remnant. The available $\gamma$-ray, X-ray and radio data provide a clear probe for the nature of the relativistic ejecta and the non-thermal processes occurring within, while the ultraviolet, optical and infrared emission are shown to probe material torn during the merger and subsequently heated by the decay of freshly synthesized $r$-process material. The simplest hypothesis that the non-thermal emission is due to a low-luminosity short $\gamma$-ray burst (sGRB) seems to agree with the present data. While low luminosity sGRBs might be common, we show here that the collective prompt and multi-wavelength observations are also consistent with a typical, powerful sGRB seen off-axis. Detailed follow-up observations are thus essential before we can place stringent constraints on the nature of the relativistic ejecta in GW170817.


title: The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. IV. Detection of Near-infrared Signatures of r-process Nucleosynthesis with Gemini-South

abstract: We present a near-infrared spectral sequence of the electromagnetic counterpart to the binary neutron star merger GW170817 detected by Advanced LIGO/Virgo. Our dataset comprises seven epochs of J+H spectra taken with FLAMINGOS-2 on Gemini-South between 1.5 and 10.5 days after the merger. In the initial epoch, the spectrum is dominated by a smooth blue continuum due to a high-velocity, lanthanide-poor blue kilonova component. Starting the following night, all subsequent spectra instead show features that are similar to those predicted in model spectra of material with a high concentration of lanthanides, including spectral peaks near 1.07 and 1.55 microns. Our fiducial model with 0.04 M_sun of ejecta, an ejection velocity of v=0.1c, and a lanthanide concentration of X_lan=1e-2 provides a good match to the spectra taken in the first five days, although it over-predicts the late-time fluxes. We also explore models with multiple fitting components, in each case finding that a significant abundance of lanthanide elements is necessary to match the broad spectral peaks that we observe starting at 2.5 d after the merger. These data provide direct evidence that binary neutron star mergers are significant production sites of even the heaviest r-process elements.


title: The Emergence of a Lanthanide-Rich Kilonova Following the Merger of Two Neutron Stars

abstract: We report the discovery and monitoring of the near-infrared counterpart (AT2017gfo) of a binary neutron-star merger event detected as a gravitational wave source by Advanced LIGO/Virgo (GW170817) and as a short gamma-ray burst by Fermi/GBM and Integral/SPI-ACS (GRB170817A). The evolution of the transient light is consistent with predictions for the behaviour of a "kilonova/macronova", powered by the radioactive decay of massive neutron-rich nuclides created via r-process nucleosynthesis in the neutron-star ejecta. In particular, evidence for this scenario is found from broad features seen in Hubble Space Telescope infrared spectroscopy, similar to those predicted for lanthanide dominated ejecta, and the much slower evolution in the near-infrared Ks-band compared to the optical. This indicates that the late-time light is dominated by high-opacity lanthanide-rich ejecta, suggesting nucleosynthesis to the 3rd r-process peak (atomic masses A~195). This discovery confirms that neutron-star mergers produce kilo-/macronovae and that they are at least a major - if not the dominant - site of rapid neutron capture nucleosynthesis in the universe.


title: The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. III. Optical and UV Spectra of a Blue Kilonova From Fast Polar Ejecta

abstract: We present optical and ultraviolet spectra of the first electromagnetic counterpart to a gravitational wave (GW) source, the binary neutron star merger GW170817. Spectra were obtained nightly between 1.5 and 9.5 days post-merger, using the SOAR and Magellan telescopes; the UV spectrum was obtained with the \textit{Hubble Space Telescope} at 5.5 days. Our data reveal a rapidly-fading blue component ($T\approx5500$ K at 1.5 days) that quickly reddens; spectra later than $\gtrsim 4.5$ days peak beyond the optical regime. The spectra are mostly featureless, although we identify a possible weak emission line at $\sim 7900$ \AA\ at $t\lesssim 4.5$ days. The colours, rapid evolution and featureless spectrum are consistent with a "blue" kilonova from polar ejecta comprised mainly of light $r$-process nuclei with atomic mass number $A\lesssim 140$. This indicates a sight-line within $\theta_{\rm obs}\lesssim 45^{\circ}$ of the orbital axis. Comparison to models suggests $\sim0.03$ M$_\odot$ of blue ejecta, with a velocity of $\sim 0.3c$. The required lanthanide fraction is $\sim 10^{-4}$, but this drops to $<10^{-5}$ in the outermost ejecta. The large velocities point to a dynamical origin, rather than a disk wind, for this blue component, suggesting that both binary constituents are neutron stars (as opposed to a binary consisting of a neutron star and a black hole). For dynamical ejecta, the high mass favors a small neutron star radius of $\lesssim 12$ km. This mass also supports the idea that neutron star mergers are a major contributor to $r$-process nucleosynthesis.


title: The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VI. Radio Constraints on a Relativistic Jet and Predictions for Late-Time Emission from the Kilonova Ejecta

abstract: We present Very Large Array (VLA) and Atacama Large Millimeter/sub-millimeter Array ALMA radio observations of GW\,170817, the first Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo gravitational wave (GW) event from a binary neutron star merger and the first GW event with an electromagnetic (EM) counterpart. Our data include the first observations following the discovery of the optical transient at both the centimeter ($13.7$ hours post merger) and millimeter ($2.41$ days post merger) bands. We detect faint emission at 6 GHz at 19.47 and 39.23 days after the merger, but not in an earlier observation at 2.46 d. We do not detect cm/mm emission at the position of the optical counterpart at frequencies of 10-97.5 GHz at times ranging from 0.6 to 30 days post merger, ruling out an on-axis short gamma-ray burst (SGRB) for energies $\gtrsim 10^{48}$ erg. For fiducial SGRB parameters, our limits require an observer viewer angle of $\gtrsim 20^{\circ}$. The radio and X-ray data can be jointly explained as the afterglow emission from an SGRB with a jet energy of $\sim 10^{49}-10^{50}$ erg that exploded in a uniform density environment with $n\sim 10^{-4}-10^{-2}$ cm$^{-3}$, viewed at an angle of $\sim 20^{\circ}-40^{\circ}$ from the jet axis. Using the results of our light curve and spectral modeling, in conjunction with the inference of the circumbinary density, we predict the emergence of late-time radio emission from the deceleration of the kilonova (KN) ejecta on a timescale of $\sim 5-10$ years that will remain detectable for decades with next-generation radio facilities, making GW\,170817 a compelling target for long-term radio monitoring.


title: The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/VIRGO GW170817. VII. Properties of the Host Galaxy and Constraints on the Merger Timescale

abstract: We present the properties of NGC 4993, the host galaxy of GW170817, the first gravitational wave (GW) event from the merger of a binary neutron star (BNS) system and the first with an electromagnetic (EM) counterpart. We use both archival photometry and new optical/near-IR imaging and spectroscopy, together with stellar population synthesis models to infer the global properties of the host galaxy. We infer a star formation history peaked at $\gtrsim 10$ Gyr ago, with subsequent exponential decline leading to a low current star formation rate of 0.01 M$_{\odot}$ yr$^{-1}$, which we convert into a binary merger timescale probability distribution. We find a median merger timescale of $11.2^{+0.7}_{-1.4}$ Gyr, with a 90% confidence range of $6.8-13.6$ Gyr. This in turn indicates an initial binary separation of $\approx 4.5$ R$_{\odot}$, comparable to the inferred values for Galactic BNS systems. We also use new and archival $Hubble$ $Space$ $Telescope$ images to measure a projected offset of the optical counterpart of $2.1$ kpc (0.64$r_{e}$) from the center of NGC 4993 and to place a limit of $M_{r} \gtrsim -7.2$ mag on any pre-existing emission, which rules out the brighter half of the globular cluster luminosity function. Finally, the age and offset of the system indicates it experienced a modest natal kick with an upper limit of $\sim 200$ km s$^{-1}$. Future GW$-$EM observations of BNS mergers will enable measurement of their population delay time distribution, which will directly inform their viability as the dominant source of $r$-process enrichment in the Universe.


title: The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. I. Dark Energy Camera Discovery of the Optical Counterpart

abstract: We present the Dark Energy Camera (DECam) discovery of the optical counterpart of the first binary neutron star merger detected through gravitational wave emission, GW170817. Our observations commenced 10.5 hours post-merger, as soon as the localization region became accessible from Chile. We imaged 70 deg$^2$ in the $i$ and $z$ bands, covering 93\% of the initial integrated localization probability, to a depth necessary to identify likely optical counterparts (e.g., a kilonova). At 11.4 hours post-merger we detected a bright optical transient located $10.6''$ from the nucleus of NGC\,4993 at redshift $z=0.0098$, consistent (for $H_0 = 70$\, km s$^{-1}$ Mpc$^{-1}$) with the distance of $40 \pm 8$\, Mpc reported by the LIGO Scientific Collaboration and the Virgo Collaboration (LVC). At detection the transient had magnitudes $i\approx 17.30$ and $z\approx 17.45$, and thus an absolute magnitude of $M_i = -15.7$, in the luminosity range expected for a kilonova. We identified 1,500 potential transient candidates. Applying simple selection criteria aimed at rejecting background events such as supernovae, we find the transient associated with NGC\,4993 as the only remaining plausible counterpart, and reject chance coincidence at the 99.5\% confidence level. We therefore conclude that the optical counterpart we have identified near NGC\,4993 is associated with GW170817. This discovery ushers in the era of multi-messenger astronomy with gravitational waves, and demonstrates the power of DECam to identify the optical counterparts of gravitational-wave sources.


title: AGILE Observations of the Gravitational-wave Source GW170817: Constraining Gamma-Ray Emission from a NS-NS Coalescence

abstract: The LIGO-Virgo Collaboration (LVC) detected, on 2017 August 17, an exceptional gravitational-wave (GW) event temporally consistent within $\sim\,1.7 \, \rm s$ with the GRB 1708117A observed by Fermi-GBM and INTEGRAL. The event turns out to be compatible with a neutron star-neutron star (NS-NS) coalescence that subsequently produced a radio/optical/X-ray transient detected at later times. We report the main results of the observations by the AGILE satellite of the GW170817 localization region (LR) and its electromagnetic (e.m.) counterpart. At the LVC detection time $T_0$, the GW170817 LR was occulted by the Earth. The AGILE instrument collected useful data before and after the GW-GRB event because in its spinning observation mode it can scan a given source many times per hour. The earliest exposure of the GW170817 LR by the gamma-ray imaging detector (GRID) started about 935 s after $T_0$. No significant X-ray or gamma-ray emission was detected from the LR that was repeatedly exposed over timescales of minutes, hours, and days before and after GW170817, also considering Mini-calorimeter and Super-AGILE data. Our measurements are among the earliest ones obtained by space satellites on GW170817 and provide useful constraints on the precursor and delayed emission properties of the NS-NS coalescence event. We can exclude with high confidence the existence of an X-ray/gamma-ray emitting magnetar-like object with a large magnetic field of $10^{15} \, \rm G$. Our data are particularly significant during the early stage of evolution of the e.m. remnant.


title: MASTER optical detection of the first LIGO/Virgo neutron stars merging GW170817

abstract: Following the reported discovery of the gravitational-wave pulse GW170817/ G298048 by three LIGO/Virgo antennae (Abbott et al., 2017a), the MASTER Global Robotic Net telescopes obtained the first image of the NGC 4993 galaxy after the NS+NS merging. The optical transient MASTER OTJ130948.10-232253.3/SSS17a was later found, which appears to be a kilonova resulting from a merger of two neutron stars. In this paper we report the independent detection and photometry of the kilonova made in white light and in B, V, and R filters. We note that luminosity of the discovered kilonova NGC 4993 is very close to another possible kilonova proposed early GRB 130603 and GRB 080503.


title: Optical Observations of LIGO Source GW 170817 by the Antarctic Survey Telescopes at Dome A, Antarctica

abstract: The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW~170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the $i$-band magnitude declining from $17.23\pm0.13$ magnitude to $17.72\pm0.09$ magnitude in $\sim 1.8$ hour. The brightness and time evolution of the optical transient associated with GW~170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about $\sim 10^{-2}$ solar mass of radioactive material at a speed of up to $30\%$ the speed of light.


title: Origin of the heavy elements in binary neutron-star mergers from a gravitational wave event

abstract: The cosmic origin of the elements heavier than iron has long been uncertain. Theoretical modelling shows that the matter that is expelled in the violent merger of two neutron stars can assemble into heavy elements such as gold and platinum in a process known as rapid neutron capture (r-process) nucleosynthesis. The radioactive decay of isotopes of the heavy elements is predicted to power a distinctive thermal glow (a 'kilonova'). The discovery of an electromagnetic counterpart to the gravitational-wave source GW170817 represents the first opportunity to detect and scrutinize a sample of freshly synthesized r-process elements. Here we report models that predict the detailed electromagnetic emission of kilonovae and enable the mass, velocity and composition of ejecta to be derived from the observations. We compare the models to the optical and infrared radiation associated with GW170817 event to argue that the observed source is a kilonova. We infer the presence of two distinct components of ejecta, one composed primarily of light (atomic mass number less than 140) and one of heavy (atomic mass number greater than 140) r-process elements. Inferring the ejected mass and a merger rate from GW170817 implies that such mergers are a dominant mode of r-process production in the Universe.


title: Learning the seasonality of disease incidences from empirical data

abstract: Investigating the seasonality of disease incidences is very important in disease surveillance in regions with periodical climatic patterns. In lieu of the paradigm about disease incidences varying seasonally in line with meteorology, this work seeks to determine how well standard epidemic models can capture such seasonality for better forecasts and optimal futuristic interventions. Once incidence data are assimilated by a periodic model, asymptotic analysis in relation to the long-term behavior of the disease occurrences can be performed using the classical Floquet theory, which explains the stability of the existing periodic solutions. For a test case, we employed an IR model to assimilate weekly dengue incidence data from the city of Jakarta, Indonesia, which we present in their raw and moving-average-filtered versions. To estimate a periodic parameter toward performing the asymptotic analysis, eight optimization schemes were assigned returning magnitudes of the parameter that vary insignificantly across schemes. Furthermore, the computation results combined with the analytical results indicate that if the disease surveillance in the city does not improve, then the incidence will raise to a certain positive orbit and remain cyclical.


title: Flow: A Modular Learning Framework for Mixed Autonomy Traffic

abstract: The rapid development of autonomous vehicles (AVs) holds vast potential for transportation systems through improved safety, efficiency, and access to mobility. However, the progression of these impacts, as AVs are adopted, is not well understood. Numerous technical challenges arise from the goal of analyzing the partial adoption of autonomy: partial control and observation, multi-vehicle interactions, and the sheer variety of scenarios represented by real-world networks. To shed light into near-term AV impacts, this article studies the suitability of deep reinforcement learning (RL) for overcoming these challenges in a low AV-adoption regime. A modular learning framework is presented, which leverages deep RL to address complex traffic dynamics. Modules are composed to capture common traffic phenomena (stop-and-go traffic jams, lane changing, intersections). Learned control laws are found to improve upon human driving performance, in terms of system-level velocity, by up to 57% with only 4-7% adoption of AVs. Furthermore, in single-lane traffic, a small neural network control law with only local observation is found to eliminate stop-and-go traffic - surpassing all known model-based controllers to achieve near-optimal performance - and generalize to out-of-distribution traffic densities.


title: Large phonon time-of-flight fluctuations in expanding flat condensates of cold fermi gases

abstract: We reexamine how quantum density fluctuations in condensates of ultra-cold fermi gases lead to fluctuations in phonon times-of-flight, an effect that increases as density is reduced. We suggest that these effects should be measurable in pancake-like (two-dimensional) condensates on their release from their confining optical traps, providing their initial (width/thickness) aspect ratio is suitably large.


title: A novel study of temperature effects on the viscoelastic behavior of articular cartilage

abstract: This paper presents a new approach to study the effects of temperature on the poro- elastic and viscoelastic behavior of articular cartilage. Biphasic solid-fluid mixture theory is applied to study the poro-mechancial behavior of articular cartilage in a fully saturated state. The balance of linear momentum, mass, and energy are considered to describe deformation of the solid skeleton, pore fluid pressure, and temperature distribution in the mixture. The mechanical model assumes both linear elastic and viscoelastic isotropic materials, infinitesimal strain theory, and a time-dependent response. The influence of temperature on the mixture behavior is modeled through temperature dependent mass density and volumetric thermal strain. The fluid flow through the porous medium is described by the Darcy's law. The stress-strain relation for time-dependent viscoelastic deformation in the solid skeleton is described using the generalized Maxwell model. A verification example is presented to illustrate accuracy and efficiency of the developed finite element model. The influence of temperature is studied through examining the behavior of articular cartilage for confined and unconfined boundary conditions. Furthermore, articular cartilage under partial loading condition is modeled to investigate the deformation, pore fluid pressure, and temperature dissipation processes. The results suggest significant impacts of temperature on both poro- elastic and viscoelastic behavior of articular cartilage.


title: Generalization in Deep Learning

abstract: This paper provides theoretical insights into why and how deep learning can generalize well, despite its large capacity, complexity, possible algorithmic instability, nonrobustness, and sharp minima, responding to an open question in the literature. We also discuss approaches to provide non-vacuous generalization guarantees for deep learning. Based on theoretical observations, we propose new open problems and discuss the limitations of our results.


title: 800 nm pumped SPDC from chirped crystals

abstract: We consider crystal chirp effects on SPDC when pumping at 800 nm. The typical distribution produced in frequency-momentum space is a pop tab-like structure which turns out to be suitable for the implementation of versatile light sources. Our analyzes consider the effect of internal and external parameters in the process; in the former we include the crystal chirp and length, while in the latter temperature, as well as pump chirp and beam properties. We report evidence of the appropriateness of SPDC from chirped crystals to manipulate the frequency and transverse momentum properties of the light produced. We briefly comment on potential usefulness of the types of light produced, in particular for quantum information applications.


title: Asynchronous Early Output Section-Carry Based Carry Lookahead Adder with Alias Carry Logic

abstract: A new asynchronous early output section-carry based carry lookahead adder (SCBCLA) with alias carry output logic is presented in this paper. To evaluate the proposed SCBCLA with alias carry logic and to make a comparison with other CLAs, a 32-bit addition operation is considered. Compared to the weak-indication SCBCLA with alias logic, the proposed early output SCBCLA with alias logic reports a 13% reduction in area without any increases in latency and power dissipation. On the other hand, in comparison with the early output recursive CLA (RCLA), the proposed early output SCBCLA with alias logic reports a 16% reduction in latency while occupying almost the same area and dissipating almost the same average power. All the asynchronous CLAs are quasi-delay-insensitive designs which incorporate the delay-insensitive dual-rail data encoding and adhere to the 4-phase return-to-zero handshaking. The adders were realized and the simulations were performed based on a 32/28nm CMOS process.


title: Topologically protected elastic waves in one-dimensional phononic crystals of continuous media

abstract: We report the design of silica-based 1D phononic crystals (PnCs) with topologically distinct complete phononic bandgaps (PnBGs) and the observation of a topologically protected state of elastic waves at their interface. By choosing different structural parameters of unit cells, two PnCs can possess a common PnBG with different topological nature. At the interface between the two PnCs, a topological interface mode with a quality factor of ~5,650 is observed in the PnBG. Spatial confinement of the interface mode is also confirmed by using photoelastic imaging technique. Such topologically protected elastic states are potentially applicable for constructing novel phononic devices.


title: Safe Learning of Quadrotor Dynamics Using Barrier Certificates

abstract: To effectively control complex dynamical systems, accurate nonlinear models are typically needed. However, these models are not always known. In this paper, we present a data-driven approach based on Gaussian processes that learns models of quadrotors operating in partially unknown environments. What makes this challenging is that if the learning process is not carefully controlled, the system will go unstable, i.e., the quadcopter will crash. To this end, barrier certificates are employed for safe learning. The barrier certificates establish a non-conservative forward invariant safe region, in which high probability safety guarantees are provided based on the statistics of the Gaussian Process. A learning controller is designed to efficiently explore those uncertain states and expand the barrier certified safe region based on an adaptive sampling scheme. In addition, a recursive Gaussian Process prediction method is developed to learn the complex quadrotor dynamics in real-time. Simulation results are provided to demonstrate the effectiveness of the proposed approach.


title: On the decomposition numbers of $\mathrm{SO}_8^+(2^f)$

abstract: Let $q=2^f$, and let $G=\mathrm{SO}_8^+(q)$ and $U$ be a Sylow $2$-subgroup of $G$. We first describe the fusion of the conjugacy classes of $U$ in $G$. We then use this information to prove the unitriangularity of the $\ell$-decomposition matrices of $G$ for all $\ell \ne 2$ by inducing certain irreducible characters of $U$ to $G$; the characters of $U$ of degree $q^3/2$ play here a major role. We then determine the $\ell$-decomposition matrix of $G$ in the case $\ell \mid q+1$, when $\ell \ge 5$ and $(q+1)_\ell>5$, up to two non-negative indeterminates in one column.


title: Approximate Ripple Carry and Carry Lookahead Adders - A Comparative Analysis

abstract: Approximate ripple carry adders (RCAs) and carry lookahead adders (CLAs) are presented which are compared with accurate RCAs and CLAs for performing a 32-bit addition. The accurate and approximate RCAs and CLAs are implemented using a 32/28nm CMOS process. Approximations ranging from 4- to 20-bits are considered for the less significant adder bit positions. The simulation results show that approximate RCAs report reductions in the power-delay product (PDP) ranging from 19.5% to 82% than the accurate RCA for approximation sizes varying from 4- to 20-bits. Also, approximate CLAs report reductions in PDP ranging from 16.7% to 74.2% than the accurate CLA for approximation sizes varying from 4- to 20-bits. On average, for the approximation sizes considered, it is observed that approximate CLAs achieve a 46.5% reduction in PDP compared to the approximate RCAs. Hence, approximate CLAs are preferable over approximate RCAs for the low power implementation of approximate computer arithmetic.


title: Note on the gravitational and electromagnetic radiation for the Einstein-Maxwell equations with cosmological constant

abstract: In the middle of last century, Bondi and his coworkers proposed an out going boundary condition for the Einstein equations. Based on such boundary condition the authors theoretically solved the puzzle of the existence problem of gravitational wave. Currently many works on gravitational wave source modeling use Bondi's result to do the analysis including the gravitational wave calculation in numerical relativity. Recently more and more observations imply that the Einstein equations should be modified with an nonzero cosmological constant. In this work we use Bondi's original method to treat the Einstein equations with cosmological constant for theoretical curiosity. We find that the gravitational wave does not essentially exist if Bondi's original out going boundary condition is imposed. We find this unphysical result is due to the non-consistency between Bondi's original out going boundary condition and the Einstein equations with cosmological constant. Based on theoretical analysis, we propose an new Bondi-type out going boundary condition for the Einstein equations with cosmological constant. With this new boundary condition, the gravitational wave behavior for the Einstein equations with cosmological constant is similar to the Einstein equations without cosmological constant. In particular, when the cosmological constant goes to zero, the behavior of the Einstein equations without cosmological constant is recovered. Together with gravitational wave, electromagnetic wave is also analyzed. The boundary conditions for electromagnetic wave does not depend on the cosmological constant.


title: Calibrated Boosting-Forest

abstract: Excellent ranking power along with well calibrated probability estimates are needed in many classification tasks. In this paper, we introduce a technique, Calibrated Boosting-Forest that captures both. This novel technique is an ensemble of gradient boosting machines that can support both continuous and binary labels. While offering superior ranking power over any individual regression or classification model, Calibrated Boosting-Forest is able to preserve well calibrated posterior probabilities. Along with these benefits, we provide an alternative to the tedious step of tuning gradient boosting machines. We demonstrate that tuning Calibrated Boosting-Forest can be reduced to a simple hyper-parameter selection. We further establish that increasing this hyper-parameter improves the ranking performance under a diminishing return. We examine the effectiveness of Calibrated Boosting-Forest on ligand-based virtual screening where both continuous and binary labels are available and compare the performance of Calibrated Boosting-Forest with logistic regression, gradient boosting machine and deep learning. Calibrated Boosting-Forest achieved an approximately 48% improvement compared to a state-of-art deep learning model. Moreover, it achieved around 95% improvement on probability quality measurement compared to the best individual gradient boosting machine. Calibrated Boosting-Forest offers a benchmark demonstration that in the field of ligand-based virtual screening, deep learning is not the universally dominant machine learning model and good calibrated probabilities can better facilitate virtual screening process.


title: A multi-branch convolutional neural network for detecting double JPEG compression

abstract: Detection of double JPEG compression is important to forensics analysis. A few methods were proposed based on convolutional neural networks (CNNs). These methods only accept inputs from pre-processed data, such as histogram features and/or decompressed images. In this paper, we present a CNN solution by using raw DCT (discrete cosine transformation) coefficients from JPEG images as input. Considering the DCT sub-band nature in JPEG, a multiple-branch CNN structure has been designed to reveal whether a JPEG format image has been doubly compressed. Comparing with previous methods, the proposed method provides end-to-end detection capability. Extensive experiments have been carried out to demonstrate the effectiveness of the proposed network.


title: On estimates of transition density for subordinate Brownian motions with Gaussian components in $C^{1,1}$-open sets

abstract: We consider a subordinate Brownian motion $X$ with Gaussian components when the scaling order of purely discontinuous part is between $0$ and $2$ including $2$. In this paper we establish sharp two-sided bounds for transition density of $X$ in ${\mathbb R}^d$ and $C^{1,1}$-open sets. As a corollary, we obtain a sharp Green function estimates.


title: Novel Fiber Design for Wideband Conversion and Amplification in Multimode Fibers

abstract: We propose an operating principle to achieve broadband and highly tunable mode conversion and amplification exploiting inter-modal four wave mixing in a multimode fiber. A bandwidth of 30 nanometers is demonstrated by properly designing a simple step-index silica fiber.


title: Fraction-variant beam orientation optimization for non-coplanar IMRT

abstract: Conventional beam orientation optimization (BOO) algorithms for IMRT assume that the same set of beam angles is used for all treatment fractions. In this paper we present a BOO formulation based on group sparsity that simultaneously optimizes non-coplanar beam angles for all fractions, yielding a fraction-variant (FV) treatment plan. Beam angles are selected by solving a multi-fraction FMO problem involving 500-700 candidate beams per fraction, with an additional group sparsity term that encourages most candidate beams to be inactive. The optimization problem is solved using the Fast Iterative Shrinkage-Thresholding Algorithm. Our FV BOO algorithm is used to create non-coplanar, five-fraction treatment plans for prostate and lung cases, as well as a non-coplanar 30-fraction plan for a head and neck case. A homogeneous PTV dose coverage is maintained in all fractions. The treatment plans are compared with fraction-invariant plans that use a fixed set of beam angles for all fractions. The FV plans reduced mean and max OAR dose on average by 3.3% and 3.7% of the prescription dose, respectively. Notably, mean OAR dose was reduced by 14.3% of prescription dose (rectum), 11.6% (penile bulb), 10.7% (seminal vesicle), 5.5% (right femur), 3.5% (bladder), 4.0% (normal left lung), 15.5% (cochleas), and 5.2% (chiasm). Max OAR dose was reduced by 14.9% of prescription dose (right femur), 8.2% (penile bulb), 12.7% (prox. bronchus), 4.1% (normal left lung), 15.2% (cochleas), 10.1% (orbits), 9.1% (chiasm), 8.7% (brainstem), and 7.1% (parotids). Meanwhile, PTV homogeneity defined as D95/D5 improved from .95 to .98 (prostate case) and from .94 to .97 (lung case), and remained constant for the head and neck case. Moreover, the FV plans are dosimetrically similar to conventional plans that use twice as many beams per fraction. Thus, FV BOO offers the potential to reduce delivery time for non-coplanar IMRT.


title: Small-depth Multilinear Formula Lower Bounds for Iterated Matrix Multiplication, with Applications

abstract: In this paper, we study the algebraic formula complexity of multiplying $d$ many $2\times 2$ matrices, denoted $\mathrm{IMM}_{d}$, and show that the well-known divide-and-conquer algorithm cannot be significantly improved at any depth, as long as the formulas are multilinear. Formally, for each depth $\Delta \leq \log d$, we show that any product-depth $\Delta$ multilinear formula for $\mathrm{IMM}_d$ must have size $\exp(\Omega(\Delta d^{1/\Delta})).$ It also follows from this that any multilinear circuit of product-depth $\Delta$ for the same polynomial of the above form must have a size of $\exp(\Omega(d^{1/\Delta})).$ In particular, any polynomial-sized multilinear formula for $\mathrm{IMM}_d$ must have depth $\Omega(\log d)$, and any polynomial-sized multilinear circuit for $\mathrm{IMM}_d$ must have depth $\Omega(\log d/\log \log d).$ Both these bounds are tight up to constant factors. 1. Depth-reduction: A well-known result of Brent (JACM 1974) implies that any formula of size $s$ can be converted to one of size $s^{O(1)}$ and depth $O(\log s)$; further, this reduction continues to hold for multilinear formulas. Our lower bound implies that any depth-reduction in the multilinear setting cannot reduce the depth to $o(\log s)$ without a superpolynomial blow-up in size. 2. Separations from general formulas: Our result, along with a non-trivial upper bound for $\mathrm{IMM}_{d}$ implied by a result of Gupta, Kamath, Kayal and Saptharishi (SICOMP 2016), shows that for any size $s$ and product-depth $\Delta = o(\log s),$ general formulas of size $s$ and product-depth $\Delta$ cannot be converted to multilinear formulas of size $s^{\omega(1)}$ and product-depth $\Delta,$ when the underlying field has characteristic zero.


title: Spreading with two speeds and mass segregation in a diffusive competition system with free boundaries

abstract: We investigate the spreading behavior of two invasive species modeled by a Lotka-Volterra diffusive competition system with two free boundaries in a spherically symmetric setting. We show that, for the weak-strong competition case, under suitable assumptions, both species in the system can successfully spread into the available environment, but their spreading speeds are different, and their population masses tend to segregate, with the slower spreading competitor having its population concentrating on an expanding ball, say Bt, and the faster spreading competitor concentrating on a spherical shell outside Bt that disappears to infinity as time goes to infinity.


title: Using Deep Learning and Satellite Imagery to Quantify the Impact of the Built Environment on Neighborhood Crime Rates

abstract: The built environment has been postulated to have an impact on neighborhood crime rates, however, measures of the built environment can be subjective and differ across studies leading to varying observations on its association with crime rates. Here, we illustrate an accurate and straightforward approach to quantify the impact of the built environment on neighborhood crime rates from high-resolution satellite imagery. Using geo-referenced crime reports and satellite images for three United States cities, we demonstrate how image features consistently identified using a convolutional neural network can explain up to 82% of the variation in neighborhood crime rates. Our results suggest the built environment is a strong predictor of crime rates, and this can lead to structural interventions shown to reduce crime incidence in urban settings.


title: Global Lagrangian solutions of the Camassa-Holm equation

abstract: In this paper we construct global weak conservative solutions of the Camassa-Holm equation on the periodic domain. We first express the equation in Lagrangian flow variable $\eta$ and then transform it using a change of variable $\rho=\sqrt{\eta_x}$. The new variable removes the singularity of the CH equation, and we obtain both the global weak conservative solution and global spatial smoothness of the Lagrangian trajectories of the CH equation. This work is motivated by J. Lenells who proved similar results for the Hunter-Saxton equation using the geometric interpretation.


title: Asymptotic spreading speed for the weak competition system with a free boundary

abstract: This paper is concerned with a diffusive Lotka-Volterra type competition system with a free boundary in one space dimension. Such a system may be used to describe the invasion of a new species into the habitat of a native competitor. We show that the longtime dynamical behavior of the system is determined by a spreading-vanishing dichotomy, and provide sharp criteria for spreading and vanishing of the invasive species. Moreover, we determine the asymptotic spreading speed of the invasive species when its spreading is successful, which involves two systems of traveling wave type equations, and is highly nontrivial to establish.


title: Two-component instantaneous wall shear rate measurements in reversing flows

abstract: Validation experiments of the two-dimensional inverse algorithm are performed in a pulsed Poiseuille flow exposing shear reversal phases. The method is applied to the three-segment electrodiffusion (ED) probe for which a specific nondimensionalization process is suggested, allowing to better link measurements from a real ED probe to the modeled one in the inverse problem. This approach provided a two-component wall shear rate in good agreement with the one obtained from laser Doppler anemometry (LDA) measurements, thus validating the ability of ED probes to deal with high-amplitude unsteady flows. The classic linear velocity approximation ($u=sy$) in the probe vicinity is also investigated in such a flow.


title: SpooQySats: cubesats to demonstrate quantum key distribution technologies

abstract: Satellite-based QKD offers the potential to share highly secure encryption keys between optical ground stations all over the planet. SpooQySats is a programme for establishing the space worthiness of highly-miniaturised, polarization entangled, photon pair sources using CubeSat nanosatellites. The sources are being developed iteratively with an early version in orbit already and improved versions soon to be launched. Once fully developed, the photon pair sources can be deployed on more advanced satellites that are equipped with optical links. These can allow for very secure uplinks and downlinks and can be used to establish a global space-based quantum key distribution network. This would enable highly secure symmetric encryption keys to be shared between optical ground stations all over the planet.


title: A Geometric View of Optimal Transportation and Generative Model

abstract: In this work, we show the intrinsic relations between optimal transportation and convex geometry, especially the variational approach to solve Alexandrov problem: constructing a convex polytope with prescribed face normals and volumes. This leads to a geometric interpretation to generative models, and leads to a novel framework for generative models. By using the optimal transportation view of GAN model, we show that the discriminator computes the Kantorovich potential, the generator calculates the transportation map. For a large class of transportation costs, the Kantorovich potential can give the optimal transportation map by a close-form formula. Therefore, it is sufficient to solely optimize the discriminator. This shows the adversarial competition can be avoided, and the computational architecture can be simplified. Preliminary experimental results show the geometric method outperforms WGAN for approximating probability measures with multiple clusters in low dimensional space.


title: Demographics of Star-forming Galaxies since $z\sim2.5$. I. The $UVJ$ Diagram in CANDELS

abstract: This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.


title: Probabilistic cellular automata with memory two: invariant laws and multidirectional reversibility

abstract: We focus on a family of one-dimensional probabilistic cellular automata with memory two: the dynamics is such that the value of a given cell at time $t+1$ is drawn according to a distribution which is a function of the states of its two nearest neighbours at time $t$, and of its own state at time $t-1$. Such PCA naturally arise in the study of some models coming from statistical physics ($8$-vertex model, directed animals and gaz models, TASEP, etc.). We give conditions for which the invariant measure has a product form or a Markovian form, and we prove an ergodicity result holding in that context. The stationary space-time diagrams of these PCA present different forms of reversibility. We describe and study extensively this phenomenon, which provides families of Gibbs random fields on the square lattice having nice geometric and combinatorial properties.


title: Balanced Judicious Partition is Fixed-Parameter Tractable

abstract: The family of judicious partitioning problems, introduced by Bollob\'as and Scott to the field of extremal combinatorics, has been extensively studied from a structural point of view for over two decades. This rich realm of problems aims to counterbalance the objectives of classical partitioning problems such as Min Cut, Min Bisection and Max Cut. While these classical problems focus solely on the minimization/maximization of the number of edges crossing the cut, judicious (bi)partitioning problems ask the natural question of the minimization/maximization of the number of edges lying in the (two) sides of the cut. In particular, Judicious Bipartition (JB) seeks a bipartition that is "judicious" in the sense that neither side is burdened by too many edges, and Balanced JB also requires that the sizes of the sides themselves are "balanced" in the sense that neither of them is too large. Both of these problems were defined in the work by Bollob\'as and Scott, and have received notable scientific attention since then. In this paper, we shed light on the study of judicious partitioning problems from the viewpoint of algorithm design. Specifically, we prove that BJB is FPT (which also proves that JB is FPT).


title: Surjections of unit groups and semi-inverses

abstract: Given a surjective ring homomorphism, we study when the induced group homomorphism on unit groups is surjective. To this end, we introduce notions of generalized inverses and units, as well as a class of rings such that the set of closed points in the spectrum is a closed set. It is shown that any surjection out of such a ring induces a surjection on unit groups.


title: An application of Nakayama functor in representation stability theory

abstract: Using the Nakayama functor, we construct an equivalence from a Serre quotient category of a category of finitely generated modules to a category of finite-dimensional modules. We then apply this result to the categories FI$_G$ and VI$_q$, and answer positively an open question of Nagpal on representation stability theory.


title: REM sleep complicates period adding bifurcations from monophasic to polyphasic sleep behavior in a sleep-wake regulatory network model for human sleep

abstract: The structure of human sleep changes across development as it consolidates from the polyphasic sleep of infants to the single nighttime sleep period typical in adults. Across this same developmental period, time scales of the homeostatic sleep drive, the physiological drive to sleep that increases with time spent awake, also change and presumably govern the transition from polyphasic to monophasic sleep behavior. Using a physiologically-based, sleep-wake regulatory network model for human sleep, we investigated the dynamics of wake, rapid eye movement (REM) sleep, and non-REM (NREM) sleep during this transition by varying the homeostatic sleep drive time constants. Previously, we introduced an algorithm for constructing a one-dimensional circle map that represents the dynamics of the full sleep-wake network model. By tracking bifurcations in the piecewise continuous circle map as the homeostatic sleep drive time constants are varied, we establish evidence for a border collision bifurcation that results in period-adding-like behavior in the number of sleep cycles per day. Interestingly, this bifurcation is preceded by bifurcations in the number of REM bouts per sleep cycle that exhibit truncated period-adding-like behavior. The interaction of these bifurcations in numbers of sleep episodes and numbers of REM bouts per sleep episode generate non-monotonic variation in sleep cycle patterns as well as quasi-periodic patterns during the transition from polyphasic to monophasic sleep behavior. This analysis may have implications for understanding changes in sleep in early childhood when preschoolers transition from napping to non-napping behavior, and the wide interindividual variation observed during this transition.


title: Fast two-snapshot structured illumination for temporal focusing microscopy with enhanced axial resolution

abstract: We present a new two-snapshot structured light illumination (SLI) reconstruction algorithm for fast image acquisition. The new algorithm, which only requires two mutually {\pi} phase-shifted raw structured images, is implemented on a custom-built temporal focusing fluorescence microscope (TFFM) to enhance its axial resolution via a digital micromirror device (DMD). First, the orientation of the modulated sinusoidal fringe patterns is automatically identified via spatial frequency vector detection. Subsequently, the modulated in-focal-plane images are obtained via rotation and subtraction. Lastly, a parallel amplitude demodulation method, derived based on Hilbert transform, is applied to complete the decoding processes. To demonstrate the new SLI algorithm, a TFFM is custom-constructed, where a DMD replaces the generic blazed grating in the system and simultaneously functions as a diffraction grating and a programmable binary mask, generating arbitrary fringe patterns. The experimental results show promising depth-discrimination capability with an axial resolution enhancement factor of 1.25, which matches well with the theoretical estimation, i.e, 1.27. Imaging experiments on pollen grain samples have been performed. The results indicate that the two-snapshot algorithm presents comparable contrast reconstruction and optical cross-sectioning capability with those adopting the conventional root-mean-square (RMS) reconstruction method. The two-snapshot method can be readily applied to any sinusoidally modulated illumination systems to realize high-speed 3D imaging as less frames are required for each in-focal-plane image restoration, i.e., the image acquisition speed is improved by 2.5 times for any two-photon systems.


title: Infinite prime avoidance

abstract: We investigate prime avoidance for an arbitrary set of prime ideals in a commutative ring. Various necessary and/or sufficient conditions for prime avoidance are given, which yield natural classes of infinite sets of primes that satisfy prime avoidance. Examples and counterexamples are given throughout to illustrate the phenomena that can occur. As an application, we show how to use prime avoidance to construct counterexamples among rings essentially of finite type.


title: Quantum oscillations and Dirac-Landau levels in Weyl superconductors

abstract: When magnetic field is applied to metals and semimetals quantum oscillations appear as individual Landau levels cross the Fermi level. Quantum oscillations generally do not occur in superconductors (SC) because magnetic field is either expelled from the sample interior or, if strong enough, drives the material into the normal state. In addition, elementary excitations of a superconductor -- Bogoliubov quasiparticles -- do not carry a well defined electric charge and therefore do not couple in a simple way to the applied magnetic field. We predict here that in Weyl superconductors certain types of elastic strain have the ability to induce chiral pseudo-magnetic field which can reorganize the electronic states into Dirac-Landau levels with linear band crossings at low energy. The resulting quantum oscillations in the quasiparticle density of states and thermal conductivity can be experimentally observed under a bending deformation of a thin film Weyl SC and provide new insights into this fascinating family of materials.


title: Nonnegative Ricci curvature, stability at infinity, and finite generation of fundamental groups

abstract: We study the fundamental group of an open $n$-manifold $M$ of nonnegative Ricci curvature. We show that if there is an integer $k$ such that any tangent cone at infinity of the Riemannian universal cover of $M$ is a metric cone, whose maximal Euclidean factor has dimension $k$, then $\pi_1(M)$ is finitely generated. In particular, this confirms the Milnor conjecture for a manifold whose universal cover has Euclidean volume growth and the unique tangent cone at infinity.


title: Computation Offloading and Activation of Mobile Edge Computing Servers: A Minority Game

abstract: With the ever-increasing popularity of resource-intensive mobile applications, Mobile Edge Computing (MEC), e.g., offloading computationally expensive tasks to the cellular edge, has become a prominent technology for the next generation wireless networks. Despite its great performance in terms of delay and energy, MEC suffers from restricted power allowance and computational capability of the edge nodes. Therefore, it is imperative to develop distributed mechanisms for computation offloading, so that not only the computational servers are utilized at their best capacity, but also the users' latency constraints are fulfilled. In this paper, by using the theory of Minority Games, we develop a novel distributed server activation mechanism for computational offloading. Our scheme guarantees energy-efficient activation of servers as well as satisfaction of users' quality-of-experience (QoE) requirements in terms of latency.


title: Multiscale convergence properties for spectral approximations of a model kinetic equation

abstract: In this work, we prove rigorous convergence properties for a semi-discrete, moment-based approximation of a model kinetic equation in one dimension. This approximation is equivalent to a standard spectral method in the velocity variable of the kinetic distribution and, as such, is accompanied by standard algebraic estimates of the form $N^{-q}$, where $N$ is the number of modes and $q>0$ depends on the regularity of the solution. However, in the multiscale setting, the error estimate can be expressed in terms of the scaling parameter $\epsilon$, which measures the ratio of the mean-free-path to the characteristic domain length. We show that, for isotropic initial conditions, the error in the spectral approximation is $\mathcal{O}(\epsilon^{N+1})$. More surprisingly, the coefficients of the expansion satisfy super convergence properties. In particular, the error of the $\ell^{th}$ coefficient of the expansion scales like $\mathcal{O}(\epsilon^{2N})$ when $\ell =0$ and $\mathcal{O}(\epsilon^{2N+2-\ell})$ for all $1\leq \ell \leq N$. This result is significant, because the low-order coefficients correspond to physically relevant quantities of the underlying system. All the above estimates involve constants depending on $N$, the time $t$, and the initial condition. We investigate specifically the dependence on $N$, in order to assess whether increasing $N$ actually yields an additional factor of $\epsilon$ in the error. Numerical tests will also be presented to support the theoretical results.


title: Spectroscopy of excited states of unbound nuclei $^{30}$Ar and $^{29}$Cl

abstract: Several states of proton-unbound isotopes $^{30}$Ar and $^{29}$Cl were investigated by measuring their in-flight decay products, $^{28}$S+proton+proton and $^{28}$S+proton, respectively. A refined analysis of $^{28}$S-proton angular correlations indicates that the ground state of $^{30}$Ar is located at $2.45^{+0.05}_{-0.10}$ MeV above the two-proton emission threshold. The theoretical investigation of the $^{30}$Ar ground state decay demonstrates that its mechanism has the transition dynamics with a surprisingly strong sensitivity of the correlation patterns of the decay products to the two-proton decay energy of the $^{30}$Ar ground state and the one-proton decay energy as well as the one-proton decay width of the $^{29}$Cl ground state. The comparison of the experimental $^{28}$S-proton angular correlations with those resulting from Monte Carlo simulations of the detector response illustrates that other observed $^{30}$Ar excited states decay by sequential emission of protons via intermediate resonances in $^{29}$Cl. Based on the findings, the decay schemes of the observed states in $^{30}$Ar and $^{29}$Cl were constructed. For calibration purposes and for checking the performance of the experimental setup, decays of the previously-known states of a two-proton emitter $^{19}$Mg were remeasured. Evidences for one new excited state in $^{19}$Mg and two unknown states in $^{18}$Na were found.


title: Non-iterative SLAM for Warehouse Robots Using Ground Textures

abstract: We present a novel visual SLAM method for the warehouse robot with a single downward-facing camera using ground textures. Traditional methods resort to feature matching or point registration for pose optimization, which easily suffers from repetitive features and poor texture quality. In this paper, we present a robust kernel cross-correlator for robust image-level registration. Compared with the existing methods that often use iterative solutions, our method, named non-iterative visual SLAM (NI-SLAM), has a closed-form solution with a complexity of $O(n\log n)$. This allows it to run very efficiently, yet still provide better accuracy and robustness than the state-of-the-art methods. In the experiments, we demonstrate that it achieves 78% improvement over the state-of-the-art systems for indoor and outdoor localization. We have successfully tested it on warehouse robots equipped with a single downward camera, showcasing its product-ready superiority in a real operating area.


title: Toward Crowd-Sensitive Path Planning

abstract: If a robot can predict crowds in parts of its environment that are inaccessible to its sensors, then it can plan to avoid them. This paper proposes a fast, online algorithm that learns average crowd densities in different areas. It also describes how these densities can be incorporated into existing navigation architectures. In simulation across multiple challenging crowd scenarios, the robot reaches its target faster, travels less, and risks fewer collisions than if it were to plan with the traditional A* algorithm.


title: Monte Carlo population synthesis on massive star binaries: Astrophysical implications for gravitational wave sources

abstract: There are important but unresolved processes in the standard formation scenarios of double compact star binaries (DCBs; BH-BH, BH-NS, NS-NS systems), such as mass transfer and the common envelope (CE) phase. We analyze the effects of different assumptions on key physical processes and binary initial conditions on massive star binary evolution with binary population synthesis (BPS), including a survey of proposed prescriptions for the mass transfer ($\rm q_{\rm cr}$) and the binding energy parameter ($\lambda$) in the CE phase. We find that $\rm q_{\rm cr}$ clearly affects the properties of NS-NS systems while $\lambda$ has influence on the mass distributions of BH-BH systems. The merger rates of DCBs are increased by efficient CE ejection, which in our prescription is related to the binding energy parameter including all the possible budgets to the energy content. It has been suggested that the difference in the properties of GW150914 and GW151226 may reflect different metallicity. We reproduce their properties with our BPS calculations and find that the property of BH-BH systems at low metallicity is sensitive to $\lambda$; the efficient CE ejection leads to a top-heavy mass distribution both for the primary and secondary BHs, which is favored to explain the nature of GW150914. The efficient CE ejection also leads to enhancement of both the BH-BH and NS-NS merger rates to the level consistent with the observational constraints from the detected gravitational wave sources including GW170817.


title: Periodic Points in Genus Two: Holomorphic Sections over Hilbert Modular Varieties, Teichmuller Dynamics, and Billiards

abstract: We show that all GL(2, R)-equivariant point markings over orbit closures of primitive genus two translation surfaces arise from marking pairs of points exchanged by the hyperelliptic involution, Weierstrass points, or the golden points in the golden eigenform locus. As corollaries, we classify the holomorphically varying families of points over orbifold covers of genus two Hilbert modular surfaces, solve the finite blocking problem on genus two translation surfaces, and show that there is at most one nonarithmetic rank two orbit closure in the minimal stratum in genus four.


title: Multiple Lattice Tilings in Euclidean Spaces

abstract: This paper proves the following results: Besides parallelograms and centrally symmetric hexagons, there is no other convex domain which can form a two-, three- or four-fold lattice tiling in the Euclidean plane. If a centrally symmetric octagon can form a lattice multiple tiling, then the multiplicity is at least seven. However, there are decagons which can form five-fold $($or six-fold$)$ lattice tilings. Consequently, whenever $n\ge 3$, there are non-parallelohedral polytopes which can form five-fold lattice tilings in the $n$-dimensional Euclidean space.


title: Rank One Orbit Closures in H^{hyp}(g-1,g-1)

abstract: Every GL(2,R)-orbit in hyperelliptic components of strata of abelian differentials in genus greater than two is either closed, dense, or contained in a locus of branched covers.


title: Quenched local central limit theorem for random walks in a time-dependent balanced random environment

abstract: We prove a quenched local central limit theorem for continuous-time random walks in $\mathbb Z^d, d\ge 2$, in a uniformly-elliptic time-dependent balanced random environment which is ergodic under space-time shifts. We also obtain Gaussian upper and lower bounds for quenched and (positive and negative) moment estimates of the transition probabilities and asymptotics of the discrete Green function.


title: On stochastic and deterministic quasi-Newton methods for non-Strongly convex optimization: Asymptotic convergence and rate analysis

abstract: Motivated by applications arising from large scale optimization and machine learning, we consider stochastic quasi-Newton (SQN) methods for solving unconstrained convex optimization problems. The convergence analysis of the SQN methods, both full and limited-memory variants, require the objective function to be strongly convex. However, this assumption is fairly restrictive and does not hold for applications such as minimizing the logistic regression loss function. To the best of our knowledge, no rate statements currently exist for SQN methods in the absence of such an assumption. Also, among the existing first-order methods for addressing stochastic optimization problems with merely convex objectives, those equipped with provable convergence rates employ averaging. However, this averaging technique has a detrimental impact on inducing sparsity. Motivated by these gaps, the main contributions of the paper are as follows: (i) Addressing large scale stochastic optimization problems, we develop an iteratively regularized stochastic limited-memory BFGS (IRS-LBFGS) algorithm, where the stepsize, regularization parameter, and the Hessian inverse approximation matrix are updated iteratively. We establish the convergence to an optimal solution of the original problem both in an almost-sure and mean senses. We derive the convergence rate in terms of the objective function's values and show that it is of the order $\mathcal{O}\left(k^{-\left(\frac{1}{3}-\epsilon\right)}\right)$, where $\epsilon$ is an arbitrary small positive scalar; (ii) In deterministic regime, we show that the regularized limited-memory BFGS algorithm displays a rate of the order $\mathcal{O}\left(\frac{1}{k^{1 -\epsilon'}}\right)$, where $\epsilon'$ is an arbitrary small positive scalar. We present our numerical experiments performed on a large scale text classification problem.


title: Shannon entropies and Fisher information of K-shell electrons of neutral atoms

abstract: We represent the two K-shell electrons of neutral atoms by Hylleraas-type wave function which fulfils the exact behavior at the electron-electron and electron-nucleus coalescence points and, derive a simple method to construct expressions for single-particle position- and momentum-space charge densities, $\rho(\vec{r})$ and $\gamma(\vec{p})$ respectively. We make use of the results for $\rho(\vec{r})$ and $\gamma(\vec{p})$ to critically examine the effect of correlation on bare (uncorrelated) values of Shannon information entropies ($S$) and of Fisher information ($F$) for the K-shell electrons of atoms from helium to neon. Due to inter-electronic repulsion the values of the uncorrelated Shannon position-space entropies are augmented while those of the momentum-space entropies are reduced. The corresponding Fisher information are found to exhibit opposite behavior in respect of this. Attempts are made to provide some plausible explanation for the observed response of $S$ and $F$ to electronic correlation.


title: Reliable numerical key rates for quantum key distribution

abstract: In this work, we present a reliable, efficient, and tight numerical method for calculating key rates for finite-dimensional quantum key distribution (QKD) protocols. We illustrate our approach by finding higher key rates than those previously reported in the literature for several interesting scenarios (e.g., the Trojan-horse attack and the phase-coherent BB84 protocol). Our method will ultimately improve our ability to automate key rate calculations and, hence, to develop a user-friendly software package that could be used widely by QKD researchers.


title: The Feeling of Success: Does Touch Sensing Help Predict Grasp Outcomes?

abstract: A successful grasp requires careful balancing of the contact forces. Deducing whether a particular grasp will be successful from indirect measurements, such as vision, is therefore quite challenging, and direct sensing of contacts through touch sensing provides an appealing avenue toward more successful and consistent robotic grasping. However, in order to fully evaluate the value of touch sensing for grasp outcome prediction, we must understand how touch sensing can influence outcome prediction accuracy when combined with other modalities. Doing so using conventional model-based techniques is exceptionally difficult. In this work, we investigate the question of whether touch sensing aids in predicting grasp outcomes within a multimodal sensing framework that combines vision and touch. To that end, we collected more than 9,000 grasping trials using a two-finger gripper equipped with GelSight high-resolution tactile sensors on each finger, and evaluated visuo-tactile deep neural network models to directly predict grasp outcomes from either modality individually, and from both modalities together. Our experimental results indicate that incorporating tactile readings substantially improve grasping performance.


title: Robust Maximum Likelihood Estimation of Sparse Vector Error Correction Model

abstract: In econometrics and finance, the vector error correction model (VECM) is an important time series model for cointegration analysis, which is used to estimate the long-run equilibrium variable relationships. The traditional analysis and estimation methodologies assume the underlying Gaussian distribution but, in practice, heavy-tailed data and outliers can lead to the inapplicability of these methods. In this paper, we propose a robust model estimation method based on the Cauchy distribution to tackle this issue. In addition, sparse cointegration relations are considered to realize feature selection and dimension reduction. An efficient algorithm based on the majorization-minimization (MM) method is applied to solve the proposed nonconvex problem. The performance of this algorithm is shown through numerical simulations.


title: Centers and Azumaya loci of finite $W$-algebras

abstract: In this paper, we study the center $Z$ of the finite $W$-algebra $\mathcal{T}(\mathfrak{g},e)$ associated with a semi-simple Lie algebra $\mathfrak{g}$ over an algebraically closed field $\mathds{k}$ of characteristic $p\gg0$, and an arbitrarily given nilpotent element $e\in\mathfrak{g}$. We obtain an analogue of Veldkamp's theorem on the center. For the maximal spectrum $\text{Specm}(Z)$, we show that its Azumaya locus coincides with its smooth locus of smooth points. The former locus reflects irreducible representations of maximal dimension for $\mathcal{T}(\mathfrak{g},e)$.


title: Gamma-ray emission from internal shocks in novae

abstract: Gamma-ray emission at energies >100MeV has been detected from nine novae using the Fermi-LAT, and it can be explained by particle acceleration at shocks in these systems. Eight out of these nine objects are classical novae in which interaction of the ejecta with a tenuous circumbinary material is not expected to generate detectable gamma-ray emission. We examine whether particle acceleration at internal shocks can account for the gamma-ray emission from these novae. The shocks result from the interaction of a fast wind radiatively-driven by nuclear burning on the white dwarf with material ejected in the initial runaway stage of the nova outburst. We present a one-dimensional model for the dynamics of a forward and reverse shock system in a nova ejecta, and for the associated time-dependent particle acceleration and high-energy gamma-ray emission. Non-thermal proton and electron spectra are calculated by solving a time-dependent transport equation for particle injection, acceleration, losses, and escape from the shock region. The predicted emission is compared to LAT observations of V407 Cyg, V1324 Sco, V959 Mon, V339 Del, V1369 Cen, and V5668 Sgr. The 100MeV gamma-ray emission arises predominantly from particles accelerated up to ~100GeV at the reverse shock and undergoing hadronic interactions in the dense cooling layer downstream of the shock. The internal shock model can account for the gamma-ray emission of the novae detected by Fermi-LAT, including the main features in the observations of the recent gamma-ray nova ASASSN-16ma. Gamma-ray observations hold potential for probing the mechanism of mass ejection in novae, but should be combined to diagnostics of the thermal emission at lower energies to be more constraining. (abridged version)


title: The Structure of Root Data and Smooth Regular Embeddings of Reductive Groups

abstract: We investigate the structure of root data by considering their decomposition as a product of a semisimple root datum and a torus. Using this decomposition we obtain a parameterisation of the isomorphism classes of all root data. By working at the level of root data we introduce the notion of a smooth regular embedding of a connected reductive algebraic group, which is a refinement of the commonly used regular embeddings introduced by Lusztig. In the absence of Steinberg endomorphisms such embeddings were constructed by Benjamin Martin. In an unpublished manuscript Asai proved three key reduction techniques that are used for reducing statements about arbitrary connected reductive algebraic groups, equipped with a Frobenius endomorphism, to those whose derived subgroup is simple and simply connected. By using our investigations into root data we give new proofs of Asai's results and generalise them so that they are compatible with Steinberg endomorphisms. As an illustration of these ideas, we answer a question posed to us by Olivier Dudas concerning unipotent supports.


title: Anisotropies of Ultra-high Energy Cosmic Rays Dominated by a Single Source in the Presence of Deflections

abstract: This work presents a scenario of ultra-high energy cosmic ray source distribution where a nearby source is solely responsible for the anisotropies in arrival directions of cosmic rays while the rest of the sources contribute only isotropically. An analytical approach focused on large-scale anisotropies, which are influenced by deflections in a Kolmogorov-type turbulent magnetic field, is employed to provide more general results. When the recent Pierre Auger Observatory angular power spectrum above 8 EeV is used the restricted model gives, under the assumption of the small angle approximation, a solution where the RMS deflection with respect to the line of sight is $\alpha_{\rm rms} = \left(50^{+11}_{-10}\right)^\circ$, while the relative flux from the single source $\eta=0.03\pm 0.01$. Furthermore, the solution can be translated into constraints on the source distance, luminosity, and extra-galactic magnetic field strength. For Centaurus A and the Virgo cluster the required relation between the coherence length and the RMS magnetic field strength is obtained: a coherence length of $~\sim 100\,\mathrm{kpc}$ would imply the RMS field strength around $1\,\mathrm{nG}$ for iron dominated and above $10\,\mathrm{nG}$ for proton dominated composition. We also performed trajectory simulations with our publicly available code CRPropa to show that our analytical model can serve as a good approximation as long as the deflections in cosmic magnetic fields can be described as a random walk. The simulations showed that generally structured fields tend to suppress large-scale anisotropies, especially the dipole, compared to anisotropies at smaller scales described by higher multipoles.


title: Tilting Modules Under Special Base Changes

abstract: Given a non-unit, non-zero-divisor, central element $x$ of a ring $\Lambda$, it is well known that many properties or invariants of $\Lambda$ determine, and are determined by, those of $\Lambda / x \Lambda$ and $\Lambda_x$. In the present paper, we investigate how the property of "being tilting" behaves in this situation. It turns out that any tilting module over $\Lambda$ gives rise to tilting modules over $\Lambda_x$ and $\Lambda / x \Lambda$ after localization and passing to quotient respectively. On the other hand, it is proved that under some mild conditions, a module over $\Lambda$ is tilting if its corresponding localization and quotient are tilting over $\Lambda_x$ and $\Lambda / x \Lambda$ respectively.


title: BKTreebank: Building a Vietnamese Dependency Treebank

abstract: Dependency treebank is an important resource in any language. In this paper, we present our work on building BKTreebank, a dependency treebank for Vietnamese. Important points on designing POS tagset, dependency relations, and annotation guidelines are discussed. We describe experiments on POS tagging and dependency parsing on the treebank. Experimental results show that the treebank is a useful resource for Vietnamese language processing.


title: Entanglement Entropy of Target Functions for Image Classification and Convolutional Neural Network

abstract: The success of deep convolutional neural network (CNN) in computer vision especially image classification problems requests a new information theory for function of image, instead of image itself. In this article, after establishing a deep mathematical connection between image classification problem and quantum spin model, we propose to use entanglement entropy, a generalization of classical Boltzmann-Shannon entropy, as a powerful tool to characterize the information needed for representation of general function of image. We prove that there is a sub-volume-law bound for entanglement entropy of target functions of reasonable image classification problems. Therefore target functions of image classification only occupy a small subspace of the whole Hilbert space. As a result, a neural network with polynomial number of parameters is efficient for representation of such target functions of image. The concept of entanglement entropy can also be useful to characterize the expressive power of different neural networks. For example, we show that to maintain the same expressive power, number of channels $D$ in a convolutional neural network should scale with the number of convolution layers $n_c$ as $D\sim D_0^{\frac{1}{n_c}}$. Therefore, deeper CNN with large $n_c$ is more efficient than shallow ones.


title: On the Hybrid Minimum Principle

abstract: The Hybrid Minimum Principle (HMP) is established for the optimal control of deterministic hybrid systems with both autonomous and controlled switchings and jumps where state jumps at the switching instants are permitted to be accompanied by changes in the dimension of the state space. First order variational analysis is performed via the needle variation methodology and the necessary optimality conditions are established in the form of the HMP. A feature of special interest in this work is the explicit presentations of boundary conditions on the Hamiltonians and the adjoint processes before and after switchings and jumps. In addition to an analytic example, the HMP results are illustrated for the optimal control of an electric vehicle with transmission, where the modelling of the powertrain requires the consideration of both autonomous and controlled switchings accompanied by dimension changes.


title: Evidence for Weyl fermions in a canonical heavy-fermion semimetal YbPtBi

abstract: The manifestation of Weyl fermions in strongly correlated electron systems is of particular interest. We report evidence for Weyl fermions in the heavy fermion semimetal YbPtBi from electronic structure calculations, angle-resolved photoemission spectroscopy, magnetotransport and calorimetric measurements. At elevated temperatures where $4f$-electrons are localized, there are triply degenerate points, yielding Weyl nodes in applied magnetic fields. These are revealed by a contribution from the chiral anomaly in the magnetotransport, which at low temperatures becomes negligible due to the influence of electronic correlations. Instead, Weyl fermions are inferred from the topological Hall effect, which provides evidence for a Berry curvature, and a cubic temperature dependence of the specific heat, as expected from the linear dispersion near the Weyl nodes. The results suggest that YbPtBi is a Weyl heavy fermion semimetal, where the Kondo interaction renormalizes the bands hosting Weyl points. These findings open up an opportunity to explore the interplay between topology and strong electronic correlations.


title: An open window to neutrino mass hierarchy at maximal Dirac CP violation

abstract: Non-zero leptonic CP phases in the neutrino sector are clear evidence for physics beyond the Standard Model and have many implications in particle physics and cosmology. Some clues to maximal Dirac CP violation $\delta_{CP}=3\pi/2$ are reviewed. An approach to connect the Dirac CP violation phase and the neutrino mass hierarchy in terms of the measurement of CP asymmetry is proposed. Under appropriate requirements for baseline distance and beam energy, the problem of neutrino mass hierarchy can be revealed at maximal Dirac CP violation. General allowed parameter spaces are investigated, within which the two unknown neutrino issues, mass hierarchy and Dirac CP violation, can be solved simultaneously. The abilities of addressing the mass hierarchy of current long baseline neutrino experiments are also discussed.