Search for: All records

ABSTRACT One-dimensional (1D) stellar evolution models are widely used across various astrophysical fields, however they are still dominated by important uncertainties that deeply affect their predictive power. Among those, the merging of independent convective regions is a poorly understood phenomenon predicted by some 1D models but whose occurrence and impact in real stars remain very uncertain. Being an intrinsically multi-D phenomenon, it is challenging to predict the exact behaviour of shell mergers with 1D models. In this work, we conduct a detailed investigation of a multiple shell merging event in a 20 M$$_\odot$$ star using 3D hydrodynamic simulations. Making use of the active tracers for composition and the nuclear network included in the 3D model, we study the merging not only from a dynamical standpoint but also considering its nucleosynthesis and energy generation. Our simulations confirm the occurrence of the merging also in 3D, but reveal significant differences from the 1D case. Specifically, we identify entrainment and the erosion of stable regions as the main mechanisms that drive the merging, we predict much faster convective velocities compared to the mixing-length theory velocities, and observe multiple burning phases within the same merged shell, with important effects for the chemical composition of the star, which presents a strongly asymmetric (dipolar) distribution. We expect that these differences will have important effects on the final structure of massive stars and thus their final collapse dynamics and possible supernova explosion, subsequently affecting the resulting nucleosynthesis and remnant. 

ABSTRACT The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever-increasing computing power makes now possible to simulate in three-dimensional (3D) part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, and abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same time-scale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur. 

ABSTRACT Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multidimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M⊙ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution; thus, revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions, and compact remnants. 

Abstract Rotating nonaxisymmetric neutron stars (NSs) are promising sources for continuous gravitational waves (CWs). CWs may, if detected, inform us about the internal structure and equation of state of NSs. Here, we present a narrowband search for CWs from known pulsars, for which a matched-filter search can be applied. Narrowband searches are robust to mismatches between electromagnetic (EM) and gravitational emissions, in contrast to fully targeted searches where they are assumed to be phase-locked. In this work, we search for the CW counterparts emitted by 34 pulsars using data from the first and second parts of the fourth LIGO–Virgo–KAGRA observing run. We use the 5n-vector narrowband pipeline, which applies frequency-domain matched filtering. In previous searches, it covered a narrow range in the frequency—frequency time derivative (f— $\dot{f}$) space. Here, we also explore a range in the second time derivative of the frequency $\ddot{f}$around the EM observations. Additionally, for the first time, we target sources in a binary system with this kind of search. We find no evidence for CWs and therefore set upper limits on the strain amplitude emitted by each pulsar. For 20 analyses, we report an upper limit below the theoretical spin-down limit. The tightest constraint is for pulsar PSR J0534+2200 (the Crab pulsar), for which our strain upper limit on the CW amplitude corresponds to ≤0.04% of the spin-down power being radiated in the CW channel. 

Abstract We detail the population properties of binary neutron star, neutron star–black hole binary, and binary black hole mergers using 158 events from the cumulative Gravitational-Wave Transient Catalog 4.0. The black hole primary mass distribution consists of a power-law-like continuum that steepens above 35M_⊙with overdensities at 10M_⊙and 35M_⊙. Binary black holes with primary masses near 10M_⊙are more likely to have less massive secondaries, with a mass ratio distribution peaking at $q=0.71_{-0.11}^{+0.12}$, potentially a signature of stable mass transfer during binary evolution. Black hole spins are inferred to be nonextremal, with 90% of black holes havingχ < 0.6, and preferentially aligned with binary orbits, implying many merging binaries form in isolation. However, we find that a significant fraction, 0.23–0.41, of binaries have negative effective inspiral spins, suggesting many could be formed dynamically in gas-free environments. We find evidence for correlation between effective inspiral spin and mass ratio, driven either by variation in the mode of the distribution or the width. The binary black hole merger rate increases with redshift, consistent with the cosmic star formation density. While there is no evidence of the mass spectrum evolving with redshift, the distribution of effective inspiral spin is found to broaden as redshift increases out toz ≈ 1. We infer the local merger rates (atz= 0) to be 8.8–250 Gpc⁻³yr⁻¹for binary neutron stars, 9.3–86 Gpc⁻³yr⁻¹for neutron star–black hole binaries, and 13–26 Gpc⁻³yr⁻¹for binary black holes; all values reflect central 90% credible intervals. 

The angular distribution of gravitational-wave power from persistent sources may exhibit anisotropies arising from the large-scale structure of the Universe. This motivates directional searches for astrophysical and cosmological gravitational-wave backgrounds, as well as continuous-wave emitters. We present results of such a search using data from the first observing run through the first portion of the fourth observing run of the LIGO-Virgo-KAGRA Collaborations. We apply gravitational-wave radiometer techniques to generate skymaps and search for both narrowband and broadband persistent gravitational-wave sources. Additionally, we use spherical harmonic decomposition to probe spatially extended sources. No evidence of persistent gravitational-wave signals is found, and we set the most stringent constraints to date on such emissions. For narrowband point sources, our sensitivity estimate to effective strain amplitude lies in the range $\left(0.03–8.4\right)\times {10}^{-24}$across all-sky and frequency range ( $20-160$) Hz. For targeted sources—Scorpius X-1, SN 1987A, the Galactic Center, Terzan 5, and NGC 6397—we constrain the strain amplitude with best limits ranging from $\sim 1.1\times {10}^{-25}$to $6.5\times {10}^{-24}$. For persistent broadband sources, we constrain the gravitational-wave flux $F_{\alpha ,\hat{n}}^{95\%,\mathrm{UL}}\left(25\mathrm{Hz}\right)<\left(0.008–5.5\right)\times {10}^{-8}\mathrm{erg}{\mathrm{cm}}^{-2}{\mathrm{s}}^{-1}{\mathrm{Hz}}^{-1}$, depending on the sky direction $\hat{n}$and spectral index $\alpha =0,2/3$, 3. Finally, for extended sources, we place upper limits on the angular power spectrum $C_{\ell }^{1/2}<\left(0.63–17\right)\times {10}^{-10}{\mathrm{sr}}^{-1}$. 

We present a search for gravitational waves from inspiraling, planetary-mass ultracompact binaries using data from the first part of the fourth observing run of LIGO, Virgo, and KAGRA. Finding no evidence of such systems, we determine the maximum distance reach for such objects and their merger rate densities. Then, we identify classes of primordial black hole mass distributions for which these rate limits can be translated into relevant constraints on the mass distribution of primordial black holes, assuming that they are all composed of dark matter, in the mass range $[{10}^{-6},{10}^{-3}]M_{\odot }$. Our constraints for the class of primordial black hole mass distribution functions $f\left(m\right)$we consider reach down to $f\left(m\right)\sim 0.1$. They are consistent with existing microlensing results in the planetary-mass range, provide a complementary probe to subsolar mass objects, and are publicly available [1]. 

Abstract Among known neutron stars, the Vela pulsar is one of the best targets for gravitational-wave searches. It is also one of the most prolific in terms of glitches, which are sudden frequency changes in a pulsar’s rotation. Such glitches could cause a variety of transient gravitational-wave signals. Here, we search for signals associated with a Vela glitch on 2024 April 29 in data of the two LIGO detectors from the fourth LIGO–Virgo–KAGRA observing run. We search both for seconds-scale burst-like emission, primarily from fundamental (f-)mode oscillations, and for longer quasi-monochromatic transients up to 4 months in duration, primarily from quasi-static quadrupolar deformations. We find no significant detection candidates, but for the first time we set direct observational upper limits on gravitational strain amplitude that are stricter than what can be indirectly inferred from the overall glitch energy scale. We discuss the short- and long-duration observational constraints in the context of specific emission models. These results demonstrate the potential of gravitational-wave probes of glitching pulsars as detector sensitivity continues to improve. 

Abstract LIGO, Virgo, and KAGRA form a network of gravitational-wave observatories. Data and analysis results from this network are made publicly available through the Gravitational Wave Open Science Center. This paper describes open data from this network, including the addition of data from the first part of the fourth observing run and selected periods from the preceding engineering run, collected from 2023 May to 2024 January. The public dataset includes calibrated strain time series for each instrument, data from additional channels used for noise subtraction and detector characterization, and analysis data products from version 4.0 of the Gravitational-Wave Transient Catalog.