We analyse a n-dimensional Generalized Uncertainty Principle (GUP) quantization framework, characterized by a non-commutative nature of the configurational variables. First, we identify a set of states which are maximally localized only along a single direction, at the expense of being less localized in all the other ones. Subsequently, in order to recover information about localization on the whole configuration space, we use the only state of the theory which exhibits maximal localization simultaneously in every direction to construct a satisfactory quasi-position representation, by virtue of a suitable translational operator. The resultant quantum framework is then applied to model the dynamics of the Bianchi I cosmology. The corresponding Wheeler-DeWitt equation is reduced to Schr\"odinger dynamics for the two anisotropy degrees of freedom, using a WKB representation for the volume-like variable of the Universe, in accordance with the Vilenkin scenario. The main result of our cosmological implementation of the constructed quantum theory demonstrates how the dynamics of a wave packet peaked at some point in the configuration space represented in the quasi-position variables, favours as the most probable configuration exactly the initial one for a relatively long time, if compared with the ordinary quantum theory. This preference arises from the distinct behavioral dynamics exhibited by wave packets in the two quantum theories.

In recent study in Ref.~\cite{Dashi} (arXiv: 2401.12525), we have introduced a method aimed at calculating the weak-field asymptotic deflection angle. This method offers an efficient computational approach that avoids the complexities of integration and cumbersome iterative procedures typically associated with deflection angle calculations. In the present paper, we expand upon this method to encompass finite-distance deflection scenarios, wherein it is postulated that both the distances from the source to the lens and from the observer to the lens are finite. Importantly, this extension naturally encompasses the case of a lens in asymptotically non-flat spacetime. As an illustrative example, we apply this method to compute the gravitational deflection angle of massive particles in Kerr spacetime while accounting for the effects of finite distance.

Intermediate mass black holes (IMBHs, $\sim 10^2-10^5M_{\odot}$) are often dubbed as the missing link between stellar mass ($\lesssim 10^2M_{\odot}$) and super-massive ($\gtrsim 10^{5-6} M_{\odot}$) black holes. Observational signatures of these can result from tidal disruption of white dwarfs (WDs), which would otherwise be captured as a whole by super-massive black holes. Recent observations indicate that IMBHs might be rapidly spinning, while it is also known that isolated white dwarfs might have large spins, with spin periods of the order of minutes. Here, we aim to understand the effects of ``coupling'' between black hole and stellar spin, focussing on the tidal disruption of spinning WDs in the background of spinning IMBHs. Using smoothed particle hydrodynamics, we perform a suite of numerical simulations of partial tidal disruptions, where spinning WDs are in eccentric orbits about spinning IMBHs. We take a hybrid approach, where we integrate the Kerr geodesic equations while being in a regime where we can treat the internal stellar fluid dynamics in the Newtonian limit. We find substantial effects of the ``coupling'' between the black hole spin and the spin of the white dwarf, although the pericenter distance of the white dwarf is taken to be large enough so that the Newtonian limit of its fluid dynamics is a robust approximation. In particular, the core mass, the bound tail mass, and the mass difference between the two tidal tails strongly depend on such ``coupled'' spin effects. However, the late time fallback rate of the debris behaves similar to the non-spinning cases. We also compute gravitational wave amplitudes and find that while the black hole spin influences the same, there is no evidence of influence of stellar spin on such amplitudes in our regime of interest.

We present a new method and implementation to obtain Bayesian posteriors on the amplitude parameters $\{h_0, \cos \iota, \psi, \phi_0\}$ of continuous-gravitational waves emitted by known pulsars. This approach leverages the well-established $\mathcal{F}$-statistic framework and software. We further explore the benefits of employing a likelihood function that is analytically marginalized over $\phi_0$, which avoids signal degeneracy problems in the $\psi$-$\phi_0$ subspace. The method is tested on simulated signals, hardware injections in Advanced-LIGO detector data, and by performing percentile-percentile (PP) self-consistency tests of the posteriors via Monte-Carlo simulations. We apply our methodology to PSR J1526-2744, a recently discovered millisecond pulsar. We find no evidence for a signal and obtain a Bayesian upper limit $h_0^{95\%}$ on the gravitational-wave amplitude of approximately $7 \times 10^{-27}$, consistent with a previous frequentist upper limit.

The recently reported signal of common red noise between pulsars by several pulsar timing array collaborations has been thought as evidence of the stochastic gravitational wave background (SGWB) due to the Helling-Downs correlation. In this study, we will search for the non-linearities of SGWB through its effect on the overlap reduction function in NANOGrav 15-year data set. In particular, we focus on a folded non-Gaussian component to SGWB whose amplitude is quantified with a single parameter $\alpha$ in the unpolarized case. The results reveal that such a non-Gaussianity of SGWB is favored, and $\alpha=0$ (the Gaussian SGWB) is excluded at about $3\sigma$ level.

We analyse the deformations of a cylindrical elastic body resulting from displacements in a varying gravitational field.

In this work we establish the analog at background level between the Created Cold Dark Matter (CCDM) model and diverse dark energy (DE) models. We find that for a barotropic or linear equation of state (EoS) for the DE, $p=w \rho$, and standard matter sector the corresponding CCDM model is described by the same functional structure of the ratio of particle production, $\Gamma$. For a different EoS the functional form of the $\Gamma$ term is not longer maintained, however, in the case of a polytropic EoS given by the Chaplygin gas the resulting $\Gamma$ term can be written as the one obtained in the barotropic case under certain considerations.

In this research contribution, we explore the relativistic quantum dynamics of spin-0 scalar particles within a curved space-time background containing topological defects, while incorporating the effects of rainbow gravity. The Klein Gordon equation is solved in a spherically symmetrical Godel type space time, considering the presence of topological defects and the influence of rainbow gravity, described by the line element:

In the framework of the extended Einstein-aether-axion theory we study the model of a two-level aetheric control over the evolution of a spatially isotropic homogeneous Universe filled with axionic dark matter. Two guiding functions are introduced, which depend on the expansion scalar of the aether flow, equal to the tripled Hubble function. The guiding function of the first type enters the aetheric effective metric, which modifies the kinetic term of the axionic system; the guiding function of the second type predetermines the structure of the potential of the axion field. We obtained new exact solutions of the total set of master equations of the model (with and without cosmological constant), and studied in detail four analytically solvable submodels, for which both guiding functions are reconstructed and illustrations of their behavior are presented.

Simulating realistic time-domain observations of gravitational waves (GWs) and GW detector glitches can help in advancing GW data analysis. Simulated data can be used in downstream tasks by augmenting datasets for signal searches, balancing data sets for machine learning, and validating detection schemes. In this work, we present Conditional Derivative GAN (cDVGAN), a novel conditional model in the Generative Adversarial Network framework for simulating multiple classes of time-domain observations that represent gravitational waves (GWs) and detector glitches. cDVGAN can also generate generalized hybrid samples that span the variation between classes through interpolation in the conditioned class vector. cDVGAN introduces an additional player into the typical 2-player adversarial game of GANs, where an auxiliary discriminator analyzes the first-order derivative time-series. Our results show that this provides synthetic data that better captures the features of the original data. cDVGAN conditions on three classes, two denoised from LIGO blip and tomte glitch events from its 3rd observing run (O3), and the third representing binary black hole (BBH) mergers. Our proposed cDVGAN outperforms 4 different baseline GAN models in replicating the features of the three classes. Specifically, our experiments show that training convolutional neural networks (CNNs) with our cDVGAN-generated data improves the detection of samples embedded in detector noise beyond the synthetic data from other state-of-the-art GAN models. Our best synthetic dataset yields as much as a 4.2% increase in area-under-the-curve (AUC) performance compared to synthetic datasets from baseline GANs. Moreover, training the CNN with hybrid samples from our cDVGAN outperforms CNNs trained only on the standard classes, when identifying real samples embedded in LIGO detector background (4% AUC improvement for cDVGAN).

Bosonic fields (within suitable mass range) may be collectively generated by rotating black holes through the black hole superradiance process. The resulting black hole is surrounded by a ``cloud" of particles whose wave function populates the superradiant energy level of the black hole. For comparable mass ratio binary black hole systems, it has been suggested that these clouds are mostly depleted at large binary separations because of a resonant level mixing effect. As a result, these clouds may not be dynamically relevant for black hole and neutron star binaries that enter the LIGO and LISA detection frequency band. In this work, we point out that the common envelope process during a compact binary evolution may bring the binary to $\sim 0.01$AU in hundreds to thousands of years, so that the resonant level mixing of hyperfine levels are no longer important. We derive a relevant regime of binary parameters where the clouds are still present for binary entering the LISA band. When the binary separation further decreases due to gravitational wave radiation, we discuss the impact of non-resonant level mixing for cloud depletion, as well as possible cloud mass transfer between the binary objects.

Neutron stars are compact and dense celestial objects that offer the unique opportunity to explore matter and its interactions under conditions that cannot be reproduced elsewhere in the Universe. Their extreme gravitational, rotational and magnetic energy reservoirs fuel the large variety of their emission, which encompasses all available multi-messenger tracers: electromagnetic and gravitational waves, neutrinos, and cosmic rays. However, accurately measuring global neutron-star properties such as mass, radius, and moment of inertia poses significant challenges. Probing internal characteristics such as the crustal composition or superfluid physics is even more complex. This article provides a comprehensive review of the different methods employed to measure neutron-star characteristics and the level of reliance on theoretical models. Understanding these measurement techniques is crucial for advancing our knowledge of neutron-star physics. We also highlight the importance of employing independent methods and adopting a multi-messenger approach to gather complementary data from various observable phenomena as exemplified by the recent breakthroughs in gravitational-wave astronomy and the landmark detection of a binary neutron-star merger. Consolidating the current state of knowledge on neutron-star measurements will enable an accurate interpretation of the current data and errors, and better planning for future observations and experiments.

There are currently two open questions in white dwarf physics: why are massive dwarfs observed less often in astronomical surveys, and why have not any super-Chandrasekhar white dwarfs been found despite the discovery of more than a dozen peculiar, overly-luminous type Ia supernovae in about a couple of decades? According to different research, magnetic fields appear to somewhat resolve these issues, but stability remains a concern. For the first time, we investigate how modified gravity affects the specific heat of electrons and ions, the crystallization process, and the cooling mechanism in white dwarfs. We demonstrate it for the Ricci-based gravity. We show that massive white dwarfs fade faster and conclude that it could be a physical reason, apart from the presence of high magnetic fields, both for finding fewer massive white dwarfs and the lack of direct detection of super-Chandrasekhar white dwarfs.

We present a theoretical investigation of the expected experimental signals produced by freely falling atoms with time oscillating mass and transition frequency. These oscillations could be produced in a variety of models, in particular, models of scalar dark matter (DM) non universally coupled to the standard matter (SM) such as axion-like particles (ALP) and dilatons. Performing complete and rigorous calculations, we show that, on one hand, two different atomic species would accelerate at a different rate, and on the other hand, they would produce a non-zero differential phase shift in atom interferometers (AI). The former would produce observable signals in equivalence principle tests like the recent MICROSCOPE mission, and we provide a corresponding sensitivity estimate, showing that MICROSCOPE can reach beyond the best existing searches in the ALP case. We also compare the expected sensitivity of two future AI experiments, namely the AION-10 gradiometer and an isotope differential AI considered for MAGIS-100, that we will refer to as SPID. We show that the SPID setup would be more sensitive to these dark matter fields compared to the gradiometer one, assuming equivalent experimental parameters.

Entanglement can be generated through the gravitational interaction between two massive bodies that are initially in a product state. This shows that the gravitational field is quantum. When the third massive body is introduced and the gravitational interaction only between the third body with either one of the former two bodies is considered, we find that no entanglement is generated between the former two bodies up to the monopole approximation, even though the considered gravitational interaction is quantum. This resembles the behavior of two accelerating two-level atoms that is usually regarded as the Unruh-DeWitt detectors. By linking the acceleration to that generated by the gravitational field, we show that the equivalence principle is still valid even though the gravitational field is quantum.

We report an approximate, non-trivial mapping of angular (linear) momentum in gravitational waves obtained from numerical relativity (NR) and adiabatic point-particle black hole perturbation theory (BHPT) in the comparable mass regime for quasi-circular, non-spinning binary black holes. This mapping involves two time-independent scaling parameters, $\alpha_{J}$ ($\alpha_{P}$) and $\beta_{J}$ ($\beta_{P}$), that adjust the BHPT angular (linear) momentum and the BHPT time respectively such a way that it aligns with NR angular (linear) momentum. Our findings indicate that this scaling mechanism works really well until close to the merger. In addition to the comparison of $\alpha_{J}$ ($\alpha_{P}$) with relevant values obtained from the waveform and flux scalings, we explore the mass ratio dependence of the scaling parameter $\alpha_{J}$ ($\alpha_{P}$). Finally, we investigate their possible connection to the missing finite size correction for the secondary black hole within the BHPT framework and the implication of these scalings on the remnant properties of the binary.

Gravitational wave physics can probe theories that extend beyond General Relativity. Motivated by recent attention on the Kalb-Ramond field as a dark matter candidate, in this work, we study a parity violating dimension four operator which couples the dual Riemann curvature to the 2-form field. After mapping the equations of motion for the right- and left-handed gravitational wave amplitudes to the novel parameterization recently presented in \cite{Jenks:2023pmk}, we discuss various constraints on the model parameters in light of the coincident electromagnetic/gravitational wave signal of GW170817 and the GWTC-3 dataset.

Evidence for a stochastic gravitational wave (GW) background, plausibly originating from the merger of supermassive black holes (SMBHs), is accumulating with observations from pulsar timing arrays. An outstanding question is how inspiraling SMBHs get past the "final parsec" of separation, where they have a tendency to stall before GW emission alone can make the binary coalesce. We argue that dynamical friction from the dark matter (DM) spike surrounding the black holes is sufficient to resolve this puzzle, if the DM has a self-interaction cross section of order $0.2-5\,$cm$^2$/g. The same effect leads to a softening of the GW spectrum at low frequencies as suggested by the current data. For collisionless cold DM, the friction deposits so much energy that the spike is disrupted and cannot bridge the final parsec, while for self-interacting DM, the isothermal core of the halo can act as a reservoir for the energy liberated from the SMBH orbits. A realistic velocity dependence, such as generated by the exchange of a massive mediator like a dark photon, is favored to give a good fit to the GW spectrum while providing a large enough core. A similar velocity dependence has been advocated for solving the small-scale structure problems of cold DM.

We construct and parametrize solutions to the constraint equations of general relativity in a neighborhood of Minkowski spacetime with arbitrary prescribed decay properties at infinity. We thus provide a large class of initial data for the results on stability of Minkowski which include a mass term in the asymptotics. Due to the symmetries of Minkowski, a naive linear perturbation fails. Our construction is based on a simplified conformal method, a reduction to transverse traceless perturbations and a nonlinear fixed point argument where we face linear obstructions coming from the cokernels of both the linearized constraint operator and the Laplace operator. To tackle these obstructions, we introduce a well-chosen truncated black hole around which to perturb. The control of the parameters of the truncated black hole is the most technical part of the proof, since its center of mass and angular momentum could be arbitrarily large.

The stochastic gravitational wave background for pulsar timing arrays is often modeled by a Gaussian ensemble which is isotropic and unpolarized. However, the Universe has a discrete set of polarized gravitational wave sources at specific sky locations. Can we trust that the Gaussian ensemble is an accurate description? To investigate this, we explicitly construct an ensemble containing $N$ individual binary sources with circular orbits. The orbital inclination angles are randomly distributed, hence the individual sources are elliptically polarized. We then compute the first two moments of the Hellings and Downs correlation, as well as the pulsar-averaged correlation mean and (cosmic) variance. The first moments are the same as for a previously studied ensemble of circularly polarized sources. However, the second moments, and hence the variances, are different for the two ensembles. While neither discrete source model is exactly described by a Gaussian ensemble, we show that in the limit of large $N$, the differences are small.

A method is introduced which, for the first time, allows us to construct axial vector fields without which formal definitions of quasi-local angular momentum, in general, would remain empty. The introduced method is practical, it can be used to construct all such axial vector fields, and it allows the quasi-local angular momentum to be represented by a triple vector in three-dimensional Euclidean space. We also derive balance relations which allow us to monitor the variation of the magnitude and direction of this vector, and also to monitor the angular momentum transports in generic spacetimes without symmetries.

This is the second in a series of two articles introducing the Gravitas computational general relativity framework, in which we now focus upon the design and capabilities of Gravitas's numerical subsystem, including its ability to perform general 3+1 decompositions of spacetime via the ADM formalism, its support for the definition and construction of arbitrary Cauchy surfaces as initial data, its support for the definition and enforcement of arbitrary gauge and coordinate conditions, its various algorithms for ensuring the satisfaction of the ADM Hamiltonian and momentum constraints, and its unique adaptive refinement algorithms based on hypergraph rewriting via Wolfram model evolution. Particular attention is paid to the seamless integration between Gravitas's symbolic and numerical subsystems, its ability to configure, run, analyze and visualize complex numerical relativity simulations and their outputs within a single notebook environment, and its capabilities for handling generic curvilinear coordinate systems and spacetimes with general (and often highly non-trivial) topologies using its specialized and highly efficient hypergraph-based numerical algorithms. We also provide illustrations of Gravitas's functionality for the visualization of hypergraph geometries and spacetime embedding diagrams, the ability for Gravitas's symbolic and numerical subsystems to be used in concert for the extraction of gravitational wave signals and other crucial simulation data, and Gravitas's in-built library of standard initial data, matter distributions and gauge conditions. We conclude by demonstrating how the numerical subsystem can be used to set up, run, visualize and analyze a standard yet nevertheless reasonably challenging numerical relativity test case: a binary black hole collision and merger within a vacuum spacetime (including the extraction of its outgoing gravitational wave profile).

Gravastars are considered as one of the prime exotic type compact objects which may be found at the end state of gravitational collapse of massive stars with a view to resolve the complexities that are pertinent in case of a black hole \cite{Mazur}-\cite{Mazur2}. In this paper, we analyse the role of charge on the possible formation of isotropic spherically symmetric gravastar configuration in the framework of Rastall gravity. Gravastar contains three distinct layers {\it viz.} i) Interior region, ii) Thin shell and iii) Exterior region. The interior region is characterised by the equation of state $p=-\rho$ that defines the repulsive outward pressure in radial direction at all points on the thin shell. The thin shell, contains ultra-relativistic stiff fluid which is denoted by the equation of state $p=\rho$ following Zel'dovich's criteria \cite{Zeldovich,Zeldovich1} for cold baryonic universe, can withstand the repulsive pressure exerted by the interior region. The exterior region is the vacuum space-time represented by the Reissner-Nordstr$\ddot{o}$m solution. In view of the above specifications, we construct and analyse a charged gravastar model in Rastall theory of gravity which represents several salient features. The basic physical attributes, {\it viz.} proper length, energy, entropy and equation of state parameter of the shell are investigated. In this model, it is interesting to note that for large value of the radius of hyper-surface (R) the EoS parameter of the thin shell corresponds to dark energy EoS with $\mathcal{W}(R)\rightarrow-1$. However, for small value of $R$ the EoS parameter $\mathcal{W}(R)\rightarrow0$, defines a dust shell. The stability of the model is ensured through the study of gravitational surface redshift and maximisation of shell entropy within the framework of Rastall theory of gravity.

We investigate the turnaround radius of the Reissner-Nordstr\"{o}m deSitter Spacetime and how the turnaround radius changes if a test particle carries charge. We also consider the Mart\'{i}nez-Troncoso-Zanelli (MTZ) solution of conformally coupled gravity and investigate how the turnaround radius changes for a scalar test charge. In both scalar and electric interaction cases we find that the Turnaround Radius depends on the particle's energy.

Time-domain (TD) Bayesian inference is important in ringdown analysis for gravitational wave (GW) astronomy. The validity of this method has been well studied by Isi and Farr [1]. Using GW190521 as an example, we study the TD method in detail by comparing it with the frequency-domain (FD) method as a complement to previous study. We argue that the autocovariance function (ACF) should be calculated from the inverse fast Fourier transform of the power spectral density (PSD), which is usually estimated by the Welch method. In addition, the total duration of the GW data that are used to estimate the PSD and the slice duration of the truncated ACF should be long enough. Only when these conditions are fully satisfied can the TD method be considered sufficiently equivalent to the FD method.

This study deals with the formation and evolution of a strange star in the Krori Barua Rainbow spacetime from collapsing charged stringy quark matter. The dynamical variables are explored from the field equations, taking into account the effects of particle's energy on the mass density, pressure, and string tension. The electric field is also computed using the MIT Bag model. The real time data of SAX J1808.4-3658 is used to analyzed the physical properties including gradients, energy conditions, anisotropy, stability, Tolman Oppenheimer Volkoff equation, mass function, compactness, and red-shift. The graphical analysis has been made according to both the theories of rainbow gravity and general relativity. The energy conditions and anisotropy are found to be satisfied, indicating the physical existence of suggested model. Tolman Oppenheimer Volkoff equation is satisfied indicating equilibrium of forces and stability of the compact object. Overall, our model is consistent with the observational information of SAX J1808.4-3658.

We provide spinor-based proof of the positivity of the Wang-Yau \cite{yau} quasi-local mass. More precisely we prove that the gravitational mass bounded by a spacelike topological $2-$sphere is non-negative in a generic spacetime verifying dominant energy condition and vanishes only if the surface is embedded in the Minkowski space. This proof is purely quasi-local in nature and in particular does not rely on Bartanik's gluing and asymptotic extension construction \cite{bartnik1993quasi} and subsequent application of the Schoen-Yau \cite{schoen1979proof,schoen1981proof} positive mass theorem.

Spin is one fundamental property of microscopic particles. A lot of theoretical work postulated the possible coupling between spin and gravitation, which could result in the violation of equivalence principle. In a recent joint mass-and-energy test of weak equivalence principle (WEP) with a 10-meter $^{85}$Rb-$^{87}$Rb dual-species atom interferometer, the E${\rm \ddot{o}}$tv${\rm \ddot{o}}$s parameters of four $^{85}$Rb-$^{87}$Rb combinations with specific atomic spin states were measured to the $10^{-10}$-level. These experimental results are used to constrain the postulated spin-gravity coupling effects in a robust way. The bounds on the spin-independent and spin-dependent anomalous passive gravitational mass tensors are consistently set to the $10^{-10}$-level, which improve existing bounds by three orders of magnitude.

The current standard model of cosmology - the {\ensuremath{\Lambda}}CDM model - is appropriately named after its controversial foreign ingredients: a cosmological constant ({\ensuremath{\Lambda}}) that accounts for the recent accelerated expansion of the Universe and cold dark matter needed to explain the formation and dynamics of large scale structures. Together, these form the dark sector, whose nature remains a mystery. After 25 years of withstanding confirmation and support for the {\ensuremath{\Lambda}}CDM model, enough to bypass some of its unclear theoretical issues, this paradigm is facing its biggest crisis yet. The rapid advent of technology has brought cosmology to an unprecedented observational era, with increased technical precision and the emergence of independent measures, including probes of phenomena that were thought impossible to detect or even exist, such as the gravitational ripples that propagate in the spacetime. However, such precision has unveiled cracks in the porcelain of {\ensuremath{\Lambda}}CDM, with pieces that seem glued together and difficult to reconcile. Particularly worrying is the apparent lack of compatibility between measurements of the Universe's present expansion rate based on local measurements and those based on phenomena that occurred far in the early Universe and that can only be translated into present quantities through physical propagation under a cosmological model. In this dissertation, we delve into extensions to the standard model that consider alternatives to the mysterious nature of the dark sector and any possible new interactions therein. We analyse these alternative models, hoping to identify measurable observational signatures of extra degrees of freedom in the dark sector.

Intrinsic spin of fermions can generate torsion in spacetime. This torsion is a non-propagating field that can be integrated out, leaving an effective non-universal four-fermion interaction. This geometrical interaction affects fermions inside a matter distribution and can be expected to become stronger as the density grows. In this work, we investigate the role of this interaction in a gravitationally collapsing fermionic distribution. Our specific aim is to explore if this interaction can provide a repulsive contribution and prevent the final singularity formation. We consider a collapsing distribution of massive fermions, ignoring other interactions. Using simplified yet reasonable assumptions, we establish that a repulsive contribution can arise depending on how torsion couples with different chiralities. Also, the interaction starts to dominate as the collapse proceeds, accelerating or decelerating the collapse depending on the relative signs of the geometrical interaction between different species of fermions.

The search for supersymmetric partners at Large Hadron Collider revealed negative result. Though, strictly speaking, it does not exclude low energy supersymmetry, but still it leads to strong constraints of the parameter space. Therefore the search for supersymmetric particles at higher energies becomes of interest. It is shown that in $R^2$-modified cosmology heavy particles with the interaction strength typical for supersymmetry could be promising candidates for carriers of dark matter. We consider the heating of the Universe at the post-inflationary stage via particle production by oscillating curvature scalar (scalaron). The bounds on the masses of dark matter particles are obtained for different dominant decay modes of the scalaron. Possible impact of superheavy particle decays on the spectrum of ultra high energy cosmic rays is discussed.

Gravitational-wave astronomy has developed enormously over the last decade with the first detections across different frequency bands, but has yet to access $0.1-10$ $\mathrm{Hz}$ gravitational waves. Gravitational waves in this band are emitted by some of the most enigmatic sources, including intermediate-mass binary black hole mergers, early inspiralling compact binaries, and possibly cosmic inflation. To tap this exciting band, we propose the construction of a detector based on pulsar timing principles, the Artificial Precision Timing Array (APTA). We envision APTA as a solar system array of artificial "pulsars"$-$precision-clock-carrying satellites that emit pulsing electromagnetic signals towards Earth or other centrum. In this fundamental study, we estimate the clock precision needed for APTA to successfully detect gravitational waves. Our results suggest that a clock relative uncertainty of $10^{-17}$, which is currently attainable, would be sufficient for APTA to surpass LISA's sensitivity in the decihertz band and observe $10^3-10^4$ $\mathrm{M}_\odot$ black hole mergers. Future atomic clock technology realistically expected in the next decade would enable the detection of an increasingly diverse set of astrophysical sources, including stellar-mass compact binaries that merge in the LIGO-Virgo-KAGRA band, extreme-mass-ratio inspirals, and Type Ia supernovae. This work opens up a new area of research into designing and constructing artificial gravitational-wave detectors relying on the successful principles of pulsar timing.

Hybrid neutron stars, the compact objects consisting hadronic matter and strange quark matter, can be considered as the probes for the scalar tensor gravity. In this work, we explore the scalarization of hybrid neutron stars in the scalar tensor gravity. For the hadronic phase, we apply a piecewise polytropic equation of state constrained by the observational data of GW170817 and the data of six low-mass X-ray binaries with thermonuclear burst or the symmetry energy of the nuclear interaction. In addition, to describe the strange quark matter inside the hybrid neutron star, different MIT bag models are employed. We study the effects of the value of bag constant, the mass of s quark, the perturbative quantum chromodynamics correction parameter, and the density jump at the surface of quark-hadronic phase transition on the scalarization of hybrid neutron stars. Our results confirm that the scalarization is more sensitive to the value of bag constant, the mass of s quark, and the density jump compared to the perturbative quantum chromodynamics correction parameter.

Modelling the gravitational wave signal from binaries beyond comparable mass is an important open issue in gravitational wave astronomy. For non-spinning binaries and when the spins are aligned with the orbital angular momentum, some first studies concerning the transition between the comparable and extreme mass ratio regime are already available, which suggest that extreme mass ratio results at times extrapolate to comparable mass ratios with surprising precision. Here we study the case of misaligned spins: We present new NR simulations performed with the Einstein Toolkit code at mass ratios up to 18 and construct a heterogeneous dataset that spans all mass ratios, including data from NR simulations, numerical approximations to extreme mass ratio binaries, and data from the geodesic approximation. As a first application we provide fits for the remnant mass and spin magnitude in single spin precessing systems, omitting consideration of the in-plane spin orientation. These fits demonstrate accuracy comparable to the state-of-the-art NRSur7dq4EmriRemnant model, all while retaining the simplicity and efficiency inherent in previous phenomenological fits.

Our representation of the Universe is built with sequences of symbols, numbers, operators, rules and undecidable propositions defining our mathematical truths, represented either by classical, quantum and probabilistic Turing Machines containing intrinsic randomness. Each representation is at all effects a physical subset of the Universe, a metastructure of events in space and time, which actively participate to the evolution of the Universe as we are internal observers. The evolution is a deterministic sequence of local events, quantum measurements, originated from the local wavefunction collapse of the complementary set of the observers that generate the local events in the Universe. With these assumptions, the Universe and its evolution are described in terms of a semantically closed structure without a global object-environment loss of decoherence as a von Neumann's universal constructor with a semantical abstract whose structure cannot be decided deterministically a-priori from an internal observer. In a semantically closed structure the realization of a specific event writing the semantical abstract of the constructor is a problem that finds a "which way" for the evolution of the Universe in terms of a choice of the constructor's state in a metastructure, the many-world Everett scenario from the specific result of a quantum measurement, a classical G\"odel undecidable proposition for an internal observer, exposing the limits of our description and possible simulation of the Universe.

We prove under certain weak assumptions a black hole no-hair theorem in spherically symmetric spacetimes for self-gravitating time-dependent multiple scalar fields with an arbitrary target space admitting a Killing field with a non-empty axis and arbitrary non-negative potential invariant under the flow of the Killing field. It is shown that for such configurations the only spherically symmetric and asymptotically flat black hole solutions consist of the Schwarzschild metric and a constant multi-scalar map. In due course of the proof we also unveil the intrinsic connection of the time-dependence of the scalar fields with the symmetries of the target space.

We continue ongoing research work on applying the homological algebraic conceptual and technical machinery of Abstract Differential Geometry towards formulating a finitary, causal and quantal version of vacuum Einstein Lorentzian gravity and free Yang-Mills theories.

We present the early release of the atlas of continuous gravitational waves covering frequencies from 20 Hz to 1500 Hz and spindowns from -5e-10 to 5e-10 Hz/s. Compared to the previous atlas release we have greatly expanded the parameter space, and we now also provide polarization-specific data - both for signal-to-noise ratios and for the upper limits. Continuous wave searches are computationally difficult and take a long time to complete. The atlas enables new searches to be performed using modest computing power. To allow new searches to start sooner, we are releasing this data early, before our followup stages have completed.

Stellar evolution predicts the existence of a mass gap for black hole remnants produced by pair-instability supernova dynamics, whose lower and upper edges are very uncertain. We study the possibility of constraining the location of the upper end of the pair-instability mass gap, which is believed to appear around ${m_\text{min}} \sim130M_\odot$, using gravitational wave observations of compact binary mergers with next-generation ground-based detectors. While high metallicity may not allow for the formation of first-generation black holes on the "far side" beyond the gap, metal-poor environments containing Population III stars could lead to such heavy black hole mergers. We show that, even in the presence of contamination from other merger channels, next-generation detectors will measure the location of the upper end of the mass gap with a relative precision close to $\Delta {m_\text{min}}/{m_\text{min}} \simeq 4\% (N_\text{det}/100 )^{-1/2}$ at 90% C.L., where $N_\text{det} $ is the number of detected mergers with both members beyond the gap. These future observations could reduce current uncertainties in nuclear and astrophysical processes controlling the location of the gap.

We present the first application of the weighted skew-spectra to analyze non-Gaussian information in galaxy survey data. Using the tree-level galaxy skew-spectra together with the one-loop power spectrum multipoles, we analyze the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) galaxy clustering data, and target our search towards the equilateral bispectrum shape of primordial non-Gaussianity. We use the Effective Field Theory model for the galaxy power spectrum and bispectrum, and account for systematic effects, such as the survey geometry. From our likelihood analysis, we find $f_{\rm NL}^{\rm equil}=-34^{+296}_{-334}$ at $68\%$ CL, consistent with previous works, while systematic errors from our treatment of the survey geometry lead to an unreliable estimation of $f_{\rm NL}^{\rm ortho}$. We further constrain the bias and counterterm parameters, while keeping the cosmology fixed to $\textit{Planck }2018$ values. As a check, we also validate our analysis pipeline using the ${\tt Nseries}$ simulation suite.

In this review, we aim to analyze the mathematical interpretation of the Ashtekar-Barbero-Immirzi formulation of General Relativity. Along with a brief introduction to the necessary mathematical structures and tools, we illustrate some relevant physical theory quantities as geometrical objects within the framework of principal bundle theory.

We consider conservation of momentum in AQUAL, a field-theoretic extension to Modified Newtonian Dynamics (MOND). We show that while there is a sense in which momentum is conserved, it is only if momentum is attributed to the gravitational field, and thus Newton's third law fails as usually understood. We contrast this situation with that of Newtonian gravitation on a field theoretic formulation. We then briefly discuss the situation in TeVeS, a relativistic theory that has AQUAL as a classical limit.

We address in this article the criticism in a recently submitted article by Clement and Noiucer (arXiv:2312.17662 [gr-qc]) on "Cotton Gravity" (CG), a gravity theory alternative to General Relativity. These authors claim that CG is "not predictive" for producing "too many" spherically symmetric vacuum solutions, while taking the Bianchi I vacuum as test case they argue that geometric constraint on the Cotton tensor lead to an undetermined problem, concluding in the end that CG "is not a physical theory". We provide arguments showing that this critique is incorrect and misrepresents the theory.

This paper presents the study of dark-energy compact stars in the context of modified Rastall teleparallel gravity. It is the first time that dark energy celestial phenomena have been explored in this modified gravitational theory. Employing the torsion-based functions, $f(T)$ and $h(T)$, we analyzed their effects in a spherically symmetric spacetime chosen as the interior geometry, while using the Schwarzschild geometry as an outer spacetime. In this study, we explored various dark energy stellar properties, including dark energy pressure components, energy conditions, and equation of state components. Our findings reveal that the observed negative behavior of these stellar properties served as compelling evidence, validating the presence of dark energy in stellar configurations. Detailed investigations of the energy conditions, pressure profiles, sound speeds, TOV equation, adiabatic index, gradients, mass function, compactness, and redshift function forecasts a comprehensive assessment, affirming the acceptability and realism of the investigated stellar configuration.

The equations of state (EoSs) governing neutron star (NS) matter obtained for both non-relativistic and relativistic mean-field models are systematically confronted with a diverse set of terrestrial data and astrophysical observations within the Bayesian framework. The terrestrial data, spans from bulk properties of finite nuclei to the heavy-ion collisions, constrain the symmetric nuclear matter EoS and the symmetry energy up to twice the saturation density ($\rho_0$= 0.16 fm$^{-3}$). The astrophysical observations encompass the NS radius, the tidal deformability, and the lower bound on maximum mass. Three distinct posterior distributions of EoSs are generated by gradually updating the priors with different constraints: (i) only the maximum NS mass, (ii) incorporating additional terrestrial data, (iii) combining both the terrestrial data and astrophysical observations. These distributions are then compared using the Kullback-Liebler divergence which highlights the significant constraints imposed on the EoSs by the currently available lower bound of NS maximum mass and terrestrial data. The remaining astrophysical observations marginally refine the EoS within the density range $\sim$ 2-3$\rho_0$. It is observed that the non-relativistic mean field model yields stiffer EoS for the symmetric nuclear matter while softer density dependent symmetry energy compared to the relativistic one when all the data/observations are considered.

Starting from parametrization of scalar perturbations generated during inflation in terms of $e$--fold $N$ and using an approach recently developed by Starobinsky, dubbed "direct smooth reconstruction", we show that, in the slow--roll approximation, it is possible to reconstruct the inflation potential in the Einstein frame and its corresponding Lagrangian $f(R)$ model

Aims: Ab initio global particle-in-cell (PIC) simulations of compact neutron star magnetospheres in the align rotator configuration to investigate the role of GR and plasma supply on the polar cap particle acceleration efficiency - precursor of coherent radio emission.

Methods: A new module for the PIC code OSIRIS to model plasma dynamics around compact objects with full self-consistent GR effects is presented. A detailed description of the extensions and implementation methods are provided for the main sub-algorithms of the PIC loop, including the field solver, particle pusher and charge-conserving current-deposit scheme.

Results: Leptons are efficiently accelerated in the polar caps of neutron stars with force-free magnetospheres. This solution supports strong poloidal currents, which are easily turned to spacelike at any stellar compactness due to the GR frame-dragging effect. Charge-separated magnetospheric solutions, in opposition, depend on the plasma supply and compactness to activate the polar cap. The furthest away the solution is from force-free, the highest the compactness is required for the polar cap activation.

Conclusions: GR effects are crucial to explain the pulsar mechanism for low obliquity rotators. Focusing on the aligned rotator, we show that GR relaxes the minimum required poloidal magnetospheric current for the transition of the polar cap to the accelerator regime, thus justifying the observation of weak pulsars beyond the expected death line. Also, we demonstrate how the interplay between the polar cap and outer gaps might explain the intermittent behaviour of the measured spin-down luminosity and the existence of radio sub-pulse nullings for older pulsars.

The symmetry frame formalism is an effective tool for computing the symmetries of a Riemann-Cartan geometry and, in particular, in metric teleparallel geometries. In the case of non-vanishing torsion in a four dimensional RC geometry, the Minkowski geometry is the only geometry admitting ten affine frame symmetries. Excluding this almost trivial geometry, the maximal number of affine frame symmetries is seven. A natural question is to ask what four dimensional geometries admit a seven-dimensional group of affine frame symmetries. Such geometries are locally homogeneous and admit the largest isotropy group permitted, and hence are called {\it maximally isotropic}. Using the symmetry frame formalism to compute affine frame symmetries along with the additional structure of the torsion tensor, we employ the Cartan-Karlhede algorithm to determine all possible seven-dimensional symmetry groups for Riemann-Cartan geometries.

We present for the first time a study of the quasinormal modes of rapidly rotating Ellis-Bronnikov wormholes in General Relativity. We compute the spectrum of the wormholes using a spectral decomposition of the metric perturbations on a numerical background. We focus on the $M_z=2,3$ sector of the perturbations, and show that the triple isospectrality of the symmetric and static Ellis-Bronnikov wormhole is broken due to rotation, giving rise to a much richer spectrum than the spectrum of Kerr black holes. We do not find any instabilities for $M_z=2,3$ perturbations.

We quantize a closed homogeneous and isotropic universe in the modified teleparallel gravity, wherein a scalar field is non-minimally coupled to both the torsion and a boundary term. In this regard, we study exact solutions at both the classical and quantum frameworks by utilizing the corresponding Wheeler-DeWitt (WDW) equation of the model. To correspond to the comprehensive classical and quantum levels, we propose an appropriate initial condition for the wave packets and observe that they closely adhere to the classical trajectories and reach their peak. We quantify this correspondence using the de-Broglie Bohm interpretation of quantum mechanics. According to this proposal, the classical and Bohmian trajectories coincide when the quantum potential vanishes along the Bohmian paths. Furthermore, we apply the de-parameterization technique to our model in the realm of the problem of time in quantum cosmological models based on the WDW equation, utilizing the global internal time denoted as $\chi$, which represents a scalar field.

The Milky Way galaxy is estimated to be home to ten million to a billion stellar-mass black holes (BHs). Accurately determining this number and distribution of BH masses can provide crucial information about the processes involved in BH formation, the possibility of the existence of primordial BHs, and interpreting gravitational wave (GW) signals detected in LIGO-VIRGO-KAGRA. Sahu et al. recently confirmed one isolated stellar-mass BH in our galaxy using astrometric microlensing. This work proposes a novel method to identify such BHs using the gravitational analog of the Gertsenshtein-Zel'dovich (GZ) effect. We explicitly demonstrate the generation of GWs when a kilohertz(kHz) electromagnetic (EM) pulse from a pulsar is intervened by a spherically symmetric compact object situated between the pulsar and Earth. Specifically, we show that the curvature of spacetime acts as the catalyst, akin to the magnetic field in the GZ effect. Using the covariant semi-tetrad formalism, we quantify the GW generated from the EM pulse through the Regge-Wheeler tensor and express the amplitude of the generated GW in terms of the EM energy and flux. We demonstrate how GW detectors can detect stellar-mass BHs by considering known pulsars within our galaxy. This approach has a distinct advantage in detecting stellar mass BHs at larger distances since the GW amplitude falls as $1/r$.

It is currently unknown how matter behaves at the extreme densities found within the cores of neutron stars. Measurements of the neutron star equation of state probe nuclear physics that is otherwise inaccessible in a laboratory setting. Gravitational waves from binary neutron star mergers encode details about this physics, allowing the equation of state to be inferred. Planned third-generation gravitational-wave observatories, having vastly improved sensitivity, are expected to provide tight constraints on the neutron star equation of state. We combine simulated observations of binary neutron star mergers by the third-generation observatories Cosmic Explorer and Einstein Telescope to determine future constraints on the equation of state across a plausible neutron star mass range. In one year of operation, a network consisting of one Cosmic Explorer and the Einstein Telescope is expected to detect $\gtrsim 3\times 10^5$ binary neutron star mergers. By considering only the 75 loudest events, we show that such a network will be able to constrain the neutron star radius to at least $\lesssim 200$ m (90% credibility) in the mass range $1-1.97$ $M_{\odot}$ -- about ten times better than current constraints from LIGO-Virgo-KAGRA and NICER. The constraint is $\lesssim 75$ m (90% credibility) near $1.4-1.6$ $M_{\odot}$ where we assume the the binary neutron star mass distribution is peaked. This constraint is driven primarily from the loudest $\sim 20$ events.

We present an overview of recent numerical advances in the theoretical characterization of massive binary black hole (MBBH) mergers in astrophysical environments. These systems are among the loudest sources of gravitational waves (GWs) in the universe and particularly promising candidates for multimessenger astronomy. Coincident detection of GWs and electromagnetic (EM) signals from merging MBBHs is at the frontier of contemporary astrophysics. One major challenge in observational efforts searching for these systems is the scarcity of strong predictions for EM signals arising before, during, and after merger. Therefore, a great effort in theoretical work to-date has been to characterize EM counterparts emerging from MBBHs concurrently to the GW signal, aiming to determine distinctive observational features that will guide and assist EM observations. To produce sharp EM predictions of MBBH mergers it is key to model the binary inspiral down to coalescence in a full general relativistic fashion by solving Einstein's field equations coupled with the magnetohydrodynamics equations that govern the evolution of the accreting plasma in strong-gravity. We review the general relativistic numerical investigations that have explored the astrophysical manifestations of MBBH mergers in different environments and focused on predicting potentially observable smoking-gun EM signatures that accompany the gravitational signal.

We explore the prospects for identifying differences in simulated gravitational-wave signals of binary neutron star (BNS) mergers associated with the way thermal effects are incorporated in the numerical-relativity modelling. We consider a hybrid approach in which the equation of state (EoS) comprises a cold, zero temperature, piecewise-polytropic part and a thermal part described by an ideal gas, and a tabulated approach based on self-consistent, microphysical, finite-temperature EoS. We use time-domain waveforms corresponding to BNS merger simulations with four different EoS. Those are injected into Gaussian noise given by the sensitivity of the third-generation detector Einstein Telescope and reconstructed using BayesWave, a Bayesian data-analysis algorithm that recovers the signals through a model-agnostic approach. The two representations of thermal effects result in frequency shifts of the dominant peaks in the spectra of the post-merger signals, for both the quadrupole fundamental mode and the late-time inertial modes. For some of the EoS investigated those differences are large enough to be told apart, especially in the early post-merger phase when the signal amplitude is the loudest. These frequency shifts may result in differences in the inferred tidal deformability, which might be resolved by third-generation detectors up to distances of about tens of Mpc at most.

In this paper, we construct the slowly rotating case of an asymptotically flat supermassive black hole embedded in dark matter using Newman-Janis procedure. Our analysis is carried with respect to the involved parameters including the halo total mass $M$ and the galaxy's lengthscale $a_0$. Concretly, we investigate the dark matter impact on the effective potential and the photon sphere. In particular, we find that the lengthscale $a_0$ controles such potential values. Indeed, for low $a_0$ values, we find that the halo total mass $M$ decreases the potential values significantly while for high $a_0$ values such impact is diluted. Regarding the shadow aspects, we show that the shadow size is much smaller for high values of $a_0$ while the opposite effect is observed when the halo total mass $M$ is increased. By comparing our case to the slowly rotating case, we notice that the former exhibits a shadow shifted from its center to the left side. Finally, we compute the deflection angle in the weak-limit approximation and inspect the dark matter parameters influence. By ploting such quantity, we observe that one should expect lower bending angle values for black holes in galactic nuclei.

The search for subsolar mass primordial black holes (PBHs) poses a challenging problem due to the low signal-to-noise ratio, extended signal duration, and computational cost demands, compared to solar mass binary black hole events. In this paper, we explore the possibility of investigating the mass range between subsolar and planetary masses, which is not accessible using standard matched filtering and continuous wave searches. We propose a systematic approach employing the Viterbi algorithm, a dynamic programming algorithm that identifies the most likely sequence of hidden Markov states given a sequence of observations, to detect signals from small mass PBH binaries. We formulate the methodology, provide the optimal length for short-time Fourier transforms, and estimate sensitivity. Subsequently, we demonstrate the effectiveness of the Viterbi algorithm in identifying signals within mock data containing Gaussian noise. Our approach offers the primary advantage of being agnostic and computationally efficient.

In this research, we delve into the localization patterns of fermionic fields within a braneworld setting, employing a modified gravity model denoted as $f(Q)$. Our investigation revolves around two specific models, $f_1(Q)=Q+kQ^n$ and $f_2(Q)=Q+k_1Q^2+k_2Q^3$, where we systematically vary the parameters $n$ and $k_{1,2}$. Through an in-depth analysis encompassing the effective potential, massless, and massive modes, we elucidate how deviations from the conventional symmetric teleparallel equivalent of general relativity (STEGR) gravity impact the localization of fermionic fields. To ensure greater precision, our methodology integrates probabilistic measures such as Shannon entropy and relative probability. Moreover, we gauge the stability of these models employing differential configurational entropy (DCE), revealing a compelling correlation between the most stable configurations and the emergence of novel structures within the background scalar field. This work significantly contributes to our understanding of the gravitational modifications' intricate influence on fermionic field localization within braneworld scenarios. By shedding light on these dynamics, it advances the broader comprehension of the interplay between gravity modifications and fermionic field behaviors in these theoretical frameworks.

This study revisits supersymmetric (SUSY) hybrid inflation in light of CMB experiments and swampland conjectures. We first show that if one adds radiative, soft mass, and SUGRA corrections to the scalar potential, supersymmetric hybrid inflation is still consistent with Planck 2018 despite an impression that it does not. Usually, in SUSY hybrid inflation with minimal K\"ahler potential, the gauge symmetry breaking scale $M$ turns out to be ${\cal O}(10^{15})$ GeV, which causes proton decay rate problem. In this study, we present a new parameter space where the proton decay rate problem can be avoided by achieving $M$ of the order of $10^{16}$ GeV with $M_{S}^{2}<$0 and $am_{3/2}>$0. In this scenario, one requires a soft SUSY breaking scale $|M_{S}| \gtrsim 10^{6}$ GeV. Moreover, the tensor to scalar ratio $r$ is in the range $10^{-16}$ to $10^{-6}$, which is quite small. In this case, modified swampland hold, but it is difficult to satisfy trans-Planckian censorship conjecture. For this reason, we also consider non-minimal K\"ahler potential. We fixed spectral index $n_{S}=$0.9665 (central value) of Planck 2018 data and $M=2\times 10^{16}$ GeV and present our calculations. We show the canonical measure of primordial gravity waves $r$ for $M_{S}=$ 1 TeV, $m_{3/2}=$ 1 TeV, $\kappa_{S}<0$ for $\cal{N}=$1 and $\cal{N}=$2, ranges from $10^{-5}$ to $0.01$ which can be observed in Planck and next-generation experiments such as LiteBIRD, Simons Observatory, PRISM, PIXIE,CORE, CMB-S4 and CMB-HD experiments that are gearing up to measure it. In addition to it, we present the parametric space and benchmark points for a non-minimal case which is consistent with modified swampland and trans-Planckian censorship conjectures.

We study the evolution of scalar and tensor cosmological perturbations in the framework of the Einstein-Cartan theory of gravity. The value of the gravitational slip parameter which is defined as the ratio of the two scalar potentials in the Newtonian gauge, can be used to determine whether or not the gravity is modified. We calculate the value of slip parameter in the Einstein-Cartan cosmology and show that it falls within the observed range. We also discuss the evolution of the cosmic gravitational waves as another measure of the modification of gravity.

In this manuscript, we show that three fundamental building blocks are supporting the Cosmological Principle. The first of them states that there is a special frame in the universe where the spatial geometry is intrinsically homogeneous and isotropic. The second demands the existence of a fiducial observer to whom the Hubble parameter is isotropic. The last piece states that matter and radiation behave as a perfect fluid. We show that these three hypotheses give us the Friedmann-Lema\^itre-Robertson-Walker (FLRW) spacetimes, the central pillar of the standard model of Cosmology. We keep with the first of them and start to investigate the so-called intrinsically homogeneous and isotropic spacetimes. They emerge after the decoupling of the CMB with the geometric frame of reference. Furthermore, a ``$\Lambda$CDM-like" effective theory arises naturally in those backgrounds, together with some new density parameters relating to the local inhomogeneities, the internal energy density, and the local and global magnitudes of the Hubble anisotropy. All those properties make this class of inhomogeneous models, which roughly speaking, keeps "1/3" of the Cosmological Principle, worth investigating in applications to Cosmology, for it can accommodate the same ingredients of the standard model, as a geometric frame and a free-falling isotropic cosmic background radiation, and reduce to the latter when some observable parameters vanish.

As a potential candidate for the late-time accelerating expansion of the Universe, the Chaplygin gas and its generalized models have significant implications to modern cosmology. In this work we investigate the effects of dark energy on the internal structure of a neutron star composed of two phases, which leads us to wonder: Do stable neutron stars have a dark-energy core? To address this question, we focus on the radial stability of stellar configurations composed by a dark-energy core -- described by a Chaplygin-type equation of state (EoS) -- and an ordinary-matter external layer which is described by a polytropic EoS. We examine the impact of the rate of energy densities at the phase-splitting surface, defined as $\alpha= \rho_{\rm dis}^-/\rho_{\rm dis}^+$, on the radius, total gravitational mass and oscillation spectrum. The resulting mass-radius diagrams are notably different from dark energy stars without a common-matter crust. Specifically, it is found that both the mass and the radius of the maximum-mass configuration decrease as $\alpha$ becomes smaller. Furthermore, our theoretical predictions for mass-radius relations consistently describe the observational measurements of different massive millisecond pulsars as well as the central compact object within the supernova remnant HESS J1731-347. The analysis of the normal oscillation modes reveals that there are two regions of instability on the $M(\rho_c)$ curve when $\alpha$ is small enough indicating that the usual stability criterion $dM/d\rho_c>0$ still holds for rapid phase transitions. However, this is no longer true for the case of slow transitions.

Mergers of binary compact objects, accompanied with electromagnetic (EM) counterparts, offer excellent opportunities to explore varied cosmological models, since gravitational waves (GWs) and EM counterparts always carry the information of luminosity distance and redshift, respectively. $f(T)$ gravity, which alters the background evolution and provides a friction term in the propagation of GWs, can be tested by comparing the modified GW luminosity distance with the EM luminosity distance. Considering the third-generation gravitational-wave detectors, Einstein Telescope and two cosmic explorers, we simulate a series of GW events of binary neutron stars and neutron-star-black-hole binaries with EM counterparts. These simulations can be used to constrain $f(T)$ gravity [especially the power-law model $f(T)=T+\alpha(-T)^\beta$ in this work] and other cosmological parameters, such as $\beta$ and the Hubble constant. In addition, combining simulations with current observations of type Ia supernovae and baryon acoustic oscillations, we obtain tighter limitations for $f(T)$ gravity. We find that the estimated precision significantly improved when all three datasets are combined ($\Delta \beta \sim 0.03$), compared to analyzing the current observations alone ($\Delta \beta \sim 0.3$). Simultaneously, the uncertainty of the Hubble constant can be reduced to approximately $1\%$.

We propose a model of a twisted accretion disc around a Kerr black hole interacting with a secondary black hole of a smaller mass on an inclined eccentric orbit. We use parameters of the system, which may be appropriate for the so-called 'precessing massive' model of OJ 287. We calculate expressions for torque exerted on the disc by the secondary and a contribution of the secondary to the apsidal precession of disc elements by a double averaging procedure over the periods of the secondary and the disc elements. These expressions are used at all scales of interest, including the ones inside the binary orbit. We calculate numerically the evolution of the disc tilt and twist assuming a flat initial configuration. We consider the disc aspect ratio $h/r=10^{-3}$, a rather large viscosity parameter $\alpha=0.1$ and several values of the primary rotational parameter, $\chi$. We find that, after a few periods of Lense-Thirring precession of the orbit, the disc relaxes to a quasi-stationary configuration in the precessing frame with a non-trivial distribution of the disc inclination angle, $\beta$, over the radial scale. We propose an analytic model for this configuration. We show that the presence of the twisted disc leads to multiple crossings of the disc by the secondary per one orbital period, with time periods between the crossings being different from the flat disc model. Our results should be taken into account in the modelling of OJ 287. They can also be applied to similar sources.

Methods were developed in Ref. [1] for constructing reference metrics (and from them differentiable structures) on three-dimensional manifolds with topologies specified by suitable triangulations. This note generalizes those methods by expanding the class of suitable triangulations, significantly increasing the number of manifolds to which these methods apply. These new results show that this expanded class of triangulations is still a small subset of all possible triangulations. This demonstrates that fundamental changes to these methods are needed to further expand the collection of manifolds on which differentiable structures can be constructed numerically.

Utilizing the recently established connection between Palatini-like gravity and linear Generalized Uncertainty Principle (GUP) models, we have formulated an approach that facilitates the examination of Bose gases. Our primary focus is on the ideal Bose-Einstein condensate and liquid helium, chosen as illustrative examples to underscore the feasibility of tabletop experiments in assessing gravity models. The non-interacting Bose-Einstein condensate imposes constraints on linear GUP and Palatini $f(R)$ gravity (Eddington-inspired Born-Infeld gravity) within the ranges of $-10^{12}\lesssim\sigma\lesssim 3\times 10^{24}{\text{ s}}/{\text{kg m}}$ and $-10^{-1}\lesssim\bar\beta\lesssim 10^{11} \text{ m}^2$ ($-4\times10^{-1}\lesssim\epsilon\lesssim 4\times 10^{11} \text{ m}^2$), respectively. In contrast, the properties of liquid helium suggest more realistic bounds, specifically $-10^{23}\lesssim\sigma\lesssim 10^{23}{\text{ s}}/{\text{kg m}}$ and $-10^{9}\lesssim\bar\beta\lesssim 10^{9} \text{ m}^2$. Additionally, we argue that the newly developed method employing Earth seismic waves provides improved constraints for quantum and modified gravity by approximately one order of magnitude.

Stochastic gravitational-wave (GW) background (SGWB) contains information about the early Universe and astrophysical processes. The recent evidence of SGWB by pulsar timing arrays in the nanohertz band is a breakthrough in the GW astronomy. For ground-based GW detectors, while unfortunately in data analysis the SGWB can be masked by loud GW events from compact binary coalescences (CBCs). Assuming a next-generation ground-based GW detector network, we investigate the potential for detecting the astrophysical and cosmological SGWB with non-CBC origins by subtracting recovered foreground signals of loud CBC events. As an extension of the studies by Sachdev et al. (2020) and Zhou et al. (2023), we incorporate aligned spin parameters in our waveform model. Because of the inclusion of spins, we obtain significantly more pessimistic results than the previous work, where the residual energy density of foreground is even larger than the original background. The degeneracy between the spin parameters and symmetric mass ratio is strong in the parameter estimation process and it contributes most to the imperfect foreground subtraction. Our results have important implications for assessing the detectability of SGWB from non-CBC origins for ground-based GW detectors.

We show that the seemingly different methods used to derive non-Lorentzian (Galilean and Carrollian) gravitational theories from Lorentzian ones are equivalent. Specifically, the pre-nonrelativistic and the pre-ultralocal parametrizations can be constructed from the gauging of the Galilei and Carroll algebras, respectively. Also, the pre-ultralocal approach of taking the Carrollian limit is equivalent to performing the ADM decomposition and then setting the signature of the Lorentzian manifold to zero. We use this uniqueness to write a generic expansion for the curvature tensors and construct Galilean and Carrollian limits of all metric theories of gravity of finite order ranging from the $f(R)$ gravity to a completely generic higher derivative theory, the $f(g_{\mu\nu},R_{\mu\nu\sigma \rho},\nabla_{\mu})$ gravity. We present an algorithm for calculation of the $n$-th order of the Galilean and Carrollian expansions that transforms this problem into a constrained optimization problem. We also derive the condition under which a gravitational theory becomes a modification of general relativity in both limits simultaneously.

We develop a numerical approach to compute polar parity perturbations within fully relativistic models of black hole systems embedded in generic, spherically symmetric, anisotropic fluids. We apply this framework to study gravitational wave generation and propagation from extreme mass-ratio inspirals in the presence of several astrophysically relevant dark matter models, namely the Hernquist, Navarro-Frenk-White, and Einasto profiles. We also study dark matter spike profiles obtained from a fully relativistic calculation of the adiabatic growth of a BH within the Hernquist profile, and provide a closed-form analytic fit of these profiles. Our analysis completes prior numerical work in the axial sector, yielding a fully numerical pipeline to study black hole environmental effects. We study the dependence of the fluxes on the DM halo mass and compactness. We find that, unlike the axial case, polar fluxes are not adequately described by simple gravitational-redshift effects, thus offering an exciting avenue for the study of black hole environments with gravitational waves.

We report unbiased AI measurements of the fine structure constant alpha in two proximate absorption regions in the spectrum of the quasar HE0515-4414. The data are high resolution, high signal to noise, and laser frequency comb calibrated, obtained using the ESPRESSO spectrograph on the VLT. The high quality of the data and proximity of the regions motivate a differential comparison, exploring the possibility of spatial variations of fundamental constants, as predicted in some theories. We show that if the magnesium isotopic relative abundances are terrestrial, the fine structure constants in these two systems differ at the 7-sigma level. A 3-sigma discrepancy between the two measurements persists even for the extreme non-terrestrial case of 100% ^{24}Mg, if shared by both systems. However, if Mg isotopic abundances take independent values in these two proximate systems, one terrestrial, the other with no heavy isotopes, both can be reconciled with a terrestrial alpha, and the discrepancy between the two measurements falls to 2-sigma. We discuss varying constant and varying isotope interpretations and resolutions to this conundrum for future high precision measurements.

Electrostatic force actuation is a key component of the system of geodesic reference test masses (TM) for the LISA orbiting gravitational wave observatory and in particular for performance at low frequencies, below 1 mHz, where the observatory sensitivity is limited by stray force noise. The system needs to apply forces of order 10$^{-9}$ N while limiting fluctuations in the measurement band to levels approaching 10$^{-15}$ N/Hz$^{1/2}$. We present here the LISA actuation system design, based on audio-frequency voltage carrier signals, and results of its in-flight performance test with the LISA Pathfinder test mission. In LISA, TM force actuation is used to align the otherwise free-falling TM to the spacecraft-mounted optical metrology system, without any forcing along the critical gravitational wave-sensitive interferometry axes. In LISA Pathfinder, on the other hand, the actuation was used also to stabilize the TM along the critical $x$ axis joining the two TM, with the commanded actuation force entering directly into the mission's main differential acceleration science observable. The mission allowed demonstration of the full compatibility of the electrostatic actuation system with the LISA observatory requirements, including dedicated measurement campaigns to amplify, isolate, and quantify the two main force noise contributions from the actuation system, from actuator gain noise and from low frequency ``in band'' voltage fluctuations. These campaigns have shown actuation force noise to be a relevant, but not dominant, noise source in LISA Pathfinder and have allowed performance projections for the conditions expected in the LISA mission.

We investigate the tachyonic instability of Kerr-Newman (KN) black hole with a rotation parameter $a$ in the Einstein-Chern-Simons-scalar theory coupled with a quadratic massive scalar field. This instability analysis corresponds to exploring the onset of spontaneous scalarization for KN black holes. First, we find no $a$-bound for $\alpha<0$ case by considering (1+1)-dimensional analytical method. A direct numerical method is adopted to explore (2+1)-dimensional time evolution of a massive scalar perturbation with positive and negative $\alpha$ to obtain threshold curves numerically. We obtain threshold curves $\alpha_{\rm th}(a)$ of tachyonic instability for positive $\alpha$ without any $a$-bounds. We expect to find the same threshold curves $\alpha_{\rm th}(a)$ of tachyonic instability for negative $\alpha$ without any $a$-bound because its linearized scalar theory is invariant under the transformation of $\alpha\to -\alpha $ and $\theta\to -\theta$. In addition, it is found that the scalar mass term suppresses tachyonic instability of KN black holes.