The Latest ArXiv Papers (Except Weekends) of 2018-03-18

We demonstrate that the relic gravitational wave background from a multi-step phase transition may deviate from the simple sum of the single spectra, for phase transitions with similar nucleation temperatures $T_N$. We demonstrate that the temperature range $\Delta T$ between the volume fractions $f(T)=0.1$ and $f(T)=0.9$ occupied by the vacuum bubbles can span $\sim 20$ GeV. This allows for a situation in which phase transitions overlap, such that the later bubbles may nucleate both in high temperature and intermediate temperature phases. Such scenarios may lead to more exotic gravitational wave spectra, which cannot be fitted that of a consecutive PTs. We demonstrate this explicitly in the singlet extension of the Standard Model. Finally, we comment on potential additional effects due to the more exotic dynamics of overlapping phase transitions.

The evolution of galaxies is governed by equations with chaotic solutions: gravity and compressible hydrodynamics. While this micro-scale chaos and stochasticity has been well studied, it is poorly understood how it couples to macro-scale properties examined in simulations of galaxy formation. In this paper, we show how perturbations introduced by floating-point roundoff, random number generators, and seemingly trivial differences in algorithmic behaviour can produce non-trivial differences in star formation histories, circumgalactic medium (CGM) properties, and the distribution of stellar mass. We examine the importance of stochasticity due to discreteness noise, variations in merger timings and how self-regulation moderates the effects of this stochasticity. We show that chaotic variations in stellar mass can grow until halted by feedback-driven self-regulation or gas exhaustion. We also find that galaxy mergers are critical points from which large (as much as a factor of 2) variations in quantities such as the galaxy stellar mass can grow. These variations can grow and persist for more than a Gyr before regressing towards the mean. These results show that detailed comparisons of simulations require serious consideration of the magnitude of effects compared to run-to-run chaotic variation, and may significantly complicate interpreting the impact of different physical models. Understanding the results of simulations requires us to understand that the process of simulation is not a mapping of an infinitesimal point in configuration space to another, final infinitesimal point. Instead, simulations map a point in a space of possible initial conditions points to a volume of possible final states.

It is a very well established matter nowadays that many modified gravity models can offer a sound alternative to General Relativity for the description of the accelerated expansion of the universe. But it is also equally well known that no clear and sharp discrimination between any alternative theory and the classical one has been found so far. In this work, we attempt at formulating a different approach starting from the general class of $f(R)$ theories as test probes: we try to reformulate $f(R)$ Lagrangian terms as explicit functions of the redshift, i.e., as $f(z)$. In this context, the $f(R)$ setting to the consensus cosmological model, the $\Lambda$CDM model, can be written as a polynomial including just a constant and a third-order term. Starting from this result, we propose various different polynomial parameterizations $f(z)$, including new terms which would allow for deviations from $\Lambda$CDM, and we thoroughly compare them with observational data. While on the one hand we have found no statistically preference for our proposals (even if some of them are as good as $\Lambda$CDM by using Bayesian Evidence comparison), we think that our novel approach could provide a different perspective for the development of new and observationally reliable alternative models of gravity.

We report a new fast radio burst (FRB) discovered in real-time as part of the UTMOST project at the Molonglo Observatory Synthesis Radio Telescope (MOST). FRB170827 is the first detected with our low-latency ($< 24$ s), machine-learning-based FRB detection system. The FRB discovery was accompanied by the capture of voltage data at the native time and frequency resolution of the observing system, enabling coherent dedispersion and detailed off-line analysis, which have unveiled fine temporal and frequency structure. The dispersion measure (DM) of 176.80 $\pm$ 0.04 pc cm$^{-3}$, is the lowest of the FRB population. The Milky Way contribution along the line of sight is $\sim$ 40 pc cm$^{-3}$, leaving an excess DM of $\sim$ 145 pc cm$^{-3}$. The FRB has a fluence $>$ 20 $\pm$ 7 Jy ms. The burst's width at 10$\%$ maximum is 400 $\mu$s, making it one of the narrowest FRBs. The burst shows three temporal components, the narrowest of which is $\sim$ 30 $\mu$s, and a scattering of $4.1 \pm 2.7$ $\mu$s. The FRB shows spectral modulations on frequency scales of 1.5 MHz and 0.1 MHz. Both are prominent in the dynamic spectrum, which shows a bright region of emission between 841 and 843 MHz, and weaker, patchy emission across the entire band. We show the fine spectral structure could arise in the FRB host galaxy, or its immediate vicinity.

Having an accurate calibration method is crucial for any scientific research done by a radio telescope. The next generation radio telescopes such as the Square Kilometre Array (SKA) will have a large number of receivers which will produce exabytes of data per day. In this paper we propose new direction-dependent and independent calibration algorithms that, while requiring much less storage during calibration, converge very fast. The calibration problem can be formulated as a non-linear least square optimization problem. We show that combining a block-LDU decomposition with Gauss-Newton iterations produces systems of equations with convergent matrices. This allows significant reduction in complexity per iteration and very fast converging algorithms. We also discuss extensions to direction-dependent calibration. The proposed algorithms are evaluated using simulations.

We currently know about 30 magnetars: seemingly isolated neutron stars whose properties can be (in part) comprehended only acknowledging that they are endowed with magnetic fields of complex morphology and exceptional intensity-at least in some components of the field structure. Although magnetars represent only a small percentage of the known isolated neutron stars, there are almost certainly many more of them, since most magnetars were discovered in transitory phases called outbursts, during which they are particularly noticeable. In outburst, in fact, a magnetar can be brighter in X-rays by orders of magnitude and usually emit powerful bursts of hard-X/soft-gamma-ray photons that can be detected almost everywhere in the Galaxy with all-sky monitors such as those on board the Fermi satellite or the Neil Gehrels Swift Observatory. Magnetars command great attention because the large progress that has been made in their understanding is proving fundamental to fathom the whole population of isolated neutron stars, and because, due to their extreme properties, they are relevant for a vast range of different astrophysical topics, from the study of gamma-ray bursts and superluminous supernovae, to ultraluminous X-ray sources, fast radio bursts, and even to sources of gravitational waves. Several excellent reviews with different focuses were published on magnetars in the last few years: among others, Israel and Dall'Osso (2011); Rea and Esposito (2011); Turolla and Esposito (2013); Mereghetti et al. (2015); Turolla et al. (2015); Kaspi and Beloborodov (2017). Here, we quickly recall the history of these sources and travel through the main observational facts, trying to touch some recent and sometimes little-discussed ramifications of magnetars.

Gamma-ray burst (GRB) jets are narrow, and thus typically point away from us. They are initially ultra-relativistic, causing their prompt $\gamma$-ray and early afterglow emission to be beamed away from us. However, as the jet gradually decelerates its beaming cone widens and eventually reaches our line of sight and the afterglow emission may be detected. Such orphan afterglows were not clearly detected so far. Nevertheless, they should be detected in upcoming optical or radio surveys, and it would be challenging to clearly distinguish between them and other types of transients. Therefore, we perform detailed, realistic calculations of the expected afterglow emission from GRB jets viewed at different angles from the jet's symmetry axis. The dynamics are calculated using 2D relativistic hydrodynamics simulations of jets propagating into different power-law external density profiles, $\rho_{\rm ext}\propto{}R^{-k}$ for $k=0,\,1,\,1.5,\,2$, ranging from a uniform ISM-like medium ($k=0$) to a stratified steady stellar-wind like profile ($k=2$). We calculate radio, optical and X-ray lightcurves, and the evolution of the radio afterglow image size, shape and flux centroid. This may help identify misaligned relativistic jets, whether initially ultra-relativistic and producing a GRB for observers within their beam, or (possibly intrinsically more common) moderately relativistic, in either (i) nearby supernovae Ib/c (some of which are associated with long duration GRBs), or (ii) in binary neutron star mergers, which may produce short duration GRBs, and may also be detected in gravitational waves (e.g. GW$\,$170827/GRB$\,$170817A with a weak prompt $\gamma$-ray emission may harbor an off-axis jet).

The X-ray to radio afterglow emission of GRB 170817A/GW 170817 so far scales as $F_\nu\propto\nu^{-0.6}t^{0.8}$ with observed frequency and time, consistent with a single power-law segment of the synchrotron spectrum from the external shock going into the ambient medium. This requires the effective isotropic equivalent afterglow shock energy in the visible region to increase as $\sim t^{1.7}$. The two main channels for such an energy increase are (i) \emph{radial}: more energy carried by slower material (in the visible region) gradually catches up with the afterglow shock and energizes it, and (ii) \emph{angular}: more energy in relativistic outflow moving at different angles to our line of sight, whose radiation is initially beamed away from us but its beaming cone gradually reaches our line of sight as it decelerates. One cannot distinguish between these explanations (or combinations of them) using only the X-ray to radio $F_\nu(t)$. Here we demonstrate that the most promising way to break this degeneracy is through afterglow imaging and polarization, by calculating the predicted evolution of the afterglow image (its size, shape and flux centroid) and linear polarization $\Pi(t)$ for different angular and/or radial outflow structures that fit $F_\nu(t)$. We consider two angular profiles -- a Gaussian and a narrow core with power-law wings in energy per solid angle, as well as a (cocoon motivated) (quasi-) spherical flow with radial velocity profile. For a jet viewed off-axis (and a magnetic field produced in the afterglow shock) $\Pi(t)$ peaks when the jet's core becomes visible, at $\approx2t_p$ where the lightcurve peaks at $t_p$, and the image can be elongated with aspect ratios$\;\gtrsim2$. A quasi-spherical flow has an almost circular image and a much lower $\Pi(t)$ (peaking at $\approx t_p$) and flux centroid displacement.

We present possible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next decade, with the intention of providing information to the astronomy community to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals, and study the capability of the network to determine the sky location of the source. We report our findings for gravitational-wave transients, with particular focus on gravitational-wave signals from the inspiral of binary neutron star systems, which are the most promising targets for multi-messenger astronomy. The ability to localize the sources of the detected signals depends on the geographical distribution of the detectors and their relative sensitivity, and 90% credible regions can be as large as thousands of square degrees when only two sensitive detectors are operational. Determining the sky position of a significant fraction of detected signals to areas of 5 square degrees to 20 square degrees requires at least three detectors of sensitivity within a factor of ~2 of each other and with a broad frequency bandwidth.

Recently we have witnessed the first multi-messenger detection of colliding neutron stars through gravitational waves (GWs) and Electromagnetic (EM) waves (GW170817), thanks to the joint efforts of LIGO/Virgo and Space/Ground-based telescopes. In this paper, we report on the RATIR followup observation strategies and show the results for the trigger G194575. This trigger is not of astrophysical interest; however, is of great interests to the robust design of a followup engine to explore large sky error regions. We discuss the development of an image-subtraction pipeline for the 6-color, optical/NIR imaging camera RATIR. Considering a two band ($i$ and $r$) campaign in the Fall of 2015, we find that the requirement of simultaneous detection in both bands leads to a factor $\sim$10 reduction in false alarm rate, which can be further reduced using additional bands. We also show that the performance of our proposed algorithm is robust to fluctuating observing conditions, maintaining a low false alarm rate with a modest decrease in system efficiency that can be overcome utilizing repeat visits. Expanding our pipeline to search for either optical or NIR detections (3 or more bands), considering separately the optical $riZ$ and NIR $YJH$ bands, should result in a false alarm rate $\approx 1\%$ and an efficiency $\approx 90\%$. RATIR's simultaneous optical/NIR observations are expected to yield about one candidate transient in the vast 100 $\mathrm{deg^2}$ LIGO error region for prioritized followup with larger aperture telescopes.

The distribution of elements in galaxies provides a wealth of information about their production sites and their subsequent mixing into the interstellar medium. Here we investigate the distribution of elements within stars in the IllustrisTNG simulations. In particular, we analyze the abundance ratios of magnesium and europium in Milky Way-like galaxies from the TNG100 simulation (stellar masses ${\log} (M_\star / {\rm M}_\odot) \sim 9.7 - 11.2$). As abundances of magnesium and europium for individual stars in the Milky Way are observed across a variety of spatial locations and metallicities, comparison with the stellar abundances in our more than $850$ Milky Way-like galaxies provides stringent constraints on our chemical evolutionary methods. To this end we use the magnesium to iron ratio as a proxy for the effects of our SNII and SNIa metal return prescription, and a means to compare our simulated abundances to a wide variety of galactic observations. The europium to iron ratio tracks the rare ejecta from neutron star -- neutron star mergers, the assumed primary site of europium production in our models, which in turn is a sensitive probe of the effects of metal diffusion within the gas in our simulations. We find that europium abundances in Milky Way-like galaxies show no correlation with assembly history, present day galactic properties, and average galactic stellar population age. In general, we reproduce the europium to iron spread at low metallicities observed in the Milky Way, with the level of enhancement being sensitive to gas properties during redshifts $z \approx 2-4$. We show that while the overall normalization of [Eu/Fe] is susceptible to resolution and post-processing assumptions, the relatively large spread of [Eu/Fe] at low [Fe/H] when compared to that at high [Fe/H] is very robust.

Short GRBs: opening angles, local neutron star merger rate and off-axis events for GRB/GW association. (arXiv:1708.07008v3 [astro-ph.HE] UPDATED)

The jet breaks in the afterglow lightcurves of short gamma-ray bursts (SGRBs), rarely detected so far, are crucial for estimating the half-opening angles of the ejecta ($\theta_{\rm j}$) and hence the neutron star merger rate. In this work we report the detection of jet decline behaviors in GRB 150424A and GRB 160821B and find $\theta_{\rm j}\sim 0.1$ rad. Together with five events reported before 2015 and other three "identified" recently (GRB 050709, GRB 060614 and GRB 140903A), we have a sample consisting of nine SGRBs and one long-short GRB with reasonably estimated $\theta_{\rm j}$. In particular, three {\it Swift} bursts in the sample have redshifts $z\leq 0.2$, with which we estimate the local neutron star merger rate density {to be $\sim 1109^{+1432}_{-657}~{\rm Gpc^{-3}~yr^{-1}}$ or $162^{+140}_{-83} {\rm Gpc^{-3}yr^{-1}}$ if the narrowly-beamed GRB 061201 is excluded}. Inspired by the typical $\theta_{\rm j}\sim 0.1$ rad found currently, we further investigate whether the off-beam GRBs (in the uniform jet model) or the off-axis events (in the structured jet model) can significantly enhance the GRB/GW association or not. For the former the enhancement is at most moderate, while for the latter the enhancement can be much greater and a high GRB/GW association probability of $\sim 10\%$ is possible. We also show that the data of GRB 160821B may contain a macronova/kilonova emission component with a temperature of $\sim 3100$ K at $\sim 3.6$ days after the burst and more data are needed to ultimately clarify.

We use 1169 Pan-STARRS supernovae (SNe) and 195 low-$z$ ($z < 0.1$) SNe Ia to measure cosmological parameters. Though most Pan-STARRS SNe lack spectroscopic classifications, in a previous paper (I) we demonstrated that photometrically classified SNe can be used to infer unbiased cosmological parameters by using a Bayesian methodology that marginalizes over core-collapse (CC) SN contamination. Our sample contains nearly twice as many SNe as the largest previous SN Ia compilation. Combining SNe with Cosmic Microwave Background (CMB) constraints from Planck, we measure the dark energy equation of state parameter $w$ to be -0.989$\pm$0.057 (stat$+$sys). If $w$ evolves with redshift as $w(a) = w_0 + w_a(1-a)$, we find $w_0 = -0.912 \pm 0.149$ and $w_a =$ -0.513$\pm$0.826. These results are consistent with cosmological parameters from the Joint Lightcurve Analysis and the Pantheon sample. We try four different photometric classification priors for Pan-STARRS SNe and two alternate ways of modeling CC SN contamination, finding that no variant gives a $w$ differing by more than 2% from the baseline measurement. The systematic uncertainty on $w$ due to marginalizing over CC SN contamination, $\sigma_w^{\textrm{CC}} = 0.012$, is the third-smallest source of systematic uncertainty in this work. We find limited (1.6$\sigma$) evidence for evolution of the SN color-luminosity relation with redshift, a possible systematic that could constitute a significant uncertainty in future high-$z$ analyses. Our data provide one of the best current constraints on $w$, demonstrating that samples with $\sim$5% CC SN contamination can give competitive cosmological constraints when the contaminating distribution is marginalized over in a Bayesian framework.

We present 3D global kinetic pulsar magnetosphere models, where the charged particle trajectories and the corresponding electromagnetic fields are treated self-consistently. For our study, we have developed a cartesian 3D relativistic particle-in-cell code that incorporates the radiation reaction forces. We describe our code and discuss the related technical issues, treatments, and assumptions. Injecting particles up to large distances in the magnetosphere, we apply arbitrarily low to high particle injection rates and get an entire spectrum of solutions from close to the Vacuum-Retarded-Dipole to close to the Force-Free solution, respectively. For high particle injection rates (close to FF solutions) significant accelerating electric field components are confined only near the equatorial current sheet outside the light-cylinder. A judicious interpretation of our models allows the calculation of the particle emission and consequently the derivation of the corresponding realistic high-energy sky-maps and spectra. Using model parameters that cover the entire range of spin-down powers of Fermi young and millisecond pulsars, we compare the corresponding model $\gamma$-ray light-curves, cutoff energies, and total $\gamma$-ray luminosities with those observed by Fermi to discover a dependence of the particle injection-rate, $\mathcal{F}$, on the spin-down power, $\dot{\mathcal{E}}$, indicating an increase of $\mathcal{F}$ with $\dot{\mathcal{E}}$. Our models guided by Fermi observations provide field-structures and particle distributions that are not only consistent with each other but also able to reproduce a broad range of the observed $\gamma$-ray phenomenology of both young and millisecond pulsars.

A proto-neutron star (PNS) is a newly formed compact object in a core collapse supernova. In this Letter, the neutrino emission from the cooling process of a PNS is investigated using two types of nuclear equation of state (EOS). It is found that the neutrino signal is mainly determined by the high-density EOS. The neutrino luminosity and mean energy are higher and the cooling time scale is longer for the softer EOS. Meanwhile, the neutrino mean energy and the cooling time scale are also affected by the low-density EOS because of the difference in the population of heavy nuclei. Heavy nuclei have a large scattering cross section with neutrinos owing to the coherent effects and act as thermal insulation near the surface of a PNS. The neutrino mean energy is higher and the cooling time scale is longer for an EOS with a large symmetry energy at low densities, namely a small density derivative coefficient of the symmetry energy, $L$.

The Stellar Populations of Two Ultra-Diffuse Galaxies from Optical and Near-infrared Photometry. (arXiv:1711.05272v3 [astro-ph.GA] UPDATED)

We present observational constraints on the stellar populations of two ultra-diffuse galaxies (UDGs) using optical through near-infrared (NIR) spectral energy distribution (SED) fitting. Our analysis is enabled by new $Spitzer$-IRAC 3.6 $\mu$m and 4.5 $\mu$m imaging, archival optical imaging, and the prospector fully Bayesian SED fitting framework. Our sample contains one field UDG (DGSAT I), one Virgo cluster UDG (VCC 1287), and one Virgo cluster dwarf elliptical for comparison (VCC 1122). We find that the optical--NIR colors of the three galaxies are significantly different from each other. We infer that VCC 1287 has an old ($\gtrsim7.7$ Gyr) and surprisingly metal-poor ($[Z/Z_{\odot}]\lesssim-1.0$) stellar population, even after marginalizing over uncertainties on diffuse interstellar dust. In contrast, the field UDG DGSAT I shows evidence of being younger than the Virgo UDG, with an extended star formation history and an age posterior extending down to $\sim3$ Gyr. The stellar metallicity of DGSAT I is sub-solar but higher than that of the Virgo UDG, with $[Z/Z_{\odot}]=-0.63^{+0.35}_{-0.62}$; in the case of exactly zero diffuse interstellar dust, DGSAT I may even have solar metallicity. With VCC 1287 and several Coma UDGs, a general picture is emerging where cluster UDGs may be "failed" galaxies, but the field UDG DGSAT I seems more consistent with a stellar feedback-induced expansion scenario. In the future, our approach can be applied to a large and diverse sample of UDGs down to faint surface brightness limits, with the goal of constraining their stellar ages, stellar metallicities, and circumstellar and diffuse interstellar dust content.

The search for extra dimensions is a challenging endeavor to probe physics beyond the Standard Model. The joint detection of gravitational waves (GW) and electromagnetic (EM) signals from the merging of a binary system of compact objects like neutron stars (NS), can help constrain the geometry of extra dimensions beyond our 3+1 spacetime ones. A theoretically well-motivated possibility is that our observable Universe is a 3+1-dimensional hypersurface, or brane, embedded in a higher 4+1-dimensional Anti-de Sitter (AdS$_5$) spacetime, in which gravity is the only force which propagates through the infinite bulk space, while other forces are confined to the brane. In these types of brane-world models, GW and EM signals between two points on the brane would, in general, travel different paths. This would result in a time-lag between the detection of GW and EM signals emitted simultaneously from the same source. We consider the recent near-simultaneous detection of the GW event GW170817 from the LIGO/Virgo collaboration, and its EM counterpart, the short gamma-ray burst GRB170817A detected by the Fermi Gamma-ray Burst Monitor and the INTEGRAL Anti-Coincidence Shield spectrometer. Assuming the standard $\Lambda$-Cold Dark Matter ($\Lambda$CDM) scenario and performing a likelihood analysis which takes into account astrophysical uncertainties associated to the measured time-lag, we set an upper limit of $\ell \lesssim 0.535\,$Mpc at $68\%$ confidence level on the AdS$_5$ radius of curvature $\ell$. Although the bound is not competitive with current Solar System constraints, it is the first time that data from a multi-messenger GW-EM measurement is used to constrain extra-dimensional models. Thus, our work provides a proof-of-principle for the possibility of using multi-messenger astronomy for probing the geometry of our space-time.

It was shown that the dark matter(DM) minihalo around an intermediate mass black hole(IMBH) can be redistributed into a cusp, called the DM minispike. We consider an intermediate-mass-ratio inspiral consisting of an IMBH harbored in a DM minispike with nonannihilating DM particles and a small black hole(BH) orbiting around it. We investigate gravitational waves(GWs) produced by this system and analyze the waveforms with the comprehensive consideration of gravitational pull, dynamical friction and accretion of the minispike and calculate the time difference and phase difference caused by it. We find that for a certain range of frequency, the inspiralling time of the system is dramatically reduced for smaller central IMBH and large density of DM. For the central IMBH with $10^5M_\odot$, the time of merger is ahead, which can be distinguished by LISA, Taiji and Tianqin. We focus on the effect of accretion and compare it with that of gravitational pull and friction. We find that the accretion mass is a small quantity compared to the initial mass of the small BH and the accretion effect is inconspicuous compared with friction. However, the accumulated phase shift caused by accretion is large enough to be detected by LISA, Taiji and Tianqin, which indicate that the accretion effect can not be ignored in the detection of GWs.

We present models of realistic globular clusters with post-Newtonian dynamics for black holes. By modeling the relativistic accelerations and gravitational-wave emission in isolated binaries and during three- and four-body encounters, we find that nearly half of all binary black hole mergers occur inside the cluster, with about 10% of those mergers entering the LIGO/Virgo band with eccentricities greater than 0.1. In-cluster mergers lead to the birth of a second generation of black holes with larger masses and high spins, which, depending on the black hole natal spins, can sometimes be retained in the cluster and merge again. As a result, globular clusters can produce merging binaries with detectable spins regardless of the birth spins of black holes formed from massive stars. These second-generation black holes would also populate any upper mass gap created by pair-instability supernovae.

Deviations from the predictions of general relativity due to energy-momentum squared gravity (EMSG) are expected to become pronounced in the high density cores of neutron stars. We derive the hydrostatic equilibrium equations in EMSG and solve them numerically to obtain the neutron star mass-radius relations for four different realistic equations of state. We use the existing observational measurements of the masses and radii of neutron stars to constrain the free parameter, $\alpha$, that characterises the coupling between matter and spacetime in EMSG. We show that $-10^{-38}\,\mathrm{cm^{3}/erg}<\alpha <+10^{-37}\,\mathrm{cm^{3}/erg}$. Under this constraint, we discuss what contributions EMSG can provide to the physics of neutron stars, in particular, their relevance to the so called \textit{hyperon puzzle} in neutron stars. We also discuss how EMSG alters the dynamics of the early universe from the predictions of the standard cosmological model. We show that EMSG leaves the standard cosmology safely unaltered back to $t\sim 10^{-4}$ seconds at which the energy density of the universe is $\sim 10^{34}\,\mathrm{erg\,cm^{-3}}$.

We study the effects of free space nucleon radius on nuclear matter and neutron stars within the framework of quark mean field model. The nucleon radius is treated self-consistently with this model, where quark confinement is adjusted to fit different values of nucleon radius. Corrections due to center-of-mass motion, quark-pion coupling, and one gluon exchange are included to obtain the nucleon mass in vacuum. The meson coupling constants that describe the behavior of many-body nucleonic system are newly-constructed by reproducing the empirical saturation properties of nuclear matter, including the recent determinations of symmetry energy parameters. Our results show that the nucleon radius in free space have negligible effects on nuclear matter equation of state and neutron star mass-radius relations, which is different from the conclusion drawn in previous studies. We further explore that the sensitivity of star radius on the nucleon radius found in earlier publications is actually from the symmetry energy and its slope.

It is a very well established matter nowadays that many modified gravity models can offer a sound alternative to General Relativity for the description of the accelerated expansion of the universe. But it is also equally well known that no clear and sharp discrimination between any alternative theory and the classical one has been found so far. In this work, we attempt at formulating a different approach starting from the general class of $f(R)$ theories as test probes: we try to reformulate $f(R)$ Lagrangian terms as explicit functions of the redshift, i.e., as $f(z)$. In this context, the $f(R)$ setting to the consensus cosmological model, the $\Lambda$CDM model, can be written as a polynomial including just a constant and a third-order term. Starting from this result, we propose various different polynomial parameterizations $f(z)$, including new terms which would allow for deviations from $\Lambda$CDM, and we thoroughly compare them with observational data. While on the one hand we have found no statistically preference for our proposals (even if some of them are as good as $\Lambda$CDM by using Bayesian Evidence comparison), we think that our novel approach could provide a different perspective for the development of new and observationally reliable alternative models of gravity.

We present possible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next decade, with the intention of providing information to the astronomy community to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals, and study the capability of the network to determine the sky location of the source. We report our findings for gravitational-wave transients, with particular focus on gravitational-wave signals from the inspiral of binary neutron star systems, which are the most promising targets for multi-messenger astronomy. The ability to localize the sources of the detected signals depends on the geographical distribution of the detectors and their relative sensitivity, and 90% credible regions can be as large as thousands of square degrees when only two sensitive detectors are operational. Determining the sky position of a significant fraction of detected signals to areas of 5 square degrees to 20 square degrees requires at least three detectors of sensitivity within a factor of ~2 of each other and with a broad frequency bandwidth.

We present an elastic constitutive model of gravity where we identify physical space with the mid-hypersurface of an elastic hyperplate called the "cosmic fabric" and spacetime with the fabric's world volume. Using a Lagrangian formulation, we show that the fabric's behavior as derived from Hooke's Law is analogous to that of spacetime per the Field Equations of General Relativity. The study is conducted in the limit of small strains, or analogously, in the limit of weak and nearly static gravitational fields. The Fabric's Lagrangian outside of inclusions is shown to have the same form as the Einstein-Hilbert Lagrangian for free space. Properties of the fabric such as strain, stress, vibrations, and elastic moduli are related to properties of gravity and space, such as the gravitational potential, gravitational acceleration, gravitational waves, and the energy density of free space. By introducing a mechanical analogy of General Relativity, we enable the application of Solid Mechanics tools to address problems in Cosmology.

We have recently proposed a new action principle for combining Einstein equations and the Dirac equation for a point mass. We used a length scale $L_{CS}$, dubbed the Compton-Schwarzschild length, to which the Compton wavelength and Schwarzschild radius are small mass and large mass approximations, respectively. Here we write down the field equations which follow from this action. We argue that the large mass limit yields Einstein equations, provided we assume wave function collapse and locaLISAtion for large masses. The small mass limit yields the Dirac equation. We explain why the Kerr-Newman black hole has the same gyromagnetic ratio as the Dirac electron, both being twice the classical value. The small mass limit also provides compelling reasons for introducing torsion, which is sourced by the spin density of the Dirac field. There is thus a symmetry between torsion and gravity: torsion couples to quantum objects through Planck's constant $\hbar$ (but not $G$) and is important in the microscopic limit. Whereas gravity couples to classical matter, as usual, through Newton's gravitational constant $G$ (but not $\hbar$), and is important in the macroscopic limit. We construct the Einstein-Cartan-Dirac equations which include the length $L_{CS}$. We find a potentially significant change in the coupling constant of the torsion driven cubic non-linear self-interaction term in the Dirac-Hehl-Datta equation. We speculate on the possibility that gravity is not a fundamental interaction, but emerges as a consequence of wave function collapse, and that the gravitational constant maybe expressible in terms of Planck's constant and the parameters of dynamical collapse models.

Short GRBs: opening angles, local neutron star merger rate and off-axis events for GRB/GW association. (arXiv:1708.07008v3 [astro-ph.HE] UPDATED)

The jet breaks in the afterglow lightcurves of short gamma-ray bursts (SGRBs), rarely detected so far, are crucial for estimating the half-opening angles of the ejecta ($\theta_{\rm j}$) and hence the neutron star merger rate. In this work we report the detection of jet decline behaviors in GRB 150424A and GRB 160821B and find $\theta_{\rm j}\sim 0.1$ rad. Together with five events reported before 2015 and other three "identified" recently (GRB 050709, GRB 060614 and GRB 140903A), we have a sample consisting of nine SGRBs and one long-short GRB with reasonably estimated $\theta_{\rm j}$. In particular, three {\it Swift} bursts in the sample have redshifts $z\leq 0.2$, with which we estimate the local neutron star merger rate density {to be $\sim 1109^{+1432}_{-657}~{\rm Gpc^{-3}~yr^{-1}}$ or $162^{+140}_{-83} {\rm Gpc^{-3}yr^{-1}}$ if the narrowly-beamed GRB 061201 is excluded}. Inspired by the typical $\theta_{\rm j}\sim 0.1$ rad found currently, we further investigate whether the off-beam GRBs (in the uniform jet model) or the off-axis events (in the structured jet model) can significantly enhance the GRB/GW association or not. For the former the enhancement is at most moderate, while for the latter the enhancement can be much greater and a high GRB/GW association probability of $\sim 10\%$ is possible. We also show that the data of GRB 160821B may contain a macronova/kilonova emission component with a temperature of $\sim 3100$ K at $\sim 3.6$ days after the burst and more data are needed to ultimately clarify.

The search for extra dimensions is a challenging endeavor to probe physics beyond the Standard Model. The joint detection of gravitational waves (GW) and electromagnetic (EM) signals from the merging of a binary system of compact objects like neutron stars (NS), can help constrain the geometry of extra dimensions beyond our 3+1 spacetime ones. A theoretically well-motivated possibility is that our observable Universe is a 3+1-dimensional hypersurface, or brane, embedded in a higher 4+1-dimensional Anti-de Sitter (AdS$_5$) spacetime, in which gravity is the only force which propagates through the infinite bulk space, while other forces are confined to the brane. In these types of brane-world models, GW and EM signals between two points on the brane would, in general, travel different paths. This would result in a time-lag between the detection of GW and EM signals emitted simultaneously from the same source. We consider the recent near-simultaneous detection of the GW event GW170817 from the LIGO/Virgo collaboration, and its EM counterpart, the short gamma-ray burst GRB170817A detected by the Fermi Gamma-ray Burst Monitor and the INTEGRAL Anti-Coincidence Shield spectrometer. Assuming the standard $\Lambda$-Cold Dark Matter ($\Lambda$CDM) scenario and performing a likelihood analysis which takes into account astrophysical uncertainties associated to the measured time-lag, we set an upper limit of $\ell \lesssim 0.535\,$Mpc at $68\%$ confidence level on the AdS$_5$ radius of curvature $\ell$. Although the bound is not competitive with current Solar System constraints, it is the first time that data from a multi-messenger GW-EM measurement is used to constrain extra-dimensional models. Thus, our work provides a proof-of-principle for the possibility of using multi-messenger astronomy for probing the geometry of our space-time.

It was shown that the dark matter(DM) minihalo around an intermediate mass black hole(IMBH) can be redistributed into a cusp, called the DM minispike. We consider an intermediate-mass-ratio inspiral consisting of an IMBH harbored in a DM minispike with nonannihilating DM particles and a small black hole(BH) orbiting around it. We investigate gravitational waves(GWs) produced by this system and analyze the waveforms with the comprehensive consideration of gravitational pull, dynamical friction and accretion of the minispike and calculate the time difference and phase difference caused by it. We find that for a certain range of frequency, the inspiralling time of the system is dramatically reduced for smaller central IMBH and large density of DM. For the central IMBH with $10^5M_\odot$, the time of merger is ahead, which can be distinguished by LISA, Taiji and Tianqin. We focus on the effect of accretion and compare it with that of gravitational pull and friction. We find that the accretion mass is a small quantity compared to the initial mass of the small BH and the accretion effect is inconspicuous compared with friction. However, the accumulated phase shift caused by accretion is large enough to be detected by LISA, Taiji and Tianqin, which indicate that the accretion effect can not be ignored in the detection of GWs.

We present models of realistic globular clusters with post-Newtonian dynamics for black holes. By modeling the relativistic accelerations and gravitational-wave emission in isolated binaries and during three- and four-body encounters, we find that nearly half of all binary black hole mergers occur inside the cluster, with about 10% of those mergers entering the LIGO/Virgo band with eccentricities greater than 0.1. In-cluster mergers lead to the birth of a second generation of black holes with larger masses and high spins, which, depending on the black hole natal spins, can sometimes be retained in the cluster and merge again. As a result, globular clusters can produce merging binaries with detectable spins regardless of the birth spins of black holes formed from massive stars. These second-generation black holes would also populate any upper mass gap created by pair-instability supernovae.

Deviations from the predictions of general relativity due to energy-momentum squared gravity (EMSG) are expected to become pronounced in the high density cores of neutron stars. We derive the hydrostatic equilibrium equations in EMSG and solve them numerically to obtain the neutron star mass-radius relations for four different realistic equations of state. We use the existing observational measurements of the masses and radii of neutron stars to constrain the free parameter, $\alpha$, that characterises the coupling between matter and spacetime in EMSG. We show that $-10^{-38}\,\mathrm{cm^{3}/erg}<\alpha <+10^{-37}\,\mathrm{cm^{3}/erg}$. Under this constraint, we discuss what contributions EMSG can provide to the physics of neutron stars, in particular, their relevance to the so called \textit{hyperon puzzle} in neutron stars. We also discuss how EMSG alters the dynamics of the early universe from the predictions of the standard cosmological model. We show that EMSG leaves the standard cosmology safely unaltered back to $t\sim 10^{-4}$ seconds at which the energy density of the universe is $\sim 10^{34}\,\mathrm{erg\,cm^{-3}}$.

By precisely monitoring the "ticks" of Nature's most precise clocks (millisecond pulsars), scientists are trying to detect the "ripples in spacetime" (gravitational waves) produced by the inspirals of supermassive black holes in the centers of distant merging galaxies. Here we describe a relatively simple demonstration that uses two metronomes and a microphone to illustrate several techniques used by pulsar astronomers to search for gravitational waves. An adapted version of this demonstration could be used as an instructional laboratory investigation at the undergraduate level.

The High-Luminosity LHC will provide the unique opportunity to explore the nature of physics beyond the Standard Model of strong and electroweak interactions. Highly selective first level triggers are essential for the physics programme of the ATLAS experiment at the HL-LHC where the instantaneous luminosity will exceed the LHC Run 1 instantaneous luminosity by almost an order of magnitude. The ATLAS first level muon trigger rate is dominated by low momentum muons, selected due to the moderate momentum resolution of the resistive plate and thin gap trigger chambers. This limitation can be overcome by including the data of the precision muon drift tube (MDT) chambers in the first level trigger decision. This requires the fast continuous transfer of the MDT hits to the off-detector trigger logic and a fast track reconstruction algorithm performed in the trigger logic.

In order to demonstrate the feasibility of reconstructing tracks in MDT chambers within the short available first-level trigger latency of about 3~$\mu$s we implemented a seeded Hough transform on the ARM Cortex A9 microprocessor of a Xilinx Zynq FPGA and studied its performance with test-beam data recorded in CERN's Gamma Irradiation Facility. We could show that by using the ARM processor's Neon Single Instruction Multiple Data Engine to carry out 4 floating point operations in parallel the challenging latency requirement can be matched.

We present a non-iterative solver based on the Schur complement method for sparse linear systems of special form which appear in Quantum Monte-Carlo (QMC) simulations of strongly interacting fermions on the lattice. While the number of floating-point operations for this solver scales as the cube of the number of lattice sites, for practically relevant lattice sizes it is still significantly faster than iterative solvers such as the Conjugate Gradient method in the regime of strong inter-fermion interactions, for example, in the vicinity of quantum phase transitions. The speed-up is even more dramatic for the solution of multiple linear systems with different right-hand sides. We present benchmark results for QMC simulations of the tight-binding models on the hexagonal graphene lattice with on-site (Hubbard) and non-local (Coulomb) interactions, and demonstrate the potential for further speed-up using GPU.

We present a modular analysis program written in Python devoted to the estimation of autocorrelation times for Monte Carlo simulations by means of the $\Gamma$-method algorithm. We give a brief review of this method and describe the main features of the program. The latter is characterized by a user-friendly interface and an open source environment which, along with its modularity, make it a versatile tool. Finally we present a simple application as an operational test for the program.

We demonstrate that the relic gravitational wave background from a multi-step phase transition may deviate from the simple sum of the single spectra, for phase transitions with similar nucleation temperatures $T_N$. We demonstrate that the temperature range $\Delta T$ between the volume fractions $f(T)=0.1$ and $f(T)=0.9$ occupied by the vacuum bubbles can span $\sim 20$ GeV. This allows for a situation in which phase transitions overlap, such that the later bubbles may nucleate both in high temperature and intermediate temperature phases. Such scenarios may lead to more exotic gravitational wave spectra, which cannot be fitted that of a consecutive PTs. We demonstrate this explicitly in the singlet extension of the Standard Model. Finally, we comment on potential additional effects due to the more exotic dynamics of overlapping phase transitions.

A proto-neutron star (PNS) is a newly formed compact object in a core collapse supernova. In this Letter, the neutrino emission from the cooling process of a PNS is investigated using two types of nuclear equation of state (EOS). It is found that the neutrino signal is mainly determined by the high-density EOS. The neutrino luminosity and mean energy are higher and the cooling time scale is longer for the softer EOS. Meanwhile, the neutrino mean energy and the cooling time scale are also affected by the low-density EOS because of the difference in the population of heavy nuclei. Heavy nuclei have a large scattering cross section with neutrinos owing to the coherent effects and act as thermal insulation near the surface of a PNS. The neutrino mean energy is higher and the cooling time scale is longer for an EOS with a large symmetry energy at low densities, namely a small density derivative coefficient of the symmetry energy, $L$.

In this paper a complete generaLISAtion of the Catani-Seymour dipole subtraction method to next-to-leading order electroweak calculations is presented. All singularities due to photon and gluon radiation off both massless and massive partons in the presence of both massless and massive spectators are accounted for. Particular attention is paid to the simultaneous subtraction of singularities of both QCD and electroweak origin which are present in the next-to-leading order corrections to processes with more than one perturbative order contributing at Born level. Similarly, embedding non-dipole-like photon splittings in the dipole subtraction scheme discussed. The implementation of the formulated subtraction scheme in the framework of the Sherpa Monte-Carlo event generator is detailed and numerous internal consistency checks validating the obtained results are presented.

We investigate the perturbative renormaLISAtion of deformed conformal field theories from the Hamiltonian perspective. We discuss the relation with conformal perturbation theory, to which we provide an explicit match up to third order in the coupling, and show how second-order anomalous dimensions in the Wilson-Fisher fixed points are straightforwardly computed in the Hamiltonian framework. The second part of the paper focuses on the cutoff employed in the truncated conformal space approach of Yurov and Zamolodchikov. We discuss the appearance of non-covariant and non-local counterterms to second order in the cutoff, which we concretise in the $\phi^4$ theories, and find a smooth cutoff to deal with subleading oscillations.

Recently, the existence of an Amplituhedron for tree level amplitudes in the bi-adjoint scalar field theory has been proved by Arkhani-Hamed et al. We argue that hyperbolic geometry constitutes a natural framework to address the study of positive geometries in moduli spaces of Riemann surfaces, and thus to try to extend this achievement beyond tree level. In this paper we begin an exploration of these ideas starting from the simplest example of hyperbolic geometry, the hyperbolic plane. The hyperboloid model naturally guides us to re-discover the moduli space Associahedron, and a new version of its kinematical avatar. As a by-product we obtain a solution to the scattering equations which can be interpreted as a special case of the two well known solutions in terms of spinor-helicity formalism. The construction is done in $1+2$ dimensions and this makes harder to understand how to extract the amplitude from the dlog of the space time Associahedron. Nevertheless, we continue the investigation accommodating a loop momentum in the picture. By doing this we are led to another polytope called Halohedron, which was already known to mathematicians. We argue that the Halohedron fulfils many criteria that make it plausible to be understood as a 1-loop Amplituhedron for the cubic theory. Furthermore, the hyperboloid model again allows to understand that a kinematical version of the Halohedron exists and is related to the one living in moduli space by a simple generaLISAtion of the tree level map.

We have recently proposed a new action principle for combining Einstein equations and the Dirac equation for a point mass. We used a length scale $L_{CS}$, dubbed the Compton-Schwarzschild length, to which the Compton wavelength and Schwarzschild radius are small mass and large mass approximations, respectively. Here we write down the field equations which follow from this action. We argue that the large mass limit yields Einstein equations, provided we assume wave function collapse and locaLISAtion for large masses. The small mass limit yields the Dirac equation. We explain why the Kerr-Newman black hole has the same gyromagnetic ratio as the Dirac electron, both being twice the classical value. The small mass limit also provides compelling reasons for introducing torsion, which is sourced by the spin density of the Dirac field. There is thus a symmetry between torsion and gravity: torsion couples to quantum objects through Planck's constant $\hbar$ (but not $G$) and is important in the microscopic limit. Whereas gravity couples to classical matter, as usual, through Newton's gravitational constant $G$ (but not $\hbar$), and is important in the macroscopic limit. We construct the Einstein-Cartan-Dirac equations which include the length $L_{CS}$. We find a potentially significant change in the coupling constant of the torsion driven cubic non-linear self-interaction term in the Dirac-Hehl-Datta equation. We speculate on the possibility that gravity is not a fundamental interaction, but emerges as a consequence of wave function collapse, and that the gravitational constant maybe expressible in terms of Planck's constant and the parameters of dynamical collapse models.

We construct the twist operator for the Snyder space. Our starting point is a non-associative star product related to a Hermitian reaLISAtion of the noncommutative coordinates originally introduced by Snyder. The corresponding coproduct of momenta is non-coassociative. The twist is constructed using a general definition of the star product in terms of a bi-differential operator in the Hopf algebroid approach. The result is given by a closed analytical expression. We prove that this twist reproduces the correct coproducts of the momenta and the Lorentz generators. The twisted Poincar\'{e} symmetry is described by a non-associative Hopf algebra, while the twisted Lorentz symmetry is described by the undeformed Hopf algebra. This new twist might be important in the construction of different types of field theories on Snyder space.