2025-07-14 | | Total: 114
A variety of mechanisms in the early Universe lead to the generation of gravitational waves (GWs). We introduce here a novel source of GWs generated by vacuum fluctuations after inflation. Being gravitons minimally coupled particles, their quantum creation takes place during inflation, but is absent in an unperturbed Universe during the radiation-dominated epoch, since they behave as conformally invariant particles. However, the presence of inhomogeneities breaks the conformal flatness of the metric, allowing scalar metric perturbations to induce the quantum production of gravitons. We compute the resulting GW spectrum from this mechanism for different models of the primordial scalar power spectrum. We find that this GW signal peaks around the GHz frequency range, distinguishing it from other astrophysical and cosmological backgrounds and underscoring the need for detectors sensitive to these high frequencies.
Hubble tension can be alleviated by altering either early- or late-time ΛCDM. So far, extensive studies have shown that only early dark energy or ad hoc combinations of those two-fold effects can reduce the tension to 3σ or lower. In this work, we improve the later solution by considering a cascade decaying dark matter sector, where a parent particle produces relativistic dark matter in the early Universe and the dark matter subsequently decays to neutrinos in the late-time Universe. We parametrize this model with four model parameters, carry out a Markov Chain Monte Carlo fit to Planck 2018+BAO+Pantheon data sets, and compare it to Big Bang Nucleosynthesis limits, neutrino flux bounds, and Planck and SDSS data about matter power spectra. Within parameter regions reducing the Hubble tension to 3.2σ and compatible with the existing constraints, the main phenomenological features include a larger S8∼0.84, smaller ns∼0.95 and larger matter power spectra relative to the ΛCDM, which are left for future tests.
In this work, we investigated the influence of constants that characterize the motion of thin tubes of a massless scalar field (TToMSF) in a gas on the global properties of such a gas, in particular, on its ability to accelerated expansion. The peculiarity of the TToMSF gas is related to the fact that the parameter of its equation of state w depends on the degree of rarefaction of the gas and can take values in a very wide range of values −1<w<1. In this case, fluctuations in the density of the TToMSF gas, caused by its gravitational interaction with non-uniformly distributed baryonic matter, should lead to the implementation of local equations of state with different values and different laws of change of the parameter w. The possibility of implementing local equations of state can alleviate the tensions associated with the observed non-uniformity of the distribution of matter and the growth of structures in the Universe.
As ground-based gravitational-wave (GW) detectors improve in sensitivity, gravitational-wave background (GWB) signals will progressively become detectable. Currently, searches for the GWB model the signal as a power law, however deviations from this model will be relevant at increased sensitivity. Therefore, to prepare for the range of potentially detectable GWB signals, we propose an interpolation model implemented through a transdimensional reverse-jump Markov chain Monte Carlo (RJMCMC) algorithm. This interpolation model foregoes a specific physics-informed model (of which there are a great many) in favor of a flexible model that can accurately recover a broad range of potential signals. In this paper, we employ this framework for an array of GWB applications. We present three dimensionless fractional GW energy density injections and recoveries as examples of the capabilities of this spline interpolation model. We further demonstrate how our model can be implemented for hierarchical GW analysis on ΩGW.
Quantum effects are expected to modify the cosmological dynamics of the early universe while maintaining some (potentially discrete) notion of space-time structure. In one approach, loop quantum cosmology, current models are shown here to either be incompatible with a consistent space-time structure, or to have physical singularities. The latter happens in spite of a non-zero scale factor in the isotropic background dynamics. A new effective Friedmann equation shows that a bounce is obtained at sub-Planckian densities, preceded by a physical singularity at infinite scale factor that resembles a time-reversed big rip. The entire phase is accompanied by rapid changes of the Hubble radius. In addition, a new version of perturbative inhomogeneity in loop quantum cosmology is introduced that maintains a consistent space-time structure and has a non-singular background dynamics.
Accurate and efficient maneuver detection is critical for ensuring the safety and predictability of spacecraft trajectories. This paper presents a novel maneuver detection approach based on comparing the confidence levels associated with the orbital state estimation and the observation likelihood. First, a confidence-dominance maneuver indicator (CDMI) is proposed by setting a confidence level for the state estimation and computing the maximum likelihood of the observation and its confidence level. The CDMI then flag a maneuver when the observation's confidence level exceeds that of the state estimation, indicating that the observation is unlikely under the no-maneuver hypothesis while maintaining consistency with the prior state estimation confidence. To efficiently compute the maximum likelihood of the observation and obtain the CDMI, a recursive polynomial optimization method is developed, taking advantage of convex optimization and polynomial approximation. In addition, an integrated CDMI approach is developed to eliminate the need to manually select the state confidence level. The integrated CDMI approach maintains high detection accuracy while simultaneously providing an indication of maneuver likelihood, thereby enhancing robustness and practical applicability. The performance of the proposed CDMI-based maneuver detection approaches is evaluated against an optimal control distance metric and two mixture-based approaches. The simulation results demonstrate that the proposed integrated CDMI approach can achieve up to 99.33\% detection accuracy, at least 10% higher than the competing methods, while substantially reducing computational costs.
We introduce a machine learning (ML) framework called TIER for improving the sensitivity of gravitational wave search pipelines. Typically, search pipelines only use a small region of strain data in the vicinity of a candidate signal to construct the detection statistic. However, extended strain data (∼10 s) in the candidate's vicinity can also carry valuable complementary information. We show that this information can be efficiently captured by ML classifier models trained on sparse summary representation/features of the extended data. Our framework is easy to train and can be used with already existing candidates from any search pipeline, and without requiring expensive injection campaigns. Furthermore, the output of our model can be easily integrated into the detection statistic of a search pipeline. Using TIER on triggers from the IAS-HM pipeline, we find up to ∼20% improvement in sensitive volume time in LIGO-Virgo-Kagra O3 data, with improvements concentrated in regions of high masses and unequal mass ratios. Applying our framework increases the significance of several near-threshold gravitational-wave candidates, especially in the pair-instability mass gap and intermediate-mass black hole (IMBH) ranges.
Standardizing the definition of eccentricity is necessary for unambiguous inference of the orbital eccentricity of compact binaries from gravitational wave observations. In previous works, we proposed a definition of eccentricity for systems without spin-precession that relies solely on the gravitational waveform, is applicable to any waveform model, and has the correct Newtonian limit. In this work, we extend this definition to spin-precessing systems. This simple yet effective extension relies on first transforming the waveform from the inertial frame to the coprecessing frame, and then adopting an amplitude and a phase with reduced spin-induced effects. Our method includes a robust procedure for filtering out spin-induced modulations, which become non-negligible in the small eccentricity and large spin-precession regime. Finally, we apply our method to a set of Numerical Relativity and Effective One Body waveforms to showcase its robustness for generic eccentric spin-precessing binaries. We make our method public via Python implementation in \texttt{gw\_eccentricity}.
Recently, we introduced the Lorentzian-Euclidean black hole, a static and spherically symmetric solution of vacuum Einstein equations that exhibits a change in metric signature across the event horizon. In this framework, the analysis of radial trajectories of freely falling bodies proves that the central singularity can be avoided via a mechanism we interpret as atemporality, which is responsible for the shift of the time variable from real to imaginary values. In this paper, we further explore this model by first examining the behavior of null geodesics. Our investigation requires a set of signature-adaptive coordinate changes that generalize the local Lorentz transformations underlying General Relativity. We find that photon orbits, like their massive counterparts, cannot traverse the event horizon, thereby strengthening the previous result on the impossibility to reach the r=0 singularity. Additionally, we discuss the causal structure of the spacetime, provide the corresponding Penrose diagram, and analyze the process of matter accretion in the outer region of the black hole.
In warm inflation (WI), the persistent thermal bath that is sustained by dissipative interactions with the inflaton field produces a stochastic background of gravitational waves (GWs). In this paper we study the production and evolution of these GWs. Specifically, we investigate the emission of thermal gravitons (gravitons emitted by a thermal bath) from particle scattering in the bath and the evolution of the corresponding GWs. We find that the bulk of thermal graviton production in WI occurs during the transition to radiation domination after inflation. Further, the energy density of thermal gravitons is enhanced by roughly one to two orders of magnitude compared to that in a radiation-dominated scenario with the same reheating temperature. We also calculate the spectrum of the resulting stochastic GW background and find that it has a distinctive shape, consisting of a peak at high frequencies ~100 GHz and an almost flat spectrum extending to low frequencies. The peak arises from emission of sub-horizon modes that follow the temperature of the bath. The flat part of the spectrum corresponds to the modes that exit the horizon during WI and re-enter during radiation domination. We show that the detection prospects for the high-frequency peak of the GW spectrum, while improved slightly compared to the radiation-dominated case, still remain challenging. The thermal spectrum's low-frequency plateau is typically subdominant to the amplitude of the standard vacuum tensor modes from inflation, although WI models can exist where the thermal graviton plateau surpasses the vacuum contribution without exceeding current observational limits on the tensor-to-scalar ratio. Furthermore, we calculate the thermal graviton contribution from WI to dark radiation and show that WI models are generally expected to satisfy current observational bounds, including those from the cosmic microwave background.
On January 14, 2025 the LIGO interferometers detected a gravitational wave from the merger of two black holes, GW250114. Using publicly available information, we estimate that the signal-to-noise ratio (SNR) of GW250114 was ∼80. This would make it three to four times louder than any other gravitational wave detected to date. GW250114 therefore offers a unique opportunity to make precise measurements of its source parameters and to test general relativity. In anticipation of its public data release, we analyze a set of simulated signals that have parameters similar to what we estimate for GW250114 and explore what new insights may be gained from this significant event.
Deep imaging of structures from the Cosmic Dawn (CD) and the Epoch of Reionization (EoR) in five targeted fields is one of the highest priority scientific objectives for the Square Kilometre Array (SKA). Selecting 'quiet' fields, which allow deep imaging, is critical for future SKA CD/EoR observations. Pre-observations using existing radio facilities will help estimate the computational capabilities required for optimal data quality and refine data reduction techniques. In this study, we utilize data from the Murchison Widefield Array (MWA) Phase II extended array for a selected field to study the properties of foregrounds. We conduct deep imaging across two frequency bands: 72-103 MHz and 200-231 MHz. We identify up to 2,576 radio sources within a 5-degree radius of the image center (at RA (J2000) 8h , Dec (J2000) 5{\deg}), achieving approximately 80% completeness at 7.7 mJy and 90% at 10.4 mJy for 216 MHz, with a total integration time of 4.43 hours and an average RMS of 1.80 mJy. Additionally, we apply a foreground removal algorithm using Principal Component Analysis (PCA) and calculate the angular power spectra of the residual images. Our results indicate that nearly all resolved radio sources can be successfully removed using PCA, leading to a reduction in foreground power. However, the angular power spectra of the residual map remains over an order of magnitude higher than the theoretically predicted CD/EoR 21 cm signal. Further improvements in data reduction and foreground subtraction techniques will be necessary to enhance these results.
This paper is devoted to study the viability of DESI and fast radio burst to constraint free parameters of modified theories of gravity. Thus, we present a model supported in f(R) gravity involving a function of the Ricci scalar named Starobinsky-type with the peculiarity that the Non-commutative essence is intrinsic to the coefficients. Additionally, to understand the dynamics within a flat Friedmann-Lemaître-Robertson-Walker Universe, we explore the possibility of deriving a Friedman equation (in measure) that results from an adequate mathematical treatment. As we mention previously, to test the outlined model, a Monte Carlo Markov Chain analysis is implemented, using Cosmic Chronometers, type Ia supernovae, Hydrogen II galaxies, Intermediate-luminosity quasars, Baryon Acoustic oscillations and Fast Radio Bursts data, to constraint the free parameters of the model and presenting H(z), q(z) and ωeff(z). The final results are compared with the ΛCDM model and it is presented a robust discussion about the viability of DESI and fast radio burst to constraint free parameters in specific of a modified theory of gravity.
Although the field of neutrino astronomy has blossomed in the last decade, physicists have struggled to fully map the high-energy neutrino sky. TAMBO, a mountain-based neutrino observatory, aims to solve that issue -- and find clues of new physics along the way.
The identification of stellar structures in the Galactic halo, including stellar streams and merger remnants, often relies on the dynamics of their constituent stars. However, this approach has limitations due to the complex dynamical interactions between these structures and their environment. Perturbations such as tidal forces exerted by the Milky Way, the potential escape of stars, and passages through the Galactic plane can result in the loss of dynamical coherence of stars in these structures. Consequently, relying solely on dynamics may be insufficient for detecting such disrupted or dispersed remnants. We combine chemistry and dynamics, integrated through a system of neural networks, to develop a clustering method for identifying accreted structures in the Galactic halo. We developed an integrated approach combining Siamese neural networks (SNNs), graph neural networks (GNNs), autoencoders, and the OPTICS algorithm to create a comprehensive procedure named CREEK. This method is designed to uncover stellar structures in the Galactic halo. Initially, CREEK was trained on known globular clusters (GCs) and then applied to the dataset to identify stellar streams. CREEK successfully recovered 80% of the GCs present in the APOGEE dataset, re-identified several known stellar streams, and identified a potential new stream. Additionally, within highly populated stellar structures, CREEK can identify substructures that exhibit distinct chemical compositions and orbital energies. This approach provides an objective data-driven method for selecting stars associated with streams and stellar structures in general.
Cosmic inflation provides a compelling framework for explaining several observed features of our Universe, but its viability depends on an efficient reheating phase that converts the inflaton's energy into Standard Model particles. This conversion often proceeds through non-perturbative mechanisms such as parametric resonance, which is described by Hill's equation. In this work, we investigate how stochastic fluctuations in the parameters of Hill's equation can influence particle production during reheating. We show that such fluctuations can arise from couplings to light scalar fields, and can significantly alter the stability bands in the resonance structure, thereby enhancing the growth of fluctuations and broadening the region of efficient energy transfer. Using random matrix theory and stochastic differential equations, we decompose the particle growth rate into deterministic and noise-induced components and demonstrate analytically and numerically that even modest noise leads to substantial particle production in otherwise stable regimes. These results suggest that stochastic effects can robustly enhance the efficacy of reheating across a wide swath of parameter space, with implications for early Universe cosmology, UV completions involving multiple scalar fields, and the resolution of the cosmological moduli problem.
We present a multi-phase study of the gas in the circumnuclear region (~1.1x1.0 kpc^2) of the nearby Seyfert 1.8 galaxy NGC 1365, observed in the context of the Mid-IR Activity of Circumnuclear Line Emission (MIRACLE) program. We combined spatially resolved spectroscopic observations from JWST/MIRI, VLT/MUSE, and ALMA to investigate the ionized atomic gas and the warm and cold molecular phases. MIRI data revealed over 40 mid-IR emission lines from ionized and warm molecular gas. Moment maps show that both cold and warm molecular gas follow the rotation of the stellar disk along the circumnuclear ring. The ionized gas displays flux and kinematic patterns that depend on ionization potential (IP): low-IP species (<25 eV) trace the disk, while higher-IP lines (up to ~120 eV) trace outflowing material. The [O III]5700 and [Ne V]14 lines both trace the southeast nuclear outflow cone. Additionally, [Ne V]14 detects the northwest counter-cone, obscured in the optical and thus invisible in [O III]5700. Mid-IR diagnostics, unlike optical ones, clearly reveal the AGN as the primary ionization source in the nucleus. Emission from high-IP species is spatially coincident with the ionization cones and not with star-forming regions. Using the [Ne V]24/[Ne V]14 ratio, we derive an electron density of (750+-440) cm^(-3), in agreement with values from the [S II] optical doublet. For the first time, we apply a fully self-consistent approach combining advanced photoionization and kinematic models (HOMERUN+MOKA3D) to constrain intrinsic outflow properties, overcoming the limitations of simplified classical methods. Exploiting the synergy of JWST/MIRI and VLT/MUSE, HOMERUN reproduces fluxes of over 60 emission lines from optical to mid-IR, disentangling AGN and star formation contributions and yielding robust estimates of outflow mass, geometry, and energetics.
Supernova 2023ixf occurred on May 18th 2023 in the nearby galaxy Messier 101 (D ~ 6.85 Mpc), making it the closest supernova in the last decade. Following its discovery, astronomers around the world rushed to observe the explosion across the electromagnetic spectrum in order to uncover its early-time properties. Based on multi-wavelength analysis during its first year post-explosion, supernova 2023ixf is a type II supernova that interacted with dense, confined circumstellar material in its local environment -- this material being lost from its red supergiant progenitor in the final years before explosion. In this article, we will review the findings of >80 studies already published on this incredible event as well as explore how the synthesis of SN 2023ixf observations across the electromagnetic spectrum can be used to constrain both type II supernova explosion physics in addition to the the uncertain mass loss histories of red supergiant stars in their final years.
Central to our understanding of stellar evolution and its impact on processes in our Galaxy and across the Universe is the study of mass loss. While the general framework is well established, recent JWST observations of objects like NGC 3132 have revealed intricate nebular structures, suggesting complex mass-loss processes likely driven by multiple star system at its core. These findings pose new challenges for the currently available investigation tools. The primary goal of this study is the first detailed comparison of the physical properties and chemical composition obtained for NGC 3132, based on the latest detailed 3D model and observations from MUSE, JWST and Spitzer. We evaluate the reliability of the traditional empirical method and photoionization model for abundances estimations, both based on the same available high-quality, spatially resolved observations. We find that the model and empirical method yield consistent results for the integrated total properties such as Te, ne and chemical abundances. However, when applied to simulated observations from the model, the empirical method fails to recover the model input abundances, providing only an approximate estimate. This discrepancy arises in part from the loss of information when summing fluxes over regions which have complex ionisation structures. This discrepancy in the case of oxygen has been estimated to be up to 35%. Moreover, the latest IR data reveal a spatial correlation between H2, c(Hb) as well as the [8.0]/[4.5] IRAC ratio. Finally, new clumps are discovered in [Ni II] 7378 \AA, [Fe II] 8617 Åand [Fe III] 5270 Åemission lines.
Tidal disruption events (TDEs) offer a unique probe of supermassive black hole (SMBH) demographics, but their observed rates remain difficult to reconcile with standard single-SMBH models. In this work, we use simulations of SMBH binaries, including the combined effects of eccentric Kozai-Lidov oscillations and two-body relaxation, to explore how TDE rates scale with SMBH mass and redshift. We find that binary systems exhibit increasing TDE rates with mass, in contrast to the declining trend expected for single SMBHs. These binary-driven rates match those observed in post-starburst galaxies, suggesting that a subset of TDE hosts may contain SMBH binaries. TDE light curves in some massive galaxies exhibit unexpectedly short durations, suggesting that the disrupting SMBH may be less massive than implied by host galaxy scaling relations, consistent with disruptions by the less massive black hole in a binary. By convolving our mass-dependent rates with the SMBH mass function, we predict redshift-dependent TDE rates, which we show can be used to constrain the supermassive black hole binary fraction. Our results provide a testable framework for interpreting TDE demographics in upcoming wide-field surveys such as LSST and Roman.
Line broadening and variability are observational hallmarks of active galactic nuclei which allow us to measure supermassive black hole masses as well as constrain the geometry and kinematics of the emitting gas, with the most precise measurements requiring a degree of modeling. Two popular models of the broad-line region describe the emitting gas as either a distribution of puffed up clouds or a thin disk with strong velocity gradients. As we show in this work, key features in the reverberation mapping dataset obtained by the AGN STORM team in NGC 5548 cannot be accounted for by either simple model. In several emission lines the observed broad-line profile has a single peak yet the delay profile has a distinct double peak, strongly motivating a BLR with emission from multiple components. We demonstrate a few possibilities that may alleviate the tension and better represent the true nature of the broad-line emitting gas in NGC 5548 and beyond.
One of the most commonly observed solar radio sources in the metric and decametric wavelengths is the solar noise storm. These are generally associated with active regions and are believed to be powered by the plasma emission mechanism. Since plasma emission emits primarily at the fundamental and harmonic of the local plasma frequency, it is significantly affected by density inhomogeneities in the solar corona. The source can become significantly scatter-broadened due to the multi-path propagation caused by refraction from the density inhomogeneities. Past observational and theoretical estimates suggest some minimum observable source size in the solar corona. The details of this limit, however, depends on the modeling approach and details of the coronal turbulence model chosen. Hence pushing the minimum observable source size to smaller values can help constrain the plasma environment of the observed sources. In this work, we for the first time, use data from the upgraded Giant Metrewave Radio Telescope in the 250--500 MHz band, to determine multiple instances of very small-scale structures in the noise storms. We also find that these structures are stable over timescales of 15--30 minutes. By comparing the past observations of Type III radio bursts and noise storms, we hypothesize that the primary reason behind the detection of these small sources in noise storm is due to the local environment of the noise storm. We also build an illustrative model and propose some conditions under which the minimum observable source size predicted by theoretical models, can be lowered significantly.
Gamma-ray emission from the plane of the Milky Way is understood as partly originating from the interaction of cosmic rays with the interstellar medium. The same interaction is expected to produce a corresponding flux of neutrinos. In 2023, IceCube reported the first observation of this galactic neutrino flux at 4.5σ confidence level. The analysis relied on neutrino flux predictions - based on gamma ray observations - to model the expected neutrino emission from the galactic plane. Three signal hypotheses describing different possible spatial and energy distributions were tested, where the single free parameter in each test was the normalization of the neutrino flux. We present first results of an analysis that can improve the characterization of Galactic neutrino emission by dividing the galactic plane into segments in galactic longitude. An unbinned maximum-likelihood analysis is used that can fit the spectral index and the flux normalization separately in each segment. While gamma ray telescopes can not differentiate between hadronic and leptonic emission, neutrino production must come from hadronic processes. Measuring a spectral index can further help to understand the contribution of unresolved neutrino sources inside the galactic plane. This work uses a full-sky cascade dataset and provides model-independent insight into the variation of the neutrino flux and energy distribution from different regions of the galactic plane.
We report the detection of near- and mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) out to ~ 35 kpc in the Makani Galaxy, a compact massive galaxy with a record-breaking 100-kpc scale starburst-driven wind at redshift z = 0.459. The NIRCam and MIRI observations with JWST take advantage of a coincidental match between the PAH spectral features at 3.3, 7.7, and (11.3 + 12.2) microns in Makani and the bandpasses of the MIRI and NIRCam filters. The warm dust is not only detected in the cool-gas tracers of the galactic wind associated with the more recent (7 Myr) starburst episode, but also in the outer warm-ionized gas wind produced by the older (0.4 Gyr) episode. The presence of PAHs in the outer wind indicates that the PAHs have survived the long (R/v ~ 10^8 yrs) journey to the halo despite the harsh environment of the galactic wind. The measured F1800W/F1130W flux ratios in the unresolved nucleus, inner halo (R = 10 - 20 kpc), and outer halo (R = 20 - 35 kpc), tracers of the PAH (11.3 + 12.2)/7.7 ratios, indicate decreasing starlight intensity incident on the PAHs, decreasing PAH sizes, and increasing PAH ionization fractions with increasing distance from the nucleus. These data provide the strongest evidence to date that the ejected dust of galactic winds survives the long journey to the CGM, but is eroded along the way.
The recent discovery of the third interstellar object (3I/ATLAS) expands the known census from two to three and significantly improves statistical inferences regarding the underlying galactic population. 3I/ATLAS exhibits detectable activity that may increase as the object approaches perihelion. In this paper, we argue that the cometary activity likely significantly contributes to 3I/ATLAS's brightness, since the nuclear size of 3I/ATLAS assuming an asteroidal reflectance implies an untenable interstellar object mass per star. 3I/ATLAS exhibits an excess velocity of v∞=58 km/s relative to the Sun which is significantly larger than that of 1I/`Oumuamua or 2I/Borisov. This velocity implies that 3I/ATLAS is relatively old in comparison to previous interstellar objects. Here, we calculate the posterior distribution of ages implied by the kinematics of the interstellar objects and find that 3I/ATLAS is likely ∼3−11 Gyr old, assuming that the interstellar object and stellar age-velocity relations are equivalent. We also calculate the distribution of host star metallicities and find that 3I/ATLAS likely originated in a lower-metallicity system than the previous interstellar objects. These results show that interstellar object formation is likely efficient at low metallicities and early in the history of the Galaxy. Finally, we estimate and approximate interstellar object formation rate throughout Galactic history implied by these three objects.