Papers for Monday, May 05 2025

We investigate the relevance of {\sl self-organized criticality (SOC)} models in previously published empirical datasets, which includes statistical observations in astrophysics, geophysics, biophysics, sociophysics, and informatics. We study 25 interdisciplinary phenomena with five different event detection and power law fitting methods. The total number of analyzed size distributions amounts to 64 cases, of which 80\% are found to be nearly consistent ($\alpha_s=1.99\pm0.30$) with the SOC model predictions. The fractal-diffusive SOC model predicts power law slopes of $\alpha_F=(9/5)=1.80$ for the flux $F$, $\alpha_E=(5/3)\sim1.67$ for the fluence or energy $E$, and $\alpha_T=2.00$ for the avalanche duration $T$. We find that the phenomena of solar flares, earthquakes, and forest fires are consistent with the theoretical predictions, while the size distributions of other phenomena are not conclusive due to neglected background treatment, inadequacy of power law fitting range, small-number statistics, and finite-system size effects.

Significant numbers of free-floating planetary-mass objects have been discovered in nearby star-forming regions by the James Webb Space Telescope, including a substantial number (42) of Jupiter Mass Binary Objects ('JuMBOs') in the Orion Nebula Cluster. The JuMBOs have much wider separations than other populations of substellar binaries, and their existence challenges conventional theories of substellar and planetary-mass object formation. Whilst several theories have been proposed to explain their formation, there has yet to be a study that determines whether they could survive the dynamical encounters prevalent within a dense star-forming region. We place a population of planet-planet binaries in N-body simulations of dense star-forming regions and calculate their binary fraction over time. We find that between 50-90 per cent of planet-planet binaries are destroyed on timescales of a few Myr, which implies that many more must form if we are to observe them in their current numbers. Furthermore, if the ONC was much more dense at formation, the initial separation distribution of the JuMBOs must have been even wider (and less similar to other substellar binaries) than the observed distribution.

The kinematics of stars in dwarf spheroidal galaxies have been studied to understand the structure of dark matter halos. However, the kinematic information of these stars is often limited to celestial positions and line-of-sight velocities, making full phase space analysis challenging. Conventional methods rely on projected analytic phase space density models with several parameters and infer dark matter halo structures by solving the spherical Jeans equation. In this paper, we introduce an unsupervised machine learning method for solving the spherical Jeans equation in a model-independent way as a first step toward model-independent analysis of dwarf spheroidal galaxies. Using equivariant continuous normalizing flows, we demonstrate that spherically symmetric stellar phase space densities and velocity dispersions can be estimated without model assumptions. As a proof of concept, we apply our method to Gaia challenge datasets for spherical models and measure dark matter mass densities given velocity anisotropy profiles. Our method can identify halo structures accurately, even with a small number of tracer stars.

Adaptive optics systems are critical in any application where highly resolved imaging or beam control must be performed through a dynamic medium. Such applications include astronomy and free-space optical communications, where light propagates through the atmosphere, as well as medical microscopy and vision science, where light propagates through biological tissue. Recent works have demonstrated common-path wavefront sensors for adaptive optics using the photonic lantern, a slowly varying waveguide that can efficiently couple multi-moded light into single-mode fibers. We use the SCExAO astrophotonics platform at the 8-m Subaru Telescope to show that spectral dispersion of lantern outputs can improve correction fidelity, culminating with an on-sky demonstration of real-time wavefront control. To our best knowledge, this is the first such result for either a spectrally dispersed or a photonic lantern wavefront sensor. Combined with the benefits offered by lanterns in precision spectroscopy, our results suggest the future possibility of a unified wavefront sensing spectrograph using compact photonic devices.

We use a one-dimensional line-by-line radiative-convective model to simulate hot, dense terrestrial-planet atmospheres. We find that strong shortwave absorption by H2O and CO2 inhibits near-surface convection, reducing surface temperatures by up to approximately 2000 K compared to fully convective predictions. Pure CO2 atmospheres are typically 1000 K cooler than pure H2O atmospheres, with only a few percent of H2O needed to elevate surface temperatures by hundreds of kelvin for a fixed incident stellar radiation. We also show that minor greenhouse gases such as SO2 and NH3 have a limited warming effect when H2O is abundant. Even at insolation values as high as 12,500 W/m2 (about 37 times Earth's current solar flux), planets with mixed CO2-H2O envelopes have surface temperatures in the 1200 to 2000 K range, limiting surface melting. Our results highlight the critical role of shortwave heating on magma ocean planets and the need for improved high-temperature spectroscopy beyond 20,000 cm-1.

Using the RR Lyrae surveys Gaia DR3 Specific Objects Study, PanSTARRS1 and ASAS-SN-II, we determine the Milky Way's thick disc scale length and scale height as well as the radial scale length of the galaxy's inner halo. We use a Bayesian approach to estimate these values using two independent techniques: Markov chain Monte Carlo sampling, and importance nested sampling. We consider two vertical density profiles for the thick disc. In the exponential model, the scale length of the thick disc is $h_R=2.14_{-0.17}^{+0.19}$ kpc, and its scale height is $h_z=0.64_{-0.06}^{+0.06}$ kpc. In the squared hyperbolic secant profile $sech^2$, those values are correspondingly $h_R=2.10_{-0.17}^{+0.19}$ kpc and $h_z=1.02_{-0.08}^{+0.09}$ kpc. The density distribution of the inner halo can be described as a power law function with the exponent $n =-2.35_{-0.05}^{+0.05}$ and flattening $q =0.57_{-0.02}^{+0.02}$. We also estimate the halo to disc concentration ratio as $\gamma=0.19_{-0.02}^{+0.02}$ for the exponential disc and $\gamma=0.32_{-0.03}^{+0.03}$ for the $sech^2$ disc.

We present a new epoch of JWST spectroscopic variability monitoring of the benchmark binary brown dwarf WISE 1049AB, the closest, brightest brown dwarfs known. Our 8-hour MIRI low resolution spectroscopy (LRS) and 7-hour NIRSpec prism observations extended variability measurements for any brown dwarfs beyond 11 $\mu$m for the first time, reaching up to 14 $\mu$m. Combined with the previous epoch in 2023, they set the longest JWST weather monitoring baseline to date. We found that both WISE 1049AB show wavelength-dependent light curve behaviours. Using a robust k-means clustering algorithm, we identified several clusters of variability behaviours associated with three distinct pressure levels. By comparing to a general circulation model (GCM), we identified the possible mechanisms that drive the variability at these pressure levels: Patchy clouds rotating in and out of view likely shaped the dramatic light curves in the deepest layers between 1-2.5 $\mu$m, whereas hot spots arising from temperature / chemical variations of molecular species likely dominate the high-altitude levels between 2.5-3.6 $\mu$m and 4.3-8.5 $\mu$m. Small-grain silicates potentially contributed to the variability of WISE 1049A at 8.5-11 $\mu$m. While distinct atmospheric layers are governed by different mechanisms, we confirmed for the first time that each variability mechanism remains consistent within its layer over the long term. Future multi-period observations will further test the stability of variability mechanisms on this binary, and expanded JWST variability surveys across the L-T-Y sequence will allow us to trace and understand variability mechanisms across a wider population of brown dwarfs and planetary-mass objects.

Pulsar Timing Arrays (PTAs) have recently found strong evidence for low-frequency gravitational waves (GWs) in the nanohertz frequency regime. As GWs pass, they produce deviations in measured lengths and light-travel times. PTA experiments utilize the highly-consistent radio bursts from millisecond pulsars, distributed throughout the local galaxy, to identify miniscule timing deviations indicative of GWs. To distinguish GWs from noise, PTAs search for a particular correlation pattern between different pulsars called Hellings & Downs correlations. The type of GW signal that has recently been identified is a stochastic GW background (GWB), which is observed to have more power at lower GW frequencies. A GWB matching these observations has long been predicted from super-massive black-hole (SMBH) binaries. SMBHs are known to exist in the centers of galaxies, which can then form binaries when two SMBHs are brought together following the merger of galaxies. No example of an SMBH binary has confidently been identified to date, and tremendous uncertainties about their formation and evolution remain. Alternative sources of the GWB have also been proposed, based on models for new fundamental physics, particularly in the early Universe. Improved sensitivity of PTAs will eventual lead to the characterization of GWB anisotropy and constraints on GWs from individual SMBH binaries, either of which could definitively demonstrate the true origin of the GWB. If the source is SMBH binaries, a variety of electromagnetic counterparts are possible, allowing for multimessenger astrophysics with low-frequency GWs.

We report the first X-ray polarimetric results of the neutron star (NS) low-mass X-ray binary (LMXB) Z-source GX 349+2 using the Imaging X-ray Polarimetry Explorer (IXPE). We discovered that the X-ray source was polarized at PD = 1.1 +/- 0.3% (1-sigma errors) with a polarization angle of PA = 32 +/- 6 degree (1-sigma errors). Simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR) observations show that the source transitioned through the normal branch (NB), flaring branch (FB), and soft apex (SA) of the Z-track during our IXPE observations. The X-ray spectro-polarimetry results suggest a source geometry comprising an accretion disk component, a blackbody representing the emission from the NS surface, and a Comptonized component. We discuss the accretion geometry of the Z source in light of the spectro-polarimetric results.

The convectively driven, weakly magnetized regions of the solar photosphere dominate the Sun's surface at any given time, but the temporal variations of these quiet regions of the photosphere throughout the solar cycle are still not well known. To look for cycle-dependent changes in the convective properties of quiet Sun photosphere, we use high spatial and spectral resolution spectropolarimetric observations obtained by the Hinode Solar Optical Telescope (SOT) and apply the Spectropolarimetric Inversions Based on Response Functions (SIR) code to infer physical conditions in the lower solar photosphere. Using a homogeneous set of 49 datasets, all taken at disk center, we analyze the temperature stratification and the line-of-sight velocities of the granules and intergranules over a period of 15 years. We use a k-means clustering technique applied to the spectral profiles to segment the granules and intergranules based on both intensity and velocity. We also examine the profile bisectors of these different structures and compare these to past analyses. Our results show fairly constant properties over this period with no clear dependence on the solar cycle. We do, however, find a slight increase in the photospheric temperature gradient during the declining phase of the solar cycle. Our findings could have significant implications for understanding the coupling between the quiet Sun atmosphere and the global solar dynamo.

Fifty years after the very first sounding rocket measurement of cosmic X-ray polarization, the Imaging X-ray Polarimetry Explorer (IXPE) mission has effectively opened a new window into the X-ray sky. Prior to launch of IXPE, an extensive calibration campaign was carried out to fully characterize the response of this new type of instrument. Specifically, the polarization-sensitive detectors were intensively calibrated in Italy, where they were developed and built. The X-ray optics, which collect and focus X rays onto the detectors, were built and calibrated in the U.S. A key question was whether the telescope (optics + detectors) calibrations could be synthesized from the individual component calibrations, avoiding time consuming and costly end-to-end calibrations for a flight program with a fixed schedule. The data presented here are from a calibration of the flight spare telescope utilizing the flight spare detector and flight spare mirror assembly combined. These data show that the presence of the mirror module does not affect the polarization response of the detectors (within the required calibration accuracy) and that the angular resolution of the telescopes could be accurately determined. Thus, the original extensive stand-alone ground calibration data of all the flight detectors and all the flight optic can be utilized in full to derive the flight telescopes calibrations.

At visible wavelengths, Venus appears serene and pale-yellow, but since the 1920s, observers have noted high-contrast features in the ultraviolet. These features track the about 4-day superrotation of the upper cloud deck and vary widely over time and space. The identity of the UV absorber(s)-active between at least 280 and 500 nm-remains unknown, as no proposed candidate fully matches all observational data. From remote observations of Venus, and accounting for light scattering by sub-micrometer droplets, we modeled the 365-455 nm absorbance per cm of the bulk liquids forming Venus's clouds. Assuming a uniform distribution in mode 1 and 2 particles across a 6 km layer below the cloud top at 65 km, we constrain the bulk absorbance with a peak at A375 nm being 2942 per cm. This extremely high absorbance implies the presence of a highly efficient absorber, most likely conjugated organics, at relatively high concentration-e.g. about 25 g/L for porphyrin type pigments. Inorganic absorbers, with molar absorption coefficients typically in the range of 1,000-10,000 per M per cm, would either need to comprise a large portion of the aerosols or are simply not light absorbent enough, even if present in pure form. We emphasize that all candidate absorbers must be evaluated against Venus's reflectance curve using (i) known molar absorption coefficients, (ii) realistic atmospheric distributions, and (iii) appropriate particle size distributions. The upcoming Rocket Lab mission will test the hypothesis of organics in Venus's clouds.

We fit archival near-infrared spectra of 305 brown dwarfs with atmosphere models from the Sonora and Phoenix groups. Using the parameters of the best-fit models as estimates for the physical properties of the brown dwarfs in our sample, we have performed a survey of how brown dwarf atmospheres evolve with spectral type and temperature. We present the fit results and observed trends. We find that clouds have a more significant impact on near infrared spectra than disequilibrium chemistry, and that silicate clouds influence the near infrared spectrum through the late T types. We note where current atmosphere models are able to replicate the data and where the models and data conflict. We also categorize objects with similar spectral morphologies into families and discuss possible causes for their unique spectral traits. We identify two spectral families with morphologies that are likely indicative of binarity.

Forecasting future solar activity has become crucial in our modern world, where intense eruptive phenomena mostly occurring during solar maximum are likely to be strongly damaging to satellites and telecommunications. We present a 4D variational assimilation technique applied for the first time to real solar data. Our method is tested against observations of past cycles 22, 23, 24 and on the ongoing cycle 25 for which we give an estimate of the imminent maximum value and timing and also provide a first forecast of the next solar minimum. We use a variational data assimilation technique applied to a solar mean-field Babcock-Leighton flux-transport dynamo model. Ensemble predictions are produced in order to obtain uncertainties on the timing and value of the maximum of cycle $n+1$, when data on cycle $n$ is assimilated. We study in particular the influence of the phase during which data is assimilated in the model and of the weighting of various terms in the objective function. The method is validated on cycles 22, 23 and 24 with very satisfactory results. For cycle 25, predictions vary depending on the extent of the assimilation window but start converging past 2022 to a solar maximum reached between mid-2024 up to the beginning of 2025 with a sunspot number value of $143.1 \pm 15.0$. Relatively close values of the maximum are found in both hemispheres within a time lag of a few months. We also forecast a next minimum around late 2029, with still significant errorbars. The data assimilation technique presented here combining a physics-based model and real solar observations produces promising results for future solar activity forecasting.

Faraday rotation measurements of extragalactic radio sources occulted by the solar corona serve as a powerful complementary tool for probing the pre-eruption electron density and magnetic field structure. These measurements thereby allow us to refine predictions from global MHD models. In this paper, we discuss our recent study of the morphological evolution of a CME-driven shock event that occurred on August 3, 2012. Our analysis used white-light coronagraphic observations from three different vantage points in space (SOHO and STEREO A and B). Obtaining data from these spacecraft, we derived key parameters such as the radius of curvature of the driving flux rope, the shock speed, and the standoff distance from the CMEs' leading edge. A notable feature of this event was the availability of rare Faraday rotation measurements of a group of extragalactic radio sources occulted by the solar corona, which were obtained a few hours before the eruption. These observations from the VLA radio interferometer provide independent information on the integrated product of the line-of-sight magnetic field component and electron density. By modeling the shock standoff distance and using constraints from the Faraday rotation measurements, we achieve a high level of agreement between the fast-mode Mach number predicted by the Magnetohydrodynamic Algorithm outside a Sphere (MAS) code in its thermodynamic mode and the value deduced from the analysis of the 3D reconstruction of coronagraphic data, provided that appropriate correction factors (f_b = 2.4 and f_n = 0.5) are applied in advance to scale the simulated magnetic field and electron density, respectively. Our results are consistent with previous estimates and provide critical information for fine-tuning future MHD simulations.

Meteor light curves are sometimes known to display flickering: rapid, quasi-periodic variations in brightness. This effect is generally attributed to the rotational modulation of the ablation rate, which is caused by the time-varying cross section area presented by a nonspherical rotating meteoroid to the oncoming airflow. In this work we investigate the effects that the rotation of a meteoroid of given shape (spherical, cubic, or cylindrical) has on the meteor's light curve, given state-of-the-art experimental laboratory estimates of the drag and lift coefficients of hypersonic flow (Mach number > 5) around various shaped objects. The meteoroid's shape is important in determining these two forces, due to the different response of the drag and lift coefficients according to the angle of attack. As a case study, the model was applied to a fireball observed on 2018 April 17 by the PRISMA network, a system of all-sky cameras that achieves a systematic monitoring of meteors and fireballs in the skies over the Italian territory. The results show that this methodology is potentially able to yield a powerful diagnostic of the rotation rate of meteoroids prior to their encounter with the atmosphere, while also providing essential information on their pre-fall actual shapes.

Studying the atmospheres of exoplanets is one of the most promising ways to learn about distant worlds beyond our solar system. The composition of an exoplanet's atmosphere can provide critical insights into its geology and potential habitability. For instance, the presence of certain molecules such as water vapor, oxygen, or methane have been proposed to indicate the possibility of life. From an observation point of view, over the past fifteen years, significant progress has been made in characterizing exoplanetary atmospheres. This work reviews recent developments in ground-based high-resolution spectroscopic instruments that make it possible to analyze distant atmospheres in great detail. High-resolution transmission spectroscopy, one of the most effective methods used, has examined the atmospheres of Jupiter-like and is pushing towards the smaller, sub-Neptunian exoplanets. Numerous molecules have been detected using this technique, including CO,H2O,TiO,HCN,CH4,NH3,C2H2,OH. We explore the intriguing possibilities that lie ahead for future ground-based instrumentation, particularly in the context of detecting biologically relevant molecules within Earth-analog exoplanetary atmospheres including molecular oxygen (O2). With detailed exposure time calculations for detecting O2 we find that at the same exposure time spectral resolution of 300,000 reaches higher significance compared to 100,000. The exposure time and therefore the needed number of transits is reduced by a factor of 4 in challenging haze and cloud scenarios.

In the recent sixth data release (DR6) of the Atacama Cosmology Telescope (ACT) collaboration, the value of $n_{\rm s}=0.9743 \pm 0.0034$ for the scalar spectral index is reported, which excludes the Starobinsky and Higgs inflationary models at $2\sigma$ level. In this paper, we perform a Bayesian inference of the parameters of the Starobinsky or Higgs inflationary model with non-instantaneous reheating using the Markov chain Monte Carlo method. For the analysis, we use observational data on the cosmic microwave background collected by the Planck and ACT collaborations and on baryonic acoustic oscillations from the DESI collaboration. The reheating stage is modelled by a single parameter $R_{\rm reh}$. Using the modified Boltzmann code CLASS and the cobaya software with the GetDist package, we perform a direct inference of the model parameter space and obtain their posterior distributions. Using the Kullback--Leibler divergence, we estimate the information gain from the data, yielding $2.52$ bits for the reheating parameter. Inclusion of the ACT DR6 data provides $75\%$ more information about the reheating stage compared to analysis without ACT data. We draw constraints on the reheating temperature and the average equation of state. While the former can vary within $10$ orders of magnitude, values in the $95\,\%$ credible interval indicate a sufficiently low reheating temperature; for the latter there is a clear preference for values greater than $0.5$, which means that the conventional equations of state for dust $\omega=0$ and relativistic matter $\omega=1/3$ are excluded with more than $2\sigma$ level of significance. However, there still is a big part of parameter space where Starobinsky and Higgs inflationary models exhibit a high degree of consistency with the latest observational data, particularly from ACT DR6. Therefore, it is premature to reject these models.

Adding an independent estimate of the mean stellar density, $\rho_{\star}$, as a constraint in the analysis of stars that host transiting exoplanets can significantly improve the precision of the planet radius estimate in cases where the light curve is too noisy to yield an accurate value of the transit impact parameter, e.g. the light curves of Earth-size planets orbiting in the habitable zone of Sun-like stars that will be obtained by the PLATO mission. I have compiled a sample of 36 solar-type stars for which analysis of high-quality light curves together with constraints on the orbital eccentricity yield mean stellar density measurements with a median error of 2.3%. Of these, 8 are in transiting exoplanet systems and 28 in eclipsing binary systems with very low mass companions that contribute <0.1% of the total flux in the V band. A re-calibrated empirical relation for stellar mass as a function of T$_{\rm eff}$, $\rho_{\star}$ and [Fe/H] has been used to find mass estimates with a typical precision of 5.2% for the stars in this sample. Examples are given of how this sample can be used to test the accuracy and precision of $\rho_{\star}$ and log(g) estimates from catalogues of stellar parameters for solar-type stars.

Gamma ray burst (GRB) afterglow light curves have the potential to inform us about presently unobserved stages in the aftermath of a neutron star merger. Using numerical simulations of short GRB afterglows we obtain an approximate quantitative connection between key aspects of the emission mechanism and the shapes of the resulting light curves. Employing simple, but efficient, parameterizations of the light curve based on a broken power law in terms of physical parameters, fitted to a large dataset of synthetic light curves, we apply basic machine learning techniques to determine the approximate connection between key input parameters of the forward shock model and the light curve parameters. Solving then the inverse problem, we find that the strength of the central engine can be reasonably accurately estimated even with very limited information. In particular, merely the position of jet-break in the on-axis, respectively the maximum in the off-axis light curve determines the kinetic energy at the tens of percent level.

Direct detection experiments have established the most stringent constraints on potential interactions between particle candidates for relic, thermal dark matter and Standard Model particles. To surpass current exclusion limits a new generation of experiments is being developed. The upcoming upgrade of the CRESST experiment will incorporate $\mathcal{O}$(100) detectors with different masses ranging from $\sim$2g to $\sim$24g, aiming to achieve unprecedented sensitivity to sub-GeV dark matter particles with a focus on spin-independent dark matter-nucleus scattering. This paper presents a comprehensive analysis of the planned upgrade, detailed experimental strategies, anticipated challenges, and projected sensitivities. Approaches to address and mitigate low-energy excess backgrounds $-$ a key limitation in previous and current sub-GeV dark matter searches $-$ are also discussed. In addition, a long-term roadmap for the next decade is outlined, including other potential scientific applications.

The magneto-rotational instability (MRI) is widely believed to play a central role in generating large-scale, poloidal magnetic fields during binary neutron star mergers. However, the few simulations that begin with a weak seed magnetic field and capture its growth until saturation predominantly show the effects of small-scale turbulence and winding, but lack convincing evidence of MRI activity. In this work, we investigate how the MRI is affected by the complex magnetic field topologies characteristic of the post-merger phase, aiming to assess the actual feasibility of MRI in such environments. We first derive the MRI instability criterion, as well as expressions for the characteristic wavelength and growth timescale of the fastest-growing modes, under conditions that include significant magnetic field gradients. Our analysis reveals that strong radial magnetic field gradients can impact significantly on the MRI, slowing its growth or suppressing it entirely if large enough. We then apply this extended framework to both idealized analytical disk models and data from a numerical relativity simulation of a long-lived neutron star merger remnant. We find that conditions favourable to MRI growth on astrophysically relevant timescales may occur only in limited regions of the post-merger disk, and only at late times $t\gtrsim 100$ ms after the merger. These results suggest that the MRI plays a limited role in amplifying poloidal magnetic fields in post-merger environments during the first $\mathcal{O}(100)$ms.

Massive neutrinos imprint distinctive signatures on the evolution of cosmic structures, notably suppressing small-scale clustering. We investigate the impact of massive neutrinos on the galaxy bispectrum in redshift-space, adopting a spherical harmonic multipole decomposition $B_L^m(k_1, \mu, t)$, that captures the full angular dependence. We develop an analytical and numerical framework incorporating neutrino-corrected perturbation theory kernels and redshift-space distortions. Our results demonstrate that the linear triangle configurations are particularly sensitive to massive neutrinos, with deviations reaching up to $\sim 2\%$ for a total mass $\sum m_\nu = 0.12\,\mathrm{eV}$. To assess detection prospects in galaxy surveys like \textit{Euclid}, we compute the signal-to-noise ratio (SNR) for individual multipoles, including the effects of Finger-of-God damping and shot noise. The neutrino-induced signatures in $B_0^0$ and $B_2^0$ are found to be detectable with SNR $\gtrsim 5$ across a range of configurations, even after accounting for small-scale suppression. Higher-order multipoles such as $B_2^1$ and $B_2^2$ are moderately sensitive, with SNR $\gtrsim$ ($2-3$) in squeezed limits, while hexadecapole moments are more suppressed but still exhibit measurable signals at high $k_1$. Additionally, the SNR generally increases with wave number $k_1$, particularly for squeezed and stretched triangles, suggesting that access to smaller scales significantly enhances detection this http URL study highlights the potential of the redshift-space bispectrum multipoles as sensitive probes of massive neutrinos, complementing traditional power spectrum analyses, and underscores the importance of angular information and higher-order statistics for galaxy surveys.

Type-I thermonuclear bursts (TNBs) from neutron star low-mass X-ray binaries (NSLMXBs) originate on the neutron star's surface from the unstable burning of the accreted material. On the other hand, kHz quasi-periodic oscillations (QPOs) are thought to originate in the innermost regions of the in-spiraling accretion disk. Due to the violent nature of the bursts, it is anticipated that Type-I TNBs will impact the inner accretion flow, and, therefore, the kHz QPOs. In this work, we systematically study the evolution of the upper and the lower kHz QPOs immediately before and after a Type-I TNB on 4U 1636$-$536 using {\it{AstroSat}} observations in the 3--20 keV band. The analysis of the power-density-spectra show the presence of kHz QPOs within 200 seconds before the onset of the Type--I burst. However, we have not detected any prominent signature of the same within 100-200 second posterior to the burst. The kHz QPOs then re-emerges after $\approx$ 200 sec. The fractional rms variation in the 3--20 keV band drops by $\approx$ 5-6 \%, supporting the non-existence of kHz QPOs in the 200 sec post-Burst Zone. The time scale of 200 sec coincides with the viscous time scale, highlighting a scenario where the inner disk is temporarily disrupted by the intense radiation from the Type-I TNB. The kHz QPO then re-establishes as the inner disk is restored.

In this work, we explore how the size and surrounding tidal fields of dark matter protohalos at high redshift influence their angular momentum (AM) evolution. While tidal torque theory (TTT) states that AM arises from the misalignment between protohalo shape and tidal fields, it remains unclear what is the characteristic scale of the perturbations that couple with each protohalo, and its correlation with protohalo properties such as size. Moreover, although the assumptions of the TTT are assumed to hold during the linear and quasi-linear regime, cosmological simulations reveal that discrepancies between its predictions and the true AM of halos emerge earlier than expected. To address this, we analyze cosmological simulations to study tidal fields at $z=80$ using different smoothing lengths, and determine which best predicts AM under TTT. We then investigate discrepancies between predicted and actual AM across redshifts, considering the effect of evolving tidal and inertia tensors. Our results show that the early tidal field couples with the inertia tensor of protohalos on scales about half of their characteristic size and confirm that disagreements between theory and simulation arise before the non-linear regime, suggesting a systematic effect from protohalo shape interacting with the forming cosmic web.

The problem of quasar classification comes in the class of highly imbalanced classification problems since Radio-loud (RL) quasars are rare and make up only about 10% of the quasar population. In this work, we use the Sloan Digital Sky Survey-DR3 dataset and introduce a PCA-based regression pipeline designed to maximize recall for rare classes in class-imbalanced astronomical data. We demonstrate an effective methodology to identify the key features of the dataset and apply Principal Component Analysis (PCA) for dimensionality reduction. For the PCA transformed SDSS-DR3 dataset, first two components account for the 97% of the observed variance. We perform classification of Radio-Loud (RL) and Radio-Quiet (RQ) quasars with Random Forest Classifier (RFC), novel PCA based balanced linear regression classifier (PBC), Random forest integrated with SMOTE classifier and XGBoost classifier with threshold tuning. RFC achieves an overall accuracy of 92% while PBC achieves an overall accuracy of 62%. XGBoost achieves an overall accuracy of 72% and SMOTE integrated RFC achieves an accuracy of 85%. Higher precision is obtained for RQ quasars in all classification methods. For the RL class, RFC achieves a recall of 0.04, XGBoost achieves a recall of 0.39, SMOTE integrated RFC achieves a recall of 0.25 and PBC achieves a recall of 0.52 attributed to the balanced logistic regression. RFC and PBC achieve F1 score of 0.08 and 0.19 respectively for RL while XGBoost achieves an improved F1 score of 0.22 but at the cost of reduced recall of the RL class. SMOTE integrated RFC achieves a better F1 score of 0.21 over RFC and PBC. Overall results of classifiers point to extreme class imbalance between RQ and RL classes in the data set.

The gravitational wave signal of binary compact objects has two main contributions at frequencies below the characteristic merger frequency: the gravitational wave signal associated with the early inspiral stage of the binary and the non-linear gravitational wave memory. We compare the sensitivity of upcoming gravitational wave detectors to these two contributions, with a particular interest in events with a merger phase at frequencies higher than the detector's peak sensitivity. We demonstrate that for light primordial black holes, current and upcoming detectors are more sensitive to the inspiral signal. Our analysis incorporates the evolution history of primordial black hole binaries, key to accurately estimating the relevant event rates. We also discuss the waveform templates of the memory signal at ground- and space-based interferometers, and the implications for a matched filtering search. This allows us to compare the sensitivity of high-frequency gravitational wave detectors, sensitive to the merger phase, with the sensitivity of existing interferometers.

The mixing of material from stellar convective cores into their adjacent radiative layers has been a matter of long-standing debate. Pulsating subdwarf B stars offer excellent conditions to advance our understanding of this problem. In this work we use a model-independent approach to infer information about the cores of three subdwarf B stars and compare it with similar inferences from earlier analysis of red giants in the helium core-burning phase. This is achieved by fitting an analytical description of the gravity-mode pulsation periods to pulsation data collected by the Kepler satellite. From the fits we infer the reduced asymptotic period spacings and the amplitude and position of sharp structural variations associated with chemical discontinuities in the stellar interiors. Our results indicate the presence of sharp structural variations with similar properties in all three stars, located near the edge of the gravity-mode propagation cavity and likely associated with the C-O/He transition. We find that these structural variations differ systematically from those of helium core-burning red giant stars, having larger amplitudes and being located at a larger buoyancy radius. This suggests that chemical mixing beyond the adiabatically stratified core into the radiatively stratified layers may be more extensive in subdwarf B stars than in helium core-burning red giants. Alternatively, the stratification of the mixing region beyond the adiabatically stratified core may differ significantly between the two types of stars. The model-independent constraints set on the structural variations inside these three stars are the first of a kind and will be key to enhance the modelling of layers adjacent to stellar convective cores and to test non-canonical stellar evolution channels leading to the formation of hot subdwarf stars.

We perform high-resolution atmospheric flow simulations of hot and warm giant exoplanets that are tidally locked. The modeled atmospheres are representative of those on KELT-11b and WASP-39b, which possess markedly different equilibrium temperatures but reside in a similar dynamical regime: in this regime, their key dynamical numbers (e.g., Rossby and Froude numbers) are comparable. Despite their temperature difference, both planets exhibit qualitatively similar atmospheric circulation patterns, which are characterized by turbulent equatorial flows, anticyclonic polar vortices, and large-scale Rossby waves that gives rise to quasi-zonal flows in the extra-tropics (i.e., near ~20 degrees). Quantitative differences between the KELT-11b and WASP-39b atmospheres reflect their different Rossby deformation scales, which influence the horizontal length scale of wave--vortex interactions and the overall structure of the circulation.

We study the evolution of the bar fraction in disc galaxies between $0.5 < z < 4.0$ using multi-band coloured images from JWST CEERS. These images were classified by citizen scientists in a new phase of the Galaxy Zoo project called GZ CEERS. Citizen scientists were asked whether a strong or weak bar was visible in the host galaxy. After considering multiple corrections for observational biases, we find that the bar fraction decreases with redshift in our volume-limited sample (n = 398); from $25^{+6}_{-4}$% at $0.5 < z < 1.0$ to $3^{+6}_{-1}$% at $3.0 < z < 4.0$. However, we argue it is appropriate to interpret these fractions as lower limits. Disentangling real changes in the bar fraction from detection biases remains challenging. Nevertheless, we find a significant number of bars up to $z = 2.5$. This implies that discs are dynamically cool or baryon-dominated, enabling them to host bars. This also suggests that bar-driven secular evolution likely plays an important role at higher redshifts. When we distinguish between strong and weak bars, we find that the weak bar fraction decreases with increasing redshift. In contrast, the strong bar fraction is constant between $0.5 < z < 2.5$. This implies that the strong bars found in this work are robust long-lived structures, unless the rate of bar destruction is similar to the rate of bar formation. Finally, our results are consistent with disc instabilities being the dominant mode of bar formation at lower redshifts, while bar formation through interactions and mergers is more common at higher redshifts.

In this paper, we investigate the impact of charge on the stability of the pulsar star SAX J1748.9-2021 within the context of $f(\mathbb{Q}, \mathcal{T})$ theory, where $\mathbb{Q}$ and $\mathcal{T}$ represent the non-metricity scalar and the energy-momentum tensor, respectively. To achieve this, we employ the Krori-Barua metric ansatz with anisotropic fluid. We obtain exact relativistic solutions for the corresponding field equations. Additionally, we examine its physical and geometric properties using astrophysical observations of the pulsar SAX J1748.9-2021. This approach provides a basis for connecting several physical quantities, including fluid parameters, anisotropy, mass-radius relationship, compactness, redshift, \emph{energy}, the Zeldovich and causality conditions, equation of state parameter, adiabatic index and the Tolman-Oppenheimer-Volkoff equation. Our findings align with observational evidence, indicating that the pulsar SAX J1748.9-2021 remains both feasible and stable under this modified theory.

We study the purely gravitational signatures of dark matter from the ultralight to the ultraheavy mass range in proposed long-baseline atom gradiometers. Due to their exceptional acceleration sensitivity, proposed space-based experiments, such as MAGIS-space, AEDGE and AEDGE+, could detect a dark matter subcomponent which constitutes $\mathcal{O}(10\%)$ of the local dark matter energy density and is populated by compact clumps of mass between $10^6$ kg and $10^{10}$ kg ($10^{-25}~M_\odot\lesssim M \lesssim 10^{-21}~M_\odot$), in an otherwise unexplored region of dark matter model space. Furthermore, because the gravitational observable depends on the relative gravitational time delay measured by spatially separated atomic clouds, we find that atom gradiometers are parametrically more sensitive than laser interferometers, such as LIGO and LISA, to fast-oscillating spacetime perturbations sourced by energy density and pressure fluctuations of ultralight dark matter. In optimistic scenarios, space-based experiments are projected to probe a DM overdensity of $\mathcal{O}(10)$ times the local dark matter energy density for masses $m\lesssim 10^{-17}$ eV. We therefore conclude that proposed atom gradiometers may be poised to detect dark matter through purely gravitational interactions.

We present a novel type of soliton dubbed soft oscillons. In contrast with conventional oscillons the soft counterparts come in a continuum of unboundedly large sizes. They are peculiar also in that the oscillation frequency is set by their size and in their high harmonic content, that leads to a peculiar core structure and pulsating emission of radiation. Soft oscillons appear in a broad class of scalar field theories with plateau potentials and owe their existence to a simple confinement mechanism: a large field perturbation automatically acts as a cavity that traps massless modes. We study the classical evaporation rate, lifetime and stability under anisotropic deformations. Soft oscillons are longer lived in gapless models, where lifetimes can reach 100 times the initial radius or more.

We review the dynamics of spectator scalar fields non-minimally coupled to gravity in the post-inflationary Universe, with particular emphasis on scenarios where the end of inflation is followed by a period of kination. In this context, the evolution of the Ricci scalar can lead to the spontaneous breaking of discrete or continuous symmetries through tachyonic instabilities. These Hubble-induced phase transitions may cause the amplification of field fluctuations, the formation of transient topological defects, and an efficient energy transfer into relativistic degrees of freedom. We analyze this general mechanism and its cosmological consequences, including (re)heating and gravitational wave production. As a concrete realization, we discuss the postinflationary evolution of the Standard Model Higgs, examining the interplay between curvature effects, vacuum stability, and non-perturbative dynamics. This framework provides a minimal and predictive connection between high-energy physics and the early Universe, with potential observational signatures.

In the present review, we consider the status of the classification of the vacuum, stationary and asymptotically flat black holes in scalar-tensor gravity. Contrary to the similar problem in general relativity, the black hole classification in scalar-tensor theories is much more challenging due to the very complicated character of the field equations and the very complex mathematical structure of the scalar-tensor gravity as a whole. We review most of the known no-hair results, and where possible new ones, as well demonstrate some of the difficulties that appear in our attempts to classify the black holes within scalar-tensor gravity. The proofs of the theorems and the underlying mathematical techniques are given in sufficient detail. To make the review self-contained we also present the vacuum black hole uniqueness theorems in general relativity and their proofs.

We compute by means of post-Newtonian (PN) methods the innermost stable circular orbit (ISCO) of arbitrary-mass (in particular, comparable-mass) compact binaries. Two methods are used with equivalent results: dynamical perturbation of the conservative equations of motion in harmonic coordinates, and dynamical perturbation of the conservative Hamiltonian in ADM coordinates. The perturbation of the non-local tail term at 4PN order in both approaches is carefully investigated. Our final gauge invariant result for the location of the ISCO at 4PN order is close to the numerical value of the ISCO shift computed by the gravitational self-force (GSF) approach in the small mass-ratio limit, and is also in good agreement with the full numerical-relativity calculation in the case of equal masses. The PN method followed here is considered in standard Taylor-expanded form, without any resummation techniques applied. As a complement, we also compute explicitly the gauge transformation from harmonic coordinates to ADM coordinates up to 4PN order, including the tail contribution therein.

We present an exact analytical investigation of spin precession for a test gyroscope in the magnetized Kerr spacetime--an exact electrovacuum solution to the Einstein-Maxwell equations. Our approach accommodates arbitrary magnetic field strengths, enabling a unified treatment across both weak and ultra-strong field regimes. The analysis reveals distinct spin precession behaviors near rotating collapsed objects, which differ characteristically between black holes and naked singularities, offering a potential observational means to differentiate them. The external magnetic field induces a nontrivial modification of the precession frequency through its interaction with the spacetime's gravitoelectromagnetic structure. In the weak-field limit, magnetic fields generally reduce the precession rate, though the effect depends sensitively on the motion and orientation of the test gyro close to the collapsed object. As a special case, we show that in the presence of magnetic fields, the spin precession frequency due to gravitomagnetic effect acquires a long-range $1/r$ (where $r$ is the distance from the central object to the test gyro) correction in contrast to the standard $1/r^3$ falloff. In addition, we obtain the exact geodetic precession (gravitoelectric effect) frequency in magnetized Schwarzschild spacetime, showing that the magnetic field enhances ($\propto r^{1/2}$) geodetic precession in contrast to the standard $1/r^{5/2}$ falloff. Our results provide observationally testable predictions relevant for black holes in strong magnetic environments, including those possibly realized near magnetars or in the early universe. In particular, the strong-field behavior of spin precession could have important implications for transmuted black holes formed via collapse or mergers of magnetized progenitors in both astrophysical and cosmological contexts.