Papers for Friday, Apr 19 2024

Context. Most ground-based planet search direct imaging campaigns use angular differential imaging, which distorts the signal from extended sources like protoplanetary disks. In the case PDS 70, a young system with two planets found within the cavity of a protoplanetary disk, obtaining a reliable image of both planets and disk is essential to understanding planet-disk interactions. Aims. Our goals are to reveal the true intensity of the planets and disk without self-subtraction effects for the first time, search for new giant planets beyond separations of 0.1" and to study the morphology of the disk shaped by two massive planets. Methods. We present YJHK-band imaging, polarimetry, and spatially resolved spectroscopy of PDS 70 using near-simultaneous reference star differential imaging, also known as star-hopping. We created a radiative transfer model of the system to match the near-infrared imaging and polarimetric data, along with sub-millimeter imaging from ALMA. Furthermore, we extracted the spectra of the planets and the disk and compared them. Results. We find that the disk is quite flared with a scale height of ~15% at the outer edge of the disk at ~90 au, similar to some disks in the literature. The gap inside of ~50 au is estimated to have ~1% of the dust density of the outer disk. The Northeast outer disk arc seen in previous observations is likely the outer lip of the flared disk. Abundance ratios of grains estimated by the modeling indicate a shallow grain-size index > -2.7, instead of the canonical -3.5. There is both vertical and radial segregation of grains. Planet c is well separated from the disk and has a spectrum similar to planet b, clearly redder than the disk spectra. Planet c is possibly associated with the sudden flaring of the disk starting at ~50 au. No new planets > 5 Mj were found.

JWST is discovering a large population of $z>4$ supermassive black holes (SMBHs) that are overmassive with respect to the stellar content of their hosts. A previous study developed a physical model to interpret this overmassive population as the result of quasar feedback acting on a compact host galaxy. In this Note, we apply this model to JADES GN 1146115, a dormant supermassive black hole at $z=6.7$ whose mass is $\sim40\%$ of the host's mass in stars and accreting at $\sim2\%$ of the Eddington limit. The host has been forming stars at the low rate of $\sim 1 \, \rm M_\odot \,yr^{-1}$ for the past $\sim 100$ Myr. Our model suggests that this galactic system is on the verge of a resurgence of global star formation activity. This transition comes after a period of domination by the effect of its overmassive black hole, whose duration is comparable to typical quasar lifetimes.

Very massive stars (VMSs) formed via a sequence of stellar collisions in dense star clusters have been proposed as the progenitors of massive black hole seeds. VMSs could indeed collapse to form intermediate-mass black holes (IMBHs), which would then grow by accretion to become the supermassive black holes observed at the centers of galaxies and powering high-redshift quasars. Previous studies have investigated how different cluster initial conditions affect the formation of a VMS, including mass segregation, stellar collisions, and binaries, among others. In this study, we investigate the growth of VMSs with a new grid of Cluster Monte Carlo (CMC) star cluster simulations -- the most expansive to date. The simulations span a wide range of initial conditions, varying the number of stars, cluster density, stellar initial mass function (IMF), and primordial binary fraction. We find a gradual shift in the mass of the most massive collision product across the parameter space; in particular, denser clusters born with top-heavy IMFs provide strong collisional regimes that form VMSs with masses easily exceeding 1000 solar masses. Our results are used to derive a fitting formula that can predict the typical mass of a VMS formed as a function of the star cluster properties. Additionally, we study the stochasticity of this process and derive a statistical distribution for the mass of the VMS formed in one of our models, recomputing the model 50 times with different initial random seeds.

Stellar bars are common morphological structures in the local Universe: according to optical and NIR surveys, they are present in about two-thirds of disc galaxies. These elongated structures are also believed to play a crucial role in secular evolutionary processes, since they are able to efficiently redistribute gas, stars and angular momentum within their hosts, although it remains unclear whether they enhance or suppress star formation. A useful tool to investigate such an ambiguity is the Main Sequence (MS) relation, which tightly links stellar mass ($M_{\star}$) and star formation rate (SFR). The main goal of this work is to explore star formation processes in barred galaxies, in order to assess the relevance of possible bar quenching effects on the typical log-linear trend of the resolved MS. To this purpose, we carry out a spatially resolved analysis on sub-kpc scales for a sample of six nearby barred galaxies. Multi-wavelength photometric data (from far-UV to far-IR) are collected from the DustPedia database and a panchromatic Spectral Energy Distribution (SED) fitting procedure is applied on square apertures of fixed angular size (8" $\times$ 8"), making use of the magphys code. For each galaxy we obtain the distributions of stellar mass and star formation rate surface densities and relate them in the $\log \Sigma_{\star}$ - $\log \Sigma_{\rm SFR}$ plane deriving the spatially resolved MS relation. Although significant galaxy-to-galaxy variations are in place, we infer the presence of a common anti-correlation track in correspondence with the bar-hosting region, which shows systematically lower values of SFR. Such a central quiescent signature can be interpreted as the result of a bar-driven depletion of gas reservoirs and a consequent halt of star formation. This seems to point in the direction of an inside-out quenching scenario.

We present ALMA band-7 observations of J2057-0030, a multi-component merger system at $z$ $\sim$ 4.68 spanning at least 50 kpc in size, using the [CII] $\lambda$157.74 $\mu$m line and underlying far-infrared (FIR) continuum. We find two main components, the quasar (QSO) and a dusty star-forming galaxy (DSFG), both detected in [CII] and continuum emission as well as multiple neighboring clumps detected only in [CII]. Three of these clumps form a (tidal) tail that extends from the QSO in a straight direction to the west, covering a projected distance of $\sim$ 10 kpc. This perturbed morphology, added to a spatial distance of $\sim$ 20 kpc and a velocity offset of $\Delta{v}$ = 68 km s$^{-1}$ between the QSO and the DSFG, strongly supports a merging scenario. By fitting a spectral energy distribution model to the continuum data, we estimate star formation rates of $\approx$ 402 $M_{\odot}$ yr$^{-1}$ for the QSO host and $\approx$ 244 $M_{\odot}$ yr$^{-1}$ for the DSFG, which locate them on or close to the main sequence of star-forming galaxies. The J2057-0030 QSO was selected for being one of the brightest unobscured quasars at its redshift while presenting a rather modest star formation rate. Based on a commonly accepted paradigm regarding the formation of quasars, this result is expected for a quasar that has already passed an obscured phase of rapid star formation during a major merger. However, we see that the merger event in this system is far from being finished, and it is rather likely somewhere between the first pericenter and subsequent close passages. This is presumably another case of a high-$z$ quasar residing in a high-density environment with a companion obscured galaxy.

We present a sample of over 100 highly X-ray luminous galaxy clusters at $z\sim$ 0.5-0.9, discovered by the extended Massive Cluster Survey (eMACS) in ROSAT All-Sky Survey (RASS) data. Follow-up observations of a subset at higher resolution and greater depth with the Chandra X-ray Observatory are used to map the gaseous intra-cluster medium, while strong-gravitational-lensing features identified in Hubble Space Telescope imaging allow us to constrain the total mass distribution. We present evidence of the exceptional gravitational-lensing power of these massive systems, search for substructure along the line of sight by mapping the radial velocities of cluster members obtained through extensive ground-based spectroscopy, and identify dramatic cases of galaxy evolution in high-density cluster environments. The available observations of the eMACS sample presented here provide a wealth of insights into the properties of very massive clusters ($\gtrsim 10^{15} M_\odot$) at $z > 0.5$, which act as powerful lenses to study galaxies in the very distant Universe. We also discuss the evolutionary state, galaxy population, and large-scale environment of eMACS clusters and release to the community all data and science products to further the understanding of the first generation of truly massive clusters to have formed in the Universe.

As the hunt for an Earth-like exoplanets has intensified in recent years, so has the effort to characterise and model the stellar signals that can hide or mimic small planetary signals. Stellar variability arises from a number of sources, including granulation, supergranulation, oscillations and activity, all of which result in quasi-periodic or stochastic behaviour in photometric and/or radial velocity observations. Traditionally, the characterisation of these signals has mostly been done in the frequency domain. However, the recent development of scalable Gaussian process regression methods makes direct time-domain modelling of stochastic processes a feasible and arguably preferable alternative, obviating the need to estimate the power spectral density of the data before modelling it. In this paper, we compare the two approaches using a series of experiments on simulated data. We show that frequency domain modelling can lead to inaccurate results, especially when the time sampling is irregular. By contrast, Gaussian process regression results are often more precise, and systematically more accurate, in both the regular and irregular time sampling regimes. While this work was motivated by the analysis of radial velocity and photometry observations of main sequence stars in the context of planet searches, we note that our results may also have applications for the study of other types of astrophysical variability such as quasi-periodic oscillations in X-ray binaries and active galactic nuclei variability.

Though free-floating planets (FFPs) that have been ejected from their natal star systems may outpopulate their bound counterparts in the terrestrial-mass range, they remain one of the least explored exoplanet demographics. Due to their negligible electromagnetic emission at all wavelengths, the only observational technique able to detect these worlds is gravitational microlensing. Microlensing by terrestrial-mass FFPs induces rare, short-duration magnifications of background stars, requiring high-cadence, wide-field surveys to detect these events. The Transiting Exoplanet Survey Satellite (TESS), though designed to detect close-bound exoplanets via the transit technique, boasts a cadence as short as 200 seconds and has monitored hundreds of millions of stars, making it well-suited to search for short-duration microlensing events as well. We have used existing data products from the TESS Quick-Look Pipeline (QLP) to perform a preliminary search for FFP microlensing candidates in 1.3 million light curves from TESS Sector 61. We find one compelling candidate associated with TIC-107150013, a source star at $d_s = 3.194$ kpc. The event has a duration $t_E = 0.074^{+0.002}_{-0.002}$ days and shows prominent finite-source features ($\rho = 4.55^{+0.08}_{-0.07}$), making it consistent with an FFP in the terrestrial-mass range. This exciting result indicates that our ongoing search through all TESS sectors has the opportunity to shed new light on this enigmatic population of worlds.

The enrichment history of $r$-process elements has been imprinted on the stellar abundances that change in accordance with increasing metallicity in galaxies. Close examination of the [Eu/Fe] feature caused by stars in nearby galaxies, including the Large Magellanic Cloud (LMC), shows its perplexity. The decreasing trend of the [Eu/Fe] feature is followed by a nearly constant value; this trend is generally attributed to an onset of the delayed Fe release from type Ia supernovae (SNe Ia), which is the same interpretation of the [$\alpha$/Fe] feature. However, this feature appears in the LMC at [Fe/H] of approximately -0.7, which is significantly higher than that for the [alpha/Fe] case ($\approx$ -2). This result potentially indicates the presence of an overlooked property of the $r$-process site that remains unseen in the study of the Milky Way. Here, we propose that this [Eu/Fe]-knee feature is created by a fade-out of core-collapse SNe producing $r$-process elements; these elements along with neutron star mergers (NSMs) promote the $r$-process enrichment under the condition for this specific SNe such that their occurrence is limited to a low-metallicity environment. This metallicity threshold for the occurrence rate of $r$-process SNe at a subsolar is nearly identical to that for long gamma-ray bursts whose origin may be connected to fast-rotating massive stars. Moreover, we reason that the contribution of Eu from NSMs is crucial to maintain a high [Eu/Fe] at an early stage in dwarf galaxies by a balance with Fe from SNe Ia; both enrichments via NSMs and SNe Ia proceed with similar delay time distributions.

Rapid formation of supermassive black holes occurs in dense nuclear star clusters that are initially gas-dominated. Stellar-mass black hole remnants of the most massive cluster sink into the core, where a massive runaway black hole forms as a consequence of combined effects of repeated mergers and Eddington-limited gas accretion. The associated gravitational-wave signals of high-redshift extreme mass-ratio inspirals are a unique signature of the nuclear star cluster scenario.

The redshifted 21-cm signal from the epoch of reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically-averaged EoR 21-cm signal power spectrum. The goal is to study these scale-dependent features at different stages of reionization using numerical simulations and build a source model-independent framework to probe the properties of the intergalactic medium using EoR 21-cm signal power spectrum measurements. Under the assumption of high spin temperature, we modelled the redshift evolution of the ratio of EoR 21-cm brightness temperature power spectrum and the corresponding density power spectrum using an ansatz consisting of a set of redshift and scale-independent parameters. This set of eight parameters probes the redshift evolution of the average ionization fraction and the quantities related to the morphology of the ionized regions. We have tested this ansatz on different reionization scenarios generated using different simulation algorithms and found that it is able to recover the redshift evolution of the average neutral fraction within an absolute deviation $\lesssim 0.1$. Our framework allows us to interpret 21-cm signal power spectra in terms of parameters related to the state of the IGM. This source model-independent framework can efficiently constrain reionization scenarios using multi-redshift power spectrum measurements with ongoing and future radio telescopes such as LOFAR, MWA, HERA, and SKA. This will add independent information regarding the EoR IGM properties.

The presence of an abundant population of low frequency photons at high redshifts (such as a radio background) can source leading order effects on the evolution of the matter and spin temperatures through rapid free-free absorptions. This effect, known as soft photon heating, can have a dramatic impact on the differential brightness temperature, $\Delta T_{\rm b}$, a central observable in $21$cm cosmology. Here, we introduce a semi-analytic framework to describe the dynamics of soft photon heating, providing a simplified set of evolution equations and a useful numerical scheme which can be used to study this generic effect. We also perform quasi-instantaneous and continuous soft photon injections to elucidate the different regimes in which soft photon heating is expected to impart a significant contribution to the global $21$cm signal and its fluctuations. We find that soft photon backgrounds produced after recombination with spectral index $\gamma > 3.0$ undergo significant free-free absorption, and therefore this heating effect cannot be neglected. The effect becomes stronger with steeper spectral index, and in some cases the injection of a synchrotron-like spectrum ($\gamma = 3.6$) can suppress the amplitude of $\Delta T_{\rm b}$ relative to the standard model prediction, making the global $21$cm signal even more difficult to detect in these scenarios.

Andromeda XVIII is an isolated dwarf galaxy 579 kpc away from the nearest large galaxy, M31. It is a candidate ``backsplash galaxy'' that might have been affected by a close passage to M31. We present new Keck/DEIMOS spectroscopy of Andromeda XVIII to assess the likelihood that it is a backsplash galaxy. We estimated the velocities, metallicities ([Fe/H]), and $\alpha$-enhancements ([$\alpha$/Fe]) for 56 probable members. We estimated Andromeda XVIII's mean heliocentric velocity, rotation velocity, position angle of the rotation axis, and velocity dispersion using maximum likelihood coupled with a Monte Carlo Markov chain (MCMC). There is no evidence for bulk rotation, though sub-populations might be rotating. The mean heliocentric velocity is -337.2 km s$^{-1}$. The line-of-sight velocity relative to M31 is -36 km s$^{-1}$ (approaching), which is lower than the escape velocity from M31. Based on the abundances of 38 stars with low errors ($\delta [Fe/H] < 0.3$) compared to a total of 56 probable members, parameters for the simplest chemical evolution models were estimated using maximum likelihood coupled with an MCMC. The metallicity distribution is inconsistent with these models due to a sharp metal-rich cut-off. Together, the metallicity distribution and the mean velocity are consistent with a sudden interruption of star formation. One possible cause for this quenching may be rapid gas loss due to ram pressure stripping during a close passage by M31 in the past. We also consider a past major merger as another possible cause.

Since their first discovery in the late 1960s, Gamma-ray bursts have attracted an exponentially growing interest from the international community due to their central role in the most highly debated open questions of the modern research of astronomy, astrophysics, cosmology, and fundamental physics. These range from the intimate nuclear composition of high density material within the core of ultra-dense neuron stars, to stellar evolution via the collapse of massive stars, the production and propagation of gravitational waves, as well as the exploration of the early Universe by unveiling first stars and galaxies (assessing also their evolution and cosmic re-ionization). GRBs have stimulated in the past $\sim$50 years the development of cutting-edge technological instruments for observations of high energy celestial sources from space, leading to the launch and successful operations of many different scientific missions (several of them still in data taking mode nowadays). In this review, we provide a brief description of the GRB-dedicated missions from space being designed and developed for the future. The list of these projects, not meant to be exhaustive, shall serve as a reference to interested readers to understand what is likely to come next to lead the further development of GRB research and associated phenomenology.

Surfaces of carbonaceous asteroids (C-complex) have shown diverse contrasting spectral variations, which may be related to space weathering. We performed laser irradiation experiments on CI and CM simulant material under vacuum to mimic the spectral alteration induced by micrometeorite impacts. We used in situ ultraviolet-visible and near-infrared reflectance spectroscopy to analyze spectral alterations in response to pulsed laser irradiation, as well as scanning electron microscopy and x-ray photoelectron spectroscopy to search for microstructural and compositional changes. Laser irradiation causes an increase in spectral slope (reddening) and a decrease in the albedo (darkening), and these changes are stronger in the ultraviolet-visible region. These spectral changes are likely driven by the excess iron found in the altered surface region, although other factors, such as the observed structural changes, may also contribute. Additionally, while the 0.27~${\mu}m$ band appears relatively stable under laser irradiation, a broad feature at 0.6~${\mu}m$ rapidly disappears with laser irradiation, suggesting that space weathering may inhibit the detection of any feature in this spectral region, including the 0.7~${\mu}m$ band, which has typically been used an indicator of hydration. Comparing our laboratory results with optical spectrophotometry observations of C-complex asteroids, we find that the majority of objects are spectrally red and possess colors that are similar to our irradiated material rather than our fresh samples. Furthermore, we also find that ``younger'' and ``older'' C-complex families have similar colors, suggesting that the space weathering process is near-equal or faster than the time it takes to refresh the surfaces of these airless bodies.

The kinetic energy in tidal flows, when converted into heat, can affect the internal structure of a star and shift its location on a color-magnitude diagram from that of standard models. In this paper we explore the impact of injecting heat into stars with masses near the main sequence turnoff mass (1.26 $M_\odot$) of the open cluster M67. The heating rate is obtained from the tidal shear energy dissipation rate which is calculated from first principles by simultaneously solving the equations that describe orbital motion and the response of a star's layers to the gravitational, Coriolis, centrifugal, gas pressure and viscous forces. The stellar structure models are computed with MESA. We focus on the effects of injecting heat in pulses lasting 0.01 Gyr, a timeframe consistent with the synchonization timescale in binary systems. We find that the location of the tidally perturbed stars in the M67 color-magnitude diagram is shifted to significantly higher luminosities and effective temperatures than predicted by the standard model isochrone and include locations corresponding to some of the Blue Straggler Stars. Because tidal heating takes energy from the orbit causing it to shrink, Blue Straggler Stars could be merger or mass-transfer progenitors as well as products of these processes.

Cyclotron resonant scattering features (CRSFs) are the absorption features in the X-ray spectra of strongly magnetized accretion neutron stars (NSs), which are probably the most reliable probe to the surface magnetic fields of NSs. The high mass X-ray binary GX 301--2 exhibits a very wide, variable and complicated CRSF in the average spectra, which should be two absorption lines based on NuStar and Insight-HXMT observations. With the Insight-HXMT frequent observations, we performed the phase-resolved spectroscopy and confirmed two cyclotron absorption lines in the phase-resolved spectra, with their centroid energy ratio $\sim 1.6-1.7$ in the super-critical luminosity case. A major hindrance in understanding those CRSFs is the very poorly constrained magnetic inclination angle, which is also a fundamental property of a NS and key to understanding the emission characteristics of a pulsar. Comparing the phase-resolved CRSF with simulated X-ray spectra, the magnetic inclination angle is found to be $\gtrsim 70^{\circ}$, i.e., nearly orthogonal between the NS's spin and magnetic axies. The implications of an orthogonal rotator and magnetic structure evolution in the accreting X-ray binary are also discussed.

Dust grains play a significant role in several astrophysical processes, including gas/dust dynamics, chemical reactions, and radiative transfer. Replenishment of small-grain populations is mainly governed by fragmentation during pair-wise collisions between grains. The wide spectrum of fragmentation outcomes, from complete disruption to erosion and/or mass transfer, can be modelled by the general non-linear fragmentation equation. Efficiently solving this equation is crucial for an accurate treatment of the dust fragmentation in numerical modelling. However, similar to dust coagulation, numerical errors in current fragmentation algorithms employed in astrophysics are dominated by the numerical over-diffusion problem -- particularly in 3D hydrodynamic simulations where the discrete resolution of the mass density distribution tends to be highly limited. With this in mind, we have derived the first conservative form of the general non-linear fragmentation with a mass flux highlighting the mass transfer phenomenon. Then, to address cases of limited mass density resolution, we applied a high-order discontinuous Galerkin scheme to efficiently solve the conservative fragmentation equation with a reduced number of dust bins. An accuracy of 0.1 -1% is reached with 20 dust bins spanning a mass range of 9 orders of magnitude.

Principal regular satellites of gas giants are thought to be formed by the accumulation of solid materials in circumplanetary disks (CPDs). While there has been significant progress in the study of satellite formation in CPDs, details of the supply of satellite building blocks to CPDs remain unclear. We performed orbital integration of solid particles in the protoplanetary disk (PPD) approaching a planet, considering the gas drag force using the results of three-dimensional hydrodynamical simulations of a local region around the planet. We investigated planetary-mass dependence of the capture positions and capture rates of dust particles accreting onto the CPD. We also examined the degree of dust retention in accreting gas onto the CPD, which is important for determining the ratio of dust-to-gas inflow rates, a key parameter in satellite formation. We found that the degree of dust retention increases with increasing planetary mass for a given dust scale height in the PPD. In the case of a small planet ($M_{\rm p}=0.2M_{\rm Jup}$), most particles with insufficient initial altitudes in the PPD are isolated from the gas in the accreting region. On the other hand, in the case of a massive planet ($M_{\rm p}=1M_{\rm Jup}$), dust particles can be coupled to the vertically accreting gas, even when the dust scale height is about $10-30$\% of the gas scale height. The results of this study can be used for models of dust delivery and satellite formation in the CPDs of gas giants of various masses, including exoplanets.

Oscillating neutron stars are sources of continuous gravitational waves. We study analytically the excitation of stellar oscillations by the mechanical impact on the stellar surface of ''clumps'' of stochastically accreted matter. We calculate the waveform and spectrum of the gravitational wave signal emitted by the accretion-driven pulsations. Results are generated for an idealised model of a nonrotating, unmagnetised, one-component star with uniform polytropic index $n_{\rm poly}$ assuming Newtonian gravity and the Cowling approximation. We find that the excited mode amplitudes grow with increasing $n_{\rm poly}$ and mode order $n$. The gravitational wave signal forms a sequence of amplitude-modulated packets for $n_{\rm poly}=1$, lasting $\sim 10^{-3}$s after each impact. The gravitational wave strain increases with increasing $n_{\rm poly}$, but decreases with increasing $n$ and increasing multipole order $l$ for $n_{\rm poly}=1$. In the observing band of current long-baseline interferometers, $g$-modes emit higher, narrower peaks in the amplitude spectral density than $f$- and $p$-modes, with the highest peaks reaching $\sim 10^{-26}$Hz$^{-1/2}$ for modes with damping time $\tau_{nl} \sim 10^{8}$yr. The root-mean-square strain $h_{\text{rms}}$, calculated by summing over modes with $2\leq l\leq4$ and $\tau_{nl} \leq 10^{8}$yr, spans the range $10^{-33} \leq h_{\text{rms}} \leq 10^{-32}$ for $1\leq n_{\text{poly}}\leq 2$.

By using the \textit{Fermi}/GBM data of the X-ray bursts (XRBs) of SGR J1935+2154, we investigate the temporal clustering of the bursts and the cumulative distribution of the waiting time and fluence/flux. It is found that the bursts occurring in the episode hosting FRB 20200428 have obviously shorter waiting times than those in the other episodes. The general statistical properties of the XRBs further indicate they could belong to a self-organized critical (SOC) system (e.g., starquakes), making them very similar to the earthquake phenomena. Then, according to a unified scaling law between the waiting time and energy of the earthquakes as well as their aftershocks, we implement an analogy analysis on the XRBs and find that the FRB episode owns more dependent burst events than the other episodes. It is indicated that the FRB emission could be produced by the interaction between different burst events, which could correspond to a collision between different seismic/Alfven waves or different explosion outflows. Such a situation could appear when the magnetar enters into a global intensive activity period.

We performed a $^{12}$CO- and $^{13}$CO($J=1-0$)-line study of the "Brick" (G0.253+0.016) in the Galactic Centre (GC) by analyzing the archival data obtained with the Nobeyama 45-m telescope. We present kinematics and molecular gas distributions in the longitude-velocity diagrams, and suggest that the Brick is located along the GC Arm I in the central molecular zone (CMZ) in front of the GC, which yields a distance of 8 kpc and GC radius 0.2 kpc. The major and minor-axis diameters of the Brick are $D_x\times D_y=8.4 \ {\rm pc} \times 4.1 \ {\rm pc}$, and the scale radius is $r_{\rm bri}=\sqrt{D_x D_y}=2.96$ pc. The molecular mass inferred from the CO-line integrated intensity is $M_{\rm brixco} \sim 5.1\times 10^4 M_\odot$ for a conversion factor $X_{\rm gc}=1.0\times 10^{20}$ H$_2$ cm$^{-2}$ [K km s$^{-1}$]$^{-1}$, while the Virial mass for the velocity dispersion of $\sigma_v=10.0 $ km s $^{-1}$ is calculated to be $M_{\rm brivir}\sim 6.8 \times 10^4 \ M_\odot$, which yields a new conversion factor of $X_{\rm bri} =1.3\times 10^{20}$ H$_2$ cm$^{-2}$ [K km s$^{-1}$]$^{-1}$. No thermal radio emission indicative of HII region and star formation is found in radio-continuum archive. The Brick's center has a cavity surrounded by a spherical molecular bubble of radius $r_{\rm bub}=1.85$ pc and mass $\sim 1.7\times 10^4 M_\odot$ expanding at $v_{\rm exp}=10$ km s$^{-1}$ with kinetic energy of $E_0\sim 1.7\times 10^{49}$ erg. If the bubble is approximated by an adiabatic spherical shock wave, its age is $t\sim 2/5 r_{\rm bub}/v_{\rm exp}\sim 7.2\times 10^4$ y. We suggest that the bubble will be a dark supernova remnant buried in the dense molecular cloud. The Brick, therefore, experienced massive-star formation followed by a supernova explosion more than $\sim 10^5$ y ago.

Flows driven by photons have been studied for almost a century, and a quantitative description of the radiative forces on atoms and ions is important for understanding a wide variety of systems with outflows and accretion disks, such as active galactic nuclei. Quantifying the associated forces is crucial to determining how these outflows enable interactive mechanisms within these environments, such as AGN feedback. The total number of spectral lines in any given ion of the outflow material must be tabulated in order to give a complete characterization of this force. Here we provide calculations of the dimensionless line force multiplier for AGN environments. For a wide array of representative AGN sources, we explicitly calculate the photoionization balance at the proposed wind-launching region above the accretion disk, compute the strength of the line-driving force on the gas, and revisit and formalize the role of the commonly-used ionization parameter $\xi$ in ultimately determining the line-driving force. We perform these computations and analyses for a variety of AGN central source properties, such as black hole mass, initial wind velocity, and number density. We find that, while useful, the ionization parameter provides an incomplete description of the overall ionization state of the outflow material. We use these findings to provide an updated method for calculating the strength of the radiative line-driving using both the X-ray spectral index $\Gamma_X$ and the ionization parameter.

This work aims to measure the mass accretion rate, the accretion luminosity, and more generally the physical conditions of the warm emitting gas in the inner disk of the very low-mass star 2MASS-J16053215-1933159. We investigate the source mid-infrared spectrum for atomic and molecular hydrogen line emission. We present the full James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) Medium Resolution Spectrometer (MRS) spectrum of the protoplanetary disk around the very low-mass star 2MASS-J16053215-1933159 from the MINDS GTO program, previously shown to be abundant in hydrocarbon molecules. We analyzed the atomic and molecular hydrogen lines in this source by fitting one or multiple Gaussian profiles. We then built a rotational diagram for the H2 lines to constrain the rotational temperature and column density of the gas. Finally, we compared the observed atomic line fluxes to predictions from two standard emission models. We identify five molecular hydrogen pure rotational lines and 16 atomic hydrogen recombination lines. The spectrum indicates optically thin emission for both species. We use the molecular hydrogen lines to constrain the mass and temperature of the warm emitting gas. The HI (7-6) recombination line is used to measure the mass accretion rate and luminosity onto the central source. HI recombination lines can also be used to derive the physical properties of the gas using atomic recombination models. The JWST-MIRI MRS observations for the very low-mass star 2MASS-J16053215-1933159 reveal a large number of emission lines, many originating from atomic and molecular hydrogen because we are able to look into the disk warm molecular layer. Their analysis constrains the physical properties of the emitting gas and showcases the potential of JWST to deepen our understanding of the physical and chemical structure of protoplanetary disks

Growing observations of temporal, spectral, and polarization properties of fast radio bursts (FRBs) indicate that the radio emission of the majority of bursts is likely produced inside the magnetosphere of its central engine, likely a magnetar. We revisit the idea that FRBs are generated via coherent inverse Compton scattering (ICS) off low-frequency X-mode electromagnetic waves (fast magnetosonic waves) by bunches at a distance of a few hundred times of the magnetar radius. Following findings are revealed: 1. Crustal oscillations during a flaring event would excite kHz Alfv\'en waves. Fast magnetosonic waves with the same frequency can be generated directly or be converted from Alfv\'en waves at a large radius, with an amplitude large enough to power FRBs via the ICS process. 2. The cross section increases rapidly with radius and significant ICS can occur at $r \gtrsim 100 R_\star$ with emission power much greater than the curvature radiation power but still in the linear scattering regime. 3. The low-frequency fast magnetosonic waves naturally redistribute a fluctuating relativistic plasma in the charge-depleted region to form bunches with the right size to power FRBs. 4. The required bunch net charge density can be sub-Goldreich-Julian, which allows a strong parallel electric field to accelerate the charges, maintain the bunches, and continuously power FRB emission. 5. This model can account for a wide range of observed properties of repeating FRB bursts, including high degrees of linear and circular polarization and narrow spectra as observed in many bursts from repeating FRB sources.

In June 1976, a pristine meteorite stone weighing approximately 1 kg, fully covered with a fresh black fusion crust, was collected on a mountain road in the high-altitude Alpine environment. The recovery took place while clearing the remnants of a snow avalanche, 2 km northwest of Ischgl in Austria. Subsequent to its retrieval, the specimen remained in the finder's private residence without undergoing any scientific examination or identification until 2008, when it was brought to the University of Innsbruck. The sample was classified as a well-preserved LL6 chondrite, with a W0 weathering grade, implying a relatively short time between the meteorite fall and its retrieval. To investigate the potential connection between the Ischgl meteorite and a recorded fireball event, we have reviewed all documented fireballs ever photographed by German fireball camera stations. This examination led us to identify the fireball EN241170 observed in Germany by ten different European Network stations on the night of November 23/24, 1970, as the most likely candidate. We employed state-of-the-art techniques to reconstruct the fireball's trajectory, and to reproduce both its luminous and dark flight phases in detail. We find that the determined strewn field and the generated heat map closely align with the recovery location of the Ischgl meteorite. Furthermore, the measured radionuclide data reported here indicate that the pre-atmospheric size of the Ischgl meteoroid is consistent with the mass estimate inferred from our deceleration analysis along the trajectory. Our findings strongly support the conclusion that the Ischgl meteorite originated from the EN241170 fireball, effectively establishing it as a confirmed meteorite fall. This discovery enables to determine, along with the physical properties, also the heliocentric orbit and cosmic history of the Ischgl meteorite.

Recent simulations have shown that asymmetries in the ejecta distribution of supernova remnants (SNRs) may be a reflection of asymmetries left over from the initial supernova explosion. Thus, SNR studies provide a vital means for testing and constraining model predictions in relation to the distribution of heavy elements, which are key to improving our understanding of the explosion mechanisms in Type Ia supernovae. The use of a novel blind source separation method applied to the megasecond X-ray observations of the historic Kepler and Tycho supernova remnants has revealed maps of the ejecta distribution. These maps are endowed with an unprecedented level of detail and clear separations from the continuum emission. Our method also provides a three-dimensional (3D) view of the ejecta by individually disentangling red- and blueshifted spectral components associated with images of the Si, S, Ar, Ca, and Fe emission. This approach provides insights into the morphology of the ejecta distribution in those two remnants. Those mappings have allowed us to thoroughly investigate the asymmetries in the intermediate-mass elements and Fe distribution in two Type Ia supernova remnants. We also compared the results with the core-collapse Cassiopeia A remnant, which we had studied previously. The images obtained confirm, as expected for Type Ia SNRs, that the Fe distribution is mostly closer to the core than that of intermediate-mass elements. They also highlight peculiar features in the ejecta distribution, such as the Fe-rich southeastern knot in Tycho.

Solar flares are accompanied by an enhanced emission of electromagnetic waves from the radio up to the gamma-ray range. The associated hard X-ray (HXR) and microwave radiation is generated by energetic electrons, which carry a substantial part of the energy released during a flare. The flare is generally understood as a manifestation of magnetic reconnection in the corona. The so-called standard CSHKP model is one of the most widely accepted models for eruptive flares. The solar flare on September 10, 2017 offers a unique opportunity to study this model. The observations from the Expanded Owens Valley Solar Array (EOVSA) show that 1.6x10^4 electrons with energies >300 keV were generated in the flare region. There are signatures in solar radio and extreme ultraviolet observations as well as numerical simulations that a termination shock (TS) appears in the magnetic reconnection outflow region. Electrons accelerated at the TS can be considered to generate the loop-top HXR sources. In contrast to previous studies, we investigate whether the heating of the plasma at the TS provides enough relativistic electrons needed for the HXR and microwave emission observed during the X8.2 solar flare on September 10, 2017. We studied the heating of the plasma at the TS by evaluating the jump in the temperature across the shock by means of the Rankine-Hugoniot relationships under coronal circumstances measured during that event. The part of relativistic electrons was calculated in the heated downstream region. In the magnetic reconnection outflow region, the plasma is strongly heated at the TS. Thus, there are enough energetic electrons in the tail of the electron distribution function needed for the microwave and HXR emission observed during that event. The generation of relativistic electrons at the TS is a possible mechanism to explain the enhanced microwave and HXR radiation emitted during flares.

Three-dimensional magnetic nulls are the points where magnetic field vanishes and are preferential sites for magnetic reconnection: a process which converts magnetic energy into heat and accelerates charged particles along with a rearrangement of magnetic field lines. In the solar corona, the reconnections manifest as coronal transients including solar flares, coronal mass ejections and coronal jets. The nulls are generally found to be collocated with complex active regions on the solar photosphere. Extrapolation of magnetic field from corresponding photospheric magnetogram indicate an abundance of these nulls in solar atmosphere. Nevertheless, their generation is still not well understood. Recently, Maurya et al. (2023) have demonstrated magnetic reconnection to be a cause for generation and annihilation of magnetic nulls through magnetohydrodynamics simulation, where the initial magnetic field is idealized to have a single radial null. This article further extends the study in a more realistic scenario where the initial magnetic field is constructed by extrapolating photospheric magnetogram data and hence, incorporates field line complexities inherent to a complex active region. For the purpose, the active region NOAA 11977 hosting a C6.6 class flare is selected. The simulation is initiated using non-force-free extrapolated magnetic field from the photospheric vector magnetogram at around 02:48:00 UT on 17 February 2014, 16 minutes before the flare peak. The generation, annihilation and dynamics of nulls are explored by a complimentary usage of trilinear null detection technique and tracing of magnetic field line dynamics. It is found that the nulls can spontaneously generate/annihilate in pairs while preserving the topological degree and can have observational implications like footpoint brightenings. Magnetic reconnection is found to be the cause of such generation and annihilation.

Electromagnetic wave lensing, a common physical phenomenon recognized in visible light for centuries, finds extensive applications in manipulating light in optical systems such as telescopes and cameras. Magnetohydrodynamic wave is a common perturbation phenomenon in the corona. By using high spatio-temporal resolution observations from the Solar Dynamics Observatory, here, we report the observation of a magnetohydrodynamic wave lensing in the highly ionized and magnetized coronal plasma, where quasi-periodic wavefronts emanated from a flare converged at a specific point after traversing a coronal hole. The entire process resembles an electromagnetic wave lensing from the source to the focus. Meanwhile, the magnetohydrodynamic wave lensing is well reproduced through a magnetohydrodynamic numerical simulation with full spatio-temporal resolution. We further investigate potential applications for coronal seismology, as the lensing process encodes information on the Alfv\'en speed, in conjunction with favorable geometric and density variations.

We present a method for a non-linear asteroseismic inversion suitable for gravity-mode pulsators and apply it to slowly pulsating B-type (SPB) stars. Our inversion method is based on the iterative improvement of a parameterised static stellar structure model, which in turn is based on constraints from the observed oscillation periods. We present tests to demonstrate that the method is successful in recovering the properties of artificial targets both inside and outside the parameter space. We also present a test of our method on the well-studied SPB star KIC 7760680. We believe that this method is promising for carrying out detailed analyses of observations of SPB and $\gamma$ Dor stars and will provide complementary information to evolutionary models.

We investigate the impact of the Dark Energy Spectroscopic Instrument (DESI) 2024 data on dark energy scenarios. We thus analyze three typologies of models, the first in which the cosmic speed up is related to thermodynamics, the second associated with Taylor expansions of the barotropic factor, whereas the third based on \emph{ad hoc} dark energy parameterizations. In this respect, we perform Monte Carlo Markov chain analyses, adopting the Metropolis-Hastings algorithm, of 12 models. To do so, we first work at the background, inferring \emph{a posteriori} kinematic quantities associated with each model. Afterwards, we obtain early time predictions, computing departures on the growth evolution with respect to the model that better fits DESI data. We find that the best model to fit data \emph{is not} the Chevallier-Polarski-Linder (CPL) parametrization, but rather a more complicated log-corrected dark energy contribution. To check the goodness of our findings, we further directly fit the product, $r_d h_0$, concluding that $r_d h_0$ is anticorrelated with the mass. This treatment is worked out by removing a precise data point placed at $z=0.51$. Surprisingly, in this case the results again align with the $\Lambda$CDM model, \emph{indicating that the possible tension between the concordance paradigm and the CPL model can be severely alleviated}. We conclude that future data points will be essential to clarify whether dynamical dark energy is really in tension with the $\Lambda$CDM model.

$f(R)$ gravity is one of the serious alternatives of general relativity having large range of astronomical consequences. In this work, we study Big Bang Nucleosynthesis (BBN) in $f(R)$ gravity theory. We consider modification to gravity due to the existence of primordial black holes in the radiation era which introduce additional degrees of freedom known as scalarons. We calculate the light element abundances by using the BBN code PArthENoPE. It is found that for a range of scalaron mass $(2.2-3.5) \times 10^4$ eV, the abundance of lithium is lowered by $3-4$ times the value predicted by general relativistic BBN which is a level desired to address the cosmological lithium problem. For the above scalaron mass range helium abundance is within observed bound. However, deuterium abundance is found to be increased by $3-6$ times the observed primordial abundance which calls for high efficiency of stellar formation and evolution processes for destruction of the same. A novel relation between scalaron mass and black hole mass has been used to estimate that the above scalaron mass range corresponds to primordial black holes of sub-planetary mass ($\sim 10^{19}$ g) serving as one of the potential non-baryonic dark matter candidates. We infer Big Bang equivalence of power law $f(R)$ gravity with primordial black holes.

We employ Hubble data and Gaussian Processes in order to reconstruct the dynamical connection function in $f(Q)$ cosmology beyond the coincident gauge. In particular, there exist three branches of connections that satisfy the torsionless and curvatureless conditions, parameterized by a new dynamical function $\gamma$. We express the redshift dependence of $\gamma$ in terms of the $H(z)$ function and the $f(Q)$ form and parameters, and then we reconstruct it using 55 $H(z)$ observation data. Firstly, we investigate the case where ordinary conservation law holds, and we reconstruct the $f(Q)$ function, which is very well described by a quadratic correction on top of Symmetric Teleparallel Equivalent of General Relativity. Proceeding to the general case, we consider two of the most studied $f(Q)$ models of the literature, namely the square-root and the exponential one. In both cases we reconstruct $\gamma(z)$, and we show that according to AIC and BIC information criteria its inclusion is favoured compared to both $\Lambda$CDM paradigm, as well as to the same $f(Q)$ models under the coincident gauge. This feature acts as an indication that $f(Q)$ cosmology should be studied beyond the coincident gauge.

Future data provided by the \Euclid mission will allow us to better understand the cosmic history of the Universe. A metric of its performance is the figure-of-merit (FoM) of dark energy, usually estimated with Fisher forecasts. The expected FoM has previously been estimated taking into account the two main probes of \Euclid, namely the three-dimensional clustering of the spectroscopic galaxy sample, and the so-called 3$\times$2\,pt signal from the photometric sample (i.e., the weak lensing signal, the galaxy clustering, and their cross-correlation). So far, these two probes have been treated as independent. In this paper, we introduce a new observable given by the ratio of the (angular) two-point correlation function of galaxies from the two surveys. For identical (normalised) selection functions, this observable is unaffected by sampling noise, and its variance is solely controlled by Poisson noise. We present forecasts for \Euclid where this multi-tracer method is applied and is particularly relevant because the two surveys will cover the same area of the sky. This method allows for the exploitation of the combination of the spectroscopic and photometric samples. When the correlation between this new observable and the other probes is not taken into account, a significant gain is obtained in the FoM, as well as in the constraints on other cosmological parameters. The benefit is more pronounced for a commonly investigated modified gravity model, namely the $\gamma$ parametrisation of the growth factor. However, the correlation between the different probes is found to be significant and hence the actual gain is uncertain. We present various strategies for circumventing this issue and still extract useful information from the new observable.

We present an analysis of the whole 2018 outburst of the black hole X-ray binary MAXI J1820+070 with Insight-HXMT data. We focus our study on the temporal evolution of the parameters of the source. We employ two different models to fit the thermal spectrum of the disk: the Newtonian model DISKBB and the relativistic model NKBB. These two models provide different pictures of the source in the soft state. With DISKBB, we find that the inner edge of the disk is close to the innermost stable circular orbit of a fast-rotating black hole and the corona changes geometry from the hard to the soft state, probably becoming compact and close to the black hole. With NKBB, we find that the disk is truncated in the soft state and that the coronal geometry does not change significantly during the whole outburst. However, we find that the model with NKBB can predict an untruncated disk around a fast-rotating black hole if we assume that the inclination angle of the disk is around 30 degrees (instead of 60 degrees, which is the inclination angle of the jet and is usually adopted as the inclination angle of the disk in the literature) and we employ a high-density reflection model. In such a case, we measure a high value of the black hole spin parameter with observations in the soft state, in agreement with the high spin value found from the analysis of the reflection features and in disagreement with the low spin value found by previous continuum-fitting method measurements with the inclination angle of the disk set to the value of the inclination angle of the jet.

Amino acids are essential for the synthesis of protein. Amino acids contain both amine (R$-$NH$_{2}$) and carboxylic acid (R$-$COOH) functional groups, which help to understand the possible formation mechanism of life in the universe. Among the 20 types of amino acids, glycine (NH$_{2}$CH$_{2}$COOH) is known as the simplest non-essential amino acid. In the last 40 years, all surveys of NH$_{2}$CH$_{2}$COOH in the interstellar medium, especially in the star-formation regions, have failed at the millimeter and sub-millimeter wavelengths. We aimed to identify the possible precursors of NH$_{2}$CH$_{2}$COOH, because it is highly challenging to identify NH$_{2}$CH$_{2}$COOH in the interstellar medium. Many laboratory experiments have suggested that methylenimine (CH$_{2}$NH) plays a key role as a possible precursor of NH$_{2}$CH$_{2}$COOH in the star-formation regions via the Strecker synthesis reaction. After spectral analysis using the local thermodynamic equilibrium (LTE) model, we successfully identified the rotational emission lines of CH$_{2}$NH towards the hot molecular core G10.47+0.03 using the Atacama Compact Array (ACA). The estimated column density of CH$_{2}$NH towards G10.47+0.03 is (3.40$\pm$0.2)$\times$10$^{15}$ cm$^{-2}$ with a rotational temperature of 218.70$\pm$20 K, which is estimated from the rotational diagram. The fractional abundance of CH$_{2}$NH with respect to H$_{2}$ towards G10.47+0.03 is 2.61$\times$10$^{-8}$. We found that the derived abundance of CH$_{2}$NH agree fairly well with the existing two-phase warm-up chemical modelling abundance value of CH$_{2}$NH. We discuss the possible formation pathways of CH$_{2}$NH within the context of hot molecular cores, and we find that CH$_{2}$NH is likely mainly formed via neutral-neutral gas-phase reactions of CH$_{3}$ and NH radicals towards G10.47+0.03.

A near-infrared spectrum of nova V1716 Scorpii (PNV J17224490-4137160), a recent bright (V_max = 7.3 mag), Fermi-LAT detected gamma-ray source, was modeled using the photoionization code CLOUDY. Abundances were estimated for He, C, N, O, Si, Al, Mg, Fe, Ne, S, Ca, and P. Notably, P (a factor of 120) and N (a factor of 248) are highly overabundant. It was necessary to assume the ejecta consist of two components (with a cylindrical geometry): a dense component from which the bulk of the H, He, and neutral O~I and N emission arises and a more diffuse component from which most of the coronal lines arise. Some of the coronal lines are found to originate from both the dense and diffuse components. The mass of the ejecta, including neutral and ionized gas, is ~ 4.19e-4 solar masses. Our analysis indicates that in the case of V1716 Sco (which has a carbon-oxygen white dwarf), a fraction of 25% white dwarf matter rather than 50% is favored for the mixing between white dwarf and the accreted envelope before the outburst. This mixing ratio is like that found for Oxygen-Neon novae where a 25% mixing fraction is also indicated. Helium hydride -- the first molecule to form after the Big Bang -- may have formed in the ejecta of V1716 Sco based on photoionization modeling. This prediction suggests that novae may be potential formation sites of this important molecular ion.

Several stars show deep transits consistent with discs of roughly 1 Solar radius seen at moderate inclinations, likely surrounding planets on eccentric orbits. We show that this configuration arises naturally as a result of planet-planet scattering when the planets possess satellite systems. Planet-planet scattering explains the orbital eccentricities of the discs' host bodies, while the close encounters during scattering lead to the exchange of satellites between planets and/or their destabilisation. This leads to collisions between satellites and their tidal disruption close to the planet. Both of these events lead to large quantities of debris being produced, which in time will settle into a disc such as those observed. The mass of debris required is comparable to a Ceres-sized satellite. Through N-body simulations of planets with clones of the Galilean satellite system undergoing scattering, we show that 90 percent of planets undergoing scattering will possess debris from satellite destruction. Extrapolating to smaller numbers of satellites suggests that tens of percent of such planets should still possess circumplanetary debris discs. The debris trails arising from these events are often tilted at tens of degrees to the planetary orbit, consistent with the inclinations of the observed discs. Disruption of satellite systems during scattering thus simultaneously explains the existence of debris, the tilt of the discs, and the eccentricity of the planets they orbit.

In a previous paper using Gaia DR2 data, we demonstrated that the two closely situated open clusters Collinder 135 and UBC 7 might have formed together about 50 Myr ago. In this work, we performed star-by-star dynamical modelling of the evolution of the open clusters Collinder 135 and UBC 7 from their supposed initial state to their present-day state, reproducing observational distributions of members. Modelling of the Collinder 135 and UBC 7 dynamical evolution was done using the high-order parallel N-body code \phi-GPU with up-to-date stellar evolution. Membership and characteristics of the clusters were acquired based on Gaia DR3 data. The comparison of the present-day radial cumulative star count obtained from the N-body simulations with the current observational data gave us full consistency of the model with observational data, especially in the central 8 pc, where 80% of the stars reside. The proper motion velocity components obtained from the N-body simulations of the stars are also quite consistent with the observed distributions and error bars. These results show that our numerical modelling is able to reproduce the open clusters' current complex 6D observed phase-space distributions with a high level of confidence. Thus, the model demonstrates that the hypothesis of a common origin of Collinder 135 and UBC 7 complies with present-day observational data.

Broad Absorption Line Quasars (BALQSOs) displaying distinct blue-shifted broad absorption lines. These serve as invaluable probes for unraveling the intricate structure and evolution of quasars, shedding light on the profound influence exerted by supermassive black holes on galaxy formation. The proliferation of large-scale spectroscopic surveys such as LAMOST, SDSS, and DESI has exponentially expanded the repository of quasar spectra at our disposal. In this study, we present an innovative approach to streamline the identification of BALQSOs, leveraging the power of dimensionality reduction and machine learning algorithms. Our dataset is curated from the SDSS DR16, amalgamating quasar spectra with classification labels sourced from the DR16Q quasar catalog. We employ a diverse array of dimensionality reduction techniques, including Principal Component Analysis (PCA), t-Distributed Stochastic Neighbor Embedding (t-SNE), Locally Linear Embedding (LLE), and Isometric Mapping (ISOMAP), to distill the essence of the original spectral data. The resultant low-dimensional representations serve as inputs for a suite of machine learning classifiers, including XGBoost and Random Forest models. Through experimentation, we unveil PCA as the most effective dimensionality reduction methodology, adeptly navigating the intricate balance between dimensionality reduction and preservation of vital spectral information. Notably, the synergistic fusion of PCA with the XGBoost classifier emerges as the pinnacle of efficacy in the BALQSO classification endeavor, boasting impressive accuracy rates of 97.60% by 10-cross validation and 96.92% on the outer test sample. This study not only introduces a novel machine learning-based paradigm for quasar classification but also offers invaluable insights transferrable to a myriad of spectral classification challenges pervasive in the realm of astronomy.

We have developed a new computational method to explore astrophysical and heliophysical phenomena, especially those considerably influenced by non-thermal energetic particles. This novel approach considers the backreaction from these energetic particles by incorporating the non-thermal fluid pressure into Magnetohydrodynamics (MHD) equations. The pressure of the non-thermal fluid is evaluated from the energetic particle distribution evolved through Parker's transport equation, which is solved using stochastic differential equations. We implement this method in the HOW-MHD code (Seo \& Ryu 2023), which achieves 5th-order accuracy. We find that without spatial diffusion, the method accurately reproduces the Riemann solution in the hydrodynamic shock tube test when including the non-thermal pressure. Solving Parker's transport equation allows the determination of pressure terms for both relativistic and non-relativistic non-thermal fluids with adiabatic indices $\gamma_{\rm{NT}}=4/3$ and $\gamma_{\rm{NT}}=5/3$, respectively. The method also successfully replicates the Magnetohydrodynamic shock tube test with non-thermal pressure, successfully resolving the discontinuities within a few cells. Introducing spatial diffusion of non-thermal particles leads to marginal changes in the shock but smooths the contact discontinuity. Importantly, this method successfully simulates the energy spectrum of the non-thermal particles accelerated through shock, which includes feedback from the non-thermal population. These results demonstrate that this method is very powerful for studying particle acceleration when a significant portion of the plasma energy is taken by energetic particles.

This first VELOCE data release comprises 18,225 high-precision RV measurements of 258 bona fide classical Cepheids on both hemispheres collected mainly between 2010 and 2022, alongside 1161 additional observations of 164 other stars. The median per-observation RV uncertainty is 0.037 km/s, and some reach 0.002 km/s. Non-variable standard stars characterize RV zero-point stability and provide a base for future cross-calibrations. We determined zero-point differences between VELOCE and 31 literature data sets using template fitting and measured linear period changes of 146 Cepheids. Seventy six spectroscopic binary Cepheids and 14 candidates are identified using VELOCE data alone and are investigated in detail in a companion paper (VELOCE II). Several new insights into Cepheid pulsations were obtained, including: a) the most detailed description of the Hertzsprung progression by RVs; b) the identification of double-peaked bumps in the RV curve; c) clear evidence that virtually all Cepheids feature spectroscopic variability signals that lead to modulated RV variability. We identified 36 such stars, of which 4 also exhibit orbital motion. Linear radius variations depend strongly on pulsation period and a steep increase in slope of the $\Delta$R/p versus logP-relation is found near 10d, challenging the existence of a tight relation between Baade-Wesselink projection factors and pulsation periods. We investigated the accuracy of RV time series measurements, v$_\gamma$, and RV amplitudes published in Gaia's DR3 and determined an average offset of 0.65 \pm 0.11 km/s relative to VELOCE. We recommend adopting a single set of template correlation parameters for distinct classes of large-amplitude variable stars to avoid systematic offsets in v$_\gamma$ among stars belonging to the same class. Peak-to-peak amplitudes of Gaia RVs exhibit significant (16%) dispersion compared to VELOCE. [abridged]

We present a comprehensive multiwavelength investigation into flares and activity in nearby M~dwarf stars. We leverage the most extensive contemporaneous dataset obtained through the Transiting Exoplanet Sky Survey (TESS), Kepler/K2, the Neil Gehrels Swift Observatory (\textit{Swift}), and the Hubble Space Telescope (HST), spanning the optical and near-ultraviolet (NUV) regimes. In total, we observed 213 NUV flares on 24 nearby M dwarfs, with $\sim$27\% of them having detected optical counterparts, and found that all optical flares had NUV counterparts. We explore NUV/optical energy fractionation in M dwarf flares. Our findings reveal a slight decrease in the ratio of optical to NUV energies with increasing NUV energies, a trend in agreement with prior investigations on G-K stars' flares at higher energies. Our analysis yields an average NUV fraction of flaring time for M0-M3 dwarfs of 2.1\%, while for M4-M6 dwarfs, it is 5\%. We present an empirical relationship between NUV and optical flare energies and compare to predictions from radiative-hydrodynamic and blackbody models. We conducted a comparison of the flare frequency distribution (FFDs) of NUV and optical flares, revealing the FFDs of both NUV and optical flares exhibit comparable slopes across all spectral subtypes. NUV flares on stars affect the atmospheric chemistry, the radiation environment, and the overall potential to sustain life on any exoplanets they host. We find that early and mid-M dwarfs (M0-M5) have the potential to generate NUV flares capable of initiating abiogenesis.

Research has shown that many young and intermediate-age clusters (younger than $\sim$2 Gyr) have extended main sequences and main-sequence turnoffs (eMSTOs), which cannot be adequately described by a single isochrone. The reason for the extended main sequences is now known, with the most probable cause being the fast rotation of stars. However, a significant fraction of slowly rotating stars form a younger stellar population than their fast-rotating counterparts, leading to speculation that they have undergone thorough rotational mixing processes internally. One speculation is that a considerable number of slowly rotating stars reside in close binary systems, where tidal forces from companion stars are the cause of their rotational deceleration. In this work, we report a relatively old open star cluster in the Milky Way, NGC 2423 ($\sim$1 Gyrs old), which exhibits an apparent eMSTO. As anticipated, many characteristics of NGC 2423 indicate that its eMSTO is driven by stellar rotations. Our calculations indicate that if slowly rotating stars commonly have a close companion star, they should exhibit significant differences in radial velocities observationally, and binary systems that can be tidally locked within the age of NGC 2423 should have a mass ratio close to 1. However, none of these predictions align with our observations. Interestingly, among the only two equal-mass binary systems in the observed region for which spectroscopic data could be obtained, we discovered that one of them is a tidally locked binary system. This further suggests the validity of our numerical simulation results.

We present the rest-frame ultraviolet-optical spectral properties of 65 broad absorption line (BAL) quasars from the Gemini Near Infrared Spectrograph-Distant Quasar Survey (GNIRS-DQS). These properties are compared with those of 195 non-BAL quasars from GNIRS-DQS in order to identify the drivers for the appearance of BALs in quasar spectra. In particular, we compare equivalent widths and velocity widths, as well as velocity offsets from systemic redshifts, of principal emission lines. In spite of the differences between their rest-frame ultraviolet spectra, we find that luminous BAL quasars are generally indistinguishable from their non-BAL counterparts in the rest-frame optical band at redshifts $1.55 \lesssim z \lesssim 3.50$. We do not find any correlation between BAL trough properties and the H$\beta$-based supermassive black hole masses and normalized accretion rates in our sample. Considering the Sloan Digital Sky Survey quasar sample, which includes the GNIRS-DQS sample, we find that a monochromatic luminosity at rest-frame 2500 A of $\gtrsim 10^{45}$ erg s$^{-1}$ is a necessary condition for launching BAL outflows in quasars. We compare our findings with other BAL quasar samples and discuss the roles that accretion rate and orientation play in the appearance of BAL troughs in quasar spectra.

Year 1 results of the Legacy Survey of Space and Time (LSST) will provide tighter constraints on small-scale cosmology, beyond the validity of linear perturbation theory. This heightens the demand for a computationally affordable prescription that can accurately capture nonlinearities in beyond-$\Lambda$CDM models. The COmoving Lagrangian Acceleration (COLA) method, a cost-effective \textit{N}-body technique, has been proposed as a viable alternative to high-resolution \textit{N}-body simulations for training emulators of the nonlinear matter power spectrum. In this study, we evaluate this approach by employing COLA emulators to conduct a cosmic shear analysis with LSST-Y1 simulated data across three different nonlinear scale cuts. We use the $w$CDM model, for which the \textsc{EuclidEmulator2} (\textsc{ee2}) exists as a benchmark, having been trained with high-resolution \textit{N}-body simulations. We primarily utilize COLA simulations with mass resolution $M_{\rm part}\approx 8 \times 10^{10} ~h^{-1} M_{\odot}$ and force resolution $\ell_{\rm force}=0.5 ~h^{-1}$Mpc, though we also test refined settings with $M_{\rm part}\approx 1 \times 10^{10} ~h^{-1}M_{\odot}$ and force resolution $\ell_{\rm force}=0.17 ~h^{-1}$Mpc. We find the performance of the COLA emulators is sensitive to the placement of high-resolution \textit{N}-body reference samples inside the prior, which only ensure agreement in their local vicinity. However, the COLA emulators pass stringent criteria in goodness-of-fit and parameter bias throughout the prior, when $\Lambda$CDM predictions of \textsc{ee2} are computed alongside every COLA emulator prediction, suggesting a promising approach for extended models.

Detection of radio emission from Jupiter was identified quickly as being due to its planetary-scale magnetic field. Subsequent spacecraft investigations have revealed that many of the planets, and even some moons, either have or have had large-scale magnetic fields. In the case of the Earth, Jupiter, Saturn, Uranus, and Neptune, the their magnetic fields are generated by dynamo processes within these planets, and an interaction between the solar wind and their magnetic fields generates intense radio emission via the electron cyclotron maser instability. In the case of Jupiter, its magnetic field interacts with the moon Io to result in radio emission as well. Extrasolar planets reasonably may be expected to generate large-scale magnetic fields and to sustain an electron cyclotron maser instability. Not only may these radio emissions be a means for discovering extrasolar planets, because magnetic fields are tied to the properties of planetary interiors, radio emissions may be a remote sensing means of constraining extrasolar planetary properties that will be otherwise difficult to access. In the case of terrestrial planets, the presence or absence of a magnetic field may be an indicator for habitability. Since the first edition of the Handbook, there have been a number of advances, albeit there remain no unambigous detection of radio emission from extrasolar planets. New ground-based telescopes and new possibilities for space-based telescopes provide promise for the future.

Measurements of the luminosity distance of propagating gravitational waves can provide invaluable information on the geometry and content of our Universe. Due to the clustering of cosmic structures, in realistic situations we need to average the luminosity distance of events coming from patches inside a volume. In this work we evaluate, in a gauge-invariant and fully-relativistic treatment, the impact of cosmological perturbations on such averaging process. We find that clustering, lensing and peculiar velocity effects impact estimates for future detectors such as Einstein Telescope, Cosmic Explorer, the Big Bang Observer and DECIGO. The signal-to-noise ratio of the angular power spectrum of the average luminosity distance over all the redshift bins is 17 in the case of binary black holes detected by Einstein Telescope and Cosmic Explorer. We also provide fitting formulas for the corrections to the average luminosity distance due to general relativistic effects.

A population of non-stellar black holes ($\gtrsim$100 M$_{\odot}$) has been long predicted to wander the Milky Way. We aim to characterize this population by using the L-Galaxies semi-analytical model applied on top of the high resolution Millennium-II merger trees. Our results predict $\sim$10 wandering black holes with masses $\sim$2 $\times$ 10$^{3}$ M$_{\odot}$ in a typical $z$ = 0 Milky Way galaxy, accounting for $\sim$2$\%$ of the total non-stellar black hole mass budget of the galaxy. We find that the locations of these wanderers correlate with their formation scenario. While the ones concentrated at $\lesssim$1 kpc from the galactic nucleus on the disk come from past galactic mergers, the ones formed as a consequence of ejections due to gravitational recoils or the disruption of satellite galaxies are typically located at $\gtrsim$100 kpc. Such small and large distances might explain the absence of strong observational evidence for wandering black holes in the Milky Way. Our results also indicate that $\sim$67$\%$ of the wandering population is conformed by the leftovers of black hole seeds that had little to no growth since their formation. We find that wandering black holes that are leftover seeds become wanderers at an earlier time with respect to grown seeds, and also come from more metal-poor galaxies. Finally, we show that the number of wandering black holes in a Milky Way-type galaxy depends on the seeding efficiency.

General circulation models of gas giant exoplanets predict equatorial jets that drive inhomogeneities across the planetary atmosphere. We studied the transmission spectrum of the hot Jupiter WASP-127b during one transit in the K band with CRIRES+. Telluric and stellar signals were removed from the data using SYSREM and the planetary signal was investigated using the cross-correlation (CCF) technique. After detecting a spectral signal indicative of atmospheric inhomogeneities, we employed a Bayesian retrieval framework with a 2D modelling approach tailored to address this scenario. We detected strong signals of H$_2$O and CO, which exhibited not one but two distinct CCF peaks. The double peaked signal can be explained by a supersonic equatorial jet and muted signals at the poles, with the two peaks representing the signals from the planet's morning and evening terminators, respectively. We calculated an equatorial jet velocity of $7.7\pm0.2$km/s from our retrieved overall equatorial velocity and the planet's tidally locked rotation, and derive distinct atmospheric properties for the two terminators as well as the polar region. The evening terminator is found to be hotter than the morning terminator by $175^{+116}_{-133}$K and the muted signals from the poles can be explained by significantly lower temperatures or a high cloud deck. Our retrieval yields a solar C/O ratio and metallicity and challenges previous studies of WASP-127b's atmosphere. The presence of a double peaked signal highlights the importance of accounting for planetary 3D structure during interpretation of atmospheric signals. The measured supersonic jet velocity and the lack of signal from the polar regions, representing a detection of latitudinal inhomogeneity in a spatially unresolved target, showcases the power of high-resolution transmission spectroscopy for the characterization of global circulation in exoplanets.

Axion and axion-like particles (ALPs) are a prominent candidate for physics beyond the Standard Model, and can play an important role in cosmology, serving as dark matter or dark energy, or both, drawing motivation in part from the string theory axiverse. Axion-like particles (ALPs) can also arise as composite degrees of freedom following chiral symmetry breaking in a dark confining gauge theory, analogous to the Standard Model (SM) pion. A dark sector with arbitrary $N_f$ flavors of dark quarks leads to $N_f^2-1$ axion-like states, effectively a field theory axiverse (or '$\pi$-axiverse'). A portal to the visible sector can be achieved through the standard kinetic mixing between the dark photon and SM photon, generating millicharges for the dark quarks and consequently couplings, both parity-even and parity-odd, between the SM and the dark pions. This scenario has been studied for the $N_f=2$ case and more recently for a dark Standard Model with $N_f=6$. In this work, we study the spectrum of this field theory axiverse for an arbitrary number of flavors, and apply this to the example $N_f=10$. We calculate the couplings to the SM photon analogous to the conventional axion-photon coupling, including the $N_f$ and $N_c$ dependence, and compute the present and future constraints on the $N_f=10$, $N_c=3$, $\pi$-axiverse. We elucidate the accompanying 'bary-verse' of superheavy dark baryons, namely an ensemble of charged and neutral dark baryons with a mass set by the dark pion decay constant.

A semiclassical investigation of the electromagnetic radiation emitted by a charged particle in a radially freely falling motion in Schwarzschild spacetime is carried out. We use quantum field theory at tree level to obtain the one-particle-emission amplitudes. We analyze and compare the energy spectrum and total energy released, which are calculated from these amplitudes, for particles with varying initial positions and for particles originating from infinity with varying kinetic energy. We also compare the results with those due to a falling charged "string" extended in the radial direction.

We explore compactifications of the form of three tori and one circle in the framework of 11D supergravity. By imposing suitable gauge conditions and boundary conditions, we find that the four-dimensional FRW universe emerges as a solution representing cosmological D3-branes in the eleven-dimensional bulk. These specific compactification methods can produce cosmological inflation that aligns with the observational constraints set by the 2021 BICEP/Keck and Planck 2018 results. In the cosmological inflation models we construct, the inflaton can be interpreted as the conformal vibrations of extra dimensions with a size around 10^5 times the reduced Planck length. Additionally, we offer the theoretical predictions for the mass of the inflaton, and the tree-level Newton's gravity law between two massive point particles surrounded by a spherically symmetric distribution of the inflaton, which can reproduce the Tully-Fisher relation and explain the flat rotation curves of galaxies.

We have developed a set of four fully coupled Boltzmann equations to precisely determine the relic density and temperature of dark matter by including three distinct sectors: dark matter, light scalar, and standard model sectors. The intricacies of heat transfer between DM and the SM sector through a light scalar particle are explored, inspired by stringent experimental constraints on the scalar-Higgs mixing angle and the DM-scalar coupling. Three distinct sectors emerge prior to DM freeze-out, requiring fully coupled Boltzmann equations to accurately compute relic density. Investigation of forbidden, resonance, and secluded DM scenarios demonstrates significant deviations between established methods and the novel approach with fully coupled Boltzmann equations. Despite increased computational demands, this emphasizes the need for improved precision in relic density calculations, underlining the importance of incorporating these equations in comprehensive analyses.

The gravitational-wave signal GW170817 is a result of a binary neutron star coalescence event. The observations of electromagnetic counterparts suggest that the event didn't led to the prompt formation of a black-hole. In this work, we first classify the GW170817 LIGO-Virgo data sample into prompt collapse to a black-hole using the $q$-dependent threshold mass fits and then remove these cases from the data sample. We find that the cases without a prompt black-hole formation do not support radii $ <$ 10 km unlike the LIGO-Virgo data sample. This is consistent with the maximum mass constraint, based on the binary pulsar J0348+0432, imposed LIGO-Virgo data sample. Additionally, we find that the cases without the prompt collapse to a black-hole improve the uncertainty range of neutron star radii from 3.3 km to 2.6 km for the data sample without the mass constraint and from 2.8 km to 2.5 km for the data sample with the mass constraint, implying improved constraints on the neutron star radii and hence the equation-of-state.

With the rise of satellites and mankind's growing dependence on technology, there is an increasing awareness of space weather phenomena related to high-energy particles. Shock waves driven by coronal mass ejections (CMEs) and corotating interaction regions (CIRs) occasionally act as potent particle accelerators, generating hazardous solar energetic particles (SEPs) that pose risks to satellite electronics and astronauts. Numerical simulation tools capable of modelling and predicting large SEP events are thus highly demanded. We introduce the new Icarus$+$PARADISE model as an advancement of the previous EUHFORIA$+$PARADISE model. Icarus, based on the MPI-AMRVAC framework, is a three-dimensional magnetohydrodynamic code that models solar wind configurations from 0.1 au onwards, encompassing transient structures like CMEs or CIRs. Differing from EUHFORIA's uniform-only grid, Icarus incorporates solution adaptive mesh refinement (AMR) and grid stretching. The particle transport code PARADISE propagates energetic particles as test particles through these solar wind configurations by solving the focused transport equation in a stochastic manner. We validate our new model by reproducing EUHFORIA+PARADISE results. This is done by modelling the acceleration and transport of energetic particles in a synthetic solar wind configuration containing an embedded CIR. Subsequently, we illustrate how the simulation results vary with grid resolution by employing different levels of AMR. The resulting intensity profiles illustrate increased particle acceleration with higher levels of AMR in the shock region, better capturing the effects of the shock.

In this paper, we present numerical and experimental results on helicity oscillations in a liquid-metal Rayleigh-B\'enard (RB) convection cell, with an aspect ratio of 0.5. We find that helicity oscillations occur during transitions of flow states that are characterised by significant changes in the Reynolds number. Moreover, we also observe helicity oscillations at flow conditions where the temporal gradient of the change in the Reynolds number is significantly smaller than that of the helicity. Notably, the helicity oscillations observed during the transient double-roll state exhibit characteristics remarkably similar to those associated with the Tayler Instability.

We adopt general relativistic ray-tracing (GRRT) schemes to study images of Kerr-MOG black holes surrounded by geometrically thick magnetized equilibrium tori, which belong to steady-state solutions of thick accretion disks within the framework of general relativistic magnetohydrodynamics (GRMHD). The black hole possesses an extra dimensionless MOG parameter described its deviation from usual Kerr one. Our results show that the presence of the MOG parameter leads to smaller disks in size, but enhances the total flux density and peak brightness in their images. Combining with observation data of black hole M87* from the Event Horizon Telescope (EHT), we make a constraint on parameters of the Kerr-MOG black hole and find that the presence of the MOG parameter broadens the allowable range of black hole spin.

We propose a novel method (floZ), based on normalizing flows, for estimating the Bayesian evidence (and its numerical uncertainty) from a set of samples drawn from the unnormalized posterior distribution. We validate it on distributions whose evidence is known analytically, up to 15 parameter space dimensions, and compare with two state-of-the-art techniques for estimating the evidence: nested sampling (which computes the evidence as its main target) and a k-nearest-neighbors technique that produces evidence estimates from posterior samples. Provided representative samples from the target posterior are available, our method is more robust to posterior distributions with sharp features, especially in higher dimensions. It has wide applicability, e.g., to estimate the evidence from variational inference, Markov-chain Monte Carlo samples, or any other method that delivers samples from the unnormalized posterior density.

In the context of $f(R,T)$ gravity and other modified theories of gravity, the knowledge of the first order variation of the trace $T$ of the energy-momentum tensor with respect to the metric is essential for an accurate characterization of the gravitational field. In this paper, by considering a paradigmatic example of a perfect fluid whose dynamics is described by a pure k-essence matter Lagrangian in $f(R,T)=R+\mathcal F(T)$ gravity, we show that the first order variation of the trace of the energy-momentum tensor cannot in general be determined from the proper density, proper pressure and 4-velocity of the fluid alone, and that the sound speed of the fluid can directly influence the dynamics of gravity. We also confirm that the second variation of the matter Lagrangian with respect to the metric should not in general be neglected. These results can be particularly relevant for cosmological studies of $f(R,T)$ gravity in which some of the material content of the Universe is modeled as a perfect fluid.