Papers for Monday, Aug 01 2022

The tidal disruption of stars by supermassive black holes (SMBHs) probes relativistic gravity. In the coming decade, the number of observed tidal disruption events (TDEs) will grow by several orders of magnitude, allowing statistical inferences of the properties of the SMBH and stellar populations. Here we analyse the probability distribution functions of the pericentre distances of stars that encounter an SMBH in the Schwarzschild geometry, where the results are completely analytic, and the Kerr metric. From this analysis we calculate the number of observable TDEs, defined to be those that come within the tidal radius $r_{\rm t}$ but outside the direct capture radius (which is, in general, larger than the horizon radius). We find that relativistic effects result in a steep decline in the number of stars that have pericenter distances $r_{\rm p} \lesssim 10\,r_{\rm g}$, where $r_{\rm g} = GM/c^2$, and that for maximally spinning SMBHs the distribution function of $r_{\rm p}$ at such distances scales as $f_{\rm r_{\rm p}}\propto r_{\rm p}^{4/3}$, or in terms of $\beta \equiv r_{\rm t}/r_{\rm p}$ scales as $f_{\beta} \propto \beta^{-10/3}$. We find that spin has little effect on the TDE fraction until the very high-mass end, where instead of being identically zero the rate is small ($\lesssim 1\%$ of the expected rate in the absence of relativistic effects). Effectively independent of spin, if the progenitors of TDEs reflect the predominantly low-mass stellar population and thus have masses $\lesssim 1M_{\odot}$, we expect a substantial reduction in the rate of TDEs above $10^{7}M_{\odot}$.

We model satellite quenching at $z \sim 1$ by combining $14$ massive ($10^{13.8} < M_{\mathrm{halo}}/\mathrm{M}_{\odot} < 10^{15}$) clusters at $0.8 < z < 1.3$ from the GOGREEN and GCLASS surveys with accretion histories of $56$ redshift-matched analogs from the IllustrisTNG simulation. Our fiducial model, which is parameterized by the satellite quenching timescale ($\tau_{\rm quench}$), accounts for quenching in our simulated satellite population both at the time of infall by using the observed coeval field quenched fraction and after infall by tuning $\tau_{\rm quench}$ to reproduce the observed satellite quenched fraction versus stellar mass trend. This model successfully reproduces the observed satellite quenched fraction as a function of stellar mass (by construction), projected cluster-centric radius, and redshift and is consistent with the observed field and cluster stellar mass functions at $z \sim 1$. We find that the satellite quenching timescale is mass dependent, in conflict with some previous studies at low and intermediate redshift. Over the stellar mass range probed ($M_{\star}> 10^{10}~\mathrm{M}_{\odot}$), we find that the satellite quenching timescale decreases with increasing satellite stellar mass from $\sim1.6~{\rm Gyr}$ at $10^{10}~\mathrm{M}_{\odot}$ to $\sim 0.6 - 1~{\rm Gyr}$ at $10^{11}~\mathrm{M}_{\odot}$ and is roughly consistent with the total cold gas (H{\scriptsize I}+H$_{2}$) depletion timescales at intermediate $z$, suggesting that starvation may be the dominant driver of environmental quenching at $z < 2$. Finally, while environmental mechanisms are relatively efficient at quenching massive satellites, we find that the majority ($\sim65-80\%$) of ultra-massive satellites ($M_{\star} > 10^{11}~\mathrm{M}_{\odot}$) are quenched prior to infall.

We investigate the effects of stellar populations and sizes on Ly$\alpha$ escape in 27 spectroscopically confirmed and 35 photometric Lyman-Alpha Emitters (LAEs) at z $\approx$ 2.65 in seven fields of the Bo\"otes region of the NOAO Deep Wide-Field Survey. We use deep $HST$/WFC3 imaging to supplement ground-based observations and infer key galaxy properties. Compared to typical star-forming galaxies (SFGs) at similar redshifts, the LAEs are less massive ($M_{\star} \approx 10^{7} - 10^{9}~M_{\odot}$), younger (ages $\lesssim$ 1 Gyr), smaller ($r_{e} <$ 1 kpc), less dust-attenuated (E(B$-$V) $\le$ 0.26 mag), but have comparable star-formation-rates (SFRs $\approx 1 - 100~M_{\odot} {\rm yr^{-1}}$). Some of the LAEs in the sample may be very young galaxies having low nebular metallicities (${\rm Z_{neb} \lesssim 0.2 Z_{\odot}}$) and/or high ionization parameters ($\log{(\rm U)} \gtrsim -2.4$). Motivated by previous studies, we examine the effects of the concentration of star formation and gravitational potential on Ly$\alpha$ escape, by computing star-formation-rate surface density, $\Sigma_{\rm SFR}$ and specific star-formation-rate surface density, $\Sigma_{\rm sSFR}$. For a given $\Sigma_{\rm SFR}$, the Ly$\alpha$ escape fraction is higher for LAEs with lower stellar masses. LAEs have higher $\Sigma_{\rm sSFR}$ on average compared to SFGs. Our results suggest that compact star formation in a low gravitational potential yields conditions amenable to the escape of Ly$\alpha$ photons. These results have important implications for the physics of Ly$\alpha$ radiative transfer and for the type of galaxies that may contribute significantly to cosmic reionization.

The Argus Optical Array is a synoptic survey observatory, currently in development, that will have a total collecting area equivalent to a 5-meter monolithic telescope and an all-sky field of view, multiplexed from 900 commercial off-the-shelf telescopes. The Array will observe 7916 deg$^2$ every second during high-speed operations ($m_g\leq16.1$) and every 30 seconds at base cadence ($m_g\leq19.1$), producing 4.3 PB and 145 TB respectively of data per night with its 55-gigapixel mosaic of cameras. The Argus Array Hierarchical Data Processing System (Argus-HDPS) is the instrument control and analysis pipeline for the Argus Array project, able to create fully-reduced data products in real time. We pair sub-arrays of cameras with co-located compute nodes, responsible for distilling the raw 11 Tbps data rate into transient alerts, full-resolution image segments around selected targets at 30-second cadence, and full-resolution coadds of the entire field of view at $15+$-min cadences. Production of long-term light curves and transient discovery in deep coadds out to 5-day cadence ($m_g\leq24.0$) will be scheduled for daytime operations. In this paper, we describe the data reduction strategy for the Argus Optical Array and demonstrate image segmentation, coaddition, and difference image analysis using the GPU-enabled Argus-HDPS pipelines on representative data from the Argus Array Technology Demonstrator.

We test the assumption of entropy conservation between Big Bang nucleosynthesis (BBN) and recombination by considering a massive particle that decays into a mixture of photons and other relativistic species. We employ Planck temperature and polarization anisotropies, COBE/FIRAS spectral distortion bounds, and the observed primordial deuterium abundance to constrain these decay scenarios. If between $56\%$ and $71\%$ of the decaying particle's energy is transferred to photons, then $N_{\mathrm{eff}}$ at recombination is minimally altered, and Planck data alone allows for significant entropy injection. If photons are injected by the decay, the addition of spectral distortion bounds restricts the decay rate of the particle to be $\Gamma_Y > 1.91\times10^{-6} \text{s}^{-1}$ at $95\%$ confidence level (CL). We find that constraints on the energy density of the decaying particle are significantly enhanced by the inclusion of bounds on the primordial deuterium abundance, allowing the particle to contribute at most $2.35\%$ ($95\%$ CL) of the energy density of the universe before decaying.

Galaxy mergers at high redshifts trigger the activity of their central supermassive black holes, eventually also leading to their coalescence -- and a potential source of low-frequency gravitational waves detectable by the SKA Pulsar Timing Array (PTA). Two key parameters related to the fuelling of black holes are the Eddington ratio of quasar accretion $\eta_{\rm Edd}$, and the radiative efficiency of the accretion process, $\epsilon$ (which affects the so-called active lifetime of the quasar, $t_{\rm QSO}$). We forecast the regime of detectability of gravitational wave events with SKA PTA, finding the associated binaries to have orbital periods on the order of weeks to years, observable through relativistic Doppler velocity boosting and/or optical variability of their light curves. Combining the SKA regime of detectability with the latest observational constraints on high-redshift black hole mass and luminosity functions, and theoretically motivated prescriptions for the merger rates of dark matter haloes, we forecast the number of active counterparts of SKA PTA events expected as a function of primary black hole mass at $z \gtrsim 6$. We find that the quasar counterpart of the most massive black holes will be ${uniquely \ localizable}$ within the SKA PTA error ellipse at $z \gtrsim 6$. We also forecast the number of expected counterparts as a function of the quasars' Eddington ratio and active lifetime. Our results show that SKA PTA detections can place robust constraints on the seeding and growth mechanisms of the first supermassive black holes.

\emph{Wide-field Infrared Survey Explorer} all-sky survey has discovered a new population of hot dust-obscured galaxies (Hot DOGs), which has been confirmed to be dusty quasars. Previous statistical studies have found significant overdensities of sub-millimeter and mid-IR selected galaxies around Hot DOGs, indicating they may reside in dense regions. Here we present the near-infrared ($J$ and $K_s$ bands) observations over a $7.5'\times 7.5'$ field centered on a Hot DOG W1835$+$4355 at $z \sim 2.3$ using the wide-field infrared camera on the Palomar 200-inch telescope. We use the color criterion $J-K_s>2.3$ for objects with $K_s<20$, to select Distant Red Galaxies (DRGs). We find a significant excess of number density of DRGs in W1835$+$4355 field compared to three control fields, by a factor of about 2. The overdensity of red galaxies around W1835$+$4355 are consistent with the multi-wavelength environment of Hot DOGs, suggesting that Hot DOGs may be a good tracer for dense regions at high redshift. We find that W1835$+$4355 do not reside in the densest region of the dense environment traced by itself. A possible scenario is that W1835$+$4355 is undergoing merging process, which lowers the local number density of galaxies in its surrounding region.

Fast radio bursts (FRBs) are brief, energetic, extragalactic flashes of radio emission whose progenitors are largely unknown. Although studying the FRB population is essential for understanding how these astrophysical phenomena occur, such studies have been difficult to conduct without large numbers of FRBs and characterizable observational biases. Using the recently released catalog of 536 FRBs published by the Canadian Hydrogen Intensity Mapping Experiment/Fast Radio Burst (CHIME/FRB) collaboration, we present a study of the FRB population that also calibrates for selection effects. Assuming a Schechter luminosity function, we infer a characteristic energy cut-off of $E_\mathrm{char} =$ $2.38^{+5.35}_{-1.64} \times 10^{41}$ erg and a differential power-law index of $\gamma =$ $-1.3^{+0.7}_{-0.4}$. Simultaneously, we infer a volumetric rate of [$7.3^{+8.8}_{-3.8}$(stat.)$^{+2.0}_{-1.8}$(sys.)]$\times 10^4$ bursts Gpc$^{-3}$ year$^{-1}$ above a pivot energy of 10$^{39}$ erg and below a scattering timescale of 10 ms at 600 MHz, and find we cannot significantly constrain the cosmic evolution of the FRB population with star formation rate. Modeling the host dispersion measure (DM) contribution as a log-normal distribution and assuming a total Galactic contribution of 80 pc cm$^{-3}$, we find a median value of $\mathrm{DM}_\mathrm{host} =$ $84^{+69}_{-49}$ pc cm$^{-3}$, comparable with values typically used in the literature. Proposed models for FRB progenitors should be consistent with the energetics and abundances of the full FRB population predicted by our results. Finally, we infer the redshift distribution of FRBs detected with CHIME, which will be tested with the localizations and redshifts enabled by the upcoming CHIME/FRB Outriggers project.

The Argus Optical Array will be the first all-sky, arcsecond-resolution, 5-m class telescope. The 55 GPix Array, currently being prototyped, will consist of 900 telescopes with 61 MPix very-low-noise CMOS detectors enabling sub-second cadences. Argus will observe every part of the northern sky for 6-12 hours per night, achieving a simultaneously high-cadence and deep-sky survey. The array will build a two-color, million-epoch movie, reaching dark-sky depths of $m_g$=19.6 each minute and $m_g$=23.6 each week over 47% of the entire sky, enabling the most-sensitive-yet searches for high-speed transients, gravitational-wave counterparts, exoplanet microlensing events, and a host of other phenomena. In this paper we present our newly-developed array arrangement, which mounts all telescopes into the inside of a hemispherical bowl (turning the original dome design inside-out). The telescopes' beams thus converge at a single ``pseudofocal'' point. When placed along the telescope's polar axis, this point does not move as the telescope tracks, allowing every telescope to simultaneously look through a single, unmoving window in a fixed enclosure. This telescope bowl is suspended from a simple free-swinging pivot (turning the usual telescope mounting support upside-down), with polar alignment afforded by the creation of a virtual polar axis defined by a second mounting pivot. This new design, currently being prototyped with the 38-telescope Argus Pathfinder, eliminates the need for a movable external dome and thus greatly reduces the cost and complexity of the full Argus Array. Coupled with careful software scope control and the use of existing software pipelines, the Argus Array could thus become one of the deepest and fastest sky surveys, within a midscale-level budget.

Here we present a statistical study on tidal features around massive early-type galaxies (ETGs). Utilizing the imaging data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), we measure the flux fraction of tidal features ($f_{\rm tidal}$) in 2649 ETGs with stellar mass $M_*>10^{11}M_{\odot}$ and redshift $0.05<z<0.15$ using automated techniques. The Wide-layer of HSC-SSP reaches a depth of $\sim 28.5$ mag arcsec$^{-2}$ in $i$-band. Under this surface brightness limit, we find that about 28 % of these galaxies harbor prominent tidal features with $f_{\rm tidal}>1\%$, among which the number of ETGs decreases exponentially with $f_{\rm tidal}$, with a logarithmic slope of $\sim100$. Within the stellar mass range we probe, we note that $f_{\rm tidal}$ increases by a factor of 2 from $M_*\approx10^{11}M_{\odot}$ to $M_*\approx10^{12}M_{\odot}$. We also perform pair-count to estimate the merger rate of these massive ETGs. Combining the merger rates with $f_{\rm tidal}$, we estimate that the typical lifetime of tidal features is $\sim$ 3 Gyr, consistent with previous studies.

We use Molecular Dynamics simulations to study the formation and stability of single and multicomponent lattices in the outer crust of Neutron Stars. Including an improved treatment for Gaussian charge distribution of ions we obtain the expressions for the potential and forces for electron screened Coulomb interactions using the efficient Ewald sum procedure. Our findings show that for baryon densities in the outer crust a point-like ion treatment can not fully describe the crystallization behaviour thus the Coulomb theory melting parameter, $\Gamma_C$ , must be complemented with an additional parameter, $\eta$, providing information on the finite size of ions. In our approach we find that including beyond point-like approaches in screened plasmas has a robust impact on calculated lattice energetic stability decreasing crystallization energies per baryon up to $\sim$40% with respect to point-like interaction and, as a consequence, melting point resulting displaced to lower temperatures.

The eROSITA X-ray telescope, launched in 2019, is predicted to observe roughly 100,000 galaxy clusters. Follow-up observations of these clusters from Chandra, for example, will be needed to resolve outstanding questions about galaxy cluster physics. Deep Chandra cluster observations are expensive and follow-up of every eROSITA cluster is infeasible, therefore, objects chosen for follow-up must be chosen with care. To address this, we have developed an algorithm for predicting longer duration, background-free observations based on mock eROSITA observations. We make use of the hydrodynamic cosmological simulation Magneticum, have simulated eROSITA instrument conditions using SIXTE, and have applied a novel convolutional neural network to output a deep Chandra-like "super observation'' of each cluster in our simulation sample. Any follow-up merit assessment tool should be designed with a specific use case in mind; our model produces observations that accurately and precisely reproduce the cluster morphology, which is a critical ingredient for determining cluster dynamical state and core type. Our model will advance our understanding of galaxy clusters by improving follow-up selection and demonstrates that image-to-image deep learning algorithms are a viable method for simulating realistic follow-up observations.

Molecular gas flows are analyzed in 14 cluster galaxies (BCGs) centered in cooling hot atmospheres. The BCGs contain $10^{9}-10^{11}~\rm M_\odot$ of molecular gas, much of which is being moved by radio jets and lobes. The molecular flows and radio jet powers are compared to molecular outflows in 45 active galaxies within $z<0.2$. We seek to understand the relative efficacy of radio, quasar, and starburst feedback over a range of active galaxy types. Molecular flows powered by radio feedback in BCGs are $\sim$10--1000 times larger in extent compared to contemporary galaxies hosting quasar nuclei and starbursts. Radio feedback yields lower flow velocities but higher momenta compared to quasar nuclei, as the molecular gas flows in BCGs are usually $\sim$10--100 times more massive. The product of the molecular gas mass and lifting altitude divided by the AGN or starburst power -- a parameter referred to as the lifting factor -- exceeds starbursts and quasar nuclei by two to three orders of magnitude, respectively. When active, radio feedback is generally more effective at lifting gas in galaxies compared to quasars and starburst winds. The kinetic energy flux of molecular clouds generally lies below and often substantially below a few percent of the driving power. We find tentatively that star formation is suppressed in BCGs relative to other active galaxies, perhaps because these systems rarely form molecular disks that are more impervious to feedback and are better able to promote star formation.

Ionizing photons must have escaped from high-redshift galaxies, but the neutral high-redshift intergalactic medium makes it unlikely to directly detect these photons during the Epoch of Reionization. Indirect methods of studying ionizing photon escape fractions present a way to infer how the first galaxies may have reionized the universe. Here, we use HET/LRS2 observations of J0919+4906, a confirmed z$\approx$0.4 emitter of ionizing photons to achieve spatially resolved (12.5 kpc in diameter) spectroscopy of \ion{Mg}{II}$\lambda2796$, \ion{Mg}{II}$\lambda2803$, [\ion{O}{II}]$\lambda\lambda3727,3729$ , [\ion{Ne}{III}]$\lambda3869$, H$\gamma$, [\ion{O}{III}]$\lambda4363$, H$\beta$, [\ion{O}{III}]$\lambda4959$, [\ion{O}{III}]$\lambda5007$, and H$\alpha$. From these data we measure \ion{Mg}{II} emission, which is a promising indirect tracer of ionizing photons, along with nebular ionization and dust attenuation in multiple spatially-resolved apertures. We find that J0919+4906 has significant spatial variation in its \ion{Mg}{II} escape and thus ionizing photon escape fraction. Combining our observations with photoionization models, we find that the regions with the largest relative \ion{Mg}{II} emission and \ion{Mg}{II} escape fractions have the highest ionization and lowest dust attenuation. Some regions have an escape fraction that matches that required by models to reionize the early universe, while other regions do not. We observe a factor of 36 spatial variation in the inferred LyC escape fraction, which is similar to recently observed statistical samples of indirect tracers of ionizing photon escape fractions. These observations suggest that spatial variations in neutral gas properties lead to large variations in the measured LyC escape fractions. Our results suggest that single sightline observations may not trace the volume-averaged escape fraction of ionizing photons.

We consider a generic model of quadratic gravity coupled to a single scalar and investigate the effects of gravitational degrees of freedom on inflationary parameters. We find that quantum corrections arising from the massive spin-2 ghost generate significant contributions to the effective inflationary potential and allow for a realization of the spontaneous breakdown of global scale invariance without the need for additional scalar fields. We compute inflationary parameters, compare the resulting predictions to well-known inflationary models, and find that they fit well within the Planck collaboration's constraints on inflation.

Large amounts of Adaptive-Optics (AO) control loop data and telemetry are currently inaccessible to end-users. Broadening access to those data has the potential to change the AO landscape on many fronts, addressing several use-cases such as derivation of the system's PSF, turbulence characterisation and optimisation of system control. We address one of the biggest obstacles to sharing these data: the lack of standardisation, which hinders access. We propose an object-oriented Python package for AO telemetry, whose data model abstracts the user from an underlining archive-ready data exchange standard based on the Flexible Image Transport System (FITS). Its design supports data from a wide range of existing and future AO systems, either in raw format or abstracted from actual instrument details. We exemplify its usage with data from active AO systems on 10m-class observatories, of which two are currently supported (AOF and Keck), with plans for more.

UV wind line variability in OB stars appears to be universal. We review the evidence that the variability is due to large, dense, optically thick structures rooted in or near the photosphere. Using repeated bservations and a simple model we translate observed profile variations into optical depth variations and, consequently, variations in measured mass loss rates. Although global rates may be stable, measured rates vary. Consequently, profile variations infer how mass loss rates determined from UV wind lines vary. These variations quantify the intrinsic error inherent in any mass loss rate derived from a single observation. These derived rates can differ by factors of 3 or more. Our results also imply that rates from non-simultaneous observations (such as UV and ground based data) need not agree. Finally, we use our results to examine the nature of the structures responsible for the variability.

The {\em Planck} DR3 measurements of the temperature and polarization anisotropies power spectra of the cosmic microwave background (CMB) show an excess of smoothing of the acoustic peaks with respect to $\Lambda$CDM, often quantified by a phenomenological parameter $A_{\rm L}$. A specific feature superimposed to the primordial power spectrum has been suggested as a physical solution for this smoothing excess. Here, we investigate the impact of this specific localized oscillation with a frequency linear in the wavenumber, designed to mimic the smoothing of CMB temperature spectrum corresponding to $A_{\rm L} \simeq 1.1-1.2$ on the matter power spectrum. We verify the goodness of the predictions in perturbation theory at next-to-leading order with a set of N-body simulations, a necessary step to study the non-linear damping of these primordial oscillations. We show that for a large portion of the parameter space, the amplitude of this primordial oscillation can be strongly damped on the observed nonlinear matter power spectrum at $z=0$, but a larger signal is still persistent at $z \lesssim 2$ and is therefore a target for future galaxy surveys at high redshifts. From an analysis of the BOSS DR12 two-point correlation function, we find ${\cal A}_{\rm lin} < 0.26$ at 95\% CL by keeping the frequency fixed to the best-fit of {\em Planck} data.

The development of the Keck All sky Precision Adaptive optics (KAPA) project was initiated in September 2018 to upgrade the Keck I adaptive optics (AO) system to enable laser tomography adaptive optics (LTAO) with a four laser guide star (LGS) asterism. The project includes the replacement of the existing LMCT laser with a Toptica laser, the implementation of a new real-time controller (RTC) and wavefront sensor optics and camera, and a new daytime calibration and test platform to provide the required infrastructure for laser tomography. The work presented here describes the new daytime calibration infrastructure to test the performance for the KAPA tomographic algorithms. This paper outlines the hardware infrastructure for daytime calibration and performance assessment of tomographic algorithms. This includes the implementation of an asterism simulator having fiber-coupled light sources simulating four Laser Guide Stars (LGS) and two Natural Guide Stars (NGS) at the AO bench focus, as well as the upgrade of the existing TelSim on the AO bench to simulate focal anisoplanatism and wind driven atmospheric turbulence. A phase screen, that can be adjusted in effective altitude, is used to simulate wind speeds up to 10 m/s for a duration of upto 3 s.

It is an open question whether and how gravitational wave events involving neutron stars can be preceded by electromagnetic counterparts. This work shows that the collision of two neutron stars with magnetic fields well below magnetar-level strengths can produce millisecond Fast-Radio-Burst-like transients. Using global force-free electrodynamics simulations, we demonstrate that electromagnetic flares, produced by overtwisted common flux tubes in the binary magnetosphere, collide with the orbital current sheet and compress it, resulting in enhanced magnetic reconnection. As a result, the current sheet fragments into a sequence of plasmoids, which collide with each other leading to the emission of coherent electromagnetic waves. The resulting millisecond-long burst of radiation should have frequencies in the range of $10-20\,\rm GHz$ for magnetic fields of $B^{\ast}=10^{11}\, \rm G$ at the stellar surfaces.

Numerical simulation of hydrodynamic equations forms the central part of solving various modern astrophysical problems. In the case of shocks, one can have either dynamical equations or jump conditions (the conservation equations without any time evolution). The solution of the jump condition in curve space-time is derived and analyzed in detail in the present work. We also derive the Taub adiabat or combustion adiabat equation from the jump condition. We have analyzed both time-like and space-like shocks in the present work. We find that the change in entropy for the weak shocks for curved space-time is small similar to that for flat space-time. We also find that for general relativistic space-like shocks, the Chapman-Jouguet point does not necessarily correspond to the sonic point for downstream matter, unlike the relativistic case. To analyze the shock wave solution for the curved space-time, one needs the information of metric potentials describing the space-time, which for the present work is taken to be a neutron star. We assume that a shock wave is generated at the centre of the star and is propagating outward. As the shock wave is propagating outwards, it combusts nuclear matter to quark matter, and we have a combustion scenario. We find that the general relativistic treatment of shock conditions is necessary to study shocks in neutron stars so that the results are consistent with the solution of the TOV equation while calculating the maximum mass for a given equation of state. We also find that with such general relativistic treatment, the combustion process in neutron stars is always a detonation.

We have performed simulations of equal and unequal mass binary neutron star mergers in this study. We have compared how the observational gravitational signals change depending upon the equation of state, which governs matter at high density. Mainly we have compared results between the hadronic and quark equation of state. The quark matter is modelled with Gibbs formalism, where a mixed-phase appears between the pure hadronic and pure quark state. We have varied the density where quark matter first appears. It is found that when two almost equal mass binaries merge, the final star attains a stable configuration (does not collapse) irrespective of the matter properties. However, if the post-merger scenario of the binary merger can be probed, a significant difference in the gravitational wave amplitude can be seen if quark matter appears at low densities. It is found that if the matter properties with hadronic and quark degrees differ significantly, it is reflected in the stability of the final merger product. Hadronic matter can give a stable remnant star, whereas if quark matter appears, it makes the collapse possible (assuming that the hadronic equation of state is stiffer than the quark equation of state). However, when unequal mass binaries (the mass difference is significant) merge, the difference in the observational signals depending on the equation of state (hadronic or quark) is evident just from the point of first contact. It is also seen that the difference in gravitational signal (depending on the equation of state) is more significant for unequal mass binary merging than for equal mass binary having the same total baryonic mass.

Some cosmological models with non-negligible dark energy fractions in particular windows of the pre-recombination epoch are capable of alleviating the Hubble tension quite efficiently, while keeping the good description of the data that are used to build the cosmic inverse distance ladder. There has been an intensive discussion in the community on whether these models enhance the power of matter fluctuations, leading de facto to a worsening of the tension with the large-scale structure measurements. We address this pivotal question in the context of several early dark energy (EDE) models, considering also in some cases a coupling between dark energy and dark matter, and the effect of massive neutrinos. We fit them using the Planck 2018 likelihoods, the supernovae of Type Ia from the Pantheon compilation and data on baryon acoustic oscillations. We find that ultra-light axion-like (ULA) EDE can actually alleviate the $H_0$ tension without increasing the values of $\sigma_{12}$ with respect to those found in the $\Lambda$CDM, whereas EDE with an exponential potential does not have any impact on the tensions. A coupling in the dark sector tends to enhance the clustering of matter, and the data limit a lot the influence of massive neutrinos, since the upper bounds on the sum of their masses are too close to those obtained in the standard model. We find that in the best case, namely ULA, the Hubble tension is reduced to $\sim 2\sigma$.

Many spectacular polarimetric images have been obtained in recent years with adaptive optics (AO) instruments at large telescopes because they profit significantly from the high spatial resolution. This paper summarizes some basic principles for AO polarimetry, discusses challenges and limitations of these systems, and describes results which illustrate the performance of AO polarimeters for the investigation of circumstellar disks, of dusty winds from evolved stars, and for the search of reflecting extra-solar planets.

The non-singular bouncing cosmology is an alternative paradigm to inflation, wherein the background energy density vanishes at the bounce point, in the context of Einstein gravity. Therefore, the non-linear effects in the evolution of density fluctuations ($\delta \rho$) may be strong in the bounce phase, which potentially provides a mechanism to enhance the abundance of primordial black holes (PBHs). This article presents a comprehensive illustration for PBH enhancement due to the bounce phase. To calculate the non-linear evolution of $\delta \rho$, the Raychaudhuri equation is numerically solved here. Since the non-linear processes may lead to a non-Gaussian probability distribution function for $\delta \rho$ after the bounce point, the PBH abundance is calculated in a modified Press-Schechter formalism. In this case, the criterion of PBH formation is complicated, due to complicated non-linear evolutionary behavior of $\delta \rho$ during the bounce phase. Our results indicate that the bounce phase indeed has potential to enhance the PBH abundance sufficiently. Furthermore, the PBH abundance is applied to constrain the parameters of bounce phase, providing a complementary to the surveys of cosmic microwave background and large scale structure.

In a $\Lambda$CDM universe, most galaxies evolve by mergers and accretions, leaving faint and/or diffuse structures, such as tidal streams and stellar halos. Although these structures are a good indicator of galaxies' recent mass assembly history, they have the disadvantage of being difficult to observe due to their low surface brightness (LSB). To recover these LSB features by minimizing the photometric uncertainties introduced by the optical system, we developed a new optimized telescope named K-DRIFT pathfinder, adopting a linear astigmatism free-three mirror system. Thanks to the off-axis design, it is expected to avoid the loss and scattering of light on the optical path within the telescope. To assess the performance of this prototype telescope, we investigate the photometric depth and capability to identify LSB features. We find that the surface brightness limit reaches down to $\mu_{r,1\sigma}\sim28.5$ mag arcsec$^{-2}$ in $10^{\prime\prime}\times10^{\prime\prime}$ boxes, enabling us to identify a single stellar stream to the east of NGC 5907. We also examine the characteristics of the point spread function (PSF) and find that the PSF wing reaches a very low level. Still, however, some internal reflections appear within a radius of $\sim$6 arcmin from the center of sources. Despite a relatively small aperture (0.3 m) and short integration time (2 hr), this result demonstrates that our telescope is highly efficient in LSB detection.

We investigated the solar cycle dependency on the presence and periodicity of the Quasi-Biennial Oscillation (QBO). Using helioseismic techniques, we used solar oscillation frequencies from the Global Oscillations Network Group (GONG), Michelson Doppler Imager (MDI) and Helioseismic & Magnetic Imager (HMI) in the intermediate-degree range to investigate the frequency shifts over Cycles 23 and 24. We also examined two solar activity proxies, the F10.7 index and the MgII index, for the last four solar cycles to study the associated QBO. The analyses were performed using Empirical Mode Decomposition (EMD) and the Fast Fourier Transform (FFT). We found that the EMD analysis method is susceptible to detecting statistically significant Intrinsic Mode Functions (IMFs) with periodicities that are overtones of the length of the dataset under examination. Statistically significant periodicites, which were not due to overtones, were detected in the QBO range. We see a reduced presence of the QBO in Cycle 24 compared to Cycle 23. The presence of the QBO was not sensitive to the depth to which the p-mode travelled, nor the average frequency of the p-mode. The analysis further suggested that the magnetic field responsible for producing the QBO in frequency shifts of p-modes is anchored above approximately 0.95 solar radii.

The cross-correlation of optical galaxies with the neutral hydrogen (HI) radiation intensity can enhance the signal-to-noise ratio (SNR) of the HI intensity measurement. In this paper, we investigate the cross-correlation of the galaxy samples obtained by the spectroscopic survey of the China Space Station Telescope (CSST) with the HI Intensity mapping (IM) survey of the Five-hundred-meter Aperture Spherical Telescope (FAST). Using the IllusitrisTNG simulation result at redshift $0.2 \sim 0.3$, we generate mock data of the CSST survey and a FAST L-band drift scan survey. The CSST spectroscopic survey can yield a sample of galaxies with a high comoving number density of $10^{-2} (\unit{Mpc}/h)^{-3}$ at $z \sim 0.3$. We cross-correlate the foreground-removed radio intensity with the CSST galaxies, including both the whole sample, and red and blue galaxy sub-samples separately. We find that in all cases the HI and optical galaxies are well correlated. The total HI abundance can be measured with a high precision from this correlation. A relative error of $\sim 0.6\%$ for $\Omega_{\rm HI}$ could be achieved at $z\sim 0.3$ for an overlapping survey area of $10000 \unit{deg}^2$.

We present N-PDFs of 29 Galactic regions obtained from Herschel imaging at high angular resolution, covering diffuse and quiescent clouds, and those showing low-, intermediate-, and high-mass star formation (SF), and characterize the cloud structure using the Delta-variance tool. The N-PDFs are double-log-normal at low column densities, and display one or two power law tails (PLTs) at higher column densities. For diffuse, quiescent, and low-mass SF clouds, we propose that the two log-normals arise from the atomic and molecular phase, respectively. For massive clouds, we suggest that the first log-normal is built up by turbulently mixed H2 and the second one by compressed (via stellar feedback) molecular gas. Nearly all clouds have two PLTs with slopes consistent with self-gravity, where the second one can be flatter or steeper than the first one. A flatter PLT could be caused by stellar feedback or other physical processes that slow down collapse and reduce the flow of mass toward higher densities. The steeper slope could arise if the magnetic field is oriented perpendicular to the LOS column density distribution. The first deviation point (DP), where the N-PDF turns from log-normal into a PLT, shows a clustering around values of a visual extinction of AV (DP1) around 2-5. The second DP, which defines the break between the two PLTs, varies strongly. Using the Delta-variance, we observe that the AV value, where the slope changes between the first and second PLT, increases with the characteristic size scale in the variance spectrum. We conclude that at low column densities, atomic and molecular gas is turbulently mixed, while at high column densities, the gas is fully molecular and dominated by self-gravity. The best fitting model N-PDFs of molecular clouds is thus one with log-normal low column density distributions, followed by one or two PLTs.

A laboratory plasma-wall interaction-based astrophysical gravito-electrostatic sheath (GES) model is methodologically applied to study the dynamic stability of the magnetoactive bi-fluidic solar plasma system in the presence of turbulence effect. The spherically symmetric GES-model formalism couples the solar interior plasma (SIP, internally self-gravitating, bounded) and the solar wind plasma (SWP, externally point-gravitating, unbounded) through the diffused solar surface boundary (SSB). A normal spherical mode ansatz results in a generalized linear quadratic dispersion relation depicting the modal fluctuations on both the SIP and SWP scales. A constructive numerical platform reveals the evolution of both dispersive and non-dispersive modal features of the modified-GES mode excitations. The reliability of the derived non-planar dispersion laws is concretized with the help of an exact analytic shape matching the previously reported results founded on the plane-wave approximation. It is found that the thermo-statistical GES stability depends mainly on the magnetic field, equilibrium plasma density, and plasma temperature. It is speculated that the dispersive features are more pronounced in the self-gravitational domains against the electrostatic ones. The magneto-thermal interplay introduces decelerating (accelerating) and destabilizing (stabilizing) influences on the SIP (SWP), and so forth. At last, we briefly indicate the applicability of the proposed analysis to understand diverse helioseismic activities from the collective plasma dynamical viewpoint in accordance with the recent astronomical observational scenarios reported in the literature.

Spiral density waves observed in protoplanetary discs have often been used to infer the presence of embedded planets. This inference relies both on simulations as well as the linear theory of planet-disc interaction developed for planets on circular orbits to predict the morphology of the density wake. In this work we develop and implement a linear framework for calculating the structure of the density wave in a gaseous disc driven by an eccentric planet. Our approach takes into account both the essential azimuthal and temporal periodicities of the problem, allowing us to treat any periodic perturbing potential (i.e. not only that of an eccentric planet). We test our framework by calculating the morphology of the density waves excited by an eccentric, low-mass planet embedded in a globally isothermal disc and compare our results to the recent direct numerical simulations (and heuristic wavelet analysis) of the same problem by Zhu and Zhang. We find excellent agreement with the numerical simulations, capturing all the complex eccentric features including spiral bifurcations, wave crossings and planet-wave detachments, with improved accuracy and detail compared with the wavelet method. This illustrates the power of our linear framework in reproducing the morphology of complicated time-dependent density wakes, presenting it as a valuable tool for future studies of eccentric planet-disc interactions.

Current searches for gravitational waves from compact-binary objects are primarily designed to detect the dominant gravitational-wave mode and assume that the binary components have spins which are aligned with the orbital angular momentum. These choices lead to observational biases in the observed distribution of sources. Sources with significant spin-orbit precession or unequal-mass-ratios, which have non-negligible contributions from sub-dominant gravitational-wave modes, may be missed; in particular, this may significantly suppress or bias the observed neutron star -- black hole (NSBH) population. We simulate a fiducial population of NSBH mergers and determine the impact of using searches that only account for the dominant-mode and aligned spin. We compare the impact for the Advanced LIGO design, A+, LIGO Voyager, and Cosmic Explorer observatories. We find that for a fiducial population where the spin distribution is isotropic in orientation and uniform in magnitude, we will miss $\sim 25\%$ of sources with mass-ratio $q > 6$ and up to $\sim 60 \%$ of highly precessing sources $(\chi_p > 0.5)$, after accounting for the approximate increase in background. In practice, the true observational bias can be even larger due to strict signal-consistency tests applied in searches. The observation of low spin, unequal-mass-ratio sources by Advanced LIGO design and Advanced Virgo may in part be due to these selection effects. The development of a search sensitive to high mass-ratio, precessing sources may allow the detection of new binaries whose spin properties would provide key insights into the formation and astrophysics of compact objects.

The main contribution to the effective shear modulus of neutron star crust can be calculated within Coulomb solid model and can be approximated by simple analytical expression for arbitrary (even multicomponent) composition. Here I consider correction associated with electron screening within Thomas-Fermi approximation. In particular, I demonstrate that for relativistic electrons (density $\rho>10^6$ g\,cm$^{-3}$) this correction can be estimated as $\delta \mu_\mathrm{eff}^\mathrm{V}= -9.4\times 10^{-4}\sum_Z n_Z Z^{7/3} e^2/a_\mathrm{e}$, where summation is taken over ion species, $n_Z$ is number density of ions with charge $Ze$, $k_\mathrm{TF}$ is Thomas-Fermi screening wave number. Finally, $a_\mathrm{e}=(4 \pi n_\mathrm{e}/3)^{-1/3}$ is electron sphere radius. Quasineutrality condition $n_\mathrm{e}=\sum_Z Z n_Z$ is assumed. This result holds true for arbitrary (even multicomponent and amorphous) matter and can be applied for neutron star crust and (dense) cores of white dwarfs. For example, the screening correction reduces shear modulus by $\sim 9$\% for $Z\sim40$, which is typical for inner layers of neutron star crust.

Electromagnetic Star-Planet Interaction (SPI) describes the phenomenon, when a planet couples to its host star via electromagnetic forces. Alfv\'en waves can establish such a coupling by forming Alfv\'en wings. SPI allows phenomena that we do not know from the Solar System. Wing-wing interaction is such an example, where the Alfv\'en wings of two planets merge and interact non-linearly. In this paper we focus on the effects that SPI has on other planets and the stellar wind. First, we analyse the different wave structures connected to SPI. The second part then investigates wing-wing interaction. Our study applies a magnetohydrodynamic model to describe a stellar system with multiple possible planets. As an example, we chose TRAPPIST-1 and its two innermost planets. We extended the PLUTO code to simulate collisions between atmospheric neutral particles and plasma ions. Neutral gas clouds imitate the planets and move through the simulation domain. That allows the simulation of fully time-dependent stellar systems. We analysed the wave structures, which result from the interaction between stellar wind and TRAPPIST-1 b. The inward going wave structure is an Alfv\'en wing. The outward going part of the interaction consists of an Alfv\'en wing, slow mode waves, the planetary wake and a slow shock. We quantified the strength of the respective wave perturbations at the outer planets to be on the order of 10\% to 40\% of the local background values of thermal, magnetic and dynamic pressure. Wing-wing interaction occurs due to the relative position of two planets during their conjunction and shows three phases. First there is an initial, non-linear intensification of the Poynting flux by $20\%$, an intermediate phase with reduced Poynting flux and a third phase when the Alfv\'en wing of planet c goes through planet b's wave structures with another intensification of the Poynting flux.

Ionizing radiation is known to have a destructive impact on biology by causing damage to the DNA, cells, and production of Reactive Oxygen Species (ROS) among other things. While direct exposure to high radiation dose is indeed not favorable for biological activity, ionizing radiation can, and in some cases is known to produce a number of biologically useful products. One such mechanism is the production of biologically useful products via charged particle-induced radiolysis. Energetic charged particles interact with surfaces of planetary objects such as Mars, Europa and Enceladus without much shielding from their rarefied atmospheres. Depending on the energy of said particles, they can penetrate several meters deep below the surface and initiate a number of chemical reactions along the way. Some of the byproducts are impossible to produce with lower-energy radiation (such as sunlight), opening up new avenues for life to utilize them. For each of these cases, we calculate the energy deposition rate as a function of depth, and estimate the energy availability for potential metabolic activity. We discuss various mechanisms through which life could support itself utilizing the byproducts of these ionizing radiation-induced reactions, such as chemoautotrophs using solvated electrons, extracellular electron transfer, and indirect electrophy to facilitate processes like carbon fixation, nitrogen fixation and sulfate reduction, and possibly for ATP production.

With data from Pantheon, we have at our disposal a sample of more than a thousand supernovae Ia covering a wide range of redshifts with good precision. Here we make fits to the corresponding Hubble--Lema\^itre diagram with various cosmological models, with intergalactic extinction, evolution of the luminosity of supernovae, and redshift components due to partially non-cosmological factors. The data are well fitted by the standard model to include dark energy, but there is a degeneracy of solutions with several other variables. Therefore, the Hubble--Lema\^itre diagram of SNe Ia cannot be used alone to infer the existence of the accelerated expansion scenario with dark energy. Within this degeneracy, models that give good fits to the data include the following alternative solutions: Einstein--de Sitter with gray extinction $a_V=1.2\times 10^{-4}$ Mpc$^{-1}$; linear Hubble--Lema\^itre law static Euclidean with gray extinction $a_V=0.4\times 10^{-4}$ Mpc$^{-1}$; Static Euclidean with tired light and gray extinction $a_V=2.8\times 10^{-4}$ Mpc$^{-1}$; Einstein--de Sitter with absolute magnitude evolution $\alpha =-0.10$ mag Gyr$^{-1}$; Friedmann model with $\Omega _M=0.07 - 0.29$, $\Omega _\Lambda =0$ and partially non-cosmological tired-light redshifts/blueshift with attenuation/enhancement $|K_i|<2.2\times 10^{-4}$ Mpc$^{-1}$ (although requiring calibration of $M$ incompatible with local SNe measurements).

Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-body simulation with 2.1 trillion particles in a volume of $25.7~\mathrm{Gpc}^3$, that has sufficient mass resolution to resolve dwarf-size systems, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$\nu^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$\nu^2$GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at $z\sim6-7$, the cosmic variance depends weakly on AGN luminosity, in particular for small survey areas (0.01 and 0.1 deg$^2$). This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST). In particular, we predict that the Euclid deep survey will find 120--240 (16--80) AGN -- depending on the model -- with rest-frame UV absolute magnitude brighter than $-20$ ($-20.5$) at $z=6.3$ ($z=7$) in the Euclid H-band deep survey (abridged).

In this study, we investigate the period changes of eight short-period Type II Cepheids of the BL Her subtype, i.e., with periods in the 1-4 day range. The $O-C$ diagrams for these stars are constructed using all suitable observational data from ground and space surveys. This spans a time interval of over one century and includes digitized photographic plates as well as photometry from the literature. The $O-C$ diagrams show parabolic evolutionary trends, which indicate the presence of both increasing and decreasing periods for these eight short period stars. These period changes are in good agreement with the recent theoretical evolutionary framework and stellar evolution models for BL Her stars. The pulsation stability test proposed by Lombard and Koen also suggests that the changes in the periods are real.

We attempt to establish links between a summary curve, or modulus summary curve, MSC, of the solar background magnetic field (SBMF) derived from Principal Component Analysis, with the averaged sunspot numbers (SSN). The comparison of MSC with the whole set of SSN reveals rather close correspondence of cycle timings, duration and maxima times for the cycles 12- 24, 6,7 and -4,-3. Although, in 1720-1760 and 1830-1860 there are discrepancies in maximum amplitudes of the cycles, durations and shifts of the maximum times between MSC and SSN curves. The MSC curve reveals pretty regular cycles with double maxima (cycles 1-4), triple maximum amplitude distributions for cycles 0 and 1 and for cycles -1 and -2 just before Maunder minimum. The MSC cycles in 1700-1750 reveal smaller maximal magnitudes in cycles -3 to 0 and in cycle 1-4 than the amplitudes of SSN, while cycles -2 to 0 have reversed maxima with minima with SSN. Close fitting of MSC or Bayesian models to the sunspot curve distorts the occurrences of either Maunder Minimum or/and modern grand solar minimum (2020-2053). These discrepancies can be caused by poor observations and by difference in solar magnetic fields responsible for these proxies. The dynamo simulations of toroidal and poloidal magnetic field in the grand solar cycle (GSC) from 1650 until 2050 demonstrate the clear differences between their amplitude variations during the GSC. The use of eigen vectors of SBMF can provide additional information to that derived from SSN that can be useful for understanding solar activity.

Parameterised Post-Newtonian Cosmology (PPNC) is a theory-agnostic framework for testing gravity in cosmology, which connects gravitational physics on small and large scales in the Universe. It is a direct extension of the Parameterised Post-Newtonian (PPN) approach to testing gravity in isolated astrophysical systems, and therefore allows constraints on gravity from vastly different physical regimes to be compared and combined. We investigate the application of this framework to a class of example scalar-tensor theories of gravity in order to verify theoretical predictions, and to investigate for the first time the scale-dependence of the gravitational couplings that appear within its perturbation equations. In doing so, we evaluate the performance of some simple interpolating functions in the transition region between small and large cosmological scales, as well as the uncertainties that using such functions would introduce into the calculation of observables. We find that all theoretical predictions of the PPNC framework are verified to high accuracy in the relevant regimes, and that simple interpolating functions perform well (but not perfectly) between these regimes. This study is an important step towards being able to use the PPNC framework to analyse cosmological datasets, and to thereby test if/how the gravitational interaction has changed as the Universe has evolved.

We present an analysis of the red giant component of the recurrent nova V3890 Sgr, using data obtained before and after its 2019 eruption. Its effective temperature is $T_{\rm eff}=3050\pm$200 K for $\log{g}=0.7$, although there are modest changes in $T_{\rm eff}$. There is an overabundance of both carbon ($0.20\pm0.05$~dex) and sodium ($1.0\pm0.3$~dex) relative to their solar values, possibly the result of ejecta from the 1990 nova eruption being entrained into the red giant photosphere. We find $^{12}$C/$^{13}$C $=25\pm2$, a value similar to that found in red giants in other recurrent novae. The interpretation of the quiescent spectrum in the 5--38$\,mu$m region requires the presence of photospheric SiO absorption and cool ($\sim400$~K) dust in the red giant environment. The spectrum in the region of the Na{\sc i} D lines is complex, and includes at least six interstellar components, together with likely evidence for interaction between ejecta from the 2019 eruption and material accumulated in the plane of the binary. Three recurrent novae with giant secondaries have been shown to have environments with different dust content, but photospheres with similar $^{12}$C/$^{13}$C ratios. The SiO fundamental bands most likely have a photospheric origin in the all three stars.

After the amazing discoveries by the GRAVITY collaboration in the last few years on the star S2 orbiting the black hole Sgr A* in the center of the Milky Way, we present a detailed investigation of the impact of gravitational lensing on the reconstruction of stellar orbits around this massive black hole. We evaluate the lensing astrometric effects on the stars S2, S38 and S55 and how these systematically affect the derived orbital parameters. The effect is below current uncertainties, but not negligible. With the addition of more observations on these stars, it will be possible to let the astrometric shift by lensing emerge from the statistical noise and be finally detected. By repeating the analysis on a smaller semimajor axis $a$ and various inclinations $i$, we are able to quantify the lensing effects on a broader range of parameters. As expected, for smaller semimajor axes and for nearly edge-on orbits lensing effects increase by about an order of magnitude.

Over the last few years, both ALMA and Spitzer/IRAC observations have revealed a population of likely massive galaxies at $z>3$ that was too faint to be detected in HST rest-frame ultraviolet imaging. However, due to the very limited photometry for individual galaxies, the true nature of these so-called HST-dark galaxies has remained elusive. Here, we present the first sample of such galaxies observed with very deep, high-resolution NIRCam imaging from the Early Release Science Program CEERS. 33 HST-dark sources are selected based on their red colours across 1.6 $\mu$m to 4.4 $\mu$m. Their physical properties are derived from 12-band multi-wavelength photometry, including ancillary HST imaging. We find that these galaxies are generally heavily dust-obscured ($A_{V}\sim2$ mag), massive ($\log (M/M_{\odot}) \sim10$), star-forming sources at $z\sim2-8$ with an observed surface density of $\sim0.8$ arcmin$^{-2}$. This suggests that an important fraction of massive galaxies may have been missing from our cosmic census at $z>3$ all the way into the reionization epoch. The HST-dark sources lie on the main-sequence of galaxies and add an obscured star formation rate density (SFRD) of $\mathrm{1.3^{+1.6}_{-1.0} \times 10^{-3} M_{\odot}/yr/Mpc^{3}}$ at $z\sim6$, similar to previous estimates. Our analysis shows the unique power of JWST to reveal this previously missing galaxy population and to provide a complete census of galaxies at $z=2-8$ based on rest-frame optical imaging.

We consider the consequences of a matter power spectrum which rises on small scales until eventually being cutoff by microphysical processes associated with the particle nature of dark matter. Evolving the perturbations of a weakly interacting massive particle from before decoupling until deep in the nonlinear regime, we show that nonlinear structure can form abundantly at very high redshifts. In such a scenario, dark matter annihilation is substantially increased after matter-radiation equality. Furthermore, since the power spectrum can be increased over a broad range of scales, the first star forming halos may form earlier than usual as well. The next challenge is determining how early Universe observations may constrain such enhanced dark matter perturbations.

There are indications that the third known eruption of the recurrent nova T CrB is imminent, and multi-wavelength observations prior to the eruption are important to characterise the system before it erupts. T CrB is known to display the SiO fundamental vibrational feature at 8$\,\mu$m. When the anticipated eruption occurs, it is possible that the shock produced when the ejected material runs into the wind of the red giant in the system may be traced using SiO maser emission. We have used the 100m Effelsberg Radio Telescope to search for $^{28}$SiO emission in the $\upsilon=1$, $\upsilon=2$, $J=1\rightarrow0$ transitions, at 43.122 GHz and 42.820~GHz respectively, while the system is in quiescence. We find no evidence for such emission.

A renewed interest about the origin of \emph{r}-process elements has been stimulated by the multi-messenger observation of the gravitational event GW170817, with the detection of both gravitational waves and electromagnetic waves corresponding to the merger of two neutron stars. Such phenomenon has been proposed as one of the main sources of the \emph{r}-process. However, the origin of the \emph{r}-process elements at different metallicities is still under debate. We aim at investigating the origin of the \emph{r}-process elements in the Galactic thin disc population. From the sixth internal data release of the \emph{Gaia}-ESO we have collected a large sample of Milky Way thin- and thick-disc stars for which abundances of Eu, O, and Mg are available. The sample consists of members of 62 open clusters, located at a Galactocentric radius from $\sim 5$ kpc to $\sim 20$ kpc in the disc, in the metallicity range $[-0.5, 0.4]$ and covering an age interval from 0.1 to 7 Gy, and about 1300 Milky Way disc field stars in the metallicity range $[-1.5, 0.5]$. We compare the observations with the results of a chemical evolution model, in which we varied the nucleosynthesis sources for the three considered elements. Our main result is that Eu in the thin disc is predominantly produced by sources with short lifetimes, such as magneto-rotationally driven SNe. There is no strong evidence for additional sources at delayed times. Our findings do not imply that there cannot be a contribution from mergers of neutron stars in other environments, as in the halo or in dwarf spheroidal galaxies, but such a contribution is not needed to explain Eu abundances at thin disc metallicities.

The Product of Exponentials (PoE) formulation is most commonly used in the field of robotics, but has recently been adapted for use in describing orbital motion. The PoE formula for orbital mechanics is an alternate method for defining and drawing an orbit based on its orbital elements set. Currently the PoE formula for orbital mechanics has only been derived through the first derivative (velocity). This work explores the second derivative of the adapted PoE formula for orbital mechanics, which gives a more complete description of the orbital motion of a satellite in a two-body system. This comprehensive approach employs a unified approach to account for all six time-varying orbital elements, therefore broadening the scope of the research and applications.

We present multiwavelength observations and a model for flat spectrum radio quasar NVSS J141922-083830, originally classified as a blazar candidate of unknown type (BCU II object) in the Third Fermi-LAT AGN Catalog (3LAC). Relatively bright flares (>3 magnitudes) were observed on 21 February 2015 (MJD 57074) and 8 September 2018 (MJD 58369) in the optical band with the MASTER Global Robotic Net (MASTER-Net) telescopes. Optical spectra obtained with the Southern African Large Telescope (SALT) on 1 March 2015 (MJD 57082), during outburst, and on 30 May 2017 (MJD 57903), during quiescence, showed emission lines at 5325\r{A} and at $\approx$3630\r{A} that we identified as the Mg II 2798\r{A} and C III] 1909\r{A} lines, respectively, and hence derived a redshift z = 0.903. Analysis of Fermi-LAT data was performed in the quiescent regime (5 years of data) and during four prominent flaring states in February-April 2014, October-November 2014, February-March 2015 and September 2018. We present spectral and timing analysis with Fermi-LAT. We report a hardening of the gamma-ray spectrum during the last three flaring periods, with a power-law spectral index $\Gamma = 2.0$-$2.1$. The maximum gamma-ray flux level was observed on 24 October 2014 (MJD 56954) at $(7.57 \pm 1.83) \times 10^{-7}$ ph~cm$^{-2}$s$^{-1}$. The multi-wavelength spectral energy distribution during the February-March 2015 flare supports the earlier evidence of this blazar to belong to the FSRQ class. The SED can be well represented with a single-zone leptonic model with parameters typical of FSRQs, but also a hadronic origin of the high-energy emission can not be ruled out.

Astrophysical black holes (BHs) can be fully described by their mass and spin. However, producing rapidly spinning ones is extremely difficult as the stars that produce them lose most of their angular momentum before the BH is formed. Binaries where the progenitor is paired with a low-mass star in a tight orbit can produce rapidly spinning BHs (through tides), whereas those with massive companions cannot (as they do not fit in such an orbit). A few rapidly-spinning black holes (BHs) have been observed paired with very massive companion stars, defying stellar-formation paradigm. Models which reduce the stellar-core--envelope interaction (and winds) do not match observations nor theory well; I show they also miss explaining the energetics. BH spins cannot be produced during stellar collapse; using orbital spin through explosion-fallback material does not match the observations; spinning BHs up through accepted mass-transfer channels takes longer than their lifetimes, it is usually discarded. I show that fast mass-transfer mechanisms, predicted to merge the BH and star, successfully spin the BHs up and show a mechanism to avoid said mergers and main dangers of the alternatives while naturally explaining the observations. The implications are potentially paradigm-shifting and far-reaching in the high-energy, BH astrophysics context.

The classical scenario of terrestrial planet formation is characterized by a phase of giant impacts among Moon-to-Mars mass planetary embryos. While the classic model and its adaptations have produced adequate analogs of the outer three terrestrial planets, Mercury's origin remains elusive. Mercury's high-core mass fraction compared to the Earth's is particularly outstanding. Among collisional hypotheses, this feature has been long interpreted as the outcome of an energetic giant impact among two massive protoplanets. Here, we revisit the classical scenario of terrestrial planet formation with focus on the outcome of giant impacts. We have performed a large number of N-body simulations considering different initial distributions of planetary embryos and planetesimals. Our simulations tested the effects of different giant planet configurations, from virtually circular to very eccentric configurations. We compare the giant impacts produced in our simulations with those that are more likely to account for the formation of Mercury and the Moon according to smoothed hydrodynamic simulations. Impact events that could lead to Moon's formation are observed in all our simulations with up to ~20% of all giant impacts, consistent with the range of the expected Moon-forming event conditions. On the other hand, Mercury-forming events via a single giant impact are extremely rare, accounting for less than ~1% of all giant impacts. Our results suggest that producing Mercury as a remnant of a single giant impact that strips out the mantle of a differentiated planetary object with Earth-like iron-silicate ratio is challenging and alternative scenarios may be required (e.g. multiple collisions).

We analyse the BOSS DR12 galaxy power spectrum data jointly with BAO data for three models of dark energy. We use recent measurements using a windowless estimator, and an independent and fast pipeline based on EFTofLSS implemented via the FAST-PT algorithm to compute the redshift-space loop corrections. We accelerate our analysis by using the BACCO linear emulator instead of a Boltzmann solver. We perform two sets of analyses: one with $3\sigma$ Planck priors on $A_s$ and $n_s$, and another that is CMB-free, without such priors. Firstly, we study $\Lambda$CDM, reproducing previous results obtained with the same estimator. We find a low value of $A_s$ in the CMB-free case, in agreement with many previous analyses of the BOSS data. We then study $w$CDM, finding a lower value of the amplitude in the CMB-free run, coupled with a preference for phantom dark energy with $w=-1.17^{+0.12}_{-0.11}$, again in broad agreement with previous results. Finally, we investigate the dark scattering model, which we label $wA$CDM. In the CMB-free analysis, we find a large degeneracy between the interaction strength $A$ and the amplitude $A_s$, hampering measurements of those parameters. On the contrary, in our run with a CMB prior, we are able to constrain the dark energy parameters to be $w=-0.972^{+0.036}_{-0.029}$ and $A = 3.9^{+3.2}_{-3.7}$, which show a 1$\sigma$ hint of interacting dark energy. This is the first measurement of this parameter and demonstrates the ability of this model to alleviate the $\sigma_8$ tension. Our analysis can be used as a guide for any model with scale-independent growth. Finally, we study the dependence of the results on the priors of the nuisance parameters and find these priors to be informative, with their broadening creating shifts in the contours. We argue for an in depth study of this issue, which can affect current and forthcoming analyses of LSS data.

The BICEP3 Polarimeter is a small aperture, refracting telescope, dedicated to the observation of the Cosmic Microwave Background (CMB) at 95GHz. It is designed to target degree angular scale polarization patterns, in particular the very-much-sought-after primordial B-mode signal, which is a unique signature of cosmic inflation. The polarized signal from the sky is reconstructed by differencing co-localized, orthogonally polarized superconducting Transition Edge Sensor (TES) bolometers. In this work, we present absolute measurements of the polarization response of the detectors for more than $\sim 800$ functioning detector pairs of the BICEP3 experiment, out of a total of $\sim 1000$. We use a specifically designed Rotating Polarized Source (RPS) to measure the polarization response at multiple source and telescope boresight rotation angles, to fully map the response over 360 degrees. We present here polarization properties extracted from on-site calibration data taken in January 2022. A similar calibration campaign was performed in 2018, but we found that our constraint was dominated by systematics on the level of $\sim0.5^\circ$. After a number of improvements to the calibration set-up, we are now able to report a significantly lower level of systematic contamination. In the future, such precise measurements will be used to constrain physics beyond the standard cosmological model, namely cosmic birefringence.

We present the serendipitous discovery of a late T-type brown dwarf candidate in JWST NIRCam observations of the Early Release Science Abell 2744 parallel field. The discovery was enabled by the sensitivity of JWST at 4~$\mu$m wavelengths and the panchromatic 0.9--4.5~$\mu$m coverage of the spectral energy distribution. The unresolved point source has magnitudes F115W = 27.95$\pm$0.15 and F444W = 25.84$\pm$0.01 (AB), and its F115W$-$F444W and F356W$-$F444W colors match those expected for other, known T dwarfs. We can exclude it as a reddened background star, high redshift quasar, or very high redshift galaxy. Comparison with stellar atmospheric models indicates a temperature of $T_{eff}$ $\approx$ 600~K and surface gravity $\log{g}$ $\approx$ 5, implying a mass of 0.03~M$_{\odot}$ and age of 5~Gyr. We estimate the distance of this candidate to be 570--720~pc in a direction perpendicular to the Galactic plane, making it a likely thick disk or halo brown dwarf. These observations underscore the power of JWST to probe the very low-mass end of the substellar mass function in the Galactic thick disk and halo.

Accurate models of the structural evolution of dark matter subhaloes, as they orbit within larger systems, are fundamental to understanding the detailed distribution of dark matter at the present day. Numerical simulations of subhalo evolution support the idea that the mass loss associated with tidal stripping is most naturally understood in energy space, with the particles that are the least bound being removed first. Starting from this premise, we recently proposed a zero-parameter "energy-truncation model" for subhalo evolution. We tested this model with simulations of tidal stripping of satellites with initial NFW profiles, and showed that the energy-truncation model accurately predicts both the mass loss and density profiles. In this work, we apply the model to a variety of Hernquist, Einasto and King profiles. We show that it matches the simulation results quite closely in all cases, indicating that it may serve as a universal model to describe tidally stripped collisionless systems. A key prediction of the energy-truncation model is that the central density of dark matter subhaloes is conserved as they lose mass; this has important implications for dark matter annihilation calculations, and for other observational tests of dark matter.

Recent JWST observations suggest an excess of $z\gtrsim10$ galaxy candidates above most theoretical models. Here, we explore how the interplay between halo formation timescales, star formation efficiency and dust attenuation affects the properties and number densities of galaxies we can detect in the early universe. We calculate the theoretical upper limit on the UV luminosity function, assuming star formation is 100% efficient and all gas in halos is converted into stars, and that galaxies are at the peak age for UV emission ($\sim10$ Myr). This upper limit is $\sim4$ orders of magnitude greater than current observations, implying these are fully consistent with star formation in $\Lambda$CDM cosmology. In a more realistic model, we use the distribution of halo formation timescales derived from extended Press-Schechter theory as a proxy for star formation rate (SFR). We predict that the galaxies observed so far at $z\gtrsim10$ are dominated by those with the fastest formation timescales, and thus most extreme SFRs and young ages. These galaxies can be upscattered by $\sim1.5$ mag compared to the median UV magnitude vs halo mass relation. This likely introduces a selection effect at high redshift whereby only the youngest ($\lesssim 10$ Myr), most highly star forming galaxies (specific SFR$\gtrsim$30 Gyr$^{-1}$) have been detected so far. Furthermore, our modelling suggests that redshift evolution at the bright end of the UV luminosity function is substantially affected by the build-up of dust attenuation. We predict that deeper JWST observations (reaching $m\sim30$) will reveal more typical galaxies with relatively older ages ($\sim100$ Myr) and less extreme specific SFRs ($\sim 10$ Gyr$^{-1}$ for a $M_\mathrm{UV} \sim -20$ galaxy at $z\sim10$).

The inspiral-merger-ringdown (IMR) consistency test checks the consistency of the final mass and final spin of a binary black hole merger remnant, independently inferred via the inspiral and merger-ringdown parts of the waveform. As binaries are expected to be nearly circularized when entering the frequency band of ground-based detectors, tests of general relativity (GR) currently employ quasi-circular waveforms. We quantify the effect of residual orbital eccentricity on the IMR consistency test. We find that eccentricity causes a significant systematic bias in the inferred final mass and spin of the remnant black hole at an orbital eccentricity (defined at $10$ Hz) of $e_0 \gtrsim 0.1$ in the LIGO band (for a total binary mass in the range $65-200 \,M_{\odot}$). For binary black holes observed by Cosmic Explorer (CE), the systematic bias becomes significant for $e_0 \gtrsim 0.015$ (for $200-600 \,M_{\odot}$ systems). This eccentricity-induced bias on the final mass and spin leads to an apparent inconsistency in the IMR consistency test, manifesting as a false violation of GR. Hence, eccentric corrections to waveform models are important for constructing a robust test of GR, especially for 3rd-generation (3G) detectors. We also estimate the eccentric corrections to the relationship between the inspiral parameters and the final mass and final spin; they are shown to be quite small.

We consider the general class of theories in which there is a new ultralight scalar field that mediates an equivalence principle violating, long-range force. In such a framework, the sun and the earth act as sources of the scalar field, leading to potentially observable location dependent effects on atomic and nuclear spectra. We determine the sensitivity of current and next-generation atomic and nuclear clocks to these effects and compare the results against the existing laboratory and astrophysical constraints on equivalence principle violating fifth forces. We show that in the future, the annual modulation in the frequencies of atomic and nuclear clocks in the laboratory caused by the eccentricity of the earth's orbit around the sun may offer the most sensitive probe of this general class of equivalence principle violating theories. Even greater sensitivity can be obtained by placing a precision clock in an eccentric orbit around the earth and searching for time variation in the frequency, as is done in anomalous redshift experiments. In particular, an anomalous redshift experiment based on current clock technology would already have a sensitivity to fifth forces that couple primarily to electrons at about the same level as the existing limits. Our study provides well-defined sensitivity targets to aim for when designing future versions of these experiments.

In the framework of gravothermal evolution of an ideal monatomic fluid, I examine the dynamical instability of the fluid sphere in ($N$+1) dimensions by exploiting Chandrasekhar's criterion to each quasistatic equilibrium along the sequence of the evolution. Once the instability is triggered, it would probably collapse into a black hole if no other interaction halts the process. From this viewpoint, the privilege of (3+1)-dimensional spacetime is manifest, as it is the marginal dimensionality in which the ideal monatomic fluid is stable but not too stable. Moreover, it is the unique dimensionality that allows stable hydrostatic equilibrium with positive cosmological constant. While all higher dimensional ($N>3$) spheres are genuinely unstable. In contrast, in (2+1)-dimensional spacetime it is too stable either in the context of Newton's theory of gravity or Einstein's general relativity. It is well known that the role of negative cosmological constant is crucial to have the Ba\~nados-Teitelboim-Zanelli (BTZ) black hole solution and the equilibrium configurations of a fluid disk. Owing to the negativeness of the cosmological constant, there is no unstable configuration for a homogeneous fluid disk to collapse into a naked singularity, which supports the cosmic censorship conjecture. However, BTZ holes of mass $\mathcal{M}_{\rm BTZ}>0$ could emerge from collapsing fluid disks. The implications of spacetime dimensionality are briefly discussed.

Orbital eccentricity is a key signature of dynamical binary black hole formation. The gravitational waves from a coalescing binary contain information about its orbital eccentricity, which may be measured if the binary retains sufficient eccentricity near merger. Dedicated waveforms are required to measure eccentricity. Several models have been put forward, and show good agreement with numerical relativity at the level of a few percent or better. However, there are multiple ways to define eccentricity for inspiralling systems, and different models internally use different definitions of eccentricity, making it difficult to directly compare eccentricity measurements. In this work, we systematically compare two eccentric waveform models, $\texttt{SEOBNRE}$ and $\texttt{TEOBResumS}$, by developing a framework to translate between different definitions of eccentricity. This mapping is constructed by minimizing the relative mismatch between the two models over eccentricity and reference frequency, before evolving the eccentricity of one model to the same reference frequency as the other model. We show that for a given value of eccentricity passed to $\texttt{SEOBNRE}$, one must input a $20$-$50\%$ smaller value of eccentricity to $\texttt{TEOBResumS}$ in order to obtain a waveform with the same empirical eccentricity. We verify this mapping by repeating our analysis for eccentric numerical relativity simulations, demonstrating that $\texttt{TEOBResumS}$ reports a correspondingly smaller value of eccentricity than $\texttt{SEOBNRE}$.

We establish new experiments to search for dark matter based on a model of a light scalar field with a dilaton-like coupling to the electromagnetic field, which is strongly motivated by superstring theory. We estimate the power of the photon signal in the process of a non-resonant scalar-photon transition and in a cavity resonator permeated by electric and magnetic fields. We show that existing cavity resonators employed in the experiments like ADMX have a low but non-vanishing sensitivity to the scalar-photon coupling. As a result, by re-purposing the results of the ADMX experiment, we find new limits on the scalar-photon coupling in the range of the scalar field masses from 2.7 to 4.2 $\mu$eV. We discuss possible modifications of this experiment, which enhance the sensitivity to the scalar field dark matter. We also propose a broadband experiment for scalar field dark matter searches based on a high-voltage capacitor. The estimated sensitivity of this experiment exceeds by nearly two orders in magnitude the sensitivity of the experiment based on molecular spectroscopy.

When a gravitational wave or a graviton travels through an electric or magnetic background, it could convert into a photon with some probability. In this paper, a dipole magnetic field is considered as this kind of background in both the Minkowski spacetime and the curved spacetime in the near-zone of a neutron star. In the former case, we find that the graviton traveling vertically rather than parallel to the background magnetic field could be more effectively converted into an electromagnetic radiation field. In the latter case, we focus on the situation, in which the graviton travels along the radial direction near a neutron star. The radius of a neutron star is about ten kilometers, so the gravitational wave with long wavelength or low frequency may bypass neutron stars by diffraction. For high frequency gravitational wave, the conversion probability is proportional to the distance square as that in the static electric or magnetic background case. The smaller the inclination angle between the dipole field and the neutron star north pole is, the larger magnetic amplitude will be. The term that described curved spacetime will slightly enhance this kind of probability. We estimate that this value is about the order of $\sim 10^{-14}- 10^{-10}$. Therefore, it is expectable that this kind of conversion process may have a potential to open a window for observing high frequency gravitational waves.

We study the dilution of dark matter (DM) relic density caused by the electroweak first-order phase transition (FOPT) in the singlet extension models, including the singlet extension of the standard model (xSM), of the two-Higgs-doublet model (2HDM+S) and the next-to-minimal supersymmetric standard model (NMSSM). We find that in these models the entropy released by the strong electroweak FOPT can dilute the DM density to 1/3 at most. Nevertheless, in the xSM and NMSSM where the singlet field configure is relevant to the phase transition temperature, the strong FOPT always happens before the DM freeze-out, making the dilution effect negligible for the current DM density. On the other hand, in the 2HDM+S where the DM freeze-out temperature is independent of FOPT, the dilution may salvage some parameter space excluded by excessive DM relic density or by DM direct detections.