Papers for Thursday, Jun 01 2023

The interplay between dark matter (DM) and baryons has long been ignored when building galaxies semi-empirically and observationally. Here I show that baryonic gravity leads to an adiabatic contraction of DM halos, which is most significant in massive galaxies. Ignoring this effect, the derived DM halos are not guaranteed in dynamical equilibrium. I present a new approach to deriving DM halos from rotation curves, which incorporates the adiabatic contraction. The compressed halos turn out super cuspy with respect to NFW halos, which require smaller baryonic contributions and less concentrated primordial halos. I also examine the semi-empirical approach to building galaxies, and find the adiabatic contraction can shift massive galaxies from the observed radial acceleration relation dramatically. Both approaches lead to super cuspy DM halos for massive galaxies, demonstrating the importance of the baryon-driven contraction, which has to be taken into account in order to make an apple-to-apple comparison with simulations.

Context. AM CVn binaries, that is, systems in which a white dwarf accretes matter from a helium-rich (semi-)degenerate object, are potential progenitors of thermonuclear supernovae and strong sources of persistent gravitational wave radiation. For a long time, it has been believed that these systems cannot descend from cataclysmic variables (CVs), at least not in large numbers, because the initial conditions need to be fine-tuned and, even worse, the resulting surface hydrogen abundance would be high enough to be detected which contradicts a defining feature of AM CVn binaries. Aims. Here we show that both claimed weaknesses of the CV formation channel for AM CVn binaries are model-dependent and rely on poorly constrained assumptions for magnetic braking. Methods. We performed binary evolution simulations with the MESA code for different combinations of post-common-envelope white dwarf and companion masses as well as orbital periods assuming strong magnetic braking. Results. We found that AM CVn binaries with extremely-low surface hydrogen abundances are one natural outcome of CV evolution if the donor star has developed a non-negligible helium core prior to the onset of mass transfer. In this case, after hydrogen envelope exhaustion during CV evolution, the donor becomes degenerate and its surface hydrogen abundance substantially drops and becomes undetectable. Our simulations also show that the CV formation channel is able to explain the observed AM CVn binaries with very low mass and bloated donor stars (Gaia14aae and ZTF J1637+49) which has not been demonstrated for any alternative formation channel. Conclusions. CVs with evolved donors are likely the progenitors of at least a fraction of AM CVn binaries.

We present the redshift evolution of quasar luminosity functions decomposed by halo mass, galaxy mass, supermassive black hole (SMBH) mass, and Eddington ratio, as well as SMBH kinetic/radiative energy output ratios from TRINITY, a flexible empirical model that self-consistently infers the halo--galaxy--SMBH connection that match observational data. Key findings include: 1) The normalization of QLF increases by ~3-4 dex from z~10 to z~4, due to the fast mass build-up of different SMBH populations; 2) From z~4 to z~1, less massive galaxies and SMBHs make up bigger and bigger fractions of QLFs, due to the AGN downsizing effect; 3) At z~0, massive haloes/galaxies/SMBHs are responsible for most bright quasars due to low Eddington ratios among all SMBHs; 4) The bright ends of quasar luminosity functions (QLFs) are dominated by SMBHs that are at least 0.3 dex over-massive relative to the median SMBH mass-galaxy mass relation; 5) QLFs at z~6-7 are dominated by SMBHs accreting at Eddington ratios 0.1 < $\eta_\mathrm{rad}$ < 1, but super-Eddington AGNs dominate QLFs towards z~9-10.

The Galactic bulge and bar are critical to our understanding of the Milky Way. However, due to the lack of reliable stellar distances, the structure and kinematics of the bulge/bar beyond the Galactic center have remained largely unexplored. Here, we present a method to measure distances of luminous red giants using a period-amplitude-luminosity relation anchored to the Large Magellanic Cloud, with random uncertainties of 10-15% and systematic errors below 1-2%. We apply this method to data from the Optical Gravitational Lensing Experiment (OGLE) to measure distances to $190,302$ stars in the Galactic bulge and beyond out to 20 kpc. Using this sample we measure a distance to the Galactic center of $R_0$ = $8108\pm106_{\rm stat}\pm93_{\rm sys}$ pc, consistent with astrometric monitoring of stars orbiting Sgr A*. We cross-match our distance catalog with Gaia DR3 and use the subset of $39,566$ overlapping stars to provide the first constraints on the Milky Way's velocity field ($V_R,V_\phi,V_z$) beyond the Galactic center. We show that the $V_R$ quadrupole from the bar's near side is reflected with respect to the Galactic center, indicating that the bar is both bi-symmetric and aligned with the inner disk, and therefore dynamically settled along its full extent. We also find that the vertical height $V_Z$ map has no major structure in the region of the Galactic bulge, which is inconsistent with a current episode of bar buckling. Finally, we demonstrate with N-body simulations that distance uncertainty plays a major factor in the alignment of the major and kinematic axes of the bar and distribution of velocities, necessitating caution when interpreting results for distant stars.

Solar flares and coronal mass ejections, the primary space weather disturbances affecting the entire heliosphere and near-Earth environment, mainly emanate from sunspot regions harbouring high degrees of magnetic twist. However, it is not clear how magnetic helicity, the quantity for measuring the magnetic twist, is supplied to the upper solar atmosphere via the emergence of magnetic flux from the turbulent convection zone. Here, we report state-of-the-art numerical simulations of magnetic flux emergence from the deep convection zone. By controlling the twist of emerging flux, we find that with the support of convective upflow, the untwisted emerging flux can reach the solar surface without collapsing, in contrast to previous theoretical predictions, and eventually create sunspots. Because of the turbulent twisting of magnetic flux, the produced sunspots exhibit rotation and inject magnetic helicity into the upper atmosphere, amounting to a substantial fraction of injected helicity in the twisted cases that is sufficient to produce flare eruptions. This result indicates that the turbulent convection is responsible for supplying a non-negligible amount of magnetic helicity and potentially contributes to solar flares.

We present two- and three-dimensional hydrodynamic simulations of $\sim$kpc-scale AGN jets with mean jet powers in the range $1-7\times10^{45}\,$erg~s$^{-1}$, in which the jet power varies (through variation of the Lorentz factor) according to a flicker or pink noise power spectrum. We find the morphology and dynamics of the jet-cocoon system depends on the amplitude of the variability with a clear correspondence between the shape of the cocoon and the historical activity. The jet advances quickly during high-power states, whereas quiescent periods instead produce passive periods of inflation resembling Sedov-Taylor blast waves. Periods of high activity preferentially produce hotspots and create stronger backflow as they maximise the pressure gradient between the jet head and cocoon. The variability can also lead to propagating internal shock structures along the jet. Our work suggests that variability and flickering in the jet power has important implications, which we discuss, for observations of radio galaxies, ultrahigh energy cosmic ray acceleration and jet power to luminosity correlations. We explore the link between morphology and fuelling, and suggest that chaotic cold accretion should introduce a relatively small scatter in radio luminosity ($\sim0.2$ dex) and modest imprints on morphology; sources such as Hercules A and Fornax A, which show evidence for more dramatic variability, may therefore require redder power spectra, or be triggered by mergers or other discrete events. We suggest ways to search for jet flickering observationally and propose that radio galaxies may be an important diagnostic of Myr timescale AGN fuelling, due to their `long-term memory'.

[Abridged] We combine HST images from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey with JWST images from the Cosmic Evolution Early Release Science (CEERS) survey to measure the stellar and dust-obscured star formation distributions of a mass-complete ($>10^{10}M_\odot$) sample of 95 star-forming galaxies (SFGs) at $0.1<z<2.5$. Rest-mid-infrared (rest-MIR) morphologies (sizes and S\'ersic indices) are determined using their sharpest Mid-InfraRed Instrument (MIRI) images dominated by dust emission. Rest-MIR S\'ersic indices are only measured for the brightest MIRI sources ($S/N>75$; 38 galaxies). At lower $S/N$, simulations show that simultaneous measurements of the size and S\'ersic index become unreliable. We extend our study to fainter sources ($S/N>10$; 95 galaxies) by fixing their S\'ersic index to unity. The S\'ersic index of bright galaxies ($S/N>75$) has a median value of 0.7, which, together with their axis ratio distribution, suggests a disk-like morphology in the rest-MIR. Galaxies above the main sequence (MS; i.e., starbursts) have rest-MIR sizes that are a factor 2 smaller than their rest-optical sizes. The median rest-optical to rest-MIR size ratio of MS galaxies increases with stellar mass, from 1.1 at $10^{9.8}M_\odot$ to 1.6 at $10^{11}M_\odot$. This mass-dependent trend resembles the one found in the literature between the rest-optical and rest-near-infrared sizes of SFGs, suggesting that it is due to dust attenuation gradients affecting rest-optical sizes and that the sizes of the stellar and star-forming components of SFGs are, on average, consistent at all masses. There is, however, a small population of SFGs (15%) with a compact star-forming component embedded in a larger stellar structure. This could be the missing link between galaxies with an extended stellar component and those with a compact stellar component; the so-called blue nuggets.

Context: Combining high-contrast imaging with medium- or high-resolution integral field spectroscopy has the potential to boost the detection rate of exoplanets, especially at small angular separations. Furthermore, it immediately provides a spectrum of the planet that can be used to characterise its atmosphere. The achievable spectral resolution, wavelength coverage, and FOV of such an instrument are limited by the number of available detector pixels. Methods: The trade-offs are studied through end-to-end simulations of a typical high-contrast imaging instrument, analytical considerations, and atmospheric retrievals. The results are then validated with archival VLT/SINFONI data of the planet beta Pictoris b. Results: We show that molecular absorption spectra generally have decreasing power towards higher spectral resolution and that molecule mapping is already powerful for moderate resolutions (R>300). When choosing between wavelength coverage and spectral resolution for a given number of spectral bins, it is best to first increase the spectral resolution until R~2,000 and then maximise the bandwidth within an observing band. We find that T-type companions are most easily detected in the J/H band through methane and water features, while L-type companions are best observed in the H/K band through water and CO features. Such an instrument does not need to have a large FOV, as most of the gain in contrast is obtained in the speckle-limited regime close to the star. We show that the same conclusions are valid for the constraints on atmospheric parameters such as the C/O ratio, metallicity, surface gravity, and temperature, while higher spectral resolution (R~10,000) is required to constrain the radial velocity and spin of the planet.

Cold Jovian planets play an important role in sculpting the dynamical environment in which inner terrestrial planets form. The core accretion model predicts that giant planets cannot form around low-mass M dwarfs, although this idea has been challenged by recent planet discoveries. Here, we investigate the occurrence rate of giant planets around low-mass (0.1-0.3M$_\odot$) M dwarfs. We monitor a volume-complete, inactive sample of 200 such stars located within 15 parsecs, collecting four high-resolution spectra of each M dwarf over six years and performing intensive follow-up monitoring of two candidate radial-velocity variables. We use TRES on the 1.5 m telescope at the Fred Lawrence Whipple Observatory and CHIRON on the Cerro Tololo Inter-American Observatory 1.5 m telescope for our primary campaign, and MAROON-X on Gemini North for high-precision follow-up. We place a 95%-confidence upper limit of 1.5% (68%-confidence limit of 0.57%) on the occurrence of $M_{\rm P}$sin$i > $1M$_{\rm J}$ giant planets out to the water snow line and provide additional constraints on the giant planet population as a function of $M_{\rm P}$sin$i$ and period. Beyond the snow line ($100$ K $< T_{\rm eq} < 150$ K), we place 95%-confidence upper limits of 1.5%, 1.7%, and 4.4% (68%-confidence limits of 0.58%, 0.66%, and 1.7%) for 3M$_{\rm J} < M_{\rm P}$sin$i < 10$M$_{\rm J}$, 0.8M$_{\rm J} < M_{\rm P}$sin$i < 3$M$_{\rm J}$, and 0.3M$_{\rm J} < M_{\rm P}$sin$i < 0.8$M$_{\rm J}$ giant planets; i.e., Jupiter analogs are rare around low-mass M dwarfs. In contrast, surveys of Sun-like stars have found that their giant planets are most common at these Jupiter-like instellations.

A large fraction of the baryons at low redshift are undetected, and likely reside in the tenuous, hot intergalactic medium (IGM). One way to probe the missing baryons is through their absorption of bright sources. The anomalous absorption excess in the X-ray afterglows of $\gamma$-ray bursts (GRBs) has been suggested to result from the missing baryons. In order to test this hypothesis, the present paper employs IllustrisTNG simulations to compute the X-ray absorption effect on cosmological distances. The simulation shows that ionization of H and He in the IGM leaves the metals responsible for $>60\%$ of the X-ray opacity of high-$z$ sources (the ionization of He isn't available in the simulation, we used here external knowledge). The high-$z$ asymptotic optical depth at 0.5\,keV in the simulation reaches $0.15\pm0.07$, while the GRB afterglow values tend to $\approx 0.4$, implying the missing baryons can account for a significant fraction of the observed opacity. The remaining discrepancy is ascribed mainly to the low average metallicity in the simulation, which drops from 0.06 solar at $z=0$ to 0.01 at $z=3$, and which is below previously measured values.

Unexpectedly strong X-ray emission from extragalactic radio jets on kiloparsec scales has been one of the major discoveries of Chandra, the only X-ray observatory capable of sub-arcsecond-scale imaging. The origin of this X-ray emission, which appears as a second spectral component from that of the radio emission, has been debated for over two decades. The most commonly assumed mechanism is inverse Compton upscattering of the Cosmic Microwave Background (IC-CMB) by very low-energy electrons in a still highly relativistic jet. Under this mechanism, no variability in the X-ray emission is expected. Here we report the detection of X-ray variability in the large-scale jet population, using a novel statistical analysis of 53 jets with multiple Chandra observations. Taken as a population, we find that the distribution of p-values from a Poisson model is strongly inconsistent with steady emission, with a global p-value of 1.96e-4 under a Kolmogorov-Smirnov test against the expected Uniform (0,1) distribution. These results strongly imply that the dominant mechanism of X-ray production in kpc-scale jets is synchrotron emission by a second population of electrons reaching multi-TeV energies. X-ray variability on the time-scale of months to a few years implies extremely small emitting volumes much smaller than the cross-section of the jet.

Through their lifetime sunspots undergo a change in their area and shape and, as they decay, they fragment into smaller structures. Here, for the first time we analyze the spatial structure of magnetohydrodynamic (MHD) slow body and fast surface modes in observed umbrae as their cross-sectional shape changes. The Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) techniques were used to analyze 3 and 6 hours SDO/HMI time series of Doppler velocities at the photospheric level of approximately circular and elliptically shaped sunspots. Each time series were divided equally into time intervals, to evidence the change of the sunspots' shape. To identify physical wave modes, the POD/DMD modes were cross-correlated with a slow body mode model using the exact shape of the umbra, whereas the shape obtained by applying a threshold level of the mean intensity for every time interval. Our results show that the spatial structure of MHD modes are affected, even by apparently small changes of the umbral shape, especially in the case of the higher-order modes. For the datasets used in our study, the optimal time intervals to consider the influence of the change in the shape on the observed MHD modes is 37 - 60 minutes. The choice of these intervals is crucial to properly quantify the energy contribution of each wave mode to the power spectrum.

Traditionally, the solar activity cycle is thought as an interplay of the main dipole component of the solar poloidal magnetic field and the toroidal magnetic field. However, the real picture as presented in the extended solar-cycle models is much more complicated. Here, we develop the concept of the extended solar cycle clarifying what zonal harmonics are responsible for the equatorward and polarward propagating features in the surface activity tracers. We arrive at a conclusion that the zonal harmonics with L = 5 play a crucial role in separating the phenomena of both types, which are associated with the odd zonal harmonics. Another objective of our analysis is the role of even zonal harmonics, which prove to be rather associated with the North-South asymmetry of the solar activity than with its 11-year solar periodicity.

Young massive clusters (YMCs) are the most massive star clusters forming in nearby galaxies and are thought to be a young analogue to the globular clusters. Understanding the formation process of YMCs leads to looking into very efficient star formation in high-redshift galaxies suggested by recent JWST observations. We investigate possible observational signatures of their formation stage, particularly when the mass of a cluster is increasing via accretion from a natal molecular cloud. To this end, we study the broad-band continuum emission from ionized gas and dust enshrouding YMCs, whose formation is followed by recent radiation-hydrodynamics simulations. We perform post-process radiative transfer calculations using simulation snapshots and find characteristic spectral features at radio and far-infrared frequencies. We show that a striking feature is long-lasting, strong free-free emission from a $\sim$ 10pc-scale HII region with a large emission measure of $\gtrsim 10^7 \mathrm{cm}^{-6} \ \mathrm{pc}$, corresponding to the mean electron density of $\gtrsim 10^3~\mathrm{cm}^{-3}$. There is a turnover feature below $\sim$ 10 GHz, a signature of the optically-thick free-free emission, often found in Galactic ultra-compact HII regions. These features come from the peculiar YMC formation process, where the cluster's gravity effectively traps photoionized gas for a long duration and enables continuous star formation within the cluster. Such large and dense HII regions show distinct distribution on the density-size diagram, apart from the standard sequence of Galactic HII regions. This is consistent with the observational trend inferred for extragalactic HII regions associated with YMCs.

We review the theoretical underpinnings, evolutionary status, calibrations and current applications of three bright Population II extragalactic distance indicators: Tip of the Red Giant Branch (TRGB) stars, RR Lyrae variables and J-Branch Asymptotic Giant Branch (JAGB/Carbon) stars. For M_I (TRGB) = -4.05 mag the Hubble constant is determined to be Ho = 69.8 +/- 0.6 (stat) +/-1.6 (sys) km/s/Mpc.

Exoplanet atmospheres are the key to understanding the nature of exoplanets. To this end, transit spectrophotometry provides us opportunities to investigate the physical properties and chemical compositions of exoplanet atmospheres. We aim to detect potential atmospheric signatures in 12 gaseous giant exoplanets using transit spectrophotometry and we try to constrain their atmospheric properties. The targets of interest were observed using transit spectrophotometry with the GTC OSIRIS instrument. We estimated the transit parameters and obtained the optical transmission spectra of the target planets using a Bayesian framework. We analyzed the spectral features in the transmission spectra based on atmospheric retrievals. Most of the observed transmission spectra were found to be featureless, with only the spectrum of CoRoT-1b showing strong evidence for atmospheric features. However, in combination with the previously published near-infrared transmission spectrum, we found multiple interpretations for the atmosphere of CoRoT-1b due to the lack of decisive evidence for alkali metals or optical absorbers. Featureless spectra are not necessarily indicative of cloudy atmospheres if they poorly constrain the altitudes of cloud decks. Precise constraints on the models of hazes and clouds strongly depend on the significance of the observed spectral features. Further investigations on these exoplanets, especially CoRoT-1b, are required to confirm the properties of their atmospheres.

Energetic particle populations in the solar corona and in the heliosphere appear to have different characteristics even when produced in the same solar flare. It is not clear what causes this difference: properties of the acceleration region, the large-scale magnetic field configuration in the flare, or particle transport effects, such as scattering. In this study we use a combination of magnetohydrodynamic and test-particle approaches to investigate magnetic reconnection, particle acceleration and transport in two solar flares: an M-class flare on June 19th, 2013, and an X-class flare on September 6th, 2011. We show that in both events , the same regions are responsible for the acceleration of particles remaining in the coronal and being ejected towards the heliosphere. However, the magnetic field structure around the acceleration region acts as a filter, resulting in different characteristics (such as energy spectra) acquired by these two populations. We argue that this effect is an intrinsic property of particle acceleration in the current layers created by the interchange reconnection and, therefore, may be ubiquitous, particularly, in non-eruptive solar flares with substantial particle emission into the heliosphere.

We present the first detailed chemical analysis from APOGEE-2S observations of stars in six regions of recently discovered substructures in the outskirts of the Magellanic Clouds (MCs) extending to 20 degrees from the Large MC (LMC) center. We also present, for the first time, the metallicity and alpha-abundance radial gradients of the LMC and SMC out to 11 degrees and 6 degrees, respectively. Our chemical tagging includes 13 species including light, alpha, and Fe-peak elements. We find that the abundances of all of these chemical elements in stars populating two regions in the northern periphery - along the northern "stream"-like feature - show good agreement with the chemical patterns of the LMC, and thus likely have an LMC origin. For substructures located in the southern periphery of the LMC, we find more complex chemical and kinematical signatures, indicative of a mix of LMC-like and SMC-like populations. However, the southern region closest to the LMC shows better agreement with the LMC, whereas that closest to the SMC shows a much better agreement with the SMC chemical pattern. When combining this information with 3-D kinematical information for these stars, we conclude that the southern region closest to the LMC has likely an LMC origin, whereas that closest to the SMC has an SMC origin, and the other two southern regions have a mix of LMC and SMC origins. Our results add to the evidence that the southern substructures of the LMC periphery are the product of close interactions between the LMC and SMC, and thus likely hold important clues that can constrain models of their detailed dynamical histories.

Sub-Neptune type exoplanets are abundant in our galaxy yet have no solar system analogs. They exist in a broad range of stellar forcing and rotational regimes that are distinctly different from solar system planets and more commonly studied hot Jupiters. Here we present simulations that explore global atmospheric circulation of sub-Neptunes generated with a two-dimensional shallow-water model, SWAMPE. We explore the circulation regimes of synchronously rotating sub-Neptunes with a focus on the interaction of planetary rotation rate and radiative timescale in a variety of stellar insolations. In highly irradiated, short-timescale regimes, our models exhibit high day-night geopotential contrasts. As the timescales become longer, the geopotential contrasts and longitudinal variability decrease, while temporal variability increases. The transition from day-to-night flow to jet-dominated flow is primarily driven by the radiative timescale. Strong- and medium-forcing regimes exhibit transitions between day-to-night flow and jet-dominated flow at similar points in the parameter space. Weak-forcing regime differs due to comparatively stronger rotational effects. Planetary rotation period dominates in determining equator-to-pole geopotential contrast. Our simulations exhibit higher time variability when either radiative timescale or rotation period is long.

Luminous infrared galaxies (LIRGs), the most extreme star-forming galaxies in the nearby (D$<$30 Mpc) Universe, show a notable X-ray emission deficiency (up to a factor of $\sim$10) compared with predictions from scaling relations of galaxy-wide high mass X-ray binary (HMXB) luminosity with star-formation rate. In the nearby ($\approx$20 Mpc) LIRG NGC 7552, the majority of the IR emission originates in a circumnuclear starburst ring, which has been resolved into several discrete knots of star formation. We present results from recent Chandra observations of NGC 7552, which reveal significant deficits in the 2-7 keV X-ray luminosities from two of the most powerful star-forming knots. We hypothesize that the expected luminous HMXB populations in these knots are either (1) obscured by very large column densities or (2) suppressed due to the knots having relatively high metallicity and/or very young ages ($\lesssim$ 5 Myr). We distinguish between these possibilities using data from recent NuSTAR observations, whose sensitivity above 10 keV is capable of uncovering heavily obscured HMXB populations, since emission at these energies is more immune to absorption effects. We find no evidence of a heavily obscured HMXB population in the central region of NGC 7552, suggesting suppressed HMXB formation. We further show that metallicity-dependent scaling relations cannot fully account for the observed deficit from the most powerful star-forming knots or the central region as a whole. Thus, we suggest that recent bursts in local star formation activity likely drive the high $L_{\rm{IR}}$ within these regions on timescales $\lesssim$ 5 Myr, shorter than the timescale required for the formation of HMXBs.

Space-based stellar coronagraph instruments aim to directly image exoplanets that are a fraction of an arcsecond separation and ten billion times fainter than their host star. To achieve this, one or more deformable mirrors (DMs) are used in concert with coronagraph masks to control the wavefront and minimize diffracted starlight in a region of the image known as the ``dark zone" or ``dark hole." The DMs must have a high number of actuators (50 to 96 across) to allow dark holes that are large enough to image a range of desired exoplanet separations. In addition, the surfaces of the DMs must be controlled at the picometer level to enable the required contrast. Any defect in the mechanical structure of the DMs or electronic system could significantly impact the scientific potential of the mission. Thus, NASA's Exoplanet Exploration Program (ExEP) procured two 50$\times$50 microelectromechanical (MEMS) DMs manufactured by Boston Micromachines Corporation (BMC) to test their robustness to the vibrational environment that the DMs will be exposed to during launch. The DMs were subjected to a battery of functional and high-contrast imaging tests before and after exposure to flight-like random vibrations. The DMs did not show any significant functional nor performance degradation at $10^{-8}$ contrast levels.

Direct imaging is the primary technique currently used to detect young and warm exoplanets and understand their formation scenarios. The extreme flux ratio between an exoplanet and its host star requires the use of coronagraphs to attenuate the starlight and create high contrast images. However, their performance is limited by wavefront aberrations that cause stellar photons to leak through the coronagraph and on to the science detector preventing the observation of fainter extrasolar companions. The VLT/SPHERE instrument takes advantage of its efficient adaptive optics system to minimize dynamical aberrations to improve the image contrast. In good seeing conditions, the performance is limited by quasi-static aberrations caused by slowly varying aberrations and manufacturing defects in the optical components. The mitigation of these aberrations requires additional wavefront sensing and control algorithms to enhance the contrast performance of SPHERE. Dark hole algorithms initially developed for space-based application and recently performed on SPHERE calibration unit have shown significant improvement in contrast. This work presents a status update of dark hole algorithms applied on SPHERE and the results obtained during the on-sky tests performed on February 15th 2022.

The probability density function of the arrival time of \v{C}erenkov light on a photo-multiplier tube has been studied. This study covers light production, transmission and detection. The light production includes the light from a muon, the light from a shower and the light due to the energy loss of a muon. For the transmission of light, the effects of dispersion, absorption and scattering in the medium are considered. For the detection of light, the angular acceptance and the quantum efficiency of the photo-multiplier tube are taken into account.

We use the TNG300 magneto-hydrodynamic simulation and mock catalogues built using subhalo abundance matching (SHAM) to study the origin of the redshift evolution of the halo occupation distribution (HOD). We analyse stellar-mass selected galaxy samples with fixed number densities, spanning the redshift range $0 \le z \le 3$. We measure their halo occupation functions and fit the HOD parameters to study their evolution over cosmic time. The TNG300 galaxy population strongly depends on the baryonic physics implemented in the simulation. In contrast, the galaxy population predicted by a basic SHAM model without scatter is a direct result of the cosmology of the dark matter simulation. We find that the HOD evolution is similar for both models and is consistent with a previous study of the HOD evolution in semi-analytical models. Specifically, this is the case for the ratio between the characteristic halo masses for hosting central and satellite galaxies. The only HOD parameter whose evolution varies across models is $\sigma_{\rm logM}$, which contains information about the stellar mass-halo mass relation of the galaxies and does not strongly impact galaxy clustering. We also demonstrate that the dependence on the specific values of the cosmological parameters is small. We conclude that the cosmology of the galaxy sample, i.e. the cosmological hierarchical growth of structure, and not the baryonic physics prescriptions, governs the evolution of the HOD for stellar mass-selected samples. These results have important implications for populating simulated lightcones with galaxies and can facilitate the interpretation of clustering data at different redshifts.

We establish a new and cosmological-model-independent method to explore the cosmic background dynamics in this work. Utilizing the latest Pantheon+ type Ia supernova sample and the Hubble parameter measurements, we obtain the values of the Hubble parameter and the deceleration parameter at five different redshift points ranging from 0.2 to 0.6, and find that they can deviate from the predictions of the $\Lambda$CDM model at more than $2\sigma$. We further probe the equation of state of dark energy and obtain that a slightly oscillating equation of state of dark energy around the $-1$ line is favored.

Cyanides, ranging from three carbon atoms to PAHs, and alkenyl compounds are abundant in the interstellar medium (ISM). Aminoacrylonitrile (3-Amino-2-propenenitrile, H$_{2}$N-CH=CH-CN), an alkenyl cyanide, thus represents a promising candidate for new interstellar detection. A comprehensive spectroscopic laboratory investigation of aminoacrylonitrile in its rotational ground vibrational state has been herein performed. The measurements carried out up to the THz regime made it possible to generate a precise set of reliable rest frequencies for its search in space up to sub-millimetre wavelengths. The $Z$-aminoacrylonitrile ($Z$-apn) isomer spectrum has been recorded employing a source-modulated sub-millimetre spectrometer, from 80 GHz to 1 THz. A combination of Doppler and sub-Doppler measurement regimes allowed to record 600 new lines. The collected data have enabled the characterisation of a set of spectroscopic parameters up to decic centrifugal distortion constants. The catalogue generated from the improved spectral data has been used for the search of $Z$-apn in the spectral survey of the G+0.693-0.027 molecular cloud located in the central molecular zone, in the proximity of the Galactic centre.

Calculations of stellar evolution at initial abundances of helium $Y=0.28$ and heavier elements $Z=0.014$ were done for stars with masses on the main sequence $1.7M_\odot\le M_\textrm{ZAMS}\le 5.2M_\odot$. Evolutionary sequences corresponding to the AGB stage were used for modelling the pulsation period decrease observed for almost two centuries in the Mira--type variable R Hya. Diminution of the period from $\Pi\approx$ 495 d in the second half of the eighteenth century to $\Pi\approx 380$ d in the 1950s is due stellar radius decrease accompanying dissipation of the radiation--diffusion wave generated by the helium flash. For all the history of its observations R Hya was the fundamental mode pulsator. The best agreement with observations is obtained for eight evolutionary models with initial mass $M_\textrm{ZAMS}=4.8M_\odot$ and the mass loss rate parameter of the Bl\"ocker formula $0.03\le\eta_\mathrm{B}\le 0.07$. Theoretical mass estimates of R Hya are in the range $4.44M_\odot\le M\le 4.63M_\odot$, whereas the mean stellar radius ($421R_\odot\le \bar R \le 445R_\odot$) corresponding to the pulsation period $\Pi\approx 380$ agrees well with measurements of the angular diameter by methods of the optical interferometric imaging.

Currently, the $\Lambda$ Cold Dark Matter model, which relies on the existence of cold dark matter and a cosmological constant $\Lambda$, best describes the Universe. However, we lack information in the high-redshift ($z$) region between Type Ia Supernovae (SNe Ia) (up to $z=2.26$) and the Cosmic Microwave Background ($z=1100$), an interval crucial to test cosmological models and their possible evolution. We have defined a sample of 983 Quasars up to $z=7.54$ with reduced intrinsic dispersion $\delta=0.007$ which determines the matter density parameter $\Omega_M$ with the same precision of SNe Ia. Although previous analysis have been used Quasars as cosmological tools (e.g. Risaliti and Lusso 2019), this is the first time that high-redshift sources, in this case Quasars, as standalone cosmological probes yield such tight constraints on $\Omega_M$. Our results show the importance of correcting cosmological relationships for selection biases and redshift evolution and how the choice of a golden sample reduces considerably the intrinsic scatter. This proves the reliability of Quasars as standard cosmological candles.

We report the discovery of HN Lib b, a sub-Neptunian mass planet orbiting the nearby ($d \approx$ = 6.25 pc) M4.0 V star HN Lib detected by our CARMENES radial-velocity (RV) survey. We determined a planetary minimum mass of $M_\text{b}\sin i = $ 5.46 $\pm$ 0.75 $\text{M}_\oplus$ and an orbital period of $P_\text{b} = $ 36.116 $\pm$ 0.029 d, using $\sim$5 yr of CARMENES data, as well as archival RVs from HARPS and HIRES spanning more than 13 years. The flux received by the planet equals half the instellation on Earth, which places it in the middle of the conservative habitable zone (HZ) of its host star. The RV data show evidence for another planet candidate with $M_\text{[c]}\sin i = $ 9.7 $\pm$ 1.9 $\text{M}_\oplus$ and $P_\text{[c]} = $ 113.46 $\pm$ 0.20 d. The long-term stability of the signal and the fact that the best model for our data is a two-planet model with an independent activity component stand as strong arguments for establishing a planetary origin. However, we cannot rule out stellar activity due to its proximity to the rotation period of HN Lib, which we measured using CARMENES activity indicators and photometric data from a ground-based multi-site campaign as well as archival data. The discovery adds HN Lib b to the shortlist of super-Earth planets in the habitable zone of M dwarfs, but HN Lib [c] probably cannot be inhabited because, if confirmed, it would most likely be an icy giant.

Can pulsar-like compact objects release further huge free energy besides the kinematic energy of rotation? This is actually relevant to the equation of states of cold supra-nuclear matter, which is still under hot debate. Enormous energy is surely needed to understand various observations, such as $\gamma-$ray bursts, fast radio bursts and soft $\gamma-$ray repeaters. The elastic/gravitational-free energy of solid strangeon star is revisited, with two approaches to calculate in general relativity. It is found that huge free energy (> $10^{46}$ erg) could be released via starquakes, given an extremely small anisotropy ($(p_{\rm t}-p_{\rm r})/p_{\rm r} \sim 10^{-4}$, with $p_{\rm t}$/$p_{\rm r}$ the tangential/radial pressures).

I present a study of the several contributions to the single-photon angular resolution of pair telescopes in the MeV energy range. I examine some test cases, the presently active {\sl Fermi} LAT, the ``pure-silicon'' projects ASTROGAM and AMEGO-X, and the emulsion-based project GRAINE.

We reconstruct the extragalactic dispersion measure \ -- redshift relation (${\rm DM_E}-z$ relation) from well-localized fast radio bursts (FRBs) using Bayesian inference method. Then the ${\rm DM_E}-z$ relation is used to infer the redshift and energy of the first CHIME/FRB catalog. We find that the distributions of extragalactic dispersion measure and inferred redshift of the non-repeating CHIME/FRBs follow cut-off power law, but with a significant excess at the low-redshift range. We apply a set of criteria to exclude events which are susceptible to selection effect, but find that the excess at low redshift still exists in the remaining FRBs (which we call Gold sample). The cumulative distributions of fluence and energy for both the full sample and the Gold sample do not follow the simple power law, but they can be well fitted by the bent power law. The underlying physical implications remain to be further investigated.

Recent advances in modeling density distributions, so-called neural density fields, can accurately describe the density distribution of celestial bodies without, e.g., requiring a shape model - properties of great advantage when designing trajectories close to these bodies. Previous work introduced this approach, but several open questions remained. This work investigates neural density fields and their relative errors in the context of robustness to external factors like noise or constraints during training, like the maximal available gravity signal strength due to a certain distance exemplified for 433 Eros and 67P/Churyumov-Gerasimenko. It is found that both models trained on a polyhedral and mascon ground truth perform similarly, indicating that the ground truth is not the accuracy bottleneck. The impact of solar radiation pressure on a typical probe affects training neglectable, with the relative error being of the same magnitude as without noise. However, limiting the precision of measurement data by applying Gaussian noise hurts the obtainable precision. Further, pretraining is shown as practical in order to speed up network training. Hence, this work demonstrates that training neural networks for the gravity inversion problem is appropriate as long as the gravity signal is distinguishable from noise. Code and results are available at https://github.com/gomezzz/geodesyNets

In addition to the Evershed flow directed from the umbra towards the outer boundary of the sunspot, under special circumstances, a counter Evershed flow (CEF) in the opposite direction also occurs. We aim to characterize the proper motions and evolution of three CEFs observed by the Solar Optical Telescope onboard the Japanese Hinode spacecraft and the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory. We use state-of-the-art inversions of the radiative transfer equation of polarized light applied to spectropolarimetric observations of the Fe I line pair around 630 nm. The three CEFs appeared within the penumbra. Two of the CEF structures, as part of their decay process, were found to move radially outwards through the penumbra parallel to the penumbral filaments with speeds, deduced from their proper motions, ranging between 65 and 117 m/s. In these two cases, a new spot appeared in the moat of the main sunspot after the CEFs reached the outer part of the penumbra. Meanwhile, the CEFs moved away from the umbra, and their magnetic field strengths decreased. The expulsion of these two CEFs seems to be related to the normal Evershed flow. The third CEF appeared to be dragged by the rotation of a satellite spot. Chromospheric brightenings were found to be associated with the CEFs, and those CEFs that reached the umbra-penumbra boundary showed enhanced chromospheric activity. The two CEFs, for which line-of-sight velocity maps were available during their formation phase, appear as intrusions into the penumbra. They may be associated with magnetic flux emergence.

In the coming decades, the space-based gravitational-wave (GW) detectors such as Taiji, TianQin, and LISA are expected to form a network capable of detecting millihertz GWs emitted by the mergers of massive black hole binaries (MBHBs). In this work, we investigate the potential of GW standard sirens from the Taiji-TianQin-LISA network in constraining cosmological parameters. For the optimistic scenario in which electromagnetic (EM) counterparts can be detected, we predict the number of detectable bright sirens based on three different MBHB population models, i.e., pop III, Q3d, and Q3nod. Our results show that the Taiji-TianQin-LISA network alone could achieve a constraint precision of $0.9\%$ for the Hubble constant, meeting the standard of precision cosmology. Moreover, the Taiji-TianQin-LISA network could effectively break the cosmological parameter degeneracies generated by the CMB data, particularly in the dynamical dark energy models. When combined with the CMB data, the joint CMB+Taiji-TianQin-LISA data offer $\sigma(w)=0.036$ in the $w$CDM model, which is close to the latest constraint result obtained from the CMB+SN data. We also consider a conservative scenario in which EM counterparts are not available. Due to the precise sky localizations of MBHBs by the Taiji-TianQin-LISA network, the constraint precision of the Hubble constant is expected to reach 1.1\%. In conclusion, the GW standard sirens from the Taiji-TianQin-LISA network will play a critical role in helping solve the Hubble tension and shedding light on the nature of dark energy.

A near-IR high-resolution, R=80000 spectrometer has been developed at IPAG to directly characterize the atmosphere of exoplanets using adaptive optics (AO) assisted telescopes, and a single-mode fiber-injection unit. A first technical test with the 200' Hale telescope at Palomar Observatory occurred in March 2022 using the PALM3000 AO system offered by this telescope. Observations have also been made at the same time with the PARVI spectrometer so that a direct comparison can be made between the two instruments. This spectrometer uses a virtually imaged phased array (VIPA) instead of an echelle grating, resulting in a very compact optical layout that fits in a 0.25m3 cryostat. Using a quarter of an H2RG detector, the spectrometer analyses the middle part of the H-band, from 1.57 to 1.7 microns for 2 sources whose light is transferred from the telescope to the spectrometer using single-mode fibers. By design, the transmission of the spectrometer is expected to be 40-50%, which is 2-3 times higher than the transmission of current high-resolution spectrometers such as CRIRES+ and NIRSPEC. A damaged cross-disperser limited it to 21%, however. A replacement grating with a correct, twice as high efficiency has been procured after the on-sky demonstration. In addition to recalling the main specifications of the VIPA spectrometer, this paper presents the control software, the calibration process, and the reduction pipeline that have been developed for the instrument. It also presents the results of the on-sky technical test with the Hale telescope, as well as measurements of the effective resolution and transmission, along with a comparison of a spectrum of the sun obtained with the spectrometer with the BASS2000 reference spectrum. Planned modifications are also discussed. That includes the integration of a new dedicated H2RG detector, and of K-band optics.

Quasi-periodic pulsations (QPPs) are frequently observed in solar and stellar flare emission, with recent studies suggesting that an increasing instantaneous period is a common characteristic of QPPs. Determining the prevalence of non-stationarity in QPPs contributes to a better understanding of which mechanisms are responsible in QPP generation. We obtain the rate of period evolution from QPPs in 98 M- and X-class flares from Solar Cycle 24 with average periods between 8-130s and investigate the prevalence of QPP non-stationarity. We also investigate whether the presence of a Coronal Mass Ejection (CME) impacts the period evolution of QPPs. We analyse soft X-ray lightcurves obtained from GOES' X-Ray Sensor (XRS) and assess the dominant periods in the impulsive and decay phases of the flares using the Fast Fourier Transform. We relate the rate of period evolution to flare duration, peak flare energy, and average QPP period. We find evidence of non-stationarity in 81% of the flares assessed, with most QPPs exhibiting a period evolution of less than 10s between the impulsive and decay phases, of which 66% exhibited an apparent period growth and 14% showed an apparent period shrinkage. We find a positive correlation between the absolute magnitude of period evolution and the duration of the flare and no correlation between the period evolution of the QPPs and flare energy or CME presence. Furthermore, we conclude that non-stationarity is common in solar QPPs and must be accounted for in flare analysis.

We present a multi-scale and multi-wavelength study to unveil massive star formation (MSF) processes around sites AFGL 5180, and AFGL 6366S, both hosting a Class II 6.7 GHz methanol maser emission. The radio continuum map at 8.46 GHz reveals a small cluster of radio sources toward AFGL 5180. Signatures of the early stages of MSF in our target sites are spatially seen at the opposite edges of a filamentary cloud (length $\sim$5 pc), which is observed in the sub-millimeter dust continuum maps. Using the near-infrared photometric data, the spatial distribution of young stellar objects is found toward the entire filament, primarily clustered at its edges. The getsf utility on the Herschel far-infrared images reveals a hub-filament system (HFS) toward each target site. The analysis of the molecular line data, which benefits from large area coverage ($\sim$1 degree $\times$ 1 degree), detects two cloud components with a connection in both position and velocity space. This supports the scenario of a cloud-cloud collision (CCC) that occurred $\sim$1 Myr ago. The filamentary cloud, connecting AFGL 5180 and AFGL 6366S, seems spatially close to an HII region Sh2-247 excited by a massive O9.5 star. Based on the knowledge of various pressures exerted by the massive star on its surroundings, the impact of its energetic feedback on the filamentary cloud is found to be insignificant. Overall, our observational outcomes favor the possibility of the CCC scenario driving MSF and the formation of HFSs toward the target sites.

We describe the defining observations of the solar cycle that provide constraints for the dynamo processes operating within the Sun. Specifically, we report on the following topics: historical sunspot numbers and revisions; active region (AR) flux ranges and lifetimes; tilt angles; Hale and Joy's law; the impact of rogue ARs on cycle progression; the spatio-temporal emergence of ARs that creates the butterfly diagram; polar fields; large-scale flows including zonal, meridional, and AR in-flows; short-term cycle variability; and helioseismic results including mode parameter changes.

In this work, we extend the model proposed by White concerning the post-collapse evolution of density peaks while considering the role of angular momentum. On a timescale smaller than the peak collapse, $t_{0}$, the inner regions of the peak reach the equilibrium forming a cuspy profile, as in White's paper, but the power-law density profile is flatter, namely $\rho \propto r^{-1.52}$, using the specific angular momentum $J$ obtained in theoretical models of how it evolves in CDM universes, namely $J \propto M^{2/3}$. The previous result shows how angular momentum influences the slope of the density profile, and how a slightly flatter profile obtained in high-resolution numerical simulations, namely $\rho \propto r^{\alpha}$, $(\alpha \simeq -1.5)$ can be reobtained. Similarly to simulations, in our model adiabatic contraction was not taken into account. This means that more comprehensive simulations could give different values for the slope of the density profile, similar to an improvement of our model.

We have performed 3-D numerical simulations to investigate the effect of partial ionization on the process of magnetic flux emergence. In our study, we have modified the single-fluid MHD equations to include the presence of neutrals and have performed two basic experiments: one that assumes a fully ionized plasma (FI case) and one that assumes a partially ionized plasma (PI case). We find that the PI case brings less dense plasma to and above the solar surface. Furthermore, we find that partial ionization alters the emerging magnetic field structure, leading to a different shape of the polarities in the emerged bipolar regions compared to the FI case. The amount of emerging flux into the solar atmosphere is larger in the PI case, which has the same initial plasma beta as the FI case, but a larger initial magnetic field strength. The expansion of the field above the photosphere occurs relatively earlier in the PI case, and we confirm that the inclusion of partial ionization reduces cooling due to adiabatic expansion. However, it does not appear to work as a heating mechanism for the atmospheric plasma. The performance of these experiments in three dimensions shows that PI does not prevent the formation of unstable magnetic structures, which erupt into the outer solar atmosphere.

The detection of a sub-solar mass black hole could yield dramatic new insights into the nature of dark matter and early-Universe physics, as such objects lack a traditional astrophysical formation mechanism. Gravitational waves allow for the direct measurement of compact object masses during binary mergers, and we expect the gravitational-wave signal from a low-mass coalescence to remain within the LIGO frequency band for thousands of seconds. However, it is unclear whether one can confidently measure the properties of a sub-solar mass compact object and distinguish between a sub-solar mass black hole or other exotic objects. To this end, we perform Bayesian parameter estimation on simulated gravitational-wave signals from sub-solar mass black hole mergers to explore the measurability of their source properties. We find that the LIGO/Virgo detectors during the O4 observing run would be able to confidently identify sub-solar component masses at the threshold of detectability; these events would also be well-localized on the sky and may reveal some information on their binary spin geometry. Further, next-generation detectors such as Cosmic Explorer and the Einstein Telescope will allow for precision measurement of the properties of sub-solar mass mergers and tighter constraints on their compact-object nature.

We present a new constraint on the Hubble constant $H_0$ using the latest measurements of the electromagnetic counterpart to the gravitational wave (GW) event GW170817. We use the latest optical, X-ray and radio observations of the afterglow up to $\sim 3.5$ years after the GW detection, and properly take into account the impact of the host galaxy peculiar velocity. We find $75.46^{+5.34}_{-5.39}$ km s$^{-1}$ Mpc$^{-1}$ (68\% Credible Interval), a $\sim7\%$ precision measurement, which is a significant improvement compared to the $14\%$ precision of the first standard siren measurement. Our result is consistent within $1\sigma$ with the Cepheid-anchored Supernova and within $1.5 \sigma$ with the Cosmic Microwave Background measurements of the Hubble constant. We also explore the impact of the various assumptions made when fitting for the afterglow on the Hubble constant estimate.

Wide-field radio continuum observations of galaxy clusters are revealing an increasing number of spiral galaxies hosting tens of kpc-long radio tails produced by the nonthermal interstellar medium being displaced by the ram pressure. We present a semi-empirical model for the multi-frequency radio continuum emission from ram pressure stripped tails based on the pure synchrotron cooling of a radio plasma moving along the stripping direction with a uniform velocity. We combine LOFAR and uGMRT observations at 144 and 400 MHz to study the flux density and spectral index profiles of the radio tails of 7 galaxies in Abell 2255, and use the model to reproduce the flux density and spectral index profiles, and infer the stripped radio plasma velocity. For 5 out of 7 galaxies we observe monotonic decrease in both flux density and spectral index up to $~30$ kpc from their stellar disk. Our model reproduces the observed trends with a radio plasma bulk projected velocity between 160 and 430 km s$^{-1}$. This result represents the first indirect measure of the stripped, nonthermal interstellar medium velocity. The observed spectral index trends indicate that the synchrotron cooling is faster than the adiabatic expansion losses, thus suggesting that the stripped radio plasma can survive for a few tens of Myr outside of the stellar disk. This provides a lower limit for the lifetime of the stripped ISM outside of the disk. As a proof of concept, we use the best-fit velocities to constrain the galaxies' 3D velocity in the cluster to be in the 300-1300 km s$^{-1}$. We estimate the ram pressure affecting these galaxies to be between 0.1 and 2.9 $\times10^{-11}$ erg cm$^{-3}$, and measure the inclination between their stellar disk and the ram pressure wind.

Modern stellar structure and evolution theory experiences a lack of observational calibrations for the interior physics of intermediate- and high-mass stars. This leads to discrepancies between theoretical predictions and observed phenomena mostly related to angular momentum and element transport. Analyses of large samples of massive stars connecting state-of-the-art spectroscopy to asteroseismology may provide clues on how to improve our understanding of their interior structure. We aim to deliver a sample of O- and B-type stars at metallicity regimes of the Milky Way and the Large Magellanic Cloud (LMC) galaxies with accurate atmospheric parameters from high-resolution spectroscopy, along with a detailed investigation of line-profile broadening, for future asteroseismic studies. After describing the general aims of our two Large Programs, we develop dedicated methodology to fit spectral lines and deduce accurate global stellar parameters from high-resolution multi-epoch UVES and FEROS spectroscopy. We use the best available atmosphere models for three regimes covered by our global sample, given its breadth in terms of mass, effective temperature, and evolutionary stage. Aside from accurate atmospheric parameters and locations in the Hertzsprung-Russell diagram, we deliver detailed analyses of macroturbulent line broadening, including estimation of the radial and tangential components. We find that these two components are difficult to disentangle from spectra with signal-to-noise ratios below 250. Future asteroseismic modelling of the deep interior physics of the most promising stars in our sample will improve the existing dearth of such knowledge for large samples of OB stars, including those of low metallicity in the LMC.

Primordial non-Gaussianity encodes vital information of the physics of the early universe, particularly during the inflationary epoch. To explore the local-type primordial non-Gaussianity, we study the anisotropies in gravitational wave background induced by the linear cosmological scalar perturbations during radiation domination in the early universe. We provide the first complete analysis to the angular power spectrum of such scalar-induced gravitational waves. The spectrum is expressed in terms of the initial inhomogeneities, the Sachs-Wolfe effect, and their crossing. It is anticipated to have frequency dependence and multipole dependence, i.e., $C_\ell(\nu)\propto [\ell(\ell+1)]^{-1}$ with $\nu$ being a frequency and $\ell$ referring to the $\ell$-th spherical harmonic multipole. In particular, the initial inhomogeneites in this background depend on gravitational-wave frequency. These properties are potentially useful for the component separation, foreground removal, and breaking degeneracies in model parameters, making the non-Gaussian parameter $f_{\mathrm{NL}}$ measurable. Further, theoretical expectations may be tested by space-borne gravitational-wave detectors in future.

The search for ephemeral liquid water on Mars is an ongoing activity. After the recession of the seasonal polar ice cap on Mars, small water ice patches may be left behind in shady places due to the low thermal conductivity of the Martian surface and atmosphere. During late spring and early summer, these patches may be exposed to direct sunlight and warm up rapidly enough for the liquid phase to emerge. To see the spatial and temporal occurrence of such ice patches, optical images should be searched for and checked. Previously a manual image analysis was conducted on 110 images from the southern hemisphere, captured by the High Resolution Imaging Science Experiment (HiRISE) camera onboard the Mars Reconnaissance Orbiter space mission. Out of these, 37 images were identified with smaller ice patches, which were distinguishable by their brightness, colour and strong connection to local topographic shading. In this study, a convolutional neural network (CNN) is applied to find further images with potential water ice patches in the latitude band between -40{\deg} and -60{\deg}, where the seasonal retreat of the polar ice cap happens. Previously analysed HiRISE images are used to train the model, each was split into hundreds of pieces, expanding the training dataset to 6240 images. A test run conducted on 38 new HiRISE images indicates that the program can generally recognise small bright patches, however further training might be needed for more precise predictions.Using a CNN model may make it realistic to analyse all available surface images, aiding us in selecting areas for further investigation.

If the construction of the ngVLA begins in 2026, its sensitivity is expected to match that of the VLA by late 2029. At that juncture it is anticipated that open-skies observing will cease on the VLA and commence on the ngVLA. We suggest that during 2026-2029 the VLA be held in a customized final configuration encompassing portions of its standard A, B, C and D configurations. Such a final VLA configuration would (1) help minimize the cost of VLA operations and maximize the pace of ngVLA construction and commissioning; (2) help VLA users pivot to the high-resolution, high-frequency research topics that are projected to headline the ngVLA science program; and (3) help mitigate the effects of source confusion during responses to transients in the era of the Rubin Observatory and LIGO A+.

Galaxies have been observed to exhibit a level of simplicity unexpected in the complex galaxy formation scenario posited by standard cosmology. This is particularly apparent in their dynamics, where scaling relations display much regularity and little intrinsic scatter. However, the parameters responsible for this simplicity have not been identified. Using the Spitzer Photometry & Accurate Rotation Curves galaxy catalogue, we argue that the radial acceleration relation (RAR) between galaxies' baryonic and total dynamical accelerations is the fundamental correlation governing the radial (in-disk) dynamics of late-type galaxies. In particular, we show that the RAR cannot be tightened by the inclusion of any other available galaxy property, that it is the strongest projection of galaxies' radial dynamical parameter space, and that all other statistical radial dynamical correlations stem from the RAR plus the non-dynamical correlations present in our sample. We further provide evidence that the RAR's fundamentality is unique in that the second most significant dynamical relation does not possess any of these features. Our analysis reveals the root cause of the correlations present in galaxies' radial dynamics: they are nothing but facets of the RAR. These results have important ramifications for galaxy formation theory because they imply that to explain statistically late-type galaxy dynamics within the disk it is necessary and sufficient to explain the RAR and lack of any significant, partially independent correlation. While simple in some modified dynamics models, this poses a challenge to standard cosmology.

A mathematical model and software implementation developed to predict trajectories of single lunar dust particles acted on by a high velocity gas flow is discussed. The model uses output from a computation fluid dynamics (CFD) or direct simulation Monte Carlo (DSMC) simulation of a rocket nozzle hot gas jet. The gas density, velocity vector field, and temperature predicted by the CFD/DSMC simulations, provide the data necessary to compute the forces and accelerations acting on a single particle of regolith. All calculations of trajectory assume that the duration of particle flight is much shorter than the change in gas properties, i.e., the particle trajectory calculations take into account the spatial variation of the gas jet, but not the temporal variation. This is a reasonable first-order assumption. Final results are compared to photogrammetry derived estimates of dust angles form Apollo landing videos.

The physical state of the lunar soil in the permanently shadowed craters of the moon is inferred from experimental investigation. The permanently shadowed craters do not undergo the same thermal cycling experienced by other parts of the moon and therefore could be slightly less compacted. This study is significant because excavating, roving, and landing interactions, along with the energy budgets and deployment schedules for associated technology, need to be scaled and designed properly. Results indicate that the degree of compaction due to thermal cycling is a function of the depth in the soil column.

In this paper we will present findings on the detection of Magnesium II (MgII, lambda = 2796 {\AA}, 2803 {\AA}) absorption systems observed in data from the Early Data Release (EDR) of the Dark Energy Spectroscopic Instrument (DESI). DESI is projected to obtain spectroscopy of approximately 3 million quasars (QSOs), of which over 99% are anticipated to be found at redshifts greater than z < 0.3, such that DESI would be able to observe an associated or intervening Mg II absorber illuminated by the background QSO. We have developed an autonomous supplementary spectral pipeline that detects such systems through an initial line-fitting process and then confirms line properties using a Markov Chain Monte Carlo (MCMC) sampler. Based upon both a visual inspection and the reanalysis of coadded observations, we estimate this sample of absorption systems to have a completeness of 82.56% and purity of 99.08%. As the spectra in which Mg II systems are detected are the result of coadding multiple observations, we can determine the sensitivity, and therefore completeness, of the sample by searching for known Mg II systems in coadded data with fewer observations (and therefore lower signal-to-noise). From a parent catalog containing 83,207 quasars, we detect a total of 23,921 Mg II absorption systems following a series of quality cuts. Extrapolating from this occurrence rate of 28.75% implies a catalog at the completion of the five-year DESI survey that contains over eight hundred thousand Mg II absorbers. The cataloging of these systems will enable significant further research as they carry information regarding circumgalactic medium (CGM) environments, the distribution of intervening galaxies, and the growth of metallicity across the redshift range 0.3 < z < 2.5.

The 21-cm signal is the most important measurement for us to understand physics during cosmic dawn. It is the key for us to understand the expansion history of the Universe and the nature of dark energy. In this paper, we focused on the characteristic 21-cm power spectrum of a special dynamic dark energy - the Interacting Chevallier-Polarski-Linder (ICPL) model - and compared it with those of the $\Lambda$CDM and CPL models. From the expected noise of HERA, we found more precise experiments in the future can detect the features of interacting dark energy in the 21-cm power spectra. By studying the brightness temperature, we found the ICPL model is closer to the observation of EDGES compared to the $\Lambda$CDM, thus alleviating the tension between theory and experiments.

Dust is important for star formation because it is the crucial component that couples gas to stellar radiation fields, allowing radiation feedback to influence gas fragmentation and thus the stellar initial mass function (IMF). Variations in dust abundance therefore provide a potential avenue by which variation in galaxy metallicity might affect the IMF. In this paper we present a series of radiation-magnetohydrodynamic simulations in which we vary the metallicity and thus the dust abundance from 1% of Solar to 3$\times$ Solar, spanning the range from the lowest metallicity dwarfs to the most metal-rich early-type galaxies found in the local Universe. We design the simulations to keep all dimensionless parameters constant so that the interaction between feedback and star-forming environments of varying surface density and metallicity is the only factor capable of breaking the symmetry between the simulations and modifying the IMF, allowing us to cleanly isolate and understand the effects of each environmental parameter. We find that at a fixed surface density more metal-rich clouds tend to form a slightly more bottom-heavy IMF than metal-poor ones, primarily because in metal-poor gas radiation feedback is able to propagate further, heating somewhat larger volumes of gas. However, shifts in IMF with metallicity at a fixed surface density are much smaller than shifts with surface density at fixed metallicity; metallicity-induced IMF variations are too small to explain the variations in mass-to-light ratio reported in galaxies of different mass and metallicity. We, therefore, conclude that metallicity variations are much less important than variations in surface density in driving changes in the IMF and that the latter rather than the former are most likely responsible for the IMF variations found in early-type galaxies.

We present the MARDELS catalog of 8,471 X-ray selected galaxy clusters over 25,000 deg^2 of extragalactic sky. The accumulation of deep, multiband optical imaging data, the development of the optical counterpart classification algorithm MCMF, and the release of the DESI Legacy Survey DR10 catalog covering the extragalactic sky makes it possible -- for the first time, more than 30 years after the launch of the ROSAT X-ray satellite -- to identify the majority of the galaxy clusters detected in the ROSAT All-Sky-Survey source catalog (2RXS). The resulting 90% pure MARDELS catalog is the largest ICM-selected cluster sample to date. MARDELS probes a large dynamic range in cluster mass spanning from galaxy groups to the most massive clusters in the Universe. The cluster redshift distribution peaks at z~0.1 and extends to redshifts z~1. Out to z~0.4, the MARDELS sample contains more clusters per redshift interval (dN/dz) than any other ICM-selected sample. In addition to the main sample, we present two subsamples with 6,930 and 5,522 clusters, exhibiting 95% and 99% purity, respectively. We forecast the utility of the sample for a cluster cosmological study, using realistic mock catalogs that incorporate most observational effects, including the X-ray exposure time and background variations, the existence likelihood selection adopted in 2RXS and the impact of the optical cleaning with MCMF. Using realistic priors on the observable--mass relation parameters from a DES-based weak lensing analysis, we estimate the constraining power of the MARDELSxDES sample to be of 0.026, 0.033 and 0.15 ($1\sigma$) on the parameters $\Omega_\mathrm{m}$, $\sigma_8$ and $w$, respectively.

Ultra-light scalar fields and their non-interacting class, the so-called fuzzy dark matter (FDM), are candidates for dark matter, introduced to solve the small-scale problems of the standard cold dark matter. In this paper, we address whether the small-scale effects, specifically the quantum pressure, could leave sizable imprints on the large-scale statistics of the matter. For this purpose, We utilize the Effective Field Theory of Large Scale Structures (EFT of LSS) wherein small-scale physics is integrated and represented on large scales by only a set of free parameters. These parameters can be determined by fitting to the cosmological simulations. We use the \textit{Gadget-2} code to study the evolution of $512^3$ particles in a box of side length $250\,h^{-1}\,\mathrm{Mpc}$. Fitting EFT predictions to the simulation data, we determine the value of the speed of sound. We use the suppressed FDM initial conditions for the FDM case, sufficient to produce accurate -- enough for our purpose -- results on large scales. We perform three FDM simulations with different masses and compare their sound speed with the standard cold dark matter (CDM) simulation. We found that the FDM sound speed is slightly higher than CDM's. The deviation of the sound speed for FDM from CDM is larger for lower FDM masses. We conclude that the impact of the FDM is not limited to the small scales alone, and we can search for them by studying the matter on large scales. Though it is beyond the observations' scope today, it is possible to discriminate it with upcoming observations.

Close hyperbolic encounters of black holes (BHs) generate certain Burst With Memory (BWM) events in the frequency windows of the operational, planned, and proposed gravitational wave (GW) observatories. We present detailed explorations of the detectable parameter space of such events that are relevant for the LIGO-Virgo-KAGRA and the International Pulsar Timing Array (IPTA) consortia. The underlying temporally evolving GW polarization states are adapted from Cho et al. [Phys. Rev. D 98, 024039 (2018)] and therefore incorporate general relativistic effects up to the third post-Newtonian order. Further, we provide a prescription to ensure the validity of our waveform family while describing close encounters. Preliminary investigations reveal that optimally placed BWM events should be visible to megaparsec distances for the existing ground-based observatories. In contrast, maturing IPTA datasets should be able to provide constraints on the occurrences of such hyperbolic encounters of supermassive BHs to gigaparsec distances.

The LISA telescopes must exhibit an optical path length stability of $\frac{\mathrm{pm}}{\sqrt{\mathrm{Hz}}}$ in the mHz observation band to meet mission requirements. The optical truss interferometer is a proposed method to aid in the ground testing of the telescopes, as well as a risk-mitigation plan for the flight units. This consists of three Fabry-Perot cavities mounted to the telescope which are used to monitor structural displacements. We have designed and developed a fiber-based cavity injection system that integrates fiber components, mode-matching optics, and a cavity input mirror into a compact input stage. The input stages, paired with return mirror stages, can be mounted to the telescope to form the optical truss cavities. We performed a thorough sensitivity analysis using various simulation methods to support the fabrication and assembly of three first-generation prototype cavities, each of which exhibited a satisfactory performance based on our models.

Understanding the origin and structure of mean magnetic fields in astrophysical conditions is a major challenge. Shear flows often coexist in such astrophysical conditions and the role of flow shear on dynamo mechanism is only beginning to be investigated. Here, we present a direct numerical simulation (DNS) study of the effect of flow shear on dynamo instability for a variety of base flows with controllable mirror symmetry (i.e, fluid helicity). Our observations suggest that for helical base flow, the effect of shear is to suppress the small scale dynamo (SSD) action, i.e, shear helps the large scale magnetic field to manifest itself by suppressing SSD action. For non-helical base flows, flow shear has the opposite effect of amplifying the small-scale dynamo action. The magnetic energy growth rate ($\gamma$) for non-helical base flows are found to follow an algebraic nature of the form, $\gamma = - aS + bS^\frac{2}{3}$ , where a, b > 0 are real constants and S is the shear flow strength and $\gamma$ is found to be independent of scale of flow shear. Studies with different shear profiles and shear scale lengths for non-helical base flows have been performed to test the universality of our finding.

The impact of light scalars coupled conformally and disformally to matter on the geodetic and frame-dragging (FD) precessions is calculated. For larger frequencies the disformal interaction becomes increasingly relevant. We use several satellite experiments and Pulsar time of arrival (ToA) measurements to derive bounds on the couplings, combining the Gravity Probe B, LARES, LAGEOS and GRACE results with pulsar timings. Forecasts for future constraints on the conformal and the disformal couplings based on the GINGER experiment, i.e. a future measurement of the Sagnac effect on Earth, the motion of $S$-stars around the galactic centre and future pulsar timing observations are presented.

I present a new class of nonrelativistic, modified-gravity MOND theories. The three gravitational degrees of freedom of these ``TRIMOND'' theories are the MOND potential and two auxiliary potentials, one of which emerges as the Newtonian potential. Their Lagrangians involve a function of three acceleration variables -- the gradients of the potentials. So, the transition from the Newtonian to the MOND regime is rather richer than in the aquadratic-Lagrangian theory (AQUAL) and the quasilinear MOND theory (QUMOND), which are special cases of TRIMOND, each defined by a Lagrangian function of a single variable. In particular, unlike AQUAL and QUMOND whose deep-MOND limit (DML) is fully dictated by the required scale invariance, here, the scale-invariant DML still requires specifying a function of two variables. For one-dimensional (e.g., spherical) mass distributions, in all TRIMOND theories the MOND acceleration is a (theory specific, but system independent) function of the Newtonian acceleration; their variety appears in nonsymmetric situations. Also, they all make the salient, primary MOND predictions. For example, they predict the same DML virial relation as AQUAL and QUMOND, and thus the same DML $M-\sigma$ relation, and the same DML two-body force. Yet they can differ materially on secondary predictions. Such TRIMOND theories may be the nonrelativistic limits of scalar-bimetric relativistic formulations of MOND, such as BIMOND with an added scalar.

This paper outlines the centralised design and production of the Ultra-High-Vacuum sidearm and Laser-Stabilisation systems for the AION Ultra-Cold Strontium Laboratories. Commissioning data on the residual gas and steady-state pressures in the sidearm chambers, on magnetic field quality, on laser stabilisation, and on the loading rate for the 3D Magneto-Optical Trap are presented. Streamlining the design and production of the sidearm and laser stabilisation systems enabled the AION Collaboration to build and equip in parallel five state-of-the-art Ultra-Cold Strontium Laboratories within 24 months by leveraging key expertise in the collaboration. This approach could serve as a model for the development and construction of other cold atom experiments, such as atomic clock experiments and neutral atom quantum computing systems, by establishing dedicated design and production units at national laboratories.