The Latest ArXiv Papers (Except Weekends) of 2017-11-18

We present a morphological catalogue for $\sim$ 670,000 galaxies in the Sloan Digital Sky Survey in two flavours: T-Type, related to the Hubble sequence, and Galaxy Zoo 2 (GZ2 hereafter) classification scheme. By combining accurate existing visual classification catalogues with machine learning, we provide the largest and most accurate morphological catalogue up to date. The classifications are obtained with Deep Learning algorithms using Convolutional Neural Networks (CNNs).

We use two visual classification catalogues, the GZ2 and the Nair et al. 2010, for training CNNs with color images in order to obtain T-Types and a series of GZ2 type questions (disk/features, edge-on galaxies, bar signature, bulge prominence, roundness and mergers). We also provide an additional probability enabling a separation between pure elliptical (E) from S0, where the T-Type model is not so efficient. For the T-Type, our results show smaller offset and scatter than previous models trained with support vector machines. For the GZ2 type questions, our models have large accuracy (> 97\%), precision and recall values (> 90\%) when applied to a test sample with the same characteristics as the one used for training, i.e., with small uncertainties in the GZ2 classification. When we apply our models to a subset of galaxies with larger uncertainties in GZ2 (low agreement between classifiers), the resulting probability distributions show a clear bimodality. This allows us to recover a significant fraction of galaxies with a robust classification, for which the GZ2 was classification was uncertain. The catalogue is publicly released with the paper.

We present new Chandra observations of the outer halo of the giant elliptical galaxy NGC 4472 (M49) in the Virgo Cluster. The data extend to 130 kpc (28'), and have a combined exposure time of 150 ks. After identifying optical counterparts using the Next Generation Virgo Cluster Survey to remove background active galactic nuclei and globular cluster (GC) sources, and correcting for completeness, we find that the number of field low-mass X-ray binaries (LMXBs) per unit stellar V-band light increases significantly with galactocentric radius. Because the flux limit of the complete sample corresponds to the Eddington limit for neutron stars in NGC 4472, many of the ~90 field LMXBs in this sample could host black holes. The excess of field LMXBs at large galactocentric radii may be partially caused by natal kicks on black holes and neutron stars in binary systems in the inner part of the galaxy. Furthermore, since the metallicity in the halo of NGC 4472 strongly decreases towards larger galactocentric radii, the number of field LMXBs per unit stellar mass is anti-correlated with metallicity, opposite to what is observed in GCs. Another way to explain the spatial distribution of field LMXBs is therefore a reversed metallicity effect, although we have not identified a mechanism to explain this in terms of stellar and binary evolution.

It is well known that solar filaments are features in the solar atmosphere which show a hemispheric preference in their chirality. The hemispheric preference is such that the dextral chirality dominates in the northern hemisphere while the sinistral chirality dominates in the southern. Determining the strength and cyclic variation of the degree of this hemispheric preference, however, is challenging and tedious and thus needs to be automated. In this paper, we follow Dr. Pietro Bernasconi's algorithm (Bernasconi et al. 2005) to detect filament chirality with two parallel channels of application. The algorithm is applied to H-alpha images by "Advanced Automated Filament Detection and Characterization Code" (AAFDCC) (Bernasconi et al. 2005) or the algorithm is explained to a human and the human determines the chirality of the solar filament. We have conducted this exercise on data during the month of August from years 2000 to 2016 and we found that 83 % of our visually determined filaments follow the hemispheric chirality preference, while 58 % of automatically determined filaments follow it. Our visually compiled results have over 90 % agreement with those of Pevtsov et al. (2003), yet the visually determined chiralities of filaments disagree with automated determinations significantly. We find that the hemispheric preference remained the same between solar cycles 23 and 24 but the preference is very difficult to determine during the solar minimum of 2008-2010 primarily due to the absence of filaments.

EXO 1745-248 is a transient neutron-star low-mass X-ray binary that resides in the globular cluster Terzan 5. We studied the transient during its quiescent state using 18 $Chandra$ observations of the cluster acquired between 2003 and 2016. We found an extremely variable source, with a luminosity variation in the 0.5-10 keV energy range of $\sim3$ orders of magnitude (between $3\times10^{31}$ erg s$^{-1}$ and $2\times10^{34}$ erg s$^{-1}$) on time scales from years down to only a few days. Using an absorbed power-law model to fit its quiescent spectra, we obtained a typical photon index of $\sim1.4$. This indicates that the source is very hard throughout the full X-ray luminosity range. Only at the highest luminosity the spectral fit requires an additional soft component. These quiescent properties are likely related to strong variability in the low-level accretion rate in the system. However, its behavior is strikingly different from that observed for other neutron star transients that are thought to still accrete in quiescence. We compare our results to these systems. We also discuss similarities and differences between our target and the transitional millisecond pulsar IGR J18245-2452, which also has hard spectra and strong variability during quiescence.

The production site of gamma-rays in blazar jet is a matter of active debate. We present a method to locate gamma-ray emission region in the framework of one-zone emission model. The relation between the magnetic field strengths in the radio cores of jet and the distances of these radio cores from central supermassive black hole (SMBH), which is derived from measurements of core-shift-effect, is extrapolated to estimate the distance of gamma-ray emission region from SMBH, $r$. The input magnetic field strength is evaluated from optical variability by using the criteria that the observed optical variability timescale $t_{\rm var}$ should be longer or equal to the synchrotron cooling timescale of electrons radiating in optical bands. We apply the method for PSK 1510-089 and BL Lacertae, and derive $r<0.46\$pc for PSK 1510-089 with $t_{\rm var}\sim$ a few hours, and $r<0.009\$pc for BL Lacertae with $t_{\rm var}\sim$ a few minutes. The estimates of $r$ and magnetic field strength in gamma-ray emission region are independent of multi-wavelength spectra modeling. This will advance our understandings for emission mechanisms and properties of the relativistic jets.

Observations show that heavy elements such as Sr and Ba are ubiquitous in the early Galaxy with puzzling diverse abundance patterns. We propose a new neutron-capture site in the He shells of early (including primordial) stars of $\sim 20$-$30 M_\odot$ to explain the above ubiquity and diversity. Nucleosynthesis at this site can produce abundance patterns observed in many metal-poor stars, including those attributed to the slow neutron-capture process or a combination of the slow and rapid neutron-capture processes. In contrast to the existing paradigm requiring multiple sites to explain these diverse patterns, they can be largely accounted for by the single mechanism proposed here.

Millisecond pulsars (MSPs) have a great potential to set standards in timekeeping, positioning and metadata communication.

Since the advent of the space-based photometric missions such as CoRoT and NASA's Kepler, asteroseismology has acquired a central role in our understanding about stellar physics. The Kepler spacecraft, especially, is still releasing excellent photometric observations that contain a large amount of information not yet investigated. For exploiting the full potential of these data, sophisticated and robust analysis tools are now essential, so that further constraining of stellar structure and evolutionary models can be obtained. In addition, extracting detailed asteroseismic properties for many stars can yield new insights on their correlations to fundamental stellar properties and dynamics. After a brief introduction to the Bayesian notion of probability, I describe the code DIAMONDS for Bayesian parameter estimation and model comparison by means of the nested sampling Monte Carlo (NSMC) algorithm. NSMC constitutes an efficient and powerful method, in replacement to standard Markov chain Monte Carlo, very suitable for high-dimensional and multimodal problems that are typical of detailed asteroseismic analyses, such as the fitting and mode identification of individual oscillation modes in stars (known as peak-bagging). DIAMONDS is able to provide robust results for statistical inferences involving tens of individual oscillation modes, while at the same time preserving a considerable computational efficiency for identifying the solution. In the tutorial, I will present the fitting of the stellar background signal and the peak-bagging analysis of the oscillation modes in a red-giant star, providing an example to use Bayesian evidence for assessing the peak significance of the fitted oscillation peaks.

Detection of Periodicity Based on Independence Tests - III. PhDCoP: Phase Distance Correlation Periodogram. (arXiv:1711.06075v1 [astro-ph.IM])

I present PhDCoP - a new periodicity metric, based on the Distance Correlation concept of Gabor Szekely. For each trial period PhDCoP calculates the distance correlation between the data samples and their phases. PhDCoP requires adaptation of the Szekely distance correlation to circular variables (phases). The resulting periodicity metric is best suited to sparse datasets, and it performs better than other methods for sawtooth-like periodicities. These include Cepheid and RR-Lyrae light curves, as well as radial velocity curves of eccentric spectroscopic binaries. The performance of PhDCoP in other contexts is almost as good as that of the Generalized Lomb-Scargle periodogram. The concept of phase distance correlation can be adapted also to astrometric data, and it has the potential to be suitable also for large evenly-spaced datasets, after some algorithmic perfection.

We report on time-resolved photometry of the 2015 February-March superoutburst of QZ Virginis. The superoutburst consisted of a separated precursor, main superoutburst, and rebrightening. We detected superhumps with a period of 0.061181(42) d between the precursor and main superoutburst. Based on analyses of period changes and amplitudes of superhumps, the observed superhumps were identified as growing superhumps (stage A superhumps). The duration of stage A superhumps was about 5 d, unusually long for SU UMa-type dwarf novae. Using the obtained stage A superhump period, we estimated the mass ratio of QZ Vir to be 0.108(3). This value suggests that QZ Vir is an SU UMa-type dwarf nova evolving toward the period minimum. Based on the present and the previous observations regarding long-lasting stage A superhumps, a time scale of stage A superhumps is likely to be determined by the mass ratio of the system and the temperature of the accretion disk.

On the basis of the SDSS survey and spectral observations with the 6-m telescope of SAO RAS, we have peformed a detailed study of SDSS J170745+302056. By combination of its characteristics -- exponential surface brightness distribution, central surface brightness of stellar disk $\mu_0(B) = 23.25\,^m/\square"$, blue colors, low metallicity and low star formation rate -- the galaxy is a typical low surface brightness spiral galaxy. Exponential scalelength of the galaxy is $\approx$3 kpc, its optical diameter exceeds 20 kpc. SDSS J170745+302056 is a member of a group of five galaxies and probably it is in interaction with UGC 10716. The existence of a large low surface brightness galaxy in such a dense environment is very unusual.

Multi-decade observing campaigns of the globular clusters 47 Tucanae and M15 have led to an outstanding number of discoveries. Here, we report on the latest results of the long-term observations of the pulsars in these two clusters. For most of the pulsars in 47 Tucanae we have measured, among other things, their higher-order spin period derivatives, which have in turn provided stringent constraints on the physical parameters of the cluster, such as its distance and gravitational potential. For M15, we have studied the relativistic spin precession effect in PSR B2127+11C. We have used full-Stokes observations to model the precession effect, and to constrain the system geometry. We find that the visible beam of the pulsar is swiftly moving away from our line of sight and may very soon become undetectable. On the other hand, we expect to see the opposite emission beam sometime between 2041 and 2053.

gravitational wave interferometers have proved the existence of a new class of binary black holes (BBHs) weighting tens of solar masses and they have provided the first reliable measurement of the rate of coalescing black holes (BHs) in the local universe. On another side, long gamma-ray bursts (GRBs) detected with gamma-ray satellites are believed to be associated with the birth of stellar mass BHs, providing a measure of the rate of these events across the history of the universe, thanks to the measure of their cosmological redshift. These two types of sources, which are subject to different detection biases and involve BHs born in different environments with potentially different characteristics, provide complementary information on the birth rate of stellar BHs. We compare here the birth rates of BHs found in BBH mergers and in long GRBs.

We construct a simple model which makes reasonable assumptions on the history of GRB formation, and which takes into account some major uncertainties, like the beaming angle of GRBs or the delay between the formation of BBHs and their coalescence. We use this model to evaluate the ratio of the number of stellar mass BHs formed in BBH mergers to those formed in GRBs. We find that in our reference model the birth rate of stellar BHs in BBH mergers represents from 5% to 50% of the rate of long GRBs, while comparable birth rates are allowed by models with smaller beaming angles or when a stronger evolution with redshift is considered for BBHs with respect to GRBs. We briefly discuss this result in view of our understanding of the progenitors of GRBs and BBH mergers, and we emphasize that this ratio, which will be better constrained in the coming years, can be directly compared with the prediction of stellar evolution models if a single model is used to produce GRBs and of BBH mergers with the same assumptions.

We present a technique to automatically minimise the re-computation when a data processing program is iteratively changed, or added to, as is often the case in exploratory data analysis in radio astronomy. A typical example is flagging and calibration of demanding or unusual observations where visual inspection suggests improvement to the processing strategy. The technique is based on memoization and referentially transparent tasks. We describe a prototype implementation for the CASA data reduction package. This technique improves the efficiency of data analysis while reducing the possibility for user error and improving the reproducibility of the final result.

We apply four statistical learning methods to a sample of $7941$ galaxies ($z<0.06$) from the Galaxy and Mass Assembly (GAMA) survey to test the feasibility of using automated algorithms to classify galaxies. Using $10$ features measured for each galaxy (sizes, colours, shape parameters \& stellar mass) we apply the techniques of Support Vector Machines (SVM), Classification Trees (CT), Classification Trees with Random Forest (CTRF) and Neural Networks (NN), returning True Prediction Ratios (TPRs) of $75.8\%$, $69.0\%$, $76.2\%$ and $76.0\%$ respectively. Those occasions whereby all four algorithms agree with each other yet disagree with the visual classification (`unanimous disagreement') serves as a potential indicator of human error in classification, occurring in $\sim9\%$ of ellipticals, $\sim9\%$ of Little Blue Spheroids, $\sim14\%$ of early-type spirals, $\sim21\%$ of intermediate-type spirals and $\sim4\%$ of late-type spirals \& irregulars. We observe that the choice of parameters rather than that of algorithms is more crucial in determining classification accuracy. Due to its simplicity in formulation and implementation, we recommend the CTRF algorithm for classifying future galaxy datasets. Adopting the CTRF algorithm, the TPRs of the 5 galaxy types are : E, $70.1\%$; LBS, $75.6\%$; S0-Sa, $63.6\%$; Sab-Scd, $56.4\%$ and Sd-Irr, $88.9\%$. Further, we train a binary classifier using this CTRF algorithm that divides galaxies into spheroid-dominated (E, LBS \& S0-Sa) and disk-dominated (Sab-Scd \& Sd-Irr), achieving an overall accuracy of $89.8\%$. This translates into an accuracy of $84.9\%$ for spheroid-dominated systems and $92.5\%$ for disk-dominated systems.

Primordial Black Hole (PBH) is one of the leading non-particle candidates for dark matter (DM). Although several observations severely constrain the amount of PBHs, it is recently pointed out that there is an uncertainty on the microlensing constraints below $\sim 10^{-10} M_\odot$ which was ignored originally but may weaken the constraints significantly. In this paper, facing this uncertainty, we investigate the possibility that PBHs can make up all DM in a broad mass spectrum. Moreover, we propose a concrete inflation model which can simultaneously produce PBHs for all DM in a broad mass spectrum around $\mathcal O(10^{-13}) M_\odot$ and PBHs for LIGO events in a sharp mass spectrum at $\mathcal O(10) M_\odot$.

We describe the AMReX suite of astrophysics codes and their application to modeling problems in stellar astrophysics. Maestro is tuned to efficiently model subsonic convective flows while Castro models the highly compressible flows associated with stellar explosions. Both are built on the block-structured adaptive mesh refinement library AMReX. Together, these codes enable a thorough investigation of stellar phenomena, including Type Ia supernovae and X-ray bursts. We describe these science applications and the approach we are taking to make these codes performant on current and future many-core and GPU-based architectures.

Extended gamma-ray sources around pulsars constrain the origin of the positron flux at Earth. (arXiv:1711.06223v1 [astro-ph.HE])

The unexpectedly high flux of cosmic ray positrons detected at Earth may originate from nearby astrophysical sources, dark matter, or unknown processes of cosmic-ray secondary production. We report the detection, using the HighAltitude Water Cherenkov Observatory (HAWC), of extended tera-electron volt gamma-ray emission coincident with the locations of two nearby middle-aged pulsars (Geminga and PSR B0656+14). The HAWC observations demonstrate that these pulsars are indeed local sources of accelerated leptons, but the measured tera-electron volt emission profile constrains the diffusion of particles away from these sources to be much slower than previously assumed. We demonstrate that the leptons emitted by these objects are therefore unlikely to be the origin of the excess positrons, which may have a more exotic origin.

We investigate a simple holographic model for cold and dense deconfined QCD matter consisting of three quark flavors. Varying the single free parameter of the model and utilizing a Chiral Effective Theory equation of state (EoS) for nuclear matter, we find four different types of compact stars: traditional neutron stars, strange quark stars, and two novel solutions we refer to as hybrid stars of the second and third kind (HS2 and HS3, respectively). The HS2s are composed of a nuclear matter core and a crust made of stable strange quark matter, while the HS3s have both a quark mantle and a nuclear crust on top of a nuclear matter core. For all types of stars constructed, we determine not only their mass-radius relations, but also tidal deformabilities, Love numbers, as well as moments of inertia and the mass distribution. We find that there exists a range of parameter values in our model, for which the novel hybrid stars have properties in very good agreement with all existing bounds on the macroscopic properties of compact stars. In particular, the tidal deformabilities of these solutions are smaller than those of ordinary neutron stars of the same mass, implying that they provide an excellent fit to the recent gravitational wave data GW170817 of LIGO and Virgo.

In this paper we compare the performance of two likelihood ratio based detection statistics namely maximum likelihood ratio statistic and {\it hybrid} statistic designed for the detection of gravitational waves from compact binary coalescence using multiple interferometric detector networks. We perform simulations with non-spinning double neutron star binary system and neutron star-black hole binary systems with spinning as well as non-spinning black hole component masses. The binary injections are distributed uniformly in volume up to 1 Gpc. We observe that, on average, the maximum likelihood ratio statistic recovers $\sim 34.45\%$, $\sim 49.69\%$, $\sim 61.25\%$ and $\sim 69.67\%$ of injections in 2, 3, 4 and 5 detector networks respectively in the case of neutron star-black hole injections for a fixed false alarm probability of $10^{-7}$ in Gaussian noise. Further, we note that, compared to the maximum likelihood ratio statistic, the {\it hybrid} statistic recovers $\sim 7.45\%$, $\sim 4.57\%$, $\sim 2.56\%$ and $\sim 1.22\%$ more injections in 2, 3, 4 and 5 detector networks respectively for the same false alarm probability in Gaussian noise. On the other hand, among binary neutron star injections, the maximum likelihood ratio statistic recovers $\sim 5.587\%$, $\sim 9.917\%$, $\sim 14.73\%$ and $\sim 19.86\%$ of injections in 2, 3, 4 and 5 detector networks respectively and the {\it hybrid} statistic recovers $\sim 14.63\%$, $\sim 12.91\%$, $\sim 11.49\%$ and $\sim 10.29\%$ more injections compared to maximum likelihood ratio statistic in 2, 3, 4 and 5 detector networks respectively.

This work introduces a new interpretation of the gravitational N-body problem, based on the one-point probability density $\Psi$ of finding a particle at a given location of phase space $(\mathbf{x}, \mathbf{v})$ at time $t$ and the associated expected phase-space density $\bar{f}(\mathbf{x}, \mathbf{v}, t) = M \Psi(\mathbf{x}, \mathbf{v}, t)$ where $M$ is the total mass of the system. At variance with the traditional paradigm, we consider that the problem is inherently stochastic, and $\bar{f}$ corresponds to a weighted average over all possible random reaLISAtions corresponds to a weighted average over all possible random reaLISAtions of the initial probability distribution $\Psi(\mathbf{x}, \mathbf{v}, 0)$. In practice, we run several numerical experiments in one dimension where $\bar{f}(x, v, t)$, and thus $\Psi(x, v, t)$, are estimated from the average of a finite number $S$ of independent simulations with $N$ particles each. The proposed approach is extremely efficient from a computational point of view, with modest CPU and memory requirements, and it provides an alternative to traditional N-body simulations when the goal is to study the average properties of N-body systems, at the cost of abandoning the notion of well-defined trajectories for each individual particle. Our results, fully consistent with those previously reported in the literature for the standard deterministic formulation of the problem, highlight the differences between the evolution of the one-point probability density $\Psi (x, v, t)$ and the predictions of the collisionless Boltzmann (Vlasov-Poisson) equation, as well as the relatively subtle dependence on the actual finite number N of particles in the system.

Parametric resonance is among the most efficient phenomena generating gravitational waves (GWs) in the early Universe. The dynamics of parametric resonance, and hence of the GWs, depend exclusively on the resonance parameter $q$. The latter is determined by the properties of each scenario: the initial amplitude and potential curvature of the oscillating field, and its coupling to other species. Previous works have only studied the GW production for fixed value(s) of $q$. We present an analytical derivation of the GW amplitude dependence on $q$, valid for any scenario, which we confront against numerical results. By running lattice simulations in an expanding grid, we study for a wide range of $q$ values, the production of GWs in post-inflationary preheating scenarios driven by parametric resonance. We present simple fits for the final amplitude and position of the local maxima in the GW spectrum. Our parametrization allows to predict the location and amplitude of the GW background today, for an arbitrary $q$. The GW signal can be rather large, as $h^2\Omega_{\rm GW}(f_p) \lesssim 10^{-11}$, but it is always peaked at high frequencies $f_p \gtrsim 10^{7}$ Hz. We also discuss the case of spectator-field scenarios, where the oscillatory field can be e.g.~a curvaton, or the Standard Model Higgs.

We evolve stellar models to study the rotational profiles of the pre-explosion cores of single massive stars that are progenitors of core collapse supernovae (CCSNe), and find large rotational shear above the iron core that might play an important role in the jet feedback explosion mechanism by amplifying magnetic fields before and after collapse. Initial masses of 15 Mo and 30 Mo and various values of the initial rotation velocity are considered, as well as a reduced mass-loss rate along the evolution and the effect of core-envelope coupling through magnetic fields. We find that the rotation profiles just before core collapse differ between models, but share the following properties. (1) There are narrow zones of very large rotational shear adjacent to convective zones. (2) The rotation rate of the inner core is slower than required to form a Keplerian accretion disk. (3) The outer part of the core and the envelope have non-negligible specific angular momentum compared to the last stable orbit around a black hole (BH). Our results suggest the feasibility of magnetic field amplification which might aid a jet-driven explosion leaving behind a neutron star. Alternatively, if the inner core fails in exploding the star, an accretion disk from the outer parts of the core might form and lead to a jet-driven CCSN which leaves behind a BH.

The increasing complexity of the aperture geometry of the future space- and ground based-telescopes will limit the performance of the next generation of coronagraphic instruments for high contrast imaging of exoplanets. We propose here a new closed-loop optimization technique using two deformable mirrors to correct for the effects of complex apertures on coronagraph performance, alternative to the ACAD technique previously developed by our group. This technique, ACAD-OSM, allows the use of any coronagraphs designed for continuous apertures, with complex, segmented, apertures, maintaining high performance in contrast and throughput. We show the capabilities of this technique on several pupil geometries (segmented LUVOIR type aperture, WFIRST, ELTs) for which we obtained high contrast levels with several deformable mirror setups (size, number of actuators, separation between them), coronagraphs (apodized pupil Lyot and vortex coronagraphs) and spectral bandwidths, which will help us present recommendations for future coronagraphic instruments. We show that this active technique handles, without any revision to the algorithm, changing or unknown optical aberrations or discontinuities in the pupil, including optical design misalignments, missing segments and phase errors.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.

Very recently, the gravitational wave event (GW) 170817 was discovered to be associated with the short gamma-ray burst (GRB) 170817A. Multi-wavelength follow-up observations were carried out, and X-ray, optical and radio counterparts to GW170817 were detected. The observations undoubtedly indicate that GRB170817A originates from a binary neutron star (BNS) merger. However, the GRB falls into the low-luminosity class which could have a higher statistical occurrence rate and detection probability than the normal (high-luminosity) class. This implies a possibility that GRB170817A is intrinsically powerful but we are off-axis and only observe its side emission. In this paper, we provide a timely modeling of the multi-wavelength afterglow emission from this GRB and the associated kilonova signal from the merger ejecta, under the assumption of a structured jet, a two-component jet, and an intrinsically less-energetic quasi-isotropic fireball respectively. Comparing the afterglow properties with the multi-wavelength follow-up observations, we can distinguish between these three models. Furthermore, a few model parameters (e.g., the ejecta mass and velocity) can be constrained.

Super-Massive Black Holes weighing up to $\sim 10^9 \, \mathrm{M_{\odot}}$ are in place by $z \sim 7$, when the age of the Universe is $\lesssim 1 \, \mathrm{Gyr}$. This implies a time crunch for their growth since such high masses cannot be easily reached in standard accretion scenarios. Here, we explore the physical conditions that would lead to optimal growth wherein stable super-Eddington accretion would be permitted. Our analysis suggests that the preponderance of optimal conditions depends on two key parameters: the black hole mass and the host galaxy central gas density. In the high-efficiency region of this parameter space, a continuous stream of gas can accrete onto the black hole from large to small spatial scales, assuming a global isothermal profile for the host galaxy. By using analytical initial mass functions for black hole seeds, we find an enhanced probability of high-efficiency growth for seeds with initial masses $\gtrsim 10^4 \, \mathrm{M_{\odot}}$. Our picture suggests that a large population of high-$z$ lower-mass black holes that formed in the low-efficiency region, with low duty cycles and accretion rates, might remain undetectable as quasars, since we predict their bolometric luminosities to be $\lesssim 10^{41} \, \mathrm{erg \, s^{-1}}$. The presence of these sources might be revealed only via gravitational wave detections of their mergers.

This paper discusses the effects of gravitational waves on high-accuracy astrometric observations such as those delivered by Gaia. Depending on the frequency of gravitational waves, two regimes for the influence of gravitational waves on astrometric data are identified: the regime when the effects of gravitational waves directly influence the derived proper motions of astrometric sources and the regime when those effects mostly appear in the residuals of the standard astrometric solution. The paper is focused on the second regime while the known results for the first regime are briefly summarized.

The deflection of light due to a plane gravitational wave is then discussed. Starting from a model for the deflection we derive the corresponding partial derivatives and summarize some ideas for the search strategy of such signals in high-accuracy astrometric data. In order to reduce the dimensionality of the parameter space the use of vector spherical harmonics is suggested and explained. The explicit formulas for the VSH expansion of the astrometric signal of a plain gravitational wave are derived.

Finally, potential sensitivity of Gaia astrometric data is discussed. Potential astrophysical sources of gravitational waves that can be interesting for astrometric detection are identified.

In this paper, we study the gravitational wave polarization modes for some particular $f(R)$ models using Newman-Penrose formalism. We find two extra scalar modes of gravitational wave (longitudinal and transversal modes) in addition to two tensor modes of general relativity. We conclude that gravitational waves correspond to class $II_{6}$ under the Lorentz-invariant E(2) classification of plane null waves for these $f(R)$ models.

We investigate a simple holographic model for cold and dense deconfined QCD matter consisting of three quark flavors. Varying the single free parameter of the model and utilizing a Chiral Effective Theory equation of state (EoS) for nuclear matter, we find four different types of compact stars: traditional neutron stars, strange quark stars, and two novel solutions we refer to as hybrid stars of the second and third kind (HS2 and HS3, respectively). The HS2s are composed of a nuclear matter core and a crust made of stable strange quark matter, while the HS3s have both a quark mantle and a nuclear crust on top of a nuclear matter core. For all types of stars constructed, we determine not only their mass-radius relations, but also tidal deformabilities, Love numbers, as well as moments of inertia and the mass distribution. We find that there exists a range of parameter values in our model, for which the novel hybrid stars have properties in very good agreement with all existing bounds on the macroscopic properties of compact stars. In particular, the tidal deformabilities of these solutions are smaller than those of ordinary neutron stars of the same mass, implying that they provide an excellent fit to the recent gravitational wave data GW170817 of LIGO and Virgo.

In this paper we compare the performance of two likelihood ratio based detection statistics namely maximum likelihood ratio statistic and {\it hybrid} statistic designed for the detection of gravitational waves from compact binary coalescence using multiple interferometric detector networks. We perform simulations with non-spinning double neutron star binary system and neutron star-black hole binary systems with spinning as well as non-spinning black hole component masses. The binary injections are distributed uniformly in volume up to 1 Gpc. We observe that, on average, the maximum likelihood ratio statistic recovers $\sim 34.45\%$, $\sim 49.69\%$, $\sim 61.25\%$ and $\sim 69.67\%$ of injections in 2, 3, 4 and 5 detector networks respectively in the case of neutron star-black hole injections for a fixed false alarm probability of $10^{-7}$ in Gaussian noise. Further, we note that, compared to the maximum likelihood ratio statistic, the {\it hybrid} statistic recovers $\sim 7.45\%$, $\sim 4.57\%$, $\sim 2.56\%$ and $\sim 1.22\%$ more injections in 2, 3, 4 and 5 detector networks respectively for the same false alarm probability in Gaussian noise. On the other hand, among binary neutron star injections, the maximum likelihood ratio statistic recovers $\sim 5.587\%$, $\sim 9.917\%$, $\sim 14.73\%$ and $\sim 19.86\%$ of injections in 2, 3, 4 and 5 detector networks respectively and the {\it hybrid} statistic recovers $\sim 14.63\%$, $\sim 12.91\%$, $\sim 11.49\%$ and $\sim 10.29\%$ more injections compared to maximum likelihood ratio statistic in 2, 3, 4 and 5 detector networks respectively.

Extreme mass ratio in-spirals (EMRIs) are candidate events for gravitational wave detection in the millihertz range (by detectors like LISA and eLISA). These events involve a stellar-mass black hole, or a similar compact object, descending in the gravitational field of a supermassive black hole, eventually merging with it. Properties of the in-spiralling trajectory away from resonance are well known and have been studied extensively, however little is known about the behaviour of these binary systems at resonance, when the radial and lateral frequencies of the orbit become commensurate. We describe the two existing models, the instantaneous frequency approach used by Gair, Bender, and Yunes, and the standard two timescales approach implemented by Flanagan and Hinderer. In both cases, the exact treatment depends on the modelling of the gravitational self-force, which is currently not available. We extend the results in Gair, Bender and Yunes to higher order in the on-resonance flux modification, and argue that the instantaneous frequency approach is also a valid treatment of the resonance problem. The non-linear differential equations which arise in treating resonances are interesting from a mathematical view point. We present our algorithm for perturbative solutions and the results to third order in the infinitesimal parameter, and discuss the scope of this approach.

The direct detection of gravitational waves (GW) from merging binary black holes and neutron stars mark the beginning of a new era in gravitational physics, and it brings forth new opportunities to test theories of gravity. To this end, it is crucial to search for anomalous deviations from general relativity in a model-independent way, irrespective of gravity theories, GW sources, and background spacetimes. In this paper, we propose a new universal framework for testing gravity with GW, based on the generalized propagation of a GW in an effective field theory that describes modification of gravity at cosmological scales. Then we perform a parameter estimation study, showing how well the future observation of GW can constrain the model parameters in the generalized models of GW propagation.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.

Since the discovery of neutrino oscillations, the experimental progress in the last two decades has been very fast, with the precision measurements of the neutrino squared-mass differences and of the mixing angles, including the last unknown mixing angle $\theta_{13}$.

Today a very large set of oscillation results obtained with a variety of experimental configurations and techniques can be interpreted in the framework of three active massive neutrinos, whose mass and flavour eigenstates are related by a 3 $\times$ 3 unitary mixing matrix, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, parameterized by three mixing angles $\theta_{12}$, $\theta_{23}$, $\theta_{13}$ and a CP-violating phase $\delta_{CP}$. The additional parameters governing neutrino oscillations are the squared-mass differences $\Delta m^2_{ji}=m^2_j-m^2_i$, where $m_i$ is the mass of the $i$th neutrino mass eigenstate. This review covers the rise of the PMNS three-neutrino mixing paradigm and the current status of the experimental determination of its parameters.

The next years will continue to see a rich program of experimental endeavour coming to fruition and addressing the three missing pieces of the puzzle, namely the determination of the octant and precise value of the mixing angle $\theta_{23}$, the unveiling of the neutrino mass ordering (whether $m_1 < m_2 < m_3$ or $m_3 < m_1 < m_2$) and the measurement of the CP-violating phase \dcp.

We review recent attempts at dealing with the sign problem in Monte Carlo calculations by deforming the region of integration in the path integral from real to complex fields. We discuss the theoretical foundations, the algorithmic issues and present some results for low dimensional field theories in both imaginary and real time.

We present results of the eigenvalue spectrum for the staggered Dir\"ac operator obtained using a modified Lanczos algorithm. We identify zero modes and non-zero modes. We derive the chiral Ward identity derived from the conserved $U(1)_A$ symmetry, and check it numerically. This is the first step toward construction of an improved method to identify zero modes reliably with staggered fermions.

An implementation of the Monte Carlo (MC) phase space generators coupled with adaptive MC integration/simulation program FOAM is presented. The first program is a modification of the classic phase space generator GENBOD interfaced with the adaptive sampling integrator/generator FOAM. On top of this tool the algorithm suitable for generation of the phase space for an reaction with two leading particles is presented (double-peripheral process with central production of particles). At the same time it serves as an instructive example of construction of a self-adaptive phase space generator/integrator with a modular structure for specialized particle physics calculations.

We review Zimmermann's forest formula, which solves Bogoliubov's recursive $R$-operation for the subtraction of ultraviolet divergences in perturbative Quantum Field Theory. We further discuss a generaLISAtion of the $R$-operation which subtracts besides ultraviolet also Euclidean infrared divergences. This generaLISAtion, which goes under the name of the $R^*$-operation, can be used efficiently to compute renormaLISAtion constants. We will discuss several results obtained by this method with focus on the QCD beta function at five loops as well as the application to hadronic Higgs boson decay rates at N${}^4$LO. This article summarizes a talk given at the Wolfhart Zimmermann Memorial Symposium.

Primordial Black Hole (PBH) is one of the leading non-particle candidates for dark matter (DM). Although several observations severely constrain the amount of PBHs, it is recently pointed out that there is an uncertainty on the microlensing constraints below $\sim 10^{-10} M_\odot$ which was ignored originally but may weaken the constraints significantly. In this paper, facing this uncertainty, we investigate the possibility that PBHs can make up all DM in a broad mass spectrum. Moreover, we propose a concrete inflation model which can simultaneously produce PBHs for all DM in a broad mass spectrum around $\mathcal O(10^{-13}) M_\odot$ and PBHs for LIGO events in a sharp mass spectrum at $\mathcal O(10) M_\odot$.

We investigate a simple holographic model for cold and dense deconfined QCD matter consisting of three quark flavors. Varying the single free parameter of the model and utilizing a Chiral Effective Theory equation of state (EoS) for nuclear matter, we find four different types of compact stars: traditional neutron stars, strange quark stars, and two novel solutions we refer to as hybrid stars of the second and third kind (HS2 and HS3, respectively). The HS2s are composed of a nuclear matter core and a crust made of stable strange quark matter, while the HS3s have both a quark mantle and a nuclear crust on top of a nuclear matter core. For all types of stars constructed, we determine not only their mass-radius relations, but also tidal deformabilities, Love numbers, as well as moments of inertia and the mass distribution. We find that there exists a range of parameter values in our model, for which the novel hybrid stars have properties in very good agreement with all existing bounds on the macroscopic properties of compact stars. In particular, the tidal deformabilities of these solutions are smaller than those of ordinary neutron stars of the same mass, implying that they provide an excellent fit to the recent gravitational wave data GW170817 of LIGO and Virgo.

Parametric resonance is among the most efficient phenomena generating gravitational waves (GWs) in the early Universe. The dynamics of parametric resonance, and hence of the GWs, depend exclusively on the resonance parameter $q$. The latter is determined by the properties of each scenario: the initial amplitude and potential curvature of the oscillating field, and its coupling to other species. Previous works have only studied the GW production for fixed value(s) of $q$. We present an analytical derivation of the GW amplitude dependence on $q$, valid for any scenario, which we confront against numerical results. By running lattice simulations in an expanding grid, we study for a wide range of $q$ values, the production of GWs in post-inflationary preheating scenarios driven by parametric resonance. We present simple fits for the final amplitude and position of the local maxima in the GW spectrum. Our parametrization allows to predict the location and amplitude of the GW background today, for an arbitrary $q$. The GW signal can be rather large, as $h^2\Omega_{\rm GW}(f_p) \lesssim 10^{-11}$, but it is always peaked at high frequencies $f_p \gtrsim 10^{7}$ Hz. We also discuss the case of spectator-field scenarios, where the oscillatory field can be e.g.~a curvaton, or the Standard Model Higgs.

Since the discovery of neutrino oscillations, the experimental progress in the last two decades has been very fast, with the precision measurements of the neutrino squared-mass differences and of the mixing angles, including the last unknown mixing angle $\theta_{13}$.

Today a very large set of oscillation results obtained with a variety of experimental configurations and techniques can be interpreted in the framework of three active massive neutrinos, whose mass and flavour eigenstates are related by a 3 $\times$ 3 unitary mixing matrix, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, parameterized by three mixing angles $\theta_{12}$, $\theta_{23}$, $\theta_{13}$ and a CP-violating phase $\delta_{CP}$. The additional parameters governing neutrino oscillations are the squared-mass differences $\Delta m^2_{ji}=m^2_j-m^2_i$, where $m_i$ is the mass of the $i$th neutrino mass eigenstate. This review covers the rise of the PMNS three-neutrino mixing paradigm and the current status of the experimental determination of its parameters.

The next years will continue to see a rich program of experimental endeavour coming to fruition and addressing the three missing pieces of the puzzle, namely the determination of the octant and precise value of the mixing angle $\theta_{23}$, the unveiling of the neutrino mass ordering (whether $m_1 < m_2 < m_3$ or $m_3 < m_1 < m_2$) and the measurement of the CP-violating phase \dcp.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.

We review recent attempts at dealing with the sign problem in Monte Carlo calculations by deforming the region of integration in the path integral from real to complex fields. We discuss the theoretical foundations, the algorithmic issues and present some results for low dimensional field theories in both imaginary and real time.

I provide an explicit construction of spectral curves for the affine $\mathrm{E}_8$ relativistic Toda chain. Their closed form expression is obtained by determining the full set of character relations in the representation ring of $\mathrm{E}_8$ for the exterior algebra of the adjoint representation; this is in turn employed to provide an explicit construction of both integrals of motion and the action-angle map for the resulting integrable system.

I consider two main areas of applications of these constructions. On the one hand, I consider the resulting family of spectral curves in the context of the correspondences between Toda systems, 5d Seiberg-Witten theory, Gromov-Witten theory of orbifolds of the resolved conifold, and Chern-Simons theory to establish a version of the B-model Gopakumar-Vafa correspondence for the $\mathrm{sl}_N$ L\^e-Murakami-Ohtsuki invariant of the Poincar\'e integral homology sphere to all orders in $1/N$. On the other, I consider a degenerate version of the spectral curves and prove a 1-dimensional Landau-Ginzburg mirror theorem for the Frobenius manifold structure on the space of orbits of the extended affine Weyl group of type $\mathrm{E}_8$ introduced by Dubrovin-Zhang (equivalently, the orbifold quantum cohomology of the type-$\mathrm{E}_8$ polynomial $\mathbb{C} P^1$ orbifold). This leads to closed-form expressions for the flat co-ordinates of the Saito metric, the prepotential, and a higher genus mirror theorem based on the Chekhov-Eynard-Orantin recursion. I will also show how the constructions of the paper lead to a generaLISAtion of a conjecture of Norbury-Scott to ADE $\mathbb{P}^1$-orbifolds, and a mirror of the Dubrovin-Zhang construction for all Weyl groups and choices of marked roots.

We review Zimmermann's forest formula, which solves Bogoliubov's recursive $R$-operation for the subtraction of ultraviolet divergences in perturbative Quantum Field Theory. We further discuss a generaLISAtion of the $R$-operation which subtracts besides ultraviolet also Euclidean infrared divergences. This generaLISAtion, which goes under the name of the $R^*$-operation, can be used efficiently to compute renormaLISAtion constants. We will discuss several results obtained by this method with focus on the QCD beta function at five loops as well as the application to hadronic Higgs boson decay rates at N${}^4$LO. This article summarizes a talk given at the Wolfhart Zimmermann Memorial Symposium.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.