Papers for Friday, Jul 07 2023

In the standard cosmological model of galaxy evolution, mergers and interactions play a fundamental role in shaping galaxies. Galaxies that are currently isolated are thus interesting, allowing us to distinguish between internal and external processes affecting the galactic structure. However, current observational limits may obscure crucial information in the low-mass or low-brightness regime. We use optical imaging of a subsample of the AMIGA catalogue of isolated galaxies to explore the impact of different factors on the structure of these galaxies. We study the type of disc break as a function of the degree of isolation and the presence of interaction indicators like tidal streams or plumes only detectable in the low surface brightness regime. We present deep optical imaging of a sample of 25 isolated galaxies. Through careful data processing and analysis techniques, the surface brightness limits achieved are comparable to those to be obtained on the 10-year LSST coadds. The extreme depth of our imaging allows us to study the interaction signatures of 20 galaxies, given that the presence of Galactic cirrus is a strong limiting factor in the characterisation of interactions for the remaining 5 of them. We detect previously unreported interaction features in 8 (40%) galaxies in our sample. We identify 9 galaxies (36%) showing an exponential disc (Type I), 14 galaxies (56%) with down-bending (Type II) profile and only 2 galaxies (8%) with up-bending (Type III) profiles. Isolated galaxies have considerably more purely exponential discs and fewer up-bending surface brightness profiles than field or cluster galaxies. We suggest that major mergers produce up-bending profiles while a threshold in star formation probably forms down-bending profiles. Unperturbed galaxies, evolving slowly with a low star formation rate could cause the high rate of Type I discs in isolated galaxies observed.

We present an end-to-end description of the formation process of globular clusters (GCs) which combines a treatment for their formation and dynamical evolution within galaxy haloes with a state-of-the-art semi-analytic simulation of galaxy formation. Our approach allows us to obtain exquisite statistics to study the effect of the environment and assembly history of galaxies, while still allowing a very efficient exploration of the parameter space of star cluster physics. Our reference model, including both efficient cluster disruption during galaxy mergers and a model for the dynamical friction of GCs within the galactic potential, accurately reproduces the observed correlation between the total mass in GCs and the parent halo mass. A deviation from linearity is predicted at low halo masses, which is driven by a strong dependence on morphological type: bulge-dominated galaxies tend to host larger masses of GCs than their later-type counterparts of similar stellar mass. While the significance of the difference might be affected by resolution at the lowest halo masses considered, this is a robust prediction of our model and represents a natural consequence of the assumption that cluster migration from the disk to the halo is triggered by galaxy mergers. Our model requires an environmental dependence of GC radii to reproduce the observed mass distribution of GCs in our Galaxy at the low-mass end. At GC masses $>10^6\,{\rm M}_\odot$, our model predicts fewer GCs than observed due to an overly aggressive treatment of dynamical friction. The metallicity distribution measured for Galactic GCs is well reproduced by our model, even though it predicts systematically younger GCs than observed. We argue that this adds further evidence for an anomalously early formation of the stars in our Galaxy.

We present, for the first time, dark matter halo (DMH) mass measurement of quasars at $z\sim6$ based on a clustering analysis of 107 quasars. Spectroscopically identified quasars are homogeneously extracted from the HSC-SSP wide layer over $891\,\mathrm{deg^2}$. We evaluate the clustering strength by three different auto-correlation functions: projected correlation function, angular correlation function, and redshift-space correlation function. The DMH mass of quasars at $z\sim6$ is evaluated as $5.0_{-4.0}^{+7.4}\times10^{12}\,h^{-1}M_\odot$ with the bias parameter $b=20.8\pm8.7$ by the projected correlation function. The other two estimators agree with these values, though each uncertainty is large. The DMH mass of quasars is found to be nearly constant $\sim10^{12.5}\,h^{-1}M_\odot$ throughout cosmic time, suggesting that there is a characteristic DMH mass where quasars are always activated. As a result, quasars appear in the most massive halos at $z \sim 6$, but in less extreme halos thereafter. The DMH mass does not appear to exceed the upper limit of $10^{13}\,h^{-1}M_\odot$, which suggests that most quasars reside in DMHs with $M_\mathrm{halo}<10^{13}\,h^{-1}M_\odot$ across most of the cosmic time. Our results supporting a significant increasing bias with redshift are consistent with the bias evolution model with inefficient AGN feedback at $z\sim6$. The duty cycle ($f_\mathrm{duty}$) is estimated as $0.019\pm0.008$ by assuming that DMHs in some mass interval can host a quasar. The average stellar mass is evaluated from stellar-to-halo mass ratio as $M_*=6.5_{-5.2}^{+9.6}\times10^{10}\,h^{-1}M_\odot$, which is found to be consistent with [C II] observational results.

Next-generation ultra-high-energy (UHE) neutrino telescopes, presently under planning, will have the potential to probe the decay of heavy dark matter (DM) into UHE neutrinos, with energies in excess of $10^7$~GeV. Yet, this potential may be deteriorated by the presence of an unknown background of UHE neutrinos, cosmogenic or from astrophysical sources, not of DM origin and seemingly large enough to obscure the DM signature. We show that leveraging the angular and energy distributions of detected events safeguards future searches for DM decay against such backgrounds. We focus on the radio-detection of UHE neutrinos in the planned IceCube-Gen2 neutrino telescope, which we model in state-of-the-art detail. We report promising prospects for the discovery potential of DM decay into UHE neutrinos, the measurement of DM mass and lifetime, and limits on the DM lifetime, despite the presence of a large background, without prior knowledge of its size and shape.

We perform a comprehensive search for optical precursor emission at the position of SN~2023ixf using data from the DLT40, ZTF and ATLAS surveys. By comparing the current data set with precursor outburst hydrodynamical model light curves, we find that the probability of a significant outburst within five years of explosion is low, and the circumstellar material (CSM) ejected during any possible precursor outburst is likely smaller than $\sim$0.015\msun. By comparing to a set of toy models, we find that, if there was a precursor outburst, the duration must have been shorter than $\sim$100 days for a typical brightness of $M_{r}\simeq-9$ mag or shorter than 200 days for $M_{r}\simeq-8$ mag; brighter, longer outbursts would have been discovered. Precursor activity like that observed in the normal type II SN~2020tlf ($M_{r}\simeq-11.5$) can be excluded in SN~2023ixf. If the dense CSM inferred by early flash spectroscopy and other studies is related to one or more precursor outbursts, then our observations indicate that any such outburst would have to be faint and only last for days to months, or it occurred more than five years prior to the explosion. Alternatively, any dense, confined CSM may not be due to eruptive mass loss from a single red supergiant (RSG) progenitor. Taken together, the results of SN~2023ixf and SN~2020tlf indicate that there may be more than one physical mechanism behind the dense CSM inferred around some normal type II SNe.

We present the first results of the Chandra Cool Targets (CCT) survey of the Second Bologna Catalog (B2CAT) of powerful radio sources, aimed at investigating the extended X-ray emission surrounding these sources. For the first 33 sources observed in the B2CAT CCT survey, we performed both imaging and spectral X-ray analysis, producing multi-band Chandra images, and compared them with radio observations. To evaluate the presence of extended emission in the X-rays, we extracted surface flux profiles comparing them with simulated ACIS Point Spread Functions. We detected X-ray nuclear emission for 28 sources. In addition, we detected 8 regions of increased X-ray flux originating from radio hot-spots or jet knots, and a region of decreased flux, possibly associated with an X-ray cavity. We performed X-ray spectral analysis for 15 nuclei and found intrinsic absorption significantly larger than the Galactic values in four of them. We detected significant extended X-ray emission in five sources, and fitted their spectra with thermal models with gas temperatures $\sim 2 \text{ keV}$. In the case of B2.1 0742+31, the surrounding hot gas is compatible with the ICM of low luminosity clusters of galaxies, while the X-ray diffuse emission surrounding the highly disturbed WAT B2.3 2254+35 features a luminosity similar to those of relatively bright galaxy groups, although its temperature is similar to those of low luminosity galaxy clusters. These results highlight the power of the low-frequency radio selection, combined with short Chandra snapshot observations, to investigate the properties of the X-ray emission from radio sources.

This paper reports the ULTRACAM discovery of dipolar surface spots in two cool magnetic white dwarfs with Balmer emission lines, while a third system exhibits a single spot, similar to the prototype GD 356. The light curves are modeled with simple, circular, isothermal dark spots, yielding relatively large regions with minimum angular radii of 20 deg. For those stars with two light curve minima, the dual spots are likely observed at high inclination (or colatitude), however, identical and antipodal spots cannot simultaneously reproduce both the distinct minima depths and the phases of the light curve maxima. The amplitudes of the multi-band photometric variability reported here are all several times larger than that observed in the prototype GD 356; nevertheless, all DAHe stars with available data appear to have light curve amplitudes that increase toward the blue in correlated ratios. This behavior is consistent with cool spots that produce higher contrasts at shorter wavelengths, with remarkably similar spectral properties given the diversity of magnetic field strengths and rotation rates. These findings support the interpretation that some magnetic white dwarfs generate intrinsic chromospheres as they cool, and that no external source is responsible for the observed temperature inversion. Spectroscopic time-series data for DAHe stars is paramount for further characterization, where it is important to obtain well-sampled data, and consider wavelength shifts, equivalent widths, and spectropolarimetry.

It has long been known that galaxy shapes align coherently with the large-scale density field. Characterizing this effect is essential to interpreting measurements of weak gravitational lensing, the deflection of light from distant galaxies by matter overdensities along the line of sight, as it produces coherent galaxy alignments that we wish to interpret in terms of a cosmological model. Existing direct measurements of intrinsic alignments using galaxy samples with high-quality shape and redshift measurements typically use well-understood but sub-optimal projected estimators, which do not make good use of the information in the data when comparing those estimators to theoretical models. We demonstrate a more optimal estimator, based on a multipole expansion of the correlation functions or power spectra, for direct measurements of galaxy intrinsic alignments. We show that even using the lowest order multipole alone increases the significance of inferred model parameters using simulated and real data, without any additional modeling complexity. We apply this estimator to measurements of intrinsic alignments in the Sloan Digital Sky survey, demonstrating consistent results with a factor of $\sim$2 greater precision in parameter fits to intrinsic alignments models. This result is functionally equivalent to quadrupling the survey area, but without the attendant costs -- thereby demonstrating the value in using this new estimator in current and future intrinsic alignments measurements using spectroscopic galaxy samples.

The gravitational wave (GW) signals from the Galactic population of cataclysmic variables (CVs) have yet to be carefully assessed. Here we estimate these signals and evaluate their significance for LISA. First, we find that at least three known systems are expected to produce strong enough signals to be individually resolved within the first four years of LISA's operation. Second, CVs will contribute significantly to the LISA Galactic binary background, limiting the mission's sensitivity in the relevant frequency band. Third, we predict a spike in the unresolved GW background at a frequency corresponding to the CV minimum orbital period. This excess noise may impact the detection of other systems near this characteristic frequency. Fourth, we note that the excess noise spike amplitude and location associated with $P_{\rm{min}}\sim80~\mathrm{min}$ can be used to measure the CV space density and period bounce location with complementary and simple GW biases compared to the biases and selection effects plaguing samples selected from electromagnetic signals. Our results highlight the need to explicitly include the Galactic CV population in the LISA mission planning, both as individual GW sources and generators of background noise, as well as the exciting prospect of characterising the CV population through their GW emission.

The accretion disk formed after a neutron star merger is an important contributor to the total ejecta from the merger, and hence to the kilonova and the $r$-process yields of each event. Axisymmetric viscous hydrodynamic simulations of these disks can capture thermal mass ejection due to neutrino absorption and in the advective phase -- after neutrino cooling has subsided -- and are thus likely to provide a lower-limit to the total disk ejecta relative to MHD evolution. Here we present a comparison between two viscous hydrodynamic codes that have been used extensively on this problem over the past decade: ALCAR and FLASH. We choose a representative setup with a black hole at the center, and vary the treatment of viscosity and neutrino transport. We find good overall agreement ($\sim 10\%$ level) in most quantities. The average outflow velocity is sensitive to the treatment of the nuclear binding energy of heavy nuclei, showing a larger variation than other quantities. We post-process trajectories from both codes with the same nuclear network, and explore the effects of code differences on nucleosynthesis yields, heating rates, and kilonova light curves. For the latter, we also assess the effect of varying the number of tracer particles in reconstructing the spatial abundance distribution for kilonova light curve production.

Recent low-frequency radio observations suggest that some nearby M dwarfs could be interacting magnetically with undetected close-in planets, powering the emission via the electron cyclotron maser (ECM) instability. Confirmation of such a scenario could reveal the presence of close-in planets around M dwarfs, which are typically difficult to detect via other methods. ECM emission is beamed, and is generally only visible for brief windows depending on the underlying system geometry. Due to this, detection may be favoured at certain orbital phases, or from systems with specific geometric configurations. In this work, we develop a geometric model to explore these two ideas. Our model produces the visibility of the induced emission as a function of time, based on a set of key parameters that characterise magnetic star-planet interactions. Utilising our model, we find that the orbital phases where emission appears are highly dependent on the underlying parameters, and does not generally appear at the quadrature points in the orbit as is seen for the Jupiter-Io interaction. Then using non-informative priors on the system geometry, we show that untargeted radio surveys are biased towards detecting emission from systems with planets in near face-on orbits. While transiting exoplanets are still likely to be detectable, they are less likely to be seen than those in near face-on orbits. Our forward model serves to be a powerful tool for both interpreting and appropriately scheduling radio observations of exoplanetary systems, as well as inverting the system geometry from observations.

We present the discovery and extensive follow-up of a remarkable fast-evolving optical transient, AT2022aedm, detected by the Asteroid Terrestrial impact Last Alert Survey (ATLAS). AT2022aedm exhibited a rise time of $9\pm1$ days in the ATLAS $o$-band, reaching a luminous peak with $M_g\approx-22$ mag. It faded by 2 magnitudes in $g$-band during the next 15 days. These timescales are consistent with other rapidly evolving transients, though the luminosity is extreme. Most surprisingly, the host galaxy is a massive elliptical with negligible current star formation. X-ray and radio observations rule out a relativistic AT2018cow-like explosion. A spectrum in the first few days after explosion showed short-lived He II emission resembling young core-collapse supernovae, but obvious broad supernova features never developed; later spectra showed only a fast-cooling continuum and narrow, blue-shifted absorption lines, possibly arising in a wind with $v\approx2700$ km s$^{-1}$. We identify two further transients in the literature (Dougie in particular, as well as AT2020bot) that share similarities in their luminosities, timescales, colour evolution and largely featureless spectra, and propose that these may constitute a new class of transients: luminous fast-coolers (LFCs). All three events occurred in passive galaxies at offsets of $\sim4-10$ kpc from the nucleus, posing a challenge for progenitor models involving massive stars or massive black holes. The light curves and spectra appear to be consistent with shock breakout emission, though usually this mechanism is associated with core-collapse supernovae. The encounter of a star with a stellar mass black hole may provide a promising alternative explanation.

The Kepler satellite has revealed a gap between sub-Neptunes and super-Earths that atmospheric escape models had predicted as an evaporation valley. We seek to contrast results from a simple XUV-driven energy-limited (ELIM) escape model against those from a direct hydrodynamic (HYDRO) model. Besides XUV-driven escape, the latter also includes the boil-off regime. We couple the two models to an internal structure model and follow the planets' temporal evolution over Gyr. To see the population-wide imprint of the two models, we first employ a rectangular grid in initial conditions. We then study the slope of the valley also for initial conditions derived from the Kepler planets. For the rectangular grid, we find that the power-law slope of the valley with respect to orbital period is -0.18 and -0.11 in the ELIM and HYDRO model, respectively. For the initial conditions derived from the Kepler planets, the results are similar (-0.16 and -0.10). While the slope found with the ELIM model is steeper than observed, the one of the HYDRO model is in excellent agreement with observations. The reason for the shallower slope is caused by the two regimes in which the ELIM model fails: First, puffy planets at low stellar irradiation. For them, boil-off dominates mass loss. However, boil-off is absent in the ELIM model, thus it underestimates escape relative to HYDRO. Second, massive compact planets at high XUV irradiation. For them, the ELIM approximation overestimates escape relative to the HYDRO case because of cooling by thermal conduction, neglected in the ELIM model. The two effects act together in concert to yield in the HYDRO model a shallower slope of the valley that agrees very well with observations. We conclude that an escape model that includes boil-off and a more realistic treatment of cooling mechanisms can reproduce one of the most important constraints, the valley slope.

We analyse the joint distribution of dust attenuation and projected axis ratios, together with galaxy size and surface brightness profile information, to infer lessons on the dust content and star/dust geometry within star-forming galaxies at 0 < z <2.5. To do so, we make use of large observational datasets from KiDS+VIKING+HSC-SSP and extend the analysis out to redshift z = 2.5 using the HST surveys CANDELS and 3D-DASH. We construct suites of SKIRT radiative transfer models for idealized galaxies observed under random viewing angles with the aim of reproducing the aforementioned distributions, including the level and inclination dependence of dust attenuation. We find that attenuation-based dust mass estimates are at odds with constraints from far-infrared observations, especially at higher redshifts, when assuming smooth star and dust geometries of equal extent. We demonstrate that UV-to-near-IR and far-infrared constraints can be reconciled by invoking clumpier dust geometries for galaxies at higher redshifts and/or very compact dust cores. We discuss implications for the significant wavelength- and redshift-dependent differences between half-light and half-mass radii that result from spatially varying dust columns within -- especially massive -- star-forming galaxies.

We present stellar parameters and chemical abundances of 30 elements for five stars located at large radii (3.5-10.7 times the half-light radius) in the Sextans dwarf spheroidal galaxy. We selected these stars using proper motions, radial velocities, and metallicities, and we confirm them as metal-poor members of Sextans with -3.34 <= [Fe/H] <= -2.64 using high-resolution optical spectra collected with the Magellan Inamori Kyocera Echelle spectrograph. Four of the five stars exhibit normal abundances of C (-0.34 <= [C/Fe] <= +0.36), mild enhancement of the alpha elements Mg, Si, Ca, and Ti ([alpha/Fe] = +0.12 +/- 0.03), and unremarkable abundances of Na, Al, K, Sc, V, Cr, Mn, Co, Ni, and Zn. We identify three chemical signatures previously unknown among stars in Sextans. One star exhibits large overabundances ([X/Fe] > +1.2) of C, N, O, Na, Mg, Si, and K, and large deficiencies of heavy elements ([Sr/Fe] = -2.37 +/- 0.25, [Ba/Fe] = -1.45 +/- 0.20, [Eu/Fe] < +0.05), establishing it as a member of the class of carbon-enhanced metal-poor stars with no enhancement of neutron-capture elements. Three stars exhibit moderate enhancements of Eu (+0.17 <= [Eu/Fe] <= +0.70), and the abundance ratios among 12 neutron-capture elements are indicative of r-process nucleosynthesis. Another star is highly enhanced in Sr relative to heavier elements ([Sr/Ba] = +1.21 +/- 0.25). These chemical signatures can all be attributed to massive, low-metallicity stars or their end states. Our results, the first for stars at large radius in Sextans, demonstrate that these stars were formed in chemically inhomogeneous regions, such as those found in ultra-faint dwarf galaxies.

In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code, SimSpin, is an open-source package written in R, but also with an API interface such that the code can be interacted with in any coding language. Documentation and individual examples can be found at the open-source website connected to the online repository. SimSpin is already being utilised by international IFS collaborations, including SAMI and MAGPI, for generating comparable data sets from a diverse suite of cosmological hydrodynamical simulations.

Continuum reverberation mapping probes the sizescale of the optical continuum-emitting region in active galactic nuclei (AGN). Through 3 years of multiwavelength photometric monitoring in the optical with robotic observatories, we perform continuum reverberation mapping on Mrk~876. All wavebands show large amplitude variability and are well correlated. Slow variations in the light curves broaden the cross-correlation function (CCF) significantly, requiring detrending in order to robustly recover interband lags. We measure consistent interband lags using three techniques (CCF, JAVELIN, PyROA), with a lag of around 13~days from $u$ to $z$. These lags are longer than the expected radius of 12~days for the self-gravitating radius of the disk. The lags increase with wavelength roughly following $\lambda^{4/3}$, as would be expected from thin disk theory, but the lag normalization is approximately a factor of 3 longer than expected, as has also been observed in other AGN. The lag in the $i$ band shows an excess which we attribute to variable H$\alpha$ broad-line emission. A flux-flux analysis shows a variable spectrum that follows $f_\nu \propto \lambda^{-1/3}$ as expected for a disk, and an excess in the $i$ band that also points to strong variable H$\alpha$ emission in that band.

In the past years, the results obtained by the WISSH quasar project provided a novel general picture on the distinctive multi-band properties of hyper-luminous ($L_{bol}>10^{47}$ erg/s) quasars at high redshift (z$\sim$2-4), unveiling interesting relations among active galactic nuclei, winds and interstellar medium, in these powerful sources at cosmic noon. Since 2022, we are performing a systematic and statistically-significant VLA study of the radio properties of WISSH. We carried out high-resolution VLA observations aiming at: 1) identifying young radio source from the broad-band spectral shape of these objects; 2) sample an unexplored high redshift/high luminosity regime, tracking possible evolutionary effects on the radio-loud/radio-quiet dichotomy; 3) quantifying orientation effects on the observed winds/outflows properties.

CO is one of the most abundant molecules in protoplanetary disks, and optically thin emission from its isotopologues has been detected in many of them. However, several past works have argued that reproducing the relatively low emission of CO isotopologues requires a very low disk mass or significant CO depletion. Here, we present a Python code, DiskMINT, which includes gas density and temperature structures that are both consistent with the thermal pressure gradient, isotope-selective chemistry, and conversion of CO into $\mathrm{CO_2}$ ice on grain-surfaces. The code generates a self-consistent disk structure, where the gas disk distribution is obtained from a Spectral Energy Distribution (SED)-derived dust disk structure with multiple grain sizes. We use DiskMINT to study the disk of RU~Lup, a high-accreting star whose disk was previously inferred to have a gas mass of only $\sim 1.5\times10^{-3}\,M_\odot$ and gas-to-dust mass ratio of $\sim 4$. Our best-fit model to the long-wavelength continuum emission can explain the total $\mathrm{C^{18}O}$ luminosity as well as the $\mathrm{C^{18}O}$ velocity and radial intensity profiles, and obtains a gas mass of $\sim 1.2\times10^{-2}\,M_\odot$, an order of magnitude higher than previous results. A disk model with parametric Gaussian vertical distribution that better matches the IR-SED can also explain the observables above with a similarly high gas mass $\sim 2.1\times10^{-2}\,M_\odot$. We confirm the conclusions of Ruaud et al. (2022) that optically thin $\mathrm{C^{18}O}$ rotational lines provide reasonable estimates of the disk mass and can therefore be used as gas disk tracers.

We motivate the cosmological science case of measuring Type Ia supernovae with the Nancy Grace Roman Space Telescope as part of the High Latitude Time Domain Survey. We discuss previously stated requirements for the science, and a baseline survey strategy. We discuss the various areas that must still be optimized and point to the other white papers that consider these topics in detail. Overall, the baseline case should enable an exquisite measurement of dark energy using SNe Ia from z=0.1 to z>2, and further optimization should only strengthen this once-in-a-generation experiment.

This research addresses the influence of Primordial Black Holes (PBHs) on cosmic reionization, using a robust semi-analytical model. This model encapsulates cosmological theory, PBH physics, and radiative transfer, with a lognormal PBH mass function as the fulcrum, with the mean PBH mass set at 30 solar masses and the fraction of dark matter in PBHs ($f_{PBH}$) varied from 1.26 to 10. We also adopt a standard thin-disk accretion model and a one-dimensional radiative transfer code for a comprehensive depiction of cosmic reionization under the influence of PBHs. Our major findings reveal an inverse relationship between PBH mass and ionizing efficiency. For instance, when the mean PBH mass is 10 solar masses, and $f_{PBH}$ is 0.1, the resulting ionizing efficiency is approximately 0.25, which drops to 0.15 with a mean PBH mass of 50 solar masses. However, with $f_{PBH}$ nearing unity, even large PBHs (mean mass of 100 solar masses) can achieve an ionizing efficiency close to 0.2. The ionization history indicates that for a mean PBH mass of 10 solar masses and $f_{PBH}$ at 0.01, rapid ionization commences at z=12.5+0.5/z=12.5-0.5 and concludes at z=9+0.5/z=9-0.5. Conversely, with a larger PBH mass of 100 solar masses and $f_{PBH}$ at 0.01, reionization initiates later, at z=10.8+0.4/z=10.8-0.4. Our two-point correlation function analysis unveils the formation of ionized bubbles with an initial size of approximately 10 Mpc at z=7.5+0.5/z=7.5-0.5 for an $f_{PBH}$ of 0.01, expanding to 20 Mpc by the end of reionization. With $f_{PBH}$ of 0.1, larger bubbles form, starting at approximately 15 Mpc and reaching 30 Mpc by reionization's end. In conclusion, PBHs significantly influence cosmic reionization, with PBH mass and their contribution to dark matter modulating this impact.

Binary stars evolve into chemically-peculiar objects and are a major driver of the Galactic enrichment of heavy elements. During their evolution they undergo interactions, including tides, that circularize orbits and synchronize stellar spins, impacting both individual systems and stellar populations. Using Zahn's tidal theory and MESA main-sequence model grids, we derive the governing parameters $\lambda_{lm}$ and $E_2$, and implement them in the new MINT library of the stellar population code BINARY_C. Our MINT equilibrium tides are 2 to 5 times more efficient than the ubiquitous BSE prescriptions while the radiative-tide efficiency drops sharply with increasing age. We also implement precise initial distributions based on bias-corrected observations. We assess the impact of tides and initial orbital-parameter distributions on circularization and synchronization in eight open clusters, comparing synthetic populations and observations through a bootstrapping method. We find that changing the tidal prescription yields no statistically-significant improvement as both calculations typically lie within 0.5$\sigma$. The initial distribution, especially the primordial concentration of systems at $\log_{10}(P/{\rm d}) \approx 0.8, e\approx 0.05$ dominates the statistics even when artificially increasing tidal strength. This confirms the inefficiency of tides on the main sequence and shows that constraining tidal-efficiency parameters using the $e-\log_{10}(P/{\rm d})$ distribution alone is difficult or impossible. Orbital synchronization carries a more striking age-dependent signature of tidal interactions. In M35 we find twice as many synchronized rotators in our MINT calculation as with BSE. This measure of tidal efficiency is verifiable with combined measurements of orbital parameters and stellar spins.

In this paper we revisit metal-enriched pair instability supernovae (PISNe) models which undergo chemically homogeneous evolution (CHE). By calculating multiple models, we intend to clarify mass ranges for the PISNe, $^{56}$Ni masses from the PISNe, and mass loss histories of CHE-PISNe models for metallicities consistent with the Small Magellanic Cloud (SMC) and with the Large Magellanic Cloud (LMC). We show that for an initial velocity of $v_{\rm i}/ v_{\rm k}$ = 0.1, these models undergo CHE and He-rich (Type Ib) PISNe occur in a lower mass range ($M_{\rm i} \sim 110-170 M_\odot$) than for more slowly rotating models. Interestingly, bright PISNe which have $^{56}$Ni masses larger than 10 $M_\odot$ occur in a relatively small mass range, $M_{\rm i} \sim 140-170 M_\odot$. Another notable characteristic of CHE-PISNe is the large late time mass loss rates; consequently, CSM interaction may be observable in their light curves. We also show some examples of O-rich (Type Ic) CHE-PISNe produced by $v_{\rm i}/ v_{\rm k}$ = 0.2 models. We expect these models to exhibit interaction with O-rich CSM, behavior which is consistent with the observed properties of the recently discovered PISN candidate, SN2018ibb. Finally, we present a collapsing $v_{\rm i}/ v_{\rm k}$ = 0.2 model which has sufficient angular momentum to be regarded as a candidate for a Super-Kilonova.

We present an ACA search for [CI] emission at 492GHz toward large T Tauri disks (gas radii $\gtrsim 200$au) in the $\sim 1-3$Myr-old Lupus star-forming region. Combined with ALMA 12-m archival data for IM Lup, we report [CI] detections in 6 out of 10 sources, thus doubling the known detections toward T Tauri disks. We also identify four Keplerian double-peaked profiles and demonstrate that [CI] fluxes correlate with $^{13}$CO, C$^{18}$O, and $^{12}$CO(2-1) fluxes, as well as with the gas disk outer radius measured from the latter transition. These findings are in line with the expectation that atomic carbon traces the disk surface. In addition, we compare the carbon and CO line luminosities of the Lupus and literature sample with [CI] detections with predictions from the self-consistent disk thermo-chemical models of Ruaud et al. (2022). These models adopt ISM carbon and oxygen elemental abundances as input parameters. With the exception of the disk around Sz 98, we find that these models reproduce all available line luminosities and upper limits with gas masses comparable to or higher than the minimum mass solar nebula and gas-to-dust mass ratios $\geq 10$. Thus, we conclude that the majority of large Myr-old disks conform to the simple expectation that they are not significantly depleted in gas, CO, or carbon.

We use the EAGLE (Evolution and Assembly of GaLaxies and their Environments) and IllustrisTNG (The Next Generation) cosmological simulations to investigate the properties of the baryonic specific angular momentum (j), baryonic mass (M) and atomic gas fraction ($f_{\rm{atm}}$) plane for nearby galaxies. We find an excellent agreement between EAGLE and TNG, with both also matching quite well the results obtained with xGASS (eXtended GALEX Arecibo SDSS Survey) for gas fractions greater than 0.01. This implies that the disagreements previously identified between xGASS and predictions from simple analytical disc stability arguments also holds true for EAGLE and TNG. For lower gas fraction (the regime currently unconstrained by observations), both simulations deviate from the plane but still maintain good agreement with each other. Despite the challenges posed by resolution limits at low gas fractions, our findings suggest a potential disconnect between angular momentum and gas fraction in the gas-poor regime, implying that not all gas-poor galaxies have low specific angular momentum.

Pulsar timing array (PTA) experiments are expected to detect nano-Hertz gravitational waves (GWs) emitted from individual inspiralling supermassive binary black holes (SMBBHs). The GW signals from a small fraction of these SMBBHs may be diffractively lensed by intervening galaxies. In this paper, we investigate the diffractive lensing effects on the continuous GW signals from the lensed SMBBHs and estimate the detectable number of such signals by PTAs, such as the Chinese PTA (CPTA) and the Square Kilometer Array (SKA) PTA. We find that the amplitude of the lensed GW signals may be only amplified by a factor of $\sim 1.01-1.14$ ($16\%-84\%$ range) and the phase of the signals may shift somewhat due to the lensing, significantly different from those strongly lensed high frequency GW signals from compact binary mergers in the geometric optics. We estimate that $\sim 0.01\%$ of all detected nano-Hertz GW signals from individual SMBBHs by future PTA experiments are lensed by foreground galaxies (i.e., up to $\sim 106$ for CPTA and up to $\sim 289$ for SKA-PTA). However, the lensed nano-Hertz GW signals are difficult to be distinguished from those without lensing by the PTA observations only. We further discuss the possibility about the identification of the lensed nano-Hertz GW signals from SMBBHs via the electromagnetic detection of their host galaxies or active galactic nuclei.

Future ground-based gravitational wave (GW) detectors, i.e., Einstein telescope (ET) and Cosmic Explorer (CE), are expected to detect a significant number of lensed binary neutron star (BNS) mergers, which may provide a unique tool to probe cosmology. In this paper, we investigate the detectability of the optical/infrared electromagnetic (EM) counterparts (kilonovae/afterglows) from these lensed BNS mergers by future GW detectors and EM telescopes using simple kilonova, afterglow, and lens models. ET and CE are expected to detect $\sim5.32^{+26.1}_{-5.10}$ and $67.3^{+332}_{-64.7}$ lensed BNS mergers per year. We find that the EM counterparts associated with all these mergers will be detectable by an all sky-survey in the H-band with the limiting magnitude $m_{\textrm{lim}}\gtrsim27$, while the detectable fraction is $\lesssim0.4\%$ in the g-/z-band if with $m_{\textrm{lim}}\lesssim24$. Generally it is more efficient to search the lensed EM counterparts by adopting the infrared bands than the optical/UV bands with the same $m_{\textrm{lim}}$. Future telescopes like Vera C. Rubin Observatory, China Space Station Telescope, and Euclid can hardly detect the EM counterparts of even one lensed BNS merger. Roman Space Telescope (RST) and James Webb Space Telescope (JWST) have the capability to detect about a few or more such events per year. Moreover, the time delays and separations between the lensed image pairs are typically in the ranges from minutes to months and from $0.1$ to $1$\,arcsec, suggesting that both the GW and EM images of most lensed BNS mergers can be well resolved by not only CE/ET in the time domain but also RST/JWST spatially.

High-velocity atomic clouds in the Galactic center have attracted significant attention due to their enigmatic formation process, which is potentially linked to the starburst or supermassive black hole activities in the region. Further, the discovery of high-velocity molecular clouds (HVMCs) presents a greater puzzle, because they are much denser and more massive. If the HVMCs were accelerated by the strong activities in the Galactic center, they are expected to be destroyed before they reach such a high velocity. To shed light on this phenomenon, we perform three-dimensional numerical simulations to investigate the origin and hydrodynamic evolution of HVMCs during a starburst in the Galactic center. We find that the presence of a magnetic field provides effective protection and acceleration to molecular clouds (MCs) within the galactic winds. Consequently, the MCs can attain latitudes of approximately 1 kpc with velocities around 200 km/s, consistent with the observed characteristics of HVMCs. The consistency of our findings across a wide parameter space supports the conclusion that HVMCs can indeed withstand the starburst environment in the Galactic center, providing valuable insights into their survival mechanisms.

The partcle-induced hadronic de-excitation of the Hoyle state in $^{12}$C induced by inelastic scattering in a hot and dense plasma can enhance the triple-alpha reaction rate. This prevents the production of heavy nuclei within the neutrino-driven winds of core-collapse supernovae and raises a question as to the contribution of proton-rich neutrino-driven winds as the origin of $p$-nuclei in the solar system abundances. Here we study $\nu p$-process nucleosynthesis in proton-rich neutrino-driven winds relevant to the production of $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$ by considering such particle-induced de-excitation. We show that the enhancement of the triple-alpha reaction rate induced by neutron inelastic scattering hardly affects the $\nu p$-process, while the proton scattering contributes to the nucleosynthesis in proton-rich neutrino-driven winds at low temperature. The associated enhanced triple-alpha reaction rate decreases the production of $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$ in a wind model of ordinary core-collapse supernovae. On the other hand, the abundances of these $p$-nuclei increase in an energetic hypernova wind model. Hence, we calculate the galactic chemical evolution of $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$ by taking account of both contributions from core-collapse supernovae and hypernovae. We show that the hypernova $\nu p$-process can enhance the calculated solar isotopic fractions of $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$ and make a significant impact on the GCE of $p$-nuclei regardless of the particle-induced Hoyle state de-excitation.

We present the first look at star formation histories of galaxy components using ProFuse, a new technique to model the 2D distribution of light across multiple wavelengths using simultaneous spectral and spatial fitting of purely imaging data. We present a number of methods to classify galaxies structurally/morphologically, showing the similarities and discrepancies between these schemes. Rather than identifying the best-performing scheme, we use the spread of classifications to quantify uncertainty in our results. We study the cosmic star formation history (CSFH), forensically derived using ProFuse with a sample of ~7,000 galaxies from the Galaxy And Mass Assembly (GAMA) survey. Remarkably, the forensic CSFH recovered via both our method (ProFuse) and traditional SED fitting (ProSpect) are not only exactly consistent with each other over the past 8 Gyr, but also with the in-situ CSFH measured using ProSpect. Furthermore, we separate the CSFH by contributions from spheroids, bulges and disks. While the vast majority (70%) of present-day star formation takes place in the disk population, we show that 50% of the stars that formed at cosmic noon (8-12 Gyr ago) now reside in spheroids, and present-day bulges are composed of stars that were primarily formed in the very early Universe, with half their stars already formed ~12 Gyr ago.

We observed active M dwarf star AD Leo for 146 hr in photometry by GWAC-F30 and also analyzed 528-hr photometric data of the star from TESS. A total of 9 and 70 flares are detected from GWAC-F30 and TESS, respectively. Flare durations, amplitudes and energies are calculated. The distributions of the three properties and FFDs are given. Within the same energy range of flares, the FFDs of AD Leo obtained in this research and the previous study are basically consistent, which suggests that the magnetic activity of this star has not significantly changed compared to that decades ago. Comparing with the average FFD of M-type stars, AD Leo's FFD is twice higher, indicating that its magnetic activity is more active than that of the average level of the M-type. Based on TESS light curve, AD Leo's rotation period is calculated as 2.21${+0.01 \choose -0.01}$ day , supporting the result given in previous research. During the decay phase of the most energetic flare from TESS, we identified QPPs and determined a 26.5-min oscillation period, which is currently the longest period for AD Leo, suggesting that long periodic physical process existed during flare of this star.

Fast Radio Bursts (FRBs) are mysterious bursts in the millisecond timescale at radio wavelengths. Currently, there is little understanding about the classification of repeating FRBs, based on difference in physics, which is of great importance in understanding their origin. Recent works from the literature focus on using specific parameters to classify FRBs to draw inferences on the possible physical mechanisms or properties of these FRB subtypes. In this study, we use publicly available 1652 repeating FRBs from FRB121102 detected with the Five-hundred-meter Aperture Spherical Telescope (FAST), and studied them with an unsupervised machine learning model. By fine-tuning the hyperparameters of the model, we found that there is an indication for four clusters from the bursts of FRB121102 instead of the two clusters ("Classical" and "Atypical") suggested in the literature. Wherein, the "Atypical" cluster can be further classified into three sub-clusters with distinct characteristics. Our findings show that the clustering result we obtained is more comprehensive not only because our study produced results which are consistent with those in the literature but also because our work uses more physical parameters to create these clusters. Overall, our methods and analyses produced a more holistic approach in clustering the repeating FRBs of FRB121102.

Multi-line molecular observations are an ideal tool for a systematic study of the physico-chemical processes in the Interstellar Medium (ISM), given the wide range of critical densities associated with different molecules and their transitions, and the dependencies of chemical reactions on the energy budget of the system. Recently high spatial resolution of typical shock tracers - SiO, HNCO, and CH3OH - have been studied in the potentially shocked regions in two nearby galaxies: NGC 1068 (an AGN-host galaxy) (Huang et al., Astron. Astrophys., 2022, 666, A102; Huang et al., in prep.) and NGC 253 (a starburst galaxy) (K.-Y. Huang et al., arXiv, 2023, preprint, arXiv:2303.12685, DOI: 10.48550/arXiv.2303.12685). This paper is dedicated to the comparative study of these two distinctively different galaxies, with the aim of determining the differences in their energetics and understanding large-scale shocks in dfferent types of galaxies.

Devasthal observatory, established over a time span of about 5 decades, is located in central Himalayan region of Devabhumi in Nainital district of Uttarakhand state, India. Operated and maintained by the Aryabhatta Research Institute of observational sciencES (ARIES), its location was selected after an extensive site survey. The first measurements of atmospheric seeing and extinctions at Devasthal were carried out during 1997 to 2001. Since 2010, three optical telescopes of 1.3 m, 3.6 m and 4 m apertures have been successfully installed at Devasthal. Optical and near infrared observations taken with these telescopes testify to the global competitiveness of Devasthal observatory for astronomical observations. The article chronicles the collaboration with the Tata Institute of Fundamental Research, beginning around 1996, for the purpose of establishing the observatory. A brief overview of the main science results obtained so far, using these facilities, is also presented.

We propose a model explaining the origin of transient/episodic jets in black-hole X-ray binaries, in which they are caused by transitions from a collimated, strongly magnetized, jet to a wide, un-collimated, outflow. The change occurs when the accretion flow leaves the magnetically-choked state due to an increase of the accretion rate at a constant magnetic flux. The formed powerful jet then detaches from its base, and propagates as a discrete ejection. The uncollimated outflow then produces a relativistic plasma that fills surrounding of the black hole, contributing to the formation of a low-density cavity. While the pressure in the cavity is in equilibrium with the surrounding interstellar medium (ISM), its inertia is orders of magnitude lower than that of the ISM. This implies that the plasma cannot efficiently decelerate the ejecta, explaining most of the observations. The modest deceleration within the cavities observed in some cases can be then due to the presence of clouds and/or filaments, forming a wide transition zone between the cavity and the ISM.

The cosmic evolution of the barred galaxy population provides key information about the secular evolution of galaxies and the settling of rotationally dominated discs. We study the bar fraction in the SMACSJ0723.37323 (SMACS0723) cluster of galaxies at z = 0.39 using the Early Release Observations obtained with the NIRCam instrument mounted on the JWST telescope. As already found in nearby galaxy samples, we find that the bar fraction distribution of SMACS0723 is a strong function of the galaxy stellar luminosity/mass. The analogy with local clusters, such as Virgo and Coma, reveals a similar distribution among the three clusters for low-mass galaxies (log(M_star/M_sun) \leq 9.5). The comparison with a sample of local galaxies in a field environment shows a remarkable lack of bars in this low-mass regime for the SMACS0723 cluster (and therefore in Virgo and Coma) with respect to the field. At high masses (log(M_star/M_sun) \geq 10.25), galaxies in SMACS0723 show a slightly lower bar fraction than those in Coma. Our results support a scenario where cluster environment affects the formation of bars in a mass-dependent way. At high masses, the mild increase in the bar fraction of local clusters (Coma) with respect to both SMACS0723 and local field galaxies suggests a weak effect of cluster environment possibly triggering bar formation. On the other hand, low-mass galaxies show the same bar fraction in the three clusters (different redshifts) and a significant drop with respect to field galaxies at z=0, therefore suggesting that: i) the bar fraction of low-mass galaxies in clusters is not evolving during the last 4~Gyr, and ii) bar formation is severely inhibited in low-mass galaxies living in clusters (Abridged).

High-energy neutrino and $\gamma$-ray emission has been observed from the Galactic plane, which may come from individual sources and/or diffuse cosmic rays. We evaluate the contribution of these two components through the multi-messenger connection between neutrinos and $\gamma$ rays in hadronic interactions. We derive maximum fluxes of neutrino emission from the Galactic plane using $\gamma$-ray catalogs, including 4FGL, HGPS, 3HWC, and 1LHAASO, and measurements of the Galactic diffuse emission by Tibet AS$\gamma$ and LHAASO. We find that depending on model templates, the diffuse emission is brighter than the sum of resolved sources when excluding promising leptonic sources such as pulsars, pulsar wind nebulae, and TeV halos. Our result indicates that the Galactic neutrino emission observed by the IceCube Collaboration may be dominated by the Galactic diffuse emission or unresolved $\gamma$-ray sources. Future observations of neutrino telescopes and air-shower $\gamma$-ray experiments in the Southern hemisphere are needed to accurately disentangle the source and diffuse emission of the Milky Way.

The current generation of air shower radio arrays has demonstrated that the atmospheric depth of the shower maximum Xmax can be reconstructed with high accuracy. These experiments are now contributing to mass composition studies in the energy range where a transition from galactic to extragalactic cosmic-ray sources is expected. However, we are still far away from an unambiguous interpretation of the data. Here we propose to use radio measurements to derive a new type of constraint on the mass composition, by reconstructing the shower length L. The low-frequency part of the Square Kilometer Array will have an extremely high antenna density of roughly 60.000 antennas within one square kilometer, and is the perfect site for high-resolution studies of air showers. In this contribution, we discuss the impact of being able to reconstruct L, and the unique contribution that SKA can make to cosmic-ray science.

The 21 cm signal from cosmic hydrogen is one of the most propitious probes of the early Universe. The detection of this signal would reveal key information about the first stars, the nature of dark matter, and early structure formation. We explore the impact of an emissive and reflective, or `hot', horizon on the recovery of this signal for global 21 cm experiments. It is demonstrated that using physically motivated foreground models to recover the sky-averaged 21 cm signal one must accurately describe the horizon around the radiometer. We show that not accounting for the horizon will lead to a signal recovery with residuals an order of magnitude larger than the injected signal, with a log Bayesian evidence of almost 1600 lower than when one does account for the horizon. It is shown that signal recovery is sensitive to incorrect values of soil temperature and reflection coefficient in describing the horizon, with even a 10% error in reflectance causing twofold increases in the RMSE of a given fit. We also show these parameters may be fitted using Bayesian inference to mitigate for these issues without overfitting and mischaracterising a non-detection. We further demonstrate that signal recovery is sensitive to errors in measurements of the horizon projection onto the sky, but fitting for soil temperature and reflection coefficients with priors that extend beyond physical expectation can resolve these problems. We show that using an expanded prior range can reliably recover the signal even when the height of the horizon is mismeasured by up to 20%, decreasing the RMSE from the model that does not perform this fitting by a factor of 9.

Our aim is to examine the size, kinematics, and geometry of the broad-line region (BLR) in the double-lensed quasar Q 0957+561 by analyzing the impact of microlensing on various rest-frame ultraviolet broad-emission lines (BELs). We explore the influence of intrinsic variability and microlensing on the C IV, C III], and Mg II emission lines through multiple spectroscopic observations taken between April 1999 and January 2017. By utilizing the line cores as a reference for no microlensing and correcting for the long time delay between the images, we estimate the sizes of the regions emitting the broad-line wings using a Bayesian approach. Our study of the microlensing amplitudes between the lensed images of the quasar Q 0957+561 reveals differing sizes of the regions emitting the three prominent BELs C IV, C III], and Mg II. The strength of the differential microlensing indicates that the high-ionization line C IV arises from a compact inner region of the BLR with a half-light radius of $R_{1/2} \gtrsim 16.0$ lt-days, which represents a lower limit on the overall size of the BLR and is comparable to the size of the region emitting the r-band continuum in this system. A somewhat larger size of $R_{1/2}\gtrsim 44$ lt-days is obtained for the semi-forbidden line C III]. Microlensing has a weak impact on the lower-ionization line Mg II, which is emitted from a region with a half-light radius of $R_{1/2} \gtrsim 50$ lt-days. These findings suggest that the BEL regions may have distinct geometries and kinematics, with the more extended ones being spherically symmetric, and the most compact ones being nonspherical, with motions likely confined to a plane.

We study a sample of 119 Herbig Ae/Be stars in the Galactic anti-center direction using the spectroscopic data from Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) survey program. Emission lines of hydrogen belonging to the Balmer and Paschen series, and metallic lines of species such as FeII, OI, CaII triplet are identified. A moderate correlation is observed between the emission strengths of H$\alpha$ and FeII 5169 \r{A}, suggesting a possible common emission region for FeII lines and one of the components of H$\alpha$. We explored a technique for the extinction correction of the HAeBe stars using diffuse interstellar bands present in the spectrum. We estimated the stellar parameters such as age and mass of these HAeBe stars, which are found to be in the range 0.1 -- 10 Myr and 1.5 -- 10 $M_{\odot}$, respectively. We found that the mass accretion rate of the HAeBe stars in the Galactic anti-center direction follows the relation $\dot{M}_{acc}$ $\propto$ $M_{*}^{3.12^{+0.21}_{-0.34}}$, which is similar to the relation derived for HAeBe stars in other regions of the Galaxy. The mass accretion rate of HAeBe stars is found to have a functional form of $\dot{M}_{acc} \propto t^{-1.1 \pm 0.2}$ with age, in agreement with previous studies.

The radio detection technique of cosmic ray air showers has gained renewed interest in the last two decades. While the radio experiments are very cost-effective to deploy, the Monte-Carlo simulations required to analyse the data are computationally expensive. Here we present a novel way to synthesise the radio emission from extensive air showers in simulations. It is a hybrid approach which uses a single microscopic Monte-Carlo simulation, called the origin shower, to generate the radio emission from a target shower with a different longitudinal evolution, primary particle type and energy. The method employs semi-analytical relations which only depend on the shower parameters to transform the radio signals in the simulated antennas. We apply this method to vertical air showers with energies ranging from $10^{17} \; \text{eV}$ to $10^{19} \; \text{eV}$ and compare the results with CoREAS simulations in two frequency bands, namely the broad [20, 500] MHz band and a more narrow one at [30, 80] MHz . We gauge the synthesis quality using the maximal amplitude and energy fluence contained in the signal. We observe that the quality depends primarily on the difference in \textXmax\ between the origin and target shower. After applying a linear bias correction, we find template synthesis to have an accuracy of better than 10% on the broad frequency range. On the restricted [30, 80] MHz range the accuracy is better than 5%. We surmise the template synthesis approach to be limited by intrinsic shower fluctuations which were not accounted for.

We analyse the rest-frame UV spectra of 2,531 high-redshift (3.5<z<4.0) quasars from the Sloan Digital Sky Survey DR16Q catalogue. In combination with previous work, we study the redshift evolution of the rest-frame UV line properties across the entire redshift range, 1.5<z<4.0. We improve the systemic redshift estimates at z>3.5 using a cross-correlation algorithm that employs high signal-to-noise template spectra spanning the full range in UV emission line properties. We then quantify the evolution of C IV and He II emission line properties with redshift. The increase in C IV blueshifts with cosmological redshift can be fully explained by the higher luminosities of quasars observed at high redshifts. We recover broadly similar trends between the He II EW and C IV blueshift at both 1.5<z<2.65 and 3.5<z<4.0 suggesting that the blueshift depends systematically on the spectral energy density (SED) of the quasar and there is no evolution in the SED over the redshift range 1.5<z<4.0. C IV blueshifts are highest when L/LEdd > 0.2 and Mbh > 10^9 Mo for the entire 1.5<z<4.0 sample. We find that luminosity matching samples as a means to explore the evolution of their rest-frame UV emission line properties is only viable if the samples are also matched in the Mbh - L/LEdd plane. Quasars at z>6 are on average less massive and have higher Eddington-scaled accretion rates than their luminosity-matched counterparts at 1.5<z<4.0, which could explain the observed evolution in their UV line properties.

Previous studies of the exoplanet LTT 1445Ac concluded that the light curve from the Transiting Exoplanet Survey Satellite (TESS) was consistent with both grazing and non-grazing geometries. As a result, the radius and hence density of the planet remained unknown. To resolve this ambiguity, we observed the LTT 1445 system for six spacecraft orbits of the Hubble Space Telescope (HST) using WFC3/UVIS imaging in spatial scan mode, including one partial transit of LTT 1445Ac. This imaging produces resolved light curves of each of the three stars in the LTT 1445 system. We confirm that the planet transits LTT 1445A and that LTT 1445C is the source of the rotational modulation seen in the TESS light curve, and we refine the estimate of the dilution factor for the TESS data. We perform a joint fit to the TESS and HST observations, finding that the transit of LTT 1445Ac is not grazing with 97% confidence. We measure a planetary radius of 1.10$_{-0.07}^{+0.10}$ R$_\oplus$. Combined with previous radial velocity observations, our analysis yields a planetary mass of $1.36\pm0.19$ M$_\oplus$ and a planetary density of 5.6$_{-1.5}^{+1.7}$ g cm$^{-3}$. LTT 1445Ac is an Earth analog with respect to its mass and radius, albeit with a higher instellation, and is therefore an exciting target for future atmospheric studies.

Ground-based gamma-ray astronomy is a powerful tool to study cosmic-ray physics, providing a diagnostic of the high-energy processes at work in the most extreme astrophysical accelerators of the universe. Ground-based gamma-ray detectors apply a number of experimental techniques to measure the products of air showers induced by the primary gamma-rays over a wide energy range, from about 30 GeV to few PeV. These are based either on the measurement of the atmospheric Cherenkov light induced by the air showers, or the direct detection of the shower's secondary particles at ground level. Thanks to the recent development of new and highly sensitive ground-based gamma-ray detectors, important scientific results are emerging which motivate new experimental proposals, at various stages of implementation. In this chapter we will present the current expectations for future experiments in the field.

We study the cosmological effects of a Galileon term in scalar-tensor theories of gravity. The subset of scalar-tensor theories considered are characterized by a non-minimal coupling $F(\sigma) R$, a kinetic term with arbitrary sign $Z (\partial \sigma)^2$ with $Z = \pm 1$, a potential $V(\sigma)$, and a Galileon term $G_3(\sigma, (\partial \sigma)^2) \square \sigma$. In addition to the modified dynamics, the Galileon term provides a screening mechanism to potentially reconcile the models with General Relativity predictions inside a Vainshtein radius. Thanks to the Galileon term, the stability conditions, namely ghost and Laplacian instabilities, in the branch with a negative kinetic term ($Z = -1$) are fulfilled for a large volume of the parameter space. Solving numerically the background evolution and linear perturbations, we derive the constraints on the cosmological parameters in presence of a Galileon term for different combination of the cosmic microwave background (CMB) data from Planck, baryon acoustic oscillations (BAO) measurements from BOSS, and supernovae from the Pantheon compilation. We find that the Galileon term alters the dynamics of all the studied cases. For a standard kinetic term ($Z = 1$), we find that Planck data and a compilation of BAO data constrain the Galileon term to small values that allow screening very inefficiently. For a negative kinetic term ($Z = -1$), a Galileon term and a non-zero potential lead to an efficient screening in a physically viable regime of the theory, with a value for the Hubble constant today which alleviates the tension between its CMB and local determinations. For a vanishing potential, the case with $Z=-1$ and the Galileon term driving the late acceleration of the Universe is ruled out by Planck data.

GRB 230307A is an extremely bright long duration GRB with an observed gamma-ray fluence of $\gtrsim$3$\times$10$^{-3}$ erg cm$^{-2}$ (10--1000 keV), second only to GRB 221009A. Despite its long duration, it is possibly associated with a kilonova, thus resembling the case of GRB 211211A. In analogy with GRB 211211A, we distinguish three phases in the prompt gamma-ray emission of GRB 230307A: an initial short duration, spectrally soft emission; a main long duration, spectrally hard burst; a temporally extended and spectrally soft tail. We intepret the initial soft pulse as a bright precursor to the main burst and compare its properties with models of precursors from compact binary mergers. We find that to explain the brightness of GRB 230307A, a magnetar-like ($\gtrsim 10^{15}$ G) magnetic field should be retained by the progenitor neutron star. Alternatively, in the post-merger scenario, the luminous precursor could point to the formation of a rapidly rotating massive neutron star.

The purpose of this work is to demonstrate the effectiveness of ground-based support for space missions, primarily PSP, using large Ukrainian decameter radio telescopes. Another goal of the work is to carry out cross calibration of the radiometers onboard spacecraft using the calibrated data of the ground-based radio telescopes. One of the most common methods of remote diagnostics of the solar corona is the study of radio emission, the sources of which are located in the solar corona at different heliocentric altitudes. The technique of joint space terrestrial observations consists in the simultaneous observation of individual events and their analysis in the widest possible frequency band during the maximum approach of the PSP vehicle to the Sun. At the same time, observation in the common frequency band is proposed to be used for calibration of the onboard radio receivers. The methods of planning joint space terrestrial observations are substantiated. Using the data of the UTR 2, URAN 2 radio telescopes and the PSP probe, the dynamic and polarization spectra of the simultaneously observed bursts on June 9, 2020 were obtained. The identification and comparison of individual bursts was carried out. A common dynamic spectrum of the bursts in the frequency band 0.5 ... 32 MHz was obtained. Cross calibration of the HFR receiver of the FIELDS PSP module in the frequency band 10...18 MHz was made using the calibrated data of terrestrial radio telescopes. The effectiveness of ground-based support of the PSP mission by the large Ukrainian radio telescopes is shown. Examples of joint observations are given, and the method of cross calibration of the FIELD PSP module receivers is demonstrated. Prospects for further ground based support for solar space missions are presented

In this invited review I discuss the calibration and applications of the period-luminosity relation of classical Cepheid and RR Lyrae stars. After a brief introduction, starting with results from Hipparcos and discussing some post-Hipparcos era developments, I focus on recent results using Gaia Data Release 3 data. I present an overview of the most recent period-luminosity relations, a discussion and some new results on Cepheids in open clusters. I also discuss the effect of reddening and that the use of Wesenheit indices is actually an oversimplification to dealing with the problem of reddening.

The gravitational interactions between the Milky Way and in-falling satellites offer a wealth of information about the formation and evolution of our Galaxy. In this paper, we explore the high-dimensionality of the GALAH DR3 plus Gaia eDR3 data set to identify new tidally stripped candidate stars of the nearby star cluster Omega Centauri ($\omega\,\mathrm{Cen}$). We investigate both the chemical and dynamical parameter space simultaneously, and identify cluster candidates that are spatially separated from the main cluster body, in regions where contamination by halo field stars is high. Most notably, we find candidates for $\omega\,\mathrm{Cen}$ scattered in the halo extending to more than $50^{\circ}$ away from the main body of the cluster. Using a grid of simulated stellar streams generated with $\omega\,\mathrm{Cen}$ like orbital properties, we then compare the on sky distribution of these candidates to the models. The results suggest that if $\omega\,\mathrm{Cen}$ had a similar initial mass as its present day mass, then we can place a lower limit on its time of accretion at t$_{\mathrm{acc}} > 7$ Gyr ago. Alternatively, if the initial stellar mass was significantly larger, as would be expected if $\omega\,\mathrm{Cen}$ is the remnant core of a dwarf Galaxy, then we can constrain the accretion time to t$_{\mathrm{acc}} > 4$ Gyr ago. Taken together, these results are consistent with the scenario that $\omega\,\mathrm{Cen}$ is the remnant core of a disrupted dwarf galaxy.

We consider a cosmic string moving through a gas of superfluid dark matter (SFDM) particles and analyze how it affects the dark matter distribution. We look at two different cases: first, a cosmic string passing through an already condensed region, and second, through a region that is not yet condensed. In the former, the string induces a weak shock in the superfluid, and the Bose-Einstein condensate (BEC) survives. In the latter, a wake of larger density is formed behind the string, and we study under which conditions a BEC can be formed in the virialized region of the wake. By requiring the thermalization of the DM particles and the overlap of their de Broglie wavelengths inside the wake, we obtain an upper bound on the mass of the dark matter particles on the order of 10 eV, which is compatible with typical SFDM models.

The area ratios of sunspots to white light faculae in the first two years of sunspot cycles 12-21 correlate remarkably well with the peak amplitudes of those cycles between 1878-1980 (Brown and Evans, 1980). This finding could not be used to predict subsequent cycle amplitudes because the Royal Greenwich Observatory program of facular area measurements was discontinued in 1976. We use continuum images from the Michelson Doppler Imager (MDI) and the Heliospheric and Magnetic Imager (HMI) to show that the close relation holds also for cycle 24, and we predict an amplitude of approximately 185 for the current cycle 25.

We show that the baryon asymmetry of the Universe cannot be explained by a large initial value before inflation because it inevitably predicts correlated baryon isocurvature perturbations that are already excluded by cosmic microwave background observations. Similar arguments can generally be applied to some models of dark matter.

Luminous blue variables (LBVs) are evolved massive stars close to the Eddington limit, with a distinct spectroscopic and photometric variability having unsteady mass-loss rates. These stars show a considerable change in their surface temperature from quiescent to outbursts phase. The cause of irregular variability and unsteady mass-loss rate is not properly understood. Here we present the result of linear stability analysis in two LBVs AF And and R 127 during their quiescent and outburst phase. We note that several modes are unstable in the models of the considered LBVs. Mode interaction is frequent in the modal diagrams for the models of both LBVs. For AF And, number of instabilities increase in models having temperature below 15000 K. The found instabilities may be linked with the observed irregular variabilities and surface eruptions. Observational facilities of Belgo-Indian Network for Astronomy and Astrophysics (BINA) will be very beneficial to study the spectroscopic and photometric behavior of the considered LBVs.

One of the most well-studied exoplanets to date, HD 189733 b, stands out as an archetypal hot Jupiter with many observations and theoretical models aimed at characterizing its atmosphere, interior, host star, and environment. We report here on the results of an extensive campaign to observe atmospheric escape signatures in HD 189733 b using the Hubble Space Telescope and its unique ultraviolet capabilities. We have found a tentative, but repeatable in-transit absorption of singly-ionized carbon (C II, $5.2\% \pm 1.4\%$) in the epoch of June-July/2017, as well as a neutral hydrogen (H I) absorption consistent with previous observations. We model the hydrodynamic outflow of HD 189733 b using an isothermal Parker wind formulation to interpret the observations of escaping C and O nuclei at the altitudes probed by our observations. Our forward models indicate that the outflow of HD 189733 b is mostly neutral within an altitude of $\sim 2$ R$_\mathrm{p}$ and singly ionized beyond that point. The measured in-transit absorption of C II at 133.57 nm is consistent with an escape rate of $\sim 1.1 \times 10^{11}$ g$\,$s$^{-1}$, assuming solar C abundance and outflow temperature of $12\,100$ K. Although we find a marginal neutral oxygen (O I) in-transit absorption, our models predict an in-transit depth that is only comparable to the size of measurement uncertainties. A comparison between the observed Lyman-$\alpha$ transit depths and hydrodynamics models suggests that the exosphere of this planet interacts with a stellar wind at least one order of magnitude stronger than solar.

We investigate the production of primordial black holes (PBHs) in a mixed inflaton--curvaton scenario with a quadratic curvaton potential, assuming the curvaton is in de Sitter equilibrium during inflation with $\langle \chi\rangle =0$. In this setup, the curvature perturbation sourced by the curvaton is strongly non-Gaussian, containing no leading Gaussian term. We show that for $m^2/H^2\gtrsim 0.3$, the curvaton contribution to the spectrum of primordial perturbations on CMB scales can be kept negligible but on small scales the curvaton can source PBHs. In particular, PBHs in the asteroid mass range $10^{-16}M_{\odot}\lesssim M\lesssim 10^{-10}M_{\odot}$ with an abundance reaching $f_{\rm PBH} = 1$ can be produced when the inflationary Hubble scale $H\gtrsim 10^{12}$ GeV and the curvaton decay occurs in the window from slightly before the electroweak transition to around the QCD transition.

It has recently been speculated that the NANOGrav observations point towards a first-order phase transition in the dark sector at the GeV scale [1]. Here we show that such a high energy sibling of the New Early Dark Energy (NEDE) phase transition might already have been predicted in the Hot New Early Dark Energy model (Hot NEDE) [2, 3]. There, it was argued that the NEDE phase transition is a signature of neutrino mass generation through the inverse seesaw mechanism. Specifically, it serves a double purpose by resolving the Hubble tension through an energy injection and generating the Majorana mass entry in the inverse seesaw mixing matrix. In addition, the usual NEDE phase transition is accompanied by a UV counterpart, which generates the heavy Dirac mass entry in the inverse seesaw mass matrix of a right-handed neutrino. Here, we investigate if the UV phase transition of the Hot NEDE model can occur at the GeV scale in view of the recent NANOGrav observations.

Assuming the evidence for gravitational wave background from recent data release of pulsar timing arrays to be interpreted as the scalar-induced gravitational waves (SIGWs), we study the second and third order gravitational waves simultaneously, by jointly analyzing a combination of PTA, big-bang nucleosynthesis (BBN), and cosmic microwave background (CMB) datasets. We obtain the primordial curvature spectral amplitude $0.014<A_\zeta<0.058$ and the spectral peak frequency $10^{-7.4}\ \mathrm{Hz}<f_\ast<10^{-6.4}\ \mathrm{Hz}$ at 95\% confidence level, indicating a mass range of primordial black holes (PBHs) $10^{-4.3}M_\odot<m_{\mathrm{pbh}}<10^{-2.3}M_\odot$. We further find that the third order gravitational waves have greater contributions to the integrated energy density than the second order gravitational waves when $A_\zeta\gtrsim0.06$. In addition, we anticipate that future PTA projects can not only test the above results, but also have powerful abilities to explore the origin and evolution of the universe, particularly, the inflation.

We describe the ensemble properties of the $1.9 < z < 3.5$ Lyman Alpha Emitters (LAEs) found in the HETDEX survey's first public data release, HETDEX Public Source Catalog 1 (Mentuch Cooper et al. 2023). Stacking the low-resolution ($R \sim$ 800) spectra greatly increases the signal-to-noise ratio, revealing spectral features otherwise hidden by noise, and we show that the stacked spectrum is representative of an average member of the set. The flux limited, Ly$\alpha$ signal-to-noise ratio restricted stack of 50K HETDEX LAEs shows the ensemble biweight ``average" $z \sim 2.6$ LAE to be a blue (UV continuum slope $\sim -2.4$ and E(B-V) $< 0.1$), moderately bright (M$_{\text{UV}} \sim -19.7$) star forming galaxy with strong Ly$\alpha$ emission (log $L_{Ly\alpha}$ $\sim$ 42.8 and $W_{\lambda}$(Ly$\alpha$) $\sim$ 114\AA), and potentially significant leakage of ionizing radiation. The restframe UV light is dominated by a young, metal poor stellar population with an average age 5-15 Myr and metallicity of 0.2-0.3 Z$_{\odot}$.

Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single stars, focusing on the identification of final outcomes and their long-term evolution and observational properties, using the moving-mesh hydrodynamics code AREPO. This type of encounters produces five types of outcomes: stellar disruption, stellar collision, weak perturbation of the original binary, binary member exchange, and triple formation. The two decisive parameters are the binary phase angle, which determines which two objects meet at the first closest approach, and the impact parameter, which sets the boundary between violent and non-violent interactions. When the impact parameter is smaller than the semimajor axis of the binary, tidal disruptions and star-BH collisions frequently occur when the BH and the incoming star first meet, while the two stars mostly merge when the two stars meet first instead. In both cases, the BHs accrete from an accretion disk at super-Eddington rates, possibly generating flares luminous enough to be observed. The stellar collision products either form a binary with the BH or remain unbound to the BH. Upon collision, the merged stars are hotter and larger than main sequence stars of the same mass at similar age. Even after recovering their thermal equilibrium state, stellar collision products, if isolated, would remain hotter and brighter than main sequence stars until becoming giants.

The Surface Brightness Fluctuation (SBF) method is a powerful tool for determining distances to early-type galaxies. The method measures the intrinsic variance in a galaxy's surface brightness distribution to determine its distance with an accuracy of about 5%. Here, we discuss the mathematical formalism behind the SBF technique, its calibration, and the practicalities of how measurements are performed. We review the various sources of uncertainties that affect the method and discuss how they can be minimized or controlled through careful observations and data analysis. The SBF technique has already been successfully applied to a large number of galaxies and used for deriving accurate constraints on the Hubble-Lema\^itre constant $H_0$. An approved JWST program will greatly reduce the systematic uncertainties by establishing a firm zero-point calibration using tip of the red giant branch (TRGB) distances. We summarize the existing results and discuss the excellent potential of the SBF method for improving the current constraints on $H_0$.

One of the goals of gravitational-wave astrophysics is to infer the number and properties of the formation channels of binary black holes (BBHs); to do so, one must be able to connect various models with the data. We explore benefits and potential issues with analyses using models informed by population synthesis. We consider 5 possible formation channels of BBHs, as in Zevin et al. (2021b). First, we confirm with the GWTC-3 catalog what Zevin et al. (2021b) found in the GWTC-2 catalog, i.e. that the data are not consistent with the totality of observed BBHs forming in any single channel. Next, using simulated detections, we show that the uncertainties in the estimation of the branching ratios can shrink by up to a factor of $\sim 1.7$ as the catalog size increases from $50$ to $250$, within the expected number of BBH detections in LIGO-Virgo-KAGRA's fourth observing run. Finally, we show that this type of analysis is prone to significant biases. By simulating universes where all sources originate from a single channel, we show that the influence of the Bayesian prior can make it challenging to conclude that one channel produces all signals. Furthermore, by simulating universes where all 5 channels contribute but only a subset of channels are used in the analysis, we show that biases in the branching ratios can be as large as $\sim 50\%$ with $250$ detections. This suggests that caution should be used when interpreting the results of analyses based on strongly modeled astrophysical sub-populations.

The pulsar timing array (PTA) collaborations have recently provided compelling evidence for the presence of a stochastic signal consistent with a gravitational-wave background. In this letter, we combine the latest data sets from NANOGrav, PPTA and EPTA collaborations to explore the cosmological interpretations for the detected signal from first-order phase transitions, domain walls and cosmic strings separately. We find that the first-order phase transitions and cosmic strings can give comparable interpretations compared to supermassive black hole binaries (SMBHBs) characterized by a power-law spectrum, but the domain wall model is strongly disfavored with the Bayes factor compared to the SMBHBs model being 0.009. Furthermore, the constraints on the parameter spaces indicate that: 1) a strong phase transition at temperatures below the electroweak scale is favored and the bubble collisions make the dominant contribution to the energy spectrum; 2) the cosmic string tension is $G \mu \in [1.46, 15.3]\times 10^{-12}$ at $90\%$ confidence interval and a small reconnection probability $p<6.68\times 10^{-2}$ is preferred at $95\%$ confidence level, implying that the strings in (super)string theory are strongly preferred than the classical field strings.

We present the discovery of Type II supernova (SN) 2023ixf in M101, among the closest core-collapse SNe in the last several decades, and follow-up photometric and spectroscopic observations in the first month of its evolution. The light curve is characterized by a rapid rise ($\approx5$ days) to a luminous peak ($M_V\approx-18$ mag) and plateau ($M_V\approx-17.6$ mag) extending to $30$ days with a smooth decline rate of $\approx0.03$ mag day$^{-1}$. During the rising phase, $U-V$ color shows blueward evolution, followed by redward evolution in the plateau phase. Prominent flash features of hydrogen, helium, carbon, and nitrogen dominate the spectra up to $\approx5$ days after first light, with a transition to a higher ionization state in the first $\approx2$ days. Both the $U-V$ color and flash ionization states suggest a rise in the temperature, indicative of a delayed shock-breakout inside dense circumstellar material (CSM). From the timescales of CSM interaction, we estimate its compact radial extent of $\sim(3-7)\times10^{14}$ cm. We then construct numerical light-curve models based on both continuous and eruptive mass-loss scenarios shortly before explosion. For the continuous mass-loss scenario, we infer a range of mass-loss history with $0.1-1.0$ $M_\odot {\rm yr}^{-1}$ in the final $2-1$ years before explosion, with a potentially decreasing mass loss of $0.01-0.1$ $M_\odot {\rm yr}^{-1}$ in $\sim0.7-0.4$ years towards the explosion. For the eruptive mass-loss scenario, we favor eruptions releasing $0.3-1$ $M_\odot$ of the envelope at about a year before explosion, which result in CSM with mass and extent similar to the continuous scenario. We discuss the implications of the available multi-wavelength constraints obtained thus far on the progenitor candidate and SN 2023ixf to our variable CSM models.

ADITYA-L1 is India's first dedicated mission to observe the sun and its atmosphere from a halo orbit around L1 point. Visible emission line coronagraph (VELC) is the prime payload on board at Aditya-L1 to observe the sun's corona. VELC is designed as an internally occulted reflective coronagraph to meet the observational requirements of wide wavelength band and close to the solar limb (1.05 Ro). Images of the solar corona in continuum and spectra in three emission lines 5303{\AA} [Fe xiv], 7892{\AA} [Fe xi] and 10747 [Fe xiii] obtained with high cadence to be analyzed using software algorithms automatically. A reasonable part of observations will be made in synoptic mode, those, need to be analyzed and results made available for public use. The procedure involves the calibration of instrument and detectors, converting the images into fits format, correcting the images and spectra for the instrumental effects, align the images etc. Then, develop image processing algorithms to detect the occurrence of energetic events using continuum images. Also derive physical parameters, such as temperature and velocity structure of solar corona using emission line observations. Here, we describe the calibration of detectors and the development of software algorithms to detect the occurrence of CMEs and analyze the spectroscopic data.

We present new spectroscopic and photometric follow-up observations of the known sample of extreme coronal line emitting galaxies (ECLEs) identified in the Sloan Digital Sky Survey (SDSS). With these new data, observations of the ECLE sample now span a period of two decades following their initial SDSS detections. We confirm the nonrecurrence of the iron coronal line signatures in five of the seven objects, further supporting their identification as the transient light echoes of tidal disruption events (TDEs). Photometric observations of these objects in optical bands show little overall evolution. In contrast, mid-infrared (MIR) observations show ongoing long-term declines. The remaining two objects had been classified as active galactic nuclei (AGN) with unusually strong coronal lines rather than being TDE related, given the persistence of the coronal lines in earlier follow-up spectra. We confirm this classification, with our spectra continuing to show the presence of strong, unchanged coronal-line features and AGN-like MIR colours and behaviour. We have constructed spectral templates of both subtypes of ECLE to aid in distinguishing the likely origin of newly discovered ECLEs. We highlight the need for higher cadence, and more rapid, follow-up observations of such objects to better constrain their properties and evolution. We also discuss the relationships between ECLEs, TDEs, and other identified transients having significant MIR variability.

The Cosmic Neutrino Background (CNB) encodes a wealth of information, but has not yet been observed directly. To determine the prospects of detection and to study its information content, we reconstruct the phase-space distribution of local relic neutrinos from the three-dimensional distribution of matter within 200 Mpc/h of the Milky Way. Our analysis relies on constrained realization simulations and forward modelling of the 2M++ galaxy catalogue. We find that the angular distribution of neutrinos is anti-correlated with the projected matter density, due to the capture and deflection of neutrinos by massive structures along the line of sight. Of relevance to tritium capture experiments, we find that the gravitational clustering effect of the large-scale structure on the local number density of neutrinos is more important than that of the Milky Way for neutrino masses less than 0.1 eV. Nevertheless, we predict that the density of relic neutrinos is close to the cosmic average, with a suppression or enhancement over the mean of (-0.3%, +7%, +27%) for masses of (0.01, 0.05, 0.1) eV. This implies no more than a marginal increase in the event rate for tritium capture experiments like PTOLEMY. We also predict that the CNB and CMB rest frames coincide for 0.01 eV neutrinos, but that neutrino velocities are significantly perturbed for masses larger than 0.05 eV. Regardless of mass, we find that the angle between the neutrino dipole and the ecliptic plane is small, implying a near-maximal annual modulation in the bulk velocity. Along with this paper, we publicly release our simulation data, comprising more than 100 simulations for six different neutrino masses.

We demonstrate that existing gravitational wave data from LIGO already places constraints on well motivated Pati-Salam models that allow the Standard Model to be embedded within grand unified theories. For the first time in these models we also constrain the parameter space by requiring that the phase transition completes, with the resulting constraint being competitive with the limits from LIGO data. Both constraints are complementary to the LHC constraints and can exclude scenarios that are much heavier than can be probed in colliders. Finally we show that results from future LIGO runs, and the planned Einstein telescope, will substantially increase the limits we place on the parameter space.

Recently, the piercing of a mini boson star by a black hole was studied, with tidal capture and the discovery of a "gravitational atom" being reported ( arXiv:2206.00021 [gr-qc] ). Building on this research, we extend the study by including a hexic solitonic potential and explore the piercing of a solitonic boson star by a black hole. Notably, the solitonic boson star can reach higher compactness, which one might expect could alter the dynamics in this context. Our findings suggest that even when the black hole's size approaches the test particle limit, the solitonic boson star is easily captured by the black hole due to an extreme tidal capture process. Regardless of the black hole initial mass and velocity, our results indicate that over 85% of the boson star material is accreted. Thus, the self-interaction does not alter the qualitative behavior of the system.

Ultracold atomic gases can undergo phase transitions that mimic relativistic vacuum decay, allowing us to empirically test early-Universe physics in tabletop experiments. We investigate the physics of these analogue systems, going beyond previous analyses of the classical equations of motion to study quantum fluctuations in the cold-atom false vacuum. We show that the fluctuation spectrum of this vacuum state agrees with the usual relativistic result in the regime where the classical analogy holds, providing further evidence for the suitability of these systems for studying vacuum decay. Using a suite of semiclassical lattice simulations, we simulate bubble nucleation from this analogue vacuum state in a 1D homonuclear potassium-41 mixture, finding qualitative agreement with instanton predictions. We identify realistic parameters for this system that will allow us to study vacuum decay with current experimental capabilities, including a prescription for efficiently scanning over decay rates, and show that this setup will probe the quantum (rather than thermal) decay regime at temperatures $T\lesssim10\,\mathrm{nK}$. Our results help lay the groundwork for using upcoming cold-atom experiments as a new probe of nonperturbative early-Universe physics.

This dissertation examines the impact of quantum gravity on electromagnetism and its backreaction, using perturbative general relativity as an effective field theory. Our analysis involves quantum-correcting Maxwell's equations to obtain a gauge-independent, real, and causal effective field equation that describes quantum gravitational effects on electromagnetism. Additionally, we present a perturbative mechanism through which quantum gravity induces a dimension six coupling between a massive scalar and electromagnetism. To investigate the effects of electromagnetism on the gravitational sector, we derive an exact, dimensionally regulated, Fourier mode sum for the Lorentz gauge propagator of a massive photon on an arbitrary cosmological background supported by a scalar inflaton. This allows us to calculate the effective potential induced by photons. Finally, we use a similar Fourier mode sum for a time-dependent mass to study the effective force on the inflaton 0-mode and its impact on reheating.

A prominent bottleneck for advancing our understanding of astrophysical turbulence is the limited resolution of numerical simulations, which inhibits fully sampling scales in the inertial range. Machine learning (ML) techniques have demonstrated promise in up-scaling resolution in both image analysis and numerical simulations (i.e., superresolution). Here we employ and further develop a physics-constrained convolutional neural network (CNN) ML model called "MeshFreeFlowNet'' for superresolution studies of turbulent systems. The MeshFreeFlowNet CNN is trained both on the simulation images as well as the evaluated PDEs, making it sensitive to the underlying physics of a particular fluid system. In particular, we aim to generate a superresolution framework for 2D turbulent Rayleigh-B\'enard convection (RBC) generated with the Dedalus code. We modify the MeshFreeFlowNet architecture to include the full set of simulation PDEs and the boundary conditions. Our training set includes fully developed turbulence sampling Rayleigh numbers (Ra) of Ra=10^6-10^{10}. We evaluate the success of the learned simulations by comparing the direct Dedalus simulation power spectra to the predicted CNN output power spectra. We compare both ground truth and predicted power spectral inertial range scalings to theoretical predictions. We find that the network performs well at all Ra studied here in recovering large-scale information, including the inertial range slopes. We find that our updated architecture performs well on laminar and turbulent flows, but the results towards smaller-scales are significantly better as the flow transitions to more turbulent regimes. This is likely because more turbulent systems have a rich variety of structures at many length scales compared to laminar flows. We also find that the superresolution prediction is overly dissipative at smaller scales than that of the inertial range.

I investigate the sensitivity of gravitational-wave searches by analyzing the response of matched filters in stationary Gaussian noise. In particular, I focus on the ability to analytically model the distribution of observed filter responses maximized over coalescence phase and/or a template bank as well as the response of statistics defined for a network of detectors. Semianalytic sensitivity estimates derived assuming stationary Gaussian noise are compared to sensitivity estimates obtained from real searches processing real noise, which is neither perfectly stationary nor perfectly Gaussian. I find that semianalytic estimates are able to reproduce real search sensitivity for the LIGO-Virgo-KAGRA Collaboration's third observing run with high fidelity. I also discuss how to select computational speed-ups (hopeless signal-to-noise ratio cuts) and make predictions for the fourth observing run using projected detector sensitivities.

Palatini $F(R)$ gravity proved to be a powerful tool in order to realize asymptotically flat inflaton potentials. Unfortunately, it also inevitably implies higher-order inflaton kinetic terms in the Einstein frame that might jeopardize the evolution of the system out of the slow-roll regime. We prove that a $F(R-X)$ gravity, where $X$ is the inflaton kinetic term, solves the issue. Moreover, when $F$ is a quadratic function such a choice easily leads to a new class of inflationary attractors, fractional attractors, that generalizes the already well-known polynomial $\alpha$-attractors.

The recent observation of the object HESS J1731-347 suggests the existence of a very light and very compact neutron star being a challenge for commonly used equation of state for dense matter. In this work we present a relativistic mean field model enriched with meson crossing terms among isovector and isoscalar mesons. Such interactions particularly dominate the behavior of the symmetry energy and accounts for small size of compact star radius. The proposed model fulfill the recent constraints concerning the symmetry energy slope and state-of-the-art compact stars constraints derived from the NICER measurements of PSR J0030+0451 and PSR J0740+6620 pulsars as well as from the GW170817 event and its associated electromagnetic counterparts AT2017gfo/GRB170817A.

In this paper, we study the magnetic reconnection process of energy extraction from a rapidly rotating Kerr-Newman-MOG black hole by investigating the combined effect of black hole charge and the MOG parameter. We explore the energy efficiency of energy extraction and power. Based on an attractive gravitational charge of the MOG parameter $\alpha$ that physically manifests to strengthen black hole gravity we show that the combined effect of the MOG parameter and black hole charge can play an increasingly important role and accordingly lead to high energy efficiency and power for the energy extraction via the magnetic reconnection. Further, we study to estimate the rate of energy extraction under the fast magnetic reconnection by comparing the power of the magnetic reconnection and Blandford-Znajek (BZ) mechanisms. We show that the rate of energy extraction increases as a consequence of the combined effect of black hole charge and MOG parameter. It suggests that magnetic reconnection is significantly more efficient than BZ. In fact, the magnetic reconnection is fueled by magnetic field energy due to the twisting of magnetic field lines around the black hole for the plasma acceleration, and thus MOG parameter gives rise to even more fast spin that can strongly change the magnetic field reconfiguration due to the frame dragging effect. This is how energy extraction is strongly enhanced through the magnetic reconnection, thus making the energy extraction surprisingly more efficient for the Kerr-Newman-MOG black hole than Kerr black hole under the combined effect of black hole charge and MOG parameter.

Testing gravity theories and their parameters using observations is an important issue in relativistic astrophysics. In this context, we investigate the motion of test particles and their harmonic oscillations in the spacetime of non-rotating hairy black holes (BHs) in Hordeski gravity, together with astrophysical applications of quasiperiodic oscillations (QPOs). We show possible values of upper and lower frequencies of twin-peak QPOs which may occur in the orbits from innermost stable circular orbits to infinity for various values of the Horndeski parameter $q$ in relativistic precession, warped disk models, and three different sub-models of the epicyclic resonant model. We also study the behaviour of the QPO orbits and their position relative to innermost stable circular orbits (ISCOs) with respect to different values of the parameter $q$. {It is obtained that at a critical value of the Horndeski parameter ISCO radius takes $6M$ which has been in the pure Schwarzschild case.} Finally, we obtain mass constraints of the central BH of microquasars GRS 1915+105 and XTE 1550-564 at the GR limit and the possible value of the Horndeski parameter in the frame of the above-mentioned QPO models. The analysis of orbits of twin peak QPOs with the ratio of upper and lower frequencies 3:2, around the BHs in the frame of relativistic precession (RP) and epicyclic resonance (ER4) QPO models have shown that the orbits locate close to the ISCO. The distance between QPO orbits and ISCO is obtained to be less than the error of the observations.

We study the inflationary model with a spectator scalar field $\chi$ coupled to both the inflaton and Ricci scalar. The interaction between the $\chi$ field and the gravity, denoted by $\xi R\chi^2$, can trigger the tachyonic instability of certain modes of the $\chi$ field. As a result, the $\chi$ field perturbations are amplified and serve as a gravitational wave (GW) source. When considering the backreaction of the $\chi$ field, an upper bound on the coupling parameter $\xi$ must be imposed to ensure that inflation does not end prematurely. In this case, we find that the inflaton's evolution experiences a sudden slowdown due to the production of $\chi$ particles, resulting in a unique oscillating structure in the power spectrum of curvature perturbations at specific scales. Moreover, the GW signal induced by the $\chi$ field is more significant than primordial GWs at around its peak scale, leading to a noticeable bump in the overall energy spectrum of GWs.

Models that produce Axion-Like-Particles (ALP) after cosmological inflation due to spontaneous $U(1)$ symmetry breaking also produce cosmic string networks. Those axionic strings lose energy through gravitational wave emission during the whole cosmological history, generating a stochastic background of gravitational waves that spans many decades in frequency. We can therefore constrain the axion decay constant and axion mass from limits on the gravitational wave spectrum and compatibility with dark matter abundance as well as dark radiation. We derive such limits from analyzing the most recent NANOGrav data from Pulsar Timing Arrays (PTA). The limits are compatible with the slightly stronger $N_{\rm eff}$ bounds on dark radiation for ALP masses $m_a \lesssim 10^{-10}$ eV. On the other hand, for heavy ALPs with $m_a\gtrsim 0.1$ GeV and $N_{\rm DW}\neq 1$, new regions of parameter space can be probed by PTA data due to the dominant Domain-Wall contribution to the gravitational wave background.

Low energy imprints of modifications to general relativity are often found in pressure balance equations inside stars. These modifications are then amenable to tests via astrophysical phenomena, using observational effects in stellar astrophysics that crucially depend on such equations. One such effect is tidal disruption of stars in the vicinity of black holes. In this paper, using a numerical scheme modelled with smoothed particle hydrodynamics, we study real time tidal disruption of a class of white dwarfs by intermediate-mass black holes, in the low energy limit of a theory of modified gravity that alters the internal physics of white dwarfs, namely the Eddington inspired Born-Infeld theory. In this single parameter extension of general relativity, the mass-radius relation of white dwarfs as well as their tidal disruption radius depend on the modified gravity parameter, and these capture the effect of modifications to general relativity. Our numerical simulations incorporating these show that departure from general relativity in these scenarios might be observationally significant, and should therefore be contrasted with data. In particular, we study observationally relevant physical quantities, i.e., tidal kick velocity and trajectory deviation of the remnant core and fallback rates of the tidal debris in this theory and compare them to the Newtonian limit of general relativity. We also comment on the qualitative differences between the modified gravity theory and one with stellar rotation.

Proto-neutron stars formed during core-collapse supernovae are hot and dense environments that contain a sizable population of muons. If these interact with new long-lived particles with masses up to roughly 100 MeV, the latter can be produced and escape from the stellar plasma, causing an excessive energy loss constrained by observations of SN 1987A. In this article we calculate the emission of light dark fermions that are coupled to leptons via a new massive vector boson, and determine the resulting constraints on the general parameter space. We apply these limits to the gauged $L_\mu-L_\tau$ model with dark fermions, and show that the SN 1987A constraints exclude a significant portion of the parameter space targeted by future experiments. We also extend our analysis to generic effective four-fermion operators that couple dark fermions to muons, electrons, or neutrinos. We find that SN 1987A cooling probes a new-physics scale up to $\sim7$ TeV, which is an order of magnitude larger than current bounds from laboratory experiments.

The purpose of this review it to present a renewed perspective of the problem of self-gravitating elastic bodies under spherical symmetry. It is also a companion to the papers [Phys. Rev. D105, 044025 (2022)], [Phys. Rev. D106, L041502 (2022)], and [arXiv:2306.16584 [gr-qc]], where we introduced a new definition of spherically symmetric elastic bodies in general relativity, and applied it to investigate the existence and physical viability, including radial stability, of static self-gravitating elastic balls. We focus on elastic materials that generalize fluids with polytropic, linear, and affine equations of state, and discuss the symmetries of the energy density function, including homogeneity and the resulting scale invariance of the TOV equations. By introducing invariant characterizations of physical admissible initial data, we numerically construct mass-radius-compactness diagrams, and conjecture about the maximum compactness of stable physically admissible elastic balls.

Epicyclic frequencies are usually observed in X-ray binaries and constitute a powerful astrophysical mean to probe the strong gravitational field around a compact object. We consider them in the equatorial plane around a general stationary and axially symmetric wormhole. We first search for the wormholes' existence, distinguishing them from a Kerr black hole. Once there will be available observational data on wormholes, we present a strategy to reconstruct the related metrics. Finally, we discuss the implications of our approach and outline possible future perspectives.

We investigate the interplay between numerical relativity (NR) and point-particle black hole perturbation theory (ppBHPT) in the comparable mass regime. Specifically, we reassess the $\alpha$-$\beta$ scaling technique, previously introduced by Islam et al, as a means to effectively match ppBHPT waveforms to NR waveforms within this regime. Utilizing publicly available long NR data for a mass ratio of $q=3$ (where $q:=m_1/m_2$ represents the mass ratio of the binary, with m_1 and m_2 denoting the masses of the primary and secondary black holes, respectively), encompassing the final $\sim 65$ orbital cycles of the binary evolution, we examine the range of applicability of such scalings. We observe that the scaling technique remains effective even during the earlier stages of the inspiral. Additionally, we provide commentary on the temporal evolution of the $\alpha$ and $\beta$ parameters and discuss whether they can be approximated as constant values. Consequently, we derive the $\alpha$-$\beta$ scaling as a function of orbital frequencies and demonstrate that it is equivalent to a frequency-dependent correction. We further provide a brief comparison between Post-Newtonian waveform and the rescaled ppBHPT waveform at $q=3$ and comment on their regime of validity. Finally, we explore the possibility of using PN to obtain the $\alpha$-$\beta$ calibration parameters and still provide a rescaled ppBHPT waveform that matches NR.

The symmetry can be broken at high temperature and then restored at low temperature, which is the so-called \emph{high temperature symmetry breaking}. It often appears in some theories such as the high scale electroweak baryogenesis mechanism. In this paper, we probe the high temperature $\mathbb{Z}_2$ symmetry breaking with gravitational waves (GWs) from domain wall annihilation. We first introduce a scalar with $\mathbb{Z}_2$ symmetry and few of singlet fermions that interact with scalar through a five-dimension operator. This can lead to the scalar potential has a non-zero minimum at high temperature. At the early stage, the scalar is pinned at symmetric phase due to the large Hubble fraction. When the scalar thermal mass becomes comparable to the Hubble parameter, it can quickly roll down to the minimum of potential. Then the $\mathbb{Z}_2$ symmetry is spontaneously broken and the domain walls will form. With the decrease of temperature, $\mathbb{Z}_2$ symmetry will be restored. We find that if domain walls are formed at $\mathcal{O}(10^{9})~ \rm GeV$, the GW produced by domain wall annihilation is expected to be observed by BBO, CE and ET. In addition, we also discuss the relationships between this scenario and NANOGrav signal.

The stochastic gravitational wave background (SGWB) detected recently by the pulsar timing arrays (PTAs) observations may have cosmological origins. In this work we consider a model of single field inflation containing an intermediate phase of ultra slow-roll. Fixing the amplitude of the peak of curvature perturbations by the PBHs bounds we calculate the gravitational waves (GWs) induced from the curvature perturbations enhanced during USR. The spectrum of the induced GWs depends on the sharpness of the transition from the USR phase to the final attractor phase as well as to the duration of the USR period. While the model can accommodate the current PTAs data but it has non-trivial predictions for the induced GWs on higher frequency ranges which can be tested by future observations.

In the light of the evidence of a gravitational wave background from the NANOGrav 15yr data set, we reconsider the split majoron model as a new physics extension of the standard model able to generate a needed contribution to solve the current tension between the data and the standard interpretation in terms of inspiraling supermassive black hole massive binaries. In the split majoron model the seesaw right-handed neutrinos acquire Majorana masses from spontaneous symmetry breaking of global $U(1)_{B-L}$ in a strong first order phase transition of a complex scalar field occurring above the electroweak scale. The final vacuum expectation value couples to a second complex scalar field undergoing a low scale phase transition occurring after neutrino decoupling. Such a coupling enhances the strength of this second low scale first order phase transition and can generate a sizeable primordial gravitational wave background contributing to the NANOGrav 15yr signal. Moreover, the free streaming length of light neutrinos can be suppressed by their interactions with the resulting Majoron background and this can mildly ameliorate existing cosmological tensions, thus providing a completely independent motivation for the model.