The Latest ArXiv Papers (Except Weekends) of 2018-10-16

PASIPHAE (the Polar-Areas Stellar Imaging in Polarization High-Accuracy Experiment) is an optopolarimetric survey aiming to measure the linear polarization from millions of stars, and use these to create a three-dimensional tomographic map of the magnetic field threading dust clouds within the Milky Way. This map will provide invaluable information for future CMB B-mode experiments searching for inflationary gravitational waves, providing unique information regarding line-of-sight integration effects. Optical polarization observations of a large number of stars at known distances, tracing the same dust that emits polarized microwaves, can map the magnetic field between them. The Gaia mission is measuring distances to a billion stars, providing an opportunity to produce a tomographic map of Galactic magnetic field directions, using optical polarization of starlight. Such a map will not only boost CMB polarization foreground removal, but it will also have a profound impact in a wide range of astrophysical research, including interstellar medium physics, high-energy astrophysics, and evolution of the Galaxy. Taking advantage of the novel technology implemented in our high-accuracy Wide-Area Linear Optical Polarimeters (WALOPs) currently under construction at IUCAA, India, we will engage in a large-scale optopolarimetric program that can meet this challenge: a survey of both northern and southern Galactic polar regions targeted by CMB experiments, covering over 10,000 square degrees, which will measure linear optical polarization of over 360 stars per square degree (over 3.5 million stars, a 1000-fold increase over the state of the art). The survey will be conducted concurrently from the South African Astronomical Observatory in Sutherland, South Africa in the southern hemisphere, and the Skinakas Observatory in Crete, Greece, in the north.

At present, we have almost as many theories to explain fast radio bursts as we have fast radio bursts observed. This landscape will be changing rapidly with CHIME/FRB, recently commissioned in Canada, and HIRAX, under construction in South Africa. This is an opportune time to review existing theories and their observational consequences, allowing us to efficiently curtail viable astrophysical models as more data becomes available. In this article we provide a currently up to date catalogue of the numerous and varied theories proposed for fast radio bursts so far. We also launch an online evolving repository for the use and benefit of the community to dynamically update our theoretical knowledge and discuss constraints and uses of fast radio bursts.

pulsar data analysis pipelines have historically been comprised of bespoke software systems, supporting the off-line analysis of data. However modern data acquisition systems are making off-line analyses impractical. They often output multiple simultaneous high volume data streams, significantly increasing data capture rates. This leads to the accumulation of large data volumes, which are prohibitively expensive to retain. To maintain processing capabilities when off-line analysis becomes infeasible due to cost, requires a shift to on-line data processing. This paper makes four contributions facilitating this shift with respect to the search for radio pulsars: i) it characterises for the modern era, the key components of a pulsar search science (not signal processing) pipeline, ii) it examines the feasibility of implementing on-line pulsar search via existing tools, iii) problems preventing an easy transition to on-line search are identified and explained, and finally iv) it provides the design for a new prototype pipeline capable of overcoming such problems. Realised using Commercial off-the-shelf (COTS) software components, the deployable system is open source, simple, scalable, and cheap to produce. It has the potential to achieve pulsar search design requirements for the Square Kilometre Array (SKA), illustrated via testing under simulated SKA loads.

We analyze the role played by vector-isovector meson interaction in dense matter present in the interior of neutron stars in the light of new measurements made during the double neutron-star merger GW170817. These concern measurements of tidal deformability from gravitational waves and electromagnetic observations. Our study includes three different equations of state that contain different physical assumptions and matter compositions, namely the NL3 family, MBF, and CMF models. Other related quantities/relations analyzed are the neutron matter pressure, symmetry energy slope, stellar masses and radii, and Urca process threshold for stellar cooling.

We present the results of a Keck/NIRSPEC follow-up survey of thirteen late-type T dwarfs (T6-T9), twelve of which have unusually red or blue J-H colors. Previous work suggests that J-H color outliers may represent the high-gravity, low-metallicity (old) and low-gravity, solar-metallicity (young) extremes of the late-type T dwarf population. We use medium-resolution Y- and H-band spectroscopy to probe regions of T dwarf atmospheres that are more sensitive to gravity and metallicity variations than the J band. We find that the spectral morphologies of our sample are largely homogeneous, with peak-normalized, Y- and H-band morphologies consistent with spectral standards. However, three objects stand out as potentially old, with overluminous Y-band spectra compared to their respective spectral standards, and a fourth object stands out as potentially young, with an underluminous Y band. Of these four objects, three have been previously identified as potential metallicity/gravity outliers, including the one object in our sample with a normal J-H color. We fit publicly available atmospheric model grids to our spectra and find that the best-fit physical parameters vary depending on the model used. As we continue to probe the characteristics of the late-T population, differences in synthetic spectra of ~10-20% in the blue wing of the Y band and ~45% at 1.65 microns, for the same physical parameters, must be reconciled. Further development and public availability of nonsolar metallicity models is also recommended. Future progress toward deciphering the impacts of gravity, metallicity, and variability in the late-type T dwarf population will also require high signal-to-noise, multiwavelength and multi-epoch photometry and spectroscopy.

Rotation measure (RM) synthesis is a widely used polarization processing algorithm for reconstructing polarized structures along the line of sight. Performing RM synthesis on large datasets produced by telescopes like LOFAR can be computationally intensive as the computational cost is proportional to the product of the number of input frequency channels, the number of output Faraday depth values to be evaluated and the number of lines of sight present in the data cube. The required computational cost is likely to get worse due to the planned large area sky surveys with telescopes like the Low Frequency Array (LOFAR), the Murchison Widefield Array (MWA), and eventually the Square Kilometre Array (SKA). The massively parallel General Purpose Graphical Processing Units (GPGPUs) can be used to execute some of the computationally intensive astronomical image processing algorithms including RM synthesis. In this paper, we present a GPU-accelerated code, called cuFFS or CUDA-accelerated Fast Faraday Synthesis, to perform Faraday rotation measure synthesis. Compared to a fast single-threaded and vectorized CPU implementation, depending on the structure and format of the data cubes, our code achieves an increase in speed of up to two orders of magnitude. During testing, we noticed that the disk I/O when using the Flexible Image Transport System (FITS) data format is a major bottleneck and to reduce the time spent on disk I/O, our code supports the faster HDFITS format in addition to the standard FITS format. The code is written in C with GPU-acceleration achieved using Nvidia's CUDA parallel computing platform. The code is available at https://github.com/sarrvesh/cuFFS.

The selection of sites for underground gravitational wave detectors based on spectral and cumulative characterisation of the low frequency seismic noise. The evaluation of the collected long term seismological data in the M\'atra Gravitational and Geophysical Laboratory revealed several drawbacks of the previously established characteristics. Here we demonstrate the problematic aspects of the recent measures and suggest more robust and more reliable methodology. In particular, we show, that the mode of the data is noisy, sensitive to the discretization and intrinsic averaging, and the $rms_{2Hz}$ is burdened by irrelevant information and not adapted to the technological changes. Therefore the use of median of the data instead of the mode and also the modification of the frequency limits of the $rms$ is preferable.

Relativistic Fireball Reprise: Radio Suppression at the Onset of Short Magnetar Bursts. (arXiv:1810.06353v1 [astro-ph.HE])

There is growing evidence that a clear distinction between magnetars and radio pulsars may not exist, implying the population of neutron stars that exhibit both radio pulsations and bursting activities could be potentially large. In this situation, new insights into the burst mechanism could be gained by combining the temporal behavior of radio pulsations. We present a general model for radio suppression by relativistic $e^{\pm}$ plasma outflows at the onset of magnetar flares. A sudden ejection of magnetic energy into the magnetosphere would generate a fireball plasma, which is promptly driven to expand at relativistic speed. This would make the plasma cutoff frequency significantly higher than the rest frame radio frequency, resulting in the suppression of radio waves. We analytically show that any GHz radio emission arising from the magnetosphere is suppressed for $\sim100\ {\rm s}$, depending on the total fireball energy. On the other hand, thermal radiation is expected from the hot spot(s) on the stellar surface created by an inflow of dense plasma, which could be the origin of short bursts. Since our hypothesis predicts radio suppression in coincidence with short bursts, this could be an indirect method to constrain the occurrence rate of short bursts at the faint end that remain undetected by X-ray detectors. Furthermore, ultra-fast gamma-ray flashes from the fireball photosphere is also expected as a smoking gun, although the onboard detection is challenging due to its extremely short duration $\sim\mu$s. Finally, our model is applied to the radio pulsar with magnetar-like activities, PSR J1119-6127 in light of recent observations. Implications for fast radio bursts and the possibility of plasma lensing are also discussed.

We present a method for calibrating the flux density scale for images generated by the Amsterdam ASTRON Radio Transient Facility And Analysis Centre (AARTFAAC). AARTFAAC produces a stream of all-sky images at a rate of one second in order to survey the Northern Hemisphere for short duration, low frequency transients, such as the prompt EM counterpart to gravitational wave events, magnetar flares, blazars, and other as of yet unobserved phenomena. Therefore, an independent flux density scaling solution per image is calculated via bootstrapping, comparing the measured apparent brightness of sources in the field to a reference catalogue. However, the lack of accurate flux density measurements of bright sources below 74 MHz necessitated the creation of the AARTFAAC source catalogue, at 60 MHz, which contains 167 sources across the Northern Hemisphere. Using this as a reference results in a sufficiently high number of detected sources in each image to calculate a stable and accurate flux scale per one second snapshot, in real-time.

The James Webb Space Telescope (JWST) Optical Simulation Testbed (JOST) is a hardware simulator for wavefront sensing and control designed to produce JWST-like images. A model of the JWST three mirror anas- tigmat is realized with three lenses in the form of a Cooke triplet, which provides JWST-like optical quality over a field equivalent to a NIRCam module. An Iris AO hexagonally segmented mirror stands in for the JWST primary. This setup successfully produces images extremely similar to expected JWST in-flight point spread functions (PSFs), and NIRCam images from cryotesting, in terms of the PSF morphology and sampling relative to the diffraction limit. The segmentation of the primary mirror into subapertures introduces complexity into wavefront sensing and control (WFS&C) of large space based telescopes like JWST. JOST provides a platform for independent analysis of WFS&C scenarios for both commissioning and maintenance activities on such ob- servatories. We present an update of the current status of the testbed including both single field and wide-field alignment results. We assess the optical quality of JOST over a wide field of view to inform the future imple- mentation of different wavefront sensing algorithms including the currently implemented Linearized Algorithm for Phase Diversity (LAPD). JOST complements other work at the Makidon Laboratory at the Space Telescope Science Institute, including the High-contrast imager for Complex Aperture Telescopes (HiCAT) testbed, that investigates coronagraphy for segmented aperture telescopes. Beyond JWST we intend to use JOST for WFS&C studies for future large segmented space telescopes such as LUVOIR.

The Kepler Space Telescope has discovered thousands of planets via the transit method. The transit timing variations of these planets allows us not only to infer the existence of other planets, transiting or not, but to characterize a number of parameters of the system. Using the transit timing variations of the planets Kepler-159b and 159c, the transit simulator TTVFast, and the Bayesian Inference tool MultiNest, we predict a new non-transiting planet, Kepler-159d, in a resonant 2:1 orbit with Kepler-159c. This configuration is dynamically stable on at least 10 Myr time scales, though we note that other less stable, higher-order resonances could also produce similar TTVs during the three-year window Kepler was in operation.

AM CVn binaries are hydrogen deficient compact binaries with an orbital period in the 5-65 min range and are predicted to be strong sources of persistent gravitational wave radiation. Using Gaia Data Release 2, we present the parallaxes and proper motions of 41 out of the 56 known systems. Compared to the parallax determined using the HST Fine Guidance Sensor we find that the archetype star, AM CVn, is significantly closer than previously thought. This resolves the high luminosity and mass accretion rate which models had difficulty in explaining. Using Pan-STARRS1 data we determine the absolute magnitude of the AM CVn stars. There is some evidence that donor stars have a higher mass and radius than expected for white dwarfs or that the donors are not white dwarfs. Using the distances to the known AM CVn stars we find strong evidence that a large population of AM CVn stars have still to be discovered. As this value sets the background to the gravitational wave signal of LISA, this is of wide interest. We determine the mass transfer rate for 15 AM CVn stars and find that the majority have a rate significantly greater than expected from standard models. This is further evidence that the donor star has a greater size than expected.

One of largest uncertainties in nuclear physics is the relation between the pressure and density of supranuclear matter: the equation of state. Some of this uncertainty may be removed through future gravitational wave observations of neutron star binaries by extracting the tidal deformabilities (or Love numbers) of neutron stars. Previous studies showed that only a certain combination of the individual deformabilities of each body (chirp tidal deformability) can be measured with second-generation gravitational wave interferometers, such as Adv. LIGO, due to correlations between the individual deformabilities. To overcome this, we search for approximately universal (or equation-of-state independent) relations between two combinations of the individual tidal deformabilities, such that once one of them has been measured, the other can be automatically obtained and the individual ones decoupled through these relations. We find an approximately universal relation between the symmetric and the anti-symmetric combination of the individual tidal deformabilities that is equation-of-state-insensitive to $20\%$ for binaries with masses less than $1.7M_\odot$. We show that these relations can be used to eliminate a combination of the tidal parameters from the list of model parameters, thus breaking degeneracies and improving the accuracy in parameter estimation. A simple study shows that the universal relations can improve the accuracy in the extraction of the symmetric combination of tidal parameters by as much as an order of magnitude, making the overall accuracy in the extraction of this parameter slightly better than that of the chirp tidal deformability. These new universal relations and the improved measurement accuracy on tidal parameters not only are important to astrophysics and nuclear physics, but also impact our ability to probe extreme gravity with gravitational waves and cosmology.

fast radio bursts (FRBs), bright transients with millisecond durations at $\sim$ GHz and typical redshifts probably $>0.8$, are likely to be gravitationally lensed by intervening galaxies. Since the time delay between images of strongly lensed FRB can be measured to extremely high precision because of the large ratio $\sim10^9$ between the typical galaxy-lensing delay time $\sim\mathcal{O}$(10 days) and the width of bursts $\sim\mathcal{O}$(ms), we propose strongly lensed FRBs as precision probes of the universe. We show that, within the flat $\Lambda$CDM model, the Hubble constant $H_0$ can be constrained with a $\sim0.91\%$ uncertainty from 10 such systems probably observed with the Square Kilometer Array (SKA) in $<$ 30 years. More importantly, the cosmic curvature can be model-independently constrained to a precision of $\sim0.076$. This constraint can directly test the validity of the cosmological principle and break the intractable degeneracy between the cosmic curvature and dark energy.

We present a new, fully generative model for constructing astronomical catalogs from optical telescope image sets. Each pixel intensity is treated as a Poisson random variable with a rate parameter that depends on the latent properties of stars and galaxies. These latent properties are themselves random, with prior distributions fitted by empirical Bayes. We compare two procedures for posterior inference. One procedure is based on Markov chain Monte Carlo (MCMC) while the other is based on variational inference (VI). We demonstrate that the MCMC procedure excels at quantifying uncertainty while the VI procedure is 1000x faster. For the error metric we consider, both procedures outperform the current state-of-the-art method for measuring the colors, shapes, and morphologies of stars and galaxies. On a supercomputer, the VI procedure efficiently uses 665,000 CPU cores (1.3 million hardware threads) to construct an astronomical catalog from 50 terabytes of images

The Laser Interferometer Space Antenna (LISA) will open the mHz band of the gravitational wave spectrum for exploration. Sensitivity curves are a useful tool for surveying the types of sources that can be detected by the LISA mission. Here we describe how the sensitivity curve is constructed, and how it can be used to compute the signal-to-noise ratio for a wide range of binary systems. We adopt the 2018 LISA Phase-0 reference design parameters. We consider both sky-averaged sensitivities, and the sensitivity to sources at particular sky locations. The calculations are included in a publicly available {\em Python} notebook.

On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function $p(\rho)$ of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as $R_1=10.8^{+2.0}_{-1.7}$ km for the heavier star and $R_2= 10.7^{+2.1}_{-1.5}$ km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than $1.97 \,M_\odot$ as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain $R_1= 11.9^{+1.4}_{-1.4}$ km and $R_2= 11.9^{+1.4}_{-1.4}$ km at the 90% credible level. Finally, we obtain constraints on $p(\rho)$ at supranuclear densities, with pressure at twice nuclear saturation density measured at $3.5^{+2.7}_{-1.7}\times 10^{34} \,\mathrm{dyn}/\mathrm{cm}^{2}$ at the 90% level.

Motivated by the recent anomalies in $b\rightarrow s$ transitions that emerged at LHCb, we consider a model with an $L_\mu - L_\tau$ gauge symmetry and additional vector-like fermions. We find that by introducing supersymmetry the model can be made consistent with the long-standing deviation in the measured value of the anomalous magnetic moment of the muon, $g-2$, and neutralino dark matter of broad mass ranges and properties. In particular, dark matter candidates include the well-known 1 TeV higgsino, which in the MSSM is typically not compatible with solutions to the $g-2$ puzzle. Moreover, its spin-independent cross section could be at the origin of the recent small excess in XENON-1T data. We apply to the model constraints arising from flavor precision measurements and direct searches at the Large Hadron Collider and show that they do not currently exclude the relevant parameter space regions.

We forecast the prospective of detection for a stochastic gravitational wave background sourced by cosmological first-order phase transitions. We focus on first-order phase transitions with negligible plasma effects, and consider the experimental infrastructures built by the end of the LISA mission. We make manifest the synergy among LISA, pulsar time array experiments, and ground-based interferometers. For phase transitions above the TeV scale or below the electroweak scale, LISA can detect the corresponding gravitational wave signal together with Einstein Telescope, SKA or even aLIGO-aVIRGO-KAGRA. For phase transitions at the electroweak scale, instead, LISA can be the only experiment observing the gravitational wave signal. In case of detection, by using a parameter reconstruction method that we anticipate in this work, we show that LISA on its own has the potential to determine when the phase transition occurs and, consequently, the energy scale above which the standard model of particle physics needs to be modified. The result may likely guide the collider community in the post-LHC era.

Neutrinos play a crucial role in the core-collapse supernova (CCSN) explosion mechanism. The requirement of accurately calculating the transport of neutrinos makes simulations of the CCSN mechanism extremely challenging and computationally expensive. Historically, this stiff challenge has been met by making approximations to the full transport equation. In this work, we compare CCSN simulations in one- and two-dimensions with three approximate neutrino transport schemes, each implemented in the {\tt FLASH} simulation framework. We compare a two moment M1 scheme with an analytic closure (M1), the Isotropic Diffusion Source Approximation (IDSA), and the Advanced Spectral Leakage (ASL) method. We identify and discuss the advantages and disadvantages of each scheme. For each approximate transport scheme, we use identical grid setups, hydrodynamics, and gravity solvers to investigate the transport effects on supernova shock dynamics and neutrino quantities. We find that the transport scheme has a small effect on the evolution of protoneutron star (PNS) radius, PNS mass, and the mass accretion rate. The neutrino luminosities, mean energies, and shock radii have a $\sim$ 10-20\% quantitative difference but the overall qualitative trends are fairly consistent between all three approximations. We find larger differences in the gain region properties, including the gain region mass and the net heating rate in the gain region, as well as the strength of PNS convection in the core. We investigate the progenitor, nuclear equation of state, and stochastic perturbation dependence of our simulations and find similar magnitudes of impact on key quantities. We also compare the computational expense of the various approximations.

The Eisenhart lift of a Paul Trap used to store ions in molecular physics is a linearly polarized periodic gravitational wave. A modified version of Dehmelt's Penning Trap is in turn related to circularly polarized periodic gravitational waves, sought for in inflationary models. Similar equations rule also the Lagrange points in Celestial Mechanics. The explanation is provided by anisotropic oscillators.

Extragalactic proper motions can reveal a variety of cosmological and local phenomena over a range of angular scales. These include observer-induced proper motions, such as the secular aberration drift caused by the solar acceleration about the Galactic Center and a secular extragalactic parallax resulting from our motion with respect to the cosmic microwave background rest frame. Cosmological effects include anisotropic expansion, transverse peculiar velocities induced by large scale structure, and the real-time evolution of the baryon acoustic oscillation. Long-period gravitational waves can deflect light rays, producing an apparent quadrupolar proper motion signal. We review these effects, their imprints on global correlated extragalactic proper motions, their expected amplitudes, the current best measurements (if any), and predictions for Gaia. Finally, we describe a possible long-baseline ngVLA program to measure or constrain these proper motion signals. In most cases, the ngVLA can surpass or complement the expected end-of-mission performance of the Gaia mission.

Since 2013, the highest energy IceCube cascade showers overcame the common muon neutrino tracks. This fast flavor changes, above few tens TeV, has been indebt to the injection of the long searched astrophysical neutrino. However for what concern the recent published ICECUBE 54 neutrino, High Energy Starting Events, (HESE) in 2016, as well as the most recent ones of 82 and 103 IceCube events (2017-2018) and the several dozens of thorough-going muon tracks formed around the IceCube, none of them are pointing or clustering toward any expected x, gamma or radio sources: no one in connection to GRB, no toward active BL Lac, neither to AGN source in Fermi catalog. No clear correlation with nearby mass distribution (Local Group), nor with galactic plane. Withal, there have not been any record of the expected double bang due to the tau neutrino birth and decay among several events above 200 TeV energy (we are disregarding for a moment the most celebrated recent first correlated, but unique, muon track and the recent two tau possible identification); no any self-clustering events at tens TeV energy raised in most recent searches. Furthermore, there is a tension between the internal HESE event spectra power index and the external thorough-going muon tracks one. As we will show at the conclusions a more mundane (but a bit more abundant) prompt charmed atmospheric neutrino component may pollute and rule the data, explaining most of the present enlisted IceCube puzzles. We review the last HESE event data shown in early and in most recent papers (and talks in Neutrino 2018) making the case for the simplest conclusions. We believe that most astrophysical neutrino signals are still hidden below the ashes of these new, anyway discovered, prompt atmospheric noise.

We study cosmology in a class of minimally modified gravity (MMG) with two local gravitational degrees of freedom. We classify modified gravity theories into type-I and type-II: theories of type-I have an Einstein frame and can be recast by change of variables as general relativity (GR) with a non-minimal matter coupling, while theories of type-II have no Einstein frame. Considering a canonical transformation of the lapse, the 3-dimensional induced metric and their conjugate momenta we generate type-I MMG. We then show that phenomenological deviations from GR, such as the speed of gravitational waves $c_T$ and the effective gravitational constant for scalar perturbations $G_{\rm eff}$, are characterized by two functions of an auxiliary variable. We study the phenomenology of several models all having $c_T=1$. We obtain a scenario with $c_T=1$ in which the effective equation-of-state parameter of dark energy is different from $-1$ even though the cosmic acceleration is caused by a bare cosmological constant, and we find that it is possible to reconstruct the theory on choosing a selected time-evolution for the effective dark energy component.

We perform a global fit of the extended scalar singlet model with a fermionic dark matter (DM) candidate. Using the most up-to-date results from the $\mathit{Planck}$ measured DM relic density, direct detection limits from the XENON1T (2018) experiment, electroweak precision observables and Higgs searches at colliders, we constrain the 7-dimensional model parameter space. We also find regions in the model parameter space where a successful electroweak baryogenesis (EWBG) is viable. This allows us to compute the gravitational wave (GW) signals arising from the phase transition, and discuss the potential discovery prospects of the model at current and future GW experiments. Our global fit places a strong upper $\mathit{and}$ lower limit on the second scalar mass, the fermion DM mass and the scalar-fermion DM coupling. In agreement with previous studies, we find that our model can simultaneously yield a strong first-order phase transition and saturate the observed DM abundance. More importantly, the GW spectra of viable points can often be within reach of future GW experiments such as LISA, DECIGO and BBO.

Compact neutron star binary systems are produced from binary massive stars through stellar evolution involving up to two supernova explosions. The final stages in the formation of these systems have not been directly observed. We report the discovery of iPTF 14gqr (SN 2014ft), a Type Ic supernova with a fast evolving light curve indicating an extremely low ejecta mass ($\approx 0.2$ solar masses) and low kinetic energy ($\approx 2 \times 10^{50}$ ergs). Early photometry and spectroscopy reveal evidence of shock cooling of an extended He-rich envelope, likely ejected in an intense pre-explosion mass loss episode of the progenitor. Taken together, we interpret iPTF 14gqr as evidence for ultra-stripped supernovae that form neutron stars in compact binary systems.

gravitational waves are produced by orbiting massive binary objects, such as black holes and neutron stars, and propagate as ripples in the very fabric of spacetime. As the waves carry off orbital energy, the two bodies spiral into each other and eventually merge. They are described by Einstein's equations of General Relativity. For the early phase of the orbit, called the inspiral, Einstein equations can be linearized and solved through analytical approximations, while for the late phase, near the merger, we need to solve the fully nonlinear Einstein's equations on supercomputers. In order to recover the gravitational wave for the entire evolution of the binary, a match is required between the inspiral and the merger waveforms. Our objectives are to establish a streamlined matching method, that will allow an analytical calculation of the complete gravitational waveform, while developing a gravitational wave modeling tutorial for undergraduate physics students. We use post-Newtonian (PN) theory for the inspiral phase, which offers an excellent training ground for students, and rely on Mathematica for our calculations, a tool easily accessible to undergraduates. For the merger phase we bypass Einstein's equations by using a simple analytic toy model named the Implicit Rotating Source (IRS). After building the inspiral and merger waveforms, we construct our matching method and validate it by comparing our results with the waveforms for the first detection, GW150914, available as open-source. Several future projects can be developed based from this project: building complete waveforms for all the detected signals, extending the post-Newtonian model to take into account non-zero eccentricity, employing and testing a more realistic analytic model for the merger, building a separate model for the ringdown, and optimizing the matching technique.

In the spirit of effective field theory, the Standard-Model Extension (SME) provides a comprehensive framework to systematically probe the possibility of Lorentz/CPT violation. In the pure gravity sector, operators with mass dimension larger than 4, while in general being advantageous to short-range experiments, are hard to investigate with systems of astronomical size. However, there is exception if the leading-order effects are CPT-violating and velocity-dependent. Here we study the lowest-order operators in the pure gravity sector that violate the CPT symmetry with carefully chosen relativistic binary pulsar systems. Applying the existing analytical results to the dynamics of a binary orbit, we put constraints on various coefficients for Lorentz/CPT violation with mass dimension 5. These constraints, being derived from the post-Newtonian dynamics for the first time, are complementary to those obtained from the kinematics in the propagation of gravitational waves.

A solution of Einstein's vacuum field equation is derived that describes a general boosted rotating black hole relative to a Lorentz frame at future null infinity. The metric contains five independent parameters, mass $m$, rotation $\omega$, the boost parameter $v/c$ and the boost direction defined by $(n_1,n_2,n_3)$ satisfying $(n_1)^2+(n_1)^2+(n_1)^2=1$ and reduces to the Kerr black hole when the boost parameter is zero and $n_1=1$. It corresponds to the most general configuration that an astrophysical black hole must have. The black hole rotates about the $z$-axis with angular momentum proportional to $m \omega$ and the geometry has just one Killing vector $\partial/\partial{u}$, where $u$ is a retarded time coordinate. The ergosphere is asymmetric, with dominant lobes in the direction opposite to the boost. However the event and Cauchy horizons -- which are the core of the black hole structure -- do not alter, being independent of the boost parameters. The singularity has the topology of a circle on planes orthogonal to the boost direction. An important feature of this general boosted solution is that the rotation axis and the boost axis do not coincide so that the two axes precess one about the other relative either to a local frame of the black hole or to a Lorentz frame at future null infinity. We argue that this black hole geometry is the natural set to describe the remnants of the recently observed gravitational wave events $GW150914$, $GW151226$, $GW170814$ and $GW170817$. We discuss possible astrophysical consequences connected with processes in the asymmetric ergosphere and electromagnetic dynamical effects resulting from the rotating black hole moving at relativistic speeds together with the precession of the rotation axis about the direction of the boost.

It is noticed that the Niederer transformation can be used to find the explicit relation between time-dependent linear oscillators, including the most interesting case when one of them is harmonic. A geometric interpretation of this correspondence is provided by a certain subclass of pp-waves, strictly related to proper conformal transformations. This observation allows us to show that the pulses of plane gravitational wave exhibiting the maximal conformal symmetry are analytically solvable. Particularly interesting is the circularly polarized family for which some aspects (such as the classical cross section, velocity memory effect and impulsive limit) are discussed in more detail. Possible applications to the description of interaction of electromagnetic beams with matter are also indicated.

One of largest uncertainties in nuclear physics is the relation between the pressure and density of supranuclear matter: the equation of state. Some of this uncertainty may be removed through future gravitational wave observations of neutron star binaries by extracting the tidal deformabilities (or Love numbers) of neutron stars. Previous studies showed that only a certain combination of the individual deformabilities of each body (chirp tidal deformability) can be measured with second-generation gravitational wave interferometers, such as Adv. LIGO, due to correlations between the individual deformabilities. To overcome this, we search for approximately universal (or equation-of-state independent) relations between two combinations of the individual tidal deformabilities, such that once one of them has been measured, the other can be automatically obtained and the individual ones decoupled through these relations. We find an approximately universal relation between the symmetric and the anti-symmetric combination of the individual tidal deformabilities that is equation-of-state-insensitive to $20\%$ for binaries with masses less than $1.7M_\odot$. We show that these relations can be used to eliminate a combination of the tidal parameters from the list of model parameters, thus breaking degeneracies and improving the accuracy in parameter estimation. A simple study shows that the universal relations can improve the accuracy in the extraction of the symmetric combination of tidal parameters by as much as an order of magnitude, making the overall accuracy in the extraction of this parameter slightly better than that of the chirp tidal deformability. These new universal relations and the improved measurement accuracy on tidal parameters not only are important to astrophysics and nuclear physics, but also impact our ability to probe extreme gravity with gravitational waves and cosmology.

The Laser Interferometer Space Antenna (LISA) will open the mHz band of the gravitational wave spectrum for exploration. Sensitivity curves are a useful tool for surveying the types of sources that can be detected by the LISA mission. Here we describe how the sensitivity curve is constructed, and how it can be used to compute the signal-to-noise ratio for a wide range of binary systems. We adopt the 2018 LISA Phase-0 reference design parameters. We consider both sky-averaged sensitivities, and the sensitivity to sources at particular sky locations. The calculations are included in a publicly available {\em Python} notebook.

Modern observations based on general relativity indicate that the spatial geometry of the expanding, large-scale Universe is very nearly Euclidean. This basic empirical fact is at the core of the so-called "flatness problem", which is widely perceived to be a major outstanding problem of modern cosmology and as such forms one of the prime motivations behind inflationary models. An inspection of the literature and some further critical reflection however quickly reveals that the typical formulation of this putative problem is fraught with questionable arguments and misconceptions and that it is moreover imperative to distinguish between different varieties of problem. It is shown that the observational fact that the large-scale Universe is so nearly flat is ultimately no more puzzling than similar "anthropic coincidences", such as the specific (orders of magnitude of the) values of the gravitational and electromagnetic coupling constants. In particular, there is no fine-tuning problem in connection to flatness of the kind usually argued for. The arguments regarding flatness and particle horizons typically found in cosmological discourses in fact address a mere $single$ issue underlying the standard FLRW cosmologies, namely the $extreme$ improbability of these models with respect to any "reasonable measure" on the "space of all spacetimes". This issue may be expressed in different ways and a phase space formulation, due to Penrose, is presented here. A horizon problem only arises when additional assumptions - which are usually kept implicit and at any rate seem rather speculative - are made.

On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function $p(\rho)$ of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as $R_1=10.8^{+2.0}_{-1.7}$ km for the heavier star and $R_2= 10.7^{+2.1}_{-1.5}$ km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than $1.97 \,M_\odot$ as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain $R_1= 11.9^{+1.4}_{-1.4}$ km and $R_2= 11.9^{+1.4}_{-1.4}$ km at the 90% credible level. Finally, we obtain constraints on $p(\rho)$ at supranuclear densities, with pressure at twice nuclear saturation density measured at $3.5^{+2.7}_{-1.7}\times 10^{34} \,\mathrm{dyn}/\mathrm{cm}^{2}$ at the 90% level.

Adiabatic regularization is a method to remove UV divergences in quantum fields in curved spacetime. For relic gravitational wave generated during inflation, regularization on all $k$-modes of the power spectrum to 2nd adiabatic order, and of the energy density and pressure to 4th order, respectively, causes low-frequency distortions. To avoid these, we regularize only the short modes inside the horizon during inflation (corresponding to the present frequencies $f \gtrsim 10^{9}$Hz), and keep the long modes intact. Doing this does not violate the energy conservation since the $k$-modes of RGW are independent of each other during inflation. The resulting spectra are UV convergent and simultaneously free of low-frequency distortion, and these properties remain in the present spectra after evolution, in contrast to regularization at the present time which has some distortion or irregularities. The spectra generally exhibit quick oscillations in frequency domain, even if the initial spectra during inflation have no oscillations. This pattern is due to the interference between the positive and negative frequency modes developed during cosmic expansion, and may be probed by future RGW detections.

We study the electroweak phase transition within a 5D warped model including a scalar potential with an exponential behavior, and strong back-reaction over the metric, in the infrared. By means of a novel treatment of the superpotential formalism, we explore parameter regions that were previously inaccessible. We find that for large enough values of the t'Hooft parameter (e.g. $N\simeq 25$) the holographic phase transition occurs, and it can force the Higgs to undergo a first order electroweak phase transition, suitable for electroweak baryogenesis. The model exhibits gravitational waves and colliders signatures. It typically predicts a stochastic gravitational wave background observable both at the Laser Interferometer Space Antenna and at the Einstein Telescope. Moreover the radion tends to be heavy enough such that it evades current constraints, but may show up in future LHC runs.

The Eisenhart lift of a Paul Trap used to store ions in molecular physics is a linearly polarized periodic gravitational wave. A modified version of Dehmelt's Penning Trap is in turn related to circularly polarized periodic gravitational waves, sought for in inflationary models. Similar equations rule also the Lagrange points in Celestial Mechanics. The explanation is provided by anisotropic oscillators.

We study solutions of the Einstein-Skyrme model. Firstly we consider test field Skyrmions on the Kerr background. These configurations -- hereafter dubbed Skerrmions -- can be in equilibrium with a Kerr black hole (BH) by virtue of a synchronisation condition. We consider two sectors for Skerrmions. In the sector with non-zero baryon charge, Skerrmions are akin to the known Skyrme solutions on the Schwarzschild background. These topological' configurations reduce to flat spacetime Skyrmions in a vanishing BH mass limit; moreoever, they never become "small" perturbations on the Kerr background: the non-linearities of the Skyrme model are crucial for all such Skerrmions. In the non-topological sector, on the other hand, Skerrmions have no analogue on the Schwarzschild background. Non-topological Skerrmions carry not baryon charge and bifurcate from a subset of Kerr solutions defining an existence line. Therein the appropriate truncation of the Skyrme model yield a linear scalar field theory containing a complex plus a real field, both massive and decoupled, and the Skerrmions reduce to the known stationary scalar clouds around Kerr BHs. Moreover, non-topological Skerrmions trivialise in the vanishing BH mass limit. We then discuss the backreaction of these Skerrmions, that yield rotating BHs with synchronised Skyrme hair, which continously connect to the Kerr solution (self-gravitating Skyrmions) in the non-topological (topological) sector. In particular, the non-topological hairy BHs provide a non-linear reaLISAtion, within the Skyrme model, of the synchronous stationary scalar clouds around Kerr.

We study cosmology in a class of minimally modified gravity (MMG) with two local gravitational degrees of freedom. We classify modified gravity theories into type-I and type-II: theories of type-I have an Einstein frame and can be recast by change of variables as general relativity (GR) with a non-minimal matter coupling, while theories of type-II have no Einstein frame. Considering a canonical transformation of the lapse, the 3-dimensional induced metric and their conjugate momenta we generate type-I MMG. We then show that phenomenological deviations from GR, such as the speed of gravitational waves $c_T$ and the effective gravitational constant for scalar perturbations $G_{\rm eff}$, are characterized by two functions of an auxiliary variable. We study the phenomenology of several models all having $c_T=1$. We obtain a scenario with $c_T=1$ in which the effective equation-of-state parameter of dark energy is different from $-1$ even though the cosmic acceleration is caused by a bare cosmological constant, and we find that it is possible to reconstruct the theory on choosing a selected time-evolution for the effective dark energy component.

For the past year, the HEP.TrkX project has been investigating machine learning solutions to LHC particle track reconstruction problems. A variety of models were studied that drew inspiration from computer vision applications and operated on an image-like representation of tracking detector data. While these approaches have shown some promise, image-based methods face challenges in scaling up to realistic HL-LHC data due to high dimensionality and sparsity. In contrast, models that can operate on the spacepoint representation of track measurements ("hits") can exploit the structure of the data to solve tasks efficiently. In this paper we will show two sets of new deep learning models for reconstructing tracks using space-point data arranged as sequences or connected graphs. In the first set of models, Recurrent Neural Networks (RNNs) are used to extrapolate, build, and evaluate track candidates akin to Kalman Filter algorithms. Such models can express their own uncertainty when trained with an appropriate likelihood loss function. The second set of models use Graph Neural Networks (GNNs) for the tasks of hit classification and segment classification. These models read a graph of connected hits and compute features on the nodes and edges. They adaptively learn which hit connections are important and which are spurious. The models are scaleable with simple architecture and relatively few parameters. Results for all models will be presented on ACTS generic detector simulated data.

The differential branching fraction of the rare decay $\Lambda^{0}_{b} \rightarrow \Lambda \mu^+\mu^-$ is measured as a function of $q^{2}$, the square of the dimuon invariant mass. The analysis is performed using proton-proton collision data, corresponding to an integrated luminosity of $3.0 \mbox{ fb}^{-1}$, collected by the LHCb experiment. Evidence of signal is observed in the $q^2$ region below the square of the $J/\psi$ mass. Integrating over $15 < q^{2} < 20 \mbox{ GeV}^2/c^4$ the branching fraction is measured as $d\mathcal{B}(\Lambda^{0}_{b} \rightarrow \Lambda \mu^+\mu^-)/dq^2 = (1.18 ^{+ 0.09} _{-0.08} \pm 0.03 \pm 0.27) \times 10^{-7} ( \mbox{GeV}^{2}/c^{4})^{-1}$, where the uncertainties are statistical, systematic and due to the normaLISAtion mode, $\Lambda^{0}_{b} \rightarrow J/\psi \Lambda$, respectively. In the $q^2$ intervals where the signal is observed, angular distributions are studied and the forward-backward asymmetries in the dimuon ($A^{l}_{\rm FB}$) and hadron ($A^{h}_{\rm FB}$) systems are measured for the first time. In the range $15 < q^2 < 20 \mbox{ GeV}^2/c^4$ they are found to be $A^{l}_{\rm FB} = -0.05 \pm 0.09 \mbox{ (stat)} \pm 0.03 \mbox{ (syst)}$ and $A^{h}_{\rm FB} = -0.29 \pm 0.07 \mbox{ (stat)} \pm 0.03 \mbox{ (syst)}$.

Motivated by the recent anomalies in $b\rightarrow s$ transitions that emerged at LHCb, we consider a model with an $L_\mu - L_\tau$ gauge symmetry and additional vector-like fermions. We find that by introducing supersymmetry the model can be made consistent with the long-standing deviation in the measured value of the anomalous magnetic moment of the muon, $g-2$, and neutralino dark matter of broad mass ranges and properties. In particular, dark matter candidates include the well-known 1 TeV higgsino, which in the MSSM is typically not compatible with solutions to the $g-2$ puzzle. Moreover, its spin-independent cross section could be at the origin of the recent small excess in XENON-1T data. We apply to the model constraints arising from flavor precision measurements and direct searches at the Large Hadron Collider and show that they do not currently exclude the relevant parameter space regions.

Since 2013, the highest energy IceCube cascade showers overcame the common muon neutrino tracks. This fast flavor changes, above few tens TeV, has been indebt to the injection of the long searched astrophysical neutrino. However for what concern the recent published ICECUBE 54 neutrino, High Energy Starting Events, (HESE) in 2016, as well as the most recent ones of 82 and 103 IceCube events (2017-2018) and the several dozens of thorough-going muon tracks formed around the IceCube, none of them are pointing or clustering toward any expected x, gamma or radio sources: no one in connection to GRB, no toward active BL Lac, neither to AGN source in Fermi catalog. No clear correlation with nearby mass distribution (Local Group), nor with galactic plane. Withal, there have not been any record of the expected double bang due to the tau neutrino birth and decay among several events above 200 TeV energy (we are disregarding for a moment the most celebrated recent first correlated, but unique, muon track and the recent two tau possible identification); no any self-clustering events at tens TeV energy raised in most recent searches. Furthermore, there is a tension between the internal HESE event spectra power index and the external thorough-going muon tracks one. As we will show at the conclusions a more mundane (but a bit more abundant) prompt charmed atmospheric neutrino component may pollute and rule the data, explaining most of the present enlisted IceCube puzzles. We review the last HESE event data shown in early and in most recent papers (and talks in Neutrino 2018) making the case for the simplest conclusions. We believe that most astrophysical neutrino signals are still hidden below the ashes of these new, anyway discovered, prompt atmospheric noise.

We explore the phase structure of a four dimensional $SO(4)$ invariant lattice Higgs-Yukawa model comprising four reduced staggered fermions interacting with a real scalar field. The fermions belong to the fundamental representation of the symmetry group while the three scalar field components transform in the self-dual representation of $SO(4)$. The model is a generalization of a four fermion system with the same symmetries that has received recent attention because of its unusual phase structure comprising massless and massive symmetric phases separated by a very narrow phase in which a small bilinear condensate breaking $SO(4)$ symmetry is present. The generalization described in this paper simply consists of the addition of a scalar kinetic term. We find a region of the enlarged phase diagram which shows no sign of a fermion condensate or symmetry breaking but in which there is nevertheless evidence of a diverging correlation length. Our results in this region are consistent with the presence of a single continuous phase transition separating the massless and massive symmetric phases observed in the earlier work.

We introduce a novel approach for solving the problem of identifying regions in the framework of Method of Regions by considering singularities and the associated Landau equations given a multi-scale Feynman diagram. These equations are then analyzed by an expansion in a small threshold parameter via the Power Geometry technique. This effectively leads to the analysis of Newton Polytopes which are evaluated using a Mathematica based convex hull program. Furthermore, the elements of the Gr\"{o}bner Basis of the Landau Equations give a family of transformations, which when applied, reveal regions like potential and Glauber. Several one-loop and two-loop examples are studied and benchmarked using our algorithm which we call ASPIRE.

We evaluate the theoretical uncertainties in next-to-leading order plus parton shower predictions for top quark pair production and decay in hadronic collisions. Our work is carried out using the Herwig 7 event generator and presents an in-depth study of variations in matching schemes with two systematically different shower algorithms, the traditional angular-ordered and alternative dipole shower. We also present all of the required extensions of the Herwig dipole shower algorithm to properly take into account quark mass effects, as well as its ability to perform top quark decays. The predictions are compared at parton level as well as to LHC data, including in the boosted regime. We find that the regions where predictions with a non-top-quark-specific tune differ drastically from data are plagued by large uncertainties which are consistent between our two shower and matching algorithms.

The physics of the solar dark matter (DM) that are captured and thermalise through the DM-nucleon interaction has been extensively studied. In this work, we consider the leptophilic DM scenario where the DM particles interact exclusively with the electrons through the axial-vector coupling. We investigate relevant phenomenologies in the Sun, including its capture, evaporation and thermaLISAtion, and we calculate the equilibrium distribution using the Monte Carlo methods, rather than adopting a semi-analytic approximation. Based on the analysis, we then determine the minimum testable mass for which the DM-electron coupling strength can be probed via the neutrino observation. Compared to the case of the DM-nucleon interaction, it turns out that minimum detectable mass of the DM-electron interaction is roughly 1 GeV smaller, and a cross section about two orders of magnitude larger is required for the saturation of the annihilation signal.

We study the electroweak phase transition within a 5D warped model including a scalar potential with an exponential behavior, and strong back-reaction over the metric, in the infrared. By means of a novel treatment of the superpotential formalism, we explore parameter regions that were previously inaccessible. We find that for large enough values of the t'Hooft parameter (e.g. $N\simeq 25$) the holographic phase transition occurs, and it can force the Higgs to undergo a first order electroweak phase transition, suitable for electroweak baryogenesis. The model exhibits gravitational waves and colliders signatures. It typically predicts a stochastic gravitational wave background observable both at the Laser Interferometer Space Antenna and at the Einstein Telescope. Moreover the radion tends to be heavy enough such that it evades current constraints, but may show up in future LHC runs.

Motivated by the recent anomalies in $b\rightarrow s$ transitions that emerged at LHCb, we consider a model with an $L_\mu - L_\tau$ gauge symmetry and additional vector-like fermions. We find that by introducing supersymmetry the model can be made consistent with the long-standing deviation in the measured value of the anomalous magnetic moment of the muon, $g-2$, and neutralino dark matter of broad mass ranges and properties. In particular, dark matter candidates include the well-known 1 TeV higgsino, which in the MSSM is typically not compatible with solutions to the $g-2$ puzzle. Moreover, its spin-independent cross section could be at the origin of the recent small excess in XENON-1T data. We apply to the model constraints arising from flavor precision measurements and direct searches at the Large Hadron Collider and show that they do not currently exclude the relevant parameter space regions.

We forecast the prospective of detection for a stochastic gravitational wave background sourced by cosmological first-order phase transitions. We focus on first-order phase transitions with negligible plasma effects, and consider the experimental infrastructures built by the end of the LISA mission. We make manifest the synergy among LISA, pulsar time array experiments, and ground-based interferometers. For phase transitions above the TeV scale or below the electroweak scale, LISA can detect the corresponding gravitational wave signal together with Einstein Telescope, SKA or even aLIGO-aVIRGO-KAGRA. For phase transitions at the electroweak scale, instead, LISA can be the only experiment observing the gravitational wave signal. In case of detection, by using a parameter reconstruction method that we anticipate in this work, we show that LISA on its own has the potential to determine when the phase transition occurs and, consequently, the energy scale above which the standard model of particle physics needs to be modified. The result may likely guide the collider community in the post-LHC era.

Since 2013, the highest energy IceCube cascade showers overcame the common muon neutrino tracks. This fast flavor changes, above few tens TeV, has been indebt to the injection of the long searched astrophysical neutrino. However for what concern the recent published ICECUBE 54 neutrino, High Energy Starting Events, (HESE) in 2016, as well as the most recent ones of 82 and 103 IceCube events (2017-2018) and the several dozens of thorough-going muon tracks formed around the IceCube, none of them are pointing or clustering toward any expected x, gamma or radio sources: no one in connection to GRB, no toward active BL Lac, neither to AGN source in Fermi catalog. No clear correlation with nearby mass distribution (Local Group), nor with galactic plane. Withal, there have not been any record of the expected double bang due to the tau neutrino birth and decay among several events above 200 TeV energy (we are disregarding for a moment the most celebrated recent first correlated, but unique, muon track and the recent two tau possible identification); no any self-clustering events at tens TeV energy raised in most recent searches. Furthermore, there is a tension between the internal HESE event spectra power index and the external thorough-going muon tracks one. As we will show at the conclusions a more mundane (but a bit more abundant) prompt charmed atmospheric neutrino component may pollute and rule the data, explaining most of the present enlisted IceCube puzzles. We review the last HESE event data shown in early and in most recent papers (and talks in Neutrino 2018) making the case for the simplest conclusions. We believe that most astrophysical neutrino signals are still hidden below the ashes of these new, anyway discovered, prompt atmospheric noise.

We perform a global fit of the extended scalar singlet model with a fermionic dark matter (DM) candidate. Using the most up-to-date results from the $\mathit{Planck}$ measured DM relic density, direct detection limits from the XENON1T (2018) experiment, electroweak precision observables and Higgs searches at colliders, we constrain the 7-dimensional model parameter space. We also find regions in the model parameter space where a successful electroweak baryogenesis (EWBG) is viable. This allows us to compute the gravitational wave (GW) signals arising from the phase transition, and discuss the potential discovery prospects of the model at current and future GW experiments. Our global fit places a strong upper $\mathit{and}$ lower limit on the second scalar mass, the fermion DM mass and the scalar-fermion DM coupling. In agreement with previous studies, we find that our model can simultaneously yield a strong first-order phase transition and saturate the observed DM abundance. More importantly, the GW spectra of viable points can often be within reach of future GW experiments such as LISA, DECIGO and BBO.

We study the spectral representation of finite temperature, out of time ordered (OTO) correlators on the multi-time-fold generalised Schwinger-Keldysh contour. We write the contour-ordered correlators as a sum over time-order permutations acting on a funda- mental array of Wightman correlators. We decompose this Wightman array in a basis of column vectors, which provide a natural generaLISAtion of the familiar retarded-advanced basis in the finite temperature Schwinger-Keldysh formalism. The coefficients of this de- composition take the form of generalised spectral functions, which are Fourier transforms of nested and double commutators. Our construction extends a variety of classical results on spectral functions in the SK formalism at finite temperature to the OTO case.

We explore the phase structure of a four dimensional $SO(4)$ invariant lattice Higgs-Yukawa model comprising four reduced staggered fermions interacting with a real scalar field. The fermions belong to the fundamental representation of the symmetry group while the three scalar field components transform in the self-dual representation of $SO(4)$. The model is a generalization of a four fermion system with the same symmetries that has received recent attention because of its unusual phase structure comprising massless and massive symmetric phases separated by a very narrow phase in which a small bilinear condensate breaking $SO(4)$ symmetry is present. The generalization described in this paper simply consists of the addition of a scalar kinetic term. We find a region of the enlarged phase diagram which shows no sign of a fermion condensate or symmetry breaking but in which there is nevertheless evidence of a diverging correlation length. Our results in this region are consistent with the presence of a single continuous phase transition separating the massless and massive symmetric phases observed in the earlier work.

It is noticed that the Niederer transformation can be used to find the explicit relation between time-dependent linear oscillators, including the most interesting case when one of them is harmonic. A geometric interpretation of this correspondence is provided by a certain subclass of pp-waves, strictly related to proper conformal transformations. This observation allows us to show that the pulses of plane gravitational wave exhibiting the maximal conformal symmetry are analytically solvable. Particularly interesting is the circularly polarized family for which some aspects (such as the classical cross section, velocity memory effect and impulsive limit) are discussed in more detail. Possible applications to the description of interaction of electromagnetic beams with matter are also indicated.

Modern observations based on general relativity indicate that the spatial geometry of the expanding, large-scale Universe is very nearly Euclidean. This basic empirical fact is at the core of the so-called "flatness problem", which is widely perceived to be a major outstanding problem of modern cosmology and as such forms one of the prime motivations behind inflationary models. An inspection of the literature and some further critical reflection however quickly reveals that the typical formulation of this putative problem is fraught with questionable arguments and misconceptions and that it is moreover imperative to distinguish between different varieties of problem. It is shown that the observational fact that the large-scale Universe is so nearly flat is ultimately no more puzzling than similar "anthropic coincidences", such as the specific (orders of magnitude of the) values of the gravitational and electromagnetic coupling constants. In particular, there is no fine-tuning problem in connection to flatness of the kind usually argued for. The arguments regarding flatness and particle horizons typically found in cosmological discourses in fact address a mere $single$ issue underlying the standard FLRW cosmologies, namely the $extreme$ improbability of these models with respect to any "reasonable measure" on the "space of all spacetimes". This issue may be expressed in different ways and a phase space formulation, due to Penrose, is presented here. A horizon problem only arises when additional assumptions - which are usually kept implicit and at any rate seem rather speculative - are made.

We study the electroweak phase transition within a 5D warped model including a scalar potential with an exponential behavior, and strong back-reaction over the metric, in the infrared. By means of a novel treatment of the superpotential formalism, we explore parameter regions that were previously inaccessible. We find that for large enough values of the t'Hooft parameter (e.g. $N\simeq 25$) the holographic phase transition occurs, and it can force the Higgs to undergo a first order electroweak phase transition, suitable for electroweak baryogenesis. The model exhibits gravitational waves and colliders signatures. It typically predicts a stochastic gravitational wave background observable both at the Laser Interferometer Space Antenna and at the Einstein Telescope. Moreover the radion tends to be heavy enough such that it evades current constraints, but may show up in future LHC runs.

The Eisenhart lift of a Paul Trap used to store ions in molecular physics is a linearly polarized periodic gravitational wave. A modified version of Dehmelt's Penning Trap is in turn related to circularly polarized periodic gravitational waves, sought for in inflationary models. Similar equations rule also the Lagrange points in Celestial Mechanics. The explanation is provided by anisotropic oscillators.

We study solutions of the Einstein-Skyrme model. Firstly we consider test field Skyrmions on the Kerr background. These configurations -- hereafter dubbed Skerrmions -- can be in equilibrium with a Kerr black hole (BH) by virtue of a synchronisation condition. We consider two sectors for Skerrmions. In the sector with non-zero baryon charge, Skerrmions are akin to the known Skyrme solutions on the Schwarzschild background. These topological' configurations reduce to flat spacetime Skyrmions in a vanishing BH mass limit; moreoever, they never become "small" perturbations on the Kerr background: the non-linearities of the Skyrme model are crucial for all such Skerrmions. In the non-topological sector, on the other hand, Skerrmions have no analogue on the Schwarzschild background. Non-topological Skerrmions carry not baryon charge and bifurcate from a subset of Kerr solutions defining an existence line. Therein the appropriate truncation of the Skyrme model yield a linear scalar field theory containing a complex plus a real field, both massive and decoupled, and the Skerrmions reduce to the known stationary scalar clouds around Kerr BHs. Moreover, non-topological Skerrmions trivialise in the vanishing BH mass limit. We then discuss the backreaction of these Skerrmions, that yield rotating BHs with synchronised Skyrme hair, which continously connect to the Kerr solution (self-gravitating Skyrmions) in the non-topological (topological) sector. In particular, the non-topological hairy BHs provide a non-linear reaLISAtion, within the Skyrme model, of the synchronous stationary scalar clouds around Kerr.

We study cosmology in a class of minimally modified gravity (MMG) with two local gravitational degrees of freedom. We classify modified gravity theories into type-I and type-II: theories of type-I have an Einstein frame and can be recast by change of variables as general relativity (GR) with a non-minimal matter coupling, while theories of type-II have no Einstein frame. Considering a canonical transformation of the lapse, the 3-dimensional induced metric and their conjugate momenta we generate type-I MMG. We then show that phenomenological deviations from GR, such as the speed of gravitational waves $c_T$ and the effective gravitational constant for scalar perturbations $G_{\rm eff}$, are characterized by two functions of an auxiliary variable. We study the phenomenology of several models all having $c_T=1$. We obtain a scenario with $c_T=1$ in which the effective equation-of-state parameter of dark energy is different from $-1$ even though the cosmic acceleration is caused by a bare cosmological constant, and we find that it is possible to reconstruct the theory on choosing a selected time-evolution for the effective dark energy component.

We study the spectral representation of finite temperature, out of time ordered (OTO) correlators on the multi-time-fold generalised Schwinger-Keldysh contour. We write the contour-ordered correlators as a sum over time-order permutations acting on a funda- mental array of Wightman correlators. We decompose this Wightman array in a basis of column vectors, which provide a natural generaLISAtion of the familiar retarded-advanced basis in the finite temperature Schwinger-Keldysh formalism. The coefficients of this de- composition take the form of generalised spectral functions, which are Fourier transforms of nested and double commutators. Our construction extends a variety of classical results on spectral functions in the SK formalism at finite temperature to the OTO case.