Papers for Monday, Mar 13 2023

We are witnessing the dawn of gravitational wave (GW) astronomy. With currently available detectors, observations are restricted to GW frequencies in the range between ${\sim} 10\,\mathrm{Hz}$ and $10\,\mathrm{kHz}$, which covers the signals from mergers of compact objects. The launch of the space observatory LISA will open up a new frequency band for the detection of stellar interactions at lower frequencies. In this work, we predict the shape and strength of the GW signals associated with common-envelope interaction and merger events in binary stars, and we discuss their detectability. Previous studies estimated these characteristics based on semi-analytical models. In contrast, we used detailed three-dimensional magnetohydrodynamic simulations to compute the GW signals. We show that for the studied models, the dynamical phase of common-envelope events and mergers between main-sequence stars lies outside of the detectability band of the LISA mission. We find, however, that the final stages of common-envelope interactions leading to mergers of the stellar cores fall into the frequency band in which the sensitivity of LISA peaks, making them promising candidates for detection. These detections can constrain the enigmatic common-envelope dynamics. Furthermore, future decihertz observatories such as DECIGO or BBO would also be able to observe this final stage and the post-merger signal, through which we might be able to detect the formation of Thorne-\.Zytkow objects.

We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z=0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. Conservatively using a redshift-independent cut in stellar mass ($M_\star = 10^{11}\,{\rm M}_\odot$), and matching the stellar-mass distributions of our samples, we find 7 $\sigma$ evidence of $h_4$ increasing with cosmic time, from a median value of 0.019$\pm$0.002 at z=0.8 to 0.059$\pm$0.004 at z=0.06. Alternatively, we use a physically motivated sample selection, based on the mass distribution of the progenitors of local quiescent galaxies as inferred from numerical simulations; in this case, we find 10 $\sigma$ evidence. This evolution suggests that, over the last 7 Gyr, there has been a gradual decrease in the rotation-to-dispersion ratio and an increase in the radial anisotropy of the stellar velocity distribution, qualitatively consistent with accretion of gas-poor satellites. These findings demonstrate that massive galaxies continue to accrete mass and increase their dispersion support after becoming quiescent.

Free-floating planets (FFPs) are the lightest products of star formation and they carry important information on the initial conditions of the environment in which they were formed. They were first discovered in the 2000s but still few of them have been identified and confirmed due to observational challenges. This is a review of the last advances in the detection of these objects and the understanding of their origin. Several studies indicate that the observed fraction of FFPs outnumbers the prediction of turbulent fragmentation and suggest that many were formed in planetary systems that were later abandoned. The JWST will certainly constitute a new step further in the detection and characterisation of FFPs. To interpret these new observations, precise ages for the nearby star-forming regions in which they were formed will be necessary.

We present follow-up studies of three ultracompact hot subdwarf binaries. Using data from the Zwicky Transient Facility, we find orbital periods of 33.6, 37.3, and 36.9 minutes for ZTF 1946+3203, ZTF 0640+1738, and ZTF 0643+0318 respectively. The light curves show ellipsoidal variability of the hot subdwarf star with potential eclipses of an accretion disc. Phase-resolved spectroscopic observations with Keck were used to measure a radial velocity curve and atmospheric parameters of the hot subdwarf stars. ZTF J0643 shows evidence of accretion disc emission lines in the average spectrum. Combining light curve and spectroscopic fits will allow us to measure precise system properties such as masses, to determine the evolutionary history and future evolution of the system.

Some galaxies show little to no sign of active galactic nucleus (AGN) activity, yet exhibit strong coronal emission lines (CLs) relative to common narrow emission lines. Many of these coronal lines have ionization potentials of $\geq 100$ eV, thus requiring strong extreme UV and/or soft X-ray flux. It has long been thought that such events are powered by tidal disruption events (TDEs), but owing to a lack of detailed multi-wavelength follow-up, such a connection has not been firmly made. Here we compare coronal line emitters (CLEs) and TDEs in terms of their host-galaxy and transient properties. We find that the mid-infrared (MIR) colors of CLE hosts in quiescence are similar to TDE hosts. Additionally, many CLEs show evidence of a large dust reprocessing echo in their mid-infrared colors, a sign of significant dust in the nucleus. The stellar masses and star formation rates of the CLE hosts are consistent with TDE hosts, and both populations reside within the green valley. The blackbody properties of CLEs and TDEs are similar, with some CLEs showing hot (T $\geq 40,000$ K) blackbody temperatures. Finally, the location of CLEs on the peak-luminosity/decline-rate parameter space is much closer to TDEs than many other major classes of nuclear transients. Combined, these provide strong evidence to confirm the previous claims that CLEs are indeed TDEs in gas-rich environments. We additionally propose a stricter threshold of CL flux $\geq 1/3$ $\times$ [O III] flux to better exclude AGNs from the sample of CLEs.

Calculation of the mass transfer (MT) rate $\dot{M}_\text{d}$ of a Roche lobe overflowing star is a fundamental task in binary star evolution theory. Most of the existing MT prescriptions are based on a common set of assumptions that combine optically-thick and optically-thin regimes with different flow geometries. In this work, we develop a new model of MT based on the assumption that the Roche potential sets up a nozzle converging on the inner Lagrangian point and that the gas flows mostly along the axis connecting both stars. We derive a set of 1D hydrodynamic equations governing the gas flow with $\dot{M}_\text{d}$ determined as the eigenvalue of the system. The inner boundary condition directly relates our model to the structure of the donor obtained from 1D stellar evolution codes. We obtain algebraic solution for the polytropic equation of state (EOS), which gives $\dot{M}_\text{d}$ within a factor of 0.9 to 1.0 of existing optically-thick prescriptions and which reduces to the existing optically-thin prescription for isothermal gas. For a realistic EOS, we find that $\dot{M}_\text{d}$ differs by up to a factor of 4 from existing models. We illustrate the effects of our new MT model on $30\,M_\odot$ low-metallicity star undergoing intensive thermal time-scale MT and find that it is more likely to become unstable to L2 overflow and common-envelope evolution than for existing MT prescriptions. Our model provides a framework for including additional physics such as radiation or magnetic fields.

A variety of observational campaigns seek to test dark-matter models by measuring dark-matter subhaloes at low masses. Despite their predicted lack of stars, these subhaloes may be detectable through gravitational lensing or via their gravitational perturbations on stellar streams. To set measurable expectations for subhalo populations within LambdaCDM, we examine 11 Milky Way (MW)-mass haloes from the FIRE-2 baryonic simulations, quantifying the counts and orbital fluxes for subhaloes with properties relevant to stellar stream interactions: masses down to 10^6 Msun, distances < 50 kpc of the galactic center, across z = 0 - 1 (lookback time 0 - 8 Gyr). We provide fits to our results and their dependence on subhalo mass, distance, and lookback time, for use in (semi)analytic models. A typical MW-mass halo contains ~16 subhaloes >10^7 Msun (~1 subhalo >10^8 Msun) within 50 kpc at z = 0. We compare our results with dark-matter-only versions of the same simulations: because they lack a central galaxy potential, they overpredict subhalo counts by 2-10x, more so at smaller distances. Subhalo counts around a given MW-mass galaxy declined over time, being ~10x higher at z = 1 than at z = 0. Subhaloes have nearly isotropic orbital velocity distributions at z = 0. Across our simulations, we also identified 4 analogs of Large Magellanic Cloud satellite passages; these analogs enhance subhalo counts by 1.4-2.7 times, significantly increasing the expected subhalo population around the MW today. Our results imply an interaction rate of ~5 per Gyr for a stream like GD-1, sufficient to make subhalo-stream interactions a promising method of measuring dark subhaloes.

The discovery of giant quasar Ly$\alpha$ nebulae at $z>2$ has opened up the possibility to directly study in emission the Circumgalactic and Intergalactic Medium (CGM/IGM). However, the resonant nature of the Ly$\alpha$ line and its different emission mechanisms hamper the ability to constrain both the kinematics and physical properties of the CGM/IGM. Here, we present results of a pilot project aiming at the detection of CGM H$\alpha$ emission, a line which does not suffer from these limitations. To this end, we first used KCWI to detect Ly$\alpha$ emission around three bright quasars with $2.25<z<2.27$, a range which is free from bright IR sky lines for H$\alpha$, and then selected the most extended nebula for H$\alpha$ follow-up with MOSFIRE. Within the MOSFIRE slit, we detected H$\alpha$ emission extending up to 20 physical kpc with a total H$\alpha$ flux of F$_{ \textrm{H}\alpha}$=(9.5 $\pm$ 0.9) $\times$ 10$^{-18}$ erg s$^{-1}$ cm$^{-2}$. Considering the Ly$\alpha$ flux in the same region, we found F$_{ \textrm{Ly}\alpha}$/F$_{ \textrm{H}\alpha}$=3.7 $\pm$ 0.3 consistent with that obtained for the Slug Nebula at z$=2.275$ and with recombination radiation. This implies high densities or a very broad density distribution within the CGM of high-redshift quasars. Moreover, the H$\alpha$ line profile suggests the presence of multiple emitting components overlapping along our line-of-sight and relatively quiescent kinematics, which seems incompatible with either quasar outflows capable of escaping the potential well of the host halo or disk-like rotation in a massive halo ($>10^{12}$M$_{\odot}$).

We report on the first multi-wavelength Swift monitoring campaign performed on SDSS J164100.10+345452.7, a nearby narrow-line Seyfert 1 galaxy formerly known as radio quiet which was recently detected both in the radio (at 37 GHz) and in the $\gamma$-rays, which hints at the presence of a relativistic jet. During our 20-month Swift campaign, while pursuing the primary goal of assessing the baseline optical/UV and X-ray properties of J1641, we caught two radio flaring episodes, one each year. Our strictly simultaneous multi-wavelength data closely match the radio flare and allow us to unambiguously link the jetted radio emission of J1641. Indeed, for the X-ray spectra preceding and following the radio flare a simple absorbed power-law model is not an adequate description, and an extra absorption component is required. The average spectrum of J1641 can be best described by an absorbed power law model with a photon index $\Gamma=1.93\pm0.12$, modified by a partially covering neutral absorber with a covering fraction $f=0.91_{-0.03}^{+0.02}$. On the contrary, the X-ray spectrum closest to the radio flare does not require such extra absorber and is much harder ($\Gamma_{\rm flare} \sim 0.7\pm0.4$), thus implying the emergence of a further, harder spectral component. We interpret this as the jet emission emerging from a gap in the absorber. The fractional variability we derive in the optical/UV and X-ray bands are found to be lower than the typical values reported in the literature, since our observations of J1641 are dominated by the source being in a low state. Under the assumption that the origin of the 37 GHz radio flare is the emergence of a jet from an obscuring screen also observed in the X-rays, the derived total jet power is $P^{\rm tot}_{\rm jet}=3.5\times10^{42}$ erg s$^{-1}$, comparable to the lowest measured in the literature. [Abridged]

We forecast the ability of bispectrum estimators to constrain primordial non-Gaussianity using future photometric galaxy redshift surveys. A full-sky survey with photometric redshift resolution of $\sigma_z/(1+z)=0.05$ in the redshift range $0.2<z<2$ can provide constraints $\sigma(f^\mathrm{local}_\mathrm{NL})=3.4$, $\sigma(f^\mathrm{equil}_\mathrm{NL})=15$, and $\sigma(f^\mathrm{orth}_\mathrm{NL})=17$ for the local, equilateral, and orthogonal shapes respectively, delivering constraints on primordial non-Gaussianities competitive to those from the cosmic microwave background. We generalize these results by deriving a scaling relation for the constraints on the amplitude of primordial non-Gaussianity as a function of redshift error, depth, sky coverage, and nonlinear scale cutoff. Finally, we investigate the impact that photometric calibration errors on the largest scales will have on the constraining power of future experiments. We show that peculiar velocities reconstructed via kinetic Sunyaev Zeldovich tomography can be used to mitigate the impact of calibration errors on primordial non-Gaussianity constraints.

We present the first analysis of cosmic shear measured in DES Y3 that employs the entire range of angular scales in the data. To achieve this, we build upon recent advances in the theoretical modelling of weak lensing provided by a combination of $N$-body simulations, physical models of baryonic processes, and neural networks. Specifically, we use BACCOemu to model the linear and nonlinear matter power spectrum including baryonic physics, allowing us to robustly exploit scales smaller than those used by the DES Collaboration. We show that the additional data produce cosmological parameters that are tighter but consistent with those obtained from larger scales, while also constraining the distribution of baryons. In particular, we measure the mass scale at which haloes have lost half of their gas, $\log\,M_{\rm c}=14.38^{+0.60}_{-0.56}\log(h^{-1}{\rm M_{ \odot}})$, and a parameter that quantifies the weighted amplitudes of the present-day matter inhomogeneities, $S_8=0.799^{+0.023}_{-0.015}$. Our constraint on $S_8$ is statistically compatible with that inferred from the Planck satellite's data at the $0.9\sigma$ level. We find instead a $1.4\sigma$ shift in comparison to that from the official DES Y3 cosmic shear, because of different choices in the modelling of intrinsic alignment, non-linearities, baryons, and lensing shear ratios. We conclude that small scales in cosmic shear data contain valuable astrophysical and cosmological information and thus should be included in standard analyses.

We measured precise masses of the host and planet in OGLE-2003-BLG-235 system, when the lens and source were resolving, with 2018 Keck high resolution images. This measurement is in agreement with the observation taken in 2005 with the Hubble Space Telescope (HST). In 2005 data, the lens and sources were not resolved and the measurement was made using color-dependent centroid shift only. Nancy Grace Roman Space Telescope will measure masses using data typically taken within 3-4 years of the peak of the event which is much shorter baseline compared to most of the mass measurements to date. Hence, color dependent centroid shift will be one of the primary method of mass measurements for Roman. Yet, mass measurements of only two events (OGLE-2003-BLG-235 and OGLE-2005-BLG-071) are done using the color dependent centroid shift method so far. The accuracy of the measurements using this method are neither completely known nor well studied. The agreement of Keck and HST results, shown in this paper, is very important since this agreement confirms the accuracy of the mass measurements determined at a small lens-source separation using the color dependent centroid shift method. This also shows that with >100 high resolution images, Roman telescope will be able to use color dependent centroid shift at 3-4 years time baseline and produce mass measurements. We find that OGLE-2003-BLG-235 is a planetary system consists of a 2.34 +- 0.43M_Jup planet orbiting a 0.56 +- 0.06M_Sun K-dwarf host star at a distance of 5.26 +- 0.71 kpc from the Sun.

Voyager 1 and 2 crossed the heliopause at $\sim$122 AU in 2012 and $\sim$119 AU in 2018, respectively. It was quite a surprise because the thickness of the inner heliosheath obtained by the existing at that time models of the global heliosphere was significantly larger (by 20-40 AU). Until now, the problem of the heliosheath thickness has not been fully resolved. Earlier in the frame of an oversimplified toy model of nearly isothermal solar wind plasma it has been shown that the effect of electron thermal conduction may significantly reduce the thickness of the inner heliosheath. In this paper, we present the first results of our 3D kinetic-MHD model of the global heliosphere, where the effect of thermal electron conduction has been considered rigorously. The thermal conduction acts mainly along the magnetic field lines. Classical and saturated thermal fluxes are employed when appropriate. It is shown the effects of thermal conduction are significant. The thickness of the inner heliospheric is reduced. It is desired effect since it helps to reconcile the thickness obtained in the model with Voyager data. The other effects are the strong depletion of the heliosheath plasma temperature toward the heliopause and the increase of the plasma temperature in the supersonic solar wind upstream of the termination shock.

On 26 September 2022 the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, a satellite of the asteroid 65803 Didymos. Because it is a binary system, it is possible to determine how much the orbit of the satellite changed, as part of a test of what is necessary to deflect an asteroid that might threaten Earth with an impact. In nominal cases, pre-impact predictions of the orbital period reduction ranged from ~8.8 - 17.2 minutes. Here we report optical observations of Dimorphos before, during and after the impact, from a network of citizen science telescopes across the world. We find a maximum brightening of 2.29 $\pm$ 0.14 mag upon impact. Didymos fades back to its pre-impact brightness over the course of 23.7 $\pm$ 0.7 days. We estimate lower limits on the mass contained in the ejecta, which was 0.3 - 0.5% Dimorphos' mass depending on the dust size. We also observe a reddening of the ejecta upon impact.

We use observations with the infrared-optimized MagAO system and Clio camera in 3.9 $\mu$m light to place stringent mass constraints on possible undetected companions to Sirius A. We suppress the light from Sirius A by imaging it through a grating vector-apodizing phase plate coronagraph with 180-degree dark region (gvAPP-180). To remove residual starlight in post-processing, we apply a time-domain principal-components-analysis-based algorithm we call PCA-Temporal (PCAT), which uses eigen-time-series rather than eigen-images to subtract starlight. By casting the problem in terms of eigen-time-series, we reduce the computational cost of post-processing the data, enabling the use of the fully sampled dataset for improved contrast at small separations. We also discuss the impact of retaining fine temporal sampling of the data on final contrast limits. We achieve post-processed contrast limits of $1.5 \times 10^{-6}$ to $9.8 \times 10^{-6}$ outside of 0.75 arcsec which correspond to planet masses of 2.6 to 8.0 $M_J$. These are combined with values from the recent literature of high-contrast imaging observations of Sirius to synthesize an overall completeness fraction as a function of mass and separation. After synthesizing these recent studies and our results, the final completeness analysis rules out 99% of $\ge 9 \ M_J$ planets from 2.5-7 AU.

Neutrinos propagating in a dense neutrino gas, such as those expected in core-collapse supernovae (CCSNe) and neutron star mergers (NSMs), can experience fast flavor conversions on relatively short scales. This can happen if the neutrino electron lepton number ($\nu$ELN) angular distribution crosses zero in a certain direction. Despite this, most of the state-of-the-art CCSN and NSM simulations do not provide such detailed angular information and instead, supply only a few moments of the neutrino angular distributions. In this study we employ, for the \emph{first} time, a machine learning (ML) approach to this problem and show that it can be extremely successful in detecting $\nu$ELN crossings on the basis of its zeroth and first moments. We observe that an accuracy of $\sim95\%$ can be achieved by the ML algorithms, which almost corresponds to the Bayes error rate of our problem. Considering its remarkable efficiency and agility, the ML approach provides one with an unprecedented opportunity to evaluate the occurrence of FFCs in CCSN and NSM simulations \emph{on the fly}. We also provide our ML methodologies on \href{https://github.com/sajadabbar/ML-nu_FFI/tree/main}{GitHub}.

The long-term evolution of the solar system is chaotic. In some cases, chaotic diffusion caused by an overlap of secular resonances can increase the eccentricity of planets when they enter into a linear secular resonance, driving the system to instability. Previous work has shown that including general relativistic contributions to the planets' precession frequency is crucial when modelling the solar system. It reduces the probability that the solar system destabilizes within 5 Gyr by a factor of 60. We run 1280 additional N-body simulations of the solar system spanning 12.5 Gyr where we allow the GR precession rate to vary with time. We develop a simple, unified, Fokker-Planck advection-diffusion model that can reproduce the instability time of Mercury with, without, and with time-varying GR precession. We show that while ignoring GR precession does move Mercury's precession frequency closer to a resonance with Jupiter, this alone does not explain the increased instability rate. It is necessary that there is also a significant increase in the rate of diffusion. We find that the system responds smoothly to a change in the precession frequency: There is no critical GR precession frequency below which the solar system becomes significantly more unstable. Our results show that the long-term evolution of the solar system is well described with an advection-diffusion model.

Recent observations of high-redshift galactic disks ($z=1-3$) show a strong negative trend in the dark matter fraction $f_{DM}$ with increasing baryonic mass. For this to be true, the inner baryons must dominate over dark matter in early massive galaxies, as observed in the Milky Way today. If disks are established and dominant at early times, we show that stellar bars form promptly within these disks, leading to a high bar fraction at early times. New JWST observations provide the best evidence to date for mature stellar bars in this redshift range. The disk mass fraction $f_{disk}$ within $R_s=2.2 R_{disk}$ is the dominant factor in determining how rapidly a bar forms. Using 3D simulations of halo-disk-bulge galaxies, we confirm the "Fujii relation" for the exponential dependence of the bar formation time $\tau_{bar}$ as a function of $f_{disk}$. For $f_{disk} > 0.3$, the bar formation time declines exponentially fast with increasing $f_{disk}$. This relation is a challenge to simulators - barred models with inadequate resolution fall off this curve. We find the Fujii relation is offset if we adopt a physically motivated definition of bar formation. For the first time, we exploit the exponential growth timescale associated with a positive feedback cycle as the bar emerges from the underlying disk. For a halo mass $M_{halo}$, the same trend described by the Fujii relation is observed for the mass range relevant to systems at cosmic noon, but now the bar onset is slower for higher mass halos at a fixed $f_{disk}$. If baryons dominate over dark matter within $R = R_s$, we predict that a high fraction of bars will be found in high-redshift disks long before $z = 1$. Due to its widespread use in simulations, we investigate the Efstathiou-Lake-Negroponte criterion for bar instability: this sub-optimal parameter is inversely related to $f_{disk}$, with a secondary dependence on $M_{halo}$.

Many astronomical phenomena, including Fast Radio Bursts and Soft Gamma Repeaters, consist of brief distinct aperiodic events. The intervals between these events vary randomly, but there are periods of greater activity, with shorter mean intervals, and of lesser activity, with longer mean intervals. A single dimensionless parameter, the width of a log-normal function fitted to the distribution of waiting times between events, quantifies the variability of the activity. This parameter describes its dynamics in analogy to the critical exponents and universality classes of renormalization group theory. If the distribution of event strengths is a power law, the width of the log-normal fit is independent of the detection threshold and is a robust measure of the dynamics of the phenomenon.

We numerically study the origin of the multi-armed spiral structure observed in the circumnuclear gaseous mini-disks of nearby galaxies. We show that the presence of dust in such disks and its interaction with the gravitationally stable gaseous component leads to the development of a multi-armed spiral structure. As a particular example, we study the formation of the multi-armed spiral pattern in the mini-disk of the galaxy NGC 4736, for which the observational data for the rotation and the density distribution are available. We find that the multi-armed spiral structure grows in the stable gaseous mini-disk of NGC 4736 if the gas-to-dust ratio is about 5-20 percent. We also demonstrate that together with the dust concentration, the important factor for the development of instability is the size of the dust grains. A nonlinear multi-armed spiral pattern develops in the stable gaseous disk with sizes of grains larger than one micron. If future observations confirm the presence of a large amount of dust in the mini-disks of galaxies, this will pinpoint the mechanism of the formation of the multi-armed spiral structure in them.

Fast radio bursts (FRBs) were discovered only in 2007. However, the number of known events and sources of repeating bursts grows very rapidly. In the near future the number of events will be $\gtrsim 10^4$ and the number of repeaters $\gtrsim100$. Presently, there is a consensus that most of the sources of FRBs might be neutron stars (NSs) with large magnetic fields. These objects might have different origin as suggested by studies of their host galaxies which represent a very diverse sample: from regions of very active star formation to old globular clusters. Thus, in the following decade we expect to have a very large sample of events directly related to extragalactic magnetars of different origin. This might open new possibilities to probe various aspects of NS physics. In the review we briefly discuss the main directions of such future studies and summarize our present knowledge about FRBs and their sources.

We make use of recent estimates for the parameters of the Milky Way's halo globular clusters and study the influence of the galactic bar on the dynamics of these clusters by computing their orbits. We use both an axisymmetric and non-axisymmetric galactic potentials, which include the rotating elongated bar/bulge structure. We account for observational errors both in the positions and in the velocities of the globular clusters and explore the influence of the bar on cluster's evolution. This is contained in the angular momentum-total energy plane, (Lz,E), which is widely exploited as an indicator of the groups of globular clusters that originated from the same accretion event. Particular attention is devoted to the Gaia-Sausage/Enceladus and Pontus structures identified recently as two independent accretion events. Our study shows that it is not possible to identify GSE and Pontus as different merger events.

Recent high-$z$ quasar observations strongly indicate that super-Eddington accretion is a crucial phase to describe the existence of supermassive black holes (SMBHs) with $M_\mathrm{BH} \gtrsim 10^9 M_\odot$ at $z \gtrsim 7$. Motivated by the theoretical suggestion that the super-Eddington phase efficiently produces outflows and jets bright in radio bands, we search and find a super-Eddington radio-loud dust-obscured galaxy (DOG) J1406+0102 at $z=0.236$, through cross-matching of the infrared-bright DOGs of Noboriguchi et al. (2019) with the VLA/FIRST 1.4 GHz radio and the SDSS optical spectral catalog. DOG J1406+0102 shows broad components in the Balmer lines. Assuming those lines are from the broad line region, it gives BH mass estimation of $\log\ (M_\mathrm{BH}/M_\odot)=7.30 \pm 0.25$, and AGN luminosity of $\log (L_\mathrm{bol,[OIII]}/\mathrm{erg}~\mathrm{s}^{-1}) = 45.91\pm0.38$ estimated from the intrinsic [OIII] luminosity, resulting in super-Eddington accretion of $\lambda_\mathrm{Edd}\simeq 3$. We show that 1) DOG J1406+0102 is operating strong AGN feedback: the [OIII] outflow velocity exceeds the escape velocity of the host galaxy halo and the kinetic efficiency is obtained as $\approx$ 8% that can be sufficient to quench the host galaxy, 2) the expected future growth pathway of DOG J1406+0102 would join an over-massive BH trajectory and 3) radio-loud DOGs can provide a significant contribution to the high-energy ($\gtrsim$ 100 TeV) cosmic neutrino background if we assume DOG J1406+0102 as a representative of radio-loud DOGs.

We present new optical spectroscopic observations of U Geminorum obtained during a quiescent stage. We performed a radial velocity analysis of three Balmer emission lines yielding inconsistent results. Assuming that the radial velocity semi amplitude accurately reflects the motion of the white dwarf, we arrive at masses for the primary which are in the range of M_wd= 1.21 - 1.37 M_Sun. Based on the internal radial velocity inconsistencies and results produced from the Doppler tomography -- wherein we do not detect emission from the hot spot, but rather an intense asymmetric emission overlaying the disc, reminiscent of spiral arms -- we discuss the possibility that the overestimation of the masses may be due to variations of gas opacities and a partial truncation of the disc.

We report the first detection of hot molecular cores in the Small Magellanic Cloud, a nearby dwarf galaxy with 0.2 solar metallicity. We observed two high-mass young stellar objects in the SMC with ALMA, and detected emission lines of CO, HCO+, H13CO+, SiO, H2CO, CH3OH, SO, and SO2. Compact hot-core regions are traced by SO2, whose spatial extent is about 0.1 pc, and the gas temperature is higher than 100 K based on the rotation diagram analysis. In contrast, CH3OH, a classical hot-core tracer, is dominated by extended (0.2-0.3 pc) components in both sources, and the gas temperature is estimated to be 39+-8 K for one source. Protostellar outflows are also detected from both sources as high-velocity components of CO. The metallicity-scaled abundances of SO2 in hot cores are comparable among the SMC, LMC, and Galactic sources, suggesting that the chemical reactions leading to SO2 formation would be regulated by elemental abundances. On the other hand, CH3OH shows a large abundance variation within SMC and LMC hot cores. The diversity in the initial condition of star formation (e.g., degree of shielding, local radiation field strength) may lead to the large abundance variation of organic molecules in hot cores. This work, in conjunction with previous hot-core studies in the LMC and outer/inner Galaxy, suggests that the formation of a hot core would be a common phenomenon during high-mass star formation across the metallicity range of 0.2-1 solar metallicity. High-excitation SO2 lines will be a useful hot-core tracer in the low-metallicity environments of the SMC and LMC.

The spin-orbit obliquity of a planetary system constraints its formation history. A large obliquity may either indicate a primordial misalignment between the star and its gaseous disk or reflect the effect of different mechanisms tilting planetary systems after formation. Observations and statistical analysis suggest that system of planets with sizes between 1 and 4 R$_{\oplus}$ have a wide range of obliquities ($\sim0-30^{\circ}$), and that single- and multi-planet transiting have statistically indistinguishable obliquity distributions. Here, we revisit the ``breaking the chains'' formation model with focus in understanding the origin of spin-orbit obliquities. This model suggests that super-Earths and mini-Neptunes migrate close to their host stars via planet-disk gravitational interactions, forming chain of planets locked in mean-motion resonances. After gas-disk dispersal, about 90-99\% of these planetary systems experience dynamical instabilities, which spread the systems out. Using synthetic transit observations, we show that if planets are born in disks where the disk angular momentum is virtually aligned with the star's rotation spin, their final obliquity distributions peak at about $\sim$5 degrees or less, and the obliquity distributions of single and multi-planet transiting systems are statistically distinct. By treating the star-disk alignment as a free-parameter, we show that the obliquity distributions of single and multi-planet transiting systems only become statistically indistinguishable if planets are assumed to form in primordially misaligned natal disks with a ``tilt'' distribution peaking at $\gtrsim$10-20 deg. We discuss the origin of these misalignments in the context of star formation and potential implications of this scenario for formation models.

Protoplanetary disk around a just born young star contains a lot of cosmic dust. especially polycyclic-aromatic-hydrocarbon (PAH), which would become basic component to create biological organics. This study classified many astronomically observed infrared spectra of protoplanetary disks to three typical spectra. Type-A show well known astronomical bands of 6.2, 7.8, 8.6 and 11.3 micrometer. Whereas Type-B included unknown complex bands. Type-(A+B) was their mixed type. We tried to find specific molecule by Density Functional Theory (DFT) calculation. We found that Type-A could be explained by large PAH molecules of (C$_{23}$H$_{12}$) and (C$_{53}$H$_{18}$), which are hexagon-pentagon combined molecular structure. Background molecule of Type-B was smaller ones of (c-C$_3$H$_2$), (C$_5$H$_5$), (C$_9$H$_7$) and (C$_{12}$H$_8$). Type-(A+B) was reproduced well by mixing those molecules of A and B. Astronomical detailed observation shows that central star of Type-A has larger mass and higher temperature than that of Type-B. This suggests that at very early stage of our solar system, our protoplanetary disk had been made up by Type-B molecules. It was interesting that (C$_5$H$_5$) and (C$_9$H$_7$) of Type-B molecules has similar molecular structure with biological nucleic-acid on our earth. Type-B molecules was supposed to become the template for synthesizing biological organics and finally for creating our life.

Highly magnetized neutron stars are promising candidates to explain some of the most peculiar astronomical phenomena, for instance, fast radio bursts, gamma-ray bursts, and superluminous supernovae. Pulsations of these highly magnetized neutron stars are also speculated to produce detectable gravitational waves. In addition, pulsations are important probes of the structure and equation of state of the neutron stars. The major challenge in studying the pulsations of highly magnetized neutron stars is the demanding numerical cost of consistently solving the nonlinear Einstein and Maxwell equations under minimum assumptions. With the recent breakthroughs in numerical solvers, we investigate pulsation modes of non-rotating neutron stars which harbour strong purely toroidal magnetic fields of $10^{15-17}$ G through two-dimensional axisymmetric general-relativistic magnetohydrodynamics simulations. We show that stellar oscillations are insensitive to magnetization effects until the magnetic to binding energy ratio goes beyond 10%, where the pulsation mode frequencies are strongly suppressed. We further show that this is the direct consequence of the decrease in stellar compactness when the extreme magnetic fields introduce strong deformations of the neutron stars.

Cosmic microwave background radiation (CMB) observations are unavoidably contaminated by emission from various extra-galactic foregrounds, which must be removed to obtain reliable measurements of the cosmological signal. In this paper, we demonstrate CMB lensing reconstruction in AliCPT-1 after foreground removal, combine the two bands of AliCPT-1 (90 and 150~GHz) with Planck HFI bands (100, 143, 217 and 353~GHz) and with the WMAP-K band (23~GHz). In order to balance contamination by instrumental noise and foreground residual bias, we adopt the Needlet Internal Linear Combination (NILC) method to clean the E-map and the constrained Internal Linear Combination (cILC) method to clean the B-map. The latter utilizes additional constraints on average frequency scaling of the dust and synchrotron to remove foregrounds at the expense of somewhat noisier maps. Assuming 4 modules observing 1 season from simulation data, the resulting effective residual noise in E- and B-map are roughly $15~\mu{\rm K}\cdot{\rm arcmin}$ and $25~\mu{\rm K}\cdot{\rm arcmin}$, respectively. As a result, the CMB lensing reconstruction signal-to-noise ratio (SNR) from polarization data is about SNR$\,\approx\,$4.5. This lensing reconstruction capability is comparable to that of other stage-III small aperture millimeter CMB telescopes.

PeVatrons are the most powerful naturally occurring particle accelerators in the Universe. The identification of counterparts associated to astrophysical objects such as dying massive stars, molecular gas, star-forming regions, and star clusters is essential to clarify the underlying nature of the PeV emission, i.e., hadronic or leptonic. We present $^{12,13}$CO(J=2$\rightarrow$1) observations made with the 1.85~m radio-telescope of the Osaka Prefecture University toward the Cygnus OB7 molecular cloud, which contains the PeVatron candidate LHAASO J2108+5157. We investigate the nature of the sub-PeV (gamma-ray) emission by studying the nucleon density determined from the content of HI and H$_2$, derived from the CO observations. In addition to MML[2017]4607, detected via the observations of the optically thick $^{12}$CO(J=1$\rightarrow$0) emission, we infer the presence of an optically thin molecular cloud, named [FKT-MC]2022, whose angular size is 1.1$\pm$0.2$^{\circ}$. We propose this cloud as a new candidate to produce the sub-PeV emission observed in LHAASO J2108+5157. Considering a distance of 1.7 kpc, we estimate a nucleon (HI+H$_2$) density of 37$\pm$14 cm$^{-3}$, and a total nucleon mass(HI+H$_2$) of 1.5$\pm$0.6$\times$10$^4$ M$_{\odot}$. On the other hand, we confirm that Kronberger 82 is a molecular clump with an angular size of 0.1$^{\circ}$, a nucleon density $\sim$ 10$^3$ cm$^{-3}$, and a mass $\sim$ 10$^3$ M$_{\odot}$. Although Kronberger 82 hosts the physical conditions to produce the observed emission of LHAASO J2108+5157, [FKT-MC]2022 is located closer to it, suggesting that the latter could be the one associated to the sub-PeV emission. Under this scenario, our results favour a hadronic origin for the emission.

The source of the tension between local SN Ia based Hubble constant measurements and those from the CMB or BAO+BBN measurements is one of the most interesting unknowns of modern cosmology. Sample variance forms a key component of the error on the local measurements, and will dominate the error budget in the future as more supernovae are observed. Many methods have been proposed to estimate sample variance in many contexts, and we compared results from a number of them in Zhai \& Percival (2022), confirming that sample variance for the Pantheon supernovae sample does not solve the Hubble tension. We now extend this analysis to include a method based on analytically calculating correlations between the radial peculiar velocities of supernovae, comparing this technique with results from numerical simulations, which can be considered a non-linear Monte-Carlo solution that works similarly. We consider the dependence of these errors on the linear power spectrum and how non-linear velocities contribute to the error. Using this technique, and matching sample variance errors, we can define an effective volume for supernovae samples, finding that the Pantheon sample is equivalent to a top-hat sphere of radius $\sim220~h^{-1}$Mpc. We use this link between sample-variance errors to compute $\Delta H_{0}$ for idealised surveys with particular angular distributions of supernovae. For example, a half-sky survey at the Pantheon depth has the potential to suppress the sample variance of $H_{0}$ to $\sim0.1$ km s$^{-1}$Mpc$^{-1}$, a significant improvement compared with the current result. Finally, we consider the strength of large-scale velocity power spectrum required to explain the Hubble tension using sample variance, finding it requires an extreme model well beyond that allowed by other observations.

The accumulation of gluons inside nucleons, i.e., the gluon condensation, may lead to a characteristic broken power-law gamma-ray spectrum in high-energy nucleon collisions. Here we show that the observed spectra of at least 25 sources in the second Fermi Large Area Telescope Catalog of Gamma-ray Pulsars can be well fitted by such a broken power-law function that has only four free parameters. It strongly indicates that the gamma-ray emission from these pulsars is of hadronic origin, but with gluon condensation inside hadrons. It is well known that the quark-gluon distribution in a free nucleon is different from that in a bound nucleon. This work exposes the nuclear $A$-dependence of the gluon condensation effect, where $A$ refers to the baryon number. Our study reveals the gluon condensation under the condition of $A\to\infty$, which may open a new window for eavesdropping on the structure of compact stars on the sub-nuclear level.

Spiral arms are observed in numerous protoplanetary discs. These spiral arms can be excited by companions, either on bound or unbound orbits. We simulate a scenario where an unbound perturber, i.e. a flyby, excites spiral arms during a periastron passage. We run three-dimensional hydrodynamical simulations of a parabolic flyby encountering a gaseous protoplanetary disc. The perturber mass ranges from $10\, \rm M_J$ to $1\, \rm M_{\odot}$. The perturber excites a two-armed spiral structure, with a more prominent spiral feature for higher mass perturbers. The two arms evolve over time, eventually winding up, consistent with previous works. We focus on analysing the pattern speed and pitch angle of these spirals during the whole process. The initial pattern speed of the two arms are close to the angular velocity of the perturber at periastron, and then it decreases over time. The pitch angle also decreases over time as the spiral winds up. The spirals disappear after several local orbital times. An inclined prograde orbit flyby induces similar disc substructures as a coplanar flyby. A solar-mass flyby event causes increased eccentricity growth in the protoplanetary disc, leading to an eccentric disc structure which dampens over time. The spirals' morphology and the disc eccentricity can be used to search for potential unbound stars or planets around discs where a flyby is suspected. Future disc observations at high resolution and dedicated surveys will help to constrain the frequency of such stellar encounters in nearby star-forming regions.

Homogeneous metallicities and continuous high-precision light curves play key roles in studying the pulsation properties of RR Lyrae stars. By cross-matching with LAMOST DR6, we have determined 7 and 50 Non-Blazhko RRab stars in the Kepler and K2 fields, respectively, who have homogeneous metallicities determined from low-resolution spectra of the LAMOST-Kepler/K2 project. The Fourier Decomposition method is applied to the light curves of these stars provided by the Kepler space based telescope to determine the fundamental pulsation periods and the pulsation parameters. The calculated amplitude ratios of R21, R31 and the phase differences of {\phi}21, {\phi}31 are consistent with the parameters of the RRab stars in both the Globular Clusters and the Large Magellanic Cloud. We find a linear relationship between the phase differences {\phi}21 and {\phi}31, which is in good agreement with the results in previous literature. As far as the amplitude, we find that the amplitude of primary frequency A1 and the total amplitude Atot follow either a cubic or linear relationship. For the rise time RT, we do not find its relevance with the period of the fundamental pulsation mode P1, or Atot and {\phi}21. However, it might follow a linear relationship with R31. Based on the homogeneous metallicities, we have derived a new calibration formula for the relationship of period-{\phi}31-[Fe/H], which agrees well with the previous studies.

Images of 10 galaxies in F814W and F606W filters obtained on the Hubble Space Telescope (HST) are used to construct color-magnitude diagrams for the star population of these galaxies. The distances to the galaxies are estimated from the luminosity of the tip of the red giant branch. The galaxies examined here have radial velocities from 250 to 760 km/s relative to the centroid of the Local Group and distances ranging from 3.7 to 13.0 Mpc. Several other observed galaxies with low radial velocities are noted at distances beyond the limit of 13 Mpc.

Massive neutrinos and $f(R)$ modified gravity have degenerate observational signatures that can impact the interpretation of results in galaxy survey experiments, such as cosmological parameter estimations and gravity model tests. Because of this, it is important to investigate astrophysical observables that can break these degeneracies. Cosmic voids are sensitive to both massive neutrinos and modifications of gravity and provide a promising ground for disentangling the above mentioned degeneracies. In order to analyse cosmic voids in the context of non-$\Lambda$CDM cosmologies, we must first understand how well the current theoretical framework operates in these settings. We perform a suite of simulations with the RAMSES-based N-body code ANUBISIS, including massive neutrinos and $f(R)$ modified gravity both individually and simultaneously. The data from the simulations is compared to models of the void velocity profile and the void-halo cross-correlation function (CCF). This is done both with the real space simulation data as model input and by applying a reconstruction method to the redshift space data. In addition, we run Markov chain Monte Carlo (MCMC) fits on the data sets to assess the capability of the models to reproduce the fiducial simulation values of $f\sigma_8(z)$ and the Alcock-Paczy\`{n}ski parameter, $\epsilon$. The void modelling applied performs similarly for all simulated cosmologies, indicating that more accurate models and higher resolution simulations are needed in order to directly observe the effects of massive neutrinos and $f(R)$ modified gravity through studies of the void-galaxy CCF. The MCMC fits show that the choice of void definition plays an important role in the recovery of the correct cosmological parameters, but otherwise no clear distinction between the ability to reproduce $f\sigma_8$ and $\epsilon$ for the various simulations.

The determination of supermassive black hole (SMBH) masses is the key to understanding the host galaxy build-up and the SMBH mass assembly histories. The SMBH masses of non-local quasars are frequently estimated via the single-epoch virial black-hole mass estimators, which may suffer from significant biases. Here we demonstrate a new approach to infer the mass distribution of SMBHs in quasars by modelling quasar UV/optical variability. Our inferred black hole masses are systematically smaller than the virial ones by $0.3\sim 0.6$ dex; the $\sim 0.3$ dex offsets are roughly consistent with the expected biases of the virial black-hole mass estimators. In the upcoming time-domain astronomy era, our methodology can be used to constrain the cosmic evolution of quasar mass distributions.

We searched for potential atmospheric species in KELT-10b, focusing on sodium doublet lines (Na i; 589 nm) and the Balmer alpha line (H $\alpha$; 656 nm) in the transmission spectrum. Furthermore, we measured the planet-orbital alignment with the spin of its host star. We used the Rossiter-McLaughlin Revolutions technique to analyze the local stellar lines occulted by the planet during its transit. We used the standard transmission spectroscopy method to probe the planetary atmosphere, including the correction for telluric lines and the Rossiter-McLaughlin effect on the spectra. We analyzed two new light curves jointly with the public photometry observations. We do not detect signals in the Na i and H $\alpha$ lines within the uncertainty of our measurements. We derive the 3-sigma upper limit of excess absorption due to the planetary atmosphere corresponding to equivalent height Rp to 1.8Rp (Na i) and 1.9Rp (H $\alpha$). The analysis of the Rossiter-McLaughlin effect yields the sky-projected spin-orbit angle of the system $\lambda$ = -5.2 $\pm$ 3.4 and the stellar projected equatorial velocity $v_{eq} \sin{i_\star}$ = 2.58 $\pm$ 0.12 km/s. Photometry results are compatible within 1 -sigma with previous studies. We found no evidence of Na i and H $\alpha$, within the precision of our data, in the atmosphere of KELT-10b. Our detection limits allow us to rule out the presence of neutral sodium or excited hydrogen in an escaping extended atmosphere around KELT-10b. We cannot confirm the previous detection of Na i at lower altitudes with VLT/UVES. We note, however, that the Rossiter-McLaughlin effect impacts the transmission spectrum on a smaller scale than the previous detection with UVES. Analysis of the planet-occulted stellar lines shows the sky-projected alignment of the system, which is likely truly aligned due to tidal interactions of the planet with its cool (Teff < 6250 K) host star.

I consider a flow structure by which main sequence companions that enter a common envelope evolution (CEE) with giant stars might launch jets even when the accreted gas has a sub-Keplerian specific angular momentum. I first show that after a main sequence star enters the envelope of a giant star the specific angular momentum of the accreted gas is sub-Keplerian but still sufficiently large for the accreted gas to avoid two conical-like openings along the two opposite polar directions. I suggest that the high-pressure zone that the accreted gas builds around the main sequence equatorial plane accelerates outflows along these polar opening. Most of the inflowing gas is deflected to the polar outflows, now turning to two opposite jets. The actual mass that the main sequence star accretes is only a small fraction, ~0.1, of the inflowing gas. However, the gravitational energy that this gas releases powers the inflow-outflow streaming of gas and adds energy to the common envelope ejection. This flow structure might take place during a grazing envelope evolution if occurs, during the early CEE, and possibly in some post-CEE cases. This study increases the parameter space for main sequence stars to launch jets in shaping some planetary nebulae, in adding energy to mass removal in CEE, and in powering some intermediate luminosity optical transients.

The NASA Lucy mission is aimed at the study of the very interesting population of Jupiter Trojans, considered as time capsules from the origin of our solar system. During its journey, the mission will pass near a main belt asteroid, Donaldjohanson. Recently, NASA has announced that a new asteroid in the belt will also be visited by Lucy: 152830 Dinkinesh (1999 VD57). The main goal of this work is to characterise this newly selected target, asteroid Dinkinesh, in order to provide critical information to the mission team. To achieve it, we have obtained visible spectra, colour photometry, and time-series photometry of Dinkinesh, using several telescopes located at different observatories. For the spectra we used the 10.4m Gran Telescopio Canarias (GTC), in the island of La Palma (Spain); for the colour photometry the 4.3m Lowell Discovery Telescope (LDT), near Happy Jack, Arizona (USA) was used; and for the time-series photometry we used the 82cm IAC80 telescope located in the island of Tenerife (Spain). Both visible spectrum and reflectance values computed from colour photometry show that Dinkinesh is an S-type asteroid, i.e., it is composed mainly of silicates and some metal. According to observations done by the NEOWISE survey, S-type asteroids have typical geometric albedo of $p_V$ = 0.223 $\pm$ 0.073. From our time-series photometry, we obtained an asteroid mean magnitude $r'$ = 19.99 $\pm$ 0.05, which provides an absolute magnitude $H_{r'}$ = 17.53 $\pm$ 0.07 assuming $G$ = 0.19 $\pm$ 0.25 for S-types. Using our colour-photometry, we transformed $H_{r'}$ to $H_V$ = 17.48 $\pm$ 0.05. This value of absolute magnitude combined with the geometric albedo provides a mean diameter for Dinkinesh of $\sim$900 m, ranging between a minimum size of 542 m and a maximum size of 1309 m.

Context. Recent observations of GRB 200415A, a short and very bright pulse of $\gamma$-rays, have been claimed to be an extragalactic magnetar giant flare (MGF) whose proposed host galaxy is the nearby ${\mathrm{NGC} \, 253}$. However, as the redshift of the transient object was not measured, it is possible that the measured location of the transient on the celestial sphere and the location of the local galaxy merely coincided. Thus, its real progenitor could have been arbitrarily far away, leading possibly to a much larger luminosity of the transient, and leaving the standard model of short gamma-ray bursts (sGRBs), the merger of two compact objects, as an explanation for the observations. Aims. In this study, our aim is to compute the false-alarm rate for the misinterpretation of sGRBs as magnetars in a given observation period. Methods. We simulate synthetic surveys of sGRB observations in a time period of 14 years corresponding to the operation period of the Gamma-ray Burst Monitor (GBM) detector. For all sGRBs that align on the sky with a nearby Local Volume galaxy, we generate realistic data which is folded through the response of the GBM. To identify candidates of sGRBs that may be misinterpreted as magnetars, six selections (spatial, star formation rate, GBM trigger, duration, isotropic energy release, and fluence) are applied to the simulated surveys. Results. In a non-negligible fraction, 15.7 %, of the simulated surveys, we identify at least one sGRB that has the same characteristics as a magnetar giant flare and could be thus misinterpreted as magnetar. Thus, we conclude that the selections that were proposed in previous work to unambiguously identify an extragalactic magnetar giant flare are not sufficient.

Tight Triple Systems have stars in a hierarchical configuration with a third star orbiting the inner binary with a period of fewer than 1000 days. Such systems are important for understanding the formation and evolution of stars in multiple systems. Having a detached eclipsing binary (DEB) as one of its components allows us to obtain precise stellar and orbital parameters of these systems. We discuss the process to obtain accurate parameters of these systems using high-resolution spectroscopy, radial velocity measurements, and precise space-based photometry. This enables us to have a 3D geometrical picture as well as the metallicity, age, and evolutionary status of these systems.

We present a detailed exposition on the prospects of formation of Primordial Black Holes (PBHs) during Slow Roll (SR) to Ultra Slow Roll (USR) transitions in the framework of single-field inflation. We use effective field theory (EFT) approach in order to keep the analysis model-independent and applicable to both the canonical and non-canonical cases. We show in detail how renormalizing the power spectrum to one loop order in $P(X,\phi)$ theories severely limits the prospects for PBH formation in a single-field inflationary framework. We demonstrate that for the allowed range of effective sound speed, $1<c_s<1.17$, the consistency of one-loop corrected power spectrum leaves a small window for black hole masses, $M_{\rm PBH}\sim \mathcal{O}(10^2-10^3)$gm to have sufficient e-foldings, $\Delta {\cal N}_{\rm Total}\sim {\cal O}(54-59)$ for inflation. We confirm that adding a SR regime after USR before the end of inflation, does not significantly alter our conclusions. Our findings strictly rule out the possibility of generating large masses of PBHs from all possible models of single field inflation (canonical and non-canonical) and mature into a "no-go theorem" for the class of mentioned theories.

In their most recent observing run, the LIGO-Virgo-KAGRA (LVK) Collaboration observed gravitational waves (GWs) from compact binary mergers with highly asymmetric mass ratios for the first time, including both binary black hole (BBH) and neutron-star--black hole (NSBH) systems. It appears that NSBHs with mass ratios $q \simeq 0.2$ are more common than equally asymmetric BBHs, but the reason for this remains unclear. In this work, we use the binary population synthesis code COSMIC to investigate the evolutionary pathways leading to the formation of asymmetric compact binaries and the factors that cause them to merge. We find that within the context of isolated binary stellar evolution, most asymmetric binary mergers start off as asymmetric stellar binaries. Because of the initial asymmetry, these systems tend to first undergo a dynamically unstable mass transfer phase. However, after the first star collapses to a compact object, the mass ratio is close to unity and the second phase of mass transfer is usually stable. According to our simulations, this stable mass transfer fails to shrink the orbit enough on its own for the system to merge. Instead, we find that the natal kick received by the second-born compact object during its stellar collapse plays a key role in determining how many of these systems can merge. For the most asymmetric systems with mass ratios $q \leq 0.1$, the merging systems in our models receive an average kick magnitude of 255 km s$^{-1}$ during the second collapse, while the average kick for non-merging systems is 59 km s$^{-1}$. Because lower mass compact objects, like NSs, are expected to receive larger natal kicks than higher mass BHs, this may explain why asymmetric NSBH systems merge more frequently than asymmetric BBH systems.

The time-domain (TD) surveys of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) yield high-cadence radial velocities, paving a new avenue to study binary systems including compact objects. In this work, we explore LAMOST TD spectroscopic data of four K2 plates and present a sample of six single-lined spectroscopic binaries that may contain compact objects. We conduct analyses using phase-resolved radial velocity measurements of the visible star, to characterize each source and to infer the properties of invisible companion. By fitting the radial velocity curves for the six targets, we obtain accurate orbital periods, ranging from $\sim$ (0.6-6) days, and radial velocity semi-amplitudes, ranging from $\sim$ (50-130) km s$^{-1}$. We calculate the mass function of the unseen companions to be between 0.08 and 0.17 $M_{\odot}$. Based on the mass function and the estimated stellar parameters of the visible star, we determine the minimum mass of the hidden star. Three targets, J034813, J063350, and J064850, show ellipsoidal variability in the light curves from K2, ZTF, and TESS surveys. Therefore, we can put constraints on the mass of the invisible star using the ellipsoidal variability. We identify no X-ray counterparts for these targets except for J085120, of which the X-ray emission can be ascribed to stellar activity. We note that the nature of these six candidates is worth further characterization utilizing multi-wavelength follow-up observations.

The flux of energetic particles originating from the Sun fluctuates during the solar cycles. It depends on the number and properties of Active Regions (ARs) present in a single day and associated solar activities, such as solar flares and coronal mass ejections (CMEs). Observational records of the Space Weather Prediction Center (SWPC NOAA) enable the creation of time-indexed databases containing information about ARs and particle flux enhancements, most widely known as Solar Energetic Particle events (SEPs). In this work, we utilize the data available for Solar Cycles 21-24, and the initial phase of Cycle 25 to perform a statistical analysis of the correlation between SEPs and properties of ARs inferred from the McIntosh and Hale classifications. We find that the complexity of the magnetic field, longitudinal location, area, and penumbra type of the largest sunspot of ARs are most correlated with the production of SEPs. It is found that most SEPs ($\approx$60\%, or 108 out of 181 considered events) were generated from an AR classified with the 'k' McIntosh subclass as the second component, and some of these ARs are more likely to produce SEPs if they fall in a Hale class with $\delta$ component. It is confirmed that ARs located in the western hemisphere produced the most SEPs recorded on the Earth's orbit. The resulting database containing information about SEP events and ARs is publicly available and can be used for the development of Machine Learning (ML) models to predict the occurrence of SEPs.

The Visible Spectro-Polarimeter (ViSP) of the NSF Daniel K. Inouye Solar Telescope (DKIST) collected its Science Verification data on May 7-8, 2021. The instrument observed multiple layers of a sunspot atmosphere simultaneously, in passbands of Ca-II 397 nm (H-line), Fe-I 630 nm, and Ca-II 854 nm, scanning the region with a spatial sampling of 0.041" and average temporal cadence of 7.76 seconds, for a 38.8 minute duration. The slit moves southward across the plane-of-the-sky at 3.83 km/s. The spectropolarimetric scans exhibit prominent oscillatory 'ridge' structures which lie nearly perpendicular to the direction of slit motion (north to south). These ridges are visible in maps of line intensity, central wavelength, line width, and both linear and circular polarizations. Contemporaneous Atmospheric Imaging Assembly observations indicate these ridges are purely temporal in character and likely attributed to the familiar chromospheric 3-minute umbral oscillations. We observe in detail a steady umbral flash near the center of the sunspot umbra. Although bad seeing limited the spatial resolution, the unique high signal-to-noise enable us to estimate the shock Mach numbers (= 2), propagation speeds (= 9 km/s), and their impact on longitudinal magnetic field (delta B = 50 G), gas pressure, and temperature (delta T/T = 0.1) of the subshocks over 30 seconds. We also find evidence for rarefaction waves situated between neighboring wave-train shocks. The Ca-II 854 nm line width is steady throughout the umbral flash except for a sharp 1.5 km/s dip immediately before, and comparable spike immediately after, the passage of the shock front. This zig-zag in line width is centered on the subshock and extends over 0.4".

We present a spectro-photometric study of a mass-complete sample of quiescent galaxies at $1.0 < z < 1.3$ with $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}}) \geq 10.3$ drawn from the VANDELS survey, exploring the relationship between stellar mass, age and star-formation history. Within our sample of 114 galaxies, we derive a stellar-mass vs stellar-age relation with a slope of $1.20^{+0.28}_{-0.27}$ Gyr per decade in stellar mass. When combined with recent literature results, we find evidence that the slope of this relation remains consistent over the redshift interval $0<z<4$. The galaxies within the VANDELS quiescent display a wide range of star-formation histories, with a mean star-formation timescale of $1.5\pm{0.1}$ Gyr and a mean quenching timescale of $1.4\pm{0.1}$ Gyr. We also find a large scatter in the quenching timescales of the VANDELS quiescent galaxies, in agreement with previous evidence that galaxies at $z \sim 1$ cease star formation via multiple mechanisms. We then focus on the oldest galaxies in our sample, finding that the number density of galaxies that quenched before $z = 3$ with stellar masses $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}}) \geq 10.6$ is $ 1.12_{-0.72}^{+1.47} \times 10^{-5} \ \mathrm{Mpc}^{-3}$. Although uncertain, this estimate is in good agreement with the latest observational results at $3<z<4$, tentatively suggesting that neither rejuvenation nor merger events are playing a major role in the evolution of the oldest massive quiescent galaxies within the redshift interval $1<z<3$.

Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light curve observations. We present the oscillation properties of high resolution hydrodynamic simulations of a 25 $\mathrm{M}_\odot$ star performed with the PPMStar code. The model star has a convective core mass of $\approx\, 12\, \mathrm{M}_\odot$ and approximately half of the envelope simulated. From this simulation, we extract light curves from several directions, average them over each hemisphere, and process them as if they were real photometric observations. We show how core convection excites waves with a similar frequency as the convective time scale in addition to significant power across a forest of low and high angular degree $l$ modes. We find that the coherence of these modes is relatively low as a result of their stochastic excitation by core convection, with lifetimes on the order of 10s of days. Thanks to the still significant power at higher $l$ and this relatively low coherence, we find that integrating over a hemisphere produces a power spectrum that still contains measurable power up to the Brunt--V\"ais\"al\"a frequency. These power spectra extracted from the stable envelope are qualitatively similar to observations, with same order of magnitude yet lower characteristic frequency. This work further shows the potential of long-duration, high-resolution hydrodynamic simulations for connecting asteroseismic observations to the structure and dynamics of core convection and the convective boundary.

Mass and angular momentum are key parameters of galaxies. Their co-evolution establishes an empirical relation between the specific stellar angular momentum j* and the stellar mass M* that depends on morphology. In this work, we measure j* in a sample of 32 early type galaxies (ETGs) from the ePN.S survey, using full 2D kinematic information out to a mean 6Re. We present lambda and j* profiles. We derive the distribution of these galaxies on the total j*-M* plane and determine the ratio between the stellar j* and the specific angular momentum of the host dark matter halo. The radially extended, 2D kinematic data show that the stellar halos of ETGs do not contain large stellar mass fractions of high j*. The j*-profiles of fast-rotator ETGs are largely converged within the range of the data. For slow rotators, j* is still rising and is estimated to increase beyond 6Re by up to 40%. More than 60% of their halo angular momentum is in misaligned rotation. We find that the ePN.S ETG sample displays the well-known correlation between j*, M*, and morphology: elliptical galaxies have systematically lower j* than similar mass S0 galaxies. However, fast and slow rotators lie on the same relation within errors with the slow rotators falling at the high M* end. A power-law fit to the j*-M* relation gives a slope of 0.55+-0.17 for the S0s and 0.76+-0.23 for the ellipticals, with normalisation about 4 and 9 times lower than spirals, respectively. The estimated retained fraction of angular momentum at 10^10-10^10.5 Msun is 25% for S0s and >10% for ellipticals, and decreases by 1.5 orders of magnitude at M*~10^12 Msun. Our results show that ETGs have substantially lower j* than spiral galaxies with similar M*. Their j* must be lost during their evolution, and/or retained in the hot gas component and the satellite galaxies that have not yet merged with the central galaxy. [abridged]

We present the cosmological implications of measurements of void-galaxy and galaxy-galaxy clustering from the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample (MGS), Baryon Oscillation Spectroscopic Survey (BOSS), and extended BOSS (eBOSS) luminous red galaxy catalogues from SDSS Data Release 7, 12, and 16, covering the redshift range $0.07 < z < 1.0$. We fit a standard $\Lambda$CDM cosmological model as well as various extensions including a constant dark energy equation of state not equal to $-1$, a time-varying dark energy equation of state, and these same models allowing for spatial curvature. Results on key parameters of these models are reported for void-galaxy and galaxy-galaxy clustering alone, both of these combined, and all these combined with measurements from the cosmic microwave background (CMB) and supernovae (SN). For the combination of void-galaxy and galaxy-galaxy clustering plus CMB and SN, we find tight constraints of $\Omega_\mathrm{m} = 0.3127\pm 0.0055$ for a base $\Lambda$CDM cosmology, $\Omega_\mathrm{m} = 0.3172\pm 0.0061, w = -0.930\pm 0.039$ additionally allowing the dark energy equation of state $w$ to vary, and $\Omega_\mathrm{m} = 0.3239\pm 0.0085, w = -0.889\pm 0.052, \mathrm{and}\ \Omega_\mathrm{k} = -0.0031\pm 0.0028 $ further extending to non-flat models.

Recent atmospheric models for brown dwarfs suggest that the existence of clouds in substellar objects is not needed to reproduce their spectra, nor their rotationally-induced photometric variability, believed to be due to the heterogeneous cloud coverage of brown dwarf atmospheres. Cloud-free atmospheric models also predict that their flux should not be polarized, as polarization is produced by the light-scattering of particles in the inhomogeneous cloud layers of brown dwarf atmospheres. To shed light on this dichotomy, we monitored the linear polarization and photometric variability of the most variable brown dwarf, 2MASS J21392216+0220185. We used FORS2 at the UT1 telescope to monitor the object in the z-band for six hours, split on two consecutive nights, covering one-third of its rotation period. We obtained the Stokes parameters, and we derived its time-resolved linear polarization, for which we did not find significant linear polarization (P = 0.14+\-0.07 %). We modeled the linear polarimetric signal expected assuming a map with one or two spot-like features and two bands using a polarization-enabled radiative-transfer code. We obtained values compatible with the time-resolved polarimetry obtained for 2MASS J21392216+0220185. The lack of significant polarization might be due to photometric variability produced mostly by banded structures or small-scale vortices, which cancel out the polarimetric signal from different regions of the dwarf's disk. Alternatively, the lack of clouds in 2MASS J21392216+0220185 would also explain the lack of polarization. Further linear polarimetric monitoring of 2MASS J21392216+0220185, during at least one full rotational period, would help to confirm or discard the existence of clouds in its atmosphere.

Solar wind turbulence is often perceived as weakly compressible and the density fluctuations remain poorly understood both theoretically and observationally. Compressible magnetohydrodynamic simulations provide useful insights into the nature of density fluctuations. We discuss a few important effects related to 3D simulations of turbulence and in-situ observations. The observed quantities such as the power spectrum and variance depend on the angle between the sampling trajectory and the mean magnetic field due to anisotropy of the turbulence. The anisotropy effect is stronger at smaller scales and lower plasma beta. Additionally, in-situ measurements tend to exhibit a broad range of variations, even though they could be drawn from the same population with the defined averages, so a careful averaging may be needed to reveal the scaling relations between density variations and other turbulence quantities such as turbulent Mach number from observations.

The Markov chain Monte Carlo method (MCMC), especially the Metropolis-Hastings (MH) algorithm, is a widely used technique for sampling from a target probability distribution $P$ on a state space $\Omega$ and applied to various problems such as estimation of parameters in statistical models in the Bayesian approach. Quantum algorithms for MCMC have been proposed, yielding the quadratic speedup with respect to the spectral gap $\Delta$ compered to classical counterparts. In this paper, we consider the quantum version of the MH algorithm in the case that calculating $P$ is costly because the log-likelihood $L$ for a state $x\in\Omega$ is obtained via computing the sum of many terms $\frac{1}{M}\sum_{i=0}^{M-1} \ell(i,x)$. We propose calculating $L$ by quantum Monte Carlo integration and combine it with the existing method called quantum simulated annealing (QSA) to generate the quantum state that encodes $P$ in amplitudes. We consider not only state generation but also finding a credible interval for a parameter, a common task in Bayesian inference. In the proposed method for credible interval calculation, the number of queries to the quantum circuit to compute $\ell$ scales on $\Delta$, the required accuracy $\epsilon$ and the standard deviation $\sigma$ of $\ell$ as $\tilde{O}(\sigma/\epsilon^2\Delta^{3/2})$, in contrast to $\tilde{O}(M/\epsilon\Delta^{1/2})$ for QSA with $L$ calculated exactly. Therefore, the proposed method is advantageous if $\sigma$ scales on $M$ sublinearly. As one such example, we consider parameter estimation in a gravitational wave experiment, where $\sigma=O(M^{1/2})$.

The Kramers-Kronig relation is a well-known relation, especially in the field of optics. The key to this relation is the causality that output comes only after input. We first show that gravitational lensing obeys the causality in the sense that (electromagnetic/gravitational) waves emitted from the source arrive at an observer only after the arrival of the signal in geometrical optics. This is done by extending the previous work which is based on the thin lens approximation. We then derive the Kramers-Kronig relation in gravitational lensing, as the relation between real and imaginary parts of the amplification factor, which is the amplitude ratio of the lensed wave to the unlensed wave. As a byproduct, we find a new relation that equates integration of the square of the real part of the amplification factor over frequency to that for the imaginary part of the amplification factor. We also obtain a sum rule which relates the integral of the imaginary part of the amplification factor with the magnification of the first arrival image in geometrical optics. Finally, we argue that an incorrect separation of the observed gravitational waveform into the amplification factor and the unlensed waveform generically leads to the violation of the Kramers-Kronig relation. Our work suggests that examining the violation of the Kramers-Kronig relation may be used for correctly extracting the lensing signal in the gravitational wave observations.

We study baryogenesis in a hybrid inflation model which is embedded to the minimal supersymmetric model with right-handed neutrinos. Inflation is induced by a linear combination of the right-handed sneutrinos and its decay reheats the universe. The decay products are stored in conserved numbers, which are transported under the interactions in equilibrium as the temperature drops down. We find that at least a few percent of the initial lepton asymmetry is left under the strong wash-out due to the lighter right-handed (s)neutrinos. To account for the observed baryon number and the active neutrino masses after a successful inflation, the inflaton mass and the Majorana mass scale should be $10^{13}\,{\rm GeV}$ and ${\cal O}(10^{7}$-$10^{10})\,{\rm GeV}$, respectively.

We study the interaction among gravitational and electromagnetic plane waves by means of an analogue electromagnetic model of gravity, where the gravitational properties are encoded in the electromagnetic properties of a material in flat space-time. In this setup, the variations in the metric tensor produced by the gravitational waves are codified as space-time-varying electromagnetic properties. We used an in-housed code based in the finite-difference time domain method to conduct numerical experiments, where we found that, when a monochromatic gravitational plane wave interacts with a narrow-band electromagnetic plane wave, an infinite number of sidebands, equally separated between themselves, are induced by the gravitational wave. Finally, we discuss possible future applications of this effect as an alternative method to directly detect gravitational waves.

We use the Weakly Interacting Massive Particle (WIMP) thermal decoupling scenario to probe Cosmologies in dilatonic Einstein Gauss-Bonnet (dEGB) gravity, where the Gauss-Bonnet term is non-minimally coupled to a scalar field with vanishing potential. We put constraints on the model parameters when the ensuing modified cosmological scenario drives the WIMP annihilation cross section beyond the present bounds from DM indirect detection searches. In our analysis we assumed WIMPs that annihilate to Standard Model particles through an s-wave process. For the class of solutions that comply with WIMP indirect detection bounds, we find that dEGB typically plays a mitigating role on the scalar field dynamics at high temperature, slowing down the speed of its evolution and reducing the enhancement of the Hubble constant compared to its standard value. For such solutions, we observe that the corresponding boundary conditions at high temperature correspond asymptotically to a vanishing deceleration parameter q, so that the effect of dEGB is to add an accelerating term that exactly cancels the deceleration predicted by General Relativity. The bounds from WIMP indirect detection are nicely complementary to late-time constraints from compact binary mergers. This suggest that it could be interesting to use other Early Cosmology processes to probe the dEGB scenario.

We investigate the effects of short axion kination eras on the energy spectrum of the primordial gravitational waves corresponding to modes that re-enter the Hubble horizon at the post-electroweak symmetry breaking epoch, well within the radiation domination era. Our main assumption is the existence of an extremely weakly coupled hidden sector between the Higgs and the axion, materialized by higher order non-renormalizable dimension six and dimension eight operators, active at a scale M of the order 20-100TeV. This new physics scale M which is way higher than the electroweak scale, is motivated by the lack of new particle observations in the large hadron collider to date, beyond the electroweak scale. Once the electroweak symmetry breaking occurs at T\sim GeV, the axion potential acquires a new minimum due to the new terms generated by the electroweak breaking, and the axion oscillations at the origin are destabilized. In effect after some considerable amount of time, the axion rolls swiftly to its new minimum, experiencing a short kination epoch, where its energy density redshifts as $\rho_a\sim a^{-6}$. After it reaches the new minimum, since the latter is energetically less favorable that the Higgs minimum, it decays to the Higgs minimum and the Universe is described again by the Higgs minimum. The axion returns to the origin and commences again oscillations initiated by quantum fluctuations, redshifts as dark matter, and the same procedure is repeated perpetually. These short axion kination eras may disturb the background total equation of state parameter during the radiation domination era. As we show, the energy spectrum of the gravitational waves mainly depends on how many times the short axion kination epochs occur.

We identify a geometric symmetry on the two-flavor Bloch sphere for collective flavor oscillations of a homogeneous dense neutrino gas. Based on this symmetry, analytical solutions to the periodic bipolar flavor evolution are derived. Using numerical calculations, we show that for configurations without this symmetry, the flavor evolution displays deviations from the bipolar flavor motion or even exhibits aperiodic patterns. We also discuss the implication of our finding for more general three-flavor and inhomogeneous cases.

The CNO cycle is one of the fundamental processes of hydrogen burning in stars. The first reaction of the cycle is the radiative proton capture on 12C and the rate of this 12C(p,gamma)13N reaction is related to the 12C/13C ratio observed e.g. in the Solar System. The low-energy cross section of this reaction was measured several times in the past, however, the experimental data are scarce in a wide energy range especially around the resonance at 1.7 MeV. In the present work the 12C(p,gamma)13N cross section was measured between 300 and 1900 keV using the activation method. This method was only used several decades ago in the low-energy region. As the activation method provides the total cross section and has uncertainties different from those of the in-beam gamma-spectroscopy technique, the present results provide a largely independent data set for future low-energy extrapolations and thus for astrophysical reaction rate calculations.

This article aims at clarifying the situation about astrophysical sources that might be observed with haloscope experiments sensitive to gravitational waves in the 1-10 GHz band. The GrAHal setup is taken as a benchmark. We follow a very pedagogical path so that the full analysis can easily be used by the entire community who might not be familiar with the theoretical framework. Different relevant physical regimes are considered in details and some formulas encountered in the literature are revised. The distances that can be probed and expected event rates are carefully evaluated, taking into account degeneracies between physical parameters. We show where experimental efforts should be focused to improve the sensitivity and we conclude that any detection in the near future is extremely unlikely.

Using a heavy-mass effective field theory (HEFT), we study gravitational-wave emission in the scattering of two spinless black holes or neutron stars of arbitrary masses at next-to-leading order in the Post-Minkowskian expansion. We compute the contributions to the one-loop scattering amplitude with four scalars and one graviton which are relevant to the calculation of the waveforms, also presenting expressions of classical tree-level amplitudes with four scalars and up to two radiated gravitons. The latter are obtained using a novel on-shell recursion relation for classical amplitudes with four scalars and an arbitrary number of gravitons. Our one-loop five-point amplitude is expressed in terms of a single family of master integrals with the principal value prescription for linearised massive propagators, which we evaluate using differential equations. In our HEFT approach, soft/heavy-mass expansions of complete integrands are avoided, and all hyper-classical iterations and quantum corrections are dropped at the diagrammatic level, thereby computing directly contributions to classical physics. Our result exhibits the expected factorisation of infrared divergences, the correct soft limits, and highly nontrivial cancellations of spurious poles. Finally, using our amplitude result we compute numerically the corresponding next-to-leading corrections to the spectral waveforms and the far-field time-domain waveforms using the Newman-Penrose scalar $\Psi_4$.

We investigate recent claims that gravitomagnetic effects in linearised general relativity can explain flat and rising rotation curves, such as those observed in galaxies, without the need for dark matter. If one models a galaxy as an axisymmetric, stationary, rotating, non-relativistic and pressureless 'dust' of stars in the gravitoelectromagnetic (GEM) formalism, we show that GEM effects on the circular velocity $v$ of a star are $O(10^{-6})$ smaller than the standard Newtonian (gravitoelectric) effects. Moreover, we find that gravitomagnetic effects are $O(10^{-6})$ too small to provide the vertical support necessary to maintain the dynamical equilibrium assumed. These issues are obscured if one constructs a single equation for $v$, as considered previously. We nevertheless solve this equation for a galaxy having a Miyamoto--Nagai density profile. We show that for the values of the mass, $M$, and semi-major and semi-minor axes, $a$ and $b$, typical for a dwarf galaxy, the rotation curve depends only very weakly on $M$. Moreover, for aspect ratios $a/b > 2$, the rotation curves are concave over their entire range, which does not match observations in any galaxy. Most importantly, we show that for the poloidal gravitomagnetic flux $\psi$ to provide the necessary vertical support, it must become singular at the origin. This originates from the unwitting, but forbidden, inclusion of free-space solutions of the Poisson-like equation that determines $\psi$, hence ruling out the methodology as a means of explaining flat galaxy rotation curves. We further show that recent deliberate attempts to leverage such free-space solutions against the rotation curve problem yield no deterministic modification outside the thin disk approximation, and that, in any case, the homogeneous contributions to $\psi$ are ruled out by the boundary value problem posed by any physical axisymmetric galaxy.