Papers for Wednesday, Mar 13 2024

The Lagrangian perturbation theory (LPT) provides a simple yet powerful way of computing the nonlinear matter power spectrum, and it has been applied to biased tracers such as halos and galaxies. The number conservation of matter particles allows a simple relation between the fluctuations at the initial and the late times, which is essential in deriving the exact expression for the nonlinear matter power spectrum. Biased tracers in contrast evolve through mergers and accretion, violating the assumption of number conservation. Furthermore, a sample of halos with a narrow mass bin at initial time evolves to become halos with very different masses at late time, making it difficult to directly connect predictions of the LPT to real observations, as observed samples have similar properties. Here we use $N$-body simulations to test these core predictions of the LPT for biased tracers and demonstrate that the LPT predictions for halos overestimates the power spectrum at $z=0$ by a factor of three for the mass bin sample $\Delta \log M_h~(h^{-1}M_{\odot})= 0.5$ at~$z\simeq 3$, while the predictions can be made consistent if we impose by hand the number conservation of halos in the simulations throughout the entire evolution. In reality, LPT modeling of biased tracers involves marginalization of unknown bias parameters, alleviating the problem. We discuss the implications for field-level models based on the LPT applications.

Binary neutron star mergers are expected to produce fast dynamical ejecta, with mildly relativistic velocities extending to $\beta=v/c>0.6$. In a preceding paper, we derived an analytic description of the time-dependent radio to X-ray synchrotron flux produced by collisionless shocks driven by such fast ejecta into the interstellar medium, for spherical ejecta with broken power-law mass (or energy) distributions, $M(>\gamma\beta)\propto(\gamma\beta)^{-s}$ with $s=s_\text{KN}$ at $\gamma\beta<\gamma_0\beta_0$ and $s=s_\text{ft}$ at $\gamma\beta>\gamma_0\beta_0$ (where $\gamma$ is the Lorentz factor). Here, we extend our analysis and provide analytic expressions for the self-absorption frequency, the cooling frequency, and the observed angular size of the emitting region (which appears as a ring in the sky). For parameter values characteristic of merger calculation results -- a "shallow" mass distribution, $1<s_\text{KN}<3$, for the bulk of the ejecta (at $\gamma\beta\approx0.2$), and a steep, $s_\text{ft}>5$, "fast tail" mass distribution -- the analytic results reproduce well (to tens of percent accuracy) the results of detailed numeric calculations, a significant improvement over earlier order-of-magnitude estimates (based on extrapolations of results valid for $\gamma\beta\ll1$).

PDS 70 remains as the best laboratory to investigate the influence of giant planet formation on the structure of the parental disk. One of the most intriguing discoveries is the detection of a resolved inner disk from ALMA observations, extending up to the orbit of PDS 70b. This inner disk is challenging to explain because most of the dust particles are expected to be trapped at the outer edge of the gap open by PDS 70b and PDS 70c. By performing dust evolution models and radiative transfer simulations that match the gas disk masses obtained from recent thermochemical models of PDS 70, we find that when the minimum grain size in the models is larger than 0.1$\mu$m, there is an efficient filtration of dust particles, and the inner disk is depleted during the first million-year of dust evolution. Therefore, to maintain an inner disk, the minimum grain size in the models needs to be smaller than 0.1$\mu$m. Only when grains are that small, they are diffused and dragged along with the gas throughout the planets' gap. The small grains transported in the inner disk grow and drift therein, but the constant reservoir of dust particles that are trapped in the outer edge of the gap and that are continuously fragmenting allows refilling the inner disk on million-year timescales. Our flux predictions at millimetre wavelength of these models agree with ALMA observations. These models predict a spectral index of 3.2 in the outer disk and 3.6 in the inner disk. Our simple analytical calculations show what the inner disk water emission recently observed with JWST may originate from these ice-coated small grains that flow through the gap, grow and drift towards the innermost disk regions, reaching the water snowline. These models may mirror the history and evolution of our Solar System where Jupiter and Saturn played a crucial role shaping the architecture and properties of planets in our Solar System.

Fragmentation during the early stages of high-mass star formation is crucial for understanding the formation of high-mass clusters. We investigated fragmentation within thirty-nine high-mass star-forming clumps as part of the Atacama Large Millimeter/submillimeter Array (ALMA) Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES). Considering projection effects, we have estimated core separations for 839 cores identified from the continuum emission and found mean values between 0.08 and 0.32 pc within each clump. We find compatibility of the observed core separations and masses with the thermal Jeans length and mass, respectively. We also present sub-clump structures revealed by the 7 m-array continuum emission. Comparison of the Jeans parameters using clump and sub-clump densities with the separation and masses of gravitationally bound cores suggests that they can be explained by clump fragmentation, implying the simultaneous formation of sub-clumps and cores within rather than a step-by-step hierarchical fragmentation. The number of cores in each clump positively correlates with the clump surface density and the number expected from the thermal Jeans fragmentation. We also find that the higher the fraction of protostellar cores, the larger the dynamic range of the core mass, implying that the cores are growing in mass as the clump evolves. The ASHES sample exhibits various fragmentation patterns: aligned, scattered, clustered, and sub-clustered. Using the Q-parameter, which can help to distinguish between centrally condensed and subclustered spatial core distributions, we finally find that in the early evolutionary stages of high-mass star formation, cores tend to follow a subclustered distribution.

Environmental effects are believed to play an important yet poorly understood role in triggering accretion events onto the supermassive black holes (SMBHs) of galaxies (Active Galactic Nuclei; AGN). Massive clusters, which represent the densest structures in the Universe, provide an excellent laboratory to isolate environmental effects and study their impact on black hole growth. In this work, we critically review observational evidence for the preferential activation of SMBHs in the outskirts of galaxy clusters. We develop a semi-empirical model under the assumption that the incidence of AGN in galaxies is independent of environment. We demonstrate that the model is broadly consistent with recent observations on the AGN halo occupation at $z=0.2$, although it may overpredict satellite AGN in massive halos at that low redshift. We then use this model to interpret the projected radial distribution of X-ray sources around high redshift ($z\approx1$) massive ($>5 \times 10^{14} \, M_\odot$) clusters, which show excess counts outside their virial radius. Such an excess naturally arises in our model as a result of sample variance. Up to 20% of the simulated projected radial distributions show excess counts similar to the observations, which are however, because of background/foreground AGN and hence, not physically associated with the cluster. Our analysis emphasises the importance of projection effects and shows that current observations of $z\approx1$ clusters remain inconclusive on the activation of SMBHs during infall.

We report an updated mass and magnification model of galaxy cluster Abell 370 using new NIRCam and NIRISS data from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). Using Lenstool and a combination of archival HST and MUSE data with new JWST data as constraints, we derive an improved gravitational lensing model and extract magnifications of background galaxies with uncertainties. Using our best fit model, we perform a search for new multiply imaged systems via predicted positions. We report no new multiply imaged systems with identifiable redshifts, likely due to already very deep HST and Spitzer data, but confirm a $z\sim8$ multiply imaged system by measuring its redshift with NIRISS and NIRSpec spectra. We find that the overall shape of the critical curve for a source at $z = 9.0$ is similar to previous models of Abell 370, with small changes. We investigate the $z\sim8$ galaxy with two images observable with an apparent magnitude in the F125W band of $26.0\pm0.2$ and $25.6\pm0.1$. After correcting for the magnifications of the images, 7.2$^{+0.2}_{-1.2}$ and 8.7$^{+0.4}_{-0.4}$, we use SED fitting to find an intrinsic stellar mass of log($M^*/M_{\odot})$ = 7.35$^{+0.04}_{-0.05}$, intrinsic SFR of 3.5$^{+2.2}_{-1.4}$ M$_{\odot}$/yr, and $M_{UV}$ of -21.3$^{+0.2}_{-0.2}$, which is close to the knee of the luminosity function at that redshift. Our model, and corresponding magnification, shear, and convergence maps are available on request and will be made publicly available on MAST in a CANUCS data release (DOI: 10.17909/ph4n-6n76).

We compute in detail the absorption optical depth for astrophysical $\gamma$-ray photons interacting with solar photons to produce electron positron pairs. This effect is greatest for $\gamma$-ray sources at small angular distances from the Sun, reaching optical depths as high as $\tau_{\gamma\gamma}\sim 10^{-2}$. We also calculate this effect including modifications to the absorption cross section threshold from subluminal Lorentz invariance violation (LIV). We show for the first time that subluminal LIV can lead to increases or decreases in $\tau_{\gamma\gamma}$ compared to the non-LIV case. We show that, at least in principle, LIV can be probed with this effect with observations of $\gamma$-ray sources near the Sun at $\gtrsim20$ TeV by HAWC or LHAASO, although a measurement will be extremely difficult due to the small size of the effect.

Studies show that both radio jets from the active galactic nuclei (AGN) and the star formation (SF) activity in quasar host galaxies contribute to the quasar radio emission; yet their relative contributions across the population remain unclear. Here, we present an improved parametric model that allows us to statistically separate the SF and AGN components in observed quasar radio flux density distributions, and investigate how their relative contributions evolve with AGN bolometric luminosity ($L_\mathrm{bol}$) and redshift ($z$) using a fully Bayesian method. Based on the newest data from LOFAR Two-Metre Sky Survey Data Release 2, our model gives robust fitting results out to $z\sim4$, showing a quasar host galaxy SFR evolution that increases with bolometric luminosity and with redshift out to $z\sim4$. This differs from the global cosmic SFR density, perhaps due to the importance of galaxy mergers. The prevalence of radio AGN emissions increases with quasar luminosity, but has little dependence on redshift. Furthermore, our new methodology and large sample size allow us to subdivide our dataset to investigate the role of other parameters. Specifically, in this paper, we explore quasar colour and demonstrate that the radio excess in red quasars is due to an enhancement in AGN-related emission, since the host galaxy SF contribution to the total radio emission is independent of quasar colour. We also find evidence that this radio enhancement occurs mostly in quasars with weak or intermediate radio power.

To provide constraints on the chemical processes responsible for the observed columns of organic species, we used NOEMA to observe the sight line toward NRAO150 in the 2mm spectral window. We targeted the low excitation lines of o-H2CO 2(1,1)-1(1,0) and p-H2CO 2(0,2)-1(0,1) as well as the nearby transitions of CS(3-2) and c-C3H2. We combined these data with previous observations to determine the excitation conditions, column densities, and abundances relative to H2 in the different velocity components. We performed non-LTE radiative transfer calculations including collision cross sections with ortho and para H2 and with electrons. New collision cross sections with electrons were computed for ortho and para formaldehyde. The c-C3H2 line profiles are very similar to those of HCO+ and CCH, while the CS absorption features are narrower and mostly concentrated in two main velocity components at V = -17 and -10 km/s. H2CO absorption lines present an intermediate pattern with absorption in all velocity components but larger opacities in the two main velocity components. The ortho-to-para ratios of H2CO and c-C3H2 are consistent with the statistical value of 3. While the excitation temperature of all c-C3H2 velocity components is consistent with the CMB, the two strong components detected in CS show a clear excess over the CMB indicating that CS resides at higher densities than other species along this particular sightline, n(H2) ~ 2500 cm-3 while n(H2) < 500 cm-3 for the other velocity components. We detected faint absorption from o-H213CO and C34S allowing us to derive isotopic ratios: o-H2CO/o-H213CO = 61 and C32S/C34S = 24. The excitation of the 4.8GHz line of formaldehyde is sensitive to the electron fraction and its excitation temperature is predicted to be lower than the CMB at low and moderate electron fractions, x(e)< 6E-5, and to rise above the CMB at high electron fractions, > 1e-4.

We have performed a comprehensive study of the Orion Nebula Cluster (ONC) combining the photometric data obtained by the two \textit{HST} Treasury programs that targeted this region. To consistently analyze the rich dataset obtained in a wide variety of filters, we adopted a Bayesian approach to fit the Spectral Energy Distribution of the sources, deriving mass, age, extinction, distance, and accretion for each source in the region. The three dimensional study of mass distribution for bona-fide cluster members shows that mass segregation in the ONC extends to sub-solar masses, while the age distribution strongly supports the idea that star formation in the ONC is best described by a major episode of star formation that happened $\sim 1$ Myr ago. For masses $\gtrsim 0.1$ \Msun, our derived empirical initial mass function (IMF) is in good agreement with a Chabrier system IMF. Both the accretion luminosity (\Lacc) and mass accretion rates (\dMacc) are best described by broken power-law relations. This suggests that for the majority of young circumstellar disks in this cluster the excess emission may be dominated by X-ray-driven photoevaporation by the central star rather than external photoevaporation. If this is the case, the slopes of the power-law relations may be largely determined by the initial conditions set at the onset of the star formation process, which may be quite similar between regions that eventually form clusters of different sizes.

X-ray obscuration of active galactic nuclei (AGNs) is considered in the context of ionized winds of stratified structure launched from accretion disks. We argue that a Compton-thick layer of a large-scale disk wind can obscure continuum X-rays and also lead to broad UV absorption such as in the blue wing of Civ; the former originates from the inner wind while the latter from the outer wind as a dual role. Motivated by a number of observational evidence showing strong AGN obscuration phenomena in Seyfert 1 AGNs, we demonstrate in this work, by utilizing a physically-motivated wind model coupled to post-process radiative transfer calculations, that an extended disk wind under certain physical conditions (e.g. morphology and density) could naturally cause a sufficient obscuration qualitatively consistent with UV/X-ray observations. Predicted UV/X-ray correlation is also presented as a consequence of variable spatial size of the wind in this scenario.

The exceptional spectra of the most luminous $z>10$ sources observed so far have challenged our understanding of early galaxy evolution, requiring a new observational benchmark for meaningful interpretation. As such, we construct spectroscopic templates representative of high-redshift, star-forming populations, using 482 confirmed sources at $z=5.0-12.9$ with JWST/NIRSpec prism observations, and report on their average properties. We find $z=5-11$ galaxies are dominated by blue UV continuum slopes ($\beta=-2.3$ to $-2.7$) and inverse Balmer jumps, characteristic of dust-poor and young systems, with a shift towards bluer slopes and younger ages with redshift. The evolution is mirrored by ubiquitous CIII] detections across all redshifts (EW$_{0}=5-14$ \r{A}), which increase in strength towards early times. Rest-frame optical lines reveal elevated ratios ($O32=7-31$, $R23=5-8$, and $Ne3O2=1-2$) and subsolar metallicities (log O/H$=7.3-7.9$), typical of ionization conditions and metallicities rarely observed in $z\sim0$ populations. Within our sample, we identify 57 Ly$\alpha$-emitters which we stack and compare to a matched sample of non-emitters. The former are characterized by more extreme ionizing conditions with enhanced CIII], CIV, and HeII+[OIII] line emission, younger stellar populations from inverse Balmer jumps, and a more pristine ISM seen through bluer UV slopes and elevated rest-frame optical line ratios. The novel comparison illustrates important intrinsic differences between the two populations, with implications for Ly$\alpha$ visibility. The spectral templates derived here represent a new observational benchmark with which to interpret high-redshift sources, lifting our constraints on their global properties to unprecedented heights and extending out to the earliest of cosmic times.

We tackle covariance estimation in low-sample scenarios, employing a structured covariance matrix with shrinkage methods. These involve convexly combining a low-bias/high-variance empirical estimate with a biased regularization estimator, striking a bias-variance trade-off. Literature provides optimal settings of the regularization amount through risk minimization between the true covariance and its shrunk counterpart. Such estimators were derived for zero-mean statistics with i.i.d. diagonal regularization matrices accounting for the average sample variance solely. We extend these results to regularization matrices accounting for the sample variances both for centered and non-centered samples. In the latter case, the empirical estimate of the true mean is incorporated into our shrinkage estimators. Introducing confidence weights into the statistics also enhance estimator robustness against outliers. We compare our estimators to other shrinkage methods both on numerical simulations and on real data to solve a detection problem in astronomy.

We present a study of energetic-electron trapping and acceleration in the Kelvin-Helmholtz-induced magnetohydrodynamic (MHD) turbulence of post-flare loops in the solar corona. Using the particle-tracing capabilities of MPI-AMRVAC 3.0, we evolve ensembles of test electrons (i.e. without feedback to the underlying MHD) inside the turbulent looptop, using the guiding-center approximation. With the MHD looptop model of Ruan et al. 2018, we investigate the relation between turbulence and particle trapping inside the looptop structure, showing that better-developed turbulent cascades result in more efficient trapping primarily due to mirror effects. We then quantify the electron acceleration in the time-evolving MHD turbulence, and find that ideal-MHD processes inside the looptop can produce nonthermal particle spectra from an initial Maxwellian distribution. Electrons in this turbulence are preferentially accelerated by mirror effects in the direction perpendicular to the local magnetic field while remaining confined within small regions of space between magnetic islands. Assuming dominance of Bremsstrahlung radiation mechanisms, we employ the resulting information from accelerated electrons (combined with the MHD background) to construct HXR spectra of the post-flare loop that include nonthermal-particle contributions. Our results pave the way to constructing more realistic simulations of radiative coronal structure for comparison with current and future observations.

We adopt a Bayesian X-ray spectral approach to investigate the accretion properties of unobscured ($20<\log(N_{\rm H}/{\rm cm}^{-2}<22$) and obscured ($22< \log(N_{\rm H}/{\rm cm}^{-2}<24$) active galactic nuclei (AGN) to shed light on the orientation vs evolution scenarios for the origin of the obscuring material. For a sample of 3882 X-ray-selected AGN from the Chandra COSMOS Legacy, AEGIS and CDFS extragalactic surveys, we constrain their stellar masses, $M_\star$, intrinsic X-ray luminosities, $L_{\rm X}$, obscuring column densities, $N_{\rm H}$, and specific accretion rates $\lambda\propto L_{\rm X}/M_\star$. By combining these observables within a Bayesian non-parametric approach, we infer, for the first time, the specific accretion rate distribution (SARD) of obscured and unobscured AGN to $z\approx3$, i.e. the probability of a galaxy with mass $M_\star$ at redshift $z$ hosting an AGN with column density $N_{\rm H}$ and specific accretion rate $\lambda$. Our findings indicate that (1) both obscured and unobscured SARDs share similar shapes, shifting towards higher accretion rates with redshift, (2) unobscured SARDs exhibit a systematic offset towards higher $\lambda$ compared to obscured SARD for all redshift intervals, (3) the obscured AGN fraction declines sharply at $\log\lambda_{\rm break} \sim-2$ for $z <0.5$, but shifts to higher $\lambda$ values with increasing redshift, (4) the incidence of AGN within the theoretically unstable blow-out region of the $\lambda-N_{\rm H}$ plane increases with redshift. These observations provide compelling evidence for AGN "downsizing" and radiation-regulated nuclear-scale obscuration with an increasing host galaxy contribution towards higher redshifts.

Be X-ray binaries (Be-XRBs) are crucial in understanding high-mass X-ray binaries, featuring a rapidly rotating Be star and a neutron star companion in an eccentric orbit, intermittently accreting material from the Be star's decretion disk. Originating from binary stellar evolution, Be-XRBs are of significant interest to binary population synthesis (BPS) studies, encapsulating the physics of supernovae, common envelope, and mass transfer (MT). Using the POSYDON BPS code, employing pre-computed grids of detailed binary stellar evolution models, we investigate the Galactic Be-XRB population. POSYDON incorporates stellar rotation self-consistently during MT phases, enabling a detailed examination of the rotational distribution of Be stars. Our fiducial BPS and Be-XRB model align well with the orbital properties of Galactic Be-XRBs, emphasizing the role of rotational constraints. Our modeling reveals a bimodal rotational distribution of Be-XRB-like systems, in excellent agreement with literature values. All Be-XRBs undergo an MT phase before the first compact object forms, with over half experiencing a second MT phase from a stripped helium companion (Case BB). Computing rotationally-limited MT efficiencies and applying them to our population, we find that the majority of Be-XRBs have undergone highly non-conservative MT (beta ~ 0.15). Our study underscores the importance of detailed angular momentum modeling during MT in interpreting Be-XRB populations, emphasizing this population as a key probe for the stability and efficiency of MT in interacting binaries.

Cold quasars are a rare population of luminous, unobscured quasars associated with host galaxies that have a high star formation rate. We aimed to study the host galaxies of sixty four of these cold quasars in order to probe how the supermassive black holes and host galaxies were coevolving. We compiled data from the XXL survey and crossmatched with the VHS, WISE, and HerMES surveys to obtain multiwavelength photometry spanning the Xray to the infrared and including optical spectroscopy. From the data, we calculated the supermassive black hole masses using broad emission from the magnesium II and hydrogen beta lines. We compared this with the stellar mass of the entire galaxy and find that the black holes are significantly more massive than would be predicted by local relations, indicating that the majority of black hole growth precedes the bulk of the the stellar mass formation. In addition to this, we created a spectral energy distribution for each galaxy to calculate the star formation rate. We compared the star formation rate with the black hole accretion rate and find that the stellar mass is rapidly increasing at a relative rate faster than the black hole growth, supporting the picture where the black hole grows first.

In star clusters, close stellar encounters can strongly impact the architecture of a planetary system or even destroy it. We present a systematic study on the effects of stellar flybys on two-planet systems. When such a system (with a modest initial planet semi-major axis ratio, $a_1/a_2=0.6-0.8$) experiences stellar flybys, one or both planets can be ejected, forming free-floating planets (FFPs), captured planets (CPs) around the flyby star, and free-floating binary planets (BPs); the remaining single-surviving-planets (SSPs) can have their orbital radii and eccentricities greatly changed. Through numerical experiments, we calculate the formation fractions (or branching ratios) of FFPs, SSPs, CPs and BPs as a function of the pericenter separation of the flyby (in units of the initial planet semi-major axis), and use them to derive the analytical expressions that can be used to compute the formation rates of FFPs, SSPs, CPs and BPs in general cluster environments. We find that the production rates of FFPs and SSPs are similar, while the rate for CPs is a factor of a few smaller. The formation fraction of free-floating BPs depends strongly on $a_1/a_2$ of the initial systems and on the planet masses. For Jupiter-mass planets, the formation fraction of BPs is always less than $1\%$ (for systems with $a_1/a_2=0.8$) and typically much smaller ($\lesssim 0.2\%$ for $a_1/a_2\lesssim 0.7$). The fraction remains less than $1\%$ when considering $4M_{\rm J}$ planets. Overall, when averaging over all possible flybys, the production rate of BPs is less than $0.1\%$ of that for FFPs. We also derive the velocity distribution of FFPs produced by stellar flybys, and the orbital semi-major axis and eccentricity distributions of SSPs, CPs and (rare) free-floating BPs. These results can find applications in future studies of exotic planets and planetary systems.

The ratio of the average tangential shear signal of different weak lensing source populations around the same lens galaxies, also known as a shear ratio, provides an important test of lensing systematics and a potential source of cosmological information. In this paper we measure shear ratios of three current weak lensing surveys -- KiDS, DES, and HSC -- using overlapping data from the Baryon Oscillation Spectroscopic Survey. We apply a Bayesian method to reduce bias in shear ratio measurement, and assess the degree to which shear ratio information improves the determination of important astrophysical parameters describing the source redshift distributions and intrinsic galaxy alignments, as well as cosmological parameters, in comparison with cosmic shear and full 3x2-pt correlations (cosmic shear, galaxy-galaxy lensing, and galaxy clustering). We consider both Fisher matrix forecasts, as well as full likelihood analyses of the data. We find that the addition of shear ratio information to cosmic shear allows the mean redshifts of the source samples and intrinsic alignment parameters to be determined significantly more accurately. Although the additional constraining power enabled by the shear ratio is less than that obtained by introducing an accurate prior in the mean source redshift using photometric redshift calibration, the shear ratio allows for a useful cross-check. The inclusion of shear ratio data consistently benefits the determination of cosmological parameters such as S_8, for which we obtain improvements up to 34%. However these improvements are less significant when shear ratio is combined with the full 3x2-pt correlations. We conclude that shear ratio tests will remain a useful source of cosmological information and cross-checks for lensing systematics, whose application will be further enhanced by upcoming datasets such as the Dark Energy Spectroscopic Instrument.

Renormalization of quantum loop effects generated from large fluctuations is a hugely debatable topic of research these days which rules out the Primordial Black Hole (PBH) formation within the framework of single-field inflation. In this article, we briefly discuss that the correct implementation of regularization, renormalization, and resummation techniques in a setup described by an ultra-slow-roll phase sandwiched between two slow-roll phases in the presence of smooth or sharp transitions can lead to a stringent constraint on the PBH mass (i.e. ${\cal O}(10^{2}{\rm gm}$)), which we advertise as a new No-go theorem. Finally, we will give some of the possible way-outs using which one can evade this proposed No-go theorem and produce solar/sub-solar mass PBHs.

We investigate the relationship between a dark matter halo's mass profile and measures of the velocity dispersion of kinematic tracers within its gravitational potential. By predicting the scaling relation of the halo mass with the aperture velocity dispersion, $M_\mathrm{vir} - \sigma_\mathrm{ap}$, we present the expected form and dependence of this halo mass tracer on physical parameters within our analytic halo model: parameterised by the halo's negative inner logarithmic density slope, $\alpha$, its concentration parameter, $c$, and its velocity anisotropy parameter, $\beta$. For these idealised halos, we obtain a general solution to the Jeans equation, which is projected over the line of sight and averaged within an aperture to form the corresponding aperture velocity dispersion profile. Through dimensional analysis, the $M_\mathrm{vir} - \sigma_\mathrm{ap}$ scaling relation is devised explicitly in terms of analytical bounds for these aperture velocity dispersion profiles: allowing constraints to be placed on this relation for motivated parameter choices. We predict the $M_{200} - \sigma_\mathrm{ap}$ and $M_{500} - \sigma_\mathrm{ap}$ scaling relations, each with an uncertainty of $60.5\%$ and $56.2\%$, respectively. These halo mass estimates are found to be weakly sensitive to the halo's concentration and mass scale, and most sensitive to the size of the aperture radius in which the aperture velocity dispersion is measured, the maximum value for the halo's inner slope, and the minimum and maximum values of the velocity anisotropy. Our results show that a halo's structural and kinematic profiles impose only a minor uncertainty in estimating its mass. Consequently, spectroscopic surveys aimed at constraining the halo mass using kinematic tracers can focus on characterising other, more complex sources of uncertainty and observational systematics.

Spirals and streamers are the hallmarks of mass accretion during the early stages of star formation. We present the first observations of a large-scale spiral and a streamer towards a very young brown dwarf candidate in its early formation stages. These observations show, for the first time, the influence of external environment that results in asymmetric mass accretion via feeding filaments onto a candidate proto-brown dwarf in the making. The impact of the streamer has produced emission in warm carbon-chain species close to the candidate proto-brown dwarf. Two contrasting scenarios, a pseudo-disk twisted by core rotation and the collision of dense cores, can both explain these structures. The former argues for the presence of a strong magnetic field in brown dwarf formation while the latter suggests that a minimal magnetic field allows large-scale spirals and clumps to form far from the candidate proto-brown dwarf.

Aditya-L1 is the first Indian space mission to explore the Sun and solar atmosphere with seven multi-wavelength payloads, with Visible Emission Line Coronagraph (VELC) being the prime payload. It is an internally occulted coronagraph with four channels to image the Sun at 5000 \AA~ in the field of view 1.05 - 3 \rsun, and to pursue spectroscopy at 5303 \AA, 7892 \AA~ and 10747 \AA~ channels in the FOV (1.05 - 1.5 \rsun). In addition, spectropolarimetry is planned at 10747 \AA~ channel. Therefore, VELC has three sCMOS detectors and one InGaAs detector. In this article, we aim to describe the technical details and specifications of the detectors achieved by way of thermo-vacuum calibration at the CREST campus of the Indian Institute of Astrophysics, Bangalore, India. Furthermore, we report the estimated conversion gain, full-well capacity, and readout noise at different temperatures. Based on the numbers, it is thus concluded that it is essential to operate the sCMOS detectors and InGaAs detectors at $-5^{\circ}$ and $-17^{\circ}$ C, respectively, at the spacecraft level.

We propose a novel method to constrain the tidal quality factor, $Q'$, from a non-synchronized star-planet system consisting of a slowly rotating low-mass star and a close-in Jovian planet, taking into account the tidal interaction and the magnetic braking. On the basis of dynamical system theory, the track of the co-evolution of angular momentum for such a system exhibits the existence of a forbidden region in the $\Omega_\mathrm{orb}$ -- $\Omega_\mathrm{spin}$ plane , where $\Omega_\mathrm{spin}$ and $\Omega_\mathrm{orb}$ denote the angular velocity of the stellar spin and planetary orbit, respectively. The forbidden region is determined primarily by the strength of the tidal interaction. By comparing ($\Omega_\mathrm{orb},\Omega_\mathrm{spin}$) of a single star-planet system to the forbidden region, we can constrain the tidal quality factor regardless of the evolutionary history of the system. The application of this method to the star-planet system, NGTS-10 -- NGTS-10 b, gives $Q' \gtrsim 10^8$, leading to an tight upper bound on the tidal torque. Since this cannot be explained by previous theoretical predictions for non-synchronized star-planet systems, our result requires mechanisms that suppress the tidal interaction in such systems.

Differentially rotating stars and planets transport angular momentum internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich--Schubert--Fricke (GSF) instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing angular momentum gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat, and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or ``parasitic saturation" models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline instability-driven turbulence in certain regimes.

[abridged] We aim to constrain the stellar initial mass function (IMF) during the epoch of reionization. To this purpose, we build up a semi-empirical model for the reionization history of the Universe, based on various ingredients: the latest determination of the UV galaxy luminosity function from JWST out to redshift $z\lesssim 12$; data-inferred and simulation-driven assumptions on the redshift-dependent escape fraction of ionizing photons from primordial galaxies; a simple yet flexible parameterization of the IMF $\phi(m_\star)\sim m_\star^\xi\, e^{-m_{\star,\rm c}/m_\star}$ in terms of a high-mass end slope $\xi<0$ and of a characteristic mass $m_{\star,\rm c}$ below which a flattening or a bending sets in; the PARSEC stellar evolution code to compute the UV and ionizing emission from different star's masses as a function of age and metallicity; a few physical constraints related to stellar and galaxy formation in faint galaxies at the reionization redshifts. We compare our model outcomes with the reionization observables from different astrophysical and cosmological probes, and perform Bayesian inference on the IMF parameters. We find that the IMF slope $\xi$ is within the range from $-2.8$ to $-2.3$, while appreciably flatter slopes are excluded at great significance. However, the bestfit value of the IMF characteristic mass $m_{\star,\rm c}\sim$ a few $M_\odot$ implies a suppression in the formation of small stellar masses, at variance with the IMF in the local Universe; this may be induced by the thermal background $\sim 20-30$ K provided by CMB photons at the reionization redshifts. Finally, we investigate the implications of our reconstructed IMF on the recent JWST detections of massive galaxies at and beyond the reionization epoch, showing that any putative tension with the standard cosmological framework is substantially alleviated.

Age-dating and weighting stellar populations in galaxies at various cosmic epochs are essential steps to study galaxy formation through cosmic times. Evolutionary population synthesis models with different input physics are used towards this aim. In particular, the contribution from the thermally pulsing asymptotic-giant-branch (TP-AGB) stellar phase, which peaks for intermediate-age 0.6-2 Gyr systems, has been debated upon for decades. Here we report the detection of strong cool star signatures in the rest-frame near-infrared spectra of three young (~1 Gyr), massive (~10^10 Msun) quiescent galaxies at large look-back time, z=1-2, using JWST/NIRSpec. The co-existence of oxygen- and carbon-type absorption features, spectral edges and features from rare species such as Vanadium, and possibly Zirconium, reveal a strong contribution from TP-AGB stars. Population synthesis models with significant TP-AGB contribution reproduce the observations considerably better than those with weak TP-AGB, which are those commonly used. These findings call for revisions of published stellar population fitting results, pointing to lower masses and younger ages, with additional implications on cosmic dust production and chemical enrichment. These results will stimulate new generations of improved models informed by these and future observations.

Upflows and downflows at active region (AR) boundaries have been frequently observed with spectroscopic observations at extreme ultraviolet (EUV) passbands. In this paper, we report the coexistence of upflows and downflows at the AR boundaries with imaging observations from the Solar Upper Transition Region Imager (SUTRI) and the Atmospheric Imaging Assembly (AIA). With their observations from 2022 September 21 to 2022 September 30, we find 17 persistent opposite flows occurring along the AR coronal loops. The upflows are prominent in the AIA 193 \AA images with a velocity of 50-200 km/s, while the downflows are best seen in the SUTRI 465 \AA and AIA 131 \AA images with a slower velocity of tens of kilometers per second (characteristic temperatures (log T(K)) for 193 \AA, 465 \AA and 131 \AA are 6.2, 5.7, 5.6, respectively). We also analyze the center-to-limb variation of the velocities for both upflows and downflows. The simultaneous observations of downflows and upflows can be explained by the chromosphere-corona mass-cycling process, in which the localized chromospheric plasma is impulsively heated to coronal temperature forming a upflow and then these upflows experience radiative cooling producing a downflow with the previously heated plasma returning to the lower atmosphere. In particular, the persistent downflows seen by SUTRI provide strong evidence of the cooling process in the mass cycle. For upflows associated with open loops, part of the plasma is able to escape outward and into the heliosphere as solar wind.

The NASA Lunar Reconnaissance Orbiter (LRO) has returned petabytes of lunar high spatial resolution surface imagery over the past decade, impractical for humans to fully review manually. Here we develop an automated method using a deep generative visual model that rapidly retrieves scientifically interesting examples of LRO surface imagery representing the first planetary image anomaly detector. We give quantitative experimental evidence that our method preferentially retrieves anomalous samples such as notable geological features and known human landing and spacecraft crash sites. Our method addresses a major capability gap in planetary science and presents a novel way to unlock insights hidden in ever-increasing remote sensing data archives, with numerous applications to other science domains. We publish our code and data along with this paper.

Although there is abundant and diverse observational evidence in support of white dwarf stars hosting planets or debris disks which form in the catastrophic destruction of various planetary bodies, the key processes that explain these observations are still being intensely investigated. The study of white dwarf planetary systems offers a unique perspective on exo-solar composition, that cannot be obtained by any other means. This chapter describes the various observational techniques that are used in order to find and characterize exo-planets and debris disks around white dwarfs. In turn, it discusses how to theoretically interpret these observations by surveying an array of various research tools and models currently employed in this field.

Neutron stars cooling passively since their birth may be reheated in their late-stage evolution by a number of possible phenomena: rotochemical, vortex creep, crust cracking, magnetic field decay, or more exotic processes such as removal of neutrons from their Fermi seas (the nucleon Auger effect), baryon number-violating nucleon decay, and accretion of particle dark matter. Using Exposure Time Calculator tools, we show that reheating mechanisms imparting effective temperatures of 2000--40000 Kelvin may be uncovered with excellent sensitivities at the James Webb Space Telescope (JWST), the Extremely Large Telescope (ELT), and the Thirty Meter Telescope (TMT), with imaging instruments operating from visible-edge to near-infrared. With a day of exposure, they could constrain the reheating luminosity of a neutron star up to a distance of 500 pc, within which about $10^5$ (undiscovered) neutron stars lie. Detection in multiple filters could overconstrain a neutron star's surface temperature, distance from Earth, mass, and radius. Using publicly available catalogues of newly discovered pulsars at the FAST and CHIME radio telescopes and the Galactic electron distribution models YMW16 and NE2001, we estimate the pulsars' dispersion measure distance from Earth, and find that potentially 30$-$40 of these may be inspected for late-stage reheating within viable exposure times, in addition to a few hundred candidates already present in the ATNF catalogue. Whereas the coldest neutron star observed (PSR J2144$-$3933) has an upper limit on its effective temperature of about 33000 Kelvin with the Hubble Space Telescope, we show that the effective temperature may be constrained down to 20000 Kelvin with JWST-NIRCam, 15000 Kelvin at ELT-MICADO, and 9000 Kelvin with TMT-IRIS. Campaigns to measure thermal luminosities of old neutron stars would be transformative for astrophysics and fundamental physics.

We introduce FLAME, a machine learning algorithm designed to fit Voigt profiles to HI Lyman-alpha (Ly$\alpha$) absorption lines using deep convolutional neural networks. FLAME integrates two algorithms: the first determines the number of components required to fit Ly$\alpha$ absorption lines, and the second calculates the Doppler parameter $b$, the HI column density N$_{\rm HI}$, and the velocity separation of individual components. For the current version of FLAME, we trained it on low-redshift Ly$\alpha$ forests observed with the Far Ultraviolet gratings of the Cosmic Origin Spectrograph (COS) aboard the Hubble Space Telescope (HST). Drawing on this data, we trained FLAME on $\sim$ $10^6$ simulated Voigt profiles, forward-modeled to Ly$\alpha$ absorption lines observed with HST-COS, to classify lines as either single or double components and then determine Voigt profile fitting parameters. FLAME shows impressive accuracy on the simulated data by identifying more than 98% (90%) of single (double) component lines. It determines $b$ values within $\approx \pm{8}~(15)$ km s$^{-1}$ and log $N_{\rm HI}/ {\rm cm}^2$ values within $\approx \pm 0.3~(0.8)$ for 90% of the single (double) component lines. However, when applied to real data, FLAME's component classification accuracy drops by $\sim$ 10%. Despite this, there is a reasonable agreement between the $b$ and N$_{\rm HI}$ distributions obtained from traditional Voigt profile fitting methods and FLAME's predictions. Our mock HST-COS data analysis, designed to emulate real data parameters, demonstrated that FLAME could achieve consistent accuracy comparable to its performance with simulated data. This finding suggests that the drop in FLAME's accuracy when used on real data primarily arises from the difficulty of replicating the full complexity of real data in the training sample.

We developed a grid of stellar rotation models for low-mass and solar-type Classical T Tauri stars (CTTS) ($0.3M_{\odot}<M_{\ast}<1.2M_{\odot}$). These models incorporate the star-disk interaction and magnetospheric ejections to investigate the evolution of the stellar rotation rate as a function of the mass of the star $M_{\ast}$, the magnetic field ($B_{\ast}$), and stellar wind ($\dot{M}_{wind}$). We compiled and determined stellar parameters for 208 CTTS, such as projected rotational velocity $v\sin(i)$, mass accretion rate $\dot{M}_{acc}$, stellar mass $M_{\ast}$, ages, and estimated rotational periods using TESS data. We also estimated a representative value of the mass-loss rate for our sample using the $[\text{O}\text{ I}]$ spectral line. Our results confirm that $v\sin(i)$ measurements in CTTS agree with the rotation rates provided by our spin models in the accretion-powered stellar winds (APSW) picture. In addition, we used the Approximate Bayesian Computation (ABC) technique to explore the connection between the model parameters and the observational properties of CTTS. We find that the evolution of $v\sin(i)$ with age might be regulated by variations in (1) the intensity of $B_{\ast}$ and (2) the fraction of the accretion flow ejected in magnetic winds, removing angular momentum from these systems. The youngest stars in our sample ($\sim $1 Myr) show a median branching ratio $\dot{M}_{wind}/\dot{M}_{acc}\sim$ $0.16$ and median $B_{\ast}\sim$ 2000 G, in contrast to $\sim 0.01$ and 1000 G, respectively, for stars with ages $\gtrsim 3$ Myr.

Magnetic sounding using data collected from the Juno mission can be used to provide constraints on Jupiter's interior. However, inwards continuation of reconstructions assuming zero electrical conductivity and a representation in spherical harmonics are limited by the enhancement of noise at small scales. In this paper we describe new reconstructions of Jupiter's internal magnetic field based on physics-informed neural networks and either the first 33 (PINN33) or the first 50 (PINN50) of Juno's orbits. The method can resolve local structures, and allows for weak ambient electrical currents. Compared with other methods, our reconstructions of Jupiter's magnetic field both on and above the surface are similar, and we achieve a similar fit to the Juno data. However, our models are not hampered by noise at depth, and so offer a much clearer picture of the interior structure. We estimate that the dynamo boundary is at a fractional radius of 0.8. At this depth, the magnetic field is arranged into longitudinal bands, and the great blue spot appears to be rooted in neighbouring structures of oppositely signed flux.

The Focusing Optics X-ray Solar Imager (FOXSI) sounding rocket payload flew for the second time on 2014 December 11. To enable direct Hard X-Ray (HXR) imaging spectroscopy, FOXSI makes use of grazing-incidence replicated focusing optics combined with fine-pitch solid-state detectors. FOXSI's first flight provided the first HXR focused images of the Sun. For FOXSI's second flight several updates were made to the instrument including updating the optics and detectors as well as adding a new Solar Aspect and Alignment System (SAAS). This paper provides an overview of these updates as well as a discussion of their measured performance.

We study the impact of gravitational wave memory on the distribution of far away light sources in the sky. For the first time we compute the built up of small, but permanent tensor distortions of the metric over cosmological time-scales using realistic models of compact binary coalescences (CBCs) whose rate of occurrence is extrapolated at $z\sim {\cal O}(1)$. This allows for a consistent computation of the random-walk like evolution of gravitational wave memory which, in turn, is used to estimate the overall shape and magnitude of astrometric deflections of far away sources of light. We find that for pulsar or quasar proper motions, the near-Earth contribution to the astrometric deflections dominates the result and the deflection is analogous to a stochastic gravitational wave memory background that is generally subdominant to the primary stochastic gravitational wave background. We find that this contribution can be within the reach of future surveys such as Theia. Finally, we also study the deviation of the presently observed angular distribution of quasars from perfect isotropy, which arises from the slow build-up of gravitational wave memory over the entire history of the universe. In this case, we find that astrometric deflections depend on the entire light trajectory from the source to the Earth, yielding a quadruple pattern whose magnitude is unlikely to be within reach of the next generation of astrometric surveys due to shot noise and cosmic variance limitations.

We identify a point-symmetric morphology of the supernova remnant (SNR) Cassiopeia A compatible with shaping by at least two, and more likely more than four, pairs of opposite jets, as expected in the jittering jets explosion mechanism (JJEM) of core-collapse supernovae. Using an old X-ray map of argon, we identify seven pairs of opposite morphological features that we connect with lines that cross each other at the same point on the plane of the sky. The opposite morphological features include protrusions, clumps, filaments, and funnels in the main SNR shell. In addition to these seven symmetry axes, we find two tentative symmetry axes (lines). These lines form a point-symmetric wind-rose. We place this point-symmetric wind-rose on a new JWST and X-ray images of Cassiopeia A. We find other morphological features and one more symmetry axis that strengthen the identified point-symmetric morphology. Not all symmetry axes correspond to jets; e.g., some clumps are formed by the compression of ejecta between two jet-inflated lobes (bubbles). The robust point-symmetric morphology in the iconic Cassiopeia A SNR strongly supports the JJEM and poses a severe challenge to the neutrino-driven explosion mechanism.

AstroSat observed transient neutron star low-mass X-ray binary XTE J1701-462 for a total duration of 135 ks during its 2022 outburst. The source traced a complete 'Z' shaped structure in the hardness intensity diagram (HID) during the observation. The source exhibited an extended horizontal branch and a short-dipping flaring branch in the HID. We find that most suitable spectral model comprises emission from a standard multi-color accretion disk and Comptonized radiation from a hot central corona. The observed disk component is cool having temperature in the range of 0.3-0.4 keV and truncated far ( ~100 - 520 km) from the compact object. The Compton corona has optical depth in the range of 3.3-5 and temperature in the range of 3.3-4.9 keV. The disk and corona flux as well as truncation radius varies significantly along the HID. We discuss the possible scenarios to explain the relationship between spectral evolution and motion of the source along the HID. The timing analysis revealed horizontal branch oscillations (HBOs) in the frequency range 34-60 Hz. The frequency and rms strength of HBO vary systematically as the source moves along the horizontal branch (HB). The observed correlation of HBO properties with the position on the HB is similar to that reported previously in this source using \textit{RXTE} data during the 2006 outburst of the source. The source also showed NBOs with frequency $\sim$ 6.7 Hz in the middle and lower normal branch. The energy dependent study of the HBO properties suggests that the HBO is stronger in the higher energy band. We also observed very-low frequency noise (VLFN) and band limited noise (BLN) components in the power-density spectra. The break frequency of BLN component was found to be tightly correlated with the HBO frequency. We discuss possible models to explain the origin and nature of the observed features in the PDS.

The goal of our work is to study the mode changes of the radio-quiet gamma-ray pulsar PSR J2021+4026 with improved detail. By accurately characterizing variations in the gamma-ray spectrum and pulse profile, we aim to relate the Fermi-LAT observations to theoretical models and interpret the mode changes in terms of variations in the structure of a multipolar dissipative magnetosphere. We continually monitored the rotational evolution and the gamma-ray flux of PSR J2021+4026 using more than 13 years of Fermi-LAT data with a binned likelihood approach. We clearly detect the previous mode changes and confirm a more recent mode change that occurred around June 2020. We investigated the features of the phase-resolved spectrum and pulse profile, and we inferred the macroscopic conductivity, the electric field parallel to the magnetic field, and the curvature radiation cutoff energy. These physical quantities are related to the spin-down rate and the gamma-ray flux and therefore are relevant to the theoretical interpretation of the mode changes. We computed the relative variations in the best-fit parameters, finding typical flux changes between 13% and 20%. Correlations appear between the gamma-ray flux and the spectral parameters, as the peak of the spectrum shifts by about 10% toward lower energies when the flux decreases. The analysis of the pulse profile reveals that the pulsed fraction of the light curve is larger when the flux is low. We introduced a simple magnetosphere model that combines a dipole field with a strong quadrupole component. We simulated magnetic field configurations to determine the positions of the polar caps for different sets of parameters, and we conclude that some configurations could explain the observed multiwavelength variability.

Several hundreds of thousands of eclipsing binaries (EBs) are expected to be detected in the Transiting Exoplanet Survey Satellite (TESS) full frame images (FFIs). This represents a significant increase in the number of EBs available for eclipse timing variation studies. In this paper, we investigate the feasibility of performing precise eclipse timing of TESS EBs using the FFIs. To this end, we developed a fast, automated method and applied it to a sample of $\sim$100 EBs selected from the Villanova TESS EB catalog. Our timing analysis resulted in the detection of ten new triple candidates with outer periods shorter than $\sim$1300$\,$d. For five of them, we were able to constrain the outer orbit by analyzing independently the short-cadence (SC) and FFI data and to derive the minimum mass of the third body with a precision better than 4 per cent for SC and 12 per cent for FFI data. We then compared the results obtained from the two datasets and found that using the FFI data leads to (1) a degradation of both the accuracy and precision of the tertiary mass determination for the tightest EBs and (2) an overall underestimation of the third component's mass. However, we stress that our main conclusions on the nature of the detected signals do not depend on which dataset is used. This confirms the great potential of TESS FFIs, which will allow us to search for rare objects such as substellar circumbinary companions and compact triple stellar systems.

In a recent preprint, Fernando et al. (2024) used public data from infrasound stations to constrain the localization of the fireball of the CNEOS 2014-01-08 (IM1) bolide. The analysis inferred a 90-percent-confidence ellipse with semi-minor and semi-major axes of 186 and 388 km, respectively. This large error ellipse includes the much better localization box derived by sensors aboard U.S. Government satellites which detected the fireball light. At the fireball's peak brightness, the CNEOS localization box documented by NASA/JPL measures 11.112km on a side and is centered on a latitude of 1.3S and a longitude of 147.6E. Here, we point out that the recent expedition to retrieve materials from IM1's site (Loeb et al. 2024a,b,c) surveyed a region of tens of km around the CNEOS box center, and was not dictated by the data studied by Fernando et al. (2024) because of its larger uncertainties.

Compact objects observed via gravitational waves are classified as black holes or neutron stars primarily based on their inferred mass with respect to stellar evolution expectations. However, astrophysical expectations for the lowest mass range, $\lesssim 1.2 \,M_\odot$, are uncertain. If such low-mass compact objects exist, ground-based gravitational wave detectors may observe them in binary mergers. Lacking astrophysical expectations for classifying such observations, we go beyond the mass and explore the role of tidal effects. We evaluate how combined mass and tidal inference can inform whether each binary component is a black hole or a neutron star based on consistency with the supranuclear-density equation of state. Low-mass neutron stars experience a large tidal deformation; its observational identification (or lack thereof) can therefore aid in determining the nature of the binary components. Using simulated data, we find that the presence of a sub-solar mass neutron star (black hole) can be established with odds $\sim 100:1$ when two neutron stars (black holes) merge and emit gravitational waves at signal-to-noise ratio $\sim 20$. For the same systems, the absence of a black hole (neutron star) can be established with odds $\sim 10:1$. For mixed neutron star-black hole binaries, we can establish that the system contains a neutron star with odds $\gtrsim 5:1$. Establishing the presence of a black hole in mixed neutron star-black hole binaries is more challenging, except for the case of a $\lesssim 1\,M_{\odot}$ black hole with a $\gtrsim 1\,M_{\odot}$ neutron star companion. On the other hand, classifying each individual binary component suffers from an inherent labeling ambiguity.

The accretion flow within the magnetospheric radius of bright X-ray pulsars can form an optically thick envelope, concealing the central neutron star from the distant observer. Most photons are emitted at the surface of a neutron star and leave the system after multiple reflections by the accretion material covering the magnetosphere. Reflections cause momentum to be transferred between photons and the accretion flow, which contributes to the radiative force and should thus influence the dynamics of accretion. We employ Monte Carlo simulations and estimate the acceleration along magnetic field lines due to the radiative force as well as the radiation pressure across magnetic field lines. We demonstrate that the radiative acceleration can exceed gravitational acceleration along the field lines, and similarly, radiation pressure can exceed magnetic field pressure. Multiple reflections of X-ray photons back into the envelope tend to amplify both radiative force along the field lines and radiative pressure. We analyze the average photon escape time from the magnetosphere of a star and show that its absolute value is weakly dependent on the magnetic field strength of a star and roughly linearly dependent on the mass accretion rate being $\sim 0.1\,{\rm s}$ at $\dot{M}\sim 10^{20}\,{\rm g\,s^{-1}}$. At high mass accretion rates, the escape time can be longer than free-fall time from the inner disc radius.

The easiest exoplanets to detect are those that orbit very close to their hoststars. As a result, even though these planets are quite rare, they represent amajor fraction of the current exoplanet population. A side-effect of theproximity between the planet and the star is that the two have strong mutualinteractions through a number of physical processes. One of the most importantof these processes is tides. Tides are thought to shape the orbits of close-inexoplanets, heat the planet making its radius expand, and even drive someplanets to spiral into their host stars. This chapter briefly introduces thebasics of tidal physics and describes the various fingerprints tides leavewithin the observed exoplanet population.

Galaxy mergers are expected to play a key role in the evolution of galaxies and their central supermassive black holes (SMBHs). An observational signature of this phenomenon is the detection of dual active galactic nuclei (AGNs) amongst merging systems, as predicted by cosmological models of structure formation. Dual AGNs at sub-kiloparsec-scale separation are the precursors of merging black hole binaries, an important source of gravitational waves, but a paucity of such systems have been confirmed to date by direct imaging, while other similar claims have been strongly disputed. Here we report the serendipitous multiwavelength discovery of a dual black hole system with a separation of ~100 pc, in the gas-rich luminous infrared galaxy MCG-03-34-64 (z=0.016). Chandra/ACIS imaging shows two spatially-resolved peaks of equal intensity in the neutral Fe Ka emission line, consistent with a dual SMBH, which is supported by Hubble Space Telescope (HST), and Very Large Array (VLA) observations. The separation of ~100 pc is the closest dual AGN separation reported to date with spatially-resolved, multiwavelength observations.

Tidal dissipation in a celestial body can be used to probe its internal structure. Tides govern the orbital evolution of binary systems and therefore constraints on the interior of binary system members can be derived by knowing the age and tidal state of the binary system. For asteroids, age estimates are challenging due to a lack of direct observation of their surface. However, the age of asteroid pairs formed by rotational fission of a parent body can be derived from dynamical modeling, and as such can be used to constrain the age of binary systems existing within asteroid pairs.We study 13 binary asteroid systems existing in asteroid pairs by modeling their tidal locking and eccentricity damping timescales from tidal dissipation in the primaries and secondaries. We consider the impact of thermal torques on these timescales from the YORP and BYORP effects. The resulting constraints on the tidal dissipation ratio Q/k2 are compared to monolithic and rubble pile asteroid theories, showing that all secondaries are consistent with rubble piles with regolith layers greater than 3m and suggest that Q/k2 for rubble piles increases with radius. A particular case is the first bound secondary of asteroid (3749) Balam, whose Q/k2 is constrained to be between 2.7x10^4 and 1.4x10^6, consistent with a rubble-pile with a regolith thickness between 30m and 100m.

The current interpretation of the observed late time cooling of transiently accreting neutron stars in low-mass X-ray binaries during quiescence requires the suppression of neutron superfluidity in their crust at variance with recent ab initio many-body calculations of dense matter. Focusing on the two emblematic sources KS~1731$-$260 and MXB~1659$-$29, we show that their thermal evolution can be naturally explained by considering the existence of a neutron superflow driven by the pinning of quantized vortices. Under such circumstances, we find that the neutron superfluid can be in a gapless state in which the specific heat is dramatically increased compared to that in the classical BCS state assumed so far, thus delaying the thermal relaxation of the crust. We have performed neutron-star cooling simulations taking into account gapless superfluidity and we have obtained excellent fits to the data thus reconciling astrophysical observations with microscopic theories. The imprint of gapless superfluidity on other observable phenomena is briefly discussed.

We investigate the role that dense environments have on the quenching of star formation and the transformation of morphology for a sample of galaxies selected from the Sloan Digital Sky Survey. We make a distinction between galaxies falling into groups $(13 \leq \log{(M_{\text{halo}}/M_{\odot})} < 14)$ and clusters $(\log{(M_{\text{halo}}/M_{\odot})} \geq 14)$, and compare to a large sample of field galaxies. Using galaxy position in projected phase space as a proxy for time since infall, we study how galaxy specific star formation rate (sSFR) and morphology, parameterized by the bulge-to-total light ratio (B/T), change over time. After controlling for stellar mass, we find clear trends of increasing quenched and elliptical fractions as functions of infall time for galaxies falling into both groups and clusters. The trends are strongest for low mass galaxies falling into clusters. By computing quenching and morphological transformation timescales, we find evidence that star formation quenching occurs faster than morphological transformation in both environments. Comparing field galaxies to recently infalling galaxies, we determine there is pre-processing of both star formation and morphology, with pre-processing affecting star formation rates more strongly. Our analysis favours quenching mechanisms that act quickly to suppress star formation, while other mechanisms that act on longer timescales transform morphology through bulge growth and disc fading.

iCOMs are species commonly found in the interstellar medium. They are believed to be crucial seed species for the build-up of chemical complexity in star forming regions as well as our own Solar System. Thus, understanding how their abundances evolve during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We use data from the ALMA Large Program FAUST to study the compact line emission towards the [BHB2007] 11 proto-binary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has previously been detected. More than 45 CH3OCHO lines are clearly detected, as well as 8 CH3OCH3 transitions , 1 H2CCO transition and 4 t-HCOOH transitions. We compute the abundance ratios with respect to CH3OH for CH3OCHO, CH3OCH3, H2CCO, t-HCOOH (as well as an upper limit for CH3CHO) through a radiative transfer analysis. We also report the upper limits on the column densities of nitrogen bearing iCOMs, N(C2H5CN) and N(C2H3CN). The emission from the detected iCOMs and their precursors is compact and encompasses both protostars, which are separated by only 0.2" (~ 28 au). The integrated intensities tend to align with the Southern filament, revealed by the high spatial resolution observations of the dust emission at 1.3 mm. A PV and 2D analysis are performed on the strongest and uncontaminated CH3OCH3 transition and show three different spatial and velocity regions, two of them being close to 11B (Southern filament) and the third one near 11A. All our observations suggest that the detected methanol, as well as the other iCOMs, are generated by the shocked gas from the incoming filaments streaming towards [BHB2007] 11A and 11B, respectively, making this source one of the few where chemical enrichment of the gas caused by the streaming material is observed.

The majority of stars are in binary/multiple systems. How such systems form in turbulent, magnetized cores of molecular clouds in the presence of non-ideal MHD effects remains relatively under-explored. Through ATHENA++-based non-ideal MHD AMR simulations with ambipolar diffusion, we show that the collapsing protostellar envelope is dominated by dense gravomagneto-sheetlets, a turbulence-warped version of the classic pseuodisk produced by anisotropic magnetic resistance to the gravitational collapse, in agreement with previous simulations of turbulent, magnetized single-star formation. The sheetlets feed mass, magnetic fields, and angular momentum to a Dense ROtation-Dominated (DROD) structure, which fragments into binary/multiple systems. This DROD fragmentation scenario is a more dynamic variant of the traditional disk fragmentation scenario for binary/multiple formation, with dense spiral filaments created by inhomogeneous feeding from the highly structured larger-scale sheetlets rather than the need for angular momentum transport, which is dominated by magnetic braking. Collisions between the dense spiraling filaments play a key role in pushing the local magnetic Toomre parameter $Q_\mathrm{m}$ below unity, leading to gravitational collapse and stellar companion formation provided that the local material is sufficiently demagnetized, with a plasma-$\beta$ of order unity or more. This mechanism can naturally produce {\it in situ} misaligned systems on the 100-au scale, often detected with high-resolution ALMA observations. Our simulations also highlight the importance of non-ideal MHD effects, which affect whether fragmentation occurs and, if so, the masses and orbital parameters of the stellar companions formed.

Exoplanet transmission spectra provide rich information about the chemical composition, clouds and temperature structure of exoplanet atmospheres. Most exoplanet transmission spectra only span infrared wavelengths ($\gtrsim$ 1 $\rm{\mu m}$), which can preclude crucial atmospheric information from shorter wavelengths. Here, we explore how retrieved atmospheric parameters from exoplanet transmission spectra change with the addition of optical data. From a sample of 14 giant planets with transit spectra from 0.3-4.5 $\rm{\mu m}$, primarily from the Hubble and Spitzer space telescopes, we apply a free chemistry retrieval to planetary spectra for wavelength ranges of 0.3-4.5 $\rm{\mu m}$, 0.6-4.5 $\rm{\mu m}$, and 1.1-4.5 $\rm{\mu m}$. We analyse the posterior distributions of these retrievals and perform an information content analysis, finding wavelengths below 0.6 $\rm{\mu m}$ are necessary to constrain cloud scattering slope parameters ($\log{a}$ and $\gamma$) and alkali species Na and K. There is limited improvement in the constraints on the remaining atmospheric parameters. Across the population, we find limb temperatures are retrieved colder than planetary equilibrium temperatures but have an overall good agreement with Global Circulation Models. As JWST extends to a minimum wavelength of 0.6 $\rm{\mu m}$, we demonstrate that exploration into complementing JWST observations with optical HST data is important to further our understanding of aerosol properties and alkali abundances in exoplanet atmospheres.

Solar and interplanetary radio bursts can reflect the existence and motion of energetic electrons and are therefore a kind of vital phenomenon in solar activities. The present study reported a solar radio burst (SRB) event observed by Parker Solar Probe (PSP) in its 8th orbital encounter phase, and it lasted about 20 hours in a frequency range of 0.5-15 MHz, called the type IV-like SRB. This type IV-like SRB consists of a series of numerous spikes with the center-frequency drifting slowly from ~5 MHz to ~1 MHz, and each individual spike appears a much faster frequency drifting and has a narrow frequency range of a few MHz and short duration of a few minutes. Based on the empirical models of the solar atmosphere adopted commonly, combining the in-situ measurement by PSP, we propose that these small-scale spikes were generated by a group of solitary kinetic Alfv\'en waves (SKAWs) in a magnetic loop accompanying coronal mass ejection (CME) and moving outwards, in which the frequency drifting of individual spike is caused by the SKAW's propagation and the center-frequency drifting may be attributed to the motion of the magnetic loop.

We present the first fully simulation-based hierarchical analysis of the light curves of a population of low-redshift type Ia supernovae (SNae Ia). Our hardware-accelerated forward model, released in the Python package slicsim, includes stochastic variations of each SN's spectral flux distribution (based on the pre-trained BayeSN model), extinction from dust in the host and in the Milky Way, redshift, and realistic instrumental noise. By utilising truncated marginal neural ratio estimation (TMNRE), a neural network-enabled simulation-based inference technique, we implicitly marginalise over 4000 latent variables (for a set of $\approx 100$ SNae Ia) to efficiently infer SN Ia absolute magnitudes and host-galaxy dust properties at the population level while also constraining the parameters of individual objects. Amortisation of the inference procedure allows us to obtain coverage guarantees for our results through Bayesian validation and frequentist calibration. Furthermore, we show a detailed comparison to full likelihood-based inference, implemented through Hamiltonian Monte Carlo, on simulated data and then apply TMNRE to the light curves of 86 SNae Ia from the Carnegie Supernova Project, deriving marginal posteriors in excellent agreement with previous work. Given its ability to accommodate arbitrarily complex extensions to the forward model -- e.g. different populations based on host properties, redshift evolution, complicated photometric redshift estimates, selection effects, and non-Ia contamination -- without significant modifications to the inference procedure, TMNRE has the potential to become the tool of choice for cosmological parameter inference from future, large SN Ia samples.

The Lightning Launch Commit Criteria (LLCC) are a set of complex rules to avoid natural and rocket-triggered lightning strikes to in-flight space launch vehicles. The LLCC are the leading source of scrubs and delays to space launches from Cape Canaveral Air Force Station (CCAFS) and NASA Kennedy Space Center (KSC). An LLCC climatology would be useful for designing launch concept of operations, mission planning, long-range forecasting, training, and setting LLCC improvement priorities. Unfortunately, an LLCC climatology has not been available for CCAFS/KSC. Attempts have been made to develop such a climatology, but they have not been entirely successful. The main shortfall has been the lack of a long continuous record of LLCC evaluations. Even though CCAFS/KSC is the world's busiest spaceport, the record of LLCC evaluations is not detailed enough to create the climatology. As a potential solution, the research in this study developed a proxy climatology of LLCC violations by using the long continuous record of surface electric field mills at CCAFS/KSC.

We explicitly construct the analogue of the \v{d}Alembert solution to the 1+1 wave equation in an hyperboloidal setting. This hyperboloidal \v{d}Alembert solution is used, in turn, to gain intuition into the behaviour of solutions to the wave equation in a hyperboloidal foliation and to explain some apparently anomalous behaviour observed in numerically constructed solutions discussed in the literature.

The correlators of large-scale fluctuations belong to the most important observables in modern cosmology. Recently, there have been considerable efforts in analytically understanding the cosmological correlators and the related wavefunction coefficients, which we collectively call cosmological amplitudes. In this work, we provide a set of simple rules to directly write down analytical answers for arbitrary tree-level amplitudes of conformal scalars with time-dependent interactions in power-law FRW universe. With the recently proposed family-tree decomposition method, we identify an over-complete set of multivariate hypergeometric functions, called family trees, to which all tree-level conformal scalar amplitudes can be easily reduced. Our method yields series expansions and monodromies of family trees in various kinematic limits, together with a large number of functional identities. The family trees are in a sense generalizations of polylogarithms and do reduce to polylogarithmic expressions for the cubic coupling in inflationary limit. We further show that all family trees can be decomposed into linear chains by taking shuffle products of all subfamilies, with which we find simple connection between bulk time integrals and boundary energy integrals.

We train a deep neural network (DNN) to output rates of dark matter (DM) induced electron excitations in silicon and germanium detectors. Our DNN provides a massive speedup of around $5$ orders of magnitude relative to existing methods (i.e. QEdark-EFT), allowing for extensive parameter scans in the event of an observed DM signal. The network is also lighter and simpler to use than alternative computational frameworks based on a direct calculation of the DM-induced excitation rate. The DNN can be downloaded $\href{https://github.com/urdshals/DEDD}{\text{here}}$.

We scrutinize the Sommerfeld enhancement in dark matter pair annihilation for $p$-wave and higher-$\ell$ partial waves. For the Yukawa potential these feature a super-resonant Breit-Wigner peak in their velocity-dependence close to Sommerfeld resonances as well as a universal scaling with velocity for all $\ell\geq 1$ that differs from the $s$-wave case. We provide a quantum mechanical explanation for these phenomena in terms of quasi-bound states sustained by the centrifugal barrier of the partial-wave potential, and give approximate WKB expressions capturing the main effects. The impact of quasi-bound states is exemplified for wino dark matter and models with light mediators, with a focus on indirect detection signals. We note that quasi-bound states can also explain similar peaks in the bound-state formation and self-scattering cross sections.

We study scalar cosmological perturbations in $f(R, T)$ modified gravity theories being $T$ the trace of the energy-momentum tensor. We provide detailed equations for the matter energy density contrast. We solve then numerically to promote a comparison with available large scale structure (LSS) formation observational data on $f \sigma_8$ and also addressing the $S_8$ tension. We identify $f(R,T)$ models that lead either to growth enhancement or suppression. Since recent results in the literature indicate a preference for the latter feature, this type of analysis is quite useful to select viable modifications of gravity. We studied class of such $f(R,T)$ models are either ruled out or severely restricted.

In this work, we explore the impact of higher dimensional spacetime on the stellar structure and thermodynamic properties of neutron stars. Utilizing the density-dependent relativistic hadron field theory, we introduce modifications to incorporate the influence of higher dimensionality, a novel approach not explored in existing literature to our best knowledge. Our methodology involves solving the essential stellar structure equations in D-dimensional spacetime ($D \geq 4$), starting with the modification of the Einstein-Hilbert action, derivation of the Einstein field equation in D dimensions, and application of the resulting exterior Schwarzschild spacetime metric for D-dimension. Our findings reveal that with incremental dimensions, the central density $\rho_{c} G_D$ and central pressure $p_c G_D$ gradually increase, leading to progressively stiffer neutron matter. Incremental dimensionality also results in a gradual increase in the maximum mass attained, limited to our study between $D=4$ and $D=6$, as no maximum mass value is obtained for $D>6$. We consistently observe the criteria $dM/d\rho_c>0$ fulfilled up to the maximum mass point, supported by stability analysis against infinitesimal radial pulsations. The validity of our solution is confirmed through causality conditions, ensuring that the matter sound speed remains within the speed of light for all cases. Additionally, our examination indicates that the total mass-to-radius ratio for all discussed D-dimensional cases comfortably resides within the modified Buchdahl limit, which exhibits the physical validity of achieved results.

We propose the $\alpha$-generalized no-scale supergravity, and study the corresponding inflationary models. With a new parameter $0<\alpha\leq 1$, the $\alpha$-generalized no-scale supergravity provides the continuous connections among the generic no-scale supergravity from string theory compactifications. The resulting prediction of the CMB, spectrum index $n_s$, and tensor-to-scalar ratio $r$ can be highly consistent with the latest Planck/BICEP/Keck Array observations. Notably, the models with $\alpha\neq 1$ give a smaller ratio $r\leq 10^{-3}$, which is flexible even under the anticipated tighter observational constraints at the future experiments. Additionally, these models have the potential to generate a broad-band stochastic gravitational wave background, and thus explain the NANOGrav 15yr signal. Furthermore, they predict the formation of primordial black holes (PBHs) with various mass scales, which could account for an significant portion of dark matter relic density in the Universe.

The radiation emitted from radon is a critical background in rare event search experiments conducted at the Yemi Underground Laboratory (Yemilab) in Jeongseon, Korea. A Radon Reduction System(RRS) has been developed and installed in Yemilab to reduce radon concentration in the air. The RRS primarily provides a purified air of 50 m3/h to the cleanroom used to assemble crystal detectors in the AMoRE, a neutrinoless double beta decay search experiment. RRS can reduce the radon level by a factor of 300, so a high-sensitivity radon detector was required. A highly sensitive radon detector was constructed using a 70 L chamber with a large PIN photodiode to measure radon concentration in the purified air. The radon detector shows an excellent resolution of 72 keV (FWHM) for 6.003 MeV alphas from 218Po decay and a sensitivity down to 23.8 +- 2.1 mBq/m3 with a boil-off N2 gas sample. The radon concentration level from the RRS measured by the radon detector was below 0.29 Bq/m3 with a reduction factor of about 300.

A semi--numerical approach proposed many years ago for describing gravitational collapse in the post--quasi--static approximation, is modified in order to avoid the numerical integration of the basic differential equations the approach is based upon. For doing that we have to impose some restrictions on the fluid distribution. More specifically, we shall assume the vanishing complexity factor condition, which allows for analytical integration of the pertinent differential equations and leads to physically interesting models. Instead, we show that neither the homologous nor the quasi--homologous evolution are acceptable since they lead to geodesic fluids, which are unsuitable for being described in the post--quasi--static approximation. Also, we prove that, within this approximation, adiabatic evolution also leads to geodesic fluids and therefore we shall consider exclusively dissipative systems. Besides the vanishing complexity factor condition, additional information is required for a full description of models. We shall propose different strategies for obtaining such an information, which are based on observables quantities (e.g. luminosity and redshift), and/or heuristic mathematical ansatz. To illustrate the method, we present two models. One model is inspired in the well known Schwarzschild interior solution, and another one is inspired in Tolman VI solution.

The interior of mature neutron stars is expected to contain superfluid neutrons and superconducting protons. The influence of temperature and currents on superfluid properties is studied within the self-consistent time-dependent nuclear energy-density functional theory. We find that this theory predicts the existence of a regime in which nucleons are superfluid (the order parameter remains finite) even though the energy spectrum of quasiparticle excitations exhibits no gap. We show that the disappearance of the gap leads to a specific heat that is not exponentially suppressed at low temperatures as in the BCS regime but can be comparable to that in the normal phase. Introducing some dimensionless effective superfluid velocity, we show that the behavior of the specific heat is essentially universal and we derive general approximate analytical formulas for applications to neutron-star cooling simulations.

The first loading of gadolinium (Gd) into Super-Kamiokande in 2020 was successful, and the neutron capture efficiency on Gd reached 50\%. To further increase the Gd neutron capture efficiency to 75\%, 26.1 tons of $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ was additionally loaded into Super-Kamiokande (SK) from May 31 to July 4, 2022. As the amount of loaded $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ was doubled compared to the first loading, the capacity of the powder dissolving system was doubled. We also developed new batches of gadolinium sulfate with even further reduced radioactive impurities. In addition, a more efficient screening method was devised and implemented to evaluate these new batches of $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$. Following the second loading, the Gd concentration in SK was measured to be $333.5\pm2.5$ ppm via an Atomic Absorption Spectrometer (AAS). From the mean neutron capture time constant of neutrons from an Am/Be calibration source, the Gd concentration was independently measured to be 332.7 $\pm$ 6.8(sys.) $\pm$ 1.1(stat.) ppm, consistent with the AAS result. Furthermore, during the loading the Gd concentration was monitored continually using the capture time constant of each spallation neutron produced by cosmic-ray muons,and the final neutron capture efficiency was shown to become 1.5 times higher than that of the first loaded phase, as expected.