Papers for Wednesday, Jun 22 2022

We present a mock image catalogue of ~100,000 MUV=-22.5 to -19.6 mag galaxies at z=7-12 from the BlueTides cosmological simulation. We create mock images of each galaxy with the James Webb (JWST), Hubble, Roman, and Euclid Space Telescopes, as well as Subaru, and VISTA, with a range of near- and mid-infrared filters. We perform photometry on the mock images to estimate the success of these instruments for detecting high-z galaxies. We predict that JWST will have unprecedented power in detecting high-z galaxies, with a 95% completeness limit at least 2.5 magnitudes fainter than VISTA and Subaru, 1.1 magnitudes fainter than Hubble and 0.9 magnitudes fainter than Roman, for the same wavelength and exposure time. Focusing on JWST, we consider a range of exposure times and filters, and find that the NIRCam F356W and F277W filters will detect the faintest galaxies, with 95% completeness at m=27.4 mag in 10ks exposures. We also predict the number of high-z galaxies that will be discovered by upcoming JWST imaging surveys. We predict that the COSMOS-Web survey will detect ~1000 MUV<-20.1 mag galaxies at 6.5<z<7.5, by virtue of its large survey area. JADES-Medium will detect almost 100% of MUV<-20 mag galaxies at z<8.5 due to its significant depth, however with its smaller survey area it will detect only ~100 of these galaxies at 6.5<z<7.5. Cosmic variance results in a large range in the number of predicted galaxies each survey will detect, which is more evident in smaller surveys such as CEERS and the PEARLS NEP and GOODS-S fields.

We explore the nature of low surface brightness galaxies (LSBGs) in the hydrodynamic cosmological simulation TNG100 of the IllustrisTNG project, selecting a sample of LSBGs ($r$-band effective surface brightness $\mu_r > 22.0$ mag arcsec$^{-2}$) at $z=0$ over a wide range of stellar masses ($M_{\ast} = 10^{9}$-$10^{12}$ M$_\odot$). We find LSBGs of all stellar masses, although they are particularly prevalent at $M_{\ast} < 10^{10}$ M$_\odot$. We show that the specific star formation rates of LSBGs are not significantly different from those of high surface brightness galaxies (HSBGs) but, as a population, LSBGs are systematically less massive and more extended than HSBGs, and tend to display late-type morphologies according to a kinematic criterion. At fixed stellar mass, we find that haloes hosting LSBGs are systematically more massive and have a higher baryonic fraction than those hosting HSBGs. We find that LSBGs have higher stellar specific angular momentum and halo spin parameter values compared to HSBGs, as suggested by previous works. We track the evolution of these quantities back in time, finding that the spin parameters of the haloes hosting LSBGs and HSBGs exhibit a clear bifurcation at $z \sim 2$, which causes a similar separation in the evolutionary tracks of other properties such as galactic angular momentum and effective radius, ultimately resulting in the values observed at $z =$ 0. The higher values of specific stellar angular momentum and halo spin in LSBGs seem to be responsible for their extended nature, preventing material from collapsing into the central regions of the galaxies, also causing LSBGs to host less massive black holes at their centres.

We present a detailed analysis of short GRB 201221D lying at redshift $z= 1.045$. We analyse the high-energy data of the burst and compare it with the sample of short gamma-ray bursts (SGRBs). The prompt emission characteristics are typical of those seen in the case of other SGRBs except for the peak energy ($E_{\rm p}$), which lies at the softer end (generally observed in the case of long bursts). We estimate the host galaxy properties by utilising the {\sc Python}-based software {\sc Prospector} to fit the spectral energy distribution of the host. The burst lies at a high redshift relative to the SGRB sample with a median redshift of $z=0.47$. We compare the burst characteristics with other SGRBs with known redshifts along with GRB 200826A (SGRB originated from a collapsar). A careful examination of the characteristics of SGRBs at different redshifts reveals that some of the SGRBs lying at high redshifts have properties similar to long GRBs indicating they might have originated from collapsars. Further study of these GRBs can help to explore the broad picture of progenitor systems of SGRBs.

Fast Radio Bursts (FRBs), bright millisecond-duration radio transients are happening thousands of times per day. FRBs' astrophysical mechanisms are still puzzling. Bustling Universe Radio Survey Telescope for Taiwan (BURSTT) is optimized to discover and localize a large sample of bright and nearby FRBs. BURSTT will have a large field-of-view (FoV) of $\sim$10$^{4}$ deg$^{2}$ for monitoring the whole visible sky all the time, a 400 MHz effective bandwidth between 300-800 MHz, and the sub-arcsecond localization capability with several outrigger stations hundreds to thousands of km away. Initially, BURSTT will equip with 256 antennas, which we will test different designs and improve the system performance. Through the scalable features, BURSTT could equip with more antennas and eventually optimized designs. We expect that BURSTT initially would detect and localize $\sim$100 bright ($\geq$100 Jy ms) and nearby FRBs per year to sub-arcsecond precision. Besides, the large FoV yields monitoring FRBs with high cadence, which is crucial to understanding the repetition of FRBs. Multi-wavelength/multi-messenger observations of the BURSTT localized bright samples would be the key to understanding the nature as well as the local environment of FRBs.

Current and future high contrast imaging instruments aim to detect exoplanets at closer orbital separations, lower masses, and/or older ages than their predecessors. However, continually evolving speckles in the coronagraphic science image limit contrasts of state-of-the-art ground-based exoplanet imaging instruments. For ground-based adaptive optics (AO) instruments it remains challenging for most speckle suppression techniques to attenuate both the dynamic atmospheric as well as quasi-static instrumental speckles on-sky. We have proposed a focal plane wavefront sensing and control algorithm to address this challenge, called the Fast Atmospheric Self-coherent camera (SCC) Technique (FAST), which in theory enables the SCC to operate down to millisecond timescales even when only a few photons are detected per speckle. Here we present the first experimental results of FAST on the Santa Cruz Extreme AO Laboratory (SEAL) testbed. In particular, we illustrate the benefit of ``second stage'' AO-based focal plane wavefront control, demonstrating up to 5x contrast improvement with FAST closed-loop compensation of evolving residual atmospheric turbulence -- both for low and high order spatial modes -- down to 20 millisecond-timescales.

Through a homogeneous analysis of spectroscopic literature data of red giant stars, we determine the radial metallicity profiles of 30 dwarf galaxies in the Local Group. We explore correlations between the calculated metallicity gradients and stellar mass, star formation history and environment, delivering the largest compilation to date of this type. The dwarf galaxies in our sample mostly show metallicity profiles decreasing with radius, with some exhibiting rather steep profiles. The derived metallicity gradients as a function of the half-light radius, $\nabla_{\rm [Fe/H]} (R/R_e)$, show no statistical differences when compared with the galaxies' morphological type, nor with their distance from the Milky Way or M31. No correlations are found with either stellar mass or star formation timescales. In particular, we do not find the linear relationship between $\nabla_{\rm [Fe/H]} (R/R_e)$ and the galaxies' median age $t_{50}$, as instead shown in the literature for a set of simulated systems. The presence of high angular momentum in some of our galaxies does not seem to have an impact on the gradient values. The strongest gradients in our sample are observed in systems that are likely to have experienced a past merger event. By excluding them, the analysed dwarf galaxies show mild gradients ($\sim -0.1$ dex $R_e^{-1}$) with little scatter between them, regardless of their stellar mass, dynamical state, and star formation history. These results are in good agreement with different sets of simulations presented in the literature and analysed using the same method as for the observed sample. The interplay between the multitude of factors that could drive the formation of metallicity gradients in dwarf galaxies likely combine in complex ways to produce in general comparable values.

We present a catalog of 1.4 million photometrically-selected quasar candidates in the southern hemisphere over the $\sim 5000\,{\rm deg^2}$ Dark Energy Survey (DES) wide survey area. We combine optical photometry from the DES second data release (DR2) with available near-infrared (NIR) and the all-sky unWISE mid-infrared photometry in the selection. We build models of quasars, galaxies, and stars with multivariate Skew-t distributions in the multi-dimensional space of relative fluxes as functions of redshift (or color for stars) and magnitude. Our selection algorithm assigns probabilities for quasars, galaxies, and stars, and simultaneously calculates photometric redshifts (photo-$z$) for quasar and galaxy candidates. Benchmarking on spectroscopically confirmed objects, we successfully classify (with photometry) 94.7% of quasars, 99.3% of galaxies, and 96.3% of stars when all IR bands (NIR $YJHK$ and WISE $W1W2$) are available. The classification and photo-$z$ regression success rates decrease when fewer bands are available. Our quasar (galaxy) photo-$z$ quality, defined as the fraction of objects with the difference between the photo-$z$ $z_p$ and the spectroscopic redshift $z_s$, $|\Delta z| = |z_s - z_p|/(1 + z_s)\le 0.1$, is 92.2% (98.1%) when all IR bands are available, decreasing to 72.2% (90.0%) using optical DES data only. Our photometric quasar catalog achieves estimated completeness of 89% and purity of 79% at $r<21.5$ (0.68 million quasar candidates), with reduced completeness and purity at $21.5<r\lesssim 24$. Among the 1.4 million quasar candidates, 87,857 have existing spectra and 84,978 (96.7%) of them are spectroscopically confirmed quasars. Finally, we provide quasar, galaxy, and star probabilities for all (0.69 billion) photometric sources in the DES DR2 coadded photometric catalog.

Using the QUIJOTE line survey in the 32.0-50.4 GHz range, we report the discovery of the molecule CH2CCHC4H towards the prestellar cold core TMC-1 in the Taurus region. We also present a rigorous detection of CH2CCHC3N, along with its detailed analysis. We identified a total of twenty rotational transitions for each one of these molecules. The rotational quantum numbers range from Ju=17 up to 24 and Ka<=3. The column density for CH2CCHC4H is N=(2.2+/-0.2)x 1E12 cm-2, while for CH2CCHC3N, we derived N=(1.2+/-0.15) x 1E11 cm-2. The rotational temperature is 9.0+/-0.5 K for both species. The abundance ratio between CH2CCHC4H and CH2CCHC3N is 18+/-4. We also compared the column densities of these species with those of their isomers CH3C6H and CH3C5N, derived from their J=20-19 up to J=30-29 rotational transitions observed with the QUIJOTE line survey. The observed abundances for all these species are reasonably well explained by state-of-the-art chemical models of TMC-1. The observed astronomical frequencies were merged with laboratory frequencies from the literature to derive improved spectroscopic parameters.

This work presents the medium-resolution ($R \sim 1,500$) spectroscopic follow-up of 522 low-metallicity star candidates from the Southern Photometric Local Universe Survey (S-PLUS). The objects were selected from narrow-band photometry, taking advantage of the metallicity-sensitive S-PLUS colors. The follow-up observations were conducted with the Blanco and Gemini South telescopes, using the COSMOS and GMOS spectrographs, respectively. The stellar atmospheric parameters (T$_{\rm eff}$, log$g$, and [Fe/H]), as well as carbon and $\alpha$-element abundances, were calculated for the program stars in order to assess the efficacy of the color selection. Results show that $92^{+2}_{-3}\%$ of the observed stars have [Fe/H]$\leq -1.0$, $83^{+3}_{-3}\%$ have [Fe/H]$\leq -2.0$, and $15^{+3}_{-3}\%$ have [Fe/H]$\leq -3.0$, including two ultra metal-poor stars ([Fe/H]$\leq -4.0$). The 80th percentile for the metallicity cumulative distribution function of the observed sample is [Fe/H]$= -2.04$. The sample also includes 68 Carbon-Enhanced Metal-Poor (CEMP) stars. Based on the calculated metallicities, further S-PLUS, color cuts are proposed, which can increase the fractions of stars with [Fe/H]$\leq -1.0$ and $\leq -2.0$ to $98\%$ and $88\%$, respectively. Such high success rates enable targeted high-resolution spectroscopic follow-up efforts, as well as provide selection criteria for fiber-fed multiplex spectroscopic surveys.

Arcus is a proposed soft x-ray grating spectrometer Explorer. It aims to explore cosmic feedback by mapping hot gases within and between galaxies and galaxy clusters and characterizing jets and winds from supermassive black holes, and to investigate the dynamics of protoplanetary discs and stellar accretion. Arcus features 12 m-focal-length grazing-incidence silicon pore optics (SPO) developed for the Athena mission. Critical-angle transmission (CAT) gratings efficiently disperse high diffraction orders onto CCDs. We report new and improved x-ray performance results for Arcus-like CAT gratings, including record resolving power for two co-aligned CAT gratings. Multiple Arcus prototype grating facets were illuminated by an SPO at the PANTER facility. The facets consist of $32\times32.5$ mm$^2$ patterned silicon membranes, bonded to metal frames. The bonding angle is adjusted according to the measured average tilt angle of the grating bars in the membrane. Two simultaneously illuminated facets show minor broadening of the Al-K$_{\alpha}$ doublet in 18$^{\rm th}$ and 21$^{\rm st}$ orders with a best fit record effective resolving power of $R_G \approx 1.3^{+\infty}_{-0.5}\times10^4$ ($3\sigma$), about 3-4 times the Arcus requirement. We measured the diffraction efficiency of quasi-fully illuminated gratings at O-K wavelengths in orders 4-7 in an Arcus-like configuration and compare results with synchrotron spot measurements. After corrections for geometrical effects and bremsstrahlung continuum we find agreement between full and spot illumination at the two different facilities, as well as with the models used for Arcus effective area predictions. We find that these flight-like gratings meet diffraction efficiency and greatly exceed resolving power Arcus requirements.

This paper is the second in a series aimed at examining the gaseous environments of z$\le$0.3 quasars and ultraluminous infrared galaxies (ULIRGs) as a function of AGN/host galaxy properties across the merger sequence. This second paper focuses on the Ly$\alpha$ emission and O VI and N V absorption features, tracers of highly ionized gas outflows, in ULIRGs observed with HST/COS. Ly$\alpha$ emission is detected in 15 out of 19 ULIRGs, and 12 of the 14 clear Ly$\alpha$ detections show emission with blueshifted velocity centroids and/or wings. The equivalent widths of the Ly$\alpha$ emission increase with increasing AGN luminosities and AGN bolometric fractions. The blueshifts of the Ly$\alpha$ emission correlate positively with those of the [O III] emission, where the latter traces the ionized gas outflows. The Ly$\alpha$ escape fractions tend to be slightly larger in objects with stronger AGN and larger outflow velocities, but they do not correlate with nebular line reddening. Among the 12 ULIRGs with good continuum signal-to-noise ratios, O VI and/or N V absorption features are robustly detected in 6 of them, all of which are blueshifted, indicative of outflows. In the combined ULIRG $+$ quasar sample, the outflows are more frequently detected in the X-ray weak or absorbed sources. The absorption equivalent widths, velocities and velocity dispersions of the outflows are also higher in the X-ray weak sources. No other strong correlations are visible between the properties of the outflows and those of the AGN or host galaxies.

We present a technique to identify spectrophotometrically variable L7-T3 brown dwarfs with single-epoch, low-resolution, near-infrared SpeX spectra. We calculated spectral indices on known variable brown dwarfs and used them to select 11 index-index parameter spaces where known variables can be distinguished from the rest of the general population of brown dwarfs. We find 62 candidate variables, 12 of which show significant variability amplitude in independent photometric monitoring surveys. This technique constitutes the first formal method to identify a time-dependent effect such as variability from peculiarities in their integrated light spectra. This technique will be a useful tool to prioritize targets for future photometric and spectroscopic monitoring in the era of JWST and 30m-class telescopes.

The LCO Outbursting Objects Key (LOOK) Project uses the telescopes of the Las Cumbres Observatory (LCO) Network to: (1) to systematically monitor a sample of Dynamically New Comets over the whole sky, and (2) use alerts from existing sky surveys to rapidly respond to and characterize detected outburst activity in all small bodies. The data gathered on outbursts helps to characterize each outburst's evolution with time, assess the frequency and magnitude distribution of outbursts in general, and contributes to the understanding of outburst processes and volatile distribution in the Solar System. The LOOK Project exploits the synergy between current and future wide-field surveys such as ZTF, PanSTARRS, and LSST as well as rapid-response telescope networks such as LCO, and serves as an excellent testbed for what will be needed the much larger number of objects coming from Rubin Observatory. We will describe the LOOK Project goals, the planning and target selection (including the use of NEOexchange as a Target and Observation Manager or "TOM"), and results from the first phase of observations, including the detection of activity and outbursts on the giant comet C/2014 UN271 (Bernardinelli-Bernstein) and the discovery and follow-up of outbursts on comets. Within these outburst discoveries, we present a high cadence of 7P/Pons-Winnecke with days, a large outburst on 57P/duToit-Neujmin-Delporte, and evidence that comet P/2020 X1 (ATLAS) was in outburst when discovered.

In the context of Be stars, we restudied the viscous transonic decretion disk model of these stars. This model is driven by a radiative force due to an ensemble of optically-thin lines and viscosity considering the Shakura Sunyaev prescription. The non-linear equation of motion presents a singularity (sonic point) and an eigenvalue, which is also the initial condition at the stellar surface. Then, to obtain this eigenvalue, we set it as a radial quantity and perform a detailed topological analysis. Thereafter, we describe a numerical method for solving either Nodal and Saddle transonic solutions. The value of the viscosity,"alpha", barely determine the location of the sonic point, but it determines the topology of the solution. We found two Nodal solutions, which are almost indistinguishable between them. Saddle solutions are founded for lower values of "alpha" than the required of the Nodal solutions. In addition, rotational velocity do not play a determine role in the velocity (and density) profile, because viscosity effects collapse all the solutions to almost a unique one in a small region above the stellar surface. A suitable combination of line force parameters and/or disk temperature, give location of the sonic point lower than 50 stellar radii, describing a truncated disk. This could explain the SED turndown observed in Be stars without needing a binary companion.

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 2 observations of CO(2-1) emission from the circumnuclear disks in two early-type galaxies, NGC 1380 and NGC 6861. The disk in each galaxy is highly inclined ($i\,{\sim}\,75^{\circ}$), and the projected velocities of the molecular gas near the galaxy centers are ${\sim}300\,\mathrm{km \, s^{-1}}$ in NGC 1380 and ${\sim}500\,\mathrm{km \, s^{-1}}$ in NGC 6861. We fit thin disk dynamical models to the ALMA data cubes to constrain the masses of the central black holes (BHs). We created host galaxy models using Hubble Space Telescope images for the extended stellar mass distributions and incorporated a range of plausible central dust extinction values. For NGC 1380, our best-fit model yields $M_{\mathrm{BH}} = 1.47 \times 10^8\,M_{\odot}$ with a ${\sim}40\%$ uncertainty. For NGC 6861, the lack of dynamical tracers within the BH's sphere of influence due to a central hole in the gas distribution precludes a precise measurement of $M_{\mathrm{BH}}$. However, our model fits require a value for $M_{\mathrm{BH}}$ in the range of $(1-3) \times 10^9\,M_{\odot}$ in NGC 6861 to reproduce the observations. The BH masses are generally consistent with predictions from local BH-host galaxy scaling relations. Systematic uncertainties associated with dust extinction of the host galaxy light and choice of host galaxy mass model dominate the error budget of both measurements. Despite these limitations, the measurements demonstrate ALMA's ability to provide constraints on BH masses in cases where the BH's projected radius of influence is marginally resolved or the gas distribution has a central hole.

Accretion of interplanetary dust onto gas giant exoplanets is considered. Poynting-Robertson drag causes dust particles from distant reservoirs to slowly inspiral toward the star. Orbital simulations for the three-body system of the star, planet, and dust particle show that a significant fraction of the dust may accrete onto massive planets in close orbits. The deceleration of the supersonic dust in the planet's atmosphere is modeled, including ablation by thermal evaporation and sputtering. The fraction of the accreted dust mass deposited as gas-phase atoms is found to be large for close-in orbits and massive planets. If mass outflow and vertical mixing are sufficiently weak, the accreted dust produces a constant mixing ratio of atoms and remnant dust grains below the stopping layer. When vertical mixing is included along with settling, the solutions interpolate between the mixing ratio due to the meteoric source above the homopause, and that of the well-mixed deeper atmosphere below the homopause. The line opacity from atoms and continuum opacity from remnant dust may be observable in transmission spectra for sufficiently large dust accretion rates, a grain size distribution tilted toward the blowout size, and sufficiently weak vertical mixing. If mixing is strong, the meteoric source may still act to augment heavy elements mixed up from the deep atmosphere as well as provide nucleation sites for the formation of larger particles. The possible role of the Lorentz drag force in limiting the flow speeds and mixing coefficient for pressures $\la 1\, \rm mbar$ is discussed.

To understand how well galaxies, gas and intracluster stars trace dark matter in and around galaxy clusters, we use the IllustrisTNG cosmological hydrodynamical simulation and compare the spatial distribution of dark matter with those of baryonic components in clusters. To quantify the global morphology of the density distribution of each component in clusters, we fit an ellipse to the density contour of each component and derive shape parameters at different radii. We find that ellipticity of dark matter is better correlated with that of galaxy mass-weighted number density, rather than with that of galaxy number density or galaxy velocity dispersion. We thus use the galaxy mass-weighted number density map as a representative of the galaxy maps. Among three different density maps from galaxies, gas, and intracluster stars, the ellipticity of dark matter is best reproduced by that of the galaxy map over the entire radii. The 'virialized' galaxy clusters show a better correlation of spatial distribution between dark matter and other components than the 'unvirialized' clusters, suggesting that it requires some time for each component to follow the spatial distribution of dark matter after merging events. Our results demonstrate that galaxies are still good tracers of dark matter distribution even in the non-linear regime corresponding to the scales in and around galaxy clusters, being consistent with the case where galaxies trace well the matter distribution in cosmologically large scales.

In this Letter, we give a detailed analysis to the M3.3 class flare that occurred on August 17, 2013 (SOL2013-08-17T18:16). It presents a clear picture of mutual magnetic interaction initially from the photosphere to the corona via the abrupt rapid shearing motion of a small sunspot before the flare, and then suddenly from the corona back to the photosphere via the sudden retraction motion of the same sunspot during the flare impulsive phase. About 10 hours before the flare, a small sunspot in the active region NOAA 11818 started to move northeast along a magnetic polarity inversion line (PIL), creating a shearing motion that changed the quasi-static state of the active region. A filament right above the PIL was activated following the movement of the sunspot and then got partially erupted. The eruption eventually led to the M3.3 flare. The sunspot was then suddenly pulled back to the opposite direction upon the flare onset. During the backward motion, the Lorentz force underwent a simultaneous impulsive change both in magnitude and direction. Its directional change is found to be conformable with the retraction motion. The observation provides direct evidence for the role of the shearing motion of the sunspot in powering and triggering the flare. It especially confirms that the abrupt motion of a sunspot during a solar flare is the result of a back reaction caused by the reconfiguration of the coronal magnetic field.

Winds are commonly observed in luminous active galactic nuclei (AGNs). A plausible model of those winds is magnetohydrodynamic (MHD) disc winds. In the case of disc winds from a thin accretion disc, isothermal or adiabatic assumption is usually adopted in such MHD models. In this work we perform two-dimensional MHD simulations implementing different thermal treatments (isothermal, adiabatic and radiative) to study their effects on winds from a thin accretion disc. We find that both the isothermal model and the adiabatic model overestimate the temperature, underestimate the power of disc winds, and cannot predict the local structure of the winds, compared to the results obtained by solving the energy equation with radiative cooling and heating. Based on the model with radiative cooling and heating, the ionization parameter, the column density and the velocity of the disc winds have been compared to the observed ultrafast outflows (UFOs). We find that in our simulations the UFOs can only be produced inside hundreds of Schwarzschild radius. At much larger radii, no UFOs are found. Thus, the pure MHD winds cannot interpret all the observed UFOs.

Stars do not simply pop up on the main sequence. Before the stars arrive on the zero-age main sequence, they form in the collapses of molecular clouds, gain matter through accretion processes, and compress their cores until hydrogen can burn in full equilibrium. Although this evolutionary phase lasts a relatively short time, it is the imprint of these important physical processes that is often ignored by simplified assumptions. While asteroseismology offers a great tool to investigate these physical processes, studying pre-MS oscillations in turn has the potential to further advance the field. Asteroseismology of pre-main sequence stars faces observational and theoretical challenges. The remnants of their birth environment which is often still surrounding the young stars causes variability that can interfere with the signal of pulsations. The lack of long time-base satellite observations in addition limits the applications of the method. Theoretical models of pre-main sequence stars include several assumptions and simplifications that influence the calculation of pulsation frequencies and excitation properties of pulsation modes. Keeping all this in mind, the prospects for pre-main sequence asteroseismology are manifold. An improved understanding of the structure of young stellar objects has the potential to answer some of the open questions of stellar evolution, including angular momentum transport and the formation of magnetic fields. While gyrochronology, for example, struggles to determine the ages of the youngest clusters, pulsations in pre-main sequence stars can function as an independent age indicator yielding higher precision for single stars. The increasing interest of stellar astrophysics in general to investigate the formation and early evolution of stars and planets illustrates the growing importance of pre-main sequence asteroseismology.

We present Tricolour, a package for Radio Frequency Interference mitigation of wideband finely channelized MeerKAT correlation data. The MeerKAT passband is heavily affected by interference from satellite, mobile, aircraft and terrestrial transponders. Coupled with typical data rates in excess of 100 GiB/hr at 208kHz channelization resolution, mitigation poses a significant processing challenge. Our flagger is highly configurable, parallel and optimized, employing Dask and Numba technologies to implement the widely used SumThreshold and MAD interference detection algorithms. We find that typical 208kHz channelized datasets can be processed at rates in excess of 400 GiB/hr for a typical L-band flagging strategy on a modern dual-socket Intel Xeon server.

Observational and theoretical evidence suggests that a substantial population of molecular clouds (MCs) appear to be unbound, dominated by turbulent motions. However, these estimations are made typically via the so-called viral parameter $\alpha_{\rm vir}^{\rm class}$, which is an observational proxy to the virial ratio between the kinetic and the gravitational energy. This parameter intrinsically assumes that MCs are isolated, spherical, and with constant density. However, MCs are embedded in their parent galaxy and thus are subject to compressive and disruptive tidal forces from their galaxy, exhibit irregular shapes, and show substantial substructure. We, therefore, compare the typical estimations of $\alpha_{\rm vir}^{\rm class}$ to a more precise definition of the virial parameter, $\alpha_{\rm vir}^{\rm full}$, which accounts not only for the self-gravity (as $\alpha_{\rm vir}^{\rm class}$), but also for the tidal stresses, and thus, it can take negative (self-gravity) and positive (tides) values. While we recover the classical result that most of the clouds appear to be unbound, having $\alpha_{\rm vir}^{\rm class} > 2$, we show that, with the more detailed definition considering the full gravitational energy, (i) 50\%\ of the total population is gravitationally bound, however, (ii) another 20\%\ is gravitationally dominated, but with tides tearing them apart; (iii) the source of those tides does not come from the galactic structure (bulge, halo, spiral arms), but from the molecular cloud complexes in which clouds reside, and probably (iv) from massive young stellar complexes, if they were present. (v) Finally, our results also suggest that, interstellar turbulence can have, at least partially, a gravitational origin.

Mergers of galaxy clusters are the most energetic events in the universe, driving shock and cold fronts, generating turbulence, and accelerating particles that create radio halos and relics. The galaxy cluster CL 0217+70 is a remarkable late stage merger, with a double peripheral radio relic and a giant radio halo. A Chandra study detects surface brightness edges that correspond to radio features within the halo. In this work, we present a study of this cluster with NuSTAR and Chandra data using spectro-imaging methods. The global temperature is found to be kT = 9.1 keV. We set an upper limit for the IC flux of ~2.7x10^(-12) erg s^(-1) cm^(-2), and a lower limit to the magnetic field of 0.08 microG. Our local IC search revealed a possibility that IC emission may be dominant at the outskirts of a radio halo emission and on/near shock regions within ~0.6 r500 of clusters. We detected a "hot spot" feature in our temperature map coincident a surface brightness edge, but our investigation on its origin is inconclusive. If the "hot spot" is the downstream of a shock, we set a lower limit of kT > 21 keV to the plasma, that corresponds to M~2. We found three shock fronts within 0.5 r500. Multiple weak shocks within the cluster center hint at an ongoing merger activity and continued feeding of the giant radio halo. CL 0217+70 is the only example hosting these secondary shocks in multiple form.

We identify and characterise a Milky Way-like realisation from the Auriga simulations with two consecutive massive mergers $\sim2\,$Gyr apart at high redshift, comparable to the reported Kraken and Gaia-Sausage-Enceladus. The Kraken-like merger ($z=1.6$, $M_{\rm Tot} = 8\times10^{10}\,$M$_{\odot}$) is gas-rich, deposits most of its mass in the inner $10\,$kpc, and is largely isotropic. The Sausage-like merger ($z=1.14$, $M_{\rm Tot} = 1\times10^{11}\,$M$_{\odot}$) leaves a more extended mass distribution at higher energies, and has a radially anisotropic distribution. For the higher redshift merger, the stellar mass ratio of the satellite to host galaxy is 1:3. As a result, the chemistry of the remnant is indistinguishable from contemporaneous in-situ populations, making it challenging to identify this component through chemical abundances. This naturally explains why all abundance patterns attributed so far to Kraken are in fact fully consistent with the metal-poor in-situ so-called Aurora population and thick disc. However, our model makes a falsifiable prediction: if the Milky Way underwent a gas-rich double merger at high redshift, then this should be imprinted on its star formation history with bursts about $\sim2\,$Gyrs apart. This may offer constraining power on the highest-redshift major mergers.

The gravitational lens system PS J0147+4630 (Andromeda's Parachute) consists of four quasar images ABCD and a lensing galaxy. We obtained $r$-band light curves of ABCD in the 2017$-$2021 period from a monitoring with two 2-m class telescopes. These curves and state-of-the-art curve-shifting algorithms led to three independent time delays relative to image A, one of which is accurate enough (uncertainty of about 4%) to be used in cosmological studies. Our finely sampled light curves and some additional fluxes in the years 2010$-$2013 also demonstrated the presence of significant microlensing variations. This paper also focused on new near-IR spectra of ABCD in 2018$-$2019 that were derived from archive data of two 10-m class telescopes. We analysed the spectral region including the MgII, H$\beta$, [OIII], and H$\alpha$ emission lines (0.9$-$2.4 $\mu$m), measuring image flux ratios and a reliable quasar redshift of 2.357 $\pm$ 0.002, and finding evidence of an outflow in the H$\alpha$ emission. In addition, we updated the lens mass model of the system and estimated a quasar black-hole logarithmic mass ${\log \left[ M_{\rm{BH}}/\rm{M_{\odot}} \right]}$ = 9.34 $\pm$ 0.30.

OB stars powering stellar bowshock nebulae (SBNe) have been presumed to have large peculiar velocities. We measured peculiar velocities of SBN central stars to assess their kinematics relative to the general O star population using $Gaia$ EDR3 data for 267 SBN central stars and a sample of 455 Galactic O stars to derive projected velocities $v_{\rm 2D}$. For a subset of each sample we obtained new optical spectroscopy to measure radial velocities and identify multiple-star systems. We find a minimum multiplicity fraction of 36$\pm$6% among SBN central stars, consistent with $>$28% among runaway Galactic O stars. The large multiplicity fraction among runaways implicates very efficient dynamical ejection rather than binary-supernova origins. The median $v_{\rm 2D}$ of SBN central stars is $v_{\rm 2D}$=14.6 km s$^{-1}$, larger than the median $v_{\rm 2D}$=11.4 km s$^{-1}$ for non-bowshock O stars. Central stars of SBNe have a runaway ($v_{\rm 2D}$$>$25 km s$^{-1}$) fraction of 24$^{+9}_{-7}$%, consistent with the 22$^{+3}_{-3}$% for control-sample O stars. Most (76%) of SBNe central stars are not runaways. Our analysis of alignment ($\Delta_{\rm PA}$) between the nebular morphological and $v_{\rm 2D}$ kinematic position angles reveals two populations: a highly aligned ($\sigma_{PA}$=25$^\circ$) population that includes stars with the largest $v_{\rm 2D}$ (31% of the sample) and a random (non-aligned) population (69% of the sample). SBNe that lie within or near HII regions comprise a larger fraction of this latter component than SBNe in isolated environments, implicating localized ISM flows as a factor shaping their orientations and morphologies. We outline a new conceptual approach to computing the Solar LSR motion, yielding [U$_\odot$, V$_\odot$, W$_\odot$]= [5.5, 7.5, 4.5] km s$^{-1}$.

We present the pulse arrival times and high-precision dispersion measure estimates for 14 millisecond pulsars observed simultaneously in the 300-500 MHz and 1260-1460 MHz frequency bands using the upgraded Giant Metrewave Radio Telescope (uGMRT). The data spans over a baseline of 3.5 years (2018-2021), and is the first official data release made available by the Indian Pulsar Timing Array collaboration. This data release presents a unique opportunity for investigating the interstellar medium effects at low radio frequencies and their impact on the timing precision of pulsar timing arrary experiments. In addition to the dispersion measure time series and pulse arrival times obtained using both narrowband and wideband timing techniques, we also present the dispersion measure structure function analysis for selected pulsars. Our ongoing investigations regarding the frequency dependence of dispersion measures have been discussed. Based on the preliminary analysis for five millisecond pulsars, we do not find any conclusive evidence of chromaticity in dispersion measures. Data from regular simultaneous two frequency observations are presented for the first time in this work. This distinctive feature leads us to the highest precision dispersion measure estimates obtained so far for a subset of our sample. Simultaneous multi-band uGMRT observations in Band 3 and Band 5 are crucial for high-precision dispersion measure estimation and for the prospect of expanding the overall frequency coverage upon the combination of data from the various Pulsar Timing Array consortia in the near future. Parts of the data presented in this work are expected to be incorporated into the upcoming third data release of the International Pulsar Timing Array.

We decompose the monthly cosmic-ray data, using several neutron monitor count rates, of Cycles 19-24 with principal component analysis (PCA). We show using different cycle limits that the first and second PC of cosmic-ray (CR) data explain 77-79% and 13-15% of the total variation of the Oulu CR Cycles 20-24 (C20- C24), 73-77% and 13-17% of the variation of Hermanus C20-C24, and 74-78% and 17-21% of the Climax C19-C22, respectively. The PC1 time series of the CR Cycles 19-24 has only one peak in its power spectrum at the period 10.95 years, which is the average solar cycle period for the interval SC19-SC24. The PC2 time series of the same cycles has a clear peak at period 21.90 (Hale cycle) and another peak at 1/3 of that period with no peak at the solar cycle period. We show that the PC2 of the CR is essential in explaining the differences in the intensities of the even and odd cycles of the CR. The odd cycles have positive phase in the first half and negative phase in the second half of their PC2. This leads to slow decrease of the intensity in the beginning of the cycle and at minimum for the odd cycles. On the contrary, for the even cycles the phases are vice versa and this leads to faster decrease and more rapid recovery in the CR intensity of the cycle. As a consequence the even cycles have more peak-like structure. The only exceptions of this rule are Cycles 23 and 24 such the former has almost zero line PC2, and the latter has similar PC2 than the earlier odd cycles. These results are confirmed with skewness-kurtosis (S-K) analysis. Furthermore, S-K shows that other even and odd cycles, except Cycle 21, are on the regression line with correlation coefficient 0.85. The Cycles 21 of all calculated eight stations are compactly located in the S -K coordinate system and have smaller skewnesses and higher kurtoses than the odd Cycles 23.

Homogeneous coronal data set (HCDS) of the green corona (Fe XIV) and coronal index of the solar activity (CI) have been used to study time-latitudinal distribution in solar cycles 18-24 and compared with similar distribution of sunspots, the magnetic fields and the solar radio flux 10.7 cm. The most important results are: (a) distribution of coronal intensities related to the cycle maximum are different for individual cycles, (b) the poleward migration of the HCDS from mid latitudes in each cycle exists, even in extremely weak Cycle 24, and the same is valid for the equatorward migration (c) the overall values of HCDS are slightly stronger for the northern hemisphere than for the southern one, (d) distribution of the HDCS are in coincidence with strongest photospheric magnetic fields (B>50 Gauss) and histogram of the sunspot groups, (e) Gnevyshev gap was confirmed with at least 95% confidence in the CI, however, with different behavior for odd and even cycles. Principal component analysis (PCA) showed that the first and second component account for 87.7% and 7.3% of the total variation of the CI. Furthermore, the PC2 of the green corona was quite different for cycle 21, compared with other cycles.

PKS 1413+135, a typical compact symmetric object (CSO) with a two-side pc-scale structure in its miniature radio morphology, is spatially associated with the source 4FGL J1416.1+1320 and recently is detected in the TeV band with the MAGIC telescopes. We present the analysis of its X-ray and gamma-ray observations obtained with Swift, XMM-Newton, Chandra, and Fermi-LAT for revealing its high energy radiation physics. No significant variation trend is observed in the X-ray band. Its gamma-ray light curve derived from the Fermi-LAT 13.5-year observations shows that it is in a low gamma-ray flux stage before MJD 58500 and experiences violent outbursts after MJD 58500. The confidence level of the flux variability is much higher than 5\sigma, and the flux at 10 GeV varies ~3 orders of magnitude. The flux variation is accompanied by the clearly spectral variation. The spectral shape displays a log-parabola spectrum in low-flux state and a hard power-law spectrum in the gamma-ray flares. The highest energy of the detected photons by the Fermi-LAT is ~236 GeV, which is detectable with the MAGIC telescopes. Furthermore, we compile the broadband SED during an GeV gamma-ray outburst and fit it with a two-zone leptonic model and emission in the X-ray-gamma-ray band is attributed to the inverse Compton scattering processes. The result shows that the gamma-rays of PKS 1413+135 would be detectable with the MAGIC telescopes, whether the source is located at z=0.247 or z=0.5. Based on our analysis and its CSO radio morphology, we suspect that the nuclear jet activity of PKS 1413+135 is episodic, the weak gamma-ray emission before MJD 58500 may be from its pc-scale jet structures powered by previous activities, and the violent outbursts with short timescale variability after MJD 58500 could be attributed to the recently re-started jet activity.

In this paper, we investigate the black hole (BH) spin contribution to jet power, especially for the magnetic arrested disk (MAD), where only inner accretion disk luminosity is closely coupled with the spin-jet power, and try to explain the `outliers' track of the radio $L_{\rm R}$ to X-ray luminosity $L_{\rm X}$ in two black hole X-ray binaries (BHXBs). Our results suggest that the BZ-jet and the inner-disk coupling could account for the `outliers' track of the radio/X-ray correlation in two BHXBs, H1743-322 and MAXI J1348-630. Although the accretion disk of H1743-322 in the outburst could be in the MAD state, there is a lower probability that MAXI J1348-630 is in the MAD state due to its low jet production efficiency. The difference in the inner-disk bolometric luminosity ratio of the two sources implies that these two BHXBs are in different inner-disk accretion states. We further investigate the phase-changing regime of MAXI J1348-630 and find that there is a phase transition around $L_{\rm X}/L_{\rm Edd}\sim 10^{-3}$. The assumption of sub-MAD is discussed as well.

Atmospheric studies at high spectral resolution have shown the presence of molecules, neutral and ionised metals, and hydrogen in the transmission spectrum of ultra-hot Jupiters, and have started to probe the dynamics of their atmospheres. We analyse the transmission spectrum of MASCARA-1b, one of the densest ultra-hot Jupiters orbiting a bright (V=8.3) star. We focus on the CaII H&K, NaI, LiI, H$\alpha$, and KI D1 spectral lines and on the cross-correlated FeI, FeII, CaI, YI, VI, VII, CaH, and TiO lines. For those species that are not present in the stellar spectrum, no detections are reported, but we measure upper limits with an excellent precision ($\sim10$ ppm for particular species). For those species that are present in the stellar spectrum and whose planet-occulted spectral lines induce spurious features in the planetary transmission spectrum, an accurate modelling of the Rossiter-McLaughlin effect (RM) and centre-to-limb variations (CLV) is necessary to recover possible atmospheric signals. In the case of MASCARA-1b, this is difficult due to the overlap between the radial velocities of the stellar surface regions occulted by MASCARA-1b and the orbital track along which the planet atmospheric signal is expected to be found. To try to disentangle a possible planetary signal, we compare our results with models of the RM and CLV effects, and estimate the uncertainties of our models depending on the different system parameters. Unfortunately, more precise measurements of the spin-orbit angle are necessary to better constrain the planet-occulted track and correct for the transit effects in the transmission spectrum with enough precision to be able to detect or discard possible planetary absorptions. Finally, we discuss the possibility that non-detections are related to the low absorption expected for a high surface gravity planet such as MASCARA-1b.

Employing both the number density and energy density, the equation of state of fermions for each microscopic region of a neutron star is derived in the presence of various statistical effects. It is computed as functions of temperature (T), chemical potential (\mu) and magnetic field (B) based on the statistical properties and corresponding phases of neutron stars. It is revealed that the measurable properties describing the dynamic behavior of a neutron star are determined by the unusual statistical conditions. It is demonstrated that, in a strong magnetic field (B larger than\mu), the pressure (P) shows a logarithmic dependence on B (P is proportional to logarithm of(\mu/B). It is found that the regions with a higher (lower) magnetic field (density) have a shorter (longer) life time. Moreover, the relationship between the magnetic field and chemical potential confirms that the large magnetic field and high density are short lived as compared to low magnetic field and low density at the same temperature.

Reconstructing images from very long baseline interferometry (VLBI) data with sparse sampling of the Fourier domain (uv-coverage) constitutes an ill-posed deconvolution problem. It requires application of robust algorithms maximizing the information extraction from all of the sampled spatial scales and minimizing the influence of the unsampled scales on image quality. We develop a new multiscale wavelet deconvolution algorithm DoG-HiT for imaging sparsely sampled interferometric data which combines the difference of Gaussian (DoG) wavelets and hard image thresholding (HiT). Based on DoG-HiT, we propose a multi-step imaging pipeline for analysis of interferometric data. DoG-HiT applies the compressed sensing approach to imaging by employing a flexible DoG wavelet dictionary which is designed to adapt smoothly to the uv-coverage. It uses closure properties as data fidelity terms only initially and perform non-convex, non-smooth optimization by an amplitude conserving and total flux conserving hard thresholding splitting. DoG-HiT calculates a multiresolution support as a side product. The final reconstruction is refined through self-calibration loops and imaging with amplitude and phase information applied for the multiresolution support only. We demonstrate the stability of DoG-HiT and benchmark its performance against image reconstructions made with CLEAN and Regularized Maximum-Likelihood (RML) methods using synthetic data. The comparison shows that DoG-HiT matches the superresolution achieved by the RML reconstructions and surpasses the sensitivity to extended emission reached by CLEAN. Application of regularized maximum likelihood methods outfitted with flexible multiscale wavelet dictionaries to imaging of interferometric data matches the performance of state-of-the art convex optimization imaging algorithms and requires fewer prior and user defined constraints.

We present the results from a systematic X-ray analysis of Subaru optically-selected clusters with high cluster richness in the eROSITA Final Equatorial-Depth Survey (eFEDS) field. By a joint analysis of SRG(Spectrum Roentgen Gamma)/eROSITA and Subaru/HSC surveys, we aim to investigate the dynamical status of the optically-selected clusters and derive the cluster scaling relations. The sample consists of 43 optically-selected galaxy clusters with richness $>40$ in the redshift range of 0.16--0.89. We systematically analyzed X-ray images and emission spectra using the eROSITA data. We identified the brightest cluster galaxy (BCG) using the optical and far-infrared databases and measured the offset between X-ray peak and BCG position as an indicator of the dynamical status. We investigated the luminosity-temperature and mass-luminosity relations based on eROSITA X-ray spectra and HSC weak-lensing data analyses. The fraction of relaxed clusters is estimated to be $2(<39)$\%, which is smaller than that of the X-ray-selected cluster samples. After correcting for a selection bias due to the richness cut, we obtained a shallow $L-T$ slope of $2.0\pm0.4$, which is consistent with the predictions of the self-similar model and the new baseline model by Fujita & Aung (2019). The $L-M$ slope of $1.5\pm0.3$ agrees with the above theoretical models and that of the shear-selected clusters in the eFEDs field. Our analysis of the high-richness optical clusters yielded a small fraction of relaxed clusters and a shallow slope of the luminosity-temperature relation. This suggests that the average X-ray properties of the optical clusters are likely to be different from those observed in the X-ray samples. Thus, the joint eROSITA and HSC observations are powerful in extending the analysis to a larger sample and understanding the selection effect towards establishing the cluster scaling relations.

The highly tangled magnetic field of a magnetar supports shear waves similar to Alfv\'en waves in an ordered magnetic field. Here we explore if torsional modes excited in the stellar interior and magnetosphere can explain the quasi-periodic oscillations (QPOs) observed in the tail of the giant flare of SGR 1900+14. We solve the initial value problem for a tangled magnetic field that couples interior shear waves to relativistic Alfv\'en shear waves in the magnetosphere. Assuming stellar oscillations arise from the sudden release of magnetic energy, we obtain constraints on the energetics and geometry of the process. If the flare energy is deposited initially inside the star, the wave energy propagates relatively slowly to the magnetosphere which is at odds with the observed rise time of the radiative event of $\lesssim 10$ ms. Nor can the flare energy be deposited entirely outside the star, as most of the energy reflects off the stellar surface, giving surface oscillations of insufficient magnitude to produce detectable modulations of magnetospheric currents. Energy deposition in a volume that straddles the stellar surface gives agreement with the observed rise time and excites a range of modes with substantial amplitude at observed QPO frequencies. In general, localized energy deposition excites a broad range of modes that encompasses the observed QPOs, though many more modes are excited than the number of observed QPOs. If the flare energy is deposited axisymmetrically, as is possible for a certain class of MHD instabilities, the number of modes that is excited is considerably reduced.

Context. Level population inversion of hydrogen atoms in ionized gas may lead to stimulated emission of hydrogen recombination lines, and the level populations can in turn be affected by powerful stimulated emissions. Aims. In this work the interaction of the radiation fields and the level population inversion of hydrogen atoms is studied. The effect of the stimulated emissions on the line profiles is also investigated. Methods. Our previous nl-model for calculating level populations of hydrogen atoms and hydrogen recombination lines is improved. The effects of line and continuum radiation fields on the level populations are considered in the improved model. By using this method the properties of simulated hydrogen recombination lines and level populations are used in analyses. Results. The simulations show that hydrogen radio recombination lines are often emitted from the energy level with an inverted population. The widths of Hn$\alpha$ lines can be significantly narrowed by strong stimulated emissions to be even less than 10 km s$^{-1}$. The amplification of hydrogen recombination lines is more affected by the line optical depth than by the total optical depth. The influence of stimulated emission on the estimates of electron temperature and density of ionized gas is evaluated. We find that comparing multiple line-to-continuum ratios is a reliable method for estimating the electron temperature, while the effectiveness of the estimation of electron density is determined by the relative significance of stimulated emission.

We present an analytical model, RAiSE HD, for the evolution of extended active galactic nuclei (AGNs) throughout their full lifecycle, including the initial jet expansion, lobe formation, and eventual remnant phases. A particular focus of our contribution is on the early jet expansion phase, which is traditionally not well captured in analytical models. We implement this model within the \emph{Radio AGN in Semi-Analytic Environments} (RAiSE) framework, and find that the predicted radio source dynamics are in good agreement with hydrodynamic simulations of both Fanaroff-Riley Type-I and -II radio lobes over several orders of magnitude in jet kinetic power. We construct synthetic synchrotron and inverse-Compton surface brightness images by complementing the RAiSE model with the magnetic field and shock-acceleration histories of a set of Lagrangian tracer particles taken from an existing hydrodynamic simulation. We show that a single set of particles is sufficient for an accurate description of the dynamics and observable features of Fanaroff-Riley Type-II radio lobes with very different jet and environment parameters. Our new model, RAiSE HD (RAiSE+hydrodynamics), predicts that the lobes of young ($\lesssim 10$ Myr) sources will be both longer and brighter than expected at the same age from existing analytical models which lack a jet-dominated expansion phase; this finding has important implications for interpretation of radio galaxy observations. The RAiSE HD code, written in Python, is publicly available on GitHub and PyPI.

The Seyfert 1 galaxy NGC 7469 possesses a prominent nuclear starburst ring and a luminous active galactic nucleus (AGN). Evidence of an outflow in the innermost nuclear region has been found in previous works. We detect the ionized gas outflow on a larger scale in the galaxy using the archival VLT/MUSE and {\em Chandra} observations. The optical emission lines are modeled using two Gaussian components, and a non-parametric approach is applied to measure the kinematics of [O III] and $\rm H\alpha$ emitting gas. Line ratio diagnostics and spatially resolved maps are derived to examine the origin of the outflow. The kpc-scale kinematics of [O III] is dominated by a blueshifted component whereas velocity map of $\rm H\alpha$ shows a rotational disk with complex non-rotational substructure. The starburst wind around the circumnuclear ring is confirmed, and we find evidence of an AGN-driven outflow extending to a radial distance of $\rm \sim2$ kpc from the nucleus, with a morphology consistent with a nearly face-on ionization cone. The previously reported circumnuclear outflow resembles part of the bright base. We derive mass and energy outflow rates for both the starburst wind and the AGN-driven outflow. The estimated kinetic coupling efficiency of the kpc-scale AGN outflow is $\dot{E}_{\rm out}/L_{\rm bol}\sim 0.1\%$, lower than the threshold predicted by the ``two-stage'' theoretical model for effective feedback. Our results reinforce the importance of spatially resolved study to disentangle feedback where AGN and starburst coexist, which may be common during the cosmic noon of black hole and galaxy growth.

Recording a stellar occultation is one powerful method that gives direct information about the physical properties of the occulting solar system object. In order to obtain reliable and accurate results, simultaneous observations from different locations across-track of the projected path are of great importance. However, organising all the observing stations, aggregating, and analysing the data is time-consuming and not that easy. We have developed a web portal named Occultation Portal (OP) to manage all those occultation observation campaigns from a central server. With this portal, the instrumental and observational information of all observers participating in a stellar occultation campaign and the concerned data are archived systematically in a standard format. The researchers can then visualise the archived data on an event basis. The investigators can also extract the light curve for each data-set with the added reduction pipeline to the portal base. This paper describes in detail the portal structure and the developed features.

The solar corona is host to a continuous flow of propagating disturbances (PD). These are continuous and ubiquitous across broad regions of the corona, including the quiet Sun. The aim of this paper is to present an improved, efficient method to create velocity vector field maps, based on the direction and magnitude of the PD as observed in time series of extreme ultraviolet (EUV) images. The method is presented here for use with the Atmospheric Imaging Assembly (AIA)/Solar Dynamics Observatory (SDO) EUV channels, and takes as input \app2 hours of images at the highest 12s cadence. Data from a region near disk center is extracted, and a process called time normalization applied to the co-aligned data. Following noise reduction using \atrous\ decomposition, the PD are effectively revealed. A modified Lucas Kanade algorithm is then used to map the velocity field. The method described here runs comfortably on a desktop computer in a few minutes, and offers an order of magnitude improvement in efficiency compared to a previous implementation. Applied to a region of the quiet Sun, we find that the velocity field describes a mosaic of cells of coherent outwardly-diverging PD flows, of typical size 50 to 100\arcsec\ (36 to 72Mm). The flows originate from points and narrow corridors in the cell centres, and end in the narrow boundaries between cells. Visual comparison with ultraviolet AIA images shows that the flow sources are correlated with the bright photospheric supergranular network boundaries. Assuming that the PD follow the local magnetic field, the velocity flow field is a proxy for the plane-of-sky distribution of the coronal magnetic field, and therefore the maps offer a unique insight into the topology of the corona. These are particularly valuable for quiet Sun regions where the appearance of structures in EUV images is hard to interpret.

The Beamforming Elevated Array for COsmic Neutrinos (BEACON) is a planned neutrino telescope designed to detect radio emission from upgoing air showers generated by ultrahigh energy tau neutrino interactions in the Earth. This detection mechanism provides a measurement of the tau flux of cosmic neutrinos. We have installed an 8-channel prototype instrument at high elevation at Barcroft Field Station, which has been running since 2018, and consists of 4 dual-polarized antennas sensitive between 30-80 MHz, whose signals are filtered, amplified, digitized, and saved to disk using a custom data acquisition system (DAQ). The BEACON prototype is at high elevation to maximize effective volume and uses a directional beamforming trigger to improve rejection of anthropogenic background noise at the trigger level. Here we discuss the design, construction, and calibration of the BEACON prototype instrument. We also discuss the radio frequency environment observed by the instrument, and categorize the types of events seen by the instrument, including a likely cosmic ray candidate event.

After their initial formation, disk galaxies are observed to be rotationally stable over periods of >6 Gyr, implying that any large velocity disturbances of stars and gas clouds are damped rapidly on the timescale of their rotation. However, it is also known that despite this damping, there must be a degree of random local motion to stabilize the orbits against degenerate collapse. A mechanism for such damping is proposed by a combination of inter-stellar gravitational interactions, and interactions with the Oort clouds and exo-Oort objects associated with each star. These mechanisms may produce rapid damping of large perturbations within a time period that is short on the scale of observational look-back time, but long on the scale of the disk rotational period for stars with small perturbations. This mechanism may also account for the locally observed mean perturbations in the Milky Way of 8-15 km/s for younger stars and 20-30 km/s for older stars.

In the Solar System short-lived radioisotopes, such as 26Al, played a crucial role during the formation planetary bodies by providing a significant additional source of heat. Notably, this led to early and large-scale melting and iron core formation in planetesimals and their loss of volatile elements, such as hydrogen and carbon. In the context exoplanetary systems, therefore, the prevalence of short-lived radioisotopes is key to interpreting the observed bulk volatile budget and atmospheric diversity among low-mass exoplanets. White dwarfs that have accreted planetary material provide a unique means to infer the frequency of iron core formation in extrasolar planetesimals, and hence the ubiquity of planetary systems forming with high short-lived radioisotope abundances. Here, we devise a quantitative method to infer the fraction of planetary systems enriched with shortlived radionuclides upon planetesimal formation from white dwarf data. We argue that the current evidence from white dwarfs point towards a significant fraction of exo-planetesimals having formed an iron core. Although the data may be explained by the accretion of exo-moon or Pluto-sized bodies that were able to differentiate due to gravitational potential energy release, our results suggest that the most likely explanation for the prevalence of differentiated material among polluted white dwarfs is that the Solar System is not unusual in being enriched in 26Al. The models presented here suggest a ubiquitous pathway for the enrichment of exoplanetary systems by short-lived radioisotopes, disfavouring short-lived radioisotope enrichment scenarios relying on rare chance encounters with single nearby supernovae, Wolf-Rayet or AGB stars.

The detection of gravitational waves from binary neuron star merger GW170817 and electromagnetic counterparts GRB170817 and AT2017gfo kick-started the field of gravitational wave multimessenger astronomy. The optically red to near infra-red emission (`red' component) of AT2017gfo was readily explained as produced by the decay of newly created nuclei produced by rapid neutron capture (a kilonova). However, the ultra-violet to optically blue emission (`blue' component) that was dominant at early times (up to 1.5 days) received no consensus regarding its driving physics. Among many explanations, two leading contenders are kilonova radiation from a lanthanide-poor ejecta component or shock interaction (cocoon emission). In this work, we simulate AT2017gfo-like light curves and perform a Bayesian analysis to study whether an ultra-violet satellite capable of rapid gravitational wave follow-up, could distinguish between physical processes driving the early `blue' component. We find that a Dorado-like ultra-violet satellite, with a 50 sq. deg. field of view and a limiting magnitude (AB) of 20.5 for a 10 minute exposure is able to distinguish radiation components up to at least 160 Mpc if data collection starts within 3.2 or 5.2 hours for two possible AT2017gfo-like light curve scenarios. We also study the degree to which parameters can be constrained with the obtained photometry. We find that, while ultra-violet data alone constrains parameters governing the outer ejecta properties, the combination of both ground-based optical and space-based ultra-violet data allows for tight constraints for all but one parameter of the kilonova model up to 160 Mpc. These results imply that an ultra-violet mission like Dorado would provide unique insights into the early evolution of the post-merger system and its driving emission physics.

Astroaccesible is an outreach project hosted by the Instituto de Astrof\'I sica de Andaluc\'I a - CSIC and leaded by a blind astronomer aimed at the teaching and popularisation of astronomy and astrophysics among all kind of disabled and non-disabled people. Among the different strategies followed to transmit information to blind and partially sighted people, audio description is one of the most accesible and popular in the case of films and museums, but it has not been yet widely incorporated for the description of astronomical images. In this contribution we introduce {"The Universe in words", which are a series of videos describing images of some of the most popular objects in the Messier catalogue. These audio descriptions do not only have a clear inclusive aspect, but also imply a better and deeper understanding of the represented images for everybody. This is one of the most important aspects of using inclusive resources, as they also clearly improve the efficiency of the transmission process for all kind of public. These videos can also be used as supplementary material in of in-person activities and as a complement to other kind of materials, such as sonifications or models of the same or similar type of astronomical objects.

We present an analysis of asteroid (3200) Phaethon using coronagraphic observations from 2008 to 2022 by the COR2 cameras onboard the twin Solar TErrestrial RElations Observatory (STEREO) spacecraft. Although the asteroid cannot be confidently detected in individual images, we managed to spot it in image stacks combined from the same sets of perihelion observations, yet only when observed at low phase angles ($\lesssim$30\deg) but not at large phase angles ($\gtrsim$150\deg). The lack of a strong forward-scattering enhancement that is expected for dust grains having sizes comparable to transmitted wavelengths thereby implies that the perihelion activity of Phaethon is highly unlikely to be relevant to the ejection of dust grains as suggested by earlier studies based on STEREO's HI-1 observations. Assuming the observed activity of Phaethon is caused by dust ejection will lead to an insurmountable discrepancy in the inferred amount of dust no less than an order of magnitude between the HI-1 and COR2 observations. Rather, we speculate that the perihelion activity is caused by sodium and/or iron emissions, the former of which may have become transmittable due to an ageing effect of the HI-1 cameras. The modelled emission flux is qualitatively similar to the HI-1 observations in the case where the peak of the atomic production rate is delayed by $\sim$1 day from perihelion. We encourage future observations of Phaethon at small heliocentric distances to verify our conjecture.

Aims. To create a retrieval framework which encapsulates the 3D nature of exoplanet atmospheres, and to apply it to observed emission phase curve and transmission spectra of hot Jupiter exoplanet WASP-43b. Methods. We present our 3D framework, which is freely available as a stand-alone module from GitHub. We use the atmospheric modelling and Bayesian retrieval package ARCiS (ARtful modelling Code for exoplanet Science) to perform 8 3D retrievals on simultaneous transmission (HST/WFC3) and phase-dependent emission (HST/WFC3, Spitzer/IRAC) observations of WASP-43b as a case study. We assess how input assumptions affect our retrieval outcomes. In particular we look at constraining equilibrium chemistry vs a free molecular retrieval, the case of no clouds vs parametrised clouds, and using Spitzer phase data that have been reduced from two different literature sources. For the free chemistry retrievals, we retrieve abundances of H2O, CH4, CO, CO2, AlO, and NH3 as a function of phase, with many more species considered for the equilibrium chemistry retrievals. Results. We find consistent super-solar C/O (0.6-0.9) and super-solar metallicities (1.7-2.9 dex) for all retrieval setups that assume equilibrium chemistry. Atmospheric heat distribution, hotspot shift (15.6 vs 4.5 for different Spitzer datasets), and temperature structure are influenced by the Spitzer emission phase data. Comparisons are made with other studies of WASP-43b, including available GCM simulations. Conclusions. The parametrised 3D setup we have developed provides a valuable tool to analyse extensive observational datasets such as spectroscopic phase curves. Near-future observations with missions such as the James Webb Space Telescope (JWST) and Ariel will greatly improve our understanding of the atmospheres of exoplanets such as WASP-43b.

We study pre-supernova feedback in a sample of $\sim$ 4700 HII regions in the nearby spiral galaxy M83, identified on their H$\alpha$ emission. We pectroscopically identify Wolf-Rayet (WR) stars populating the star-forming regions. For each HII region, we compute the pressure of ionised gas ($P_{\rm ion}$) and the direct radiation pressure ($P_{\rm dir}$) acting in the region, and investigate how they vary with galactocentric distance, with the physical properties of the region, and with the pressure of the galactic environment ($P_\mathrm{DE}$). For a subset of $\sim$ 500 regions, we also investigate the link between the pressure terms and the properties of the cluster population (age, mass, and LyC flux). We find that $P_{\rm ion}$ dominates over $P_{\rm dir}$ by at least a factor of 10 on average over the disk. Both pressure terms are strongly enhanced and become almost comparable in the central starburst region. In the disk ($R \geq 0.15~R_e$), we observe that $P_{\rm dir}$ stays approximately constant with galactocentric distance. $P_{\rm dir}$ is positively correlated with an increase in radiation field strength (linked to the negative metallicity gradient in the galaxy), while it decreases in low extinction regions. $P_{\rm ion}$ decreases constantly for increasing galactocentric distances. In general, we observe that HII regions near the center are underpressured with respect to their surroundings, whereas regions in the disk are overpressured and hence expanding. We find that regions hosting younger clusters or having more mass in young star clusters have a higher internal pressure, indicating that clustered star formation is likely playing a dominant role in setting the pressure. Finally, we estimate that only 13 % of HII regions hosting young clusters and WR stars have $f_{\rm esc} \geq 0$.[Abridged]

Context: It has been proposed that H$_2$near-infrared lines may be excited by cosmic rays and allow for a determination of the cosmic-ray ionization rate in dense gas. One-dimensional models show that measuring both the H$_2$gas column density and H$_2$line intensity enables a constraint on the cosmic-ray ionization rate as well as the spectral slope of low-energy cosmic-ray protons in the interstellar medium (ISM). Aims: We aim to investigate the impact of certain assumptions regarding the H$_2$chemical models and ISM density distributions on the emission of cosmic-ray induced H$_2$emission lines. This is of particular importance for utilizing observations of these lines with the James Webb Space Telescope to constrain the cosmic-ray ionization rate. Methods: We compare the predicted emission from cosmic-ray induced, ro-vibrationally excited H$_2$emission lines for different one- and three-dimensional models with varying assumptions on the gas chemistry and density distribution. Results: We find that the model predictions of the H$_2$line intensities for the (1-0)S(0), (1-0)Q(2), (1-0)O(2) and (1-0)O(4) transitions at 2.22, 2.41, 2.63 and 3.00 $\mu$m, respectively, are relatively independent of the astro-chemical model and the gas density distribution when compared against the H$_2$column density, making them robust tracer of the cosmic-ray ionization rate. Conclusions: We recommend the use of ro-vibrational H$_2$line emission in combination with estimation of the cloud's H$_2$column density, to constrain the ionization rate and the spectrum of low energy cosmic-rays.

The existence of possible early oceans in the northern hemisphere of Mars has been researched and debated for decades. The nature of the early martian climate is still somewhat mysterious, but evidence for one or more early oceans implies long-lasting periods of habitability. The primary evidence supporting early oceans is a set of proposed remnant shorelines circling large fractions of the planet. The features are thought to be older than 3.6 Ga and possibly as old as 4 Ga, which would make them some of the oldest large-scale features still identifiable on the surface of Mars. One question that has not been thoroughly addressed, however, is whether shorelines this old could survive modification and destruction processes like impact craters, tectonics, volcanism, and hydrology in recognizable form. Here we address one of these processes -- impact cratering -- in detail. We use standard crater counting age models to generate synthetic, global populations of craters and intersect them with hypothetical shorelines, tracking portions of the shoreline that are directly impacted. The oldest shorelines (>= 4 Ga) are at least 70 % destroyed by direct impacts. Shorelines of any age >3.6 Ga are dissected into relatively short, discontinuous segments no larger than about 40 km when including the effects of craters larger than 100 m in radius. When craters smaller than 500 m in radius are excluded, surviving segment lengths can be as large as ~1000 km. The oldest shorelines exhibit fractal structure after impacts, presenting as a discontinuous collection of lines over a range of scales. If the features are truly shorelines, high-resolution studies should find similar levels of destruction and discontinuity. However, our results indicate that observing shorelines as old as 4 Ga, should they exist, is a significant challenge and raises questions about prior mapping efforts.

Miller Range (MIL) 07687 is a peculiar carbonaceous chondrite officially classified as a CO3. However, it has been found to display unique petrographic properties that are atypical of this group. Moreover, Raman spectra of its polyaromatic carbonaceous matter does not reflect a structural order consistent with the metamorphic history of a type 3 chondrite. As a result, it has been suggested to be an ungrouped C2 chondrite with CO affinities, although it has not been fully excluded as a CO chondrite. The ambiguity of the meteorite classification is the motivation behind the present study. We conclude that MIL 07687 is a unique carbonaceous chondrite with possible affinities to CO, CM and/or some ungrouped carbonaceous chondrites. The difficulty in classifying this meteorite stems from (i) its heavily weathered nature, which interferes with the interpretation of our oxygen (O-)isotopic measurements (ii) the overlap in the petrographic and O-isotopic descriptions of various COs, CMs and ungrouped meteorites in the Meteoritical society database. Optical and infrared spectra are consistent with the meteorite unequilibrated nature and indicate that it is probably mildly aqueously altered. Despite traces of aqueous alteration having previously been described in MIL 07687, this is the first time that the presence of hydrated amorphous silicates is reported. In fact, our results show that its present hydration is beyond that of most CO3s, less than most CM2s, and comparable to primitive CR2s. Consequently, we support the meteorite s C2 ung label, although a CO2 or CM2 classification cannot be fully excluded.

Microlensing is proving to be one of the best techniques to detect distant, low-mass planets around the most common stars in the Galaxy. In principle, Earth's microlensing signal could offer the chance for other technological civilisations to find the Earth across Galactic distances. We consider the photometric microlensing signal of Earth to other potential technological civilisations and dub the regions of our Galaxy from which Earth's photometric microlensing signal is most readily observable as the "Earth Microlensing Zone" (EMZ). The EMZ can be thought of as the microlensing analogue of the Earth Transit Zone (ETZ) from where observers see Earth transit the Sun. Just as for the ETZ, the EMZ could represent a game-theoretic Schelling point for targeted searches for extra-terrestrial intelligence (SETI). To compute the EMZ, we use the Gaia DR2 catalogue with magnitude G<20 to generate Earth microlensing probability and detection rate maps to other observers. Whilst our Solar system is a multi-planet system, we show that Earth's photometric microlensing signature is almost always well approximated by a binary lens assumption. We then show that the Earth is in fact well-hidden to observers with technology comparable to our own. Specifically, even if observers are located around every Gaia DR2 star with G<20, we expect photometric microlensing signatures from the Earth to be observable on average only tens per year by any of them. In addition, the EMZs overlap with the ETZ near the Galactic centre which could be the main areas for future SETI searches.

As one of Stage IV space-based telescopes, the China Space Station Telescope (CSST) can perform photometric and spectroscopic surveys simultaneously to efficiently explore the Universe with extreme precision. In this work, we investigate several powerful CSST cosmological probes, including cosmic shear, galaxy-galaxy lensing, photometric and spectroscopic galaxy clustering, and number counts of galaxy clusters, and study the capability of these probes by forecasting the results of joint constraints on the cosmological parameters. By referring to relevant real observational results, we generate mock data and estimate the measured errors based on CSST observational and instrumental designs. To study the systematical effects on the results, we also consider a number of systematics in CSST photometric and spectroscopic surveys, such as the intrinsic alignment, shear calibration uncertainties, photometric redshift uncertainties, galaxy bias, non-linear effects, instrumental effects, etc. The Fisher matrix method is used to derive the constraint results from individual or joint surveys on the cosmological and systematical parameters. We find that the joint constraints by including all these CSST cosmological probes can significantly improve the results from current observations by one order of magnitude at least, which gives $\Omega_m$ and $\sigma_8$ <1% accuracy, and $w_0$ and $w_a$ <5% and 20% accuracies, respectively. This indicates that the CSST photometric and spectroscopic multi-probe surveys could provide powerful tools to explore the Universe and greatly improve the studies of relevant cosmological problems.

We present analysis of comprehensive radio observations of the black hole V404 Cyg during its 2015 outburst. These data represent the best ever coverage of jet production and particle acceleration from any black hole. We report for the first time a clear and near-linear flux-rms correlation in the radio flux densities. Investigation of individual flares reveals in nearly all cases the peak corresponds to the transition from optically thick to thin to synchrotron emission, but an extended phase of particle acceleration is required in contrast to simple impulsive injection models. The largest radio flare is preceded by a phase of optical oscillations and followed one day later by a smaller but optically thin flare, likely due to ejecta interacting with the interstellar medium. Comparing the radio emission to contemporaneous X-ray and optical data, we find that the X-ray and radio measurements are correlated on all timescales from seconds to one day. Correlation with the optical flux densities is weak at short timescales, but becomes significant on timescales greater than a few hours. We evaluate the physical conditions (size, magnetic field and internal energy) associated with 86 individual radio flares, which in turn allows us to place a lower limit on the kinetic feedback over the 15 days of intense activity. If this energy was deposited locally to the source, as implied by the failure to detect jets on angular scales larger than milliarcsec, then we predict that a nova-like shell could have been formed.

For the first time, the helium abundance relative to hydrogen (He/H), which relied on direct measurements of the electron temperature, has been derived in the narrow line regions (NLRs) from a local sample of Seyfert 2 nuclei. In view of this, optical emission line intensities [$3000 \: < \lambda \: < \: 7000$] of 65 local Seyfert 2 nuclei ($z \: < \: 0.2$), taken from Sloan Digital Sky Survey Data Release 15 and additional compilation from the literature, were considered. We used photoionization model grid to derive an Ionization Correction Factor (ICF) for the neutral helium. The application of this ICF indicates that the NLRs of Seyfert 2 present a neutral helium fraction of $\sim$50 per cent in relation to the total helium abundance. We find that Seyfert 2 nuclei present helium abundance ranging from 0.60 to 2.50 times the solar value, while $\sim$85 per cent of the sample present over-solar abundance values. The derived (He/H)-(O/H) abundance relation from the Seyfert 2 is stepper than that of star-forming regions (SFs) and this difference could be due to excess of helium injected into the Interstellar Medium by the winds of Wolf Rayet stars. From a regression to zero metallicity, by using Seyfert 2 estimates combined with SFs estimates, we obtained a primordial helium mass fraction $Y_{\rm p}$=0.2441$\pm$0.0037, a value in good agreement with the one inferred from the temperature fluctuations of the cosmic microwave background by the Planck Collaboration, i.e. $Y_{\rm p}$=0.2471$\pm$0.0003

We present the design and measured performance of an LED module for spatially mapping kinetic inductance detector arrays in the laboratory. Our novel approach uses a multiplexing (MUX) scheme that only requires seven wires to control 480 red LEDs, and the number of LEDs can be scaled up without adding any additional wires. This multiplexing approach relies on active surface mount components that can operate at cryogenic temperatures down to 10 K. Cryogenic tests in liquid nitrogen and inside our cryostat demonstrate that the MUX circuit works at 77 K and 10 K, respectively. The LED module presented here is tailored for millimeter-wave detector modules we are developing, but the approach could be broadly useful for other KID-based detector systems.

We have built a new model for the lens SDSS J1004+4112 including the recently measured time delay of the fourth quasar image. This time delay has a strong influence on the inner mass distribution of the lensing cluster ($\rho \propto r^{-\alpha}$) allowing us to determine $\alpha=1.18^{+0.02(+0.11)}_{-0.03(-0.18)}$ at 68% (95%) confidence level in agreement with hydrodynamical simulations of massive galaxy clusters. We find an offset between the brightest cluster galaxy (BCG) and the dark matter halo of $3.8^{+0.6(+1.4)}_{-0.7(-1.3)}$ kpc at 68% (95%) confidence which is compatible with other galaxy cluster measurements. As an observational challenge, the estimated time delay between the leading image C and the faint (I=24.7) fifth image E is roughly 8 years.

The processes of star formation and feedback, which regulate the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; $\sim$100pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterise the specific physical mechanisms that drive this transition. By applying a robust statistical method to high-resolution CO and narrow-band H$\alpha$ imaging from the PHANGS survey, we systematically measure the evolutionary timeline from molecular clouds to exposed young stellar regions on GMC scales, across an unprecedented sample of 54 star-forming main-sequence galaxies. We find that clouds live for about $1{-}3$ GMC turbulence crossing times ($5{-}30$Myr) and are efficiently dispersed by stellar feedback within $1{-}5$Myr once the star-forming region becomes partially exposed, resulting in integrated star formation efficiencies of $1{-}8$%. These ranges reflect physical galaxy-to-galaxy variation. In order to evaluate whether galactic environment influences GMC evolution, we correlate our measurements with average properties of the GMCs and their local galactic environment. We find several strong correlations that can be physically understood, revealing a quantitative link between galactic-scale environmental properties and the small-scale GMC evolution. Notably, the measured CO-visible cloud lifetimes become shorter with decreasing galaxy mass, most likely due to the increasing presence of CO-dark molecular gas in such environment. Our results represent a first step towards a comprehensive picture of cloud assembly and dispersal, which will further require extension and refinement with tracers of the atomic gas, dust, and deeply-embedded stellar populations.

In this paper we develop a differential astrometric framework that is appropriate for a scanning space satellite such as Gaia. We apply it to the first of the GAREQ fields, which is the fruit of dedicated efforts within the Gaia project focused on measuring the relativistic deflection of light close to Jupiter's limb. This provides a preliminary assessment of: a) the observability of the relativistic deflection of light close to Jupiter, and, b) Gaia's astrometric capabilities under extremely difficult conditions such as those around a bright extended object. Inputs to our differential astrometric model are the CCD transit times as measured by the IDU system, transformed to field angles via AGIS geometric calibrations, and the commanded/nominal spacecraft attitude. [abridged] To illustrate the model we analyze Gaia astrometric measurements after their calibration through the latest cyclic EDR3/DR3 processing of the GAREQ event of Feb 2017. We use observations of the closest bright target star successfully observed several times by Gaia in close proximity to Jupiter and surrounding reference stars brighter than G<13 mag in transits leading up to the time of closest approach and on subsequent transits. The relativistic deflection signal is detected with a S/N of 50 at closest approach by the target star. This signal is the combined effect due to Jupiter and the Sun, mainly dominated by Jupiter's monopole, demonstrating Gaia's scientific performance under extreme observational conditions. It is an unprecedented detection for the following reasons: a) closest ever to Jupiter's limb (~7") in the optical, and, b) highest S/N at any wavelength. Finally, this work sets the stage for investigations into disentangling the relativistic quadrupole deflection due to Jupiter with future Gaia astrometric measurements.

The structure of the solar corona is made of magnetic flux tubes or loops. Due to the lack of contrast with their environment, observing and studying coronal loops in the quiet Sun is extremely difficult. In this work we use a differential emission measure tomographic (DEMT) technique to reconstruct, from a series of EUV images covering an entire solar rotation, the average 3D distribution of the thermal properties of the coronal plasma. By combining the DEMT products with extrapolations of the global coronal magnetic field, we reconstruct coronal loops and obtain the energy input required to keep them at the typical million-degree temperatures of the corona. We statistically study a large number of reconstructed loops for Carrington rotation (CR) 2082 obtaining a series of typical average loops of different lengths. We look for relations between the thermal properties and the lengths of the constructed typical loops and find similar results to those found in a previous work \citep{maccormack_2020}. We also analyze the typical loop properties by comparing them with the zero-dimensional (0D) hydrodynamic model Enthalpy-Based Thermal Evolution of Loops \citep[EBTEL, ][]{klimchuk_2008}. We explore two heating scenarios. In the first one, we apply a constant heating rate assuming that typical loops are in quasi-static equilibrium. In the second scenario we heat the plasma in the loops using short impulsive events. We find that the reconstructed typical loops are overdense with respect to quasi-static equilibrium solutions of the hydrodynamic model. Impulsive heating, on the other hand, reproduces better the observed densities and temperatures for the shorter and approximately semicircular loops.

We spectroscopically detected candidate emission-lines of 8 likely, 17 probable, and 69 possible strong galaxy-galaxy gravitational lens candidates found within the spectra of ~10,000 galaxy targets contained within the completed Mapping of Nearby Galaxies at Apache Point Observatory (MaNGA) survey. This search is based upon the methodology of the Spectroscopic Identification of Lensing Objects (SILO) project, which extends the spectroscopic detection methods of the BOSS Emission-Line Lensing Survey (BELLS) and the Sloan Lens ACS Survey (SLACS). We scanned the co-added residuals that we constructed from stacks of foreground subtracted row-stacked-spectra (RSS) so a sigma-clipping method can be used to reject cosmic-rays and other forms of transients that impact only a small fraction of the combined exposures. We also constructed narrow-band images from the signal-to-noise of the co-added residuals to observe signs of lensed source images. We also use several methods to compute the probable strong lensing regime for each candidate lens to determine which candidate background galaxies may reside sufficiently near the galaxy centre for strong lensing to occur. We present the spectroscopic redshifts within a value-added catalogue (VAC) for data release 17 (DR17) of SDSS-IV. We also present the lens candidates, spectroscopic data, and narrow-band images within a VAC for DR17. High resolution follow-up imaging of these lens candidates are expected to yield a sample of confirmed grade-A lenses with sufficient angular size to probe possible discrepancies between the mass derived from a best-fitting lens model, and the dynamical mass derived from the observed stellar velocities.

We study a series of relations between physical parameters in coronal loops of the quiet Sun reconstructed by combining tomographic techniques and modeling of the coronal magnetic field. We use differential emission measure tomography (DEMT) to determine the three-dimensional distribution of the electron density and temperature in the corona, and we model the magnetic field with a potential-field source-surface (PFSS) extrapolation of a synoptic magnetogram. By tracing the DEMT products along the extrapolated magnetic field lines, we obtain loop-averaged electron density and temperature. Also, loop-integrated energy-related quantities are computed for each closed magnetic field line. We apply the procedure to Carrington rotation 2082, during the activity minimum between Solar Cycles 23 and 24. We find a scaling law between the loop-average density $N$ and loop length $L$, $N_m \sim L^{-0.35}$, but we do not find a significant relation between loop-average temperature and loop length. We confirm though the previously found result that loop-average temperatures at the equatorial latitudes are lower than at higher latitudes. We associate this behavior with the presence at the equatorial latitudes of loops with decreasing temperatures along their length (``down'' loops), which are in general colder than loops with increasing temperatures (``up'' loops). We find that the obtained scalings for quiet-Sun loops do not generally agree with those found in the case of AR loops from previous observational and theoretical studies. We suggest that to better understand the relations found, it is necessary to forward model the reconstructed loops using hydrodynamic codes working under the physical conditions of the quiet-Sun corona.

False vacuum remnants in first-order phase transitions in the early Universe can form compact objects which may constitute dark matter. Such remnants form because particles develop large mass gaps between the two phases and become trapped in the old phase. We focus on remnants generated in a class of models with trapped dark sector particles, trace their development, and determine their ultimate fate. Depending on model and phase transition parameters, the evolutionary endpoint of these remnants can be primordial black holes, Fermi-balls, Q-balls, or thermal balls, and they all have the potential to constitute some portion or the whole of dark matter within a broad mass range. Notably, dark sector thermal balls can remain at high temperatures until the present day and are a new compact dark matter candidate which derives its energy from the thermal energy of internal particles instead of their mass or quantum pressure.

A number of studies based on data collected by the $\textit{Hubble Space Telescope}$ ($\textit{HST}$) GO-13297 program "HST Legacy Survey of Galactic Globular Clusters: Shedding UV Light on Their Populations and Formation'' have investigated the photometric properties of a large sample of Galactic globular clusters and revolutionized our understanding of their stellar populations. In this paper, we expand previous studies by focusing our attention on the stellar clusters' internal kinematics. We computed proper motions for stars in 56 globular and one open clusters by combining the GO-13297 images with archival $\textit{HST}$ data. The astro-photometric catalogs released with this paper represent the most complete and homogeneous collection of proper motions of stars in the cores of stellar clusters to date, and expand the information provided by the current (and future) $\textit{Gaia}$ data releases to much fainter stars and into the crowded central regions. We also census the general kinematic properties of stellar clusters by computing the velocity-dispersion and anisotropy radial profiles of their bright members. We study the dependence on concentration and relaxation time, and derive dynamical distances. Finally, we present an in-depth kinematic analysis of the globular cluster NGC 5904.

The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly-resolved $L^*$ galaxies from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM ($\Delta x \lesssim 1$ kpc $h^{-1}$ and $M \simeq 200M_\odot$). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of $\lesssim5$ kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond $\sim 0.25R_{200}$, or $\sim50$ kpc at $z=0$. Within $\sim0.25R_{200}$, turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived from a global equilibrium, but an emergent equilibrium balancing radially inward and outward forces is obtained when averaging over large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies even when averaging out small-scale variations.

Large scale surface convection on red supergiants (RSGs) can lead to shifts in the photocenter of the star which might be measured by Gaia and used as a new probe of the surface dynamics of these rare but important stars. Unlike brightness variations, photocenter motions would provide information on the physical scale of the convective cells. The signal would be that RSGs show an excess astrometric noise at the level of a few percent of the stellar radius. Unfortunately, we find that the excess astrometric noise level of Gaia EDR3 is roughly an order of magnitude too large to detect the predicted motions and that RSGs have excess astrometric noise indistinguishable from other stars of similar magnitude and parallax. The typical excess astrometric noise steadily decreases with G magnitude (for G<11 mag), so it is crucial to compare stars of similar brightness.

We present an updated catalog of StrayCats (a catalog of NuSTAR stray light observations of X-ray sources) that includes nearly 18 additional months of observations. StrayCats v2 has an added 53 sequence IDs, 106 rows, and 3 new identified stray light (SL) sources in comparison to the original catalog. The total catalog now has 489 unique sequence IDs, 862 entries, and 83 confirmed StrayCats sources. Additionally, we provide new resources for the community to gauge the utility and spectral state of the source in a given observation. We have created long term light curves for each identified SL source using MAXI and Swift/BAT data when available. Further, source extraction regions for 632 identified SL observations were created and are available to the public. In this paper we present an overview of the updated catalog and new resources for each identified StrayCats SL source.

The transmission spectra of exoplanet atmospheres observed with the Hubble Space Telescope in the near-infrared range (1.1-1.65$\mu$m) frequently show evidence for some combination of clouds and hazes. Identification of systematic trends in exoplanet clouds and hazes is potentially important for understanding atmospheric composition and temperature structure. Here we report on the analysis of spectral modulation using the largest uniformly processed sample of exoplanet transit spectra currently available. The spectral retrieval includes the capability to detect and represent atmospheres in which the composition departs from thermochemical equilibrium. Using this catalog, we find that the impact of clouds/hazes, measured in terms of damping of spectral features, follows a trend from mostly cloudy/hazy around 500K to mostly clear around 1650K. For planets at lower temperatures ($<$ 270K), we observe a potential rise towards reduced cloudy/hazy atmospheres. We also find that a partially transparent aerosol component is frequently present and that it is typically vertically distributed throughout the atmospheric column. Our findings also suggest that while clouds and hazes are common in exoplanet atmospheres, the majority of planets have detectable spectral modulation that is sufficient for characterizing their atmospheres. Additionally, the empirical trend that clouds and hazes are minimized at 1650$^{+874}_{-620}$K revealed in our catalog has predictive utility for modelling the performance of large-scale transiting exoplanets survey, such as planned with the Ariel mission. This trend can also be used for making a probability-based forecast of spectral modulation for a given source in the context of future JWST observations.

Statistical anisotropy in the nanohertz-frequency gravitational-wave background (GWB) is expected to be detected by pulsar timing arrays (PTAs) in the near future. By developing a frequentist statistical framework that intrinsically restricts the GWB power to be positive, we establish scaling relations for multipole-dependent anisotropy decision thresholds that are a function of the noise properties, timing baselines, and cadences of the pulsars in a PTA. We verify that $(i)$ a larger number of pulsars, and $(ii)$ factors that lead to lower uncertainty on the cross-correlation measurements between pulsars, lead to a higher overall GWB signal-to-noise ratio, and lower anisotropy decision thresholds with which to reject the null hypothesis of isotropy. Using conservative simulations of realistic NANOGrav datasets, we predict that an anisotropic GWB with angular power $C_{l=1} > 0.3\,C_{l=0}$ may be sufficient to produce tension with isotropy at the $p = 3\times10^{-3}$ ($\sim3\sigma$) level in near-future NANOGrav data with a $20$~yr baseline. We present ready-to-use scaling relationships that can map these thresholds to any number of pulsars, configuration of pulsar noise properties, and sky coverage. We discuss how PTAs can improve the detection prospects for anisotropy, as well as how our methods can be adapted for more versatile searches.

We present MeerKAT Absorption Line Survey (MALS) observations of the HI gas in the Klemola31 galaxy group ($z=0.029$), located along the line of sight to the radio-loud quasar PKS2020-370 ($z=1.048$). Four galaxies of the group are detected in HI emission, and HI absorption is also detected in front of PKS2020-370 in Klemola31A. The emission and absorption are somewhat compensating on the line of sight of the quasar, and the derived column density of the absorption appears under-estimated, with respect to the neighbouring emission. A symmetric tilted-ring model of Klemola31A, assuming the absorbing gas in regular rotation in the plane, yields a rather high spin temperature of 530K. An alternative interpretation is that the absorbing gas is extra-planar, which will also account for its non-circular motion. The NaI/CaII ratio also suggests that the absorbing gas is unrelated to cold HI disk. Two of the galaxies in the Klemola group are interacting with a small companion, and reveal typical tidal tails, and velocity perturbations. Only one of the galaxies, ESO400-13, reveals a strong HI deficiency, and a characteristic ram-pressure stripping, with a total asymmetry in the distribution of its gas. Since a small galaxy group as Klemola31 is not expected to host a dense intra-group gas, this galaxy must be crossing the group at a very high velocity, mostly in the sky plane.

We report the rotational lightcurves of 21 trans-Neptunian objects (TNOs) in Neptune's 2:1 mean motion resonance obtained with the 6.5 m Magellan-Baade telescope and the 4.3 m Lowell Discovery Telescope. The main survey's goal is to find objects displaying a large lightcurve amplitude which is indicative of contact binaries or highly elongated objects. In our sample, two 2:1 resonant TNOs showed a significant short-term lightcurve amplitude: 2002 VD$_{130}$ and (531074) 2012 DX$_{98}$. The full lightcurve of 2012 DX$_{98}$ infers a periodicity of 20.80$\pm$0.06h and amplitude of 0.56$\pm$0.03mag whereas 2002 VD$_{130}$ rotates in 9.85$\pm$0.07h with a 0.31$\pm$0.04mag lightcurve amplitude. Based on lightcurve morphology, we classify (531074) 2012 DX$_{98}$ as a likely contact binary, but 2002 VD$_{130}$ as a likely single elongated object. Based on our sample and the lightcurves reported in the literature, we estimate the lower percentage of nearly equal-sized contact binaries at only 7-14$\%$ in the 2:1 resonance, which is comparable to the low fraction reported for the dynamically Cold Classical trans-Neptunian objects. This low contact binary fraction in the 2:1 Neptune resonance is consistent with the lower estimate of the recent numerical modeling. We report the Sloan g', r', i' surface colors of 2002 VD$_{130}$ which is an ultra-red TNO whereas 2012 DX$_{98}$ is a very red object based on published surface colors.

Type Ia supernovae (SNe Ia) are powerful tools for measuring the expansion history of the universe, but the impact of dust around SNe Ia remains unknown and is a critical systematic uncertainty. One way to improve our empirical description of dust is to analyse highly reddened SNe Ia ($E(B-V)>0.4$, roughly equivalent to the fitted SALT2 light-curve parameter $c>0.3$). With the recently released Pantheon+ sample, there are 57 SNe Ia that were removed because of their high colour alone (with colours up to $c=1.61$), which can provide enormous leverage on understanding line-of-sight $R_V$. Previous studies have claimed that $R_V$ decreases with colour, though it is unclear if this is due to limited statistics, selection effects, or an alternative explanation. To test this claim, we fit two separate colour-luminosity relationships, one for the main cosmological sample ($c<0.3$) and one for highly reddened ($c>0.3$) SNe Ia. We find the change in the colour-luminosity coefficient to be consistent with zero. Additionally, we compare the data to simulations with different colour models, and find that the data prefers a model with a flat dependence of $R_V$ on colour over a declining dependence. Finally, our results strongly support that line-of-sight $R_V$ to SNe Ia is not a single value, but forms a distribution.

Meteorite matrices from primitive chondrites are an interplay of ingredients at the sub-micron scale, which requires analytical techniques with the nanometer spatial resolution to decipher the composition of individual components in their petrographic context. Infrared spectroscopy is an effective method that enables to probe of vibrations at the molecule-atomic scale of organic and inorganic compounds but is often limited to a few micrometers in spatial resolution. To efficiently distinguish spectral signatures of the different constituents, we apply here nano-IR spectroscopy (AFM-IR), based on the combination of infrared and atomic force microscopy, having a spatial resolution beyond the diffraction limits. Our study aims to characterize two chosen meteorite samples to investigate primitive material in terms of bulk chemistry (the CI chondrite Orgueil) and organic composition (the CR chondrite EET 92042). We confirm that this technique allows unmixing the IR signatures of organics and minerals to assess the variability of organic structure within these samples. We report an investigation of the impact of the widely used chemical HF/HCl (Hydrogen Fluoride/Hydrochloric) extraction on the nature of refractory organics (Insoluble Organic Matter, IOM) and provide insights on the mineralogy of meteorites matrices from these two samples by comparing to reference (extra)terrestrial materials. These findings are discussed with a perspective toward understanding the impact of post-accretional aqueous alteration and thermal metamorphism on the composition of chondrites. Last, we highlight that the heterogeneity of organic matter within meteoritic materials extends down to the nanoscale, and by comparison with IOMs, oxygenated chemical groups are not affected by acid extractions.

The growing number of multi-epoch optical and infrared sky surveys are uncovering unprecedented numbers of new variable stars, of an increasing number of types. The short interval between observations in adjacent near infrared filters in the UKIDSS Galactic Plane Survey (UGPS) allows for the discovery of variability on the timescale of minutes. We report on the nature of one such object, through the use of optical spectroscopy, time-series photometry and targeted X-ray observations. We propose that UGPS J194310.32+183851.8 is a magnetic cataclysmic variable star of novel character, probably featuring a longer than average spin period and an orbital period likely to be shorter than the period gap (i.e. P$_{\text{orb}}$<2 hours). We reason that the star is likely a member of the short period Intermediate-Polar subclass that exist below this period boundary, but with the additional feature that system's SED is fainter and redder than other members of the group.

In a black hole X-ray binary (XRB) system, gas accreted from a normal star onto a black hole glows brightly in X-rays. We report on an observation of the XRB Cygnus X-1 (Cyg X-1) by the Imaging X-ray Polarimetry Explorer IXPE) yielding the first highly significant detection of X-ray polarization from an accreting black hole. The electric vector polarization angle aligns with the outflowing radio jet, supporting the hypothesis that the jet is launched from the inner X-ray emitting region. The higher than expected 2-8 keV polarization degree of 4.0+-0.2% implies that the accretion disk is viewed more edge-on than inferred from the orbital parameters. The spectropolarimetric data reveal that the hot X-ray emitting plasma is extended in the plane of the accretion disk rather than along the jet axis.

[abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$\alpha\,330$, using SPHERE/VLT and ALMA respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45\,au distance from the star in addition to two spiral arms with similar radial launching points at around 90\,au. The millimeter observations from ALMA (resolution of 0.06''$\times$0.04'') show mainly an asymmetric ring located at 110\,au distance from the star. In addition to this asymmetry, there are two faint symmetric rings at 60\,au and 200\,au. The $^{12}$CO, $^{13}$CO and C$^{18}$O lines seem to be less abundant in the inner disk (these observations have a resolution of 0.16''$\times$0.11''). The $^{13}$CO peaks at a location similar to the inner ring observed with SPHERE, suggesting that this line is optically thick and traces variations of disk temperature instead of gas surface density variations, while the C$^{18}$O peaks slightly further away at around 60\,au. We compare our observations with hydrodynamical simulations that include gas and dust evolution, and conclude that a 10\,$M_{\rm{Jup}}$ mass planet at 60\,au and in an eccentric orbit ($e=0.1$) can qualitatively explain most of the observed structures. A planet in a circular orbit leads to a much narrower concentration in the millimeter emission, while a planet in a more eccentric orbit leads to a very eccentric cavity as well. In addition, the outer spiral arm launched by the planet changes its pitch angle along the spiral due to the eccentricity and when it interacts with the vortex, potentially appearing in observations as two distinct spirals. Our observations and models show that LkH$\alpha\,330$, is an exciting target to search for (eccentric-) planets while they are still embedded in their parental disk, making it an excellent candidate for studies on planet-disk interaction.

Gravitational-wave (GW) radiation from a coalescing compact binary is a standard siren as the luminosity distance of each event can be directly measured from the amplitude of the signal. One possibility to constrain cosmology using the GW siren is to perform statistical inference on a population of binary black hole (BBH) events. In essence, this statistical method can be viewed as follows. We can modify the shape of the distribution of observed BBH events by changing cosmological parameters until it eventually matches the distribution constructed from an astrophysical population model, thereby allowing us to determine the cosmological parameters. In this work, we derive the Cram\'er-Rao bound for both cosmological parameters and those governing the astrophysical population model from this statistical dark siren method by examining the Fisher information contained in the event distribution. Our study provides analytical insights and enables fast yet accurate estimations of the statistical accuracy of dark siren cosmology. Furthermore, we consider the bias in cosmology due to unmodeled substructures in the merger rate and the mass distribution. We find a $1\%$ deviation in the astrophysical model can lead to a more than $1\%$ error in the Hubble constant. This could limit the accuracy of dark siren cosmology when there are more than $10^4$ BBH events detected.

In a previous paper, we applied the Gaia DR2 catalog to improve the astrometric accuracy of about 1.7 billion objects in Pan-STARRS1 Data Release 2 (PS1 DR2). We report here on further improvements made by utilizing Gaia EDR3 and correcting for the effects of differential chromatic refraction (DCR) in declination. We extend the correction algorithm in Paper 1 by iteratively subtracting color- and declination-dependent PS1/Gaia EDR3 declination residuals. We determine the astrometric improvement for ~440 million reference objects that are point-like and cross-match to Gaia EDR3. For this set of objects, Gaia EDR3 provides a ~3% improvement in PS1 astrometry over Gaia DR2, and DCR corrections provide an additional ~5% improvement. DCR corrections increase substantially for objects observed away from the zenith. DCR corrections lead to an astrometric improvement of ~30% for blue objects (0<g-i<1) that are 50{\deg} away from the zenith. The amplitude of systematic astrometric errors from these effects is substantially reduced to less than 1 mas for objects with PS1 colors in the range 0 < g-i < 4.5, which makes this a useful astrometric reference catalog in fields where there are few Gaia stars. The improved astrometric data will be available through the Mikulski Archive for Space Telescopes PS1 catalog interfaces.

The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R_Terra and 3.0 M_Terra that is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets. To prepare for future studies, we collected light curves of seven new transits observed with the CHEOPS space mission and new radial velocities obtained with MAROON-X/Gemini North and CARMENES/Calar Alto telescopes, together with previously published spectroscopic and photometric data from the two spectrographs and TESS. We also performed interferometric observations with the CHARA Array and new photometric monitoring with a suite of smaller telescopes. From interferometry, we measure a limb-darkened disc angular size of the star Gl 486. Together with a corrected Gaia EDR3 parallax, we obtain a stellar radius. We also measure a stellar rotation period at P_rot ~ 49.9 d, an upper limit to its XUV (5-920 AA) flux with new Hubble/STIS data, and, for the first time, a variety of element abundances (Fe, Mg, Si, V, Sr, Zr, Rb) and C/O ratio. Besides, we impose restrictive constraints on the presence of additional components, either stellar or substellar, in the system. With the input stellar parameters and the radial-velocity and transit data, we determine the radius and mass of the planet Gl 486 b at R_p = 1.343+/0.063 R_Terra and M_p = 3.00+/-0.13 M_Terra. From the planet parameters and the stellar element abundances, we infer the most probable models of planet internal structure and composition, which are consistent with a relatively small metallic core with respect to the Earth, a deep silicate mantle, and a thin volatile upper layer. With all these ingredients, we outline prospects for Gl 486 b atmospheric studies, especially with forthcoming James Webb Space Telescope observations (abridged).

The formation and evolutionary history of M31 are closely related to its dynamical structures, which remain unclear due to its high inclination. Gas kinematics could provide crucial evidence for the existence of a rotating bar in M31. Using the position-velocity diagram of [OIII] and HI, we are able to identify clear sharp velocity jump (shock) features with a typical amplitude over 100 km/s in the central region of M31 (4.6 kpc X 2.3 kpc, or 20 arcmin X 10 arcmin). We also simulate gas morphology and kinematics in barred M31 potentials and find that the bar-induced shocks can produce velocity jumps similar to those in [OIII]. The identified shock features in both [OIII] and HI are broadly consistent, and they are found mainly on the leading sides of the bar/bulge, following a hallmark pattern expected from the bar-driven gas inflow. Shock features on the far side of the disk are clearer than those on the near side, possibly due to limited data coverage on the near side, as well as obscuration by the warped gas and dust layers. Further hydrodynamical simulations with more sophisticated physics are desired to fully understand the observed gas features and to better constrain the parameters of the bar in M31.

Recently, a new class of white dwarfs (``slowly cooling WDs'') has been identified in the globular cluster M13. The cooling time of these stars is increased by stable thermonuclear hydrogen burning in their residual envelope. These WDs are thought to be originated by horizontal branch (HB) stars populating the HB blue tail, which skipped the asymptotic giant branch phase. To further explore this phenomenon, we took advantage of deep photometric data acquired with the Hubble Space Telescope in the near-ultraviolet and investigate the bright portion of the WD cooling sequence in NGC 6752, another Galactic globular cluster with metallicity, age and HB morphology similar to M13. The normalized WD luminosity function derived in NGC 6752 turns out to be impressively similar to that observed in M13, in agreement with the fact that the stellar mass distribution along the HB of these two systems is almost identical. As in the case of M13, the comparison with theoretical predictions is consistent with $\sim 70\%$ of the investigated WDs evolving at slower rates than standard, purely cooling WDs. Thanks to its relatively short distance from Earth, NGC 6752 photometry reaches a luminosity one order of a magnitude fainter than the case of M13, allowing us to sample a regime where the cooling time delay, with respect to standard WD models, reaches $\sim 300$ Myr. The results presented in this paper provide new evidence for the existence of slowly cooling WDs and further support to the scenario proposing a direct causal connection between this phenomenon and the horizontal branch morphology of the host stellar cluster.

We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation analysis, we measure the periodicity and width of the micro-structure in these three pulsars. The detected micro-structure from PSRs J1022+1001 and J1744-1134 is often highly linearly polarised. In PSR J1022+1001, the linear polarisation position angles of micro-structure pulses are in general flat with a small degree of variation. Using these results, we further examine the frequency and rotational period dependency of micro-structure properties established in previous work, along with the angular beaming and temporal modulation models that explains the appearance of micro-structure. We also discuss a possible link of micro-structure to the properties of some of the recently discovered fast radio bursts which exhibit a very similar emission morphology.

We present the Waves in Nearby Disk galaxies Survey (WiNDS) consisting of 40 nearby low inclination disk galaxies observed through H$\alpha$ high-resolution Fabry Perot interferometry. WiNDS consists of 12 new galaxy observations and 28 data archived observations obtained from different galaxy surveys. We derive two-dimensional line-of-sight velocity fields that are analyzed to identify the possible presence of vertical velocity flows in the galactic disks of these low-inclination late-type galaxies using velocity residual maps, derived from the subtraction of an axisymmetric rotation model to rotational velocity map. Large and globally coherent flows in the line-of-sight velocity of nearly face-on galaxies can be associated with large vertical displacement of the disk with respect to its mid-plane. Our goal is to characterize how frequent vertical perturbations, such as those observed in the Milky Way, arise in the Local Universe. Our currently available data have allowed us to identify 20$\%$ of WiNDS galaxies with strong velocity perturbations that are consistent with vertically perturbed galactic disks.

A large fraction of X-ray sources in our Galaxy are low-mass X-ray binaries, containing a black hole or a neutron star accreting from a gravitationally bound low-mass ($\leq$1 M$_\odot$) companion star. These systems are among the older population of stars and accreting systems in the Galaxy, and typically have long accretion histories. Low-mass X-ray binaries are categorized into various sub-classes based on their observed properties such as X-ray variability and brightness, nature of the companion star and/or the compact object, and binary configuration. In this Chapter, we review the phenomenology of sub-classes of these systems and summarize observational finding regarding their characteristics, populations, and their distribution in the Galaxy.

This paper is the first in a series of preparing and analyzing spectral and other properties for a database of already discovered Changing-Look Active Galactic Nuclei (CL AGNs). Here, we focus on the spectral fitting and analysis of broad emission lines in an exhaustive sample of 112 CL AGNs available in the literature with existing SDSS/BOSS/eBOSS spectroscopy. Additionally, we have gathered older/newer spectral epochs from all the available SDSS data releases to make the database more complete. We use PyQSOFit and perform a homogeneous spectral decomposition of all of our SDSS spectra and tabulate the AGN continuum and emission-line properties per epoch per source, chronologically. This further allows us to categorize the sources in our sample as Turn-On or Turn-Off and subsequently check for repeated occurrences of such phases. We then estimate the black hole mass (M$_{\rm BH}$) and the Eddington ratio ($\lambda_{\rm Edd}$) per epoch per source where the required parameters are available and well-estimated. We realize the movement of the source in the M$_{\rm BH}$-$\lambda_{\rm Edd}$ plane allowing us to check for systematic changes in the source's fundamental properties. We then track their transition along the optical plane of the Eigenvector 1 (EV1) schema and categorize sources that either stay within the same Population (A or B) or make an inter-population movement as a function of spectral epoch. We also test the Balmer decrement (H$\alpha$/H$\beta$) of a subset of our sample of CL AGNs as a function of time and AGN luminosity.

The most massive protoclusters virialize to become clusters at $z\sim 2$, which is also a critical epoch for the evolution of their member galaxies. XLSSC 122 is a $z=1.98$ galaxy cluster with 37 spectroscopically confirmed members. We aim to characterize their star formation histories and to put them in the context of the cluster accretion history. We measure their photometry in 12 bands and create a PSF-matched catalogue of the cluster members. We employ BAGPIPES to fit star formation histories characterized by exponentially decreasing star-forming rates. Stellar masses, metal and dust contents are also treated as free parameters. The oldest stars in the red-sequence galaxies display a range of ages, from 0.5 Gyr to over $\sim$3 Gyrs. Characteristic times are between $\sim$0.1 and $\sim$0.3 Gyr, and the oldest members present the longest times. Using MultiDark Planck 2 dark matter simulations, we calculate the assembly of XLSSC 122-like haloes, weighted by the age posteriors of the oldest members. We found that 74% of these haloes were less than 10% assembled at the onset of star formation, declining to 67% of haloes when such galaxies had formed 50% of their z=1.98 stellar masses. When 90% of their stellar masses were formed, 75% of the haloes were less than 30% assembled. The star formation histories of the red-sequence galaxies seem consistent with episodes of star formation with short characteristic times. Onset and cessation of star formation in the oldest galaxies are both likely to precede XLSSC 122 virialization.

We present a complementary methodology to constrain the total neutrino mass, $\sum m_{\nu}$, based on the diffusion coefficient of the splashback mass function of dark matter halos. Analyzing the snapshot data from the Massive Neutrino Simulations, we numerically obtain the number densities of distinct halos identified via the SPARTA code as a function of their splashback masses at various redshifts for two different cases of $\sum m_{\nu}=0.0$ eV and $0.1$ eV. Then, we fit the numerical results to the recently developed analytic formula characterized by the diffusion coefficient that quantifies the degree of ambiguity in the identification of the splashback boundaries. Our analysis confirms that the analytic formula works excellently even in the presence of neutrinos and that the decrement of its diffusion coefficient with redshift is well described by a linear fit, $B(z-z_{c})$, in the redshift range of $0.2\le z\le 2$. It turns out that the massive neutrino case yields significantly lower value of $B$ and substantially higher value of $z_{c}$ than the massless neutrino case, which indicates that the higher masses the neutrinos have, the more severely the splashback boundaries become disturbed by the surroundings. Given our result, we conclude that the total neutrino mass can in principle be constrained by measuring how rapidly the diffusion coefficient of the splashback mass function diminishes with redshifts at $z\ge 0.2$. We also discuss the anomalous behavior of the diffusion coefficient found at lower redshifts for both of the $\sum m_{\nu}$ cases, and ascribe it to the fundamental limitation of the SPARTA code at $z\le 0.13$.

We developed two different point spread function (PSF) analysis techniques for discovering wide separation binary asteroids in wide field surveys. We then applied these techniques to images of main belt asteroids in the 4 to 60 km size range captured by Pan-STARRS1. Johnston (2019) lists fewer than 10 known binaries in this size range with separations greater than 10% of the primary's Hill radius, so discovering more wide binary asteroids is crucial for understanding the limits of binary stability and improving our knowledge of asteroid masses. We analyzed each image by: i) comparing the major axis orientation of the asteroid's elliptical PSF to its non-sidereal rate on the sky, and ii) comparing the one-dimensional median profile created by collapsing the image along the asteroid's direction of motion to that of nearby field stars. For both methods, we flagged any results that deviated significantly from the expected measurements of single asteroids, and those targets with the most flags were identified as binary candidates for confirmation with high-acuity imaging.

This paper constructs an analytic description for the late stages of giant planet formation. During this phase of evolution, the planet gains the majority of its final mass through gas accretion at a rapid rate. This work determines the density and velocity fields for material falling onto the central planet and its circumplanetary disk, and finds the corresponding column density of this infalling envelope. We derive a steady-state solution for the surface density of the disk as a function of its viscosity (including the limiting case where no disk accretion occurs). Planetary magnetic fields truncate the inner edge of the disk and determine the boundary conditions for mass accretion onto the planet from both direct infall and from the disk. The properties of the forming planet and its circumplanetary disk are determined, including the luminosity contributions from infall onto the planet and disk surfaces, and from disk viscosity. The radiative signature of the planet formation process is explored using a quasi-spherical treatment of the emergent spectral energy distributions. The analytic solutions developed herein show how the protoplanet properties (envelope density distribution, velocity field, column density, disk surface density, luminosity, and radiative signatures) vary with input parameters (instantaneous mass, orbital location, accretion rate, and planetary magnetic field strength).

We present the spatially resolved absolute brightness of the Fe X, Fe XI and Fe XIV visible coronal emission lines from 1.08 to 3.4 $R_\odot$, observed during the 2019 July 2 total solar eclipse (TSE). The morphology of the corona was typical of solar minimum, with a dipole field dominance showcased by large polar coronal holes and a broad equatorial streamer belt. The Fe XI line is found to be the brightest, followed by Fe X and Fe XIV (in disk $B_\odot$ units). All lines had brightness variations between streamers and coronal holes, where Fe XIV exhibited the largest variation. However, Fe X remained surprisingly uniform with latitude. The Fe line brightnesses are used to infer the relative ionic abundances and line of sight averaged electron temperature ($T_e$) throughout the corona, yielding values from 1.25 - 1.4 MK in coronal holes up to 1.65 MK in the core of streamers. The line brightnesses and inferred $T_e$ values are then quantitatively compared to the PSI Magnetohydrodynamic model prediction for this TSE. The MHD model predicted the Fe lines rather well in general, while the forward modeled line ratios slightly underestimated the observationally inferred $T_e$ within 5 to 10 % averaged over the entire corona. Larger discrepancies in the polar coronal holes may point to insufficient heating and/or other limitations in the approach. These comparisons highlight the importance of TSE observations for constraining models of the corona and solar wind formation.

We analyze spectra of a slipping flare kernel observed during the 2015 June 22 M6.5-class flare by the Interface Region Imaging Spectrograph (IRIS). During the impulsive and peak phases of the flare, loops exhibiting an apparent slipping motion along the ribbons were observed in the 131{\AA} channel of SDO/AIA. The IRIS spectrograph slit observed a portion of the ribbons, including a moving kernel corresponding to a flare loop footpoint in Si IV, C II, and Mg II at a very-high 1 s cadence. The spectra observed in the kernel were mostly redshifted and exhibited pronounced red wings, as typically observed in large flares. However, in a small region in one of the ribbons, the Si IV 1402.77{\AA} line was partially blueshifted, with the corresponding Doppler velocity |v_{D}| exceeding 50 km s$^{-1}$. In the same region, the C II 1334.53{\AA}, 1335.66{\AA} and 1335.71{\AA} lines were weakly blueshifted (|v_{D}| < 20 km s$^{-1}$) and showed pronounced blue wings, which were observed also in the Mg II k 2796.35{\AA} as well as the Mg II triplet 2798.75{\AA} and 2798.82{\AA} lines. Using high-cadence AIA observations we found that the region where the blueshifts occurred corresponds to the accelerating kernel front as it moved through a weak-field region. The IRIS observations with high resolution allowed us to capture the acceleration of the kernel under the slit for the first time. The unique observations of blueshifted chromospheric and TR lines provide new constrains for current models of flares.

We report photometric observations of the hidden symbiotic star SU Lyn in the optical bands. In four nights we detect a weak flickering in U band with amplitude of about 0.05 magnitudes. No intranight variations are found in B, V, g' and r' bands. This is one more indication that the secondary component is a white dwarf accreting at low accretion rate. We also searched for intranight variability of a dozen related object (RR Boo, RT Boo, AM Cyg, AG Peg, BF Cyg, NQ Gem, StHa190, V627 Cas, XX Oph, FS Cet and Y Gem) in which no variability above the observational errors is detected.

Recovering credible cosmological parameter constraints in a weak lensing shear analysis requires an accurate model that can be used to marginalize over nuisance parameters describing potential sources of systematic uncertainty, such as the uncertainties on the sample redshift distribution $n(z)$. Due to the challenge of running Markov Chain Monte-Carlo (MCMC) in the high dimensional parameter spaces in which the $n(z)$ uncertainties may be parameterized, it is common practice to simplify the $n(z)$ parameterization or combine MCMC chains that each have a fixed $n(z)$ resampled from the $n(z)$ uncertainties. In this work, we propose a statistically-principled Bayesian resampling approach for marginalizing over the $n(z)$ uncertainty using multiple MCMC chains. We self-consistently compare the new method to existing ones from the literature in the context of a forecasted cosmic shear analysis for the HSC three-year shape catalog, and find that these methods recover similar cosmological parameter constraints, implying that using the most computationally efficient of the approaches is appropriate. However, we find that for datasets with the constraining power of the full HSC survey dataset (and, by implication, those upcoming surveys with even tighter constraints), the choice of method for marginalizing over $n(z)$ uncertainty among the several methods from the literature may significantly impact the statistical uncertainties on cosmological parameters, and a careful model selection is needed to ensure credible parameter intervals.

We present an in-depth study of the reflectance spectra of 39 equilibrated and 41 unequilibrated ordinary chondrites. We demonstrate that consistent measuring conditions are vital for the direct comparison of spectral features between chondrites, otherwise hampering any conclusions. We include a comparison with a total of 466 S-type asteroid reflectance spectra. We analyze (i) if a difference between EOCs and UOCs as well as between H, L and LL can be seen, (ii) if it is possible to identify unequilibrated and equilibrated S-type asteroid surfaces and (iii) if we can further constrain the match between OCs and S-type asteroids all based on reflectance spectra. We checked the classification of the 31 Antarctic UOCs analyzed in the present work, using petrography and magnetic measurements, and evidenced that 74% of them were misclassified. Reflectance spectra were compared between EOCs and UOCs as well as between H, L and LL chondrites using a set of spectral features including band depths and positions, peak reflectance values, spectral slopes and the Ol/(Ol + Px) ratio. UOCs and EOCs reflectance spectra show no clear-cut dichotomy, but a continuum with some EOCs showing stronger absorption bands and peak reflectance values, while others are comparable to UOCs. We show by the example of 6 EOCs that their band depths decrease with decreasing grain size. Based on reflectance spectra alone, it is thus highly challenging to objectively identify an unequilibrated from an equilibrated S-type surface. There is no clear distinction of the chemical groups: only LL EOCs of petrographic type > 4 can be distinguished from H and L through less deep 2-mic band depths and 1-mic band positions at longer wavelengths. Lastly, an anti-correlation between the diameter of the asteroids and their 1000 nm band depth is found indicating that larger sized S-type asteroids show finer grained surfaces.

Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program FAUST (Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars) to study the methanol line emission towards the [BHB2007] 11 protobinary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has been previously detected. Twelve methanol lines have been detected with upper energies in the range [45-537] K along with one 13CH3OH transition. The methanol emission is compact and encompasses both protostars, separated by only 28 au and presents three velocity components, not spatially resolved by our observations, associated with three different spatial regions, with two of them close to 11B and the third one associated with 11A. A non-LTE radiative transfer analysis of the methanol lines concludes that the gas is hot and dense and highly enriched in methanol with an abundance as high as 1e-5. Using previous continuum data, we show that dust opacity can potentially completely absorb the methanol line emission from the two binary objects. Although we cannot firmly exclude other possibilities, we suggest that the detected hot methanol is resulting from the shocked gas from the incoming filaments streaming towards [BHB2007] 11 A and B, respectively. Higher spatial resolution observations are necessary to confirm this hypothesis.

Born-again planetary nebulae (PNe) allow investigating stellar evolution, dust production, and nebular shocks in human timescales. Here we present an analysis of multi-epoch optical spectroscopic observations of the born-again PN A 58 around V605 Aql, which experienced a very late thermal pulse (VLTP) about a century ago. The H-deficient ejecta has experienced a considerable brightening in the time period considered, from 1996 to 2021, with notable changes also in many emission line ratios. Neither the reduction of the extinction caused by the dilution of the ejecta nor the increase of the ionizing photon flux from the central star seem capable to produce these spectral changes, which are instead attributed to shocks in the bipolar H-poor outflow dissociating molecular material and propagating through the outer nebula.

Holographic dark energy with the Hubble radius as infrared cutoff has been considered as a candidate to explain the late-time cosmic acceleration and it can solve the coincidence problem. In this scenario, a non-zero equation of state is only possible if there is an interaction between dark energy and cold dark matter. In this paper, a set of phenomenological interactions are assumed and a detailed analysis of the possible values of the coupling constants is carried out, however the resulting matter power spectrum and cosmic microwave background temperature and polarization power spectra have a shape very far from the observed ones. These results rule out any value for the free parameters and it seems to indicate that the assumed interacting holographic dark energy with a Hubble-scale cutoff is not viable to explain the accelerated expansion of the Universe, when cosmological data are taken into account.

Aims. We constrain the role of different Type Ia supernovae (SN Ia) sub-classes in the chemical enrichment of the Galaxy by studying the abundances of iron and manganese in Galactic stars. We investigate four different SN Ia sub-classes, including the classical single-degenerate near-Chandrasekhar mass SN Ia, the fainter SN Iax systems associated with He accretion from the companion, as well as two sub-Chandrasekhar mass SN Ia models. The latter include the double-detonation of a white dwarf accreting helium-rich matter and violent white dwarf mergers. Methods. The chemical abundances in Galactic stars are determined using Gaia eDR3 astrometry and photometry, and the publicly released spectra obtained within the Gaia-ESO large spectroscopic survey. Non-local thermodynamic equilibrium (NLTE) models are used in the spectroscopic analysis. The GCE models have been updated to include detailed delay time distributions arising from binary population synthesis simulations and the different SN Ia channels, as well as recent yields for core-collapse supernovae and AGB stars. The data-model comparison is performed using a Markov chain Monte Carlo framework that allows us to explore the entire parameter space allowed by the diversity of explosion mechanisms and the Galactic SN Ia rate, taking into account the uncertainties of the observed data. Results. Comparison of the new data with GCE models suggests that the observations can only be explained if the fraction of sub-Mch SNe in the models varies between 50 % and 73 %. The standard Mch SNe are not the dominant channel, but are still needed to account for the elevated [Mn/Fe] ratio in the thin disc. Our results only weakly depends on the assumptions on AGB and core collapse SN yields, providing a strong evidence that sub-Mch SNe play a major role in the chemical evolution of our Galaxy.

Infall of cold dark matter on a galaxy may result in caustic rings where the particle density is enhanced. They may be searched for as features in the galactic rotation curves. Previous studies suggested the evidence for these caustic rings with universal, that is common for different galaxies, parameters. Here we test this hypothesis with a large independent set of rotation curves by means of an improved statistical method. No evidence for universal caustic rings is found in the new analysis.

Based on the long-term light curves collected from the Catalina Sky Survey (CSS) (from 2005 to 2013) and the All-Sky Automated Survey for Supernovae (ASAS-SN) (from 2014 to 2018), optical quasi-periodic oscillations (QPOs) about 300 days can be well determined in the well-known blazar PKS~2155-304 through four different methods: the generalized Lomb-Scargle periodogram (GLSP) method, the weighted wavelet Z-transform (WWZ) technique, the epoch-folded method and redfit method. The GLSP determined significance level for the periodicity is higher than 99.9999\% based on a false alarm probability. The redfit provided confidence level for the periodicity is higher than 99\% in ASAS-SN light curve, after considering the effects of red noise. Based on continuous autoregressive (CAR) process created artificial light curves, the probability of detecting fake QPOs is lower than 0.8\%. The determined optical periodicity of 300 days from CSS and ASAS-SN light curves is well consistent with the reported optical periodicity in the literature. Moreover, three possible models are discussed to explain the optical QPOs in PKS 2155-304: the relativistic Frame-dragging effect, the binary black hole (BBH) model and the jet precession model.

Outbursts from cataclysmic variables, such as dwarf novae (DNe), are prevalent throughout the galaxy and are known to emit strongly in the ultraviolet (UV). However, the UV emission of DNe has not been studied extensively compared with the optical. Characterising in detail the physical processes responsible for DN outburst behaviour requires further UV data. Here we report, as part of our recently launched Transient UV Objects Project (TUVO), the discovery of TUVO-21acq, a new transient which we detected in the UV using data from the Ultraviolet Optical Telescope aboard Swift. We detected two separate outbursts and used the UV data to constrain source properties, focusing on the amplitudes and timescales of the outbursts. During the first outburst the source increased in brightness by >4.1, >2.4, and >3.5 magnitudes and during the second outburst by >4.4, >3.4, and >3.6 magnitudes in the UVW1, UVM2, and UVW2 bands, respectively. The outburst durations were in the range 6-21 days and 11-46 days, and the recurrence time is <316 days. We additionally obtained an optical spectrum during quiescence with SALT. The spectrum exhibited hydrogen Balmer series and helium emission lines, and a flat overall spectral shape. We confirm the nature of the source as an accreting white dwarf which underwent DN outburst based on photometric and spectroscopic properties. This study serves as a proof-of-concept for the TUVO project strategy, demonstrating its capability of discovering and classifying new, interesting UV transients. We discuss the implications of our findings for our understanding of the physics underlying DN outbursts, in particular with respect to the UV emission. We examine the need for simultaneous UV and optical observations during the onset of DN outbursts in order to help answer remaining questions, for example in the characteristics and implications of the UV delay.

Using NOEMA and ALMA 3mm line scans, we measure spectroscopic redshifts of six new dusty galaxies at 3.5<z<4.2 by solidly detecting [CI](1-0) and CO transitions. The sample was selected from the COSMOS and GOODS-North super-deblended catalogs with FIR photometric redshifts >6, based on template IR spectrum energy distribution (SED) from known submillimeter galaxies at z=4--6. Dust SED analyses explain the photo-z overestimate from seemingly cold dust temperatures (Td) and steep Rayleigh-Jeans (RJ) slopes, providing additional examples of cold dusty galaxies impacted by the Cosmic Microwave Background (CMB) as found in Jin et al. (2019). We thus study the general properties of the enlarged sample of 10 ``cold" dusty galaxies over 3.5<z<6. We conclude that these galaxies are deceivingly cold at the surface but actually warm in their starbursting cores. Several lines of evidence support this scenario: (1) The high infrared surface density and cold Td from optically thin models appear to violate the Stefan-Boltzmann law; (2) the gas masses derived from optically thin dust masses are inconsistent with estimates from dynamics and CI luminosities; (3) the implied high star formation efficiencies would conflict with cold Td; (4) high FIR optical depth is implied even using the lower, optically-thick dust masses. This work confirms the existence of a substantial population of deceivingly cold, compact dusty starburst galaxies at z>~4, together with the severe impact of the CMB on their RJ observables, paving the way for the diagnostics of optically thick dust in the early universe. Conventional gas mass estimates based on RJ dust continuum luminosities implicitly assume an optically thin case, overestimating gas masses by a factor of 2--3 on average in compact dusty star-forming galaxies.

We present the discovery of a wide retrograde moving group in the disk plane of the Milky Way using action-angle coordinates derived from the \textit{Gaia} DR3 catalog. The structure is identified from a sample of its members that are currently almost at the pericenter of their orbit and are passing through the Solar neighborhood. The motions of the stars in this group are highly correlated, indicating that the system is probably not phase mixed. With a width of at least 1.5 kpc and with a probable intrinsic spread in metallicity, this structure is most likely the wide remnant of a tidal stream of a disrupted ancient dwarf galaxy (age $\sim 12$ Gyr, $\langle {\rm [Fe/H]} \rangle \sim -1.74$). The structure presents many similarities (e.g. in energy, angular momentum, metallicity, and eccentricity) with the Sequoia merging event. However, it possesses extremely low vertical action $J_z$ which makes it unique even amongst Sequoia dynamical groups. As the low $J_z$ may be attributable to dynamical friction, we speculate that the these stars may be the remnants of the dense core of the Sequoia progenitor.

We report on the discovery in the Gaia DR3 astrometric and spectroscopic catalog of a new polar stream that is found as an over-density in action space. This structure is unique as it has an extremely large apocenter distance, reaching beyond 100 kpc, and yet is detected as a coherent moving structure in the Solar neighborhood with a width of $\sim 4$ kpc. A sub-sample of these stars that was fortuitously observed by LAMOST has a mean spectroscopic metallicity of $\langle {\rm [Fe/H]}\rangle = -1.60^{+0.15}_{-0.16}$ dex and possesses a resolved metallicity dispersion of $\sigma({\rm [Fe/H]}) = 0.32^{+0.17}_{-0.06}$ dex. The physical width of the stream, the metallicity dispersion and the vertical action spread indicate that the progenitor was a dwarf galaxy. The existence of such a coherent and highly radial structure at their pericenters in the vicinity of the Sun suggests that many other dwarf galaxy fragments may be lurking in the outer halo.

Habitability of an exoplanet is believed to be profoundly affected by activities of the host stars, although the related coronal mass ejections (CMEs) are still rarely detected in solar-like and late-type stars. We here report an observational study on flares of two M-dwarfs triggered by the high-cadence survey performed by the Ground Wide-angle Camera system. In both events, the fast, time-resolved spectroscopy enables us to identify symmetric broad H$\alpha$ emission with not only a nearly zero bulk velocity, but also a large projected maximum velocity as high as $\sim700-800\ \mathrm{km\ s^{-1}}$. This broadening could be resulted from either Stark (pressure) effect or a flaring-associated CME at stellar limb. In the context of the CME scenario, the CME mass is estimated to be $\sim4\times10^{18}$ g and $2\times10^{19}$ g. In addition, our spectral analysis reveals a temporal variation of the line center of the narrow H$\alpha$ emission in both events. The variation amplitudes are at tens of $\mathrm{km\ s^{-1}}$, which could be ascribed to the chromospheric evaporation in one event, and to a binary scenario in the other one. With the total flaring energy determined from our photometric monitor, we show a reinforced trend in which larger the flaring energy, higher the CME mass is.

Young pulsars deviate from a perfectly regular spin-down by two non-deterministic phenomena: impulsive glitches and timing noise. Both phenomena are interesting per se, and may provide insights into the superfluid properties of neutron stars, but they also act as a barrier to high-precision pulsar timing and gravitational wave experiments. We study a minimal stochastic model to describe the spin-down of a multicomponent neutron star, with fluctuations in both the internal and external torques. The power spectral density and timing noise strength of this kind of model can be obtained analytically, and compared with known results from pulsar timing observational campaigns. In particular, the presence of flat regions of the power spectral density can be interpreted as a signature of the presence of internal superfluid components. We also derive the expected scaling of the timing noise strength with the pulsar's rotational parameters (or characteristic age). Therefore, the present framework offers a theoretical guideline to interpret the observed features of timing noise in both single pulsars and across pulsar population.

A wide-frequency radio study of Active Galactic Nuclei (AGN) is crucial to evaluate the intervening radiative mechanisms responsible for the observed emission and relate them with the underlying accretion physics. We present wide frequency (5-45 GHz), high-sensitivity (few microJy/beam), (sub)-kpc JVLA observations of a sample of 30 nearby (0.003 < z < 0.3) AGN detected by INTEGRAL/IBIS at hard-X-ray. We find a high detection fraction of radio emission at all frequencies, i.e. > 95 per cent at 5, 10 and 15 GHz and > 80 per cent at 22 and 45 GHz. Two sources out of 30 remain undetected at our high sensitivities. The nuclear radio morphology is predominantly compact, sometimes accompanied by extended jet-like structures, or more complex features. The radio Spectral Energy Distributions (SEDs) of the radio cores appear either as a single or as a broken power-law, a minority of them exhibits a peaked component. The spectral slopes are either flat/inverted or steep, up to a break/peak or over the whole range. The sample mean SED shows a flat slope up to 15 GHz which steepens between 15 and 22 GHz and becomes again flat above 22 GHz. Significant radio-X-ray correlations are observed at all frequencies. About half of the sample features extended emission, clearly resolved by the JVLA, indicating low-power jets or large scale outflows. The unresolved cores, which often dominate the radio power, may be of jet, outflow, and/or coronal origin, depending on the observed frequency.

The formation of hydrides by gas-phase reactions between H2 and a heavy element atom is a very selective process. Reactions with ground-state neutral carbon, oxygen, nitrogen, and sulfur atoms are very endoergic and have high energy barriers because the H2 molecule has to be fragmented before a hydride bond is formed. In cold interstellar clouds, these barriers exclude the formation of CH, OH, NH, and SH radicals through hydrogen abstraction reactions. Here we study a very energetically unfavorable process, the reaction of N(4S) atoms with H2 molecules. We calculated the reaction rate coefficient for H2 in different vibrational levels, using quantum methods for v=0-7 and quasi-classical methods up to v=12. Owing to the high energy barrier, these rate coefficients increase with v and also with the gas temperature. We implemented the new rates in the Meudon PDR code and studied their effect on models with different ultraviolet (UV) illumination conditions. In strongly UV-irradiated dense gas (Orion Bar conditions), the presence of H2 in highly vibrationally excited levels (v>7) enhances the NH abundance by two orders of magnitude (at the PDR surface) compared to models that use the thermal rate coefficient for reaction N(4S) + H2 -> NH + H. The increase in NH column density across the PDR is a factor of ~25. We explore existing Herschel/HIFI observations of the Orion Bar and Horsehead PDRs. We report a 3-sigma emission feature at the ~974 GHz frequency of the NH N_J=1_2-0_1 line toward the Bar. The emission level implies N(NH)~10^13 cm^-2, which is consistent with PDR models using the new rate coefficients for reactions between N and UV-pumped H2. This formation route dominates over hydrogenation reactions involving the less abundant N+ ion. JWST observations will soon quantify the amount and reactivity of UV-pumped H2 in many interstellar and circumstellar environments.

In the pursuit of primordial non-Gaussianities, we hope to access smaller scales across larger comoving volumes. At low redshift, the search for primordial non-Gaussianities is hindered by gravitational collapse, to which we often associate a scale $k_{\rm NL}$. Beyond these scales, it will be hard to reconstruct the modes sensitive to the primordial distribution. When forecasting future constraints on the amplitude of primordial non-Gaussianity, $f_{\rm NL}$, off-diagonal components are usually neglected in the covariance because these are small compared to the diagonal. We show that the induced non-Gaussian off-diagonal components in the covariance degrade forecast constraints on primordial non-Gaussianity, even when all modes are well within what is usually considered the linear regime. As a testing ground, we examine the effects of these off-diagonal components on the constraining power of the matter bispectrum on $f_{\rm NL}$ as a function of $k_{\rm max}$ and redshift, confirming our results against N-body simulations out to redshift $z=10$. We then consider these effects on the hydrogen bispectrum as observed from a PUMA-like 21-cm intensity mapping survey at redshifts $2<z<6$ and show that not including off-diagonal covariance over-predicts the constraining power on $f_{\rm NL}$ by up to a factor of $5$. For future surveys targeting even higher redshifts, such as Cosmic Dawn and the Dark Ages, which are considered ultimate surveys for primordial non-Gaussianity, we predict that non-Gaussian covariance would severely limit prospects to constrain $f_{\rm NL}$ from the bispectrum.

The astronomical detection of formamide (NH$_2$CHO) toward various star-forming regions and in cometary material implies that the simplest amide might have an early origin in dark molecular clouds at low temperatures. Laboratory studies have proven the efficient NH$_2$CHO formation in interstellar CO:NH$_3$ ice analogs upon energetic processing. However, it is still under debate, whether the proposed radical-radical recombination reactions forming complex organic molecules remain valid in an abundant H$_2$O environment. The aim of this work was to investigate the formation of NH$_2$CHO in H$_2$O- and CO-rich ices under conditions prevailing in molecular clouds. Therefore, different ice mixtures composed of H$_2$O:CO:NH$_3$ (10:5:1), CO:NH$_3$ (4:1), and CO:NH$_3$ (0.6:1) were exposed to vacuum ultraviolet photons in an ultra-high vacuum chamber at 10 K. Fourier-transform infrared spectroscopy was utilized to monitor in situ the initial and newly formed species as a function of photon fluence. The infrared spectral identifications are complementarily secured by a temperature-programmed desorption experiment combined with a quadrupole mass spectrometer. The energetic processing of CO:NH$_3$ ice mixtures mainly leads to the NH$_2$CHO formation, along with its chemical derivatives such as isocyanic acid (HNCO) and cyanate ion (OCN$^-$). The formation kinetics of NH$_2$CHO shows an explicit dependency on ice ratios and compositions; the highest yield is found in H$_2$O-rich ice. The astronomical relevance of the resulting reaction network is discussed.

Black hole based tests of general relativity have proliferated in recent times with new and improved detectors and telescopes. Modelling of the black hole neighborhood, where most of the radiation carrying strong-field signature originates, is of utmost importance for robust and accurate constraints on possible violations of general relativity. As a first step, this paper presents the extension of general relativistic magnetohydrodynamic simulations of thin accretion disks to parametrically deformed black holes that generalize the Kerr solution. The extension is based on \textsc{harmpi}, a publicly available member of the \textsc{harm} family of codes, and uses a phenomenological metric to study parametric deviations away from Kerr. The extended model is used to study the disk structure, stability, and radiative efficiency. We also compute the Fe K$\alpha$ profiles in simplified scenarios and present an outlook for the future.

Supernovae feedback driven expansion has proven to be a viable mechanism to explain the average properties of Ultra Diffuse Galaxies (UDGs) such as the sizes, colors, mass and internal kinematics. Here, we explore the origin of stellar metallicity gradients in feedback driven simulated UDGs from the NIHAO project and compare them with the observed distribution of metallicity gradients of both Local Group dwarfs as well as of the recently observed UDG DF44. Simulated UDGs display a large variety of metallicity profiles, showing flat to negative gradients, similarly to what is observed in LG dwarfs, while DF44 data suggest a flat to positive gradient. The variety of metallicity gradients in simulations is set by the interplay between the radius at which star formation occurs and the subsequent supernovae feedback driven stellar redistribution: rotation supported systems tend to have flat metallicity profiles while dispersion supported galaxies show negative and steep profiles. Our results suggest that UDGs are not peculiar in what regards their metallicity gradients, when compared to regular dwarfs. Desirably, a larger observational sample of UDGs' gradients shall be available in the future, in order to test our predictions.

We present the analysis of Multi-Espectr\'ografo en GTC de Alta Resoluci\'on para Astronom\'ia (MEGARA) high-dispersion integral field spectroscopic observations of the bipolar planetary nebula (PN) M 3-38. These observations unveil the presence of a fast outflow aligned with the symmetry axis of M 3-38 that expands with a velocity up to $\pm$225 km s$^{-1}$. The deprojected space velocity of this feature can be estimated to be $\approx$320$^{+130}_{-60}$ km s$^{-1}$, which together with its highly collimated morphology suggests that it is one of the fastest jet detected in a PN. We have also used Kepler observations of the central star of M 3-38 to unveil variability associated with a dominant period of 17.7 days. We attribute this to the presence of a low-mass star with an orbital separation of $\approx$0.12-0.16 au. The fast and collimated ejection and the close binary system point towards a common envelope formation scenario for M 3-38.

Roche lobe overflow from a donor star onto a black hole or neutron star binary companion can evolve to a phase of unstable runaway mass-transfer, lasting as short as hundreds of orbits ($\lesssim 10^{2}$ yr for a giant donor), and eventually culminating in a common envelope event. The highly super-Eddington accretion rates achieved during this brief phase ($\dot{M} \gtrsim 10^{5}\dot{M}_{\rm Edd})$ are accompanied by intense mass-loss in disk winds, analogous but even more extreme than ultra-luminous X-ray (ULX) sources in the nearby universe. Also in analogy with observed ULX, this expanding outflow will inflate an energetic `bubble' of plasma into the circumbinary medium. Embedded within this bubble is a nebula of relativistic electrons heated at the termination shock of the faster $v \gtrsim 0.1 c$ wind/jet from the inner accretion flow. We present a time-dependent, one-zone model for the synchrotron radio emission and other observable properties of such ULX `hyper-nebulae'. If ULX jets are sources of repeating fast radio bursts (FRB), as recently proposed, such hyper-nebulae could generate persistent radio emission and contribute large and time-variable rotation measure to the bursts, consistent with those seen from FRB 20121102 and FRB 190520B. ULX hyper-nebulae can be discovered independent of an FRB association in radio surveys such as VLASS, as off-nuclear point-sources whose fluxes can evolve significantly on timescales as short as years, possibly presaging energetic transients from common envelope mergers.

We combine spectroscopic, photometric, and astrometric information from APOGEE data release 17 and $Gaia$ early data release 3 to perform a self-consistent characterization of $Gaia$-Sausage/Enceladus (GSE), the remnant of the last major merger experienced by the Milky Way, considering stars and globular clusters (GCs) altogether. Our novel set of chemodynamical criteria to select genuine stars of GSE yields a metallicity distribution function with a median [Fe/H] of $-1.22$ dex and $0.23$ dex dispersion. Stars from GSE present an excess of [Al/Fe] and [Mg/Mn] (also [Mg/Fe]) in comparison to surviving Milky Way dwarf satellites, which can be explained by differences in star-formation efficiencies and timescales between these systems. However, stars from Sequoia, another proposed accreted halo substructure, essentially overlap the GSE footprint in all analyzed chemical-abundance spaces, but present lower metallicities. Among probable GCs of GSE with APOGEE observations available, we find no evidence for atypical [Fe/H] spreads with the exception of $\omega$ Centauri ($\omega$Cen). Under the assumption that $\omega$Cen is a stripped nuclear star cluster, we estimate the stellar mass of its progenitor to be $M_\star \approx 1.3 \times 10^9 M_\odot$, well-within literature expectations for GSE. This leads us to envision GSE as the best available candidate for the original host galaxy of $\omega$Cen. We also take advantage of $Gaia$'s photometry and APOGEE metallicities as priors to determine fundamental parameters for eight high-probability ($>$70%) GC members of GSE via statistical isochrone fitting. Finally, the newly determined ages and APOGEE [Fe/H] values are utilized to model the age-metallicity relation of GSE.

This chapter presents deep neural network based methods for enhancing the sensitivity of X-ray telescopic observations with imaging polarimeters. Deep neural networks can be used to determine photoelectron emission directions, photon absorptions points, and photon energies from 2D photoelectron track images, with estimates for both the statistical and model uncertainties. Deep neural network predictive uncertainties can be incorporated into a weighted maximum likelihood to estimate source polarization parameters. Events converting outside of the fiducial gas volume, whose tracks have little polarization sensitivity, complicate polarization estimation. Deep neural network based classifiers can be used to select against these events to improve energy resolution and polarization sensitivity. The performance of deep neural network methods is compared against standard data analysis methods, revealing a < 0.75x improvement in minimum detectable polarization for IXPE-specific simulations. Potential future developments and improvements to these methods are discussed.

We review the multiple population (MP) phenomenon of globular clusters (GCs): i.e., the evidence that GCs typically host groups of stars with different elemental abundances and/or distinct sequences in photometric diagrams. Most Galactic and extragalactic clusters exhibit internal variations of He, C, N, O, Na, and Al. They host two distinct stellar populations: the first population of stars, which resemble field stars with similar metallicities, and one or more second stellar populations that show the signature of high-temperature H-burning. In addition, a sub-sample of clusters hosts stellar populations with different heavy-element abundances. The MP origin remains one of the most puzzling, open issues of stellar astrophysics. We summarize the scenarios for the MP formation and depict the modern picture of GCs and their stellar populations along with the main evolutionary phases. We show that the MP behavior dramatically changes from one cluster to another and investigate their complexity to define common properties. We investigate relations with the host galaxy, the parameters of the host clusters (e.g., GC's mass, age, orbit), and stellar mass. We summarize results on spatial distribution and internal kinematics of MPs. Finally, we review the relation between MPs and the so-called second-parameter problem of the horizontal-branch morphology of GCs and summarize the main findings on the extended main-sequence phenomenon in young clusters.

In recent years, an increasing amount of attention is being paid to the gravitational few-body problem and its applications to astrophysical scenarios. Among the main reasons for this renewed interest there is large number of newly discovered exoplanets and the detection of gravitational waves. Here, we present two numerical codes to model three- and few-body systems, called TSUNAMI and OKINAMI. The TSUNAMI code is a direct few-body code with algorithmic regularization, tidal forces and post-Newtonian corrections. OKINAMI is a secular, double-averaged code for stable hierarchical triples. We describe the main methods implemented in our codes, and review our recent results and applications to gravitational-wave astronomy, planetary science and statistical escape theories.

We investigate the out-of-equilibrium production of scalar dark matter (DM) from the inflaton condensate during inflation and reheating. We assume that this scalar couples only to the inflaton via a direct quartic coupling and is minimally coupled to gravity. We consider all possible production regimes: purely gravitational, weak direct coupling (perturbative), and strong direct coupling (non-perturbative). For each regime, we use different approaches to determine the dark matter phase space distribution and the corresponding relic abundance. For the purely gravitational regime, scalar dark matter quanta are copiously excited during inflation resulting in an infrared (IR) dominated distribution function and a relic abundance which overcloses the universe for a reheating temperature $T_\text{reh}>34 ~\text{GeV}$. A non-vanishing direct coupling induces an effective DM mass and suppresses the large IR modes in favor of ultraviolet (UV) modes and a minimal scalar abundance is generated when the interference between the direct and gravitational couplings is maximal. For large direct couplings, backreaction on the inflaton condensate is accounted for by using the Hartree approximation and lattice simulation techniques. Since scalar DM candidates can behave as non-cold dark matter, we estimate the impact of such species on the matter power spectrum and derive the corresponding constraints from the Lyman-$\alpha$ measurements. We find that they correspond to a lower bound on the DM mass of $\gtrsim 3\times 10^{-4} \, \rm{eV}$ for purely gravitational production, and $\gtrsim 20 \, \rm {eV}$ for direct coupling production. We discuss the implications of these results.

The detection of very high-energy neutrinos by IceCube experiment supports the existence of a comparable gamma-ray counterpart from the same cosmic accelerators. Under the likely assumption that the sources of these particles are of extragalactic origin, the emitted photon flux would be significantly absorbed during its propagation over cosmic distances. However, in the presence of photon mixing with ultra-light axion-like-particles (ALPs), this expectation would be strongly modified. Notably, photon-ALP conversions in the host galaxy would produce an ALP flux which propagates unimpeded in the extragalactic space. Then, the back-conversion of ALPs in the Galactic magnetic field leads to a diffuse high-energy photon flux. In this context, the recent detection of the diffuse high-energy photon flux by the Large High Altitude Air Shower Observatory (LHAASO) allows us to exclude at the $95\%$ CL an ALP-photon coupling $g_{a\gamma}\gtrsim 3.9-7.8 \times 10^{-11}~\mathrm{GeV^{-1}}$ for $m_{a}\lesssim 4\times10^{-7}~\mathrm{eV}$, depending on the assumptions on the magnetic fields and on the original gamma-ray spectrum. This new bound is complementary with other ALP constraints from very-high-energy gamma-ray experiments and sensitivities of future experiments.

We study inflation and the generation of primordial fluctuations in $f(R)$ gravity's rainbow. We calculate the cosmological perturbations and then the scalar and tensor primordial power spectrum. We contrast the predictions of the model with the current observational data from PLANCK and BICEP/Keck. Particularly, we found new results for the scalar spectral index $n_s$ and the tensor-to-scalar ratio $r$ along with new observational constraints on the rainbow functions.

The microwave SQUID multiplexer (umux) has enabled higher bandwidth or higher channel counts across a wide range of experiments in particle physics, astronomy, and spectroscopy. The large multiplexing factor coupled with recent commercial availability of microwave components and warm electronics readout systems make it an attractive candidate for systems requiring large cryogenic detector counts. Since the multiplexer is considered for both bolometric and calorimetric applications across several orders of magnitude of signal frequencies, understanding the bandwidth of the device and its interaction with readout electronics is key to appropriately designing and engineering systems. Here we discuss several important factors contributing to the bandwidth properties of umux systems, including the intrinsic device bandwidth, interactions with warm electronics readout systems, and aliasing. We present simulations and measurements of umux devices coupled with SLAC Microresonator RF (SMuRF) tone-tracking electronics and discuss several implications for future experimental design.

We construct a model for dense matter based on low-density nuclear matter properties that exhibits a chiral phase transition and that includes strangeness through hyperonic degrees of freedom. Empirical constraints from nuclear matter alone allow for various scenarios, from a strong first-order chiral transition at relatively low densities through a weaker transition at higher densities, even up to a smooth crossover not far beyond the edge of the allowed range. The model parameters can be chosen such that at asymptotically large densities the chirally restored phase contains strangeness and the speed of sound approaches the conformal limit, resulting in a high-density phase that resembles deconfined quark matter. Additionally, if the model is required to reproduce sufficiently massive compact stars, the allowed parameter range is significantly narrowed down, resulting for instance in a very narrow range for the poorly known slope parameter of the symmetry energy, $L\simeq (88-92)\, {\rm MeV}$. We also find that for the allowed parameter range strangeness does not appear in the form of hyperons in the chirally broken phase and the chiral transition is of first order. Due to its unified approach and relative simplicity - here we restrict ourselves to zero temperature and the mean-field approximation - the model can be used in the future to study dense matter under compact star conditions in the vicinity of the chiral phase transition, for instance to compute the surface tension or to investigate spatially inhomogeneous phases.

Several recent numerical studies have examined the effects of neutrino neutral-current scattering on fast flavor conversion (FFC). Those studies are in apparent conflict, with some finding enhancement of flavor conversion and others finding suppression. The ones that report enhancement all use homogeneous models, and we discuss in detail the self-consistency issues that they face as a result. We reproduce scattering-enhanced FFC in our own homogeneous calculations, showing that it occurs for both neutral- and charged-current scattering and that it is due to isotropization of the angular distributions over time. Because this is the very feature of the calculations that is not self-consistent, we conclude that the enhancement effect is of unclear astrophysical relevance and may not occur in natural environments.

Simpson and Visser recently proposed a phenomenological way to avoid some kinds of space-time singularities by replacing a parameter whose zero value corresponds to a singularity (say, $r$) with the manifestly nonzero expression $r(u) = \sqrt{u^2 + b^2}$, where $u$ is a new coordinate, and $b =$ \const $>0$. This trick, generically leading to a regular minimum of $r$ beyond a black hole horizon (called a ``black bounce''), may hopefully mimic some expected results of quantum gravity, and was previously applied to regularize the Schwarzschild, Reissner-Nordstrom, Kerr and some other metrics. In this paper it is applied to regularize the Fisher solution with a massless canonical scalar field in general relativity (resulting in a traversable wormhole) and a family of static, spherically symmetric dilatonic black holes (resulting in regular black holes and wormholes). These new regular metrics represent exact solutions of general relativity with a sum of stress-energy tensors of a scalar field with a nonzero self-interaction potential and a magnetic field in the framework of nonlinear electrodynamics with a Lagrangian function $L(F)$, $F = F_{\mu\nu} F^{\mu\nu}$. A novel feature in the present study is that the scalar fields involved have ``trapped ghost'' properties, that is, are phantom in a strong-field region and canonical outside it, with a smooth transition between the regions. It is also shown that any static, spherically symmetric metric can be btained as an exact solution to the Einstein equations with the stress-energy tensor of the above field combination.

In this work, we have made a systematic study of how the gravitational wave frequency of the fundamental mode from compact stars is affected by anisotropic effects using realistic equations of state. Our study is an extension of the seminal research performed by Doneva [Phys. Rev. D 85 (2012) 124023], where a polytropic equation of state was used. To achieve our objective, we considered compact stars which were built by using equations of state in the framework of a relativistic mean field theory for the case of hadronic stars and in the framework of the MIT model for the case of quark stars. In order to obtain some pertinent information that could give us the possibility to detect the anisotropy in compact stars, we also studied and analized the behaviour of various global stellar quantities, e.g., gravitational redshift, stellar mass, radius, among others. We concluded that the anisotropic effects can have important consequences, which are strongly related to the anisotropic parameter and the equation of state of high density matter. Additionally, a comparison with observational data has been made and we have shown that the anisotropic parameter $\lambda$ can be used as a tuning parameter to reproduce mass and radius observational data of neutron stars.

In this work, we investigate the equilibrium configurations of massive white dwarfs (MWD) in the context of modified gravity, namely $f(R,L_m)$ gravity, where $R$ stands for the Ricci scalar and $L_m$ is the Lagrangian matter density. We focused on the specific case $f(R,L_m) = R/2 + L_m + \sigma RL_m$, i.e., we have considered a non-minimal coupling between the gravity field and the matter field, with $\sigma$ being the coupling constant. For the first time, the theory is applied to white dwarfs, in particular to study massive white dwarfs, which is a topic of great interest in the last years. The equilibrium configurations predict maximum masses which are above the Chandrasekhar mass limit. The most important effect of the theory is to increase significantly the mass for stars with radius < 2000 km. We found that the theory can accommodate the super-Chandrasekhar white dwarfs for different star compositions. Apart from this, the theory recovers the General Relativity results for stars with radii larger than 3000 km, independent of the value of $\sigma$.

The strong interaction at low energy scales determines the equation of state (EOS) of supranuclear matters in neutron stars (NSs). It is conjectured that the bulk dense matter may be composed of strangeons, which are quark clusters with nearly equal numbers of $u$, $d$, and $s$ quarks. To characterize the strong-repulsive interaction at short distance and the nonrelativistic nature of strangeons, a phenomenological Lennard-Jones model with two parameters is used to describe the EOS of strangeon stars (SSs). For the first time, we investigate the oscillation modes of non-rotating SSs and obtain their frequencies for various parameterizations of the EOS. We find that the properties of radial oscillations of SSs are different from those of NSs, especially for stars with relatively low central energy densities. Moreover, we calculate the $f$-mode frequency of nonradial oscillations of SSs within the relativistic Cowling approximation. The frequencies of the $f$-mode of SSs are found to be in the range from $6.7\,$kHz to $ 8.7\,\rm{kHz}$. Finally, we study the universal relations between the $f$-mode frequency and global properties of SSs, such as the compactness and the tidal deformability. The results we obtained are relevant to pulsar timing and gravitational waves, and will help to probe NSs' EOSs and infer nonperturbative behaviours in quantum chromodynamics.

We investigate nuclear pasta structures at high temperatures in the framework of relativistic mean field model with Thomas-Fermi approximation. Typical pasta structures (droplet, rod, slab, tube, and bubble) are obtained adopting spherical and cylindrical approximations for Wigner-Seitz cells, where the optimum configurations are fixed by minimizing the free energy. The crystalline configurations comprised of those objects are then examined in a three-dimensional geometry with reflection symmetry. It is found that different crystalline structures can evolve into each other via volume conserving deformations, overcoming rather small barriers. For fixed densities and temperatures, the differences for the free energies per baryon of nuclear pasta in various shapes and lattice structures are typically on the order of tens of keV, suggesting the possible coexistence of those structures. As temperature increases, the thermodynamic fluctuations are expected to disrupt the long-range ordering in nuclear pasta structures. We then estimate the critical conditions for nuclear pasta to become disordered and behave like liquid, which are found to be sensitive to the densities, temperatures, proton fractions, and nuclear shapes. If we further increase temperature, eventually the nonuniform structures of nuclear pasta become unstable and are converted into uniform nuclear matter. The phase diagrams of nuclear matter are then estimated, which should be useful for understanding the evolutions of neutron stars, supernova dynamics, and binary neutron star mergers.

We study the problem of heat conduction in general relativity by using Carter's variational formulation. We write the creation rates of the entropy and the particle as combinations of the vorticities of temperature and chemical potential. We pay attention to the fact that there are two additional degrees of freedom in choosing the relativistic analog of Cattaneo equation for the parts binormal to the caloric and the number flows. Including the contributions from the binormal parts, we find a $\textit{new}$ heat-flow equations and discover their dynamical role in thermodynamic systems. The benefit of introducing the binormal parts is that it allows room for a physical ansatz for describing the whole evolution of the thermodynamic system. Taking advantage of this platform, we propose a proper ansatz that deals with the binormal contributions starting from the physical properties of thermal equilibrium systems. We also consider the stability of a thermodynamic system in a flat background. We find that $\textit{new}$ "Klein" modes exist in addition to the known ones. We also find that the stability requirement is less stringent than those in the literature.

If the spatial sections of the Universe are positively curved, extrapolating the inflationary stage backward in time inevitably leads to a classical bounce. This simple scenario, non-singular and free of exotic physics, deserves to be investigated in details. The background dynamics exhibits interesting features and is shown to be mostly insensitive to initial conditions as long as observational consequences are considered. The primordial scalar power spectrum is explicitly computed, for different inflaton potentials, and the subsequent CMB temperature anisotropies are calculated. The results are compatible with current measurements. Some deviations with respect to the standard paradigm can however appear at large scales and we carefully disentangle what is associated with the vacuum choice with what is more fundamentally due to the bounce itself.

We investigate the equatorial deflection angle of light rays propagating in Kerr-Newman black-bounce spacetime. Furthermore, we analyze the light ray trajectories and derive a closed-form formula for deflection angle in terms of elliptic integrals. The deflection angle increases with the decrease of charge and regularisation parameter for a particular impact parameter. We also study the strong field limit of the deflection angle. Using this strong deflection angle formula and lens equation, we find the radius of the first Einstein ring and study its dependence on the charge and the regularisation parameter. We demonstrate that the charge has a robust effect on the size of the Einstein rings, but the effect of the regularization parameter on the ring size is negligible. These results directly affect the observational appearance of the Kerr-Newman black-bounce.

A formalism to describe the false-vacuum decay in non-perturbative regimes was proposed recently. Here, we extend it to the presence of Einstein gravity and calculate the corresponding effective potential and decay rate for a $\lambda \phi^4$ scalar field theory. A comparison with the usual perturbative decay rate shows that the higher the coupling $\lambda$, the greater the decay probability. From the running of the self-interaction coupling, we conclude that the theory becomes weakly coupled in the infrared limit, which proves that Einstein gravity made the weak coupling approximation even more reliable as the universe cooled down. We comment on future applications of these results to cosmological phase transitions, gravitational-wave astronomy, and condensed matter systems.

The information loss paradox associated with black hole Hawking evaporation is an unresolved problem in modern theoretical physics. In a recent brief essay, we revisited the the evolution of the black hole entanglement entropy via the Euclidean path integral (EPI) of the quantum state and allow for the branching of semi-classical histories along the Lorentzian evolution. We posited that there exist at least two histories that contribute to EPI, where one is an information-losing history while the other is information-preserving. At early times, the former dominates EPI, while at late times the latter becomes dominant. By so doing we recovered the essence of the Page curve and thus the unitarity, albeit with the turning point, i.e., the Page time, much shifted toward the late time. In this full-length paper, we fill in the details of our arguments and calculations to strengthen our notion. One implication of this modified Page curve is that the entropy bound may thus be violated. We comment on the similarity and difference between our approach and that of the replica wormholes and the islands conjectures.

We investigate resonance crossings of a charged body moving around a Kerr black hole immersed in an external homogeneous magnetic field. This system can serve as an electromagnetic analogue of a weakly non-integrable extreme mass ratio inspiral (EMRI). In particular, the presence of the magnetic field renders the conservative part of the system non-integrable in the Liouville sense, while the electromagnetic self-force causes the charged body to inspiral. By studying the conservative dynamics, we show the existence of an approximate Carter-like constant and discuss how resonances grow as a function of the perturbation parameter. Then, we apply the electromagnetic self-force to investigate crossings of these resonances during an inspiral. Averaging the energy and angular momentum losses during crossings allows us to employ an adiabatic approximation for them. We demonstrate that such adiabatic approximation provides results qualitatively equivalent to the instantaneous self-force evolution, which indicates that the adiabatic approximation may describe the resonance crossing with sufficient accuracy in EMRIs.

We demonstrate that the direction of energy cascade in compressible planar turbulent flows, such as those in astrophysical disks, must reverse at higher Mach numbers. This leads to supersonic turbulence which, similarly to its three-dimensional counterpart, generates multifractal density fields. In the presence of weak initial fluctuations of the magnetic field, its horizontal components decay as in the Zeldovich antidynamo theorem for two-dimensional incompressible flows. However the vertical component may reach a strongly fluctuating steady state, if the mean field is non-zero. In such cases, the magnetic field energy is enhanced by a large factor as compared to estimates based on the mean field. We consider possible implications for the problem of accretion.

The electrons are an essential particle species in the solar wind. They often exhibit non-equilibrium features in their velocity distribution function. These include temperature anisotropies, tails (kurtosis), and reflectional asymmetries (skewness), which contribute a significant heat flux to the solar wind. If these non-equilibrium features are sufficiently strong, they drive kinetic micro-instabilities. We develop a semi-graphical framework based on the equations of quasi-linear theory to describe electron-driven instabilities in the solar wind. We apply our framework to resonant instabilities driven by temperature anisotropies. These include the electron whistler anisotropy instability and the propagating electron firehose instability. We then describe resonant instabilities driven by reflectional asymmetries in the electron distribution function. These include the electron/ion-acoustic, kinetic Alfv\'en heat-flux, Langmuir, electron-beam, electron/ion-cyclotron, electron/electron-acoustic, whistler heat-flux, oblique fast-magnetosonic/whistler, lower-hybrid fan, and electron-deficit whistler instability. We briefly comment on non-resonant instabilities driven by electron temperature anisotropies such as the mirror-mode and the non-propagating firehose instability. We conclude our review with a list of open research topics in the field of electron-driven instabilities in the solar wind.