Papers for Wednesday, Mar 10 2021

X-ray reverberation mapping has emerged as a new tool to probe accretion in AGN, providing a potentially powerful probe of accretion at the black hole scale. The lags, along with relativistic spectral signatures are often interpreted in light of the lamp-post model. Focusing specifically on testing the prediction of the relativistic reverberation model, we have targeted several of the brightest Seyfert Galaxies in X-rays with different observing programs. Here, we report the results from two large campaigns with NuSATR targeting MCG-5-23-16 and SWIFT J2127.4+5654 to test the model predictions in the 3-50 keV band. These are two of three sources that showed indications of a delayed Compton hump in early data. With triple the previously analyzed exposures, we find no evidence for relativistic reverberation in MCG-5-23-16, and the energy-dependent lags are consistent with a log-linear continuum. In SWIFT J2127.4+5654, although a continuum-only model explains the data, the relativistic reverberation model provides a significant improvement to the energy and frequency-dependent lags, but with parameters that are not consistent with the time-averaged spectrum. This adds to mounting evidence showing that the lag data is not consistent with a static lamp-post model.

The nature of dark matter (DM) is unknown. One compelling possibility is DM being composed of primordial black holes (PBHs), given the tight limits on some types of elementary particles as DM. There is only one remaining window of masses available for PBHs to constitute the entire DM density, $10^{17} - 10^{23} \mathrm{\; g}$. Here, we show that the kernel population in the cold Kuiper belt rules out this window, arguing in favor of a particle nature for DM.

Testing 3D hydrodynamic models of stellar atmospheres is feasible by retrieving spectral line shapes across stellar disks, using differential spectroscopy during exoplanet transits. From synthetic data at hyper-high spectral resolution, characteristic patterns for FeI and FeII lines were identified in Paper IV from 3D models spanning T=3964-6726K (spectral types approx. K8V-F3V). The observability of patterns among lines of different strength, excitation potential and ionization level are now examined, as observed at ordinary spectral resolutions and in the presence of noise. Time variability in 3D atmospheres induces changes in spectral-line parameters, some of which are correlated. An adequate calibration could identify proxies for the jitter in apparent radial velocity to enable adjustments to actual stellar radial motion. We also examined the center-to-limb temporal variability. Recovery of spatially resolved line profiles with fitted widths and depths is shown for various noise levels. Signals during exoplanet transit are simulated. In addition to Rossiter-McLaughlin type signatures in apparent radial velocity, analogous effects are shown for line depths and widths. From exoplanet transits, overall stellar line parameters of width, depth and wavelength position can be retrieved already with moderate efforts, but a very good signal-to-noise ratio is required to reveal the more subtle signatures between subgroups of spectral lines, where finer details of atmospheric structure are encoded. In a solar model, temporal variability shows correlations between jittering in apparent radial velocity and fluctuations in line depth. Since both fluctuations in line depth and jittering in wavelength can be measured from the ground, searches for low-mass exoplanets should explore these to adjust apparent radial velocities to actual stellar motion.

We study star cluster formation in a low-metallicity environment using three dimensional hydrodynamic simulations. Starting from a turbulent cloud core, we follow the formation and growth of protostellar systems with different metallicities ranging from $10^{-6}$ to $0.1~Z_{\odot}$. The cooling induced by dust grains promotes fragmentation at small scales and the formation of low-mass stars with $M_{*} \sim 0.01$--$0.1~M_{\odot}$ when $Z/Z_{\odot} \gtrsim 10^{-5}$. While the number of low-mass stars increases with metallicity, the stellar mass distribution is still top-heavy for $Z/Z_{\odot} \lesssim 10^{-2}$ compared to the Chabrier initial mass function (IMF). In these cases, star formation begins after the turbulent motion decays and a single massive cloud core monolithically collapses to form a central massive stellar system. The circumstellar disk preferentially feeds the mass to the central massive stars, making the mass distribution top-heavy. When $Z/Z_{\odot}=0.1$, collisions of the turbulent flows promote the onset of the star formation and a highly filamentary structure develops owing to efficient fine-structure line cooling. In this case, the mass supply to the massive stars is limited by the local gas reservoir and the mass is shared among the stars, leading to a Chabrier-like IMF. We conclude that cooling at the scales of the turbulent motion promotes the development of the filamentary structure and works as an important factor leading to the present-day IMF.

Early dark energy, as a proposed solution to the Hubble tension, faces an additional "why now" problem. Why should dark energy emerge just prior to recombination, billions of years before the onset of cosmic acceleration? Assisted quintessence explains this connection by positing that multiple scaling fields build up over time to drive the present-day cosmic acceleration. In this framework, early dark energy is inevitable. Yet, we show that scaling also leads to the demise of the scenario: the same feature that solves the coincidence problem then spoils a concordance of the Hubble constant inferred from the cosmic microwave background with that from the local distance ladder. The failure of the model offers a novel lesson on the ability of new physics to resolve the Hubble tension.

We present an approach to measure the Milky Way (MW) potential using the angular accelerations of stars in aggregate as measured by astrometric surveys like Gaia. Accelerations directly probe the gradient of the MW potential, as opposed to indirect methods using e.g. stellar velocities. We show that end-of-mission Gaia stellar acceleration data may be used to measure the potential of the MW disk at approximately 3$\sigma$ significance and, if recent measurements of the solar acceleration are included, the local dark matter density at ~2$\sigma$ significance. Since the significance of detection scales steeply as $t^{5/2}$ for observing time $t$, future surveys that include angular accelerations in the astrometric solutions may be combined with Gaia to precisely measure the local dark matter density and shape of the density profile.

Galaxies are the basic structural element of the universe; galaxy formation theory seeks to explain how these structures came to be. I trace some of the foundational ideas in galaxy formation, with emphasis on the need for non-baryonic cold dark matter. Many elements of early theory did not survive contact with observations of low surface brightness galaxies, leading to the need for auxiliary hypotheses like feedback. The failure points often trace to the surprising predictive successes of an alternative to dark matter, the Modified Newtonian Dynamics (MOND). While dark matter models are flexible in accommodating observations, they do not provide the predictive capacity of MOND. If the universe is made of cold dark matter, why does MOND get any predictions right?

Dynamical friction is typically regarded a secular process, in which the subject ('perturber') evolves very slowly (secular approximation), and has been introduced to the host over a long time (adiabatic approximation). These assumptions imply that dynamical friction arises from the LBK torque with non-zero contribution only from pure resonance orbits. However, dynamical friction is only of astrophysical interest if its timescale is shorter than the age of the Universe. In this paper we therefore relax the adiabatic and secular approximations. We first derive a generalized LBK torque, which reduces to the LBK torque in the adiabatic limit, and show that it gives rise to transient oscillations due to non-resonant orbits that slowly damp out, giving way to the LBK torque. This is analogous to how a forced, damped oscillator undergoes transients before settling to a steady state, except that here the damping is due to phase mixing rather than dissipation. Next, we present a self-consistent treatment, that properly accounts for time-dependence of the perturber potential and circular frequency (memory effect), which we use to examine orbital decay in a cored galaxy. We find that the memory effect results in a phase of accelerated, super-Chandrasekhar friction before the perturber stalls at a critical radius, $R_{\mathrm{crit}}$, in the core (core-stalling). Inside of $R_{\mathrm{crit}}$ the torque flips sign, giving rise to dynamical buoyancy, which counteracts friction and causes the perturber to stall. This phenomenology is consistent with $N$-body simulations, but has thus far eluded proper explanation.

A small fraction of early-type galaxies (ETGs) show prolate rotation, i.e. they rotate around their long photometric axis. In simulations, certain configurations of galaxy mergers are known to produce this type of rotation. We investigate the association of prolate rotation and signs of galaxy interactions among the observed galaxies. We collected a sample of 19 nearby ETGs with distinct prolate rotation from the literature and inspected their ground-based deep optical images for interaction signs - 18 in archival images and one in a new image obtained with the Milankovi\'c telescope. Tidal tails, shells, asymmetric/disturbed stellar halos, or on-going interactions are present in all the 19 prolate rotators. Comparing this with the frequency of tidal disturbance among the general sample of ETGs of a roughly similar mass range and surface-brightness limit, we estimate that the chance probability of such an observation is only 0.00087. We also found a significant overabundance of prolate rotators that are hosting multiple stellar shells. The visible tidal features imply a relatively recent galaxy interaction. That agrees with the Illustris large-scale cosmological hydrodynamical simulation, where prolate rotators are predominantly formed in major mergers during the last 6 Gyr. In the appendix, we present the properties of an additional galaxy, NGC 7052, a prolate rotator for which no deep images are available, but for which an HST image revealed the presence of a prominent shell, which had not been reported before.

A significant fraction of baryons in galaxies are in the form of diffuse gas of the circumgalactic medium (CGM). One critical component of the multi-phases of CGM, the so-called "coronal" warm-hot phase gas ($\rm 10^{5}-10^{6}$ K) traced by O VI 1032, 1038 \r{A} resonance lines, has rarely been detected in emission from galaxy halos other than Milky Way. Here we report four additional detections of O VI emission gas in the halos of nearby edge-on galaxies, NGC 4631 and NGC 891, using archival Far Ultraviolet Spectroscopic Explorer data and an updated data pipeline. We find the most intense O VI emission to be from fields forming a vertical line near the center of NGC 4631, despite the close proximity to the disk of two other fields. We theorize that the vertically-oriented fields trace a warm-hot phase gas filament that is possibly infalling onto the galaxy. The O VI kinematics of the fields closer to the edges of the disk suggest that those fields sample ejected gas that has shock heated and is cooling before returning to the disk. These results point to the existence of filamentary structures of more intense O VI emission superimposed within a more diffuse and faint O VI halo in star-forming galaxies. Volume-filled O VI emission mapping is greatly needed to determine the structure and prevalence of warm-hot gas and the role it plays in the cycling of gas between the galaxy disk and the halo. Finally, we present the sensitivity of future funded and proposed UV missions (LUVOIR-A, LUVOIR-B, CETUS, and Aspera) to the detection of diffuse and faint O VI emission in nearby galaxy halos.

We present a panchromatic study of 11 (sub-)millimetre selected DSFGs with spectroscopically confirmed redshift ($1.5< z_{\rm spec}<3$) in the GOODS-S field, with the aim of constraining their astrophysical properties (e.g., age, stellar mass, dust and gas content) and characterizing their role in the context of galaxy evolution. The multi-wavelength coverage of GOODS-S, from X-rays to radio band, allow us to model galaxy SED by using CIGALE with a novel approach, based on a physical motivated modelling of stellar light attenuation by dust. Median stellar mass ($\simeq6.5\times10^{10}$ M$_\odot$) and SFR ($\simeq241$ M$_\odot$ yr$^{-1}$) are consistent with galaxy main-sequence at $z\sim2$. The galaxies are experiencing an intense and dusty burst of star formation (median L$_{\rm IR}\simeq2\times10^{12}$ L$_\odot$), with a median age of $750$ Myr. The high median content of interstellar dust (M$_{\rm dust}\simeq5\times10^8$ M$_\odot$) suggests a rapid enrichment of the ISM (on timescales $\sim10^8$ yr). We derived galaxy total and molecular gas content from CO spectroscopy and/or Rayleigh-Jeans dust continuum ($10^{10}\lesssim$ M$_{\rm gas}/$M$_\odot\lesssim10^{11}$), depleted over a typical timescale $\tau_{\rm depl}\sim200$ Myr. X-ray and radio luminosities suggest that most of the galaxies hosts an accreting radio silent/quiet SMBH. This evidence, along with their compact multi-wavelength sizes (median r$_{\rm ALMA}\sim$ r$_{\rm VLA}=1.8$ kpc, r$_{\rm HST}=2.3$ kpc) measured from high-resolution imaging ($\theta_{\rm res}\lesssim$ 1 arcsec), indicates these objects as the high-z star-forming counterparts of massive quiescent galaxies, as predicted e.g., by the in-situ scenario. Four objects show some signatures of a forthcoming/ongoing AGN feedback, that is thought to trigger the morphological transition from star-forming disks to ETGs.

We explore hierarchical black hole (BH) mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs), accounting for both original and dynamically assembled binary BHs (BBHs). We find that the median mass of both first- and nth-generation dynamical mergers is larger in GCs and YSCs with respect to NSCs, because the lighter BHs are ejected by supernova kicks from the lower-mass clusters. Also, first- and nth-generation BH masses are strongly affected by the metallicity of the progenitor stars: the median mass of the primary BH of a nth-generation merger is $\sim{}24-38$ M$_\odot$ ($\sim{}9-15$ M$_\odot$) in metal-poor (metal-rich) NSCs. The maximum BH mass mainly depends on the escape velocity: BHs with mass up to several thousand M$_\odot$ form in NSCs, while YSCs and GCs host BHs with mass up to several hundred M$_\odot$. Furthermore, we calculate the fraction of mergers with at least one component in the pair-instability mass gap ($f_{\rm PI}$) and in the intermediate-mass BH regime ($f_{\rm IMBH}$). In the fiducial model for dynamical BBHs with metallicity $Z=0.002$, we find $f_{\rm PI}\approx{}0.05$, $0.02$ and $0.007$ ($f_{\rm IMBH}\approx{}0.012$, $0.002$ and $0.001$) in NSCs, GCs and YSCs, respectively. Both $f_{\rm PI}$ and $f_{\rm IMBH}$ drop by at least one order of magnitude at solar metallicity. Finally, we investigate the formation of GW190521 by assuming that it is either a nearly equal-mass BBH or an intermediate-mass ratio inspiral.

We consider the application of a machine learning technique to growth and expansion rate data in the context of teleparallel gravity (TG). We do this by using a combined approach of Hubble data together with redshift-space-distortion $f\sigma_8$ data which together are used to reconstruct the TG Lagrangian via Gaussian processes (GP), where the Hubble data mainly comes from cosmic chronometer and supernova type Ia data from the Pantheon release. In this work, we consider two main GP covariance functions, namely the squared-exponential and Cauchy kernels in order to show consistency (to within 1$\sigma$ uncertainties). The core results of this work are the numerical constructions of the TG Lagrangian from GP reconstructed Hubble and growth data. We take different possible combinations of the datasets and kernels to show any potential differences in this regard. We show that nontrivial cosmology beyond $\Lambda$CDM falls within the uncertainties of the reconstruction from growth data.

This paper presents the main features of a new and updated version of the program PArthENoPE, which the community has been using for many years for computing the abundances of light elements produced during Big Bang Nucleosynthesis. This is the third release of the PArthENoPE code, after the 2008 and the 2018 ones, and will be distributed from the code's website, this http URL Apart from minor changes, the main improvements in this new version include a revisited implementation of the nuclear rates for the most important reactions of deuterium destruction, H2(p,gamma)He3, H2(d, n)He3 and H2(d, p)H3, and a re-designed GUI, which extends the functionality of the previous one. The new GUI, in particular, supersedes the previous tools for running over grids of parameters with a better management of parallel runs, and it offers a brand-new set of functions for plotting the results.

As ancient, gravitationally bound stellar populations, globular clusters are abundant, vibrant laboratories characterized by high frequencies of dynamical interactions coupled to complex stellar evolution. Using surface brightness and velocity dispersion profiles from the literature, we fit $59$ Milky Way globular clusters to dynamical models from the \texttt{CMC Cluster Catalog}. Without doing any interpolation, and without any directed effort to fit any particular cluster, $26$ globular clusters are well-matched by at least one of our models. We discuss in particular the core-collapsed clusters NGC 6293, NGC 6397, NGC 6681, and NGC 6624, and the non-core-collapsed clusters NGC 288, NGC 4372, and NGC 5897. As NGC 6624 lacks well-fitting snapshots on the main \texttt{CMC Cluster Catalog}, we run six additional models in order to refine the fit. We calculate metrics for mass segregation, explore the production of compact object sources such as millisecond pulsars, cataclysmic variables, low-mass X-ray binaries, and stellar-mass black holes, finding reasonable agreement with observations. Additionally, closely mimicking observational cuts, we extract the binary fraction from our models, finding good agreement except in the dense core regions of core-collapsed clusters. Accompanying this paper are a number of \textsf{python} methods for examining the publicly accessible \texttt{CMC Cluster Catalog}, as well as any other models generated using \texttt{CMC}.

Model-independent constraints on modified gravity models hitherto exist mainly on linear scales. A recently developed formalism presented a consistent parameterisation that is valid on all scales. Using this approach, we perform model-independent modified gravity $N$-body simulations on all cosmological scales with a time-dependent $\mu$. We present convergence tests of our simulations, and we examine how well existing fitting functions reproduce the non-linear matter power spectrum of the simulations. We find that although there is a significant variation in the accuracy of all of the fitting functions over the parameter space of our simulations, the ReACT framework delivers the most consistent performance for the matter power spectrum. We comment on how this might be improved to the level required for future surveys such as Euclid and the Vera Rubin Telescope (LSST). We also show how to compute weak-lensing observables consistently from the simulated matter power spectra in our approach, and show that ReACT also performs best when fitting the weak-lensing observables. This paves the way for a full model-independent test of modified gravity using all of the data from such upcoming surveys.

WR 25 is a colliding-wind binary star system comprised of a very massive O2.5If*/WN6 primary and an O-star secondary in a 208-day period eccentric orbit. These hot stars have strong, highly-supersonic winds which interact to form a bright X-ray source from wind-collision-shocks whose conditions change with stellar separation. Different views through the WR and O star winds are afforded with orbital phase as the stars move about their orbits, allowing for exploration of wind structure in ways not easy or even possible for single stars. We have analyzed an on-axis Chandra/HETGS spectrum of WR 25 obtained shortly before periastron when the X-rays emanating from the system are the brightest. From the on-axis observations, we constrain the line fluxes, centroids, and widths of various emission lines, including He-triplets of Si XIII and Mg XI. We have also been able to include several serendipitous off-axis HETG spectra from the archive and study their flux variation with phase. This is the first report on high-resolution spectral studies of WR 25 in X-rays.

We examine the impact of baryonic physics on the halo distribution in hydrodynamic simulations (Illustris, IllustrisTNG, and EAGLE), particularly with regards to how it differs from that in dark matter only (DMO) simulations. We find that, in general, DMO simulations produce halo mass functions (HMF) that are shifted to higher halo masses than their hydrodynamic counterparts, due to the lack of baryonic physics. However, the exact nature of this mass shift is a complex function of mass, halo definition, redshift, and larger-scale environment, and it also depends on the specifics of the baryonic physics implemented in the simulation. We present fitting formulae for the corrections one would need to apply to each DMO halo catalogue in order to reproduce the HMF found in its hydrodynamic counterpart. We provide these formulae for all three simulations, for five different halo definitions at redshifts 0, 1, and 2. Additionally, we explore the dependence on environment of this HMF discrepancy, and find that, in most cases, halos in low density environments are slightly more impacted by baryonic physics than halos in high density environments. We thus also provide environment-dependent mass correction formulae that can reproduce the conditional, as well as global, HMF. We show that our mass corrections also repair the large-scale clustering of halos, though the environment-dependent corrections are required to achieve an accuracy better than 2%. Finally, we examine the impact of baryonic physics on the halo mass - concentration relation, and find that its slope in hydrodynamic simulations is consistent with that in DMO simulations. Ultimately, we recommend that any future work relying on DMO halo catalogues incorporate our mass corrections to test the robustness of their results to baryonic effects.

One of the last unexplored windows to the cosmos, the Dark Ages and Cosmic Dawn, can be opened using a simple low frequency radio telescope from the stable, quiet lunar farside to measure the Global 21-cm spectrum. This frontier remains an enormous gap in our knowledge of the Universe. Standard models of physics and cosmology are untested during this critical epoch. The messenger of information about this period is the 1420 MHz (21-cm) radiation from the hyperfine transition of neutral hydrogen, Doppler-shifted to low radio astronomy frequencies by the expansion of the Universe. The Global 21-cm spectrum uniquely probes the cosmological model during the Dark Ages plus the evolving astrophysics during Cosmic Dawn, yielding constraints on the first stars, on accreting black holes, and on exotic physics such as dark matter-baryon interactions. A single low frequency radio telescope can measure the Global spectrum between ~10-110 MHz because of the ubiquity of neutral hydrogen. Precise characterizations of the telescope and its surroundings are required to detect this weak, isotropic emission of hydrogen amidst the bright foreground Galactic radiation. We describe how two antennas will permit observations over the full frequency band: a pair of orthogonal wire antennas and a 0.3-m$^3$ patch antenna. A four-channel correlation spectropolarimeter forms the core of the detector electronics. Technology challenges include advanced calibration techniques to disentangle covariances between a bright foreground and a weak 21-cm signal, using techniques similar to those for the CMB, thermal management for temperature swings of >250C, and efficient power to allow operations through a two-week lunar night. This simple telescope sets the stage for a lunar farside interferometric array to measure the Dark Ages power spectrum.

The D68 ringlet is the innermost feature in Saturn's rings. Four clumps that appeared in D68 around 2014 remained evenly spaced about 30 degrees apart and moved very slowly relative to each other from 2014 up until the last measurements were taken in 2017. D68's narrowness and the distribution of clumps could either indicate that we have a collection of source bodies in a co-orbital configuration or imply that an outside force confines the observed dust and any source bodies. In this paper we explore the possibility that these four clumps arose from four source bodies in a co-orbital configuration. We find that there are no solutions with four masses that produce the observed spacings. We therefore consider whether an unseen fifth co-orbital object could account for the discrepancies in the angular separations and approach a stable stationary configuration. We find a range of solutions for five co-orbital objects where their mass ratios depend on the assumed location of the fifth mass. Numerical simulations of five co-orbitals are highly sensitive to initial conditions, especially for the range of masses we would expect the D68 clumps to have. The fragility of our D68 co-orbital system model implies that there is probably some outside force confining the material in this ringlet.

We investigate the structure and the tidal deformability of strange stars (SSs) with a mirror-dark-matter (MDM) core for the standard MIT bag model. We find that to explain the observations of PSR J0740+6620, PSR J0030+0451 and GW170817 simultaneously, SSs in GW170817 should have a MDM core although it is unnecessary for PSR J0740+6620 and PSR J0030+0451 to contain a MDM core. Our study leads to the result that for the standard MIT bag model, the observations of compact stars mentioned above confirm the existence of a dark-matter core inside SSs.

We use the concept of co-added rotation curves of Salucci et al. to investigate the properties of axi-symmetric multistate Scalar Field Dark Matter halos in low surface brightness galaxies and dwarf disc galaxies. We fit their rotation curves in two-state configurations and we find that all of these can be well fitted with a particle mass $\mu \sim (10^{-23} - 10^{-24})\rm{eV}/c^2$. Comparing our results with the standard cosmological model, the well-known $\Lambda$-cold dark matter, by using the Bayesian information criterion and the Akaike information criterion, we found that our two-state model seemed to be preferred.

We introduce the Australia Telescope Compact Array (ATCA) rapid-response mode by presenting the first successful trigger on the short-duration gamma-ray burst (GRB) 181123B. Early-time radio observations of short GRBs may provide vital insights into the radio afterglow properties of Advanced LIGO- and Virgo-detected gravitational wave events, which will in turn inform follow-up strategies to search for counterparts within their large positional uncertainties. The ATCA was on target within 12.6 hr post-burst, when the source had risen above the horizon. While no radio afterglow was detected during the 8.3 hr observation, we obtained force-fitted flux densities of $7 \pm 12$ and $15 \pm 11~\mu$Jy at 5.5 and 9 GHz, respectively. Afterglow modelling of GRB 181123B showed that the addition of the ATCA force-fitted radio flux densities to the Swift X-ray Telescope detections provided more stringent constraints on the fraction of thermal energy in the electrons (log$\epsilon_e = -0.75^{+0.39}_{-0.40}$ rather than log$\epsilon_e = -1.13^{+0.82}_{-1.2}$ derived without the inclusion of the ATCA values), which is consistent with the range of typical $\epsilon_e$ derived from GRB afterglow modelling. This allowed us to predict that the forward shock may have peaked in the radio band $\sim10$ days post-burst, producing detectable radio emission $\gtrsim3-4$ days post-burst. Overall, we demonstrate the potential for extremely rapid radio follow-up of transients and the importance of triggered radio observations for constraining GRB blast wave properties, regardless of whether there is a detection, via the inclusion of force-fitted radio flux densities in afterglow modelling efforts.

In this paper we analyze the implications of gravitational waves (GWs) as standard sirens on the modified gravity models by using the third-generation gravitational wave detector, i.e., the Einstein Telescope. Two viable models in $f(R)$ theories within the Palatini formalism are considered in our analysis ($f_{1}(\mathcal{R})=\mathcal{R}-\frac{\beta}{\mathcal{R}^{n}}$ and $f_{2}(\mathcal{R})=\mathcal{R}+\alpha\ln{\mathcal{R}}-\beta$), with the combination of simulated GW data and the latest electromagnetic (EM) observational data (including the recently released Pantheon type Ia supernovae sample, the cosmic chronometer data, and baryon acoustic oscillation distance measurements). Our analysis reveals that the standard sirens GWs, which provide an independent and complementary alternative to current experiments, could effectively eliminate the degeneracies among parameters in the two modified gravity models. In addition, we thoroughly investigate the nature of geometrical dark energy in the modified gravity theories with the assistance of $Om(z)$ and statefinder diagnostic analysis. The present analysis makes it clear-cut that the simplest cosmological constant model is still the most preferred by the current data. However, the combination of future naturally improved GW data most recent EM observations will reveal the consistency or acknowledge the tension between the $\Lambda$CDM model and modified gravity theories.

Jayasinghe et al. (2021) identified a dark $\approx 3\,M_\odot$ companion on a nearly edge-on $\approx 60\,\mathrm{day}$ orbit around the red giant star V723 Monoceros as a black hole candidate in the mass gap. This scenario was shown to explain most of the data presented by Jayasinghe et al. (2021), except for periodic radial velocity (RV) residuals from the circular Keplerian model. Here we show that the RV residuals are explained by orbital phase-dependent distortion of the absorption line profile associated with changing visible fractions of the approaching and receding sides of the red giant star, whose surface is tidally deformed by and rotating synchronously with the dark companion. Our RV model constrains the companion mass $M_\bullet = 2.95\pm0.17\,M_\odot$ and orbital inclination $i=82.9^{+7.0}_{-3.3}\,\mathrm{deg}$ (medians and 68.3% highest density intervals of the marginal posteriors) adopting the radius of the red giant $24.0\pm0.9\,R_\odot$ as constrained from its SED and distance. The analysis provides independent support for the companion mass from ellipsoidal variations and the limits on the companion's luminosity from the absence of eclipses both derived by Jayasinghe et al. (2021). We also show that a common scheme to evaluate the tidal RV signal as the flux-weighted mean of the surface velocity field can significantly underestimate its amplitude for RVs measured with a cross-correlation technique, and present a modified prescription that directly models the distorted line profile and its effects on the measured RVs. The formulation will be useful for estimating the component masses and inclinations in other similar binaries.

Understanding how massive stars die as supernovae is a crucial question in modern astrophysics. Supernovae are powerful stellar explosions and key drivers in the cosmic baryonic cycles by injecting their explosion energy and heavy elements to the interstellar medium that forms new stars. After decades of effort, astrophysicists have built up a stand model for the explosion mechanism of massive stars. However, this model is challenged by new kinds of stellar explosions discovered in the recent transit surveys. In particular, the new population called superluminous supernovae, which are a hundred times brighter than typical supernovae, is revolutionizing our understanding of supernovae. New studies suggest the superluminous supernovae are associated with the unusual demise of very massive stars and their extreme supernovae powered by the radioactive isotopes or compact objects formed after the explosion. Studying these supernovae fills a gap of knowledge between the death of massive stars and their explosions; furthermore, we may apply their intense luminosity to light up the distant universe. This paper aims to provide a timely review of superluminous supernovae physics, focusing on the latest development of their theoretical models.

The bright blazar OJ 287 is the best-known candidate for hosting a supermassive black hole binary system. It inspirals due to the emission of nanohertz gravitational waves (GWs). Observations of historical and predicted quasi-periodic high-brightness flares in its century-long optical lightcurve, allow us to determine the orbital parameters associated with the binary black hole (BBH) central engine. In contrast, the radio jet of OJ 287 has been covered with Very Long Baseline Interferometry (VLBI) observations for only about $30$ years and these observations reveal that the position angle (PA) of the jet exhibits temporal variations at both millimetre and centimetre wavelengths. Here we associate the observed PA variations in OJ 287 with the precession of its radio jet. In our model, the evolution of the jet direction can be associated either with the primary black hole (BH) spin evolution or with the precession of the angular momentum direction of the inner region of the accretion disc. Our Bayesian analysis shows that the BBH central engine model, primarily developed from optical observations, can also broadly explain the observed temporal variations in the radio jet of OJ 287 at frequencies of 86, 43, and 15 GHz. Ongoing Global mm-VLBI Array (GMVA) observations of OJ 287 have the potential to verify our predictions for the evolution of its $86$ GHz PA values. Additionally, thanks to the extremely high angular resolution that the Event Horizon Telescope (EHT) can provide, we explore the possibility to test our BBH model through the detection of the jet in the secondary black hole.

We generalize the stochastic theory of hierarchical clustering presented in paper I by Lapi & Danese (2020) to derive the (conditional) halo progenitor mass function and the related large-scale bias. Specifically, we present a stochastic differential equation that describes fluctuations in the mass growth of progenitor halos of given descendant mass and redshift, as driven by a multiplicative Gaussian white noise involving the power spectrum and the spherical collapse threshold of density perturbations. We demonstrate that, as cosmic time passes, the noise yields an average drift of the progenitors toward larger masses, that quantitatively renders the expectation from the standard extended Press & Schechter (EPS) theory. We solve the Fokker-Planck equation associated to the stochastic dynamics, and obtain as an exact, stationary solution the EPS progenitor mass function. Then we introduce a modification of the stochastic equation in terms of a mass-dependent collapse threshold modulating the noise, and solve analytically the associated Fokker-Planck equation for the progenitor mass function. The latter is found to be in excellent agreement with the outcomes of $N-$body simulations; even more remarkably, this is achieved with the same shape of the collapse threshold used in paper I to reproduce the halo mass function. Finally, we exploit the above results to compute the large-scale halo bias, and find it in pleasing agreement with the $N-$body outcomes. All in all, the present paper illustrates that the stochastic theory of hierarchical clustering introduced in paper I can describe effectively not only halos' abundance, but also their progenitor distribution and their correlation with the large-scale environment across cosmic times.

In the search for small exoplanets orbiting cool stars whose spectral energy distributions peak in the near infrared, the strong absorption of radiation in this region due to water vapour in the atmosphere is a particularly adverse effect for the ground-based observations of cool stars. To achieve the photometric precision required to detect exoplanets in the near infrared, it is necessary to mitigate the impact of variable precipitable water vapour (PWV) on radial-velocity and photometric measurements. The aim is to enable global PWV correction by monitoring the amount of precipitable water vapour at zenith and along the line of sight of any visible target. We developed an open source Python package that uses Geostationary Operational Environmental Satellites (GOES) imagery data, which provides temperature and relative humidity at different pressure levels to compute near real-time PWV above any ground-based observatory covered by GOES every 5 minutes or 10 minutes depending on the location. We computed PWV values on selected days above Cerro Paranal (Chile) and San Pedro M\'artir (Mexico) to benchmark the procedure. We also simulated different pointing at test targets as observed from the sites to compute the PWV along the line of sight. To asses the accuracy of our method, we compared our results with the on-site radiometer measurements obtained from Cerro Paranal. Our results show that our publicly-available code proves to be a good supporting tool for measuring the local PWV for any ground-based facility within the GOES coverage, which will help in reducing correlated noise contributions in near-infrared ground-based observations that do not benefit from on-site PWV measurements.

Variations in the grain size distribution are to be expected in the interstellar medium (ISM) due to grain growth and destruction. In this work, we present a dust collision model to be implemented inside a magnetohydrodynamical (MHD) code that takes into account grain growth and shattering of charged dust grains of a given composition (silicate or graphite). We integrate this model in the MHD code Athena, and builds on a previous implementation of the dynamics of charged dust grains in the same code. To demonstrate the performance of this coagulation model, we study the variations in the grain size distribution of a single-sized population of dust with radius 0.05 $\mu$m inside several dust filaments formed during a 2D MHD simulation. We also consider a realistic dust distribution with sizes ranging from 50 \AA~to 0.25 $\mu$m and analyze both the variations in the size distribution for graphite and silicates, as well as of the far ultraviolet extinction curve. From the obtained results, we conclude that the methodology here presented, based on the MHD evolution of the equation of motion for a charged particle, is optimal for studying the coagulation of charged dust grains in a diffuse regime such as a molecular cloud envelope. Observationally, these variations in the dust size distribution are translated into variations in the far ultraviolet extinction curve, and they are mainly caused by small graphite dust grains.

The focus of this article is a re-count of Richard Carrington's original sunspot observations from his book drawings (Carrington,1863) by an observer from the World data Center-SILSO network, Thomas H. Teague(UK). This modern recount will enable the use of Carrington's observations in the re-computation of the entire Sunspot Number series in a way Carrington's original counts (Casas and Vaquero,2014) did not. here we present comparison studies of the new recounted series with contemporary observations, new data extracted from the Journals of the Zurich Observatory and other sources of Carrington's own observations and conclude that Carrington's group counting is very close to the modern way of counting while his method for counting individual spots lags significantly behind modern counts. We also test the quality and robustness of the new re-count with methods developed in Mathieu et al.,2019.

Based on the SDSS data, the properties of galaxies with quenched star formation (QGs) within the "splashback"-radius of galaxy clusters $R_{\rm sp}$ and beyond it have been studied. We used a sample of 40 groups and galaxy clusters and a sample of field galaxies at $0.02<z<0.045$. The radii $R_{\rm sp}$ were defined from the observed integrated distribution of the number of galaxies as a function of the squared distance from the center of the galaxy systems. We show that in galaxy clusters 72% of the QGs we have found are within $R_{\rm sp}$. About 40% of these galaxies are late-type ones with ${frac}DeV < 0.8$. Approximately 80% of galaxies with quenched star formation have stellar masses in the range of $\log M_*/M_{\odot} = [10; 11]$. We found that QGs of late types and of early types in a less degree have maximum angular radii $R_{90,r}$ and $R_{50,r}$ near the "splashback"-radius of groups and clusters of galaxies. Our results confirm the assumption that in the filaments directed toward clusters, the quenched galaxies are more massive near the boundaries of clusters of galaxies than at the outskirts.

The perplexing mystery of what maintains the solar coronal temperature at about a million K, while the visible disc of the Sun is only at 5800 K, has been a long standing problem in solar physics. A recent study by Mondal(2020) has provided the first evidence for the presence of numerous ubiquitous impulsive emissions at low radio frequencies from the quiet sun regions, which could hold the key to solving this mystery. These features occur at rates of about five hundred events per minute, and their strength is only a few percent of the background steady emission. One of the next steps for exploring the feasibility of this resolution to the coronal heating problem is to understand the morphology of these emissions. To meet this objective we have developed a technique based on an unsupervised machine learning approach for characterising the morphology of these impulsive emissions. Here we present the results of application of this technique to over 8000 images spanning 70 minutes of data in which about 34,500 features could robustly be characterised as 2D elliptical Gaussians.

We study the outskirts ($R<3R_{200c}$) of 40 groups and clusters of galaxies of the local Universe ($0.02<z<0.045$) with 300~km~s$^{-1}<\sigma<950$~km~s$^{-1}$. Using the SDSS DR10 catalog data, we measured the stellar mass of galaxy clusters in accordance with the previously determined $K_s$-luminosity (2MASX data) and found their correlation in the form $M_*/M_{\odot} \propto (L_K/L_{\odot})^{1.010\pm0.004}$ ($M_K<-21\,.\!\!^{\rm m}5$, $R<R_{200c}$). We also found the dependence of the galaxy cluster stellar mass on halo mass: $M_*/M_{\odot} \propto (M_{200c}/M_{\odot})^{0.77\pm0.01}$}. Our results show that the fraction of galaxies with quenched star formation ($M_K<-21\,.\!\!^{\rm m}5$) is maximal in the central regions of the galaxy clusters and equals, on the average, $0.81\pm0.02$; it decreases to $0.44\pm0.02$ outside of the projected radius $R_{sp}$ ($2<R/R_{200c}<3$), which we found from the observed profile, but still remains higher than that in the field by 27\%. The fraction of early-type "red sequence" galaxies decreases from $0.54\pm0.02$ in the center to $0.24\pm0.01$ beyond $R_{\rm sp}$, reaching its field value.

The detection of a wide range of substructures such as rings, cavities and spirals has become a common outcome of high spatial resolution imaging of protoplanetary disks, both in the near-infrared scattered light and in the thermal millimetre continuum emission. The most frequent interpretation of their origin is the presence of planetary-mass companions perturbing the gas and dust distribution in the disk (perturbers), but so far the only bona-fide detection has been the two giant planets around PDS 70. Here, we collect a sample of 15 protoplanetary disks showing substructures in SPHERE scattered light images and present a homogeneous derivation of planet detection limits in these systems. We also estimate the mass of these perturbers through a Hill radius prescription and a comparison to ALMA data. Assuming that one single planet carves each substructure in scattered light, we find that more massive perturbers are needed to create gaps within cavities than rings, and that we might be close to a detection in the cavities of RX J1604, RX J1615, Sz Cha, HD 135344B and HD 34282. We reach typical mass limits in these cavities of 3-10 Mjup. For planets in the gaps between rings, we find that the detection limits of SPHERE are about an order of magnitude away in mass, and that the gaps of PDS 66 and HD 97048 seem to be the most promising structures for planet searches. The proposed presence of massive planets causing spiral features in HD 135344B and HD 36112 are also within SPHERE's reach assuming hot-start models.These results suggest that current detection limits are able to detect hot-start planets in cavities, under the assumption that they are formed by a single perturber located at the centre of the cavity. More realistic planet mass constraints would help to clarify whether this is actually the case, which might point to perturbers not being the only way of creating substructures.

The galaxy cluster MS 0735.6+7421 hosts two large X-ray cavities, filled with radio emission, where a decrease of the Sunyaev-Zel'dovich (SZ) effect has been detected, without establishing if its origin is thermal (from a gas with very high temperature) or non-thermal. In this paper we study how thermal and non-thermal contributions to the SZ effect in the cavities are related; in fact, Coulomb interactions with the thermal gas modify the spectrum of low energy non-thermal electrons, which dominate the non-thermal SZ effect; as a consequence, the intensity of the non-thermal SZ effect is stronger for lower density of the thermal gas inside the cavity. We calculate the non-thermal SZ effect in the cavities as a function of the thermal density, and compare the SZ effects produced by thermal and non-thermal components, and with the one from the external Intra Cluster Medium (ICM), searching for the best frequency range where it is possible to disentangle the different contributions. We find that for temperatures inside the cavities higher than $\sim1500$ keV the non-thermal SZ effect is expected to dominate on the thermal one, particularly at high frequencies ($\nu>500$ GHz), where it can also be a non-negligible fraction of the SZ effect from the external ICM. We also discuss the possible sources of astrophysical bias (as kinetic SZ effect and foreground emission from Galactic dust) and possible ways to address them, as well as necessary improvements in the modeling of the properties of cavities and the ICM.

We investigate the chaotic spin-down behavior seen from some pulsars in terms of the nonlinear superfluid dynamics. To this end, we numerically solve the set of equations for the superfluid-normal matter system whose coupling is mediated by creep of the vortex lines. We show that glitch perturbations which introduce a time-delay to the steady-state dynamics leave behind a remnant in the third time derivative of the rotational phase. This time-delay induces a hyper-chaotic spin-down for pulsars. We find that glitch-induced changes in the rotational parameters lead to non-closing cyclic patterns in the time-delayed phase difference diagram. We observe that the number of cycles, $N$, in the diagram results from $N+1$ glitches occurred in total observation time.

The properties of the first-discovered interstellar object (ISO), 1I/2017 (`Oumuamua), differ from both Solar System asteroids and comets, casting doubt on a protoplanetary disk origin. In this study, we investigate the possibility that it formed with a substantial H2 ice component in the starless core of a giant molecular cloud. While interstellar solid hydrogen has yet to be detected, this constituent would explain a number of the ISO's properties. We consider the relevant processes required to build decameter-sized, solid hydrogen bodies and assess the plausibility of growth in various size regimes. Via an energy balance argument, we find that the most severe barrier to formation is the extremely low temperature required for the favorability of molecular hydrogen ice. However, if deposition occurs, we find that the turbulence within starless cores is conducive for growth into kilometer-sized bodies on sufficiently short timescales. Then, we analyze mass loss in the interstellar medium and determine the necessary size for a hydrogen object to survive a journey to the Solar System as a function of ISO age. Finally, we discuss the implications if the H2 explanation is correct, and we assess the future prospects of ISO science. If hydrogen ice ISOs do exist, our hypothesized formation pathway would require a small population of porous, 100 micron dust in a starless core region that has cooled to 2.8K via adiabatic expansion of the surrounding gas and excellent shielding from electromagnetic radiation and cosmic rays.

Context. Turbulent transport in stellar radiative zones is a key ingredient of stellar evolution theory, but the anisotropy of the transport due to the stable stratification and the rotation of these regions is poorly understood. The assumption of shellular rotation, which is a cornerstone of the so-called rotational mixing, relies on an efficient horizontal transport. However, this transport is included in many stellar evolution codes through phenomenological models that have never been tested. Aims. We investigate the impact of horizontal shear on the anisotropy of turbulent transport. Methods. We used a relaxation approximation (also known as {\tau} approximation) to describe the anisotropising effect of stratification, rotation, and shear on a background turbulent flow by computing velocity correlations. Results. We obtain new theoretical scalings for velocity correlations that include the effect of horizontal shear. These scalings show an enhancement of turbulent motions, which would lead to a more efficient transport of chemicals and angular momentum, in better agreement with helio- and asteroseismic observations of rotation in the whole Hertzsprung-Russell diagram. Moreover, we propose a new choice for the non-linear time used in the relaxation approximation, which characterises the source of the turbulence. Conclusions. For the first time, we describe the effect of stratification, rotation, and vertical and horizontal shear on the anisotropy of turbulent transport in stellar radiative zones. The new prescriptions need to be implemented in stellar evolution calculations. To do so, it may be necessary to implement non-diffusive transport.

The ESO workshop "Ground-based thermal infrared astronomy" was held on-line October 12-16, 2020. Originally planned as a traditional in-person meeting at ESO in Garching in April 2020, it was rescheduled and transformed into a fully on-line event due to the COVID-19 pandemic. With 337 participants from 36 countries the workshop was a resounding success, demonstrating the wide interest of the astronomical community in the science goals and the toolkit of ground-based thermal infrared astronomy.

Self-gravitating disks are believed to play an important role in astrophysics in particular regarding the star and planet formation process. In this context, disks subject to an idealized cooling process, characterized by a cooling timescale $\beta$ expressed in unit of orbital timescale, have been extensively studied. We take advantage of the Riemann solver and the 3D Godunov scheme implemented in the code Ramses to perform high resolution simulations, complementing previous studies that have used Smoothed Particle Hydrodynamics (SPH) or 2D grid codes. We observe that the critical value of $\beta$ for which the disk fragments is consistent with most previous results, and is not well converged with resolution. By studying the probability density function of the fluctuations of the column density ($\Sigma$-PDF), we argue that there is no strict separation between the fragmented and the unfragmented regimes but rather a smooth transition with the probability of apparition of fragments steadily diminishing as the cooling becames less effective. We find that the high column density part of the $\Sigma$-PDF follows a simple power law whose slope turns out to be proportional to $\beta$ and we propose an explanation based on the balance between cooling and heating through gravitational stress. Our explanation suggests that a more efficient cooling requires more heating implying a larger fraction of dense material which, in the absence of characteristic scales, results in a shallower scale-free power law. We propose that the gravitational cascade proceeds in two steps, first the formation of a dense filamentary spiral pattern through a sequence of quasi-static equilibrium triggered by the viscous transport of angular momentum, and second the collapse alongside these filaments that eventually results in the formation of bounded fragments.

The radio nebula W50 is a unique object interacting with the jets of the microquasar SS433. The SS433/W50 system is a good target for investigating the energy of cosmic-ray particles accelerated by galactic jets. We report observations of radio nebula W50 conducted with the NSF's Karl G. Jansky Very Large Array (VLA) in the L band (1.0 -- 2.0 GHz). We investigate the secular change of W50 on the basis of the observations in 1984, 1996, and 2017, and find that most of its structures were stable for 33 years. We revise the upper limit velocity of the eastern terminal filament by half to 0.023$c$ assuming a distance of 5.5 kpc. We also analyze the observational data of the Arecibo Observatory 305-m telescope and identify the HI cavity around W50 in the velocity range 33.77 km s$^{-1}$ -- 55.85 km s$^{-1}$. From this result, we estimate the maximum energy of the cosmic-ray protons accelerated by the jet terminal region to be above 10$^{15.5}$ eV. We also use the luminosity of the gamma-rays in the range 0.5 -- 10 GeV to estimate the total energy of accelerated protons below 5.2 $\times$ 10$^{48}$ erg.

Upcoming measurements of the small-scale primary cosmic microwave background (CMB) temperature and polarization power spectra ($TT$/$TE$/$EE$) are anticipated to yield transformative constraints on new physics, including the effective number of relativistic species in the early universe ($N_{\rm eff}$). However, at multipoles $\ell \gtrsim 3000$, the primary CMB power spectra receive significant contributions from gravitational lensing. While these modes still carry primordial information, their theoretical modeling requires knowledge of the CMB lensing convergence power spectrum, $C_L^{\kappa\kappa}$, including on small scales where it is affected by nonlinear gravitational evolution and baryonic feedback processes. Thus, the high-$\ell$ primary CMB is sensitive to these late-time, nonlinear effects. Here, we show that inaccuracies in the modeling of $C_L^{\kappa\kappa}$ can yield surprisingly large biases on cosmological parameters inferred from the primary CMB power spectra measured by the upcoming Simons Observatory and CMB-S4 experiments. For CMB-S4, the biases can be as large as $1.6\sigma$ on the Hubble constant $H_0$ in a fit to $\Lambda$CDM and $1.2\sigma$ on $N_{\rm eff}$ in a fit to $\Lambda$CDM+$N_{\rm eff}$. We show that these biases can be mitigated by explicitly discarding all $TT$ data at $\ell>3000$ or by marginalizing over parameters describing baryonic feedback processes, both at the cost of slightly larger error bars. We also discuss an alternative, data-driven mitigation strategy based on delensing the CMB $T$ and $E$-mode maps. Finally, we show that analyses of upcoming data will require Einstein-Boltzmann codes to be run with much higher numerical precision settings than is currently standard, so as to avoid similar -- or larger -- parameter biases due to inaccurate theoretical predictions.

The $\sim$4 km diameter main belt asteroid 6478 Gault has ejected dust intermittently since at least 2013. The character of the emission, including its episodic nature and the low speed of the ejected particles ($V \sim $ 0.15 m s$^{-1}$), is most consistent with mass loss from a body rotating near rotational breakup. Owing to dust contamination of the nucleus signal, this conclusion has not yet been confirmed. To test this idea, we have obtained new images of Gault in August 2020, in the absence of dust. Our photometry shows a lightcurve having a very small amplitude (maximum $\sim 0.05$ mag) and a periodicity of $ 2.55 \pm 0.10$ hours. The new observations are consistent with a model in which Gault is rotating near breakup, with centrifugal forces responsible for its episodic mass loss. Approximated as a strengthless (fluid) spherical body, the implied density is $\rho$ = 1700 kg m$^{-3}$. We use the Froude number $Fr$, defined here as the ratio between centrifugal force and gravitational force, as a way to investigate mass loss regimes in fast spinning asteroids and find that mass shedding starts at $Fr \sim 0.5$.

We analyze 5 epochs of NICER data of the black hole X-ray binary MAXI J1820+070 during the bright hard-to-soft state transition in its 2018 outburst with both reflection spectroscopy and Fourier-resolved timing analysis. We confirm the previous discovery of reverberation lags in the hard state, and find that the frequency range where the (soft) reverberation lag dominates decreases with the reverberation lag amplitude increasing during the transition, suggesting an increasing X-ray emitting region, possibly due to an expanding corona. By jointly fitting the lag-energy spectra in a number of broad frequency ranges with the reverberation model reltrans, we find the increase in reverberation lag is best described by an increase in the X-ray coronal height. This result, along with the finding that the corona contracts in the hard state, suggests a close relationship between spatial extent of the X-ray corona and the radio jet. We find the corona expansion (as probed by reverberation) precedes a radio flare by ~5 days, which may suggest that the hard-to-soft transition is marked by the corona expanding vertically and launching a jet knot that propagates along the jet stream at relativistic velocities.

We propose a framework where a phase transition associated with a gauge symmetry breaking that occurs (not far) above the electroweak scale sets a stage for baryogenesis similar to the electroweak baryogenesis in the Standard Model. A concrete realization utilizes the breaking of $SU(2)_R \times U(1)_X \rightarrow U(1)_Y$. New chiral fermions charged under the extended gauge symmetry have nonzero lepton numbers, which makes the $B-L$ symmetry anomalous. The new lepton sector contains a large flavor-dependent CP violation, similar to the Cabibbo-Kobayashi-Maskawa phase, without inducing sizable electric dipole moments of the Standard Model particles. A bubble wall dynamics associated with the first-order phase transition and $SU(2)_R$ sphaleron processes generate a lepton asymmetry, which is transferred into a baryon asymmetry via the ordinary electroweak sphaleron process. Unlike the Standard Model electroweak baryogenesis, the new phase transition can be of the strong first order and the new CP violation is not significantly suppressed by Yukawa couplings, so that the observed asymmetry can be produced. The model can be probed by collider searches for new particles and the observation of gravitational waves. One of the new leptons becomes a dark matter candidate. The model can be also embedded into a left-right symmetric theory to solve the strong CP problem.

The first detection of gravitational waves from the binary neutron star merger GW170817 by the LIGO-Virgo Collaboration has provided fundamental new insights into the astrophysical site for r-process nucleosynthesis and on the nature of dense neutron-star matter. The detected gravitational wave signal depends upon the tidal distortion of the neutron stars as they approach merger. We report on relativistic numerical simulations of the approach to binary merger in the conformally flat, quasi-circular orbit approximation. We show that this event serves as a calibration to the quasi-circular approximation and a confirmation of the validity of the conformally flat approximation to the three-metric. We then examine how the detected chirp depends upon the adopted equation of state. This establishes a new efficient means to constrain the nuclear equation of state in binary neutron star mergers.

A general equation of state is considered for analysing the possible behaviors for (Early) dark energy that alleviates the Hubble parameter tension problem. By departing from the possible evolution for the (Early) dark energy density and the corresponding dynamical equations, the equation of state is obtained, which allow us to analyze qualitatively the cosmological evolution and the dominance of each term in the equation of state along the cosmic expansion, which show some interesting consequences as the occurrence of (past) future singularities. Then, by considering two general models, their free parameters are fit with different sources of data, showing the goodness of the fits in comparison to more standard models. Results might be considered as a promising starting point to get a better understanding of the cosmological evolution as a whole.

The spectrum of energy density fluctuations, baryon asymmetry, and coherent large-scale magnetic fields are the three observables that provide crucial information on physics at very high energies. Inflation can only provide a mechanism to explain the density perturbations, and the origin of primordial magnetic fields and baryon asymmetry require physics beyond the standard models of cosmology and particle physics. In this work, we show that the mechanism that leads to primordial helical fields also leads to baryogenesis at the beginning of the radiation-dominated epoch. The model we consider here consists of mass dimension 6 operators that include Riemann coupling between gravity and electromagnetic field without extending the Standard Model of particle physics. We explicitly show that the generation of helical magnetic fields leads to baryogenesis. We further show that the model predicts the observed amount of baryon asymmetry of the Universe for a range of reheating temperatures consistent with the observations.

We analyze the Post-Newtonian orbit of stars around a deformed Kerr black hole. The deformation we consider is a class of disformal transformations of a non-trivial Kerr solution in scalar-tensor theory which are labeled via the disformal parameter $D$. We study different limits of the disformal parameter, and compare the trajectories of stars orbiting a black hole to the case of the Kerr spacetime in general relativity, up to 2PN order. Our findings show that for generic non-zero $D$, the no-hair theorem of general relativity is violated, in the sense that the black hole's quadrupole Q is not determined by its mass $M$ and angular momentum $J$ through the relation $Q=-J^2/M$. Limiting values of $D$ provide examples of simple and exact non-circular metric solutions, whereas in a particular limit, where $1+D$ is small but finite, we obtain a leading correction to the Schwarzschild precession due to disformality. In this case, the disformal parameter is constrained using the recent measurement of the pericenter precession of the star S2 by the GRAVITY collaboration.

We present novel realizations of Higgs inflation within Supergravity which are largely tied to the existence of a pole of order two in the kinetic term of the inflaton field. This pole arises due to the selected Kaehler potentials which parameterize the (SU(1,1)/U(1))^2 or SU(2,1)/(SU(2)xU(1)) manifolds with scalar curvatures R_{(11)^2}=-4/N or R_{21}=-3/N respectively. The associated superpotential includes, in addition to the Higgs superfields, a stabilizer superfield, respects the gauge and an R symmetries and contains the first allowed nonrenormalizable term. If the coefficient of this term is almost equal to that of the renormalizable terms within about 10^-5 and N=1, the inflationary observables can be done compatible with the present data and the scale M of gauge-symmetry breaking may assume its value within MSSM. Increasing M beyond this value, though, inflation may be attained with less tuning. Modifications to the Kaehler potentials associated with the manifolds above allow for inflation, realized with just renormalizable terms, resulting to higher tensor-to-scalar ratios as N approaches its maximum at N=40.

Oumuamua, the first known object of extrasolar origin seen to enter our Solar System, has multiple unusual characteristics that, taken together, are very difficult to explain with conventional astronomical entities like asteroids and comets. Consequently, it has been hypothesized that Oumuamua is an interstellar probe that was constructed by an alien civilization. We demonstrate that the accomplishments that can be achieved with large space telescopes/interferometers in the alien's planetary system will completely quench any motivation for construction and launch of an Oumuamua-like probe. The absence of any such motivation proves that Oumuamua is not an alien creation.

"Warp drive" spacetimes and wormhole geometries are useful as "gedanken-experiments" that force us to confront the foundations of general relativity, and among other issues, to precisely formulate the notion of "superluminal" travel and communication. Here we will consider the basic definition and properties of warp drive spacetimes. In particular, we will discuss the violation of the energy conditions associated with these spacetimes, as well as some other interesting properties such as the appearance of horizons for the superluminal case, and the possibility of using a warp drive to create closed timelike curves. Furthermore, due to the horizon problem, an observer in a spaceship cannot create nor control on demand a warp bubble. To contour this difficulty, we discuss a metric introduced by Krasnikov, which also possesses the interesting property in that the time for a round trip, as measured by clocks at the starting point, can be made arbitrarily short.