Papers for Monday, Apr 24 2023

We study the $N$-point function of the density contrast to quadratic order in the squeezed limit during the matter-dominated (MD) and radiation-dominated (RD) eras in synchronous gauge. Since synchronous gauge follows the free-fall frame of observers, the equivalence principle dictates that in the gradient approximation for the long-wavelength mode there is only a single, manifestly time-independent consistency relation for the $N$-point function. This simple form is dictated by the initial mapping between synchronous and local coordinates, unlike Newtonian gauge and its correspondingly separate dilation and Newtonian consistency relations. Dynamical effects only appear at quadratic order in the squeezed limit and are again characterized by a change in the local background, also known as the separate universe approach. We show that for the 3-point function the compatibility between these squeezed-limit relations and second-order perturbation theory requires both the initial and dynamical contributions to match, as they do in single-field inflation. This clarifies the role of evolution or late-time projection effects in establishing the consistency relation for observable bispectra, which is especially important for radiation acoustic oscillations and for establishing consistency below the matter-radiation equality scale in the MD era. Defining an appropriate angle and time average of these oscillations is also important for making separate universe predictions of spatially varying local observables during the RD era, which can be useful for a wider range of cosmological predictions beyond $N$-point functions.

Models of active galactic nuclei (AGN) suggest that their circumnuclear media are complex with clumps and filaments, while recent observations hint towards polar extended structures of gas and dust, as opposed to the classical torus paradigm. The X-ray band forms an interesting window to study these circumnuclear media in great detail. In this work, we extended the radiative transfer code SKIRT with the X-ray processes that govern the broadband X-ray spectra of obscured AGN, to study the structure of AGN circumnuclear media in full 3D, based on their reflected X-ray emission. We extended the SKIRT code with Compton scattering on free electrons, photo-absorption and fluorescence by cold atomic gas, scattering on bound electrons, and extinction by dust. This includes a novel treatment of extreme-forward scattering by dust, and a detailed description of anomalous Rayleigh scattering. To verify our X-ray implementation, we performed the first dedicated benchmark of X-ray torus models, comparing five X-ray radiative transfer codes. Finally, we illustrated the 3D nature of the code by producing synthetic X-ray images and spectra of clumpy torus models. SKIRT forms a powerful new tool to model AGN circumnuclear media in full 3D from X-ray to millimetre wavelengths, and is now publicly available. In the X-ray regime, we find an excellent agreement with the simulation results of the MYTorus and RefleX codes, which validates our X-ray implementation. We find some discrepancies with other codes, which motivates the need for a robust framework that can handle non-linear 3D radiative transfer effects. The new X-ray functionalities of the SKIRT code allow for uncomplicated access to a broad suite of 3D X-ray models for AGN that can easily be tested and modified. This will be particularly useful with the advent of X-ray microcalorimeter observations.

Beyond standard summary statistics are necessary to summarize the rich information on non-linear scales in the era of precision galaxy clustering measurements. For the first time, we introduce the 2D k-th nearest neighbor (kNN) statistics as a summary statistic for discrete galaxy fields. This is a direct generalization of the standard 1D kNN by disentangling the projected galaxy distribution from the redshift-space distortion signature along the line-of-sight. We further introduce two different flavors of 2D $k$NNs that trace different aspects of the galaxy field: the standard flavor which tabulates the distances between galaxies and random query points, and a ''DD'' flavor that tabulates the distances between galaxies and galaxies. We showcase the 2D kNNs' strong constraining power both through theoretical arguments and by testing on realistic galaxy mocks. Theoretically, we show that 2D kNNs are computationally efficient and directly generate other statistics such as the popular 2-point correlation function, voids probability function, and counts-in-cell statistics. In a more practical test, we apply the 2D kNN statistics to simulated galaxy mocks that fold in a large range of observational realism and recover parameters of the underlying extended halo occupation distribution (HOD) model that includes velocity bias and galaxy assembly bias. We find unbiased and significantly tighter constraints on all aspects of the HOD model with the 2D kNNs, both compared to the standard 1D kNN, and the classical redshift-space 2-point correlation functions.

Relativistic magnetic reconnection is a non-ideal plasma process that is a source of non-thermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a non-standard finite difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than Yee, resulting in a ~40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss's law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4-3 times over the standard Yee step still yields accurate results, but only translates to modest performance improvements over CKC due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD.

Active galactic nuclei (AGNs) have been proposed as plausible sites hosting a sizable fraction of the binary black hole (BBH) mergers measured through gravitational waves (GWs) by the LIGO-Virgo-Kagra (LVK) experiment. These GWs could be accompanied by radiation feedback due to the interaction of the BBH merger remnant with the AGN disc. We present a new predicted radiation signature driven by the passage of a kicked BBH remnant throughout a thin AGN disc. We analyse the situation of a merger occurring outside the thin disc, where the merger is of second or higher generation in a merging hierarchical sequence. The coalescence produces a kicked BH remnant that eventually plunges into the disc, accretes material, and inflates jet cocoons. We consider the case of a jet cocoon propagating quasi-parallel to the disc plane and study the outflow that results when the cocoon emerges from the disc. Here we focus on the long time-scale emission produced after the disc outflow expands and becomes optically thin. The bolometric luminosity of such disc outflow evolves as $L\propto t^{-7/2}$. Depending on the parameter configuration, the flare produced by the disc outflow could be comparable to or exceed the AGN background emission at near-infrared, optical, and extreme ultraviolet wavelengths appearing $\sim$[20-500] days after the GW event and lasting for $\sim$[1-200] days, accordingly.

We present a new magnetohydrodynamic-particle-in-cell (MHD-PIC) code integrated into the Athena++ framework. It treats energetic particles as in conventional PIC codes while the rest of thermal plasmas are treated as background fluid described by MHD, thus primarily targeting at multi-scale astrophysical problems involving the kinetic physics of the cosmic-rays (CRs). The code is optimized toward efficient vectorization in interpolation and particle deposits, with excellent parallel scaling. The code is also compatible with static/adaptive mesh refinement, with dynamic load balancing to further enhance multi-scale simulations. In addition, we have implemented a compressing/expanding box framework which allows adiabatic driving of CR pressure anisotropy, as well as the $\delta f$ method that can dramatically reduce Poisson noise in problems where distribution function $f$ is only expected to slightly deviate from the background. The code performance is demonstrated over a series of benchmark test problems including particle acceleration in non-relativistic parallel shocks. In particular, we reproduce the linear growth of the CR gyro-resonant (streaming and pressure anisotropy) instabilities, under both the periodic and expanding/compressing box setting. We anticipate the code to open up the avenue for a wide range of astrophysical and plasma physics applications.

Like most galaxies, the Milky Way harbors a supermassive black hole (SMBH) at its center, surrounded by a nuclear star cluster. In this dense star cluster, direct collisions can occur between stars before they evolve off the main-sequence. Using a statistical approach, we characterize the outcomes of these stellar collisions within the inner parsec of the Galactic Center (GC). Close to the SMBH, where the velocity dispersion is larger than the escape speed from a Sun-like star, collisions lead to mass loss. We find that the stellar population within $0.01$ pc is halved within about a Gyr because of destructive collisions. Additionally, we predict a diffuse population of peculiar low-mass stars in the GC. These stars have been divested of their outer layers in the inner $0.01$ pc before migrating to larger distances from the SMBH. Between $0.01$ and $0.1$ pc from the SMBH, collisions can result in mergers. Our results suggest that repeated collisions between lower mass stars can produce massive ($\gtrsim 10$ M$_\odot$) stars, and there may be $\sim 100$ of them residing in this region. We provide predictions on the number of G objects, dust and gas enshrouded stellar objects, that may result from main-sequence stellar collisions. Lastly, we comment on uncertainties in our model and possible connections between stellar collisions and the missing red giants in the GC.

In recent years, there has been ample evidence for the existence of multiple progenitor pathways that can result in Type Ia supernova (SNe Ia), including SNe Ia of sub-Chandrasekhar mass origin best distinguished by their redder colors and higher Si II velocities near peak brightness. These SNe can contaminate the population of normal events used for cosmological analyses, creating unwanted biases in the final analyses. Given that many current and future surveys using SNe Ia as cosmological probes will not have the resources to take a spectrum of all the events, likely only getting host redshifts long after the SNe Ia have faded, we need to turn to methods that could separate these populations based purely on photometry or host properties. Here, we present a study of a sample of well observed, nearby SNe Ia and their hosts to determine if there are significant enough difference between these populations that can be discerned only from the stellar population properties of their hosts. Our results indicate that the global host properties, including star formation, stellar mass, stellar population age, and dust attenuation, of sub-Chandrasekhar mass explosions do not differ significantly from those of normal mass origin. However, we do find evidence using Na I D equivalent widths that the local environments of sub-Chandrasekhar mass explosions are more dust-affected than normal SNe Ia. Future work requires strengthening photometric probes of sub-Chandrasekhar SNe and their local environments to distinguish these events.

We present a generalised estimator for the autocorrelation function, S-ACF, which is an extended version of the standard estimator of the autocorrelation function (ACF). S-ACF is a versatile definition that can robustly and efficiently extract periodicity and signal shape information from a time series, independent of the time sampling and with minimal assumptions about the underlying process. Calculating the autocorrelation of irregularly sampled time series becomes possible by generalising the lag of the standard estimator of the ACF to a real parameter and introducing the notion of selection and weight functions. We show that the S-ACF reduces to the standard ACF estimator for regularly sampled time series. Using a large number of synthetic time series we demonstrate that the performance of the S-ACF is as good or better than commonly used Gaussian and rectangular kernel estimators, and is comparable to a combination of interpolation and the standard estimator. We apply the S-ACF to astrophysical data by extracting rotation periods for the spotted star KIC 5110407, and compare our results to Gaussian process (GP) regression and Lomb-Scargle (LS) periodograms. We find that the S-ACF periods typically agree better with those from GP regression than from LS periodograms, especially in cases where there is evolution in the signal shape. The S-ACF has a wide range of potential applications and should be useful in quantitative science disciplines where irregularly sampled time series occur. A Python implementation of the S-ACF is available under the MIT license.

We present a catalog of fundamental stellar properties for 7,673 evolved stars, including stellar radii and masses, determined from the combination of spectroscopic observations from the Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey IV (SDSS), and asteroseismology from K2. The resulting APO-K2 catalog provides spectroscopically derived temperatures and metallicities, asteroseismic global parameters, evolutionary states, and asteroseismically-derived masses and radii. Additionally, we include kinematic information from \textit{Gaia}. We investigate the multi-dimensional space of abundance, stellar mass, and velocity with an eye toward applications in Galactic archaeology. The APO-K2 sample has a large population of low metallicity stars ($\sim$288 at [M/H] $\leq$ $-$1), and their asteroseismic masses are larger than astrophysical estimates. We argue that this may reflect offsets in the adopted fundamental temperature scale for metal-poor stars rather than metallicity-dependent issues with interpreting asteroseismic data. We characterize the kinematic properties of the population as a function of $\alpha$-enhancement and position in the disk and identify those stars in the sample that are candidate components of the \textit{Gaia-Enceladus} merger. Importantly, we characterize the selection function for the APO-K2 sample as a function of metallicity, radius, mass, $\nu_{\mathrm{max}}$, color, and magnitude referencing Galactic simulations and target selection criteria to enable robust statistical inferences with the catalog.

This article presents an overview of the US National Gemini Office (US NGO) and its role within the International Gemini Observatory user community. Throughout the years, the US NGO charter changed considerably to accommodate the evolving needs of astronomers and the observatory. The current landscape of observational astronomy requires effective communication between stakeholders and reliable/accessible data reduction tools and products, which minimize the time between data gathering and publication of scientific results. Because of that, the US NGO heavily invests in producing data reduction tutorials and cookbooks. Recently, the US NGO started engaging with the Gemini user community through social media, and the results have been encouraging, increasing the observatory's visibility. The US NGO staff developed tools to assess whether the support provided to the user community is sufficient and effective, through website analytics and social media engagement numbers. These quantitative metrics serve as the baseline for internal reporting and directing efforts to new or current products. In the era of the NSF's National Optical-Infrared Astronomy Research Laboratory (NOIRLab), the US NGO is well-positioned to be the liaison between the US user base and the Gemini Observatory. Furthermore, collaborations within NOIRLab programs, such as the Astro Data Lab and the Time Allocation Committee, enhance the US NGO outreach to attract users and develop new products. The future landscape laid out by the Astro 2020 report confirms the need to establish such synergies and provide more integrated user support services to the astronomical community at large.

The Davis-Chandrasekhar-Fermi (DCF) method is widely employed to estimate the mean magnetic field strength in astrophysical plasmas. In this study, we present a numerical investigation using the DCF method in conjunction with a promising new diagnostic tool for studying magnetic fields: the polarization of spectral lines resulting from the atomic alignment effect. We obtain synthetic spectro-polarimetry observations from 3D magnetohydrodynamic (MHD) turbulence simulations and estimate the mean magnetic field projected onto the plane of the sky using the DCF method with GSA polarization maps and a modification to account for the driving scale of turbulence. We also compare the method to the classical DCF approach using dust polarization observations. Our observations indicate that the modified DCF method correctly estimates the plane-of-sky projected magnetic field strengths for sub-Alfv\'enic turbulence using a newly proposed correction factor of $\xi' \in 0.35 - 0.75$. We find that the field strengths are accurately obtained for all magnetic field inclination and azimuth angles. We also observe a minimum threshold for the mean magnetic field inclination angle with respect to the line of sight, $\theta_B \sim 16^\circ$, for the method. The magnetic field dispersion traced by the polarization from the spectral lines is comparable in accuracy to dust polarization, while mitigating some of the uncertainties associated with dust observations. The measurements of the DCF observables from the same atomic/ionic line targets ensure the same origin for the magnetic field and velocity fluctuations and offer a possibility of tracing the 3D direction of the magnetic field.

Type Ib and Ic supernovae (SNe Ib/Ic) originate from hydrogen-deficient massive star progenitors, of which the exact properties are still much debated. Using the SN data in the literature, we investigate the optical $B-V$ color of SNe Ib/Ic at the $V-$band peak and show that SNe Ib are systematically bluer than SNe Ic. We construct SN models from helium-rich and helium-poor progenitors of various masses using the radiation hydrodynamics code STELLA and discuss how the $B-V$ color at the $V-$band peak is affected by $^{56}$Ni to ejecta mass ratios, $^{56}$Ni mixing and presence/absence of the helium envelope. We argue that the dichotomy in the amounts of helium in the progenitors plays the primary role in making the observed systematic color difference at the optical peak, in favor of the most commonly invoked SN scenario that SNe Ib and SNe Ic progenitors are helium-rich and helium-poor, respectively.

Within the next decade, atmospheric O$_2$ on Earth-like M dwarf planets may be accessible with visible--near-infrared, high spectral resolution extremely large ground-based telescope (ELT) instruments. However, the prospects for using ELTs to detect environmental properties that provide context for O$_2$ have not been thoroughly explored. Additional molecules may help indicate planetary habitability, rule out abiotically generated O$_2$, or reveal alternative biosignatures. To understand the accessibility of environmental context using ELT spectra, we simulate high-resolution transit transmission spectra of previously-generated evolved terrestrial atmospheres. We consider inhabited pre-industrial and Archean Earth-like atmospheres, and lifeless worlds with abiotic O$_2$ buildup from CO$_2$ and H$_2$O photolysis. All atmospheres are self-consistent with M2V--M8V dwarf host stars. Our simulations include explicit treatment of systematic and telluric effects to model high-resolution spectra for GMT, TMT, and E-ELT configurations for systems 5 and 12 pc from Earth. Using the cross-correlation technique, we determine the detectability of major species in these atmospheres: O$_2$, O$_3$, CH$_4$, CO$_2$, CO, H$_2$O, and C$_2$H$_6$. Our results suggest that CH$_4$ and CO$_2$ are the most accessible molecules for terrestrial planets transiting a range of M dwarf hosts using an E-ELT, TMT, or GMT sized telescope, and that the O$_2$ NIR and H$_2$O 0.9 $\mu$m bands may also be accessible with more observation time. Although this technique still faces considerable challenges, the ELTs will provide access to the atmospheres of terrestrial planets transiting earlier-type M-dwarf hosts that may not be possible using JWST.

In the collisionless plasmas of radiatively inefficient accretion flows, heating and acceleration of ions and electrons is not well understood. Recent studies in the gyrokinetic limit revealed the importance of incorporating both the compressive and Alfvenic cascades when calculating the partition of dissipated energy between the plasma species. In this paper, we use a covariant analytic model of the accretion flow to explore the impact of compressive and Alfvenic heating, Coulomb collisions, compressional heating, and radiative cooling on the radial temperature profiles of ions and electrons. We show that, independent of the partition of heat between the plasma species, even a small fraction of turbulent energy dissipated to the electrons makes their temperature scale with a virial profile and the ion-to-electron temperature ratio smaller than in the case of pure Coulomb heating. In contrast, the presence of compressive cascades makes this ratio larger because compressive turbulent energy is channeled primarily into the ions. We calculate the ion-to-electron temperature in the inner accretion flow for a broad range of plasma properties, mass accretion rates, and black hole spins and show that it ranges between $5 \lesssim T_i/T_e \lesssim 40$. We provide a physically motivated expression for this ratio that can be used to calculate observables from simulations of black hole accretion flows for a wide range of conditions.

Ultra high-energy (UHE) cosmic rays (CRs) from distant sources interact with intergalactic radiation fields, leading to their spallation and attenuation. They are also deflected in intergalactic magnetic fields (IGMFs), particularly those associated with Mpc-scale structures. These deflections extend the propagation times of CR particles, forming a magnetic horizon for each CR species. The cumulative cooling and interactions of a CR ensemble also modifies their spectral shape and composition observed on Earth. We construct a transport formulation to calculate the observed UHE CR spectral composition for 4 classes of source population. The effects on CR propagation brought about by IGMFs are modeled as scattering processes during transport, by centers associated with cosmic filaments. Our calculations demonstrate that IGMFs can have a marked effect on observed UHE CRs, and that source population models are degenerate with IGMF properties. Interpretation of observations, including the endorsement or rejection of any particular source classes, thus needs careful consideration of the structural properties and evolution of IGMFs. Future observations providing tighter constraints on IGMF properties will significantly improve confidence in assessing UHE CR sources and their intrinsic CR production properties.

In spite of much work, the overall spiral structure morphology of the Milky Way remains somewhat uncertain. In the last two decades, accurate distance measurements have provided us with an opportunity to solve this issue. Using the precise locations of very young objects, for the first time, we propose that our galaxy has a multiple-arm morphology that consists of two-arm symmetry (the Perseus and Norma Arms) in the inner parts and that extends to the outer parts, where there are several long, irregular arms (the Centaurus, Sagittarius, Carina, Outer, and Local Arms).

Long-lived massive magnetars are expected to be remnants of some binary neutron star (BNS) mergers. In this paper, we argue that the magnetic powered flaring activities of these merged magnetars would occur dominantly in their early millisecond-period-spin phase, which is in the timescale of days. Such flares endure significant absorption by the ejecta from the BNS collision, and their detectable energy range is from 0.1-10 MeV, in a time-lag of $\sim$ days after the merger events indicated by the gravitational wave chirps. We estimate the rate of such flares in different energy ranges, and find that there could have been ~0.1-10 cases detected by Fermi/GBM. A careful search for $\sim10$ milliseconds spin period modulation in weak short gamma-ray bursts (GRBs) may identify them from the archival data. Future MeV detectors can detect them at a rate from a few to tens per year. The recent report on the Quasi-Period-Oscillation found in two BASTE GRBs should not be considered as cases of such flares, for they were detection in a lower energy range and with a much shorter period spin modulation.

We use the calibrations by Calamida et al. and by Hilker et al., and the standardised synthetic photometry in the v, b, and y Stromgren passbands from Gaia DR3 BP/RP spectra, to obtain photometric metallicities for a selected sample of 694233 old Galactic giant stars having |b|>20.0 and parallax uncertainties lower than 10%. The zero point of both sets of photometric metallicities has been shifted to to ensure optimal match with the spectroscopic [Fe/H] values for 44785 stars in common with APOGEE DR17, focusing on the metallicity range where they provide the highest accuracy. The metallicities derived in this way from the Calamida et al. calibration display a typical accuracy of ~0.1 dex and 1 sigma precision ~0.2 dex in the range -2.2 <=[Fe/H]<= -0.4, while they show a systematic trend with [Fe/H] at higher metallicity, beyond the applicability range of the relation. Those derived from the Hilker et al. calibration display, in general, worse precision, and lower accuracy in the metal-poor regime, but have a median accuracy < 0.05 dex for [Fe/H]>= -0.8. These results are confirmed and, consequently, the metallicities validated, by comparison with large sets of spectroscopic metallicities from various surveys. The newly obtained metallicities are used to derive metallicity distributions for several previously identified sub-structures in the Galactic halo with an unprecedented number of stars. The catalogue including both sets of metallicities and the associated uncertainties is made publicly available.

Meteoroids of a low-inclination stream hit the Earth arriving from a direction near the ecliptic. The radiant area of stream like this is often divided into two parts: one is situated northward and the other southward of the ecliptic. In other words, two showers are caused by such a stream. Well-known examples of such showers are the Northern Taurids, #17, and Southern Taurids, #2, or the Northern $\delta$-Aquariids, #26, and Southern $\delta$-Aquariids, #5. While the meteoroids of the northern shower collide with the Earth in the descending node, those of the southern shower collide with our planet in the ascending node of their orbits. Because of this circumstance and tradition, the northern and southern showers must be distinguished. Unfortunately, this is not always the case with meteor showers listed in the IAU Meteor Data Center (MDC). For the same shower, some authors reported a set of its mean parameters corresponding to the northern shower and other authors to the southern shower. We found eleven such cases in the MDC. In this paper, we propose corrections of these mis-identifications.

The relative high fluxes of the Galactic ultra-luminous X-ray pulsar Swift J0243 allow a detailed study of its spin-down regime in quiescence state, for a first time. After the 2017 giant outburst, its spin frequencies show a linear decreasing trend with some variations due to minor outbursts. The linear spin-down rate is $\sim-1.9\times10^{-12}$ Hz/s during the period of lowest luminosity, from which one can infer a dipole field $\sim1.75\times10^{13}$ G. The $\dot{\nu}-L$ relation during the spin-down regime is complex, and the $\dot{\nu}$ is close to 0 when the luminosity reaches both the high end ($L_{38}\sim0.3$) and the lowest value ($L_{38}\sim0.03$). The luminosity of zero-torque is different for the giant outburst and other minor outbursts. It is likely due to different accretion flows for different types of outburst, as evidenced by the differences of the spectra and pulse profiles at a similar luminosity for different types of outburst (giant or not). The pulse profile changes from double peaks in the spin-up state to a single broad peak in the low spin-down regime, indicating the emission beam/region is larger in the low spin-down regime. These results show that accretion is still ongoing in the low spin-down regime for which the neutron star is supposed to be in a propeller state.

Ejection mechanism of transient relativistic jets from black hole binaries is studied. Based on the observations of the limit-cycle behaviors of the superluminal jet source, GRS 1915+105, we infer that the transient jet ejections could happen just after the slim disk emerging from the standard disk at some distance runs over the standard disk and reaches the vicinity of the central black hole. The standard disk releases about a half of the gravitational energy in the course of the accretion, but the released radiative energy could be absorbed by the optically thick slim disk covering the standard disk in this situation. Then, since the inward motion of the slim disk is much faster than that of the standard disk, a quantity of energy released by an amount of gas in the standard disk is received by the much smaller amount of gas in the slim disk. As the result, the energy per mass received by the slim disk is expected to be largely amplified and is estimated to get highly relativistic. Since the energy is much larger than the gravitational energy, the height of the slim disk could significantly increase. Hence, the innermost part of the slim disk from which almost all the angular momentum has been transferred outward could have a much larger height than the black hole size and collide with one another around the central axis of the disk, turning to an outward flow along the axis normal to the disk plane. The flow in this direction can be approximated to be that through the de Laval nozzle and could become supersonic near the distance where the flow has the smallest cross section.

The initial and boundary conditions of the Galactic star formation in molecular clouds are not well understood. In an effort to shed new light on this long-standing problem, we measured properties of dense cores and filamentary structures in the Vela C molecular cloud, observed with Herschel. We applied the getsf extraction method to separate the components of sources and filaments from each other and their backgrounds, before detecting, measuring, and cataloging the structures. The cores and filamentary structures constitute 40% of the total mass of Vela C, most of the material is in the low-density molecular background cloud. We selected 570 reliable cores, of which 149 are the protostellar cores and 421 are the starless cores. Almost 78% of the starless cores were identified with the gravitationally bound prestellar cores. The exponent of the CMF (alpha = 1.35) is identical to that of the Salpeter IMF. We selected 68 filaments with at least one side that appeared not blended with adjacent structures. The filament widths are in the range of 0.15 pc to 0.63 pc, and have a median value of W = 0.3(0.11) pc. The surface densities of filaments are well correlated with their contrasts and linear densities. Within uncertainties of the filament instability criterion, many filaments may well be both supercritical and subcritical. A large fraction of filaments may definitely be considered supercritical, in which are found 94 prestellar cores, 83 protostellar cores, and only 1 unbound starless core. Taking into account the uncertainties, the supercritical filaments contain only prestellar and protostellar cores. Our findings support the idea that there exists a direct relationship between the CMF and IMF and that filaments play a key role in the formation of prestellar cores, which is consistent with the previous Herschel results.

Magnetars are believed as neutron stars (NSs) with strong magnetic fields. X-ray flares and fast radio bursts (FRBs) have been observed from the magnetar (soft gamma-ray repeater, SGR J1935+2154). We propose that the phase transition of the NS can power the FRBs and SGRs.Based on the equation of state provided by the MIT bag model and the mean field approximation, we solve the Tolman-Oppenheimer-Volkoff equations to get the NS structure. With spin-down of the NS, the hadronic shell gradually transfers to the quark shell.The gravitational potential energy released by one time of the phase transition can be achieved. The released energy, time interval between two successive phase transitions, and glitch are all consistent with the observations of the FRBs and the X-ray flares from SGR J1935+2154. We conclude that the phase transition of an NS is a plausible mechanism to power the SGRs as well as the repeating FRBs.

To accurately link the radio and optical Celestial Reference Frames (CRFs) at optical bright end, i.e., with Gaia G band magnitude < 13, increasing number and improving sky distribution of radio stars with accurate astrometric parameters from both Very Long Baseline Interferometry (VLBI) and Gaia measurements are mandatory. We selected two radio stars HD 199178 and AR Lacertae as the target for a pilot program for the frame link, using the Very Long Baseline Array (VLBA) at 15 GHz at six epochs spanning about 1 year, to measure their astrometric parameters. The measured parallax of HD 199178 is $8.949 \pm 0.059$ mas and the proper motion is $\mu_\alpha cos \delta = 26.393 \pm 0.093$, $\mu_\delta = -0.950 \pm 0.083~mas~yr^{-1}$, while the parallax of AR Lac is $23.459 \pm 0.094$ mas and the proper motion is $\mu_\alpha cos \delta = -51.906 \pm 0.138$, $\mu_\delta = 46.732 \pm 0.131~mas~yr^{-1}$. Our VLBI measured astrometric parameters have accuracies about 4-5 times better than the corresponding historic VLBI measurements and comparable accuracies with those from Gaia, validating the feasibility of frame link using radio stars. With the updated astrometric parameters for these two stars, there is a 25% reduction of the uncertainties on the Y axis for both orientation and spin parameters.

Stochastic inflation resolves primordial perturbations non-linearly, probing their probability distribution deep into its non-Gaussian tail. The strongest perturbations collapse into primordial black holes. In typical black-hole-producing single-field inflation, the strongest stochastic kicks occur during a period of constant roll. In this paper, I solve the stochastic constant-roll system, drawing the stochastic kicks from a numerically computed power spectrum, beyond the usual de Sitter approximation. The perturbation probability distribution is an analytical function of the integrated power spectrum $\sigma_k^2$ and the second slow-roll parameter $\epsilon_2$. With a large $\epsilon_2$, stochastic effects can reduce the height of the curvature power spectrum required to form asteroid mass black holes from $10^{-2}$ to $10^{-3}$. I compare these results to studies with the non-stochastic $\Delta N$ formalism.

With the coincident detections of electromagnetic radiation together with gravitational waves (GW170817) or neutrinos (TXS 0506+056), the new era of multimessenger astrophysics has begun. Of particular interest are the searches for correlation between the high-energy astrophysical neutrinos detected by the IceCube Observatory and gamma-ray photons detected by the Fermi Large Area Telescope (LAT). So far, only sources detected by the LAT have been considered in correlation with IceCube neutrinos, neglecting any emission from sources too faint to be resolved individually. Here, we present the first cross-correlation analysis considering the unresolved gamma-ray background (UGRB) and IceCube events. We perform a thorough sensitivity study and, given the lack of identified correlation, we place upper limits on the fraction of the observed neutrinos that would be produced in proton-proton (p-p) or proton-gamma (p-gamma) interactions from the population of sources contributing to the UGRB emission and dominating its spatial anisotropy (aka blazars). Our analysis suggests that, under the assumption that there is no intrinsic cutoff and/or hardening of the spectrum above Fermi-LAT energies, and that all gamma-rays from the unresolved blazars dominating the UGRB fluctuation field are produced by neutral pions from p-p (p-gamma) interactions, up to 60% (30%) of such population may contribute to the total neutrino events observed by IceCube. This translates into a O(1%) maximum contribution to the astrophysical high-energy neutrino flux observed by IceCube at 100 TeV.

The circumnuclear material around Active Galactic Nuclei (AGN) is one of the essential components of the obscuration-based unification model. However, our understanding of the circumnuclear material in terms of its geometrical shape, structure and its dependence on accretion rate is still debated. In this paper, we present the multi-epoch broadband X-ray spectral modelling of a nearby Compton-thick AGN in Circinus galaxy. We utilise all the available hard X-ray ($> 10$ keV) observations taken from different telescopes, $i.e.,$ $BeppoSAX$, $Suzaku$, $NuSTAR$ and $AstroSat$, at ten different epochs across $22$ years from $1998$ to $2020$. The $3.0-79$ keV broadband X-ray spectral modelling using physically-motivated models, namely MYTORUS, BORUS02 and UXCLUMPY, infers the presence of a torus with a low covering factor of $0.28$, an inclination angle of $77^{\circ}$ $-$ $81^{\circ}$ and Compton-thick line-of-sight column densities ($N_{\rm H,LOS} = 4.13~-~9.26~\times~10^{24}$ cm$^{-2}$) in all the epochs. The joint multi-epoch spectral modelling suggests that the overall structure of the torus is likely to remain unchanged. However, we find tentative evidence for the variable line-of-sight column density on timescales ranging from one day to one week to a few years, suggesting a clumpy circumnuclear material located at sub-parsec to tens of parsec scales.

We analyse spectroscopic and photometric transits of the hot Jupiters WASP-52b and HAT-P30b obtained with ESPRESSO, Eulercam and NGTS for both targets, and additional TESS data for HAT-P-30. Our goal is to update the system parameters and refine our knowledge of the host star surfaces. For WASP-52, the companion planet has occulted starspots in the past, and as such our aim was to use the reloaded Rossiter-McLaughlin technique to directly probe its starspot properties. Unfortunately, we find no evidence for starspot occultations in the datasets herein. Additionally, we searched for stellar surface differential rotation (DR) and any centre-to-limb variation (CLV) due to convection, but return a null detection of both. This is unsurprising for WASP-52, given its relatively cool temperature, high magnetic activity (which leads to lower CLV), and projected obliquity near 0 degrees (meaning the transit chord is less likely to cross several stellar latitudes). For HAT-P-30, this result was more surprising given its hotter effective temperature, lower magnetic field, and high projected obliquity (near 70 degrees). To explore the reasons behind the null DR and CLV detection for HAT-P-30, we simulated a variety of scenarios. We find that either the CLV present on HAT-P-30 is below the solar level or the presence of DR prevents a CLV detection given the precision of the data herein. A careful treatment of both DR and CLV is required, especially for systems with high impact factors, due to potential degeneracies between the two. Future observations and/or a sophisticated treatment of the red noise present in the data (likely due to granulation) is required to refine the DR and CLV for these particular systems; such observations would also present another opportunity to try to examine starspots on WASP-52.

Precise measurements of the spectra of secondary and primary cosmic rays are crucial for understanding the origin and propagation of those energetic particles. The High Energy cosmic-Radiation Detection (HERD) facility onboard China`s Space Station, which is expected to operate in 2027, will push the direct measurements of cosmic ray fluxes precisely up to PeV energies. In this work, we investigate the potential of HERD on studying the propagation of cosmic rays using the measurements of boron, carbon, and oxygen spectra. We find that, compared with the current results, the new HERD measurements can improve the accuracy of the propagation parameters by 8\% to 40\%. The constraints on the injection spectra at high energies will also be improved.

Numerical simulations have now shown that Population III (Pop III) stars can form in binaries and small clusters and that these stars can be in close proximity to each other. If so, they could be subject to binary interactions such as mass exchange that could profoundly alter their evolution, ionizing UV and Lyman-Werner (LW) photon emission and explosion yields, with important consequences for early cosmological reionization and chemical enrichment. Here we investigate the evolution of Pop III and extremely metal-poor binary stars with the MESA code. We find that interactions ranging from stable mass transfer to common envelope evolution can occur in these binaries for a wide range of mass ratios and initial separations. Mass transfer can nearly double UV photon yields in some of these binaries with respect to their individual stars by extending the life of the companion star, which in turn can enhance early cosmological reionization but also suppress the formation of later generations of primordial stars. Binary interactions can also have large effects on the nucleosynthetic yields of the stars by promoting or removing them into or out of mass ranges for specific SN types. We provide fits to total photon yields for the binaries in our study for use in cosmological simulations.

We present data from the Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) covering the 15-year span from April 2006 through the end of 2020. APOLLO measures the earth-moon separation by recording the round-trip travel time of photons from the Apache Point Observatory to five retro-reflector arrays on the moon. The APOLLO data set, combined with the 50-year archive of measurements from other lunar laser ranging (LLR) stations, can be used to probe fundamental physics such as gravity and Lorentz symmetry, as well as properties of the moon itself. We show that range measurements performed by APOLLO since 2006 have a median nightly accuracy of 1.7 mm, which is significantly better than other LLR stations.

We investigate the late-time neutrino emission powered by fallback mass accretion onto proto-neutron star (PNS), using neutrino radiation-hydrodynamic simulations with full Boltzmann neutrino transport. We follow the time evolution of accretion flow onto PNS until the system reaches a steady state. A standing shock wave is commonly formed in the accretion flow, whereas the shock radius varies depending on mass accretion rate and PNS mass. A sharp increase in temperature emerges in the vicinity of PNS ($\sim 10$ km), which characterizes neutrino emission. Both neutrino luminosity and average energy become higher with increasing mass accretion rate and PNS mass. The mean energy of emitted neutrinos is in the range of $10\lesssim\epsilon\lesssim20\,\mathrm{MeV}$, which is higher than that estimated from PNS cooling models ($\lesssim10\,\mathrm{MeV}$). Assuming a distance to core-collapse supernova of $10\,\mathrm{kpc}$, we quantify neutrino event rates for Super-Kamiokande (Super-K) and DUNE. The estimated detection rates are well above the background, and their energy-dependent features are qualitatively different from those expected from PNS cooling models. Another notable feature is that the neutrino emission is strongly flavor dependent, exhibiting that the neutrino event rate hinges on the neutrino oscillation model. We estimate them in the case with adiabatic Mikheev-Smirnov-Wolfenstein model, and show that the normal- and inverted mass hierarchy offer the large number of neutrino detection in Super-K and DUNE, respectively. Hence the simultaneous observation with Super-K and DUNE of the fallback neutrinos will provide a strong constraint on neutrino mass hierarchy.

A substantial body of phenomenological and theoretical work over the last few years strengthens the possibility that the vacuum energy density (VED) of the universe is dynamical, and in particular that it adopts the `running vacuum model' (RVM) form, in which the VED evolves mildly as $\delta \rho_{\rm vac}(H)\sim \nu_{\rm eff} m_{\rm Pl}^2{\cal O}\left(H^2\right)$, where $H$ is the Hubble rate and $\nu_{\rm eff}$ is a (small) free parameter. This dynamical scenario is grounded on recent studies of quantum field theory (QFT) in curved spacetime and also on string theory. It turns out that what we call the `cosmological constant', $\Lambda$, is no longer a rigid parameter but the nearly sustained value of $8\pi G(H)\rho_{\rm vac}(H)$ around (any) given epoch $H(t)$, where $G(H)$ is the gravitational coupling, which can also be very mildly running (logarithmically). Of particular interest is the possibility suggested in past works that such a running may help to cure the cosmological tensions afflicting the $\Lambda$CDM. In the current study, we reanalyze it in full and we find it becomes further buttressed. Using the modern cosmological data, namely a compilation of the latest $SNIa+BAO+$H(z)$+LSS+CMB$ observations, we probe to which extent the RVM provides a quality fit better than the concordance $\Lambda$CDM model, paying particular emphasis on its impact on the $\sigma_8$ and $H_0$ tensions. We utilize the Einstein-Boltzmann system solver $CLASS$ and the Monte Carlo sampler $MontePython$ for the statistical analysis, as well as the statistical $DIC$ criterion to compare the running vacuum against the rigid vacuum ($\nu_{\rm eff} = 0$). We show that with a tiny amount of vacuum dynamics ($|\nu_{\rm eff}|\ll 1$) the global fit can improve significantly with respect to the $\Lambda$CDM and the mentioned tensions may subside to inconspicuous levels.

Using the power of numerical relativity, we show that, beginning from generic initial conditions that are far from flat, homogeneous and isotropic and have a large Weyl curvature, a period of slow contraction rapidly drives spacetime towards vanishingly small Weyl curvature as the total energy density grows, thus providing a dynamical mechanism that satisfies the Weyl Curvature Hypothesis. We also demonstrate a tight correlation between the Weyl Curvature Hypothesis and ultralocal behavior for canonical scalar fields with a sufficiently steep negative potential energy density.

High-frequency gravitational waves can be detected by observing the frequency modulation they impart on photons. We discuss fundamental limitations to this method related to the fact that it is impossible to construct a perfectly rigid detector. We then propose several novel methods to search for O(MHz-GHz) gravitational waves based on the frequency modulation induced in the spectrum of an intense laser beam, by applying optical frequency demodulation techniques, or by using optical atomic clock technology. We find promising sensitivities across a broad frequency range.

Wavelets are widely used in various disciplines to analyse signals both in space and scale. Whilst many fields measure data on manifolds (i.e., the sphere), often data are only observed on a partial region of the manifold. Wavelets are a typical approach to data of this form, but the wavelet coefficients that overlap with the boundary become contaminated and must be removed for accurate analysis. Another approach is to estimate the region of missing data and to use existing whole-manifold methods for analysis. However, both approaches introduce uncertainty into any analysis. Slepian wavelets enable one to work directly with only the data present, thus avoiding the problems discussed above. Applications of Slepian wavelets to areas of research measuring data on the partial sphere include gravitational/magnetic fields in geodesy, ground-based measurements in astronomy, measurements of whole-planet properties in planetary science, geomagnetism of the Earth, and cosmic microwave background analyses.

We investigate the relation between quantum nonlocality and gravity at the astrophysical scale, both in the classical and quantum regimes. Considering particle pairs orbiting in the strong gravitational field of ultra-compact objects, we find that the violation of Bell inequality acquires an angular modulation factor that strongly depends on the nature of the gravitational source. We show how such gravitationally-induced modulation of quantum nonlocality readily discriminates between black holes (both classical and inclusive of quantum corrections) and string fuzzballs, i.e., the true quantum description of ultra-compact objects according to string theory. These findings promote Bell nonlocality as a potentially key tool in comparing different models of classical and quantum gravity and putting them to the test.

The increasing availability of machines relying on non-GPU architectures, such as ARM A64FX in high-performance computing, provides a set of interesting challenges to application developers. In addition to requiring code portability across different parallelization schemes, programs targeting these architectures have to be highly adaptable in terms of compute kernel sizes to accommodate different execution characteristics for various heterogeneous workloads. In this paper, we demonstrate an approach to code and performance portability that is based entirely on established standards in the industry. In addition to applying Kokkos as an abstraction over the execution of compute kernels on different heterogeneous execution environments, we show that the use of standard C++ constructs as exposed by the HPX runtime system enables superb portability in terms of code and performance based on the real-world Octo-Tiger astrophysics application. We report our experience with porting Octo-Tiger to the ARM A64FX architecture provided by Stony Brook's Ookami and Riken's Supercomputer Fugaku and compare the resulting performance with that achieved on well established GPU-oriented HPC machines such as ORNL's Summit, NERSC's Perlmutter and CSCS's Piz Daint systems. Octo-Tiger scaled well on Supercomputer Fugaku without any major code changes due to the abstraction levels provided by HPX and Kokkos. Adding vectorization support for ARM's SVE to Octo-Tiger was trivial thanks to using standard C++

We describe the design of a radio interferometer composed of a Global Navigation Satellite Systems (GNSS) antenna and a Very Long Baseline Interferometry (VLBI) radio telescope. Our eventual goal is to use this interferometer for geodetic applications including local tie measurements. The GNSS element of the interferometer uses a unique software-defined receiving system and modified commercial geodetic-quality GNSS antenna. We ran three observing sessions in 2022 between a 25 m radio telescope in Fort Davis, Texas (FD-VLBA), a transportable GNSS antenna placed within 100 meters, and a GNSS antenna placed at a distance of about 9 km. We have detected a strong interferometric response with a Signal-to-Noise Ratio (SNR) of over 1000 from Global Positioning System (GPS) and Galileo satellites. We also observed natural radio sources including Galactic supernova remnants and Active Galactic Nuclei (AGN) located as far as one gigaparsec, thus extending the range of sources that can be referenced to a GNSS antenna by 18 orders of magnitude. These detections represent the first observations made with a GNSS antenna to radio telescope interferometer. We have developed a novel technique based on a Precise Point Positioning (PPP) solution of the recorded GNSS signal that allows us to extend integration time at 1.5 GHz to at least 20 minutes without any noticeable SNR degradation when a rubidium frequency standard is used.

Symmetries play an important role in fundamental physics. In gravity and field theories, particular attention has been paid to Weyl (or conformal) symmetry. However, once the theory contains a scalar field, conformal transformations of the metric can be considered a subclass of a more general type of transformation, so-called disformal transformation. Here, we investigate the implications of pure disformal symmetry in the Universe. We derive the form of general disformal invariant tensors from which we build the most general disformal invariant action. We argue that, in cosmology, disformal symmetry amounts to require that the lapse function is fully replaced by a (time-like) scalar field at the level of the action. We then show that disformal symmetry is in general an exactly equivalent formulation of general mimetic gravity. Lastly, we go beyond mimetic gravity and find that a particular class of invariance leads to seemingly Ostrogradski-like (with higher derivatives) Lagrangians, which are nevertheless absent of Ostrogradski ghosts in a cosmological background, despite having an additional degree of freedom. We also propose an application of our formalism to find new invertible disformal transformations, where the coefficient involves higher derivatives and curvature, further expanding the theory space of scalar-tensor theories.