Papers for Thursday, Mar 30 2023

To mimic protoplanetary discs (PPDs) evolution, 2D simulations with self-gravity must introduce a softening prescription of the gravitational potential. When the disc is only made of gas the smoothing length is proportional to the gas scale height. On the other hand when a dust component is included, the question arises as whether the smoothing length approach can still be used to quantify not only the dust self-gravity but also its gravitational interaction with gas. We identified grey areas in the standard smoothing length formalism for computing self-gravity in PPDs uniquely made of gas. We revisit the smoothing length approach which is then generalised to two phases when the dust component can be considered as a pressureless fluid. Analytical developments are used to approximate the vertically averaged self-gravity when the smoothing length is not assumed to be constant but rather a spatial function. We obtained an analytical expression for the space varying smoothing length, which strongly improves the accuracy of the self-gravity computation. For the first time, this method is generalised to address bi-fluid interactions in a PPD: two additional smoothing lengths are proposed for featuring an isolated dusty disc and gas-dust self-gravity interactions. We checked that our method remains compatible with standard fast Fourier transform algorithms and evaluated computational costs. Our space varying smoothing length permits (i) to solve the contradictions inherent to a constant smoothing length hypothesis, (ii) to fit accurately the 3D vertically averaged self-gravity and (iii) is applicable to a bi-fluid description of PPDs with the use of two additional smoothing lengths. Such results are crucial to enable realistic 2D numerical simulations accounting for self-gravity and are important to deepen our understanding of planetesimals formation and type I migration.

We use the TNG50 simulation of the IllustrisTNG project to study cold, dense clouds of gas in the circumgalactic media (CGM) of Milky Way-like galaxies. We find that their CGM is typically filled with of order one hundred (thousand) reasonably (marginally) resolved clouds, possible analogs of high-velocity clouds (HVCs). There is a large variation in cloud abundance from galaxy to galaxy, and the physical properties of clouds that we explore -- mass, size, metallicity, pressure, and kinematics -- are also diverse. We quantify the distributions of cloud properties and cloud-background contrasts, providing cosmological inputs for idealized simulations. Clouds characteristically have sub-solar metallicities, diverse shapes, small overdensities ($\chi = n_{\rm cold} / n_{\rm hot} \lesssim 10$), are mostly inflowing, and have sub-virial rotation. At TNG50 resolution, resolved clouds have median masses of $\sim 10^6\,\rm{M_\odot}$ and sizes of $\sim 10$ kpc. Larger clouds are well converged numerically, while the abundance of the smallest clouds increases with resolution, as expected. In TNG50 MW-like haloes, clouds are slightly (severely) under-pressurised relative to their surroundings with respect to total (thermal) pressure, implying that magnetic fields may be important. Clouds are not distributed uniformly throughout the CGM, but are clustered around other clouds, often near baryon-rich satellite galaxies. This suggests that at least some clouds originate from satellites, via direct ram-pressure stripping or otherwise. Finally, we compare with observations of intermediate and high velocity clouds from the real Milky Way halo. TNG50 shows a similar cloud velocity distribution as observations, and predicts a significant population of currently difficult-to-detect low velocity clouds.

We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest-resolution run of the IllustrisTNG suite of $\Lambda$CDM magneto-hydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogs, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogs, i.e. galaxies with stellar disky morphology, with a stellar mass in the range of $M_* = 10^{10.5 - 11.2}$ Msun, and within a MW-like Mpc-scale environment at z=0. These are resolved with baryonic (dark matter) mass resolution of $8.5\times10^4$ Msun ($4.5\times10^5$ Msun) and $\sim150$ pc of average spatial resolution in the star-forming regions: we therefore expand by many factors (2 orders of magnitude) the sample size of cosmologically-simulated analogs with similar ($\times 10$ better) numerical resolution. The majority of TNG50 MW/M31 analogs at $z=0$ exhibit a bar, 60 per cent are star-forming, the sample includes 3 Local Group (LG)-like systems, and a number of galaxies host one or more satellites as massive as e.g. the Magellanic Clouds. Even within such a relatively narrow selection, TNG50 reveals a great diversity in galaxy and halo properties, as well as in past histories. Within the TNG50 sample, it is possible to identify several simulated galaxies whose integral and structural properties are consistent, one or more at a time, with those measured for the Galaxy and Andromeda. With this paper, we document and release a series of broadly applicable data products that build upon the IllustrisTNG public release and aim to facilitate easy access and analysis by public users. These include datacubes across snapshots ($0 \le z \le 7$) for each TNG50 MW/M31-like galaxy, and a series of value-added catalogs that will be continually expanded to provide a convenient and up to date community resource.

Normal mode seismology is a promising means of measuring rotation in gas giant interiors, and ring seismology presents a singular opportunity to do so at Saturn. We calculate Saturn's normal modes of oscillation and zonal gravity field, using nonperturbative methods for normal modes in the rigidly rotating approximation, and perturbative methods for the shifts that Saturn's deep winds induce in the mode frequencies and zonal gravity harmonics. The latter are calculated by solving the thermo-gravitational wind equation in an oblate geometry. Comparing many such models to gravity data and the frequencies of ring patterns excited by Saturn normal modes, we use statistical methods to estimate that Saturn's cloud-level winds extend inward along cylinders before decaying at a depth 0.125-0.138 times Saturn's equatorial radius, or 7,530-8,320 km, consistent with analyses of Cassini gravity and magnetic field data. The seismology is especially useful for pinning down Saturn's poorly constrained deep rotation period, which we estimate at 634.7 min (median) with a 5/95% quantile range 633.8-635.5 min. Outstanding residuals in mode frequencies at low angular degree suggest a more complicated deep interior than has been considered to date. Smaller but still significant residuals at high angular degree also show that our picture for the thermal, composition, and/or rotation profile in Saturn's envelope is not yet complete.

The ubiquitous interstellar turbulence regulates star formation and the scaling relations between the initial velocity differences and the initial separations of stars. We propose that the formation of wide binaries with initial separations $r$ in the range $\sim 10^3~\text{AU} \lesssim r \lesssim 10^5$ AU is a natural consequence of star formation in the turbulent interstellar medium. With the decrease of $r$, the mean turbulent relative velocity $v_\text{tur}$ between a pair of stars decreases, while the largest velocity $v_\text{bon}$ at which they still may be gravitationally bound increases. When $v_\text{tur} < v_\text{bon}$, a wide binary can form. In this formation scenario, we derive the eccentricity distribution $p(e)$ of wide binaries for an arbitrary relative velocity distribution. By adopting a turbulent velocity distribution, we find that wide binaries at a given initial separation generally exhibit a superthermal $p(e)$. This provides a natural explanation for the observed superthermal $p(e)$ of the wide binaries in the Solar neighborhood.

A major event in cosmic history is the genesis of the first starlight in our Universe, ending the ''Dark Ages''. During this epoch, the earliest luminous sources were enshrouded in neutral and pristine gas, which was gradually ionised in a process called ''reionisation''. Hence, one of the brightest emission lines in star-forming galaxies, Lyman-$\alpha$ (Ly-$\alpha$), was predicted to emerge only towards the end of the epoch of reionisation, about one billion years after the Big Bang. However, this picture has been challenged over the past decade by the surprising detection of Ly-$\alpha$ in galaxies less than 500 million years old. Here we show, by taking advantage of both high-resolution and high-sensitivity images from the James Webb Space Telescope programs PRIMER, CEERS and FRESCO, that all galaxies in our sample of Ly-$\alpha$ emitters deep in the epoch of reionisation have close companions. To understand the physical processes that lead to the observed Ly-$\alpha$ emission in our sample, we take advantage of novel on-the-fly radiative transfer magnetohydrodynamical simulations with cosmic ray feedback. We find that in the early Universe, the rapid build up of mass through frequent galactic mergers leads to very bursty star formation which in turn drives episodes of high intrinsic Ly-$\alpha$ emission and facilitates the escape of Ly-$\alpha$ photons along channels cleared of neutral gas. These merging galaxies reside in clustered environments thus creating sufficiently large ionised bubbles. This presents a solution to the long-standing puzzle of the detection of Ly-$\alpha$ emission deep into the epoch of reionisation.

Galaxy clusters are the most massive gravitationally bound structures in the Universe, comprising thousands of galaxies and pervaded by a diffuse, hot ``intracluster medium'' (ICM) that dominates the baryonic content of these systems. The formation and evolution of the ICM across cosmic time is thought to be driven by the continuous accretion of matter from the large-scale filamentary surroundings and dramatic merger events with other clusters or groups. Until now, however, direct observations of the intracluster gas have been limited only to mature clusters in the latter three-quarters of the history of the Universe, and we have been lacking a direct view of the hot, thermalized cluster atmosphere at the epoch when the first massive clusters formed. Here we report the detection (about $6\sigma$) of the thermal Sunyaev-Zeldovich (SZ) effect in the direction of a protocluster. In fact, the SZ signal reveals the ICM thermal energy in a way that is insensitive to cosmological dimming, making it ideal for tracing the thermal history of cosmic structures. This result indicates the presence of a nascent ICM within the Spiderweb protocluster at redshift $z=2.156$, around 10 billion years ago. The amplitude and morphology of the detected signal show that the SZ effect from the protocluster is lower than expected from dynamical considerations and comparable with that of lower-redshift group-scale systems, consistent with expectations for a dynamically active progenitor of a local galaxy cluster.

Spatially-resolved studies of the kinematics of galaxies provide crucial insights into their assembly and evolution, enabling to infer the properties of the dark matter halos, derive the impact of feedback on the ISM, characterize the outflow motions. To date, most of the kinematic studies at z=0.5-3.5 were obtained using emission lines tracing the warm, ionized gas. However, whether these provide an exhaustive or only a partial view of the dynamics of galaxies and of the properties of the ISM is still debated. Complementary insights on the cold gas kinematics are therefore needed. We present ALPAKA, a project aimed at gathering high-resolution observations of CO and [CI] emission lines of star-forming galaxies at z=0.5-3.5 from the ALMA public archive. With 147 hours of total integration time, ALPAKA assembles ~0.25'' observations for 28 star-forming galaxies, the largest sample with spatially-resolved cold gas kinematics as traced by either CO or [CI] at z>0.5. By combining multi-wavelength ancillary data, we derive the stellar masses ($M_{\star}$) and star-formation rates (SFR) for our targets, finding values of $M_{\star}\gtrsim 10^{10}$ M$_{\odot}$ and SFR of 10-3000 M$_{\odot}$/yr. A large fraction of ALPAKA galaxies (19/28) lie in overdense regions (clusters, groups, and protoclusters). We exploit the ALMA data to infer their dynamical state and we find that 19/28 ALPAKA galaxies are rotating disks, 2 are interacting systems, while for the remaining 7 sources the classification is uncertain. The disks have velocity dispersion values that are typically larger in the innermost regions than in the outskirts, with a median value for the entire disk sample of 35$^{+11}_{-9}$ km/s. Despite the bias of our sample towards galaxies hosting very energetic mechanisms, the ALPAKA disks have high ratios of ordered-to-random motion ($V/\sigma$) with a median value of 9$^{+7}_{-2}$.

We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of $M_* = 10^{10.5 - 11.2}~\rm{M_{\odot}}$, and a MW-like Mpc-scale environment at $z=0$. The stellar disks of such TNG50 MW/M31 analogs exhibit a wide diversity of structural properties, including a number of galaxies with disk scalelength and thin and thick disk scaleheights that are comparable to those measured or inferred for the Galaxy and Andromeda. With one set of physical ingredients, TNG50 returns a large variety of flaring flavours and amounts, also for mono-age stellar populations. With this paper, we hence propose a non-parametric characterization of flaring. The typical MW/M31 analogs exhibit disk scaleheights that are $1.5-2$ times larger in the outer than in the inner regions of the disk for both old and young stellar populations, but with a large galaxy-to-galaxy variation. Which stellar population flares more, and by how much, also varies from galaxy to galaxy. TNG50 de facto brackets existing observational constraints for the Galaxy and all previous numerical findings. A link between the amount of flaring and the $z=0$ global galaxy structural properties or merger history is complex. However, a connection between the scaleheights and the local stellar vertical kinematics and gravitational potential is clearly in place.

The dozens of rocky exoplanets discovered in the Circumstellar Habitable Zone (CHZ) currently represent the most suitable places to host life as we know it outside the Solar System. However, the presumed presence of liquid water on the CHZ planets does not guarantee suitable environments for the emergence of life. According to experimental studies, the building blocks of life are most likely produced photochemically in presence of a minimum ultraviolet (UV) flux. On the other hand, high UV flux can be life-threatening, leading to atmospheric erosion and damaging biomolecules essential to life. These arguments raise questions about the actual habitability of CHZ planets around stars other than Solar-type ones, with different UV to bolometric luminosity ratios. By combining the "principle of mediocricy" and recent experimental studies, we define UV boundary conditions (UV-habitable Zone, UHZ) within which life can possibly emerge and evolve. We investigate whether exoplanets discovered in CHZs do indeed experience such conditions. By analysing Swift-UV/Optical Telescope data, we measure the near ultraviolet (NUV) luminosities of 17 stars harbouring 23 planets in their CHZ. We derive an empirical relation between NUV luminosity and stellar effective temperature. We find that eighteen of the CHZ exoplanets actually orbit outside the UHZ, i.e., the NUV luminosity of their M-dwarf hosts is decisively too low to trigger abiogenesis - through cyanosulfidic chemistry - on them. Only stars with effective temperature >3900 K illuminate their CHZ planets with enough NUV radiation to trigger abiogenesis. Alternatively, colder stars would require a high-energy flaring activity.

Tidal disruption events occur when a star is disrupted by a supermassive black hole, resulting in an elongated stream of gas that partly falls back to pericenter. Due to apsidal precession, the returning stream may collide with itself, leading to a self-crossing shock that launches an outflow. If the black hole spins, this collision may additionally be affected by Lense-Thirring precession that can cause an offset between the two stream components. We study the impact of this effect on the outflow properties by carrying out local simulations of collisions between offset streams. As the offset increases, we find that the geometry of the outflow becomes less spherical and more collimated along the directions of the incoming streams, with less gas getting unbound by the interaction. However, even the most grazing collisions we consider significantly affect the trajectories of the colliding gas, likely promoting subsequent strong interactions near the black hole and rapid disc formation. We analytically compute the dependence of the offset to stream width ratio, finding that even slowly spinning black holes can cause both strong and grazing collisions. We propose that the deviation from outflow sphericity may enhance the self-crossing shock luminosity due to a reduction of adiabatic losses, and cause significant variations of the efficiency at which X-ray radiation from the disc is reprocessed to the optical band depending on the viewing angle. These potentially observable features hold the promise of constraining the black hole spin from tidal disruption events.

Roughly half of the quasi-periodic eruption (QPE) sources in galactic nuclei exhibit a remarkably regular alternating "long-short'' pattern of recurrence times between consecutive flares. In support of previous suggestions, we show that a main-sequence star (brought into the nucleus as an extreme mass-ratio inspiral; EMRI) which passes twice per orbit through the accretion disk of the supermassive black hole (SMBH) on a mildly eccentric inclined orbit, each time shocking and ejecting optically-thick gas clouds above and below the midplane (akin to dual ``supernovae''), naturally reproduces the luminosities, durations, and spectral temperatures of observed QPE flares. Similar to supernova shock break-out from a compact star, inefficient photon production in the collision debris renders the QPE emission much harder than the blackbody temperature, enabling the Comptonized spectrum to stick out from the softer quiescent disk spectrum. Destruction of the star via mass ablation typically limits the lifetime of QPE emission to decades, precluding a long-lived AGN as the source of the gaseous disk. By contrast, a tidal disruption event (TDE) naturally provides a transient gaseous disk on the requisite radial scale, with a rate sufficiently high relative to the EMRI inward migration rate, that a significant fraction of TDEs should host a QPE. This picture is consistent with the X-ray TDE observed several years prior to the QPE emission from GSN 069. Remarkably, a second TDE-like flare was observed from this event, starting immediately after detectable QPE activity ceased; this accretion event could plausibly result from the (partial or complete) destruction of the QPE-generating star triggered by runaway mass-loss. Our model can also be applied to black hole-disk collisions, such as those proposed to generate quasi-periodic flares in the candidate SMBH binary OJ 287.

The baryonic Tully-Fisher relation (BTFR), which connects the baryonic mass of galaxies with their circular velocities, has been validated across a wide range of galaxies, from dwarf galaxies to massive galaxies. Recent studies have found that several ultra-diffuse galaxies (UDGs) deviate significantly from the BTFR, indicating a galaxy population with abnormal dynamical properties. However, such studies were still confined within a small sample size. In this study, we used the 100% complete Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) to investigate the BTFR of 88 HI-rich UDGs (HUDGs), which is the largest UDG sample with dynamical information. We found that the HUDGs form a continuous distribution in the BTFR diagram, with high-velocity galaxies consistent with normal dwarf galaxies at 1 $\sigma$ level, and low-velocity galaxies deviating from the BTFR, in line with that reported in the literature. We point out that the observed deviation may be subject to various selection effects or systemic biases. Nevertheless, we found that the significance of the deviation of HUDGs from the BTFR and TFR are different, i.e., they either deviate from the BTFR or from the TFR. Our result indicates that a high-gas fraction may play an important role in explaining the deviation of HUDGs from BTFR.

We investigate the star formation main sequence (MS) (SFR-M$_{\star}$) down to 10$^{8-9}\mathrm{M}_\odot$ using a sample of 34,061 newly-discovered ultra-faint ($27\lesssim i \lesssim 30$ mag) galaxies at $1<z<3$ detected in the GOODS-N field. Virtually these galaxies are not contained in previous public catalogs, effectively doubling the number of known sources in the field. The sample was constructed by stacking the optical broad-band observations taken by the HST/GOODS-CANDELS surveys as well as the 25 ultra-deep medium-band images gathered by the GTC/SHARDS project. Our sources are faint (average observed magnitudes $<i>\sim28.2$ mag, $<H>\sim27.9$ mag), blue (UV-slope $<\beta>\sim-1.9$), star-forming (rest-frame colors $<U-V>\sim0.10$ mag, $<V-J>\sim0.17$ mag) galaxies. These observational characteristics are identified with young (mass-weighted age $<\mathrm{t_{M-w}}>\sim0.014$ Gyr) stellar populations subject to low attenuations ($<\mathrm{A(V)}>\sim0.30$ mag). Our sample allows us to probe the MS down to $10^{8.0}\,\mathrm{M}_\odot$ at $z=1$ and $10^{8.5}\,\mathrm{M}_\odot$ at $z=3$, around 0.6 dex deeper than previous analysis. In the low-mass galaxy regime, we find an average value for the slope of 0.97 at $1<z<2$ and 1.12 at $2<z<3$. Nearly $\sim$60% of our sample presents stellar masses in the range $10^{6-8}$ M$_\odot$ between $1<z<3$. If the slope of the MS remained constant in this regime, the sources populating the low-mass tail of our sample would qualify as starburst galaxies.

We introduce ASSIST, a software package for ephemeris-quality integrations of test particles. ASSIST is an extension of the REBOUND framework and makes use of its IAS15 integrator to integrate test particle trajectories in the field of the Sun, Moon, planets, and 16 massive asteroids, with the positions of the masses coming from the JPL DE441 ephemeris and its associated asteroid perturber file. The package incorporates the most significant gravitational harmonics and general relativistic corrections. ASSIST also accounts for position- and velocity-dependent non-gravitational effects. The first order variational equations are included for all terms to support orbit fitting and covariance mapping. This new framework is meant to provide an open-source package written in a modern language to enable high-precision orbital analysis and science by the small body community. ASSIST is open source, freely distributed under the GNU General Public license, version 3.

More than 50% of solar-mass stars form in multiple systems. It is therefore crucial to investigate how multiplicity affects the star and planet formation processes at the protostellar stage. We report continuum and C$^{18}$O (2-1) observations of the VLA 1623-2417 protostellar system at 50 au angular resolution as part of the ALMA Large Program FAUST. The 1.3 mm continuum probes the disks of VLA 1623A, B, and W, and the circumbinary disk of the A1+A2 binary. The C$^{18}$O emission reveals, for the first time, the gas in the disk-envelope of VLA 1623W. We estimate the dynamical mass of VLA 1623W, $M_{\rm dyn}=0.45\pm0.08$ M$_{\odot}$, and the mass of its disk, $M_{\rm disk}\sim6\times10^{-3}$ M$_{\odot}$. C$^{18}$O also reveals streamers that extend up to 1000 au, spatially and kinematically connecting the envelope and outflow cavities of the A1+A2+B system with the disk of VLA 1623W. The presence of the streamers, as well as the spatial ($\sim$1300 au) and velocity ($\sim$2.2 km/s) offset of VLA 1623W suggest that either sources W and A+B formed in different cores, interacting between them, or that source W has been ejected from the VLA 1623 multiple system during its formation. In the latter case, the streamers may funnel material from the envelope and cavities of VLA 1623AB onto VLA 1623W, thus concurring to set its final mass and chemical content.

The mid-infrared spectral region provides a unique window into the atmospheric temperature, chemistry, and dynamics of the giant planets. From more than a century of mid-infrared remote sensing, progressively clearer pictures of the composition and thermal structure of these atmospheres have emerged, along with a greater insight into the processes that shape them. Our knowledge of Jupiter and Saturn has benefitted from their proximity and relatively warm temperatures, while the details of colder and more distant Uranus and Neptune are limited, as these planets remain challenging targets. As the timeline of observations continues to grow, an understanding of the temporal and seasonal variability of the giant planets is beginning to develop, with promising new observations on the horizon.

Multi-scale deconvolution is an ill-posed inverse problem in imaging, with applications ranging from microscopy, through medical imaging, to astronomical remote sensing. In the case of high-energy space telescopes, multi-scale deconvolution algorithms need to account for the peculiar property of native measurements, which are sparse samples of the Fourier transform of the incoming radiation. The present paper proposes a multi-scale version of CLEAN, which is the most popular iterative deconvolution method in Fourier space imaging. Using synthetic data generated according to a simulated but realistic source configuration, we show that this multi-scale version of CLEAN performs better than the original one in terms of accuracy, photometry, and regularization. Further, the application to a data set measured by the NASA Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) shows the ability of multi-scale CLEAN to reconstruct rather complex topographies, characteristic of a real flaring event.

Excess infrared flux from white dwarf stars is likely to arise from a dusty debris disk or a cool companion. In this work, we present near-infrared spectroscopic observations with Keck/MOSFIRE, Gemini/GNIRS, and Gemini/Flamingos-2 of seven white dwarfs with infrared excesses identified in previous studies. We confirmed the presence of dust disks around four white dwarfs (Gaia J0611-6931, Gaia J0006+2858, Gaia J2100+2122, and WD 0145+234) as well as two new white dwarf brown dwarf pairs (Gaia J0052+4505 and Gaia J0603+4518). In three of the dust disk systems, we detected for the first time near-infrared metal emissions (Mg I, Fe I, and Si I) from a gaseous component of the disk. We developed a new Markov Chain Monte Carlo framework to constrain the geometric properties of each dust disk. In three systems, the dust disk and the gas disk appear to coincide spatially. For the two brown dwarf white dwarf pairs, we identified broad molecular absorption features typically seen in L dwarfs. The origin of the infrared excess around Gaia J0723+6301 remains a mystery. Our study underlines how near-infrared spectroscopy can be used to determine sources of infrared excess around white dwarfs, which has now been detected in hundreds of systems photometrically.

We test for an anisotropy in the mass of arriving cosmic-ray primaries associated with the galactic plane. The sensitivity to primary mass is obtained through the depth of shower maximum, $X_{\rm max}$, extracted from hybrid events measured over a 14-year period at the Pierre Auger Observatory. The sky is split into distinct on- and off-plane regions using the galactic latitude of each arriving cosmic ray to form two distributions of $X_{\rm max}$, which are compared using an Anderson-Darling 2-samples test. A scan over roughly half of the data is used to select a lower threshold energy of $10^{18.7}\,$eV and a galactic latitude splitting at $|b| = 30^\circ$, which are set as a prescription for the remaining data. With these thresholds, the distribution of $X_{\rm max}$ from the on-plane region is found to have a $9.1 \pm 1.6^{+2.1}_{-2.2}\,$g$\,$cm$^-2$ shallower mean and a $5.9\pm2.1^{+3.5}_{-2.5}\,$g$\,$cm$^-2$ narrower width than that of the off-plane region and is observed in all telescope sites independently. These differences indicate that the mean mass of primary particles arriving from the on-plane region is greater than that of those from the off-plane region. Monte Carlo studies yield a $5.9\times10^{-6}$ random chance probability for the result in the independent data, lowering to a $6.0\times10^{-7}$ post-penalization random chance probability when the scanned data is included. Accounting for systematic uncertainties leads to an indication for anisotropy in mass composition above $10^{18.7}\,$eV with a $3.3\,\sigma$ significance. Furthermore, the result has been newly tested using additional FD data recovered from the selection process. This test independently disfavors the on- and off-plane regions being uniform in composition at the $2.2\,\sigma$ level, which is in good agreement with the expected sensitivity of the dataset used for this test.

Density fluctuations in the ionised interstellar medium have a profound effect on radio pulsar observations, through angular scattering, intensity scintillations, and small changes in time delays from dispersion. Here we show that it is possible to recover the variations in dispersive delays that originate from a dominant scattering region using measurements of the dynamic spectrum of intensity scintillations, provided that the pulsar velocity and scattering region location are known. We provide a theoretical framework for the technique, which involves estimating the phase gradient from the dynamic spectra and integrating that gradient to obtain phase variations. It can be used to search for "extreme scattering events" (ESEs) in pulsars for which precision dispersion delay measurements are not otherwise possible, or to separate true dispersion variations from apparent variability caused by frequency-dependent pulse shape changes. We demonstrate that it works in practice by recovering an ESE in PSR J1603$-$7202, which is known from precision dispersion delay measurements from pulsar timing. For this pulsar, we find that the phase gradients also track the long-term variations in electron column density observed by pulsar timing, indicating that the column density variations and the scattering are dominated by the same thin scattering screen. We identify a sudden increase in the scintillation strength and magnitude of phase gradients over $\sim$days in 2010, indicating a compact structure. A decrease in the electron density in 2012 was associated with persistent phase gradients and preceded a period of decreased scintillation strength and an absence of scintillation arcs.

The progenitor systems and explosion mechanisms responsible for the thermonuclear events observationally classified as Type Ia supernovae are uncertain and difficult to uniquely constrain using traditional observations of Type Ia supernova host galaxies, progenitors, light curves, and remnants. For the subset of thermonuclear events that are prolific producers of iron, we use published theoretical nucleosynthetic yields to identify a set of elemental abundance ratios infrequently observed in metal-poor stars but shared across a range of progenitor systems and explosion mechanisms: [Na,Mg,Co/Fe]<0. We label stars with this abundance signature ``iron-rich metal-poor'' or IRMP stars. We suggest that IRMP stars formed in environments dominated by thermonuclear nucleosynthesis and consequently that their elemental abundances can be used to constrain both the progenitor systems and explosion mechanisms responsible for thermonuclear explosions. We identify three IRMP in the literature and homogeneously infer their elemental abundances. We find that the elemental abundances of BD+80 245, HE 0533--5340, and SMSS J034249.53--284216.0 are best explained by the (double) detonations of sub-Chandrasekhar mass CO white dwarfs. If our interpretation of IRMP stars is accurate, then they should be very rare in globular clusters and more common in the Magellanic Clouds and dwarf spheroidal galaxies than in the Milky Way's halo. We propose that future studies of IRMP stars will quantify the relative occurrences of different thermonuclear event progenitor systems and explosion mechanisms.

Through two-dimensional radiation-hydrodynamical simulations, we investigate the spectral properties of ionized gas irradiated by an active galactic nucleus with a supermassive black hole of 10^7 Msun. For the gas inside the dust-sublimation radius (r ~ 10^{-2} pc), we conduct post-process pseudo-three-dimensional calculations utilizing the spectral synthesis code Cloudy. We show that we can reproduce various broad emission lines in optical and ultraviolet wavelengths. The line profiles change depending on the viewing angles even for a small range from the rational axis, i.e., 5-30 degrees; most lines, such as Halpha, are characterized by a double-peaked profile, reflecting that the emissions are originated in the surface of the rotating disk. By contrast, high-ionization emission lines such as CIV 1549 show a double-peaked profile for a nearly face-on view, as these lines derive from the fast outflowing gas from the disk surface. Our results suggest that some properties of the bright UV-optical emission lines observed in Seyfert-like AGNs can be caused by the radiation-driven fountain flow inside the dust sublimation radius.

Prompted by the many controversial claims involving the NGC 1052 group, including that it hosts two dark matter-free galaxies with overluminous and monochromatic globular cluster (GC) systems, here we map out the large-scale structure (LSS) of GCs over the entire group. To recover the LSS, we use archival optical CFHT imaging data. We recover two GC density maps, one based on universal photometric properties of GCs from simple stellar population models, and one based on the properties of spectroscopically confirmed GCs in DF2 and DF4 (the two dwarf galaxies with overluminous GC populations). Both selection methods reveal overdensities around the massive galaxies in the group, as well as around NGC 1052 itself, that are coincident with the positions of previously identified stellar streams and tidal features. No intragroup GCs are found connecting these structures to any of the dwarf galaxies. We find, however, two other dwarfs in the group hosting GC systems. These include RCP32 with 2 GCs with ages equivalent to the GCs around NGC 1052, and DF9 with 3 GCs with ages similar to the GCs around DF2 and DF4. We conclude that the GC distribution in the group does not strongly support any formation scenario in particular. It favours, nonetheless, scenarios relying on galaxy-galaxy interactions and on the coeval formation of GCs around the DM-free dwarf galaxies. These may include the recently proposed bullet-dwarf formation, as well as high-redshift tidal dwarf galaxy models.

Launch of the Gaia space observatory started a new era in astrometry when the accuracy of star coordinates increased by thousands of times. Significant improvement of accuracy was also expected for the coordinates of the Solar system bodies. Gaia DR3 provided us with the data which could be used to test our expectations. In this work, we refine the orbits of a number of outer planetary satellites using both ground-based and Gaia observations. From thirteen outer satellites observed by Gaia, we chose six to obtain their orbits. Some specific moments in using observations of outer satellites made by Gaia are demonstrated. These pecularities stem from scanning motion of Gaia, in particular from the fact that the accuracy of observations is significantly different along and across the scanning direction. As expected, Gaia observations proved to be more precise than those made from Earth, which results in more accurate satellite ephemerides. We estimate accuracy of the ephemerides of considered satellites for the interval between 1996 and 2030. As astrometric positions published in Gaia DR3 were not corrected for the relativistic light deflection by the Sun, we took into account this effect, which slightly diminished the rms residuals. In addition, relativistic light deflection by the giant planets was estimated, which, as it turned out, can be neglected with the given accuracy of Gaia observations.

The study of blazar jet has been performed for several decades via the VLBI technique, while its generation and propagation stay unclear. In the present work, we compiled a sample of 407 VLBI detected \textit{Fermi} blazars (VFBs) and studied the correlations between apparent velocity (${\rm log}\,\beta_{\rm app}$) and jet/accretion disk properties. We found a positive correlation between $\gamma$-ray luminosity (${\rm log}\,L_{\rm \gamma}$) and ${\rm log}\,\beta_{\rm app}$, the correlation suggests that the apparent motion of jet knot is related to the jet power.

If, after primordial inflation, the universe undergoes a relatively long reheating period, it could present a phase of matter domination supported by the oscillating inflaton field. During this epoch, small perturbations from the inflaton that reenter the cosmological horizon could virialize to form \textit{inflaton} structures. If the primordial overdensities are large enough, their associated inflaton structures could collapse to form primordial black holes (PBHs) [L.E.Padilla, J.C.Hidalgo and K.A.Malik, Phys.Rev.D, vol.106, p.023519, Jul 2022; hereinafter P1]. For this to happen at a considerable rate, the primordial power spectrum should be enhanced at small scales, a feature typically induced in single-field inflation through an ultra-slow roll phase (produced by a nearly-inflection point in the inflationary potential). In this article we consider two specific inflationary potentials that present this nearly-inflection point and we look at the PBH formation rate through the mechanism proposed in P1. We report on constraints to these two specific models from the bounds to PBH abundances. This serves as an illustration of the usefulness of the PBH formation mechanism proposed in P1.

NGC 6302 (The ''Butterfly Nebula'') is an extremely energetic bipolar nebula whose central star is among the most massive, hottest, and presumably rapidly evolving of all central stars of planetary nebulae. Our proper-motion study of NGC 6302, based on excellent HST WFC3 images spanning 11 yr, has uncovered at least four different pairs of expanding internal lobes that were ejected at various times over the past two millennia at speeds ranging from 10 to 600 km s^-1. In addition, we find a pair of off-axis flows in constant motion at 760 +/- 100 km s^-1 within which bright [Fe II] feathers are conspicuous. Combining our results with those previously published, we find that the ensemble of flows has an ionized mass > 0.1 M_sun. The kinetic energy of the ensemble, 10^46 - 10^48 ergs, lies at the upper end of gravity-powered processes such as stellar mergers or mass accretion and is too large to be explained by stellar radiation pressure or convective ejections. The structure and dynamics of the Butterfly Nebula suggests that its central engine has had a remarkable history, and the highly unusual patterns of growth within its wings challenge our current understanding of late stellar mass ejection.

We present a sample of 330 blue edge-on low surface brightness galaxies (ELSBGs). To understand the chemical evolution of LSBGs, we derived the gas-phase abundance and the [$\alpha$/Fe] ratio. Compared with star-forming galaxies, ELSBGs show a flatter trend in the mass-metallicity ($M_*-Z$) relation, suggesting that the oxygen abundance enhancement is inefficient. We focus on 77 ELSBGs with HI data and found the closed-box model can not explain their gas fraction and metallicity relation, implying that infall and/or outflow is needed. We derived the [$\alpha$/Fe] ratio of normal ELSBG ($<$ 10$^{9.5}$M$\odot$) and massive ELSBG ($>=$ 10$^{9.5}$M$\odot$) using single stellar population grids from MILES stellar library. The mean [$\alpha$/Fe] ratios are 0.18 and 0.4 for normal ELSBG and massive ELSBG, respectively. We discussed that the long time-scale of star-formation, and/or metal-rich gas outflow event caused by SNe Ia winds are likely responsible for the $\alpha$-enhancement of massive ELSBGs.

Fast neutrino flavor conversions (FFCs) and collisional flavor instabilities (CFIs) potentially affect the dynamics of core-collapse supernovae (CCSNe) and binary neutron star mergers (BNSMs). Under the assumption of homogeneous neutrinos, we investigate effects of neutrino emission and absorption (EA) by matters through both single and multi-energy numerical simulations with physically motivated setup. In our models, FFCs dominate over CFIs in the early phase, while EA secularly and significantly give impacts on flavor conversions. They facilitate angular swaps, or the full exchange between electron neutrinos ($\nu_e$) and heavy-leptonic neutrinos ($\nu_x$). As a result, the number density of $\nu_x$ becomes more abundant than the case without EA, despite the fact that the isotropization by EA terminates the FFCs earlier. In the later phase, the system approaches new asymptotic states characterized by EA and CFIs, in which rich energy-dependent structures also emerge. Multi-energy effects sustain FFCs and the time evolution of the flavor conversion becomes energy dependent, which are essentially in line with effects of the isoenergetic scattering studied in our previous paper. We also find that $\nu_x$ in the high-energy region convert into $\nu_e$ via flavor conversions and then they are absorbed through charged current reactions, exhibiting the possibility of new path of heating matters.

We report on the result of He abundance analysis of dense circumstellar clumps in the young supernova remnant Cassiopeia A. These clumps, which are called quasi-stationary flocculi (QSFs), are known from previous optical studies to be enriched in He along with N, but the degree of He overabundance relative to H has remained uncertain. For several QSFs with near-infrared spectroscopic data, we have analyzed their He I 1.083 $\mu$m/Pa$\gamma$ ratios together with the ratios of [Fe II] lines by using the Raymond shock code. According to our analysis, He is overabundant relative to H by a factor of $\lesssim 3$ in most of these QSFs. This He abundance of QSFs is consistent with the previous conclusion from the N overabundance that QSFs were ejected when a substantial amount of the H envelope of the progenitor star had been stripped off. We discuss the mass-loss history of the progenitor star and the origin of QSFs.

To understand the evolution of the Moon, we numerically modeled mantle convection and magmatism in a two-dimensional polar rectangular mantle. Magmatism occurs as an upward permeable flow of magma generated by decompression melting through the convecting matrix. The mantle is assumed to be initially enriched in heat-producing elements (HPEs) and compositionally dense ilmenite-bearing cumulates (IBC) at its base. Here, we newly show that magma generation and migration play a crucial role in the calculated volcanic and radial expansion/contraction history. Magma is generated in the deep mantle by internal heating for the first several hundred million years. A large volume of the generated magma ascends to the surface as partially molten fingers and plumes driven by melt-buoyancy to cause a volcanic activity and radial expansion of the planet with the peak at 3.5-4 Gyr ago. Eventually, however, the planet begins to radially contract when the mantle solidifies by cooling from the surface boundary. As the mantle is cooled, the activity of partially molten plumes declines but continues for billions of years after the peak because some basal materials enriched in the dense IBC components hold HPEs. The calculated volcanic and radial expansion/contraction history is consistent with the observed history of the Moon. Our simulations suggest a substantial fraction of the mantle was solid, and there was a basal layer enriched in HPEs and the IBC components at the beginning of the history of the Moon.

In this paper the dependence of the spectra of cosmic ray nuclei on the charges of nuclei was studied, according to the data of the NUCLEON space experiment. First, we studied the dependence of the spectral index of magnetic rigidity spectra on the charge for abundant nuclei. Secondly, for the charge range $Z=9\div20$, the differences in the total spectra of rare odd and abundant even nuclei were studied. Using the GALPROP package, the inverse problem of CR propagation from a source (near supernova) to an observer was solved, a component-by-component spectrum in the source was reconstructed, and it was shown that a systematic change in the spectral index in the source exist. It is supposed that this change may be interpreted as incomplete ionization of cosmic rays at the stage of acceleration in the supernova remnant shock. The ratio of the total spectra of magnetic rigidity for low-abundance odd and abundant even nuclei from the charge range $Z=9\div20$ is obtained, and it was shown that the spectra of odd rare nuclei are harder than the stpectra of abundat even nuclei in the rigidity range 300--10000~GV.

The accurate thermoelastic properties and thermal conductivity are crucial in understanding the thermal evolution of the Martian core. A fitting method based on the ab initio calculated pressure-volume-temperature data is proposed in the formulation of the equation of state with high accuracy, by which the pressure and temperature dependent thermoelastic properties can be directly calculated by definitions. The ab initio results show that the liquid Fe0.75S0.25 under Martian core condition is thoroughly in the ferromagnetic state, without existing spin crossover. The liquid Fe0.75S0.25 in magnetic calculation owns a low thermal conductivity (21~23 W/m/K) when compared with non-magnetic calculation at the same state. Based on the Insight estimated and ab initio calculated properties of the Martian core, the iron snow model is verified when the current temperature at the core-mantle boundary is below the core melting temperature, and the simply secular cooling model is verified on the contrary.

Proto-neutron stars are born when a highly evolved and massive star collapses under gravity. In this paper, we investigate the spontaneous scalarization in proto-neutron stars. Based on the scalar tensor theory of gravity as well as the physical conditions in proto-neutron star, we examine the structure of proto-neutron star. To describe the fluid in proto-neutron star, we utilize $SU(2)$ chiral sigma model and the finite temperature extension of the Brueckner-Bethe-Goldstone quantum many-body theory in the Brueckner-Hartree-Fock approximation. Here, we apply the equation of state of proto-neutron stars considering different cases i.e. hot pure neutron matter and hot $\beta$-stable neutron star matter without neutrino trapping as well as with neutrino trapping. The effects of temperature and entropy of proto-neutron stars on the star structure are also studied. Our results confirm that the spontaneous scalarization is affected by different physical conditions in proto-neutron stars.

The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with two- and three-dimensional (2D and 3D) particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field $B$ and plasma density~$n$, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the $n\mbox{-}B$ space are developed and used to visualize and quantify the effects of compression. The $n\mbox{-}B$ histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic (QED) effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

We present the results of our study of the X-ray pulsar IGR J21343+4738 based on NuSTAR, Swift, and SRG observations in the wide energy range 0.3 - 79 keV. The absence of absorption features in the energy spectra of the source, both averaged and phase-resolved ones, has allowed us to estimate the upper and lower limits on the magnetic field of the neutron star in the binary system, $B<2.5\times10^{11}$G and $B>3.4 \times 10^{12}$G, respectively. The spectral and timing analyses have shown that IGR J21343+4738 has all properties of a quasi-persistent X-ray pulsar with a pulsation period of $322.71\pm{0.04}$s and a luminosity $L_{x} \simeq3.3$ $\times10^{35}$erg s$^{-1}$. The analysis of the long-term variability of the object in X-rays has confirmed the possible orbital period of the binary system $\sim 34.3$ days previously detected in the optical range.

In this paper, we use SPARC's HSBs, and LSBs galaxies to verify two issues. The first one is related to one claim of \citep{Donato} D09, namely: is the DM surface density (DMsd) a constant universal quantity, equal to $\log{(\rm \Sigma/M_\odot pc^{-2})}=2.15 \pm 0.2$, or does it depend on the baryon surface density of the system? The second one, is based on a MOND prediction that for HSBs the DMsd is constant, and equal to $\log{(\rm \Sigma/M_\odot pc^{-2})}=2.14$, while for LSBs the surface density is not constant and takes values that are smaller than for HSBs and D09 prediction \citep{Milgrom2009}. We find that HSBs shows a constant DMsd vs magnitude as in D09, and a constant DMsd vs $\Sigma_{\rm eff}$ as in MOND prediction, for HSBs with $\Sigma_{\rm eff}>200 L_\odot/pc^2$, and $\Sigma_{\rm eff}>300 L_\odot/pc^2$. However, the value of the DMsd is larger, $\Sigma \simeq 2.61$ (in the case of the DMsd-magnitude with $\Sigma_{\rm eff}>300 L_\odot/pc^2$), and $\Sigma \simeq 2.54$ (in the case of the surface DMsd-surface brightness with $\Sigma_{\rm eff}>200 L_\odot/pc^2$). This value slightly depends on the threshold to determine wheter a galaxy is HSB. In the case of LSBs, for $\Sigma_{\rm eff}<100 L_\odot/pc^2$, and $\Sigma_{\rm eff}<25 L_\odot/pc^2$, the surface density vs magnitude, for lower magnitudes, is approximately equal to that predicted by D09, but several galaxies, for magnitude $M>-17$, have smaller values than those predicted by D09. The DMsd vs $\Sigma_{\rm eff}$ shows a similar behavior in qualitative, but not quantitative, agreement with MOND predictions. In summary, in the case of HSBs both D09 and MOND are in qualitative, but not quantitative, agreement with the data. In the case of LSBs D09 is mainly in disagreement with the data, and MOND only in qualitative agreement with them.

Ultrahigh-energy cosmic rays are almost exclusively detected through extensive air showers, which they initiate upon interaction with the atmosphere. The longitudinal development of these air showers can be directly observed using fluorescence detector telescopes, such as those employed at the Pierre Auger Observatory or the Telescope Array. In this article, we discuss the properties of the Greisen function, which was initially derived as an approximate solution to the electromagnetic cascade equations, and its ability to describe the longitudinal shower profiles. We demonstrate that the Greisen function can be used to describe longitudinal air-shower profiles, even for hadronic air showers. Furthermore we discuss the possibility to discriminate between hadrons and photons from the shape of air-shower profiles using the Greisen function.

The flatness of galaxy rotation curves at large radii is generally considered to be a significant piece of evidence in support of the existence of dark matter. Several studies have claimed that post-Newtonian corrections to the Newtonian equations of galaxy dynamics remove the need for dark matter. A few recent studies have examined these claims, and identified errors in their reasoning. We add to this critique by giving what we consider to be particularly simple and transparent description of the errors made in these post-Newtonian calculations, some of which were of a rather technical nature, others more fundamental, e.g. the loss of the correct relativistic scaling, promoting small corrections to order unity changes. Our work reinforces the orthodoxy that post-Newtonian effects are indeed too small to significantly alter galactic rotation curves, and will hopefully serve as a useful guide for others, pointing out subtle errors that one might inadvertently make in such calculations.

Context: The search for extrasolar planets around the nearest M-dwarfs is a crucial step towards identifying the nearest Earth-like planets. One of the main challenges in this search is that M-dwarfs can be magnetically active and stellar activity can produce radial velocity (RV) signals that could mimic those of a planet. Aims: We aim to investigate whether the 2.2 day period observed in optical RVs of the nearby active M-dwarf star Gl 388 (AD Leo) is due to stellar activity or to a planet which co-rotates with the star as suggested in the past. Methods: We obtained quasi-simultaneous optical RVs of Gl 388 from 2019 to 2021 with SOPHIE (R$\sim$75k) at the OHP in France, and near-IR RV and Stokes V measurements with SPIRou at the CFHT (R$\sim$70k). Results: The SOPHIE RV time-series displays a periodic signal with 2.23$\pm$0.01 days period and 23.6$\pm$0.5 m/s amplitude, which is consistent with previous HARPS observations obtained in 2005-2006. The SPIRou RV time-series is flat at 5 m/s rms and displays no periodic signals. RV signals of amplitude higher than 5.3 m/s at a period of 2.23 days can be excluded with a confidence level higher than 99%. Using the modulation of the longitudinal magnetic field (Bl) measured with SPIRou, we derive a stellar rotation period of 2.2305$\pm$0.0016 days. Conclusions: SPIRou RV measurements provide solid evidence that the periodic variability of the optical RVs of Gl 388 is due to stellar activity rather than to a co-rotating planet. The magnetic activity nature of the optical RV signal is further confirmed by the modulation of Bl with the same period. The SPIRou campaign on Gl 388 demonstrates the power of near-IR RV to confirm or infirm planet candidates discovered in the optical around active stars. SPIRou observations reiterate how effective spectropolarimetry is at determining the stellar rotation period.

A historical planet-search on a sample of 1647 nearby southern main sequence stars has been ongoing since 1998 with the CORALIE spectrograph at La Silla Observatory, with a backup subprogram dedicated to the monitoring of binary stars. We review 25 years of CORALIE measurements and search for Doppler signals consistent with stellar or brown dwarf companions to produce an updated catalog of both known and previously unpublished binary stars in the planet-search sample, assessing the binarity fraction of the stellar population and providing perspective for more precise planet-search in the binary sample. We perform new analysis on the CORALIE planet-search sample radial velocity measurements, searching for stellar companions and obtaining orbital solutions for both known and new binary systems. We perform simultaneous radial velocity and proper motion anomaly fits on the subset of these systems for which Hipparcos and Gaia astrometry measurements are available, obtaining accurate estimates of true mass for the companions. We find 218 stars in the CORALIE sample to have at least one stellar companion, 130 of which are not yet published in the literature and for which we present orbital solutions. The use of proper motion anomaly allow us to derive true masses for the stellar companions in 132 systems, which we additionally use to estimate stability regions for possible planetary companions on circumprimary or circumbinary orbits. Finally, we produce detection limit maps for each star in the sample and obtain occurrence rates of $0.43^{+0.23}_{-0.11}\%$ and $12.69^{+0.87}_{-0.77}\%$ for brown dwarf and stellar companions respectively in the CORALIE sample.

Machine learning can be a powerful tool to discover new signal types in astronomical data. We here apply it to search for long-duration transient gravitational waves triggered by pulsar glitches, which could yield physical insight into the mostly unknown depths of the pulsar. Current methods to search for such signals rely on matched filtering and a brute-force grid search over possible signal durations, which is sensitive but can become very computationally expensive. We develop a method to search for post-glitch signals on combining matched filtering with convolutional neural networks, which reaches similar sensitivities to the standard method at false-alarm probabilities relevant for practical searches, while being significantly faster. We specialize to the Vela glitch during the LIGO-Virgo O2 run, and set upper limits on the gravitational-wave strain amplitude from the data of the two LIGO detectors for both constant-amplitude and exponentially decaying signals.

Recent observations of a small sample of repeating Fast Radio Bursts (FRBs) have revealed a periodicity in their bursting activity that may be suggestive of a binary origin for the modulation. We set out to explore the scenario where a subset of repeating FRBs originates in binary systems hosting a highly energetic neutron star and a massive companion star, akin to $\gamma$-ray binaries and young High-Mass X-ray Binaries. Firstly, we specifically focus on the host galaxy properties and binary formation rates. Subsequently, we investigate the expected evolution of the rotation and dispersion measure in this scenario, the predicted birth-site offsets, and the origin of the persistent radio emission observed in a subset of these systems. The host galaxies for repeating FRBs favour the formation of neutron star-massive star binary systems but any conclusive evidence will require future discoveries and localizations of FRBs. The birth rate of high-mass X-ray binaries, used as a proxy for all considered binaries, significantly exceeds the estimated rate of FRBs, which can be explained if only a small subset of these systems produce FRBs. We show that under simple assumptions, we can reproduce the DM and RM evolution that is seen in a subset of repeating FRBs. We also discuss the possibility of detecting a persistent radio source associated with the FRB due to an intra-binary shock between companion star wind and either the pulsar wind or giant magnetar flares. The observed long-term luminosity stability of the Persistent Radio Sources is most consistent with a giant flare-powered scenario. However, this explanation is highly dependent on the magnetic field properties of the neutron star. With these explorations, we have aimed to provide a framework to discuss future FRB observations in the context of neutron star-massive star binary scenarios.

The High Resolution Telescope (HRT) of the Polarimetric and Helioseismic Imager on board the Solar Orbiter spacecraft (SO/PHI) and the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) both infer the photospheric magnetic field from polarised light images. SO/PHI is the first magnetograph to move out of the Sun--Earth line and will provide unprecedented access to the Sun's poles. This provides excellent opportunities for new research wherein the magnetic field maps from both instruments are used simultaneously. We aim to compare the magnetic field maps from these two instruments and discuss any possible differences between them. We used data from both instruments obtained during Solar Orbiter's inferior conjunction on 7 March 2022. The HRT data were additionally treated for geometric distortion and degraded to the same resolution as HMI. The HMI data were re-projected to correct for the $3^{\circ}$ separation between the two observatories. SO/PHI-HRT and HMI produce remarkably similar line-of-sight magnetograms, with a slope coefficient of $0.97$, an offset below $1$ G, and a Pearson correlation coefficient of $0.97$. However, SO/PHI-HRT infers weaker line-of-sight fields for the strongest fields. As for the vector magnetic field, SO/PHI-HRT was compared to both the $720$-second and $90$-second HMI vector magnetic field: SO/PHI-HRT has a closer alignment with the $90$-second HMI vector. In the weak signal regime ($< 600$ G), SO/PHI-HRT measures stronger and more horizontal fields than HMI, very likely due to the greater noise in the SO/PHI-HRT data. In the strong field regime ($\gtrsim 600$ G), HRT infers lower field strengths but with similar inclinations (a slope of $0.92$) and azimuths (a slope of $1.02$). The slope values are from the comparison with the HMI $90$-second vector.

Fast radio bursts (FRBs) have been found in great numbers but the physical mechanism of these sources is still a mystery. The redshift evolutions of the FRB energy distribution function and the volumetric rate shed light on revealing the origin of the FRBs. However, such estimations rely on the dispersion measurement (DM)-redshift ($z$) relationship. A few of FRBs detected recently show large excess DM beyond the expectation from the cosmological and Milky Way contributions, which indicates large spread of DM from their host galaxies. In this work, we adopt the lognormal distributed $\rm DM_{host}$ model and estimate the energy function using the non-repeating FRBs selected from the Canadian Hydrogen Intensity Mapping Experiment (CHIME)/FRB Catalog 1. By comparing the lognormal distributed $\rm DM_{host}$ model to the constant $\rm DM_{host}$ model, the FRB energy function results are consistent within the measurement uncertainty. We also estimate the volumetric rate of the non-repeating FRBs in three different redshift bins. The volumetric rate shows that the trend is consistent with the stellar-mass density redshift evolution. Since the lognormal distributed $\rm DM_{host}$ model increases the measurement errors, the inference of FRBs tracking the stellar-mass density is nonetheless undermined.

The population of black holes inferred from the detection of gravitational waves by the LIGO-Virgo-KAGRA collaboration has revealed interesting features in the properties of black holes in the universe. We analyze the GWTC-3 dataset assuming the detected black holes have both astrophysical and primordial origins. In particular, we consider primordial black holes forming from vacuum bubbles that nucleate during inflation, with their mass distribution described by a broken power law. We find that more than half of the events could come from primordial black hole mergers. Astrophysical black holes are mainly responsible for the peak in mass distribution at $\sim 10M_\odot$ indicated by GWTC-3; whereas primordial black holes are responsible for the massive black holes, as well as the peak at $\sim 30M_\odot$. We also discuss the implications on the primordial black hole formation mechanism and the underlying inflationary model.

Strongly magnetized neutron stars are popular candidates for producing detectable electromagnetic and gravitational-wave signals. A rapid density increase in a neutron star core could also trigger the phase transition from hadrons to deconfined quarks and form a hybrid star. This formation process could release a considerable amount of energy in the form of gravitational waves and neutrinos. Hence, the formation of a magnetized hybrid star is an interesting scenario for detecting all these signals. These detections may provide essential probes for the magnetic field and composition of such stars. Thus far, a dynamical study of the formation of a magnetized hybrid star has yet to be realized. Here, we investigate the formation dynamics and the properties of a magnetized hybrid star through dynamical simulations. We find that the maximum values of rest-mass density and magnetic field strength increase slightly and these two quantities are coupled in phase during the formation. We then demonstrate that all microscopic and macroscopic quantities of the resulting hybrid star vary dramatically when the maximum magnetic field strength goes beyond a threshold of $\sim 5 \times 10^{17}$ G but they are insensitive to the magnetic field below this threshold. Specifically, the magnetic deformation makes the rest-mass density drop significantly, suppressing the matter fraction in the mixed phase. Therefore, this work provides a solid support for the magnetic effects on a hybrid star, so it is possible to link observational signals from the star to its magnetic field configuration.

We study the colour and star formation rates of paired galaxies in filaments and sheets using the EAGLE simulations. We find that the major pairs with pair separation $<50$ kpc are bluer and more star forming in filamentary environments compared to those hosted in sheet-like environments. This trend reverses beyond a pair separation of $\sim 50$ kpc. The interacting pairs with larger separations ($>50$ kpc) in filaments are on average redder and low-star forming when compared to those embedded in sheets. The galaxies in filaments and sheets may have different stellar mass and cold gas mass distributions. Using a Kolmogorov-Smirnov test, we find that for paired galaxies with pair separation $<50$ kpc, there are no significant differences in these properties in sheets and filaments. The filaments transport gas towards the cluster of galaxies. Some earlier studies find preferential alignment of galaxy pairs with filament axis. Such alignment of galaxy pairs may lead to different gas accretion efficiency in galaxies residing in filaments and sheets. We propose that the enhancement of star formation rate at smaller pair separation in filaments is caused by the alignment of galaxy pairs. A recent study with the SDSS data (Das et al., 2023) reports the same findings. The confirmation of these results by the EAGLE simulations suggests that the hydrodynamical simulations are powerful theoretical tools for studying the galaxy formation and evolution in the cosmic web.

Past years have seen various publications attempting to explain the apparent clustering features of trans-Neptunian objects, the most popular explanation being an unconfirmed "Planet 9". The recently proposed Uranus Orbiter and Probe mission by NASA's Planetary Science and Astrobiology Decadal Survey could offer the opportunity to precisely determine Planet 9's sky location and mass by carefully monitoring ranging data during the interplanetary cruise. We use Monte Carlo-Markov Chain methods to reconstruct simulated spacecraft trajectories in a simplified solar system model containing Planet 9, providing an estimate of the mission's localisation capacity depending on sky location, Earth-spacecraft Doppler link noise level and data collection rate. We characterise the noise via the Allan deviation $\sigma_{\rm A}$, scaled to the Cassini-era value $\sigma_{\rm A}^{\rm \scriptscriptstyle Cass} = 3 \times 10^{-15}$, finding that daily measurements of the spacecraft position can lead to $\sim$0.2 deg$^2$ localisation of Planet 9 (assuming $M_9 = 6.3 M_{\oplus}$, $d_9 = 460$AU). As little as a 3-fold improvement in $\sigma_{\rm A}$ drastically decreases the sky localisation area size to $\sim$0.01 deg$^2$. Thus, we showcase that a future Uranus mission carries a significant potential also for non-Uranian science.

In order to make accurate inferences about the solar interior using helioseismology, it is essential to understand all the relevant physical effects on the observations. One effect to understand is the (complex-valued) ratio of the horizontal to vertical displacement of the p- and f-modes at the height at which they are observed. Unfortunately, it is impossible to measure this ratio directly from a single vantage point, and it has been difficult to disentangle observationally from other effects. In this paper we attempt to measure the ratio directly using 7.5 hours of simultaneous observations from the Polarimetric and Helioseismic Imager on board Solar Orbiter and the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. While image geometry problems make it difficult to determine the exact ratio, it appears to agree well with that expected from adiabatic oscillations in a standard solar model. On the other hand it does not agree with a commonly used approximation, indicating that this approximation should not be used in helioseismic analyses. In addition, the ratio appears to be real-valued.

We present the discovery of 48 new and the analysis of 55, including 7 previously discovered, dwarf galaxy candidates around the giant spiral galaxy NGC~2997 using deep $BVI$ images from the KMTNet Supernova Program. Their $V$-band central surface brightness and total absolute magnitudes range between 20.3--26.7 mag arcsec$^{-2}$ and --(8.02--17.69) mag, respectively, while the $I$-band effective radii are between 0.14 and 2.97 kpc. We obtain $\alpha$ $\simeq$ --1.43 $\pm$ 0.02 for the faint-end slope of their luminosity function, comparable to previously measured values but shallower than theoretical predictions based on $\Lambda$CDM models. The distance-independent distributions of their mass and color from the host galaxy NGC~2997 suggest that the group could be dynamically young, prior to the development of significant mass segregation or radial color gradients. The systematically bluer colors of the brighter candidates than the fainter ones indicate higher star formation activities in brighter members. We suggest that the higher-mass dwarf galaxies in the group have maintained star-formation activities by effectively retaining gas content, while environmental quenching is only effective for the lower-mass galaxies. The interpretation of early evolutionary stage of this group is also consistent with the overall morphological distribution of the dwarf galaxy candidates showing a lack of morphologically evolved candidates but a plethora of irregularly shaped ones. Our detection rate of dwarf galaxy candidates in the NGC~2997 group and their inferred star formation activities are largely comparable to those found in Milky Way analog systems from the SAGA survey within the magnitude limit M$_{V}$ $\lesssim$ --13 mag, as well as those found in the ELVES survey.

At the lowest masses, the distinction between brown dwarfs and giant exoplanets is often blurred and literature classifications rarely reflect the deuterium burning boundary. Atmospheric characterisation may reveal the extent to which planetary formation pathways contribute to the population of very-low mass brown dwarfs, by revealing if their abundance distributions differ from those of the local field population or, in the case of companions, their primary stars. The T8 dwarf Ross 458c is a possible planetary mass companion to a pair of M dwarfs, and previous work suggests that it is cloudy. We here present the results of the retrieval analysis of Ross 458c, using archival spectroscopic data in the 1.0 to 2.4 micron range. We test a cloud free model as well as a variety of cloudy models and find that the atmosphere of Ross 458c is best described by a cloudy model (strongly preferred). The CH4/H2O is higher than expected at 1.97 +0.13 -0.14. This value is challenging to understand in terms of equilibrium chemistry and plausible C/O ratios. Comparisons to thermochemical grid models suggest a C/O of ~ 1.35, if CH4 and H2O are quenched at 2000 K, requiring vigorous mixing. We find a [C/H] ratio of +0.18, which matches the metallicity of the primary system, suggesting that oxygen is missing from the atmosphere. Even with extreme mixing, the implied C/O is well beyond the typical stellar regime, suggesting a either non-stellar formation pathway, or the sequestration of substantial quantities of oxygen via hitherto unmodeled chemistry or condensation processes.

As a part of the ``Systematic KMTNet Planetary Anomaly Search" series, we report five new planets (namely, OGLE-2016-BLG-1635Lb, MOA-2016-BLG-532Lb, KMT-2016-BLG-0625Lb, OGLE-2016-BLG-1850Lb, and KMT-2016-BLG-1751Lb) and one planet candidate (KMT-2016-BLG-1855), which were found by searching $2016$ KMTNet prime fields. These $buried$ planets show a wide range of masses from Earth--class to Super--Jupiter--class, and are located in both the disk and the bulge. The ultimate goal of this series is to build a complete planet sample. Because our work provides a complementary sample to other planet detection methods, which have different detection sensitivities, our complete sample will help us to obtain a better understanding of planet demographics in our Galaxy.

We revisit the possibility of using cosmological observations to constrain models that involve interactions between neutrinos and dark matter. We show that small-scale measurements of the cosmic microwave background with a few per cent accuracy are critical to uncover unique signatures from models with tiny couplings that would require a much higher sensitivity at lower multipoles, such as those probed by the Planck satellite. We analyze the high-multipole data released by the Atacama Cosmology Telescope, both independently and in combination with Planck and Baryon Acoustic Oscillation measurements, finding a compelling preference for a non-vanishing coupling, $\log_{10}u_{\nu \textrm{DM}}=-5.20^{+1.2}_{-0.74}$ at 68% CL. This aligns with other CMB-independent probes, such as Lyman-$\alpha$. We illustrate how this coupling could be accounted for in the presence of dark matter interactions with a sterile neutrino.

The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where in each arm there is a `xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogenic temperature. Here, we examine the scientific perspectives under possible variations of this reference design. We perform a detailed evaluation of the science case for a single triangular geometry observatory, and we compare it with the results obtained for a network of two L-shaped detectors (either parallel or misaligned) located in Europe, considering different choices of arm-length for both the triangle and the 2L geometries. We also study how the science output changes in the absence of the low-frequency instrument, both for the triangle and the 2L configurations. We examine a broad class of simple `metrics' that quantify the science output, related to compact binary coalescences, multi-messenger astronomy and stochastic backgrounds, and we then examine the impact of different detector designs on a more specific set of scientific objectives.

We consider the relational approach to construct gauge-invariant observables in cosmological perturbation theory using synchronous coordinates. We construct dynamical synchronous coordinates as non-local scalar functionals of the metric perturbation in the fully non-linear theory in an arbitrary gauge. We show that the observables defined in this dynamical coordinate system are gauge-independent, and that the full perturbed metric has the expected form in these coordinates. Our construction generalises the familiar synchronous gauge in linearised gravity, widely used in cosmological perturbation theory, to the non-linear theory. We also work out the expressions for the gauge-invariant Einstein equations, sourced either by an ideal fluid or a scalar field up to second order in perturbation theory, and give explicit expressions for the Hubble rate -- as measured by synchronous observers or by observers co-moving with the matter field -- up to that order. Finally, we consider quantised linear perturbations around Minkowski and de Sitter backgrounds, and compute the two-point function of the gauge-invariant metric perturbation in synchronous coordinates, starting with two-point function in a general linear covariant gauge. Although the gauge-fixed two-point function contains gauge modes, we show that the resulting gauge-invariant two-point function only contains the physical tensor modes and it is thus positive, i.e., it has a spectral representation.

The detection of black-hole binaries at high redshifts is a cornerstone of the science case of third-generation gravitational-wave interferometers. The star-formation rate peaks at z~2 and decreases by orders of magnitude by z~10. Any confident detection of gravitational waves from such high redshifts would imply either the presence of stars formed from pristine material originating from cosmological nucleosynthesis (the so-called population III stars), or black holes that are the direct relics of quantum fluctuations in the early Universe (the so-called primordial black holes). Crucially, detecting sources at cosmological distances does not imply inferring that sources are located there, with the latter posing more stringent requirements. To this end, we present two figures of merit, which we refer to as "z-z plot" and "inference horizon", that quantify the largest redshift one can possibly claim a source to be beyond. We argue that such inference requirements, in addition to detection requirements, should be investigated when quantifying the scientific payoff of future gravitational-wave facilities.

This paper investigates exact solutions of cosmological interest in fractional cosmology. Given $\mu$, the order of the fractional derivative, and $w$, the matter equation of state, we present particular exact power-law solutions. We discuss the exact general solution of the system obtained by solving a Riccati Equation, where the solution for the scale factor is a combination of power-law. Using cosmological data, we estimate the free parameters $(\alpha_0, \mu)$, where $H_{0}=100\frac{\text{km/s}}{\text{Mpc}}h$, and $\alpha_0:=t_0 H_0 = \frac{1}{6} \left(9 -2 \mu +\sqrt{8 \mu (2 \mu -9)+105}\right)(1+ 2 \epsilon_0)$, is the current age parameter. The joint analysis with data from SNe Ia + OHD leads to $h=0.684_{-0.027}^{+0.031}$, $\mu=1.840_{-0.773}^{+1.446}$ and $\epsilon_0=\left(1.213_{-1.057}^{+0.482}\right)\times 10^{-2}$, where the best-fit values are calculated at $3\sigma$ CL. On the other hand, these best-fit values lead to an age of the Universe with a value of $t_0=\alpha_0/H_0=25.62_{-4.46}^{+6.89}\;\text{Gyrs}$, a current deceleration parameter of $q_{0}=-0.37_{-0.11}^{+0.08}$, both at $3\sigma$ CL, and a current matter density parameter of $\Omega_{m,0}=0.531_{-0.260}^{+0.195}$ at $1\sigma$ CL. Finding a Universe roughly twice older as the one of $\Lambda$CDM is a distinction of Fractional Cosmology. Focusing our analysis on these results, we can conclude that the region in which $\mu>2$ is not ruled out by observations. This region of a parameter is relevant because, in the absence of matter, fractional cosmology gives a power-law solution $a(t)= \left(t/t_0\right)^{\mu-1}$, which is accelerated for $\mu>2$. We present a fractional origin model that leads to an accelerated state without appealing to $\Lambda$ or Dark Energy.

The collisionless Boltzmann equation (CBE) is a fundamental equation that governs the dynamics of a broad range of astrophysical systems from space plasma to star clusters and galaxies. It is computationally expensive to integrate the CBE directly in a phase space, and thus the applications to realistic astrophysical problems have been limited so far. Recently, Todorova \& Steijl (2020) proposed an efficient quantum algorithm for solving the CBE with a significantly reduced computational complexity. We extend the method to perform quantum simulations that follow the evolution of self-gravitating systems. We first run a 1+1 dimensional test calculation of free streaming motion on 64$\times$64 grids using 13 simulated qubits and validate our method. We then perform simulations of Jeans collapse, and compare the result with analytic and linear theory calculations. We propose a direct method to generate initial conditions as well as a method to retrieve necessary information from a register of multiple qubits. Our simulation scheme achieves $\mathcal{O}(N_v^3)$ less computational complexity than the classical method, where $N_v$ is the number of discrete velocity grids per dimension. It will thus allow us to perform large-scale CBE simulations on future quantum computers.

We discuss motion of a binary system around a supermassive black hole. Using Fermi-Walker transport, we construct a local inertial reference frame and set up a Newtonian binary system. Assuming a circular geodesic observer around a Schwarzschild black hole, we write down the equations of motion of a binary. Introducing a small acceleration of the observer, we remove the interaction terms between the center of mass (CM) of a binary and its relative coordinates. The CM follows the observer's orbit, but its motion deviates from an exact circular geodesic. We first solve the relative motion of a binary system, and then find the motion of the CM by the perturbation equations with the small acceleration. We show that there appears the Kozai-Lidov (KL) oscillations when a binary is compact and the initial inclination is larger than a critical angle. In a hard binary system, KL oscillations are regular, whereas in a soft binary system, oscillations are irregular both in period and in amplitude, although stable. We find an orbital flip when the initial inclination is large. As for the motion of the CM, the radial deviations from a circular orbit become stable oscillations with very small amplitude.

It is argued, as a working hypothesis, that "normal" and dark matter interactions can only be T and CP violating. One way to implement this idea is to consider that time reversal in dark matter is implemented, not by an antiunitary operator, but by a unitary operator. It is shown how this occurs naturally in the context of complex spacetime with an extended symmetry group.

Primordial black holes (PBHs) forming out of the collapse of enhanced cosmological perturbations provide access to the early Universe through their associated observational signatures. In particular, enhanced cosmological perturbations collapsing to form PBHs are responsible for the generation of a stochastic gravitational-wave background (SGWB) induced by second-order gravitational interactions, usually called scalar induced gravitational waves (SIGWs). This SGWB is sensitive to the underlying gravitational theory; hence it can be used as a novel tool to test the standard paradigm of gravity and constrain possible deviations from general relativity. In this work, we study the aforementioned GW signal within modified teleparallel gravity theories, developing a formalism for the derivation of the GW spectral abundance within any form of gravitational action. At the end, working within viable $f(T,\phi)$ models without matter-gravity couplings, and accounting for the effect of mono-parametric $f(T)$ gravity at the level of the source and the propagation of the tensor perturbations, we show that the respective GW signal is indistinguishable from that within GR. Interestingly, we find that in order to break the degeneracy between different $f(T)$ theories through the portal of SIGWs one should necessarily consider non-minimal matter-gravity couplings at the level of the gravitational action.

We examine electron-transport coefficients in magnetized hot and dense electron-ion plasma relevant in binary neutron star merger simulation. We calculate electrical and thermal conductivities in low density, high temperature, highly magnetized plasma of binary neutron star mergers where quantum oscillatory behavior of electrons emerge. For pronounced thermodynamic effects, we consider zeroth Landau level population of electrons for the calculation of conductivity. We solve Boltzmann equation in presence of magnetic field to obtain the dissipative components of electrical and thermal conductivities. The dissipative coefficients are formulated considering frequency dependent dynamical screening in the quantized electron-ion scattering rate. Numerical estimations show that the effect of dynamical screening of photon propagator on electrical and thermal conductivities is pronounced. We observe that dynamical screening reduces the maxima of both the electrical and thermal conductivities by factors of thirty one and twenty respectively leading to a reduction in the corresponding time scales of these coefficients. The common scaling factor between electrical and thermal conductivity is also observed to follow cubic relationship with temperature violating Wiedemann-Franz law.