The Latest ArXiv Papers (Except Weekends) of 2017-09-19

We present a new flexible Bayesian framework for directly inferring the fraction of neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR, z~6-10) from detections and non-detections of Lyman Alpha (Ly$\alpha$) emission from Lyman break galaxies (LBGs). Our framework combines sophisticated reionization simulations with empirical models of the interstellar medium (ISM) radiative transfer effects on Ly$\alpha$. We assert that the Ly$\alpha$ line profile emerging from the ISM has an important impact on the resulting transmission of photons through the IGM, and that these line profiles depend on galaxy properties. We model this effect by considering the peak velocity offset of Ly$\alpha$ lines from host galaxies' systemic redshifts which are empirically correlated with UV luminosity and redshift (or halo mass at fixed redshift). We use our framework on the sample of LBGs presented in Pentericci et al. (2014) and infer a global neutral fraction at z~7 of $\overline{x}_\mathrm{HI} = 0.59_{-0.15}^{+0.11}$, consistent with other robust probes of the EoR and confirming reionization is on-going ~700 Myr after the Big Bang. We show that using the full distribution of Ly$\alpha$ equivalent width detections and upper limits from LBGs places tighter constraints on the evolving IGM than the standard Ly$\alpha$ emitter fraction, and that larger samples are within reach of deep spectroscopic surveys of gravitationally lensed fields and JWST NIRSpec.

In this paper we investigate the formation of coalescing stellar black hole binaries (BHBs) and extreme mass ratio inspirals (EMRIs) during the assembly of a massive nucleus around an SMBH with mass $10^8$ M$_\odot$. Using direct $N$-body modeling, we show that this phase of the galaxy life is characterised by strong interactions between the SMBH and a population of massive star clusters (GCs) either formed locally or orbitally segregated due to dynamical friction. We show that if the infalling clusters contain a significant population of BHBs, the GCs-SMBH interactions can reduce the merger time-scale to $\sim 1$ Gyr. Calibrating our results to the sensitivity of the LIGO experiment, we found a detection rate of $\mathcal{R}_{\rm LIGO} \simeq 0.25-1.1$ Gpc$^{-1}$ yr$^{-1}$ for a BHB similar to the observed GW150914 source. We also investigated the rate of EMRIs caused by the GCs debris deposited around the SMBH, finding a value of $\mathcal{R}_{\rm EMRI} = 0.5-1\times 10^{-6}$ yr$^{-1}$. Finally, we show that if two or more of the infalling GCs host in their centre an intermediate mass black hole (IMBH), the three-body system composed by the two IMBHs and the SMBH causes the merging between one of the IMBH and the SMBH in more than $\sim 50\%$ of the cases within $1$ Gyr. This may give an explanation for the observed lack of IMBHs signatures in galactic nuclei.

Context. A lot of photometric data is produced by surveys such as Pan-STARRS, LONEOS, WISE or Catalina. These data are a rich source of information about the physical properties of asteroids. There are several possible approaches for utilizing these data. Lightcurve inversion is a typical method that works with individual asteroids. Our approach in this paper is statistical when we focused on large groups of asteroids like dynamical families and taxonomic classes, and the data were not sufficient for individual models.

Aims. Our aim was to study the distributions of shape elongation $b/a$ and the spin axis latitude $\beta$ for various subpopulations of asteroids and to compare our results, based on Pan-STARRS1 survey, with statistics previously done using different photometric data (Lowell database, WISE data).

Methods. We use the LEADER algorithm to compare the $b/a$ and $\beta$ distributions for different subpopulations of asteroids. The algorithm creates a cumulative distributive function (CDF) of observed brightness variations, and computes the $b/a$ and $\beta$ distributions using analytical basis functions that yield the observed CDF. A variant of LEADER is used to solve the joint distributions for synthetic populations to test the validity of the method.

Results. When comparing distributions of shape elongation for groups of asteroids with different diameters $D$, we found that there are no differences for $D < 25$ km. We also constructed distributions for asteroids with different rotation periods and revealed that the fastest rotators with $P = 0 - 4$ h are more spheroidal than the population with $P = 4 - 8$ h.

After first reviewing the gravitational wave (GW) spectral classification. we discuss the sensitivities of GW detection in space aimed at low frequency band (100 nHz-100 mHz) and middle frequency band (100 mHz-10 Hz). The science goals are to detect GWs from (i) Supermassive Black Holes; (ii) Extreme-Mass-Ratio Black Hole Inspirals; (iii) Intermediate-Mass Black Holes; (iv) Galactic Compact Binaries; (v) Stellar-Size Black Hole Binaries; and (vi) Relic GW Background. The detector proposals have arm length ranging from 100 km to 1.35x109 km (9 AU) including (a) Solar orbiting detectors and (b) Earth orbiting detectors. We discuss especially the sensitivities in the frequency band 0.1-10 microHz and the middle frequency band (0.1 Hz-10 Hz). We propose and discuss AMIGO as an Astrodynamical Middle-frequency Interferometric GW Observatory.

The Gamma-ray Cherenkov Telescope (GCT) is one of the Small Size Telescopes (SSTs) proposed for the Cherenkov Telescope Array (CTA) aimed at the 1 TeV to 300 TeV energy range. GCT will be equipped with a Compact High-Energy Camera (CHEC) containing 2048 pixels of physical size about 6$\times$6~mm$^2$, leading to a field of view of over 8 degrees. Electronics based on custom TARGET ASICs and FPGAs sample incoming signals at a gigasample per second and provide a flexible triggering scheme. Waveforms for every pixel in every event are read out are on demand without loss at over 600 events per second. A GCT prototype in Meudon, Paris saw first Cherenkov light from air showers in late 2015, using the first CHEC prototype, CHEC-M. This contribution presents results from lab and field tests with CHEC-M and the progress made to a robust camera design for deployment within CTA.

We present an unsupervised machine learning technique that automatically segments and labels galaxies in astronomical imaging surveys using only pixel data. Distinct from previous unsupervised machine learning approaches used in astronomy we use no pre-selection or pre-filtering of target galaxy type to identify galaxies that are similar. We demonstrate the technique on the HST Frontier Fields. By training the algorithm using galaxies from one field (Abell 2744) and applying the result to another (MACS0416.1-2403), we show how the algorithm can cleanly separate early and late type galaxies without any form of pre-directed training for what an 'early' or 'late' type galaxy is. We then apply the technique to the HST CANDELS fields, creating a catalogue of approximately 60,000 classifications. We show how the automatic classification groups galaxies of similar morphological (and photometric) type, and make the classifications public via a catalogue, a visual catalogue and galaxy similarity search. We compare the CANDELS machine-based classifications to human-based classifications from the Galaxy Zoo: CANDELS project. Although there is not a direct mapping between Galaxy Zoo and our hierarchical labelling, we demonstrate a good level of concordance between human and machine classifications. Finally, we show how the technique can be used to identify rarer objects and present new lensed galaxy candidates from the CANDELS imaging.

Telescopes based on the imaging atmospheric Cherenkov technique (IACTs) detect images of the atmospheric showers generated by gamma rays and cosmic rays as they are absorbed by the atmosphere. The much more frequent cosmic-ray events form the main background when looking for gamma-ray sources, and therefore IACT sensitivity is significantly driven by the capability to distinguish between these two types of events. Supervised learning algorithms, like random forests and boosted decision trees, have been shown to effectively classify IACT events. In this contribution we present results from exploratory work using deep learning as an event classification method for the Cherenkov Telescope Array (CTA). CTA, conceived as an array of tens of IACTs, is an international project for a next-generation ground-based gamma-ray observatory, aiming to improve on the sensitivity of current-generation experiments by an order of magnitude and provide energy coverage from 20 GeV to more than 300 TeV.

We develop a method to forecast the outcome of the LHC Run 3 based on the hypothetical detection of $\mathcal{O}(100)$ signal events at XENONnT. Our method relies on a systematic classification of renormaLISAble single-mediator models for dark matter-quark interactions, and is valid for dark matter candidates of spin less than or equal to one. Applying our method to simulated data, we find that at the end of the LHC Run 3 only two mutually exclusive scenarios would be compatible with the detection of $\mathcal{O}(100)$ signal events at XENONnT. In a first scenario, the energy distribution of the signal events is featureless, as for canonical spin-independent interactions. In this case, if a mono-jet signal is detected at the LHC, dark matter must have spin 1/2 and interact with nucleons through a unique velocity-dependent operator. If a mono-jet signal is not detected, dark matter interacts with nucleons through canonical spin-independent interactions. In a second scenario, the spectral distribution of the signal events exhibits a bump at non zero recoil energies. In this second case, a mono-jet signal must be detected at the LHC Run 3, dark matter must have spin 1/2 and interact with nucleons through a unique momentum-dependent operator. We therefore conclude that the observation of $\mathcal{O}(100)$ signal events at XENONnT combined with the detection, or the lack of detection, of a mono-jet signal at the LHC Run 3 would significantly narrow the range of possible dark matter-nucleon interactions. As we argued above, it can also provide key information on the dark matter particle spin.

The stability properties of rotating relativistic stars against prompt gravitational collapse to a black hole are rather well understood for uniformly rotating models. This is not the case for differentially rotating neutron stars, which are expected to be produced in catastrophic events such as the merger of binary system of neutron stars or the collapse of a massive stellar core. Combining sequences of differentially rotating models with their dynamical evolution, we show that a sufficient stability criterion for differentially rotating neutron stars exists similar to the one of their uniformly rotating counterparts. Namely: along a sequence of constant angular momentum, a dynamical instability sets in for central rest-mass densities slightly below the one of the equilibrium solution at the turning point. In addition, following Breu & Rezzolla (2016), we show that "quasi-universal" relations can be found when calculating the turning-point mass. In turn, this allows us to compute the maximum mass allowed by differential rotation, $M_{\rm max,dr}$, in terms of the maximum mass of the nonrotating configuration, $M_{_{\rm TOV}}$, finding that $M_{\rm max, dr} \simeq \left(1.54 \pm 0.05\right) M_{_{\rm TOV}}$ for all the equations of state we have considered.

We present optical polarimetry in the period May-August 2017 of the enigmatic "dipping" star KIC 8462852. During that period three ~1% photometric dips were reported by other observers. We measured the average absolute polarization of the source, and find no excess or unusual polarization compared to a nearby comparison star. We place tight upper limits on any change in the degree of polarization of the source between epochs in- and out-of-dip of <0.1% (8500-Ang.) and <0.2% (7050-Ang. and 5300-Ang.). How our limits are interpreted depends on the specific model being considered. If the whole stellar disk were covered by material with an optical depth of ~0.01 then the fractional polarisation introduced by this material must be less than 10-20%. While our non-detection does not constrain the comet scenario, it predicts that even modest amounts of dust that have properties similar to Solar System comets may be detectable. We note that the sensitivity of our method scales with the depth of the dip. Should a future ~20% photometric dip be observed (as was previously detected by Kepler) our method would constrain any induced polarization associated with any occulting material to 0.5-1.0%.

We present a demonstration of delensing the observed cosmic microwave background (CMB) B-mode polarization anisotropy. This process of reducing the gravitational-lensing generated B-mode component will become increasingly important for improving searches for the B modes produced by primordial gravitational waves. In this work, we delens B-mode maps constructed from multi-frequency SPTpol observations of a 90 deg$^2$ patch of sky by subtracting a B-mode template constructed from two inputs: SPTpol E-mode maps and a lensing potential map estimated from the $\textit{Herschel}$ $500\,\mu m$ map of the CIB. We find that our delensing procedure reduces the measured B-mode power spectrum by 28% in the multipole range $300 < \ell < 2300$; this is shown to be consistent with expectations from theory and simulations and to be robust against systematics. The null hypothesis of no delensing is rejected at $6.9 \sigma$. Furthermore, we build and use a suite of realistic simulations to study the general properties of the delensing process and find that the delensing efficiency achieved in this work is limited primarily by the noise in the lensing potential map. We demonstrate the importance of including realistic experimental non-idealities in the delensing forecasts used to inform instrument and survey-strategy planning of upcoming lower-noise experiments, such as CMB-S4.

pulsar timing and gravitational-wave (GW) detectors are superb laboratories to study gravity theories in the strong-field regime. Here we combine those tools to test the mono-scalar-tensor theory of Damour and Esposito-Far{\e}se (DEF), which predicts nonperturbative scalarization phenomena for neutron stars (NSs). First, applying Markov-chain Monte Carlo techniques, we use the absence of dipolar radiation in the pulsar-timing observations of five binary systems composed of a NS and a white dwarf, and eleven equations of state (EOSs) for NSs, to derive the most stringent constraints on the two free parameters of the DEF scalar-tensor theory. Since the binary-pulsar bounds depend on the NS mass and the EOS, we find that current pulsar-timing observations leave scalarization windows, i.e., regions of parameter space where scalarization can still be prominent. Then, we investigate if these scalarization windows could be closed and if pulsar-timing constraints could be improved by laser-interferometer GW detectors, when spontaneous (or dynamical) scalarization sets in during the early (or late) stages of a binary NS (BNS) evolution. For the early inspiral of a BNS carrying constant scalar charge, we employ a Fisher matrix analysis to show that Advanced LIGO can improve pulsar-timing constraints for some EOSs, and next-generation detectors, such as the Cosmic Explorer and Einstein Telescope, will be able to improve those bounds for all eleven EOSs. Using the late inspiral of a BNS, we estimate that for some of the EOSs under consideration the onset of dynamical scalarization can happen early enough to improve the constraints on the DEF parameters obtained by combining the five binary pulsars. Thus, in the near future the complementarity of pulsar timing and direct observations of GWs on the ground will be extremely valuable in probing gravity theories in the strong-field regime.

Planets in the mass range from 2 to 15 MEarth are very diverse. Some of them have low densities, while others are very dense. By measuring the masses and radii, the mean densities, structure, and composition of the planets are constrained. These parameters also give us important information about their formation and evolution, and about possible processes for atmospheric loss.We determined the masses, radii, and mean densities for the two transiting planets orbiting K2-106. The inner planet has an ultra-short period of 0.57 days. The period of the outer planet is 13.3 days.Although the two planets have similar masses, their densities are very different. For K2-106b we derive Mb=8.36-0.94+0.96 MEarh, Rb=1.52+/-0.16 REarth, and a high density of 13.1-3.6+5.4 gcm-3. For K2-106c, we find Mc=5.8-3.0+3.3 MEarth, Rc=2.50-0.26+0.27 REarth and a relatively low density of 2.0-1.1+1.6 gcm-3.Since the system contains two planets of almost the same mass, but different distances from the host star, it is an excellent laboratory to study atmospheric escape. In agreement with the theory of atmospheric-loss processes, it is likely that the outer planet has a hydrogen-dominated atmosphere. The mass and radius of the inner planet is in agreement with theoretical models predicting an iron core containing 80+20-30% of its mass. Such a high metal content is surprising, particularly given that the star has an ordinary (solar) metal abundance. We discuss various possible formation scenarios for this unusual planet.

Current and upcoming radio interferometric experiments are aiming to make a statistical characterization of the high-redshift 21cm fluctuation signal spanning the hydrogen reionization and X-ray heating epochs of the universe. However, connecting 21cm statistics to underlying physical parameters is complicated by the theoretical challenge of modeling the relevant physics at computational speeds quick enough to enable exploration of the high dimensional and weakly constrained parameter space. In this work, we use machine learning algorithms to build a fast emulator that mimics expensive simulations of the 21cm signal across a wide parameter space to high precision. We embed our emulator within a Markov-Chain Monte Carlo framework, enabling it to explore the posterior distribution over a large number of model parameters, including those that govern the Epoch of Reionization, the Epoch of X-ray Heating, and cosmology. As a worked example, we use our emulator to present an updated parameter constraint forecast for the Hydrogen Epoch of Reionization Array experiment, showing that its characterization of a fiducial 21cm power spectrum will considerably narrow the allowed parameter space of reionization and heating parameters, and could help strengthen Planck's constraints on $\sigma_8$. We provide both our generalized emulator code and its implementation specifically for 21cm parameter constraints as publicly available software.

Predicting the Presence of Companions for Stripped-Envelope Supernovae: The Case of the Broad-Lined Type Ic SN 2002ap. (arXiv:1705.07898v2 [astro-ph.HE] UPDATED)

Many young, massive stars are found in close binaries. Using population synthesis simulations we predict the likelihood of a companion star being present when these massive stars end their lives as core-collapse supernovae (SNe). We focus on stripped-envelope SNe, whose progenitors have lost their outer hydrogen and possibly helium layers before explosion. We use these results to interpret new Hubble Space Telescope observations of the site of the broad-lined Type Ic SN 2002ap, 14 years post-explosion. For a subsolar metallicity consistent with SN 2002ap, we expect a main-sequence companion present in about two thirds of all stripped-envelope SNe and a compact companion (likely a stripped helium star or a white dwarf/neutron star/black hole) in about 5% of cases. About a quarter of progenitors are single at explosion (originating from initially single stars, mergers or disrupted systems). All the latter scenarios require a massive progenitor, inconsistent with earlier studies of SN 2002ap. Our new, deeper upper limits exclude the presence of a main-sequence companion star $>8$-$10$ Msun, ruling out about 40% of all stripped-envelope SN channels. The most likely scenario for SN 2002ap includes nonconservative binary interaction of a primary star initially $\lesssim 23$ Msun. Although unlikely ($<$1% of the scenarios), we also discuss the possibility of an exotic reverse merger channel for broad-lined Type Ic events. Finally, we explore how our results depend on the metallicity and the model assumptions and discuss how additional searches for companions can constrain the physics that governs the evolution of SN progenitors.

The merger of stellar-mass black holes (BHs) is not expected to generate detectable electromagnetic (EM) emission. However, recently the gravitational wave (GW) events GW150914 and GW170104, detected by LIGO to be the result of merging, ~60 solar mass black hole binaries (BHBs), each have claimed coincident gamma-ray emission. Motivated by the intriguing possibility of an EM counterpart to BHB mergers, we construct a model that can reproduce the observed EM and GW signals for GW150914- and GW170104-like events, from a single-star progenitor. Following Loeb (2016), we envision a massive, rapidly rotating star within which a rotating bar instability fractures the core into two overdensities that fragment into clumps which merge to form BHs in a tight binary with arbitrary spin-orbit alignment. Once formed, the BHB inspirals due to gas and gravitational-wave drag until tidal forces trigger strong feeding of the BHs with the surrounding stellar-density gas about 10 seconds before merger. The resulting giga-Eddington accretion peak launches a jet that breaks out of the progenitor star and drives a powerful outflow that clears the gas from the orbit of the binary within one second, preserving the vacuum GW waveform in the LIGO band. The single-progenitor scenario predicts the existence of variability of the gamma-ray burst, modulated at the ~0.2 second chirping period of the BHB due to relativistic Doppler boost. The jet breakout should be accompanied by a low-luminosity supernova. Finally, because the BHBs of the single progenitor model do not exist at large separations, they will not be detectable in the low frequency gravitational wave band of LISA. Hence, the single-progenitor BHBs will be unambiguously discernible from BHBs formed through alternate, double-progenitor evolution scenarios.

Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most $\sim 0.02\;M_{\odot}$. We investigate DNS masses, spins and velocities, and in particular correlations between spin period, orbital period and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes which may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short $\gamma$-ray burst.

Ultraluminous X-ray sources (ULX) are generally assumed to host accreting black holes some of which are intermediate mass black holes (IMBH). There is recent understanding that some ULX host accreting neutron stars, i.e. pulsating ultraluminous X-ray sources (PULX). Another subgroup of ULX, ultra-luminous super-soft sources (ULS), have extremely soft X-ray emission ($kT<0.1$ keV) generally considered as being due to the reprocessing of X-rays in optically thick outflow of matter from the surrounding disc. This motivates a connection between ULSs and the general population of ULXs in terms of viewing angle. Here we propose that the bulk of the ULS/ULX population host rapidly spinning highly magnetized, $B\sim10^{13}$G, neutron stars in supercritical propeller stage. The spindown power of the neutron stars dominates luminosity of these sources. We argue that PULXs descend from ULX/ULS' when the neutron star spins-down to a period allowing for accretion to commence. We rekindle the earlier view that the galactic system SS433 is a supercritical propeller and hence is a ULX. We finally speculate that the supernova debris disc detected around the magnetar 4U~0142$+$61, at its earliest stage, would cause the system to appear as ULS or ULX, depending on the viewing angle, if observed from another nearby galaxy.

We investigate current and future prospects for coincident detection of high-energy neutrinos and gravitational waves (GWs). Short gamma-ray bursts (SGRBs) are believed to originate from mergers of compact star binaries involving neutron stars. We estimate high-energy neutrino fluences from prompt emission, extended emission, X-ray flares, and plateau emission, and show that neutrino signals associated with the extended emission are the most promising. Assuming that the cosmic-ray loading factor is $\sim10$ and the Lorentz factor distribution is lognormal, we calculate the probability of neutrino detection from extended emission by current and future neutrino detectors, and find that the quasi-simultaneous detection of high-energy neutrinos, gamma rays, and GWs is possible with future instruments or even with current instruments for nearby SGRBs having extended emission. We also discuss stacking analyses that will also be useful with future experiments such as IceCube-Gen2.

We review the measurements of magnetic fields of OBA stars. Based on these data we confirm that magnetic fields are distributed according to a lognormal law with a mean log(B)=-0.5 (B in kG) with a standard deviation sigma=0.5. The shape of the magnetic field distribution is similar to that for neutron stars. This finding is in favor of the hypothesis that the magnetic field of a neutron star is determined mainly by the magnetic field of its predecessor, the massive OB star. Further, we model the evolution of an ensemble of magnetic massive stars in the Galaxy. We use our own population synthesis code to obtain the distribution of stellar radii, ages, masses, temperatures, effective magnetic fields and magnetic fluxes from the pre-main sequence (PMS) via zero age main sequence (ZAMS) up to the terminal age main sequence (TAMS) stages. A comparison of the obtained in our model magnetic field distribution (MFD) with that obtained from the recent measurements of the stellar magnetic field allows us to conclude that the evolution of magnetic fields of massive stars is slow if not absent. The shape of the real MFD shows no indications of the magnetic desert proposed previously. Based on this finding we argue that the observed fraction of magnetic stars is determined by physical conditions at the PMS stage of stellar evolution.

Whilst astronomy as a science is historically founded on observations at optical wavelengths, studying the Universe in other bands has yielded remarkable discoveries, from pulsars in the radio, signatures of the Big Bang at submm wavelengths, through to high energy emission from accreting, gravitationally-compact objects and the discovery of gamma-ray bursts. Unsurprisingly, the result of combining multiple wavebands leads to an enormous increase in diagnostic power, but powerful insights can be lost when the sources studied vary on timescales shorter than the temporal separation between observations in different bands. In July 2015, the workshop "Paving the way to simultaneous multi-wavelength astronomy" was held as a concerted effort to address this at the Lorentz Center, Leiden. It was attended by 50 astronomers from diverse fields as well as the directors and staff of observatories and spaced-based missions. This community white paper has been written with the goal of disseminating the findings of that workshop by providing a concise review of the field of multi-wavelength astronomy covering a wide range of important source classes, the problems associated with their study and the solutions we believe need to be implemented for the future of observational astronomy. We hope that this paper will both stimulate further discussion and raise overall awareness within the community of the issues faced in a developing, important field.

After first reviewing the gravitational wave (GW) spectral classification. we discuss the sensitivities of GW detection in space aimed at low frequency band (100 nHz-100 mHz) and middle frequency band (100 mHz-10 Hz). The science goals are to detect GWs from (i) Supermassive Black Holes; (ii) Extreme-Mass-Ratio Black Hole Inspirals; (iii) Intermediate-Mass Black Holes; (iv) Galactic Compact Binaries; (v) Stellar-Size Black Hole Binaries; and (vi) Relic GW Background. The detector proposals have arm length ranging from 100 km to 1.35x109 km (9 AU) including (a) Solar orbiting detectors and (b) Earth orbiting detectors. We discuss especially the sensitivities in the frequency band 0.1-10 microHz and the middle frequency band (0.1 Hz-10 Hz). We propose and discuss AMIGO as an Astrodynamical Middle-frequency Interferometric GW Observatory.

We quantise and solve the dynamics of gravitational waves in a quantum Friedmann-Lemaitre-Robertson-Walker spacetime. The classical model is formulated canonically and includes a perfect fluid. The Hamiltonian constraint is de-parametrised by setting a fluid variable as the internal clock. The obtained reduced (or, physical) phase space is then quantised. Quantisation of perturbation variables is canonical, whereas quantisation of background variables is based on the affine group as the phase space symmetry. The initial singularity is removed and replaced with a quantum bounce. The quantum model depends on a free parameter that is naturally induced from quantisation and determines the scale of the bounce. We study the dynamics of the quantised gravitational waves across the bounce by three different methods ("thin-horizon", analytical and numerical) which give consistent results and we determine the primordial power spectrum for the case of radiation-dominated universe. Next, we use the instantaneous radiation-matter transition transfer function to make approximate predictions for late universe and constrain our model with the LIGO and Planck data. We also give an estimate of the quantum uncertainties in the present-day universe.

We put forward new explicit reaLISAtions of dS/CFT that relate ${\cal N}=2$ supersymmetric Euclidean vector models with reversed spin-statistics in three dimensions to specific supersymmetric Vasiliev theories in four-dimensional de Sitter space. The partition function of the free supersymmetric vector model deformed by a range of low spin deformations that preserve supersymmetry appears to specify a well-defined wave function with asymptotic de Sitter boundary conditions in the bulk. In particular we find the wave function is globally peaked at undeformed de Sitter space, with a low amplitude for strong deformations. This suggests that supersymmetric de Sitter space is stable in higher-spin gravity and in particular free from ghosts. We speculate this is a limiting case of the de Sitter realizations in exotic string theories.

Matched-filter searches for gravitational waves from coalescing compact binaries by the LIGO Scientific Collaboration use the FINDCHIRP algorithm: an implementation of the optimal filter with innovations to account for unknown signal parameters and to improve performance on detector data that has nonstationary and non-Gaussian artifacts. We provide details on the FINDCHIRP algorithm as used in the search for subsolar mass binaries, binary neutron stars, neutron star-black hole binaries, and binary black holes.

We study generaLISAtions of ghost-free bimetric theory which involve more than two spin-2 fields. The consistent interactions can enter in the form of two different couplings and in the majority of this work we concentrate on the simpler one. The corresponding action involves one metric coupled to $N$ tensor fields which do not interact with each other. We derive maximally symmetric solutions to the multimetric equations of motion and identify the mass eigenstates in the linearised theory around these backgrounds. Our results are then applied to the problem of singling out multimetric models which possess certain additional structures. In particular, we look for a relation between scale invariant background solutions, the perturbative emergence of Weyl invariance and the presence of partially massless spin-2 fields in the linear theory. Our findings generalise known results in bimetric theory and allow us to point out similarities and differences between the bi- and multimetric models.

pulsar timing and gravitational-wave (GW) detectors are superb laboratories to study gravity theories in the strong-field regime. Here we combine those tools to test the mono-scalar-tensor theory of Damour and Esposito-Far{\e}se (DEF), which predicts nonperturbative scalarization phenomena for neutron stars (NSs). First, applying Markov-chain Monte Carlo techniques, we use the absence of dipolar radiation in the pulsar-timing observations of five binary systems composed of a NS and a white dwarf, and eleven equations of state (EOSs) for NSs, to derive the most stringent constraints on the two free parameters of the DEF scalar-tensor theory. Since the binary-pulsar bounds depend on the NS mass and the EOS, we find that current pulsar-timing observations leave scalarization windows, i.e., regions of parameter space where scalarization can still be prominent. Then, we investigate if these scalarization windows could be closed and if pulsar-timing constraints could be improved by laser-interferometer GW detectors, when spontaneous (or dynamical) scalarization sets in during the early (or late) stages of a binary NS (BNS) evolution. For the early inspiral of a BNS carrying constant scalar charge, we employ a Fisher matrix analysis to show that Advanced LIGO can improve pulsar-timing constraints for some EOSs, and next-generation detectors, such as the Cosmic Explorer and Einstein Telescope, will be able to improve those bounds for all eleven EOSs. Using the late inspiral of a BNS, we estimate that for some of the EOSs under consideration the onset of dynamical scalarization can happen early enough to improve the constraints on the DEF parameters obtained by combining the five binary pulsars. Thus, in the near future the complementarity of pulsar timing and direct observations of GWs on the ground will be extremely valuable in probing gravity theories in the strong-field regime.

In grid-based codes that provide the combined solution of the Einstein equations and of relativistic hydrodynamics, the history of the fluid is not simple to track, especially when compared with particle-based codes. The use of tracers, namely massless particles that are advected with the flow, represents a simple and effective way to solve this problem. Yet, the use of tracers in numerical relativity is far from being settled and several issues, such as the impact of different placements in time and space of the tracers, or the relation between the placement and the description of the underlying fluid, have not yet been addressed. In this paper we present the first detailed discussion of the use tracers in numerical-relativity simulations focussing on both unbound material -- such as the one leading to r-process nucleosynthesis in binary mergers of neutron stars -- and on bound material -- such as the one in the core of the object produced from the merger of two neutron stars. In particular, when interested in unbound matter, we have evaluated four different placement schemes that could be used to initially distribute the tracers and how well their predictions match those obtained when using information from the actual fluid flow. Countering our naive expectations, we found that the most effective method does not rely on the rest-mass density distribution nor on the fluid that is unbound, but simply distributes tracers uniformly in rest-mass density. This prescription leads to the closest matching with the information obtained from the hydrodynamical solution. When considering bound matter, we demonstrate that tracers can provide insight into the fine details of the fluid motion as they can be used to track the evolution of fluid elements or to calculate the variation of quantities that are conserved along streamlines of adiabatic flows.

The basic tenet of the present work is the assumption of the lack of external and fixed time in the Universe. This assumption is best embodied by general relativity, which replaces the fixed space-time structure with the gravitational field, which is subject to dynamics. The lack of time does not imply the lack of evolution but rather brings to the forefront the role of internal clocks which are some largely arbitrary internal degrees of freedom with respect to which the evolution of timeless systems can be described. We take this idea seriously and try to understand what it implies for quantum mechanics when the fixed external time is replaced by an arbitrary internal clock. We put the issue in a solid, mathematically rigorous framework. We find that the dynamical interpretation of a quantum state of a timeless system depends on the employed internal clock. In particular, we find that the continuous spectra of well-known dynamical observables like the position of a free particle on the real line may turn discrete if measured in unusual clocks. We discuss the meaning of our result for attempts at quantization of global gravitational degrees of freedom.

Relativistic ultracompact objects without an event horizon may be able to form in nature and merge as binary systems, mimicking the coalescence of ordinary black holes. The postmerger phase of such processes presents characteristic signatures, which appear as repeated pulses within the emitted gravitational waveform, i.e., echoes with variable amplitudes and frequencies. Future detections of these signals can shed new light on the existence of horizonless geometries, and provide new information on the nature of gravity in a genuine strong-field regime. In this work we analyze phenomenological templates used to characterize echolike structures produced by exotic compact objects, and we investigate for the first time the ability of current and future interferometers to constrain their parameters. Using different models with an increasing level of accuracy, we determine the features that can be measured with the largest precision, and we span the parameter space to find the most favorable configurations to be detected. Our analysis shows that current detectors may already be able to extract all the parameters of the echoes with good accuracy, and that multiple interferometers can measure frequencies and damping factors of the signals at the level of percent.

We investigate current and future prospects for coincident detection of high-energy neutrinos and gravitational waves (GWs). Short gamma-ray bursts (SGRBs) are believed to originate from mergers of compact star binaries involving neutron stars. We estimate high-energy neutrino fluences from prompt emission, extended emission, X-ray flares, and plateau emission, and show that neutrino signals associated with the extended emission are the most promising. Assuming that the cosmic-ray loading factor is $\sim10$ and the Lorentz factor distribution is lognormal, we calculate the probability of neutrino detection from extended emission by current and future neutrino detectors, and find that the quasi-simultaneous detection of high-energy neutrinos, gamma rays, and GWs is possible with future instruments or even with current instruments for nearby SGRBs having extended emission. We also discuss stacking analyses that will also be useful with future experiments such as IceCube-Gen2.

Electromagnetic Radiation Accompanying Gravitational Waves from Black Hole Binaries. (arXiv:1706.05519v2 [astro-ph.HE] CROSS LISTED)

The transformation of powerful gravitational waves, created by the coalescence of massive black hole binaries, into electromagnetic radiation in external magnetic fields is revisited. In contrast to the previous calculations of the similar effect, we study the realistic case of the gravitational radiation frequency below the plasma frequency of the surrounding medium. The gravitational waves propagating in the plasma constantly create electromagnetic radiation dragging it with them, despite the low frequency. The plasma heating by the unattenuated electromagnetic wave may be significant in a hot rarefied plasma with strong magnetic field and can lead to a noticeable burst of electromagnetic radiation with higher frequency. The graviton-to-photon conversion effect in plasma is discussed in the context of possible electromagnetic counterparts of GW150914 and GW170104.

Charged track multiplicity is among the most powerful observables for discriminating quark- from gluon-initiated jets. Despite its utility, it is not infrared and collinear (IRC) safe, so perturbative calculations are limited to studying the energy evolution of multiplicity moments. While IRC-safe observables, like jet mass, are perturbatively calculable, their distributions often exhibit Casimir scaling, such that their quark/gluon discrimination power is limited by the ratio of quark to gluon color factors. In this paper, we introduce new IRC-safe counting observables whose discrimination performance exceeds that of jet mass and approaches that of track multiplicity. The key observation is that track multiplicity is approximately Poisson distributed, with more suppressed tails than the Sudakov peak structure from jet mass. By using an iterated version of the soft drop jet grooming algorithm, we can define a "soft drop multiplicity" which is Poisson distributed at leading-logarithmic accuracy. In addition, we calculate the next-to-leading-logarithmic corrections to this Poisson structure. If we allow the soft drop groomer to proceed to the end of the jet branching history, we can define a collinear-unsafe (but still infrared-safe) counting observable. Exploiting the universality of the collinear limit, we define generalized fragmentation functions to study the perturbative energy evolution of collinear-unsafe multiplicity.

In a recent paper we studied the effect of new-physics operators with different Lorentz structures on the semileptonic $\Lambda_b \to \Lambda_c \tau \bar{\nu}_{\tau}$ decay. This decay is of interest in light of the $R({D^{(*)}})$ puzzle in the semileptonic $\bar{B} \to D^{(*)} \tau {\bar\nu}_\tau$ decays. In this work we add tensor operators to extend our previous results and consider both model-independent new physics (NP) and specific classes of models proposed to address the $R({D^{(*)}})$ puzzle. We show that a measurement of $R(\Lambda_c) = {\cal B}[\Lambda_b \to \Lambda_c \tau \bar{\nu}_{\tau}] / {\cal B}[\Lambda_b \to \Lambda_c \ell \bar{\nu}_{\ell}]$ can strongly constrain the NP parameters of models discussed for the $R({D^{(*)}})$ puzzle. We use form factors from lattice QCD to calculate all $\Lambda_b \to \Lambda_c \tau \bar{\nu}_{\tau}$ observables. The $\Lambda_b \to \Lambda_c$ tensor form factors had not previously been determined in lattice QCD, and we present new lattice results for these form factors here.

Phase Transitions of a (Super) Quantum Mechanical Matrix Model with a Chemical Potential. (arXiv:1707.02898v2 [hep-th] UPDATED)

In this paper, we study the finite-temperature matrix quantum mechanics with chemical potential term linear in the single trace of U(N) matrices, via Monte Carlo simulation. In the bosonic case, we exhibit the existence of the Gross-Witten-Wadia (GWW) type third-order phase transition. We also extend our studies to the model with the fermionic degrees of freedom employing the non-lattice simulation via Fourier expansion, and explore the possibilities that there is a phase transition between the gapped and ungapped phase both in the absence and presence of the chemical potential term. We make a comparison of the phase diagram between the bosonic and fermionic cases.

Among various scenarios of baryon asymmetry of the Universe, electroweak baryogenesis is directly connected with physics of the Higgs sector. We discuss spectra of gravitational waves which are originated by the strongly first order phase transition at the electroweak symmetry breaking, which is required for a successful scenario of electroweak baryogenesis. In the $Z_3$ symmetric singlet scalar model, the significant gravitational waves are caused by the multi-step phase transition. We show that the model can be tested by measuring the characteristic spectra of the gravitational waves at future interferometers such as LISA and DECIGO.

We revisit the phase diagram of strong-interaction matter for the two-flavor quark-meson model using the Functional Renormalization Group. In contrast to standard mean-field calculations, an unusual phase structure is encountered at low temperatures and large quark chemical potentials. In particular, we identify a regime where the pressure decreases with increasing temperature and discuss possible reasons for this unphysical behavior.

We develop a method to forecast the outcome of the LHC Run 3 based on the hypothetical detection of $\mathcal{O}(100)$ signal events at XENONnT. Our method relies on a systematic classification of renormaLISAble single-mediator models for dark matter-quark interactions, and is valid for dark matter candidates of spin less than or equal to one. Applying our method to simulated data, we find that at the end of the LHC Run 3 only two mutually exclusive scenarios would be compatible with the detection of $\mathcal{O}(100)$ signal events at XENONnT. In a first scenario, the energy distribution of the signal events is featureless, as for canonical spin-independent interactions. In this case, if a mono-jet signal is detected at the LHC, dark matter must have spin 1/2 and interact with nucleons through a unique velocity-dependent operator. If a mono-jet signal is not detected, dark matter interacts with nucleons through canonical spin-independent interactions. In a second scenario, the spectral distribution of the signal events exhibits a bump at non zero recoil energies. In this second case, a mono-jet signal must be detected at the LHC Run 3, dark matter must have spin 1/2 and interact with nucleons through a unique momentum-dependent operator. We therefore conclude that the observation of $\mathcal{O}(100)$ signal events at XENONnT combined with the detection, or the lack of detection, of a mono-jet signal at the LHC Run 3 would significantly narrow the range of possible dark matter-nucleon interactions. As we argued above, it can also provide key information on the dark matter particle spin.

In a recent paper we studied the effect of new-physics operators with different Lorentz structures on the semileptonic $\Lambda_b \to \Lambda_c \tau \bar{\nu}_{\tau}$ decay. This decay is of interest in light of the $R({D^{(*)}})$ puzzle in the semileptonic $\bar{B} \to D^{(*)} \tau {\bar\nu}_\tau$ decays. In this work we add tensor operators to extend our previous results and consider both model-independent new physics (NP) and specific classes of models proposed to address the $R({D^{(*)}})$ puzzle. We show that a measurement of $R(\Lambda_c) = {\cal B}[\Lambda_b \to \Lambda_c \tau \bar{\nu}_{\tau}] / {\cal B}[\Lambda_b \to \Lambda_c \ell \bar{\nu}_{\ell}]$ can strongly constrain the NP parameters of models discussed for the $R({D^{(*)}})$ puzzle. We use form factors from lattice QCD to calculate all $\Lambda_b \to \Lambda_c \tau \bar{\nu}_{\tau}$ observables. The $\Lambda_b \to \Lambda_c$ tensor form factors had not previously been determined in lattice QCD, and we present new lattice results for these form factors here.

Charged track multiplicity is among the most powerful observables for discriminating quark- from gluon-initiated jets. Despite its utility, it is not infrared and collinear (IRC) safe, so perturbative calculations are limited to studying the energy evolution of multiplicity moments. While IRC-safe observables, like jet mass, are perturbatively calculable, their distributions often exhibit Casimir scaling, such that their quark/gluon discrimination power is limited by the ratio of quark to gluon color factors. In this paper, we introduce new IRC-safe counting observables whose discrimination performance exceeds that of jet mass and approaches that of track multiplicity. The key observation is that track multiplicity is approximately Poisson distributed, with more suppressed tails than the Sudakov peak structure from jet mass. By using an iterated version of the soft drop jet grooming algorithm, we can define a "soft drop multiplicity" which is Poisson distributed at leading-logarithmic accuracy. In addition, we calculate the next-to-leading-logarithmic corrections to this Poisson structure. If we allow the soft drop groomer to proceed to the end of the jet branching history, we can define a collinear-unsafe (but still infrared-safe) counting observable. Exploiting the universality of the collinear limit, we define generalized fragmentation functions to study the perturbative energy evolution of collinear-unsafe multiplicity.

We present a parameter-free scheme to combine fixed-order multi-jet results with parton-shower evolution. The scheme produces jet cross sections with leading-order accuracy in the complete phase space of multiple emissions, resumming large logarithms when appropriate, while not arbitrarily enforcing ordering on momentum configurations beyond the reach of the parton-shower evolution equation. This requires the development of a matrix-element correction scheme for complex phase-spaces including ordering conditions as well as a systematic scale-setting procedure for unordered phase-space points. The resulting algorithm does not require a merging-scale parameter. We implement the new method in the Vincia framework and compare to LHC data.

We investigate current and future prospects for coincident detection of high-energy neutrinos and gravitational waves (GWs). Short gamma-ray bursts (SGRBs) are believed to originate from mergers of compact star binaries involving neutron stars. We estimate high-energy neutrino fluences from prompt emission, extended emission, X-ray flares, and plateau emission, and show that neutrino signals associated with the extended emission are the most promising. Assuming that the cosmic-ray loading factor is $\sim10$ and the Lorentz factor distribution is lognormal, we calculate the probability of neutrino detection from extended emission by current and future neutrino detectors, and find that the quasi-simultaneous detection of high-energy neutrinos, gamma rays, and GWs is possible with future instruments or even with current instruments for nearby SGRBs having extended emission. We also discuss stacking analyses that will also be useful with future experiments such as IceCube-Gen2.

We treat the topic of the closures of the nilpotent orbits of the Lie algebras of Exceptional groups through their descriptions as moduli spaces, in terms of Hilbert series and the highest weight generating functions for their representation content. We extend the set of known Coulomb branch quiver theory constructions for Exceptional group minimal nilpotent orbits, or reduced single instanton moduli spaces, to include all orbits of Characteristic Height 2, drawing on extended Dynkin diagrams and the unitary monopole formula. We also present a representation theoretic formula, based on locaLISAtion methods, for the normal nilpotent orbits of the Lie algebras of any Classical or Exceptional group. We analyse lower dimensioned Exceptional group nilpotent orbits in terms of Hilbert series and the Highest Weight Generating functions for their decompositions into characters of irreducible representations and/or Hall Littlewood polynomials. We investigate the relationships between the moduli spaces describing different nilpotent orbits and propose candidates for the constructions of some non-normal nilpotent orbits of Exceptional algebras.

We revisit the phase diagram of strong-interaction matter for the two-flavor quark-meson model using the Functional Renormalization Group. In contrast to standard mean-field calculations, an unusual phase structure is encountered at low temperatures and large quark chemical potentials. In particular, we identify a regime where the pressure decreases with increasing temperature and discuss possible reasons for this unphysical behavior.

We put forward new explicit reaLISAtions of dS/CFT that relate ${\cal N}=2$ supersymmetric Euclidean vector models with reversed spin-statistics in three dimensions to specific supersymmetric Vasiliev theories in four-dimensional de Sitter space. The partition function of the free supersymmetric vector model deformed by a range of low spin deformations that preserve supersymmetry appears to specify a well-defined wave function with asymptotic de Sitter boundary conditions in the bulk. In particular we find the wave function is globally peaked at undeformed de Sitter space, with a low amplitude for strong deformations. This suggests that supersymmetric de Sitter space is stable in higher-spin gravity and in particular free from ghosts. We speculate this is a limiting case of the de Sitter realizations in exotic string theories.

We study generaLISAtions of ghost-free bimetric theory which involve more than two spin-2 fields. The consistent interactions can enter in the form of two different couplings and in the majority of this work we concentrate on the simpler one. The corresponding action involves one metric coupled to $N$ tensor fields which do not interact with each other. We derive maximally symmetric solutions to the multimetric equations of motion and identify the mass eigenstates in the linearised theory around these backgrounds. Our results are then applied to the problem of singling out multimetric models which possess certain additional structures. In particular, we look for a relation between scale invariant background solutions, the perturbative emergence of Weyl invariance and the presence of partially massless spin-2 fields in the linear theory. Our findings generalise known results in bimetric theory and allow us to point out similarities and differences between the bi- and multimetric models.

Phase Transitions of a (Super) Quantum Mechanical Matrix Model with a Chemical Potential. (arXiv:1707.02898v2 [hep-th] UPDATED)

In this paper, we study the finite-temperature matrix quantum mechanics with chemical potential term linear in the single trace of U(N) matrices, via Monte Carlo simulation. In the bosonic case, we exhibit the existence of the Gross-Witten-Wadia (GWW) type third-order phase transition. We also extend our studies to the model with the fermionic degrees of freedom employing the non-lattice simulation via Fourier expansion, and explore the possibilities that there is a phase transition between the gapped and ungapped phase both in the absence and presence of the chemical potential term. We make a comparison of the phase diagram between the bosonic and fermionic cases.

We revisit the the large spin asymptotics of 15j symbols in terms of cosines of the 4d Euclidean Regge action, as derived by Barrett and collaborators using a saddle point approximation. We bring it closer to the perspective of area-angle Regge calculus and twisted geometries, and compute explicitly the Hessian and phase offsets. We then extend it to more general SU(2) graph invariants, showing that saddle points still exist and have a similar structure. For graphs dual to 4d polytopes we find again two distinct saddle points leading to a cosine asymptotic formula, however a conformal shape-mismatch is allowed by these configurations, and the asymptotic action is thus a generaLISAtion of the Regge action. The allowed mismatch correspond to angle-matched twisted geometries, 3d polyhedral tessellations with adjacent faces matching areas and 2d angles, but not their diagonals. We study these geometries, identify the relevant subsets corresponding to 3d Regge data and flat polytope data, and discuss the corresponding Regge actions emerging in the asymptotics. Finally, we also provide the first numerical confirmation of the large spin asymptotics of the 15j symbol. We show that the agreement is accurate to the per cent level already at spins of order 10, and the next-to-leading order oscillates with the same frequency and same global phase.