Papers for Tuesday, Aug 01 2023

We report the discovery of ZTF J2020+5033, a high-mass brown dwarf (BD) transiting a low-mass star with an orbital period of 1.90 hours. Phase-resolved spectroscopy, optical and infrared light curves, and precise astrometry from Gaia allow us to constrain the masses, radii, and temperatures of both components with few-percent precision. We infer a BD mass of $M_{\rm BD} = 80.1\pm 1.6\,M_{\rm J}$, almost exactly at the stellar/substellar boundary, and a moderately inflated radius, $R_{\rm BD} = 1.05\pm 0.02\,R_{\rm J}$. The transiting object's temperature, $T_{\rm eff}\approx 1700\,\rm K$, is well-constrained by the depth of the infrared secondary eclipse and strongly suggests it is a BD. The system's high tangential velocity ($v_\perp = 98\,\rm km\,s^{-1}$) and thick disk-like Galactic orbit imply the binary is old; its close distance ($d\approx 140$ pc) suggests that BDs in short-period orbits are relatively common. ZTF J2020+5033 is the shortest-period known transiting BD by more than a factor of 7. Today, the entire binary would comfortably fit inside the Sun. However, both components must have been considerably larger in youth, implying that the orbit has shrunk by at least a factor of $\sim 5$ since formation. The simplest explanation is that magnetic braking continues to operate efficiently in at least some low-mass stars and BDs.

The cores of main sequence intermediate- and high-mass stars are convective. Mixing at the radiative-convective boundary, waves excited by the convection, and magnetic fields generated by convective dynamos all influence the main sequence and post-main sequence evolution of these stars. These effects must be understood to accurately model the structure and evolution of intermediate- and high-mass stars. Unfortunately, there are many challenges in simulating core convection due to the wide range of temporal and spatial scales, as well as many important physics effects. In this review, we describe the latest numerical strategies to address these challenges. We then describe the latest state-of-the-art simulations of core convection, summarizing their main findings. These simulations have led to important insights into many of the processes associated with core convection. Two outstanding problems with multidimensional simulations are, 1. it is not always straightforward to extrapolate from simulation parameters to the parameters of real stars; and 2. simulations using different methods sometimes appear to arrive at contradictory results. To address these issues, next generation simulations of core convection must address how their results depend on stellar luminosity, dimensionality, and turbulence intensity. Furthermore, code comparison projects will be essential to establish robust parameterizations that will become the new standard in stellar modeling.

(Abridged) Using the IRAM NOEMA interferometer, we measures the redshifts of 126 bright galaxies detected in the Herschel H-ATLAS, HeLMS, and HerS surveys. We report reliable spectroscopic redshifts for a total of 124 of the Herschel-selected galaxies. The redshifts are estimated from scans of the 3 and 2-mm bands (and, in one case, the 1-mm band) and are based on the detection of at least two emission lines. Together with the Pilot Programme (Neri et al. 2020), including spectroscopic redshifts of 11 sources, our survey has derived precise redshifts for 135 bright Herschel-selected galaxies, making it the largest sample of high-z galaxies with robust redshifts to date. Most emission lines detected are from 12CO (mainly from J=2-1 to 5-4), with some sources seen in [CI] and H2O emission lines. The spectroscopic redshifts are in the range 0.8<z<6.55 with a median value of z=2.56 +/- 0.10. The line widths of the sources are large, with a mean value for the full width at half maximum Delta(V) of 590 +/- 25 km/s and with 35% of the sources having widths of 700 km/s < Delta(V) < 1800 km/s. Most of the sources are unresolved or barely resolved on scales of 2 to 3 arcsec (or linear sizes of 15-25 kpc, unlensed). Some fields reveal double or multiple sources and, in some cases, sources at different redshifts. Taking these sources into account, there are, in total, 165 individual sources with robust spectroscopic redshifts, including lensed galaxies, binary systems, and over-densities. We present an overview of the z-GAL survey and provide the observed properties of the emission lines, the derived spectroscopic redshifts, and an atlas of the entire sample. The data presented here will serve as a foundation for the other z-GAL papers in this series reporting on the dust emission, the molecular and atomic gas properties, and a detailed analysis of the nature of the sources.

X-shaped Radio Galaxies (XRGs) develop when certain extra-galactic jets deviate from their propagation path. An asymmetric ambient medium (Back-flow model) or complex Active Galactic Nuclei activity (Jet-reorientation model) enforcing the jet direction to deviate may cause such structures. In this context, the present investigation focuses on the modeling of XRGs by performing 3D relativistic magneto-hydrodynamic simulations. We implement different jet propagation models applying an initially identical jet-ambient medium configuration to understand distinctive features. This study, the first of its kind, demonstrates that all adopted models produce XRGs with notable properties, thereby challenging the notion of a universal model. Jet reorientation naturally explains several contentious properties of XRGs, including wing alignment along the ambient medium's primary axis, development of collimated lobes, and the formation of noticeably longer wings than active lobes. Such XRGs disrupt the cluster medium by generating isotropic shocks and channeling more energy than the Back-flow scenario. Our synthetic thermal X-ray maps of the cluster medium reveal four `clear' elongated cavities associated with the wing-lobe alignment, regardless of projection effects, though affecting their age estimation. We show depth and geometric alignment of the evolved cavities may qualify as promising characteristics of XRGs, which may be used to disentangle different formation scenarios.

Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we create N-body models that reproduce the ratio of single-to-binary pulsars in Milky-Way-like GCs. We focus especially on NGC 6752, a typical core-collapsed cluster with many observed millisecond pulsars. Previous studies suggested that an increased rate of neutron star binary disruption in the densest clusters could explain the overabundance of single pulsars in these systems. Here, we demonstrate that binary disruption is ineffective and instead we propose that two additional dynamical processes play the dominant role: (1) tidal disruption of main-sequence stars by neutron stars; and (2) gravitational collapse of heavy white-dwarf-binary merger remnants. Neutron stars formed through these processes may also be associated with fast radio bursts similar to those observed recently in an extragalactic GC.

Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Galaxy. In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: $i$) are predominantly centrally concentrated, with $\sim50\%$ of the stars contained within $5-10$ kpc; $ii$) on average show weak but systematic net rotation in the plane of the host's disc at $z=0$ (i.e., 0.25$\lesssim$$\langle$$\kappa$/$\kappa_{\mathrm{disc}}\rangle\lesssim0.8$); $iii$) present [$\alpha$/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host's old disc; $iv$) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that $\sim$60$\%$ of the proto-Milky Way galaxies are comprised by 1 dominant system (1/5$\lesssim$M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim$4/5) and $4-5$ lower mass systems (M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim$1/10); the other $\sim$40$\%$ are comprised by 2 dominant systems and $3-4$ lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results suggest that large/rich chemical-kinematic-age samples of metal-poor stars in the inner Galaxy should help characterise the different mass fragments of the proto-Milky Way. These data may also help reveal if the Milky Way formed from one or two dominant systems.

(Abridged) We present the dust properties of 125 bright Herschel galaxies selected from the z-GAL survey. The large instantaneous bandwidth of NOEMA provides an exquisite sampling of the underlying dust continuum emission at 2 and 3 mm in the observed frame, with flux densities in at least four side bands for each source. Together with the available Herschel 250, 350, and 500 micron and SCUBA-2 850 micron flux densities, the spectral energy distribution of each source can be analyzed from the far-infrared to the millimeter, with a fine sampling of the Rayleigh-Jeans tail. This wealth of data provides a solid basis to derive robust dust properties, in particular the dust emissivity index, beta, and the dust temperature, T(dust). In order to demonstrate our ability to constrain the dust properties, we used a flux-generated mock catalog and analyzed the results under the assumption of an optically thin and optically thick modified black body emission. For the z-GAL sources, we report a range of dust emissivities with beta ~ 1.5 - 3 estimated up to high precision with relative uncertainties that vary in the range 7% - 15%, and an average of 2.2 +/- 0.3. We find dust temperatures varying from 20 to 50 K with an average of T(dust) ~ 30 K for the optically thin case and ~38 K in the optically thick case. For all the sources, we estimate the dust masses and apparent infrared luminosities (based on the optically thin approach). An inverse correlation is found between T(dust) and beta, which is similar to what is seen in the local Universe. Finally, we report an increasing trend in the dust temperature as a function of redshift at a rate of 6.5 +/- 0.5 K/z for this 500 micron-selected sample. Based on this study, future prospects are outlined to further explore the evolution of dust temperature across cosmic time.

The z-GAL survey observed 137 bright Herschel-selected targets with the IRAM NOrthern Extended Millimeter Array, with the aim to measure their redshift and study their properties. Several of them have been resolved into multiple sources. Consequently, robust spectroscopic redshifts have been measured for 165 individual galaxies in the range 0.8<z<6.5. In this paper we analyse the millimetre spectra of the z-GAL sources, using both their continuum and line emission to derive their physical properties. At least two spectral lines are detected for each source, including transitions of 12CO, [CI], and H2O. The observed 12CO line ratios and spectral line energy distributions of individual sources resemble those of local starbursts. In seven sources the para-H2O(2_11-2_02) transition is detected and follows the IR versus H2O luminosity relation of sub-millimetre galaxies. The molecular gas mass of the z-GAL sources is derived from their 12CO, [CI], and sub-millimetre dust continuum emission. The three tracers lead to consistent results, with the dust continuum showing the largest scatter when compared to 12CO. The gas-to-dust mass ratio of these sources was computed by combining the information derived from 12CO and the dust continuum and has a median value of 107, similar to star-forming galaxies of near-solar metallicity. The same combined analysis leads to depletion timescales in the range between 0.1 and 1.0 Gyr, which place the z-GAL sources between the `main sequence' of star formation and the locus of starbursts. Finally, we derived a first estimate of stellar masses - modulo possible gravitational magnification - by inverting known gas scaling relations: the z-GAL sample is confirmed to be mostly composed by starbursts, whereas ~25% of its members lie on the main sequence of star-forming galaxies (within +/- 0.5 dex).

Among several dark matter candidates, bosonic ultralight (sub meV) dark matter is well motivated because it could couple to the Standard Model (SM) and induce new forces. Previous MICROSCOPE and Eot Wash torsion experiments have achieved high accuracy in the sub-1 Hz region, but at higher frequencies there is still a lack of relevant experimental research. We propose an experimental scheme based on the diamagnetic levitated micromechanical oscillator, one of the most sensitive sensors for acceleration sensitivity below the kilohertz scale. In order to improve the measurement range, we used the sensor whose resonance frequency could be adjusted from 0.1Hz to 100Hz. The limits of the coupling constant are improved by more than 10 times compared to previous reports, and it may be possible to achieve higher accuracy by using the array of sensors in the future.

We model the observed charge states of the elements C, O, Mg, Si, and Fe in the coronal mass ejections (CMEs) ejecta. We concentrate on "halo" CMEs observed in situ by ACE/SWICS to measure ion charge states, and also remotely by STEREO when in near quadrature with Earth, so that the CME expansion can be accurately specified. Within this observed expansion, we integrate equations for the CME ejecta ionization balance, including electron heating parameterized as a fraction of the kinetic and gravitational energy gain of the CME. We also include the effects of non-Maxwellian electron distributions, characterized as a kappa function. Focusing first on the 2010 April 3 CME, we find a somewhat better match to observed charge states with kappa in the range 2-4, close to the theoretical minimum value of kappa = 3/2, implying a hard spectrum of non-thermal electrons. Similar, but more significant results come from the 2011 February 15 event, although it is quite different in terms of its evolution. We discuss the implications of these values, and of the heating required, in terms of the magnetic reconnection Lundquist number and anomalous resistivity associated with CME evolution close to the Sun.

The origin of Mercury's anomalous core and low FeO surface mineralogy are outstanding questions in planetary science. Mercury's composition may result from cosmochemical controls on the precursor solids that accreted to form Mercury. High temperatures and enrichment in solid condensates are likely conditions near the midplane of the inner solar protoplanetary disk. Silicate liquids similar to the liquids quenched in ferromagnesian chondrules are thermodynamically stable in oxygen-rich systems that are highly enriched in dust of CI-chondrite composition. In contrast, the solids surviving into the orbit of Mercury's accretion zone were probably similar to highly unequilibrated, anhydrous, interstellar organic- and presolar grain-bearing chondritic, porous interplanetary dust particles (C-IDPs). Chemical systems enriched in an assumed C-IDP composition dust produce condensates (solid + liquid assemblages in equilibrium with vapor) with super-chondritic atomic Fe/Si ratios at high temperatures, approaching 50% of that estimated for bulk Mercury. Sulfur behaves as a refractory element, but at lower temperatures, in these chemical systems. Stable minerals are FeO-poor, and include CaS and MgS, species found in enstatite chondrites. Disk gradients in volatile compositions of planetary and asteroidal precursors can explain Mercury's anomalous composition, as well as enstatite chondrite and aubrite parent body compositions. This model predicts high sulfur content, and very low FeO content of Mercury's surface rocks.

The recent discovery of enormous Ly$\alpha$ nebulae (ELANe), characterized by physical extents $>200$ kpc and Ly$\alpha$ luminosities $>10^{44}$ erg s$^{-1}$, provide a unique opportunity to study the intergalactic and circumgalactic medium (IGM/CGM) in distant galaxies. Many existing ELANe detections are associated with local overdensities of active galactic nuclei (AGN). We have initiated a search for ELANe around regions containing pairs of quasi-stellar objects (QSOs) using the Palomar Cosmic Web Imager (PCWI). The first study of this search, Cai et al., presented results of ELAN0101+0201 which was associated with a QSO pair at $z=2.45$. In this study, all targets residing in QSO pair environments analyzed have Ly$\alpha$ detections, but only one of the four targets meets the classification criteria of an ELANe associated with a QSO pair region (z$\sim2.87$). The other three sample detections of Ly$\alpha$ nebulae do not meet the size and luminosity criteria to be classified as ELANe. We find kinematic evidence that the ELANe J1613, is possibly powered {mostly by AGN outflows.} The analysis of circularly-averaged surface brightness profiles of emission from the Ly$\alpha$ regions show that the {Ly$\alpha$ emission around $z\sim2$ QSO pairs is consistent with emission around individual QSOs at $z\sim2$, which is fainter than that around $z\sim3$ QSOs. A larger sample of Ly$\alpha$ at z$\sim$2 will be needed to determine if there is evidence of redshift evolution when compared to nebular emissions at z$\sim$3 from other studies.

We continue the study of dipole cosmology framework put forward in \cite{Krishnan:2022qbv}, a beyond FLRW setting that has a preferred direction in the metric which may be associated with a cosmological tilt, a cosmic dipole. In this setup the shear and the tilt can be positive or negative given the dipole direction. Here we focus on the behaviour near the Big Bang (BB) in this setting. We first analyze a single fluid model with a generic constant equation of state $w$. In this case shear and the tilt have the same signs. While details of the behavior near the BB depends on $w$ and the other initial conditions, we find that when the shear is negative we have a shear dominated BB singularity, whereas for a positive shear we have a much milder singularity, the whimper singularity \cite{Ellis:1974ug}, at which the tilt blows up while curvature invariants remain finite. We also analyze dipole $\Lambda$CDM model and explore its near BB behavior, which besides the shear has two tilt parameters, one for radiation and one for the pressureless matter. For positive (negative) shear we again find whimper (curvature) singularity. Moreover, when the tilt parameters have opposite signs, the shear can change sign from negative to positive in the course of evolution of the Universe. We show that the relative tilt of the radiation and the matter generically remains sizeable at late times.

High-resolution JWST observations can test confusion-limited HST observations for a photometric bias that could affect extragalactic Cepheids and the determination of the Hubble constant. We present JWST NIRCAM observations in two epochs and three filters of >330 Cepheids in NGC4258 (which has a 1.5% maser-based geometric distance) and in NGC5584 (host of SNIa 2007af), near the median distance of the SH0ES HST SNIa host sample and with the best leverage among them to detect such a bias. JWST provides far superior source separation from line-of-sight companions than HST in the NIR to largely negate confusion or crowding noise at these wavelengths, where extinction is minimal. The result is a remarkable >2.5x reduction in the dispersion of the Cepheid P-L relations, from 0.45 to 0.17 mag, improving individual Cepheid precision from 20% to 7%. Two-epoch photometry confirmed identifications, tested JWST photometric stability, and constrained Cepheid phases. The P-L relation intercepts are in very good agreement, with differences (JWST-HST) of 0.00+/-0.03 and 0.02+/-0.03 mag for NGC4258 and NGC5584, respectively. The difference in the determination of H_0 between HST and JWST from these intercepts is 0.02+/-0.04 mag, insensitive to JWST zeropoints or count-rate non-linearity thanks to error cancellation between rungs. We explore a broad range of analysis variants (including passband combinations, phase corrections, measured detector offsets, and crowding levels) indicating robust baseline results. These observations provide the strongest evidence yet that systematic errors in HST Cepheid photometry do not play a significant role in the present Hubble Tension. Upcoming JWST observations of >12 SNIa hosts should further refine the local measurement of the Hubble constant.

In the third Gravitational-Wave Transient Catalog (GWTC-3), LIGO-Virgo-KAGRA have observed binary black hole (BBH) mergers out to redshifts $z_\mathrm{merge}\approx1$. Because gravitational waves (GWs) are inefficient at shrinking the binary orbit, some of these BBH systems likely experienced long delay times $\tau$ between the formation of their progenitor stars at $z_\mathrm{form}$ and their GW merger at $z_\mathrm{merge}$. The distribution of delay times predicted by isolated binary evolution resembles a power law $p(\tau)\propto\tau^{\alpha_\tau}$ with slope $-1\lesssim\alpha_\tau\lesssim-0.35$ and a minimum delay time of $\tau_\mathrm{min}=10$ Myr. We use these predicted delay time distributions to infer the formation redshifts of the $\sim70$ BBH events reported in GWTC-3 and measure the BBH progenitor formation rate out to $z_\mathrm{form}\approx3$. For our default $\alpha_\tau=-1$ delay time distribution, we find that GWTC-3 contains at least one system (with 90\% credibility) that formed earlier than $z_\mathrm{form}>3.2$. Comparing our inferred BBH progenitor rate to the star formation rate (SFR), we find that at $z_\mathrm{form}=3$, the number of BBH progenitor systems formed per stellar mass was $4.3^{+4.9}_{-3.6}\times10^{-6}M_\odot^{-1}$ and this yield dropped to $0.2^{+1.7}_{-0.2}\times10^{-6}M_\odot^{-1}$ by $z_\mathrm{form}=0$. We discuss implications of this measurement for the cosmic metallicity evolution, finding that for typical assumptions about the metallicity-dependence of the BBH yield, the average metallicity at $z_\mathrm{form}=3$ was $\langle\log_{10}(Z/Z_\odot)\rangle=-0.2^{+0.3}_{-0.2}$, although the inferred metallicity can vary by a factor of $\approx3$ for different assumptions about the BBH yield.

We investigate the variation in the upper end of stellar initial mass function (uIMF) in 375 young and compact star clusters in five nearby galaxies within $\sim 5$ Mpc. All the young stellar clusters (YSCs) in the sample have ages $\lesssim 4$ Myr and masses above 500 $M_{\odot}$, according to standard stellar models. The YSC catalogs were produced from Hubble Space Telescope images obtained as part of the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. They are used here to test whether the uIMF is universal or changes as a function of the cluster's stellar mass. We perform this test by measuring the H$\alpha$ luminosity of the star clusters as a proxy for their ionizing photon rate, and charting its trend as a function of cluster mass. Large cluster numbers allow us to mitigate the stochastic sampling of the uIMF. The advantage of our approach relative to previous similar attempts is the use of cluster catalogs that have been selected independently of the presence of H$\alpha$ emission, thus removing a potential sample bias. We find that the uIMF, as traced by the H$\alpha$ emission, shows no dependence on cluster mass, suggesting that the maximum stellar mass that can be produced in star clusters is universal, in agreement with previous findings.

Solar eruptions are explosive disruption of coronal magnetic fields, and often launch coronal mass ejections into the interplanetary space. Intriguingly, many solar eruptions fail to escape from the Sun, and the prevailing theory for such failed eruption is based on ideal MHD instabilities of magnetic flux rope (MFR); that is, a MFR runs into kink instability and erupts but cannot reach the height for torus instability. Here, based on numerical MHD simulation, we present a new model of failed eruption in which magnetic reconnection plays a leading role in the initiation and failure of the eruption. Initially, a core bipolar potential field is embedded in a background bipolar field, and by applying shearing and converging motions to the core field, a current sheet is formed within the core field. Then, tether-cutting reconnection is triggered at the current sheet, first slow for a while and becoming fast, driving an erupting MFR. Eventually, the rise of MFR is halted by the downward magnetic tension force of the overlying field, although the MFR apex has well exceeded the critical height of torus instability. More importantly, during the rise of the MFR, it experiences a significant rotation around the vertical axis (with a direction contrary to that predicted by kink instability), rendering the field direction at the rope apex almost inverse to the overlying field. As a result, a strong current sheet is formed between the MFR and the overlying flux, and reconnection occurring in this current sheet ruins completely the MFR.

We compute the contribution from clusters of galaxies to the diffuse neutrino and $\gamma-$ray background. Due to their unique magnetic-field configuration, cosmic rays (CRs) with energy $\leq10^{17}$ eV can be confined within these structures over cosmological time scales, and generate secondary particles, including neutrinos and $\gamma-$rays, through interactions with the background gas and photons. We employ three-dimensional (3D) cosmological magnetohydrodynamical (MHD) simulations of structure formation to model the turbulent intergalactic and intracluster media. We propagate CRs in these environments using multi-dimensional Monte Carlo simulations across different redshifts (from $z \sim 5$ to $z = 0$), considering all relevant photohadronic, photonuclear, and hadronuclear interactions. We also include the cosmological evolution of the CR sources. We find that for CRs injected with a spectral index $1.5 - 2.7$ and cutoff energy $E_\text{max} = 10^{16} - 10^{17}$~eV, clusters contribute to a substantial fraction to the diffuse fluxes observed by the IceCube and Fermi-LAT, and most of the contribution comes from clusters with $M > 10^{14} \, M_{\odot}$ and redshift $z < 0.3$. We also estimated the multimessenger contributions from the local galaxy cluster.

A new class of transient, which has been hypothesized to accompany the explosion of an aspherical compact supernova, would arise when streams of ejecta collide outside the star. However, conditions that favour the prompt release of radiation from the collision, such as a diffuse stellar envelope, disfavour the creation of non-radial ejecta in the first place. To determine whether the collision can both occur and be visible, we simulate aspherical explosions using the HUJI-RICH moving-mesh hydrodynamics code and analyze them in terms of diffusion measures defined for individual fluid elements. While our simulations are highly idealized, they connect to realistic explosions via a single dimensionless parameter. Defining two measures of the importance of diffusivity (two versions of the inverse P'eclet number), we find that one varies in a way that indicates colliding ejecta can release a photon flash, while the other does not. Examining the x-ray transient XT 080109 associated with supernova SN 2008D, we find that its fluence and duration are consistent with the properties of an ejecta collision in the aspherical model that is most likely to emit a flash. Our results give tentative evidence for the possibility of collision-induced flashes for a narrow and radius-dependent range of asphericity, and motivate future radiation hydrodynamics simulations.

Super-Earths within the same close-in, compact planetary system tend to exhibit a striking degree of uniformity in their radius, mass, and orbital spacing, and this 'peas-in-a-pod' phenomenon itself serves to provide one of the strongest constrains on planet formation at large. While it has been recently demonstrated from independent samples that such planetary uniformity occurs for both configurations near and distant from mean motion resonance (MMR), the question thus remains if the strength of this uniformity itself differs between near-resonant and non-resonant configurations such that the two modes may be astrophysically distinct in their evolution. We thus provide in this work a novel comparative size uniformity analysis for 48 near-resonant and 251 non-resonant multi-planet systems from the California Kepler Survey (CKS) catalog, evaluating uniformity both across systems and between planetary pairs within the same system. We find that while multi-planet configurations exhibit strong peas-in-a-pod size uniformity regardless of their proximity to resonance, near-resonant configurations display enhanced intra-system size uniformity as compared to their analogous non-resonant counterparts at the level of both entire systems and sub-system planetary pairs and chains. These results are broadly consistent with a variety of formation paradigms for multiple-planet systems, such as convergent migration within a turbulent protoplanetary disk or planet-planet interactions incited by postnebular dynamical instabilities. Nevertheless, further investigation is necessary to ascertain whether the non-resonant and near-resonant planetary configurations respectively evolve via a singular process or mechanisms that are dynamically distinct.

Context. Certain periodic variations of radial velocities (RV) of wobbling giants originate from exoplanets. Indeed, a number of exoplanets have been discovered around giant stars. Aims. The purpose of our study is to find low-amplitude and long-period RV variations around bright M (super) giants in the RGB (or AGB) stage, which are long-period variables (LPVs) or high-proper-motion (HPM) stars. Methods. High-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) at the Bohyunsan Optical Astronomy Observatory (BOAO) was used to record numerous spectra of nine giants. The observation period for the targets spans 16 years, from 2005 to 2022. Results. We found from the precise RV observations of nine M giants two sub-stellar companions, one with a 28.26$^{+2.05}_{-2.17}$ $M_{J}$ orbiting period of 663.87$^{+4.61}_{-4.31}$ days at a distance of 2.03$^{+0.01}_{-0.01}$ AU (HD 6860) and the other, with a 15.83$^{+2.33}_{-2.74}$ $M_{J}$ orbiting period of 466.63 $^{+1.47}_{-1.28}$ days at a distance of 1.33 $^{+0.08}_{-0.11}$ AU (HD 112300). Our estimate of the stellar parameters for HD 6860 makes it currently the largest star with a sub-stellar companion. We also found RV variations mimicking a planetary companion in HD 18884 and confirmed LPVs in two stars, HD 39801 and HD 42995. The RV variations of some stars seem to be associated with stellar activities rather than reflex orbital motion due to their companions. Such variations are also detected even for HD 6860 and HD 112300, hosting sub-stellar companions.

We present the results of searches for astrophysical neutrinos of few GeV energy from compact binary mergers detected during the first months of the fourth observing run of the LIGO, Virgo, and KAGRA interferometers. We describe our method, based on a selection of $0.5-5$ GeV neutrino events in IceCube, where we search for a statistically significant increase in the number of low-energy candidate events detected around the compact binary merger time. With these results, we constrain neutrino-emitting source populations. Finally, we compare our results with constraints set by neutrino searches at $> 10$ GeV energies and describe the complementarity of these low- and high-energy searches.

Recent studies show the possibility of the formation of fairly regular and global spiral patterns in dwarf galaxies (dS type). Our sample of observed dwarf objects of this class also includes galaxies with a central stellar bar. The analysis of the observational data provides a small rotation velocity and a small disk component mass for dS galaxies, which is in poor agreement with the spiral structure generation mechanism in isolated dwarfs due to the development of disk gravitational instability. Numerical simulation of the stellar-gaseous disks self-consistent dynamics imposes restrictions on the stellar disk thickness and the maximum gas rotation velocity, at which the gravitational mechanism of spiral formation can still be effective.

Stellar occultations currently provide the most accurate ground-based measurements of the positions of natural satellites (down to a few kilometres for the Galilean moons). However, when using these observations in the calculation of satellite ephemerides, the uncertainty in the planetary ephemerides dominates the error budget of the occultation. We quantify the local refinement in the central planet's position achievable by performing Very Long Baseline Interferometry (VLBI) tracking of an in-system spacecraft temporally close to an occultation. We demonstrate the potential of using VLBI to enhance the science return of stellar occultations for satellite ephemerides. We identified the most promising observation and tracking opportunities offered by the Juno spacecraft around Jupiter as perfect test cases, for which we ran simulations of our VLBI experiment. VLBI tracking at Juno's perijove close to a stellar occultation locally (in time) reduces the uncertainty in Jupiter's angular position in the sky to 250-400 m. This represents up to an order of magnitude improvement with respect to current solutions and is lower than the stellar occultation error, thus allowing the moon ephemeris solution to fully benefit from the observation. Our simulations showed that the proposed tracking and observation experiment can efficiently use synergies between ground- and space-based observations to enhance the science return on both ends. The reduced error budget for stellar occultations indeed helps to improve the moons' ephemerides, which in turn benefit planetary missions and their science products, such as the recently launched JUICE and upcoming Europa Clipper missions.

The upcoming JUICE and Europa Clipper missions to Jupiter's Galilean satellites will provide radio science tracking measurements of both spacecraft. Such data are expected to significantly help estimating the moons' ephemerides and related dynamical parameters. However, the two missions will yield an imbalanced dataset, with no flybys planned at Io, condensed over less than six years. Current ephemerides' solutions for the Galilean moons, on the other hand, rely on ground-based astrometry collected over more than a century which, while being less accurate, bring very valuable constraints on the long-term dynamics of the system. An improved solution for the Galilean satellites' complex dynamics could however be achieved by exploiting the existing synergies between these different observation sets. To quantify this, we merged simulated JUICE and Clipper radio science data with existing ground-based astrometric and radar observations, and performed the inversion. Our study specifically focusses on the resulting formal uncertainties in the moons' states, as well as Io's and Jupiter's tidal dissipation parameters. Adding astrometry stabilises the moons' state solution, especially beyond the missions' timelines. It furthermore reduces the uncertainties in $1/Q$ (inverse of the tidal quality factor) by a factor two to four for Jupiter, and about 30-35\% for Io. Among all data types, classical astrometry data prior to 1960 proved particularly beneficial. We also show that ground observations of Io add the most to the solution, confirming that ground observations can fill the lack of radio science data for this specific moon. We obtained a noticeable solution improvement when exploiting the complementarity between all different observation sets. These promising simulation results thus motivate future efforts to achieve a global solution from actual JUICE and Clipper radio science data.

The Chinese Pulsar Timing Array (CPTA) collaboration has recently reported the observational evidence of a stochastic gravitational wave background. In light of the latest CPTA observation, we aim at exploring the ability of CPTA in probing new physics. Specifically, we constrain the first-order cosmological phase transitions with CPTA data, and find that the constraining result is slightly tighter than that of NANOGrav's 12.5-yr data but weaker than NANOGrav's 15-yr data. Considering the possible complexity of gravitational wave sources, we give the constraint on a mixed scenario of cosmological phase transitions and astrophysical supermassive binary black holes. Our analysis suggests that CPTA has a great potential to probe fundamental physics in the near future.

We present an unbiased molecular line survey toward the carbon-rich circumstellar envelope CIT 6 carried out between 90 and 116 GHz with the Arizona Radio Observatory 12 m telescope. A total of 42 lines assigned to 10 molecular species and 4 isotopologues are detected. Despite the absence of any newly identified circumstellar molecules, several transitions are freshly reported for this object. This work is a complement to our previous line survey toward CIT 6 in the frequency ranges of 36-49 GHz and 131-268 GHz. Based on the measurements and in combination with previously published data, we perform a rotation-diagram analysis to determine the column densities, excitation temperatures, and fractional abundances of the molecules. The excitation temperature varies along the radius. The abundance pattern of CIT6 is broadly consistent with that of the prototype C-star IRC+10216 at the specified detection sensitivity, indicating that the molecular richness in IRC+10216 cannot be attributed to an interpretation of unusual chemistry. Nevertheless, subtle distinctions between the two C-stars are found. The higher abundance of carbon-chain molecules and lower 12C/13C ratio in CIT 6 compared to IRC+10216 imply that the former is more massive or in a more evolved phase than the latter.

The apparent close encounters of two satellites in the plane of the sky, called mutual approximations, have been suggested as a different type of astrometric observation to refine the moons' ephemerides. The main observables are the central instants of the close encounters, which have the advantage of being free of any scaling and orientation errors. However, no analytical formulation is available yet for the partials of these central instants, leaving numerical approaches or alternative observables (e.g. derivatives of the apparent distance) as options. Filling that gap, this paper develops an analytical method to include central instants as direct observables in the ephemerides estimation and assesses the quality of the resulting solution. To this end, we ran a covariance analysis to compare the estimated solutions obtained with the two types of observables (central instants versus alternative observables), using the Galilean moons of Jupiter as a test case. Alternative observables can be equivalent to central instants in most cases. However, obtaining consistent solutions between the two observable types requires to accurately weigh each alternative observable individually, based on each mutual approximation's characteristics. In that case, using central instants still yields a small improvement of 10-20% of the formal errors in the radial and normal directions (RSW frame), compared to the alternative observables' solution. This improvement increases for mutual approximations with low impact parameters and large impact velocities. Choosing between the two observables thus requires careful assessment based on the characteristics of the available observations. Using central instants over alternative observables ensures that the state estimation fully benefits from the information encoded in mutual approximations, which can be desirable for certain applications.

Using simultaneous multi-filter observations during the transit of an exoplanet around a K dwarf star, we determine the temperature of a starspot through modeling the radius and position with wavelength-dependent spot contrasts. We model the spot using the starspot modeling program STarSPot (STSP), which uses the transiting companion as a knife-edge probe of the stellar surface. The contrast of the spot, i.e. the ratio of the integrated flux of a darker spot region to the star's photosphere, is calculated for a range of filters and spot temperatures. We demonstrate this technique using simulated data of HAT-P-11, a K dwarf (T = 4780 K) with well-modeled starspot properties for which we obtained simultaneous multi-filter transits using LCO's MuSCAT3 instrument on the 2-meter telescope at Haleakala Observatory which allows for simultaneous, multi-filter, diffuser assisted high-precision photometry. We determine the average (i.e. a combination of penumbra and umbra) spot temperature for HAT-P-11's spot complexes is 4500 K $\pm$ 100 K using this technique. We also find for our set of filters that comparing the SDSS g' and i' filters maximizes the signal difference caused by a large spot in the transit. Thus, this technique allows for the determination of the average spot temperature using only one spot occultation in transit and can provide simultaneous information on the spot temperature and spot properties.

The gravitational waves produced by binary neutron star mergers offer a unique window into matter behavior under extreme conditions. In this context, we model analytically the effect of matter on the gravitational waves from binary neutron star mergers. We start with a binary black hole system, leveraging the post-Newtonian formalism for the inspiral and the Backwards-one-Body model for the merger. We combine the two methods to generate a baseline waveform and we validate our results against numerical relativity simulations. Next, we integrate tidal effects in phase and amplitude to account for matter and spacetime interaction, by using the NRTidal model, and test its accuracy against numerical relativity predictions, for two equations of state, finding a mismatch around the merger. Subsequently, we lift the restriction on the coefficients to be independent of the tidal deformability, and recalibrate them using the numerical relativity predictions. We obtain better fits for phase and amplitude around the merger, and are able to extend the phase modeling beyond the merger. We implement our method in a new open-source Python code, steered by a Jupyter Notebook. Our research offers new perspectives on analytically modeling the effect of tides on the gravitational waves from binary neutron star mergers.

We report an updated value for the distance to the planetary nebula NGC 6309 (the Box Nebula). The distance is found through two Kinematic Distance Methods (KDMs): the system of two equations reported in Zhu et al. 2013 and the Monte Carlo method reported by Wenger et al. 2018. We find the kinematic distance to NGC 6309 to be 4.1 kpc with an upper uncertainty of +0.29 kpc and a lower uncertainty of -0.38 kpc. We also calculate the distance to Cassiopeia A with the two KDMs and compare to the value reported by Reed et al 1995. The Zhu et al. method and Wenger et al. method yield a value within thirty percent and twenty percent of the Reed et al. method, respectively. The value reported by Reed et al 1995 was contained within the error bounds produced by the Wenger et al. method. The distance measurement to Cassiopeia A suggests that both KDMs, while imperfect, are moderately accurate methods for determining the distance to NGC objects in the plane of the Milky Way.

We investigate the clustering properties of 3<z<5 candidate passive galaxies from the Merlin et al. (2019) sample residing in the GOODS-North (35 sources) and GOODS-South (33 sources) fields. Within the large uncertainties due to the paucity of sources we do not detect clustering signal in GOODS-North, while this is present in GOODS-South, highlighting the importance of the effects of cosmic variance. The estimated correlation length in GOODS-South is r_0=12^+4_-5 Mpc, while the estimated minimum mass for a halo capable to host one of such high-redshift quenched galaxies is log10(M_min/M_sun) =13.0\pm 0.3, once also the constraints from their space density are taken into account. Both values are compatible with the results from GOODS-North. Putting the above findings in a cosmological context, these suggest no evolution of the dark matter content of the hosts of passive galaxies during the past 12.5 Gyr, i.e. during more than 90% of the age of the Universe. We discuss possible scenarios for the observed trend.

Choked gamma-ray bursts (CGRBs) are possible neutrino sources that have been proposed as capable of generating the flux detected by IceCube, since no accompanying gamma-ray signal is expected, as required by observations. We focus on obtaining the neutrino flux and flavor composition corresponding to CGRBs under different assumptions for the target photon density and the magnetic field of the emission region. We consider the injection of both electrons and protons into the internal shocks of CGRBs, and using a steady-state transport equation, we account for all the relevant cooling processes. In particular, we include the usually adopted background of soft photons, which is the thermalized emission originated at the shocked jet head. Additionally, we consider the synchrotron photons emitted by the electrons co-accelerated with the protons at the internal shocks in the jet. We also obtain the distributions of pions, kaons, and muons using the transport equation to account for the cooling effects due not only to synchrotron emission but also interactions with the soft photons in the ambient. We integrate the total diffuse flux of neutrinos of different flavors and compute the flavor ratios to be observed on Earth. As a consequence of the losses suffered mainly by pions and muons, we find these ratios to depend on the energy: for energies above ~(10^5-10^6) GeV (depending on the magnetic field, proton-to-electron ratio, and jet power), we find that the electron flavor ratio decreases and the muon flavor ratio increases, while the tau flavor ratio increases only moderately. Our results are sensitive to the mentioned key physical parameters of the emitting region of CGRBs. Hence, the obtained flavor ratios are to be contrasted with cumulative data from ongoing and future neutrino instruments in order to assess the contribution of these sources to the diffuse flux of astrophysical neutrinos.

A likely starburst galaxy (SBG), IRAS 13052-5711, which is the most distant SBG candidate discovered to date, was found by analyzing 14.4 years of data from the Fermi large-area telescope (Fermi-LAT). This SBG's significance level is approximately 6.55$\sigma$ in the 0.1-500 GeV band. Its spatial position is close to that of 4FGL J1308.9-5730, determined from the Fermi large telescope fourth-source Catalog (4FGL). Its power-law spectral index is approximately 2.1, and its light curve (LC) for 14.4 years has no significant variability. These characteristics are highly similar to those of SBGs found in the past. We calculate the SBG's star formation rate (SFR) to be 29.38 $\rm M_{\odot}\ yr^{-1}$, which is within the SFR range of SBGs found to date. Therefore, IRAS 13052-5711 is considered to be a likely SBG. In addition, its 0.1-500 GeV luminosity is (3.28 $\pm$ 0.67) $\times 10^{42}\ \rm erg\ s^{-1}$, which deviates from the empirical relationship of the $\gamma$-ray luminosity and the total infrared luminosity. We considered a hadronic model to explain the GeV spectrum of IRAS 13052-5711.

The megajansky radio burst, FRB 20200428, and other bright radio bursts detected from the Galactic source SGR J1935+2154 suggest that magnetars can make fast radio bursts (FRBs), but the emission site and mechanism of FRB-like bursts are still unidentified. Here we report the emergence of a radio pulsar phase of the magnetar five months after FRB 20200428. 795 pulses were detected in 16.5 hours over 13 days by the Five-hundred-meter Aperture Spherical Radio telescope, with luminosities about eight decades fainter than FRB 20200428. The pulses were emitted in a narrow phase window anti-aligned with the X-ray pulsation profile observed by the X-ray telescopes. The bursts, conversely, appear in random phases. This dichotomy suggests that radio pulses originate from a fixed region within the magnetosphere, but bursts occur in random locations and are possibly associated with explosive events in a dynamically evolving magnetosphere. This picture reconciles the lack of periodicity in cosmological repeating FRBs within the magnetar engine model.

Astrophysical transient events like Gamma Ray Bursts (GRBs) have always been promising candidates for multi-messenger astronomy, with electromagnetic and gravitational wave signals having already been observed in GRBs such as GRB 170817A. The neutrino signatures of these bursts have been long-awaited as well, with many models predicting different spectra. Most of these searches have been in the hundreds of GeV to PeV range. However, as different models indicate a possible lower energy neutrino signal, we intend to expand this search to the lowest limits of IceCube (0.5-5 GeV) as well. With the plan to look at more transient events, we present the result of the first IceCube search for < 5 GeV astrophysical neutrinos emitted from a GRB, for GRB 221009A; the brightest GRB ever observed. Furthermore, we present plans to improve the observations of < 5 GeV neutrinos in IceCube, with which we plan to probe more transient events in the future. These improvements include the addition of direction reconstruction at these energies, and optimization of the noise rejection. With these improvements, GRB 221009A is just the start of the low-energy neutrino search from transient events with IceCube.

Recent pulsar timing array (PTA) experiments have reported strong evidence of the stochastic gravitational wave background (SGWB). If interpreted as primordial gravitational waves (GWs), the signal favors a strongly blue-tilted spectrum. Consequently, the nonsingular cosmology, which is able to predict a strongly blue-tilted GW spectrum with $n_T \simeq 2$ on certain scales, offers a potential explanation for the observed SGWB signal. In this paper, we present a Genesis-inflation model capable of explaining the SGWB signal observed by the PTA collaborations while also overcoming the initial singularity problem associated with the inflationary cosmology. Furthermore, our model predicts distinctive features in the SGWB spectrum, which might be examined by forthcoming space-based gravitational wave experiments.

The launch of NASA's Kepler space telescope in 2009 revolutionized the quality and quantity of observational data available for asteroseismic analysis. While Kepler was able to detect solar-like oscillations in hundreds of main-sequence and subgiant stars, the Transiting Exoplanet Survey Satellite (TESS) is now making similar observations for thousands of the brightest stars in the sky. The Asteroseismic Modeling Portal (AMP) is an automated and objective stellar model-fitting pipeline for asteroseismic data, which was originally developed to use models from the Aarhus Stellar Evolution Code (ASTEC). We briefly summarize an updated version of the AMP pipeline that uses Modules for Experiments in Stellar Astrophysics (MESA), and we present initial modeling results for the Sun and several solar analogs to validate the precision and accuracy of the inferred stellar properties.

The physics-based method of forecasting the high energy solar flares (SFs) with the help of the neutrino detectors utilizing the coherent elastic neutrino-nucleus scattering (CE$\nu$NS) is considered. The behavior of the neutrino beams traveling through the coupled sunspots (CSs) being the sources of the future SFs is investigated. Two neutrino beams are included into consideration, namely, the $\nu_{eL}$ beam and the $\nu_{\mu L}$ beam which have been produced after the passage of the Micheev-Smirnov-Wolfenstein resonance. It is assumed that the neutrinos possess the charge radius, the magnetic and anapole moments while the CS magnetic field is vortex, nonhomogeneous and has twisting. Estimations of the weakening of the neutrino beams after traversing the resonant layers are given. It is demonstrated, that this weakening could be registered by the detectors employing CE$\nu$NS only when the neutrinos have the Dirac nature.

The study of planetary atmospheres is critical to our understanding of the origin and evolution of the Solar System. The combined effect of various physical and chemical processes over billions of years have resulted in a variety of planetary atmospheres across the Solar System. This paper performs a comparative study of planetary atmospheres and their engineering implications for future entry and aerocapture missions. The thick Venusian atmosphere results in high deceleration and heating rates and presents a demanding environment for both atmospheric entry and aerocapture. The thin Martian atmosphere allows low aerodynamic heating, but itself is not enough to decelerate a lander to sufficiently low speeds for a soft landing. With their enormous gravity wells, Jupiter and Saturn entry result in very high entry speeds, deceleration, and heating making them the most demanding destinations for atmospheric entry and impractical for aerocapture. Titan is a unique destination, with its low gravity and greatly extended thick atmosphere enabling low deceleration and heating loads for entry and aerocapture. Uranus and Neptune also have large gravity wells, resulting in high entry speeds, high deceleration and heating compared to the inner planets, but are still less demanding than Jupiter or Saturn.

We present a common unifying macroscopic framework for precursors in relativistic shock waves. These precursors transfer energy and momentum from the hot downstream to the cold upstream, modifying the shock structure. Derishev & Piran (2016} have shown that in a steady state, there is a maximal fraction of the downstream energy flux that the precursor can carry. We show that at this critical value, the shock disappears, and the flow passes through a sonic point. This behavior resembles the classical Newtonian Rayleigh flow problem. At the critical value, the transition is unstable as perturbations in the upstream accumulate at the sonic point. Thus, if such a point is reached, the shock structure is drastically modified, and the flow becomes turbulent.

Probability distributions become non-Gaussian when the flat $\Lambda$CDM model is fitted to redshift binned data in the late Universe. We explain mathematically why this non-Gaussianity arises and confirm that Markov Chain Monte Carlo (MCMC) marginalisation leads to biased inferences in observational Hubble data (OHD). In particular, in high redshift bins we find that $\chi^2$ minima, as identified from both least squares fitting and the MCMC chain, fall outside of the $1 \sigma$ confidence intervals. We resort to profile distributions to correct this bias. Doing so, we observe that $z \gtrsim 1$ cosmic chronometer (CC) data currently prefers a non-evolving (constant) Hubble parameter over a Planck-$\Lambda$CDM cosmology at $\sim 2 \sigma$. We confirm that both mock simulations and profile distributions agree on this significance. Moreover, on the assumption that the Planck-$\Lambda$CDM cosmological model is correct, using profile distributions we confirm a $> 2 \sigma$ discrepancy with Planck-$\Lambda$CDM in a combination of CC and baryon acoustic oscillations (BAO) data beyond $ z \sim 1.5$ that was noted earlier through comparison of least square fits of observed and mock data.

(Abridged) High-mass star formation is far less understood than low-mass star formation. It entails molecular outflows, which disturb the protostellar clump. Studying these outflows and the shocked gas they cause is key for a better understanding of this process. This study aims to characterise the behaviour of molecular outflows in the most massive protostellar sources in the Southern Galaxy by looking for evolutionary trends and associating shocked gas with outflow activity. We present APEX SEPIA180 observations (beamwidth $\sim$36") of SiO outflow candidates of a sample of 32 luminous and dense clumps, candidates to harbouring Hot Molecular Cores. We study the SiO(4-3) line emission, an unambiguous tracer of shocked gas and recent outflow activity, the HCO$^+$(2-1) and H$^{13}$CO$^+$(2-1) lines. 78% of our sample present SiO emission. Nine of these also have wings in the HCO$^+$ line, indicating outflow activity. The SiO emission of these 9 sources is more intense and wider than the rest, suggesting that the outflows in this group are faster and more energetic. Three positive correlations between the outflow properties were found, which suggest that more energetic outflows bear to mobilise more material. No correlation was found between the evolutionary stage indicator $L/M$ and SiO outflow properties, supporting that outflows happen throughout the whole high-mass star formation process. We conclude that sources with both SiO emission and HCO$^+$ wings and sources with only SiO emission are in virtually the same advanced stage of evolution in the high-mass star formation process. The former present more massive and more powerful SiO outflows than the latter. Thus, looking for more outflow signatures such as HCO$^+$ wings could help identify more massive and active massive star-forming regions in samples of similarly evolved sources, as well as sources with older outflow activity.

GRB 221009A is the brightest Gamma Ray Burst (GRB) ever observed. The observed extremely high flux of high and very-high-energy photons provide a unique opportunity to probe the predicted neutrino counterpart to the electromagnetic emission. We have used a variety of methods to search for neutrinos in coincidence with the GRB over several time windows during the precursor, prompt and afterglow phases of the GRB. MeV scale neutrinos are studied using photo-multiplier rate scalers which are normally used to search for galactic core-collapse supernovae neutrinos. GeV neutrinos are searched starting with DeepCore triggers. These events don't have directional localization, but instead can indicate an excess in the rate of events. 10 GeV - 1 TeV and >TeV neutrinos are searched using traditional neutrino point source methods which take into account the direction and time of events with DeepCore and the entire IceCube detector respectively. The >TeV results include both a fast-response analysis conducted by IceCube in real-time with time windows of T$_0 - 1$ to T$_0 + 2$ hours and T$_0 \pm 1$ day around the time of GRB 221009A, as well as an offline analysis with 3 new time windows up to a time window of T$_0 - 1$ to T$_0 + 14$ days, the longest time period we consider. The combination of observations by IceCube covers 9 orders of magnitude in neutrino energy, from MeV to PeV, placing upper limits across the range for predicted neutrino emission.

Among more than 800 known fast radio bursts (FRBs), only two, namely FRB 20121102A and FRB 20190520B, are confirmed to be associated with a persistent radio sources (PRS). Here we report evidence of apparent temporal variability in the PRS associated with the bursting FRB 20190520B based on the Karl G. Jansky Very Large Array (VLA) observations taken in 2020 and 2021. Based on the analysis of epoch-to-epoch variability of the PRS at L, S, C, and X band in 1-12 GHz, we detected not only overall marginal variability but also a likely radio flux decrease between the observations taken in 2020 and 2021 at 3 GHz. Assuming no spectral variation in the PRS during these observations, we found the evidence for an overall broadband radio flux decrease by about 20 percent between the 2020 and the 2021 observations, suggesting that the PRS probably evolves on the yearly time scale. If we attribute the marginal variability at 3 GHz as intrinsic or due to scintillation, the size of potential variable component of the PRS is constrained to be sub-parsec. On the other hand, the size of the PRS can be also constrained to be larger than about 0.22 parsec from the averaged radio spectrum and the integrated radio luminosity in the 1-12 GHz band based on equipartition and self-absorption arguments. We discuss potential origins of the PRS and suggest that a scenario of accreting compact objects might be able to explain the temporal and spectral properties of the PRS. Confirmation of variability or flux decline of the PRS would be critical to our understanding of the PRS and its relation to the bursting source.

We analyzed four epochs of beamformed EVN data of the Crab Pulsar at 1658.49 MHz. With the high sensitivity resulting from resolving out the Crab Nebula, we are able to detect even the faint high-frequency components in the folded profile. We also detect a total of 65951 giant pulses, which we use to investigate the rates, fluence, phase, and arrival time distributions. We find that for the main pulse component, our giant pulses represent about 80% of the total flux. This suggests we have a nearly complete giant pulse energy distribution, although it is not obvious how the observed distribution could be extended to cover the remaining 20% of the flux without invoking large numbers of faint bursts for every rotation. Looking at the difference in arrival time between subsequent bursts in single rotations, we confirm that the likelihood of finding giant pulses close to each other is increased beyond that expected for randomly occurring bursts - some giant pulses consist of causally related microbursts, with typical separations of $\sim\!30{\rm\;\mu s}$ - but also find evidence that at separations $\gtrsim\!100{\rm\;\mu s}$ the likelihood of finding another giant pulse is suppressed. In addition, our high sensitivity enabled us to detect weak echo features in the brightest pulses (at $\sim\!0.4\%$ of the peak giant pulse flux), which are delayed by up to $\sim\!300{\rm\;\mu s}$.

IceCube DeepCore is an extension of the IceCube Neutrino Observatory designed to measure GeV scale atmospheric neutrino interactions for the purpose of neutrino oscillation studies. Distinguishing muon neutrinos from other flavors and reconstructing inelasticity are especially difficult tasks at GeV scale energies in IceCube DeepCore due to sparse instrumentation. Convolutional neural networks (CNNs) have been found to have better success at neutrino event reconstruction than conventional likelihood-based methods. In this contribution, we present a new CNN model that exploits time and depth translational symmetry in IceCube DeepCore data and present the model's performance, specifically for flavor identification and inelasticity reconstruction.

Based on the \textit{Gaia} DR3 RR Lyrae catalog, we use two methods to fit the density profiles with an improved broken power law, and find that there are two break radii coinciding with the two apocenter pile-ups of high-eccentricity Gaia-Sausage-Enceladus (GSE) merger. Also, there is a break caused by the Sagittarius (Sgr) stream. Combining the positions of all breaks, we briefly analyze the metallicity and its dispersion as a function of $r$ as well as its distribution in cylindrical coordinates. For the clean sample, the $z\text{-to-}x$ ellipsoid axial ratio $q$ in $36\,{\rm kpc}\,\textless\,r\,\textless\,96\,{\rm kpc}$ becomes much smaller than that of the inner halo $(r\,\textless\,36\,{\rm kpc})$, while the major axis has a large uncertainty in the region of $36-66\,{\rm kpc}$ and the one in the region of $66-96\,{\rm kpc}$ is obviously different from that dominated by the Hercules-Aquila Cloud (HAC) and the Virgo Overdensity (VOD) in the inner halo, which indicates that there is an over-density structure distributed at low zenithal angles. Finally, we found that the over-density structure in the outer halo ($r\,\textgreater\,50\,{\rm kpc}$) is shell-shaped and relatively metal-rich compared to the outer background halo. We conclude that the shells could be the apocenter pile-ups of the high-eccentricity GSE merger, which is supported by previous numerical simulations.

Blue metal-poor stars are main-sequence stars that are bluer and brighter than typical turn-off stars in metal-poor globular clusters. They are thought to have either evolved through post-mass transfer mechanisms as field blue straggler stars or have accreted from Milky Way dwarf satellite galaxies. It has been found that a considerable fraction of blue metal poor stars are binaries, possibly with a compact companion. We observed 27 blue metal poor stars using UV imaging telescope of AstroSat in two far-UV filters, F148W and F169M. In this work, we explain the possible formation channels of two stars, BMP17 and BMP37. We fit BMP17 with a single-component spectral energy distribution whereas BMP37 with a binary-component spectral energy distribution. As both of them are known SB1s, we suggest that the WD companion of BMP17 may have cooled down so that it is out of UV imaging telescope detection limit. On the other hand, we discover a normal mass white dwarf as the hot companion of BMP37, indicating mass transfer as the possible formation channel

Validating the accelerating expansion of the universe is an important issue for understanding the evolution of the universe. By constraining the cosmic acceleration parameter $X_H$, we can discriminate between the $\Lambda \mathrm{CDM}$ (cosmological constant plus cold dark matter) model and LTB (the Lema\^itre-Tolman-Bondi) model. In this paper, we explore the possibility of constraining the cosmic acceleration parameter with the inspiral gravitational waveform of neutron star binaries (NSBs) in the frequency range of 0.1Hz-10Hz, which can be detected by the second-generation space-based gravitational wave detector DECIGO. We use a convolutional neural network (CNN), a long short-term memory (LSTM) network combined with a gated recurrent unit (GRU), and Fisher information matrix to derive constraints on the cosmic acceleration parameter $X_H$. Based on the simulated gravitational wave data with a time duration of 1 month, we conclude that CNN can limit the relative error to 14.09%, while LSTM network combined with GRU can limit the relative error to 13.53%. Additionally, using Fisher information matrix for gravitational wave data with a 5-year observation can limit the relative error to 32.94%. Compared with the Fisher information matrix method, deep learning techniques will significantly improve the constraints on the cosmic acceleration parameters at different redshifts. Therefore, DECIGO is expected to provide direct measurements of the acceleration of the universe, by observing the chirp signals of coalescing binary neutron stars.

In this paper, we present a cosmology-independent method to constrain cosmological models from the latest 221 gamma-ray bursts (GRBs) sample, including 49 GRBs from Fermi catalog with the Amati relation (the $E_{\rm p}$-${E}_{\rm iso}$ correlation), which calibrated by using a Gaussian process from the Pantheon+ type Ia supernovae (SNe Ia) sample. With 182 GRBs at $0.8\le z\le8.2$ in the Hubble diagram and the latest observational Hubble data (OHD) by the Markov Chain Monte Carlo (MCMC) method, we obtained $\Omega_{\rm m}$ = $0.348^{+0.048}_{-0.066}$ and $h$ = $0.680^{+0.029}_{-0.029}$ for the flat $\Lambda$CDM model, and $\Omega_{\rm m}$ = $0.318^{+0.067}_{-0.059}$, $h$ = $0.704^{+0.055}_{-0.068}$, $w$ = $-1.21^{+0.32}_{-0.67}$ for the flat $w$CDM model. These results are consistent with those in which the coefficients of the Amati relation and the cosmological parameters fitted simultaneously.

Phenomenological models are widely used in cosmology in relation to constraining different cosmological models, with two common examples being cosmographic expansions and modeling the equation-of-state parameter of dark energy. This work presents a study of how using different phenomenological expressions for observables and physical quantities versus using physically motivated, derived expressions affects cosmological parameter constraints. The study includes the redshift-distance relation and Hubble parameter as observables, and the dark energy equation-of-state parameter as a physical quantity, and focuses on constraining the cosmological parameter $\Omega_{\Lambda}$. The observables and equation-of-state parameter are all modeled both using the physical, derived expressions and a variety of phenomenological models with different levels of accuracy and complexity. The results suggest that the complexity of phenomenological expressions only has minor impact on the parameter constraints unless the complexity is very high. The results also indicate that statistically significantly different results can be expected from parameter constraints using different phenomenological models if the models do not have very similar accuracy. This suggests that a good practice is to use multiple phenomenological models when possible, in order to assess the model dependence of results. Straightforward examples of this is that results obtained using cosmographic expansions should always be checked against similar results obtained with expansions of other order, and when using phenomenological models such as for the equation-of-state parameters, robustness of results could be assessed using fitted models from symbolic regression, similar to what is done in this study.

We derive general expressions for how the Alcock-Paczynski distortions affect the power spectrum and the bispectrum of cosmological fields. We compute explicit formulas for the mixing coefficients of bispectrum multipoles in the linear approximation. The leading-order effect for the bispectrum is the uniform dilation of all three wavevectors. The mixing coefficients depend on the shape of the bispectrum triplet. Our results for the bispectrum multipoles are framed in terms of the "natural" basis of the lengths of three wavevectors but can be easily generalized for other bases and reduction schemes. Our validation tests confirm that the linear approximation is extremely accurate for all power spectrum multipoles. The linear approximation is accurate for the bispectrum monopole but results in sub-percent level inaccuracies for the bispectrum quadrupole and fails for the bispectrum hexadecapole. Our results can be used to simplify the analysis of the bispectrum from galaxy surveys, especially the measurement of the Baryon Acoustic Oscillation peak position. They can be used to replace numeric schemes with exact analytic formulae.

The recent extremely bright GRB 230307A from a binary neutron star merger may offer a good probe for the production of GRB-neutrinos. Within the constraint from IceCube neutrino non-detection, the limitations for key physical parameters of this burst are extracted in different scenarios including the fireball, Poynting-flux-dominated (PFD) and hybrid jet. Different from the former nearby `monsters' and due to its smaller isotropic equivalent radiated energy ($E_{\gamma,\rm iso}\sim4\times10^{52}$ erg), the constraint seems loose if non-thermal neutrinos produced from photomeson interactions are the only consideration. However, a quasi-thermal neutrino emission from hadronuclear processes is constrained in this neutron-rich post-merger environment, and the upper limit of the allowed nucleon loading factor is about a few. Based on this, a discussion is presented on the possible prompt emission mechanism and jet composition for GRB 230307A in the context of multi-messenger astrophysics.

The Sunrise Chromospheric Infrared spectroPolarimeter (SCIP) has been developed for the third flight of the SUNRISE balloon-borne stratospheric solar observatory. The aim of SCIP is to reveal the evolution of three-dimensional magnetic fields in the solar photosphere and chromosphere using spectropolarimetric measurements with a polarimetric precision of 0.03\% (1$\sigma$). Multiple lines in the 770 and 850 nm wavelength bands are simultaneously observed with two 2k$\times$2k CMOS cameras at a frame rate of 31.25 Hz. Stokes profiles are calculated onboard by accumulating the images modulated by a polarization modulation unit, and then compression processes are applied to the two-dimensional maps of the Stokes profiles. This onboard data processing effectively reduces the data rate. SCIP electronics can handle large data formats at high speed. Before the implementation into the flight SCIP electronics, a performance verification of the onboard data processing was performed with synthetic SCIP data that were produced with a numerical simulation modeling the solar atmospheres. Finally, we verified that the high-speed onboard data processing was realized on ground with the flight hardware by using images illuminated by natural sunlight or an LED.

In this population synthesis work we study a variety of possible origin channels of supernovae type Ia (SNe Ia) Among them mergers of carbon-oxygen (CO) and oxygen-neon (ONe) white dwarfs (WDs) under the influence of gravitational waves are considered as the primary channel of SNe Ia formation. We estimated frequencies of mergers of WDs with different chemical compositions and distributions of masses of merging WDs. We computed the dependence of the ratio of merger frequencies of ONe and CO WDs as primaries in corresponding binaries on time. The scatter of masses of considered sources (up to the factor $1.5-2$) of SNe Ia is important and should be carefully studied with other sophisticated methods from theoretical point of view. Our ``game of parameters'' potentially explains the increased dimming of SNe Ia in the redshift range $z\approx 0.5-1$ by the changes in the ratio of ONe and CO WDs, i.e., to describe the observed accelerated expansion of the Universe in terms of the evolution of properties of SNe Ia instead of cosmological explanations. This example shows the extreme importance of theoretical studies of problems concerning SNe Ia, because evolutionary scenario and parameter games in nature potentially lead to confusions in their empirical standardization and, therefore, they can influence on cosmological conclusions.

Classical T Tauri stars (cTTs) accrete from their circumstellar disk. The material falls onto the stellar surface, producing an accretion shock, which generates veiling in a star's spectra. In addition, the shock causes a localized accretion spot at the level of the chromosphere. Our goal is to investigate the accretion, particularly the mass accretion rates (Macc), for the cTTs DK Tau, over two periods of 17 and 29 days, using two different procedures for comparison purposes. The first method relies on the derivation of the accretion luminosity via accretion-powered emission lines. The second compares the variability of the optical veiling with accretion shock models to determine mass accretion rates. We used observations taken in 2010 and 2012 with the ESPaDOnS spectropolarimeter at the CFHT. We find peak values of the veiling (at 550 nm) ranging from 0.2 to 1.3, with a steeper trend across the wavelength range for higher peak values. When using the accretion-powered emission lines, we find mass accretion rate values ranging from log(Macc[Msol/yr]) = -8.20 to log(Macc[Msol/yr]) = -7.40. This agrees with the values found in the literature, as well as the values calculated using the accretion shock models and the veiling. In addition, we identify a power-law correlation between the values of the accretion luminosity and the optical veiling. For the 2010 observations, using the values of the filling factors (which represent the area of the star covered by an accretion spot) derived from the shock models, we infer that the accretion spot was located between +45 degrees and +75 degrees in latitude. We show that both methods of determining the mass accretion rate yield similar results. We also present a helpful means of confirming the accretion luminosity values by measuring the veiling at a single wavelength in the optical.

The accretion flow / jet correlation in neutron star (NS) low-mass X-ray binaries (LMXBs) is far less understood when compared to black hole (BH) LMXBs. In this paper we will present the results of a dense multi-wavelength observational campaign on the NS LMXB 4U 1820-30, including X-ray (Nicer, NuSTAR and AstroSAT) and quasi-simultaneous radio (ATCA) observations in 2022. 4U 1820-30 shows a peculiar 170 day super-orbital accretion modulation, during which the system evolves between "modes" of high and low X-ray flux. During our monitoring, the source did not show any transition to a full hard state. X-ray spectra were well described using a disc blackbody, a Comptonisation spectrum along with a Fe K emission line at 6.6 keV. Our results show that the observed X-ray flux modulation is almost entirely produced by changes in the size of the region providing seed photons for the Comptonisation spectrum. This region is large (about 15 km) in the high mode and likely coincides with the whole boundary layer, while it shrinks significantly (<10 km) in low mode. The electron temperature of the corona and the observed RMS variability in the hard X-rays also exhibit a slight increase in low mode. As the source moves from high to low mode, the radio emission due to the jet becomes about 5 fainter. These radio changes appear not to be strongly connected to the hard-to-soft transitions as in BH systems, while they seem to be connected mostly to variations observed in the boundary layer.

In the innermost regions of Active Galactic Nuclei (AGN), matter is understood to be flowing onto the Supermassive black hole (SMBH), which forms an accretion disk. This disk is responsible for the optical/UV continuum emission observed in the spectra of AGN. Reverberation Mapping of the accretion disk using multiple bands can yield the structure of the disk. The emission is expected to be of the black body type peaking at different wavelengths. Hence, depending on the temperature of the disk, continuous, simultaneous monitoring in multiple wavelength ranges to cover hotter inner regions and cooler outer regions can yield the structure and temperature profile of the accretion disk itself. In this study, we present initial results from our accretion disk reverberation mapping campaign targeting AGN with Super High Eddington Accreting Black Holes (SEAMBH). Our analysis on one of the sources- IRAS 04416+1215; based on the broadband observations using the Growth India telescope (GIT), reveals that the size of the accretion disk for this source, calculated by cross-correlating the continuum light curves is larger than expected from the theoretical model. We fit the light curves directly using the thin disk model available in {\sc javelin} and find that the disk sizes are approximately 4 times larger than expected from the Shakura Sunyaev (SS) disk model. Further studies are needed to understand better the structure and physics of AGN accretion disks and their role in the evolution of galaxies.

Binary post-asymptotic giant branch (post-AGB) stars are products of a poorly understood binary interaction process that occurs during the AGB phase. These systems comprise of a post-AGB primary star, a main-sequence secondary companion and a stable circumbinary disc. Studying the structure and properties of these circumbinary discs is crucial for gaining insight into the binary interaction process that governs post-AGB binaries as well as comprehending the disc's creation, evolution, and its interaction with the post-AGB binary system. We aim to use near-infrared polarimetric imaging to investigate the morphology and potential substructures of circumbinary discs around eight representative post-AGB binary stars. To achieve this, we performed polarimetric differential imaging in H and Y bands using the high-angular resolution capabilities of the European Southern Observatory-Very Large Telescope/SPHERE-Infra-Red Dual-beam Imaging and Spectroscopy instrument. We resolved the extended circumbinary disc structure for a diverse sample of eight post-AGB binary systems. Our analysis provided the first estimates of the disc scale-height for two of the systems: IW Car and IRAS 15469-5311. We also investigated the morphological differences between the full discs (with the inner rim at the dust sublimation radius) and transition discs (which are expected to have larger inner cavities), as well as similarities to protoplanetary disks around young stellar objects. We found that the transition discs displayed a more intricate and asymmetric configuration. Surprisingly, no correlation was found between the over-resolved flux in near-IR interferometric data and the polarimetric observations, suggesting that scattering of light on the disc surface may not be the primary cause of the observed over-resolved flux component.

We present a comprehensive analysis of X-ray pulsar 4U 1626-67 during its current spin-down (2SD) state, following a recent torque reversal. Since its discovery, this ultra-compact binary has experienced multiple torque states, transitioning from spin-up (1SU) during 1977-1990 to spin-down (1SD) during 1990-2008, and again spin-up (2SU) until 2023. From NuSTAR observation of May 2023, we have investigated the timing and spectral properties of this pulsar during its 2SD phase, while also comparing them with previous spin-up-down states. For energies upto 8 keV, a distinct bi-horned pulse profile was observed during the spin-up phase, while several sub-structures emerged during spin-down. Beyond 8 keV, a broad asymmetric peak was consistently observed across all torque states. The pulse fraction during the 2SD phase was higher than that during 2SU phase. A prominent ~46.8 mHz quasi-periodic oscillation has been exclusively detected during the spin-down phase. The broadband spectrum during the 2SD phase is described by empirical NPEX model, cyclotron absorption feature and its first harmonic. The spectrum during 2SU phase requires an additional blackbody component and asymmetry in the cyclotron absorption line. A significant flux drop by a factor of ~3 in the 2SD was observed.

We use 23 years of astrometric and radial velocity data on the orbit of the star S0-2 to constrain a hypothetical intermediate-mass black hole orbiting the massive black hole Sgr A* at the Galactic center. The data place upper limits on variations of the orientation of the stellar orbit (inclination, nodal angle, and pericenter) at levels between 0.02 and 0.07 degrees per year. We use a combination of analytic estimates and full numerical integrations of the orbit of S0-2 in the presence of a black-hole binary. For a companion IMBH whose semi-major axis $a_c$ is larger than that of S0-2 (1020 a.u.), we find that in the region between 1000 and 4000 a.u., a companion black hole with mass $m_c$ between $10^3$ and $10^5 M_\odot$ is excluded, with a boundary behaving as $a_c \sim m_c^{1/3}$. For a companion with $a_c < 1020$ a.u., we find that a black hole with mass between $10^3$ and $10^5 \, M_\odot$ is again excluded, with a boundary behaving as $a_c \sim m_c^{-1/2}$. These bounds arise from quadrupolar perturbations of the orbit of S0-2. However, significantly stronger bounds on the mass of an inner companion arise from the fact that the location of S0-2 is measured relative to the bright emission of Sgr A*. As a consequence, that separation is perturbed by the ``wobble'' of Sgr A* about the center of mass between it and the companion, leading to ``apparent'' perturbations of S0-2's orbit that also include a dipole component. The result is a set of bounds as small as $400 \, M_\odot$ at 200 a.u.; the numerical simulations suggest a bound from these effects varying as $a_c \sim m_c^{-1}$. We compare and contrast our results with those from a recent analysis by the GRAVITY collaboration.

Recently, three optical tidal disruption event (TDE) candidates discovered by the Zwicky Transient Facility (ZTF) have been suggested to be coincident with high-energy neutrinos. They all exhibit unusually strong dust infrared (IR) echoes, with their peak times matching the neutrino arrival time even better than the optical peaks. We hereby report on two new TDE candidates that are spatially and temporally coincident with neutrinos by matching our sample of mid-infrared outbursts in nearby galaxies (MIRONG) with Gold alerts of IceCube high-energy neutrino events up to June 2022. The two candidates show negligible optical variability according to their ZTF light curves and can therefore be classified as part of the growing population of obscured TDE candidates. The chance probability of finding two such candidates about $\sim3\%$ by redistributing the MIRONG sources randomly in the SDSS footprint, which will be as low as $\sim0.1\%$ (or $\sim0.2\%$) if we limit to sources with increased fluxes (or variability amplitudes) comparable with the matched two sources. Our findings further support the potential connection between high-energy neutrinos and TDEs in dusty environments by increasing the total number of neutrino-associated TDE and TDE candidates to five, although the underlying physics remains poorly understood.

This research is to determine at Earth the high-energy neutrino flux coming from the galactic core, and from the many other accretion disks within the galactic core. It is estimated there are 10,000 such accretion disk within the cubic parsec of the galactic core alone and many more in the galactic core halo. There are various neutrino detectors, such as IceCube, which can detect energetic neutrinos. However, the direct galactic core neutrino flux is exceptionally low, so very few neutrinos from the galactic core are measured. We created two models to simulate the galactic core neutrino flux. To better estimate the neutrino flux we randomly distributed the accretion disks and generated bodies of varying sizes. This was then used to determine the ultra-high energy neutrino flux. Since it is extremely difficult to determine neutrino direction from interactions of neutrinos, we envision an application where the energetic galactic core neutrinos are gravitationally focused by the Sun with a large light collecting power of eleven to twelve orders of magnitude. They could interact in a planet's atmosphere where the produced showers containing energetic charged particles can produce Cherenkov rings imageable by an orbiting spacecraft or upward going muons which can be observed in a cosmic ray experiment. Estimating the flux will provide a general approximation of the number of ultra-high energy neutrinos that should reach Earth, with Earth being the overall detector. Moreover, studying the neutrinos will provide more information on the conditions of the galactic region, and allow characterizations of it to be formed.

The recent Gaia third data release (DR3) has brought some new exciting data about stellar binaries. It provides new opportunities to fully characterize more stellar systems and contribute to enforce our global knowledge of stars behaviour. By combining the new Gaia non-single stars catalog with double-lined spectroscopic binaries (SB2), one can determine the individual masses and luminosities of the components. To fit an empirical mass-luminosity relation in the Gaia G band, lower mass stars need to be added. Those can be derived using Gaia resolved wide binaries combined with literature data. Using the BINARYS tool, we combine the astrometric non-single star solutions in the Gaia DR3 with SB2 data from two other catalogs : the 9th Catalogue of Spectroscopic Binary orbits (SB9) and APOGEE. We also look for low mass stars resolved in Gaia with direct imaging and Hipparcos data or literature mass fraction. The combination of Gaia astrometric non-single star solutions with double-lined spectroscopic data enabled to characterize 43 binary systems with SB9 and 13 with APOGEE. We further derive the masses of 6 low mass binaries resolved with Gaia. We then derive an empirical mass-luminosity relation in the Gaia G band down to 0.12 Msun.

The understanding of planet formation has changed recently, embracing the new idea of pebble accretion. This means that the influx of pebbles from the outer regions of planet-forming disks to their inner zones could determine the composition of planets and their atmospheres. The solid and molecular components delivered to the planet-forming region can be best characterized by mid-infrared spectroscopy. With Spitzer low-resolution (R=100, 600) spectroscopy, this approach was limited to the detection of abundant molecules such as H2O, C2H2, HCN and CO2. This contribution will present first results of the MINDS (MIRI mid-IR Disk Survey, PI: Th. Henning) project. Due do the sensitivity and spectral resolution (R~1500-3500) provided by JWST we now have a unique tool to obtain the full inventory of chemistry in the inner disks of solar-types stars and brown dwarfs, including also less abundant hydrocarbons and isotopologues. The Integral Field Unit (IFU) capabilities enable at the same time spatial studies of the continuum and line emission in extended sources such as debris disks, the flying saucer and also the search for mid-IR signatures of forming planets in systems such as PDS70. These JWST observations are complementary to ALMA and NOEMA observations of the outer disk chemistry; together these datasets provide an integral view of the processes occurring during the planet formation phase.

Galaxy cluster mass functions are a function of cosmology, but mass is not a direct observable, and systematic errors abound in all its observable proxies. Mass-free inference can bypass this challenge, but it requires large suites of simulations spanning a range of cosmologies and models for directly observable quantities. In this work, we devise a U-net - an image-to-image machine learning algorithm - to ``paint'' the IllustrisTNG model of baryons onto dark-matter-only simulations of galaxy clusters. Using 761 galaxy clusters with $M_{200c} \gtrsim 10^{14}M_\odot$ from the TNG-300 simulation at $z<1$, we train the algorithm to read in maps of projected dark matter mass and output maps of projected gas density, temperature, and X-ray flux. The models train in under an hour on two GPUs, and then predict baryonic images for $\sim2700$ dark matter maps drawn from the TNG-300 dark-matter-only (DMO) simulation in under two minutes. Despite being trained on individual images, the model reproduces the true scaling relation and scatter for the $M_{DM}-L_X$, as well as the distribution functions of the cluster X-ray luminosity and gas mass. For just one decade in cluster mass, the model reproduces three orders of magnitude in $L_X$. The model is biased slightly high when using dark matter maps from the DMO simulation, which is known to have ultra-dense clumps that in the FP simulation are smoothed out by baryonic feedback. The model performs well on inputs from TNG-300-2, whose mass resolution is 8 times coarser; further degrading the resolution biases the predicted luminosity function high. We conclude that U-net-based baryon painting is a promising technique to build large simulated cluster catalogs which can be used to improve cluster cosmology by combining existing full-physics and large $N$-body simulations.

We present high-resolution transmission spectroscopy of WASP-76b with GRACES/Gemini North obtained as part of the ExoGemS survey. With a broad spectral range of 400-1050 nm and a relatively high resolution of ~66,000, these observations are particularly well-suited to searching for atomic and molecular atmospheric species via the Doppler cross-correlation technique. We recover absorption features due to neutral iron (Fe I), sodium (Na I), and ionized calcium (Ca II) at high significance (>5$\sigma$), and investigate possible atmospheric temperatures and wind speeds. We also report tentative (>3$\sigma$) detections of Li I, K I, Cr I, and V I in the atmosphere of WASP-76b. Finally, we report non-detections of a number of other species, some of which have previously been detected with other instruments. Through model injection/recovery tests, we demonstrate that many of these species are not expected to be detected in our observations. These results allow us to place GRACES and the ExoGemS survey in context with other high-resolution optical spectrographs.

Stars form from the gravitational collapse of turbulent, magnetized molecular cloud cores. Our non-ideal MHD simulations reveal that the intrinsically anisotropic magnetic resistance to gravity during the core collapse naturally generates dense gravo-magneto-sheetlets within inner protostellar envelopes -- disrupted versions of classical sheet-like pseudodisks. They are embedded in a magnetically dominant background, where less dense materials flow along the local magnetic field lines and accumulate in the dense sheetlets. The sheetlets, which feed the disk predominantly through its upper and lower surfaces, are the primary channels for mass and angular momentum transfer from the envelope to the disk. The protostellar disk inherits a small fraction (up to 10\%) of the magnetic flux from the envelope, resulting in a disk-averaged net vertical field strength of 1-10 mG and a somewhat stronger toroidal field, potentially detectable through ALMA Zeeman observations. The inherited magnetic field from the envelope plays a dominant role in disk angular momentum evolution, enabling the formation of gravitationally stable disks in cases where the disk field is relatively well-coupled to the gas. Its influence remains significant even in marginally gravitationally unstable disks formed in the more magnetically diffusive cases, removing angular momentum at a rate comparable to or greater than that caused by spiral arms. The magnetically driven disk evolution is consistent with the apparent scarcity of prominent spirals capable of driving rapid accretion in deeply embedded protostellar disks. The dense gravo-magneto-sheetlets observed in our simulations may correspond to the ``accretion streamers" increasingly detected around protostars.

With the successful impact of the NASA DART spacecraft in the Didymos-Dimorphos binary asteroid system, we provide an initial analysis of the post-impact perturbed binary asteroid dynamics. To compare our simulation results with observations, we introduce a set of "observable elements" calculated using only the physical separation of the binary asteroid, rather than traditional Keplerian elements. Using numerical methods that treat the fully spin-orbit-coupled dynamics, we estimate the system's mass and the impact-induced changes in orbital velocity, semimajor axis, and eccentricity. We find that the changes to the mutual orbit depend strongly on the separation distance between Didymos and Dimorphos at the time of impact. If Dimorphos enters a tumbling state after the impact, this may be observable through changes in the system's eccentricity and orbit period. We also find that any DART-induced reshaping of Dimorphos would generally reduce the required change in orbital velocity to achieve the measured post-impact orbit period and will be assessed by the ESA Hera mission in 2027.

Thermal dust polarization is a powerful tool to probe magnetic fields ($\textbf{B}$), grain magnetic properties, and grain sizes. However, a systematic study of the dependence of synthetic dust polarization on grain properties in protostellar environments is not yet available. In this paper, we post-process a non-ideal MHD simulation of a collapsing protostellar core with our updated POLARIS to study in detail the effects of grain magnetic properties and grain growth on dust polarization. We found that superparamagnetic (SPM) grains can produce high polarization degree $p \sim 10-40\%$ beyond $\sim 500$ au because of their efficient magnetic alignment by magnetically enhanced Radiative Torque (MRAT) mechanism. The magnetic field tangling due to turbulence in the envelope causes the decrease of $p$ with emission intensity $I$ as $p\propto I^{\alpha}$ with the slope $\alpha \sim -0.3$. But within 500 au, SPM grains tend to have weak internal alignment and be aligned with $\textbf{B}$ by RAdiative Torque mechanism only, producing lower $p \sim 1\%$ and larger $\alpha \sim -0.6$. For paramagnetic (PM) grains, their inefficient magnetic alignment produces $p << 1\%$, and the depolarization happens with a steep slope of $\alpha \sim -0.9$ owing to the alignment loss of large grains toward the protostar. Grain growth can help to increase $p$ and weaken the depolarization effect caused by turbulence beyond $500$ au for SPM grains. But for SPM grains within $\sim 500$ au and for PM grains, increasing $a_{\rm max}$ enhances the depolarization effect due to the increasing amount of large grains with inefficient alignment. Finally, we found that the polarization angle dispersion function $S$ increases with increasing iron inclusions and $a_{\rm max}$. Our findings reveal the dependence of magnetic field strength measured using the Davis-Chandrashekhar-Fermi technique on grain alignment degree.

In active galactic nuclei, the relationship between UV and X-ray luminosity is well studied (often characterised by $\alpha_\text{ox}$) but often with heterogeneous samples. We have parametrized the intrinsic distribution of X-ray luminosity, $L_\text{X}$, for the optically-selected sample of SDSS quasars in the Stripe 82 and XXL fields across redshifts 0.5-3.5. We make use of the available XMM observations and a custom pipeline to produce Bayesian sensitivity curves that are used to derive the intrinsic X-ray distribution in a hierarchical Bayesian framework. We find that the X-ray luminosity distribution is well described by a Gaussian function in ${\log_{10}}L_\text{X}$ space with a mean that is dependent on the monochromatic 2500A UV luminosity, $L_{2500}$. We also observe some redshift dependence of the distribution. The mean of the $L_\text{X}$ distribution increases with redshift while the width decreases. This weak but significant redshift dependence leads to $L_{2500}$-$L_\text{X}$ and $L_{2500}$-$\alpha_\text{ox}$ relations that evolve with redshift, and we produce a redshift- and $L_{2500}$-dependent $\alpha_\text{ox}$ equation. The increasing average black hole mass with redshift in our sample points to black hole mass as a potential driver of the redshift evolution.

Star ADS 9173=WDS 14135+5147=Hip 69483 is a complex system. The B component has a spectroscopic companion, whose orbit with a period of 4.9 years has been known since 1986. The Gaia telescope has detected a distant faint pair over 100 arcsec away from the bright AB pair. In our article, we study the movement in a bright pair based on long-term observations with the 26-inch refractor of the Pulkovo Observatory. The AB pair orbit with a period of 6306 years was calculated using the apparent motion parameters (AMP) method. The astrometric orbit of the component B was determined on the basis of the residuals of the homogeneous CCD observations up to 2023 with the 26-inch refractor. It is in agreement with the spectroscopic one. The remaining secondary residuals show a wave with a period of approximately 20 years, the reasons for which are discussed.

The Rapid Imaging Planetary Spectrograph (RIPS) was designed as a long-slit high-resolution spectrograph for the specific application of studying atmospheres of spatially extended solar system bodies. With heritage in terrestrial airglow instruments, RIPS uses an echelle grating and order-sorting filters to obtain optical spectra at resolving powers up to R~127,000. An ultra-narrowband image from the reflective slit jaws is captured concurrently with each spectrum on the same EMCCD detector. The "rapid" portion of RIPS' moniker stems from its ability to capture high frame rate data streams, which enables the established technique known as "lucky imaging" to be extended to spatially resolved spectroscopy. Resonantly scattered emission lines of alkali metals, in particular, are sufficiently bright to be measured in short integration times. RIPS has mapped the distributions of Na and K emissions in Mercury's tenuous exosphere, which exhibit dynamic behavior coupled to the planet's plasma and meteoroid environment. An important application is daylight observations of Mercury at solar telescopes since synoptic context on the exosphere's distribution comprises valuable ground-based support for the upcoming BepiColombo orbital mission. As a conventional long slit spectrograph, RIPS has targeted the Moon's surface-bound exosphere where structure in linewidth and brightness as a function of tangent altitude are observed. At the Galilean moons, RIPS can study the plasma interaction with Io and place new constraints on the sputtered atmosphere of Europa, which in turn provides insight into the salinity of Europa's subsurface ocean. The instrumental design and construction are described herein, and these astronomical observations are presented to illustrate RIPS' performance as a visiting instrument at three different telescope facilities.

We study the perturbations up to the 2nd-order for a power-law inflation driven by a scalar field in synchronous coordinates. We solve the 1st-order perturbed Einstein equation and scalar field equation, give the 1st-order solutions for all the scalar, vector, and tensor metric perturbations, as well as the perturbed scalar field. During inflation, the 1st-order tensor perturbation is a wave and is decoupled from other perturbations, the scalar metric perturbation and the perturbed scalar field are coupled waves, propagating at the speed of light, differing from those in the dust and relativistic fluid models. The 1st-order vector perturbation is not wave and just decreases during inflation. The 2nd-order perturbed Einstein equation is similar in structure to the 1st-order one, but various products of the 1st-order perturbations occur as the effective source, among which the scalar-scalar coupling is considered in this paper. The solutions of all the 2nd-order perturbations consist of a homogeneous part similar to the 1st-order solutions, and an inhomogeneous part in a form of integrations of the effective source. The 2nd-order vector perturbation is also a wave since the effective source is composed of the 1st-order waves. We perform the residual gauge transformations between synchronous coordinates up to the 2nd-order, and identify the 1st-order and 2nd-order gauge modes. The 1st-order tensor perturbation is gauge independent, and all others are gauge dependent. We examine four 1st-order gauge invariant scalar perturbations and use the zero-point energy of the scalar field to determine the four primordial spectra.

We investigate realistic models of compact objects, focusing on neutron and strange stars, composed by dense matter and dark energy in the form of a simple fluid or scalar field interacting with matter. For the dark energy component, we use equations of state compatible with cosmological observations. This requirement strongly constrains possible deviations from the simple $\Lambda$-Cold Dark-Matter model with EoS $p_{d}=-\rho_{d}$ at least for small densities of the dark component. But we can propose that the density of dark energy interacting with matter can reach large values in relativistic stars and affects the star parameters such as the mass and radius. Simple models of dark energy are considered. Then we investigated possible effects from modified gravity choosing to study the $R^2$ model combined with dark energy. Finally, the case of dark energy as scalar field non-minimally interacting with gravity is considered.

Next generation direct dark matter (DM) detection experiments will be have unprecedented capabilities to explore coherent neutrino-nucleus scattering (CE$\nu$NS) complementary to dedicated neutrino experiments. We demonstrate that future DM experiments can effectively probe nonstandard neutrino interactions (NSI) mediated by scalar fields in the scattering of solar and atmospheric neutrinos. We set first limits on $S_1$ leptoquark models that result in sizable $\mu-d$ and $\tau-d$ sector neutrino NSI CE$\nu$NS contributions using LUX-ZEPLIN (LZ) data. As we show, near future DM experiments reaching $\sim \mathcal{O}(100)$ton-year exposure, such as argon-based ARGO and xenon-based DARWIN, can probe parameter space of leptoquarks beyond the reach of current and planned collider facilities. We also analyze for the first time prospects for testing NSI in lead-based detectors. We discuss the ability of leptoquarks in the parameter space of interest to also explain the neutrino masses and $(g-2)_\mu$ observations.

The LUX-ZEPLIN (LZ) experiment is a dark matter detector centered on a dual-phase xenon time projection chamber. We report searches for new physics appearing through few-keV-scale electron recoils, using the experiment's first exposure of 60 live days and a fiducial mass of 5.5t. The data are found to be consistent with a background-only hypothesis, and limits are set on models for new physics including solar axion electron coupling, solar neutrino magnetic moment and millicharge, and electron couplings to galactic axion-like particles and hidden photons. Similar limits are set on weakly interacting massive particle (WIMP) dark matter producing signals through ionized atomic states from the Migdal effect.

Galactic cosmic rays (GCRs) are affected by solar modulation while they propagate through the heliosphere. The study of the time variation of GCR spectra observed at Earth can shed light on the underlying physical processes, specifically diffusion and particle drifts. We combine a state-of-the art 3D numerical model of GCR transport in the heliosphere with a neural-network-accelerated Markov chain Monte Carlo to constrain the rigidity and time dependence of the global transport coefficients, using precise GCR data from the PAMELA and AMS-02 experiments between 2006 and 2019.

The IceCube South Pole Neutrino Observatory is a Cherenkov detector instrumented in a cubic kilometer of ice at the South Pole. IceCube's primary scientific goal is the detection of TeV neutrino emissions from astrophysical sources. At the lower center of the IceCube array, there is a subdetector called DeepCore, which has a denser configuration that makes it possible to lower the energy threshold of IceCube and observe GeV-scale neutrinos, opening the window to atmospheric neutrino oscillations studies. Advances in physics sensitivity have recently been achieved by employing Convolutional Neural Networks to reconstruct neutrino interactions in the DeepCore detector. In this contribution, the recent IceCube result from the atmospheric muon neutrino disappearance analysis using the CNN-reconstructed neutrino sample is presented and compared to the existing worldwide measurements.

The Advanced LIGO and Advanced Virgo detectors have enabled the confident detection of dozens of mergers of black holes and neutron stars. However, the presence of detector noise transients (glitches) hinders the search for these gravitational wave (GW) signals. We prototyped a restructuring of Gravity Spy's classification model to distinguish between glitches and astrophysical signals. Our method is able to correctly classify three-quarters of retracted candidate events in O3b as non-astrophysical and 100\% of the confirmed astrophysical events as true signals. This approach will inform candidate event validation efforts in the latest observing run.

This paper contributes to the growing body of research on deep learning methods for solar flare prediction, primarily focusing on highly overlooked near-limb flares and utilizing the attribution methods to provide a post hoc qualitative explanation of the model's predictions. We present a solar flare prediction model, which is trained using hourly full-disk line-of-sight magnetogram images and employs a binary prediction mode to forecast $\geq$M-class flares that may occur within the following 24-hour period. To address the class imbalance, we employ a fusion of data augmentation and class weighting techniques; and evaluate the overall performance of our model using the true skill statistic (TSS) and Heidke skill score (HSS). Moreover, we applied three attribution methods, namely Guided Gradient-weighted Class Activation Mapping, Integrated Gradients, and Deep Shapley Additive Explanations, to interpret and cross-validate our model's predictions with the explanations. Our analysis revealed that full-disk prediction of solar flares aligns with characteristics related to active regions (ARs). In particular, the key findings of this study are: (1) our deep learning models achieved an average TSS=0.51 and HSS=0.35, and the results further demonstrate a competent capability to predict near-limb solar flares and (2) the qualitative analysis of the model explanation indicates that our model identifies and uses features associated with ARs in central and near-limb locations from full-disk magnetograms to make corresponding predictions. In other words, our models learn the shape and texture-based characteristics of flaring ARs even at near-limb areas, which is a novel and critical capability with significant implications for operational forecasting.

We study extreme mass-ratio binary systems in which a stellar mass compact object spirals into a supermassive black hole surrounded by a scalar cloud. Scalar clouds can form through superradiant instabilities of massive scalar fields around spinning black holes and can also serve as a proxy for dark matter halos. Our framework is fully relativistic and assumes that the impact of the cloud on the geometry can be treated perturbatively. As a proof of concept, here we consider a point particle in circular, equatorial motion around a non-spinning black hole surrounded either by a spherically symmetric or a dipolar non-axisymmetric scalar cloud, but the framework can in principle be generalized to generic black hole spins and scalar cloud geometries. We compute the leading-order power lost by the point particle due to scalar radiation and show that, in some regimes, it can dominate over gravitational-wave emission. We confirm the presence of striking signatures due to the presence of a scalar cloud that had been predicted using Newtonian approximations, such as resonances that can give rise to sinking and floating orbits, as well as "sharp features" in the power lost by the particle at given orbital radii. Finally, for a spherically symmetric scalar cloud, we also compute the leading-order corrections to the black-hole geometry and to the gravitational-wave energy flux, focusing only on axial metric perturbations for the latter. We find that, for non-compact clouds, the corrections to the (axial) gravitational-wave fluxes at high frequencies can be understood in terms of a gravitational-redshift effect, in agreement with previous works.

We show that it is possible for fermion condensation of the Nambu-Jona-Lasinio type to induce a non-singular bounce that smoothly connects a phase of slow contraction to a phase of expansion. A chiral condensate -- a non-zero vacuum expectation value of the spinor bilinear $\langle \bar{\Psi}\Psi \rangle$ -- can form spontaneously after a slow contraction phase smooths and flattens the universe and the Ricci-curvature exceeds a critical value. In this approach, a high density of spin-aligned free fermions is not required, which avoids the problem of generating a large anisotropy and initiating chaotic mixmaster behavior during the bounce phase.

This paper deals with astrophysical accretion onto the magnetically charged Euler-Heisenberg black holes with scalar hair. We examine the accretion process of a variety of perfect fluids, including polytropic and isothermal fluids of the ultra-stiff, ultra-relativistic, and sub-relativistic forms, when fluid is accreting in the vicinity of the black hole. By using the Hamiltonian dynamical approach, we can find the sonic or critical points numerically for the various types of fluids that are accreting onto the black hole. Furthermore, for several types of fluids, the solution is provided in closed form, expressing phase diagram curves. We compute the mass accretion rate of a magnetically charged Euler-Heisenberg black hole with scalar hair. We observe that the maximum accretion rate is attained for small values of the black hole parameters. We may be able to understand the physical mechanism of accretion onto black holes using the outcomes of this investigation.

The main aim of this study is to reveal curved space and particle physics effects on the formation of Bose-Einstein condensate (BEC) scalar fields in cosmology and around a black hole. Cosmological scalar fields for dark energy and dark matter may be considered as a result of Bose-Einstein condensation. In this regard, our main attention will be devoted to BECs in curved space. By considering the dynamics of a BEC scalar field at a microscopic level, we first study the initial phase of the formation of condensation in cosmology. To this end, we initially introduce an effective Minkowski space formulation that enables considering only the effect of particle physics processes, excluding the effect of gravitational particle production and enabling us to see cosmological evolution more easily. Then, by using this formulation, we study a model with a trilinear coupling that induces the processes. After considering the phase evolution of the produced particles, we find that they evolve towards the formation of a BEC if some specific conditions are satisfied. In principle, the effective Minkowski space formulation introduced in this study can be applied to particle physics processes in any sufficiently smooth spacetime. In this regard, we also analyse if a BEC scalar field is realized in the spacetime around a Reissner - Nordstr\"om black hole. We find that the produced particles of particle physics processes are localized in a region around the black hole and have a tendency toward condensation if the emerged particles are much heavier than ingoing particles. We also find that such a configuration is phenomenologically viable only if the scalars and the black hole have dark electric charges. Finally, we consider gravitational collapse around Schwarzschild black holes and form a first step towards a study in future about the effects of gravitational collapse on Bose-Einstein condensation.

Magnetorotational instability (MRI) is the most likely mechanism for angular momentum transport in accretion disks. However, despite numerous efforts, the quest for an unambiguous identification of MRI in laboratory is still ongoing. For this purpose, a big Taylor-Couette experiment with liquid sodium is currently under construction within DRESDYN project. In preparation for these experiments, we have determined optimal parameters for axisymmetric standard MRI (SMRI) with an axial magnetic field and analyzed its nonlinear evolution. In this sequel paper, we investigated the linear and nonlinear dynamics of non-axisymmetric SMRI modes in a similar magnetized Taylor-Couette setup at very small magnetic Prandtl numbers $Pm \sim 10^{-5}$ of liquid sodium. We showed that the achievable magnetic Reynolds $Rm\sim 40$ and Lundquist $Lu\sim10$ numbers in this experiment are large enough for the growth of non-axisymmetric $|m|=1$ SMRI modes, but the corresponding critical $Rm_c$ is about 2-3 times higher than that of axisymmetric SMRI. Then we followed the nonlinear evolution of these modes and analyzed the structure of the saturated state and its scaling properties with $Re$. For $Re \lesssim 10^4$, the non-axisymmetric SMRI modes do not saturate and eventually decay due to the modification of the radial profile of the mean azimuthal velocity by the nonlinear saturation of the axisymmetric SMRI. By contrast, for large $Re \gtrsim 10^4$, a rapid growth and saturation of other, nonmagnetic type of non-axisymmetric modes occur, which are radially localized in the turbulent boundary layer near the inner cylinder wall. The saturation amplitude of these non-axisymmetric modes is always a few orders smaller than that of the axisymmetric one. We showed that the scaling relations for the magnetic energy and torque derived for axisymmetric SMRI in our previous study well carry over to the present 3D case.

Whether the first law of black hole mechanics is correct is an important question in black holes physics. Subjected to current limited gravitational wave events, we propose its weaker version that permits a relatively large perturbation to a black hole system and implement a simple test with the first event GW150914. Confronting the strain data with the theory, we obtain the constraint on the deviation parameter $\alpha=0.07\pm0.11$, which indicates that this weaker version is valid at the 68\% confidence level. This result implies that the first law of black hole mechanics may be correct.

We review recent trends on inflationary dynamics in the context of viable modified gravity theories. After providing a general overview of the inflationary paradigm emphasizing on what problems of hot Big Bang theory inflation solves, and a somewhat introductory presentation of single field inflationary theories with minimal and non-minimal couplings, we review how inflation can be realized in terms of several string motivated models of inflation, which involve Gauss-Bonnet couplings of the scalar field, higher order derivatives of the scalar field, and some subclasses of viable Horndeski theories. We also present and analyze inflation in the context of Chern-Simons theories of gravity, including various subcases and generalizations of string corrected modified gravities which also contain Chern-Simons correction terms, with the scalar field being identified with the invisible axion, which is the most viable to date dark matter candidate. We also provide a detailed account of vacuum $f(R)$ gravity inflation, and also inflation in $f(R,\phi)$ and kinetic-corrected $f(R,\phi)$ theories of gravity. In the end of the review we discuss the technique for calculating the overall effect of modified gravity on the waveform of the standard general relativistic gravitational wave form.

The soon-to-be-realized, global network of neutrino telescopes will allow new opportunities for collaboration between detectors. While each detector is distinct, they share the same underlying physical processes and detection principles. The full simulation chain for these telescopes is typically proprietary which limits the opportunity for joint studies. This means there is no consistent framework for simulating multiple detectors. To overcome these challenges, we introduce Prometheus, an open-source simulation tool for neutrino telescopes. Prometheus simulates neutrino injection and final state and photon propagation in both ice and water. It also supports user-supplied injection and detector specifications. In this contribution, we will introduce the software; show its runtime performance; and highlight successes in reproducing simulation results from multiple ice- and water-based observatories.

The stellar configurations of quark stars are studied using perturbative QCD (pQCD) for the equation of state (EoS). The neutron star structures with equation of state obtained from Brussels-Montreal extended Skyrme interaction are also explored. The deconfinement phase transition from quark to hadron phase in stellar interior for matter under extreme pressure is accomplished by employing the Maxwell construction. The influence of hybrid EoS on the jump in energy density has been investigated. To study hybrid stars the BSk24 hadronic model and pQCD EoS for the quark phase have been used. The properties of hybrid stars in the view of the very recent astrophysical observations have been examined. We find that the gravitational mass might exceed 2.3 $M_\odot$ in some cases, comparable with the observed mass of the pulsar PSR J0952-0607 recently detected.

We examine the anisotropy originated from a first-order vacuum phase transitions through three-dimensional numerical simulations. We apply Bianchi Type-I metric to our model that has one scalar field minimally coupled to the gravity. We calculate the time evolution of the energy density for the shear scalar and the directional Hubble parameters as well as the power spectra for the scalar field and the gravitational radiation although there are a number of caveats for the tensor perturbations in Bianchi Type-I universe. We run simulations with different mass scales of the scalar field, therefore, in addition to investigation of anisotropy via the shear scalar, we also determine at which mass scale the phase transition completes successfully, hence, neglecting the expansion of the Universe does not significantly affect the results. Finally, we showed that such an event may contribute to the total anisotropy depending on the mass scale of the scalar field and the initial population of nucleated bubbles.

We propose a minimal gauged U(1)$_X$ extension of the MSSM with R-parity conservation. In this model, U(1)$_X$ is a generalization of the well-known U(1) $B-L$. Apart from the MSSM particle content, the model includes three right-handed neutrino (RHN) chiral superfields, each carrying a unit U(1)$_X$ charge. In the presence of RHNs, the model is free from all gauge and mixed gauge-gravitational anomalies. However, there are no U(1)$_X$ Higgs chiral superfields with U(1)$_X$ charge $\pm2$ involved in the model. Two of the RHN superfields are assigned an odd R-parity, while the last one ($\Psi$) has an even parity. The U(1)$_X$ symmetry is radiatively broken by the VEV of the scalar component of $\Psi$. As a consequence of the absence of U(1)$_X$ Higgs fields and the novel R-parity assignment, the three light neutrinos consist of one massless neutrino and two Dirac neutrinos. In the early universe, the right-handed components of the Dirac neutrinos are in thermal equilibrium with the SM particles through the U(1)$_X$ gauge ($Z^\prime$) boson. The extra energy density from the RHNs is constrained to avoid disrupting the success of BBN, leading to a lower bound on the scale of U(1)$_X$ symmetry breaking. In our model, a mixture of the U(1)$_X$ gaugino and the fermionic component of $\Psi$ becomes a new dark matter (DM) candidate if it is the lightest sparticle mass eigenstate. We examine this DM phenomenology and identify a parameter region that reproduces the observed DM relic density. Furthermore, we consider constraints from the search for $Z'$ boson resonance at the LHC. The three constraints obtained from the success of BBN, the observed DM relic density, and the $Z^\prime$ resonance search at the LHC complement each other, narrowing down the allowed parameter region.

Reduced Order Quadrature (ROQ) methods can greatly reduce the computational cost of Gravitational Wave (GW) likelihood evaluations, and therefore greatly speed up parameter estimation analyses, which is a vital part to maximize the science output of advanced GW detectors. In this paper, we do an in-depth study of ROQ techniques applied to GW data analysis and present novel algorithms to enhance different aspects of the ROQ bases construction. We improve upon previous ROQ construction algorithms allowing for more efficient bases in regions of parameter space that were previously challenging. In particular, we use singular value decomposition (SVD) methods to characterize the waveform space and choose a reduced order basis close to optimal and also propose improved methods for empirical interpolation node selection, greatly reducing the error added by the empirical interpolation model. To demonstrate the effectiveness of our algorithms, we construct multiple ROQ bases ranging in duration from 4s to 256s for compact binary coalescence (CBC) waveforms including precession and higher order modes. We validate these bases by performing likelihood error tests and P-P tests and explore the speed up they induce both theoretically and empirically with positive results. Furthermore, we conduct end-to-end parameter estimation analyses on several confirmed GW events, showing the validity of our approach in real GW data.

TianQin and LISA are space-based laser interferometer gravitational wave (GW) detectors planned to be launched in the mid-2030s. Both detectors will detect low-frequency GWs around $10^{-2}\,{\rm Hz}$, however, TianQin is more sensitive to frequencies above this common sweet-spot while LISA is more sensitive to frequencies below $10^{-2}\,{\rm Hz}$. Therefore, TianQin and LISA will be able to detect the same sources but with different accuracy for different sources and their parameters. We study the detection distance and the detection accuracy with which TianQin and LISA will be able to detect some of the most important astrophysical sources -- massive black hole binaries, stellar-mass black holes binaries, double white dwarfs, extreme mass ratio inspirals, light and heavy intermediate mass ratio inspirals, as well as the stochastic gravitational background produced by binaries. We, further, study the detection distance and detection accuracy from joint detection. We compare the results obtained by the three detection scenarios highlighting the gains from joint detection as well as the contribution of TianQin and LISA to a combined study of astrophysical sources. In particular, we consider the different orientations, lifetimes, and duty cycles of the two detectors to explore how they can give a more complete picture when working together.

We evaluate several neural-network architectures, both convolutional and recurrent, for gravitational-wave time-series feature extraction by performing point parameter estimation on noisy waveforms from binary-black-hole mergers. We build datasets of 100,000 elements for each of four different waveform models (or approximants) in order to test how approximant choice affects feature extraction. Our choices include \texttt{SEOBNRv4P} and \texttt{IMRPhenomPv3}, which contain only the dominant quadrupole emission mode, alongside \texttt{IMRPhenomPv3HM} and \texttt{NRHybSur3dq8}, which also account for high-order modes. Each dataset element is injected into detector noise corresponding to the third observing run of the LIGO-Virgo-KAGRA (LVK) collaboration. We identify the Temporal Convolutional Network (TCN) architecture as the overall best performer in terms of training and validation losses and absence of overfitting to data. Comparison of results between datasets shows that the choice of waveform approximant for the creation of a dataset conditions the feature extraction ability of a trained network. Hence, care should be taken when building a dataset for the training of neural networks, as certain approximants may result in better network convergence of evaluation metrics. However, this performance does not necessarily translate to data which is more faithful to numerical relativity simulations. We also apply this network on actual signals from LVK runs, finding that its feature-extracting performance can be effective on real data.

Recent images from the Event Horizon Telescope of accreting supermassive black holes (BHs), along with upcoming observations with better sensitivity and angular resolution, offer exciting opportunities to deepen our understanding of spacetime in strong gravitational fields. A significant focus for future BH imaging observations is the direct detection of the "photon ring," a narrow band on the observer's sky that collects extremely lensed photons. The photon ring consists of self-similarly nested subrings which, in spherically-symmetric spacetimes, are neatly indexed by the maximum number of half-loops executed around the BH by the photons that arrive in them. Each subring represents an entire "higher-order" image of the horizon-scale accretion flow. Furthermore, this self-similarity is controlled by a single critical lensing exponent linked to the radial (in)stability of photon orbits near the critical (circular) photon orbit, solely determined by the spacetime geometry. However, extracting such information about the spacetime geometry can be challenging because the observed photon ring is also influenced by the structure of the emitting region. To address this, we conducted a comprehensive study by varying (a) a wide range of emission-zone morphology models and (b) families of spacetime metrics. We find that the lensing exponent can be reliably determined from future observations. This exponent can provide access to the $rr-$component of the spacetime metric, as well as significantly narrow down currently accessible BH parameter spaces. Additionally, the width of the first-order photon subring serves as yet another important discriminator of the spacetime geometry. Finally, observations of flaring events across different wavelengths might reveal time-delayed secondary images, with the delay time providing a promising new way to independently estimate the BH shadow size.