Papers for Tuesday, Dec 24 2024

The efficiency of planet formation is a fundamental question in planetary science, gaining increasing significance as observational data from planet-forming disks accumulates. Here we derive from first principles a correlation between the planet formation rate (PFR) and the gas surface density, i.e. $\rm{PFR}\propto \Sigma_g^n$. This relation serves as an analog for the well-established Kennicutt-Schmidt law for star-forming galaxies. We study the different planet formation mechanisms and the density dependence in each one of them, to finally formulate a simple relation. We find that the powerlaw ranges between $n\approx 4/3-2$, depending on the type of the forming planet, when we carry out different analyses for the formation rates of terrestrial planets, gas giants, and also planets formed by gravitational instability. We then compare our results with the available observational data. The relation we derive here aims to shed more light on the interpretation of observational data as well as analytical models, and give a new perspective on the properties of planet formation and its connection to gas.

We investigate energy release in the interacting magnetospheres of binary neutron stars (BNS) with global 3D force-free electrodynamics simulations. The system's behavior depends on the inclinations $\chi_1$ and $\chi_2$ of the stars' magnetic dipole moments relative to their orbital angular momentum. The simplest aligned configuration ($\chi_1=\chi_2=0^\circ$) has no magnetic field lines connecting the two stars. Remarkably, it still develops separatrix current sheets warping around each star and forms a dissipative region at the interface of the two magnetospheres. Dissipation at the interface is caused by a Kelvin-Helmholtz-type (KH) instability, which generates local Alfvénic turbulence and escaping fast magnetosonic waves. Binaries with inclined magnetospheres release energy in two ways: via KH instability at the interface and via magnetic reconnection flares in the twisted flux bundles connecting the two stars. The outgoing compressive waves can develop shocks and source fast radio bursts. We discuss implications for X-ray and radio precursors of BNS mergers.

Constraining the physical and chemical evolution of molecular clouds is essential to our understanding of star formation. These investigations often necessitate knowledge of some local representative number density of the gas along the line of sight. However, constraining the number density is a difficult endeavor. Robust constraints of the number density often require line observations of specific molecules along with radiation transfer modeling, which provides densities traced by that specific molecule. Column density maps of molecular clouds are more readily available, with many high-fidelity maps calculated from dust emission and extinction, in particular from surveys conducted with the Herschel Space Observatory. We introduce a new probabilistic model which is based on the assumption that the total hydrogen nuclei column density along a line of sight can be decomposed into a turbulent component and a gravitationally-dominated component. Therefore, for each pixel in a column density map, the line of sight is decomposed into characteristic diffuse (dubbed ``turbulent'') and dense (dubbed ``gravitational'') gas number densities from column density maps. The method thus exploits a physical model of turbulence to decouple the random turbulent column from gas in dense bound structures empirically using the observed column density maps. We find the model produces reasonable turbulent and gravitational densities in the Taurus L1495/B213 and Polaris Flare clouds. The model can also be used to infer an effective attenuating column density into the cloud, which is useful for astrochemical models of the clouds. We conclude by demonstrating an application of this method by predicting the emission of the [C II], [C I], and CO (J = 1-0) lines across the Taurus L1495/B213 region at the native resolution of the column density map utilizing a grid of photodissociation-region models.

Models proposing a non-gravitational interaction between dark energy (DE) and dark matter (CDM) have been extensively studied as alternatives to the standard cosmological model. A common approach to describing the DE-CDM coupling assumes it to be linearly proportional to the dark energy density. In this work, we consider the model with interaction term $Q=3H\gamma{\rho_{x}^{2}}/{(\rho_{c}+\rho_{x})}$. We show that for positive values of $\gamma$ this model predicts a future violation of the Weak Energy Condition (WEC) for the dark matter component, and for a specific range of negative values of $\gamma$ the CDM energy density can be negative in the past. We perform a parameter selection analysis for this model using data from Type Ia supernovae, Cosmic Chronometers, Baryon Acoustic Oscillations, and CMB combined with the Hubble constant $H_0$ prior. Imposing a prior to ensure that the WEC is not violated, our model is consistent with $\Lambda$CDM in 2$\sigma$ C.L.. In reality, the WEC prior shifts the constraints towards smaller values of $H_0$, highlighting an increase in the tension on the Hubble parameter. However, it significantly improves the parameter constraints, with a preference for smaller values of $\sigma_8$, alleviating the $\sigma_8$ tension between the CMB results from Planck 2018 and the weak gravitational lensing observations from the KiDS-1000 cosmic shear survey. In the case without the WEC prior, our model seems to alleviate the $H_0$ tension, which is related to the positive value of the interaction parameter $\gamma$.

We extend the co-scaling formalism of Habegger & Heitsch (2021) implemented in Athena++ to magneto-hydrodynamics. The formalism relies on flow symmetries in astrophysical problems involving expansion, contraction, and center-of-mass motion. The formalism is fully consistent with the upwind constrained transport method implemented in Athena++ and is accurate to 2nd order in space. Applying our implementation to standard magneto-hydrodynamic test cases leads to improved results and higher efficiency, compared to the fixed-grid solutions.

The Arcus Probe mission concept provides high-resolution soft X-ray and UV spectroscopy to reveal feedback-driven structure and evolution throughout the universe with an agile response capability ideal for probing the physics of time-dependent phenomena. The X-ray Spectrograph (XRS) utilizes two nearly identical CCD focal planes to detect and record X-ray photons from the dispersed spectra and zero-order of the critical angle transmission gratings. In this paper we describe the Arcus focal plane instrument and the CCDs, including laboratory performance results, which meet observatory requirements.

The R Coronae Borealis (RCB) variables are rare, hydrogen-deficient, carbon-rich supergiants known for large, erratic declines in brightness due to dust formation. Recently, the number of known RCB stars in the Milky Way and Magellanic Clouds has increased from $\sim$30 to 162. We use all-sky and targeted photometric surveys to create the longest possible light curves for all known RCB stars and systematically study their declines. Our study, the largest of its kind, includes measurements of decline activity levels, morphologies, and periodicities for nearly all RCB stars. We confirm previous predictions that cool RCB stars exhibit more declines than warm RCBs, supporting a relationship between dust formation and condensation temperatures. We also find evidence for two distinct dust production mechanisms. R CrB and SU Tau show decline onsets consistent with a Poisson process, suggesting their dust production is driven by stochastic processes, such as convection. In contrast, RY Sgr's declines correlate with its pulsation period, suggesting that its dust production is driven by pulsationally-induced shocks. Finally, we show that the dust properties of the related class of DY~Per variables differ from those of the RCB stars, suggesting differences in their evolutionary status.

Dwarf galaxies are thought of as the building blocks of large galaxies such as our Milky Way. This paper presents new high-resolution hydrodynamical simulations of dwarf galaxies and their intergalactic medium with the \texttt{GIZMO} code. Our simulations consider the key physical processes of galaxy evolution, such as gas cooling, chemistry, and stellar and black hole feedback. Unlike the previous work, the initial conditions of our simulations taking the dwarf galaxies of $2-5 \times 10^{10} \, M_\odot$ from the realistic cosmology simulations, \texttt{IllustrisTNG}. We further increase the original resolution of \texttt{IllustrisTNG} by a factor of $\sim 100$ via a particle splitting scheme. Our results show that the evolution of complex multiphase CGM and its metal content is sensitive to the redshift of dwarf galaxies. The accretion of CGM into dwarf galaxies plays a key role in providing $20 \% - 50 \%$ of the star-forming gas and replenishing $40 \% - 70 \%$ of the total mass in the galactic disk. Furthermore, the accretion history of supermassive black holes in the centers of high-$z$ dwarf galaxies shows episodic patterns with high-accreting states close to $\sim 10 \%$ of the Eddington mass accretion rate, implying the rapid growth of supermassive black holes in the early universe, which may be revealed by the coming observations from the James Webb Space Telescope (JWST).

We conduct a thorough study of the comoving curvature perturbation $\mathcal{R}$ in single-field inflation with two stages, represented by a piecewise quadratic potential, where both the first and second derivatives are allowed to be discontinuous at the transition point. We calculate the evolution of $\mathcal{R}$ by combining the perturbative and non-perturbative methods consistently, and obtain the power spectrum and the non-Gaussian features in the probability distribution function. We find that both the spectrum and the statistics of $\mathcal{R}$ depend significantly on the second derivatives of the potential at both the first and second stages. Furthermore, we find a new parameter constructed from the potential parameters, which we call $\alpha$, plays a decisive role in determining various features in the spectrum such as the amplitude, the slope, and the existence of a dip. In particular, we recover the typical $k^4$ growth of the spectrum in most cases, but the maximum growth rate of $k^5(\log k)^2$ can be obtained by fine-tuning the parameters. Then, using the $\delta N$ formalism valid on superhorizon scales, we give fully nonlinear formulas for $\cal{R}$ in terms of the scalar field perturbation $\delta\phi$ and its time derivative. In passing, we point out the importance of the nonlinear evolution of $\delta\phi$ on superhorizon scales. Finally, using the Press-Schechter formalism for simplicity, we discuss the effect of the non-Gaussian tails of the probability distribution function on the primordial black hole formation.

Thermal instability (TI) is a trigger mechanism, which can explain the formation of small condensations through some regions of the interstellar clouds. The instability criterion for flat geometry approximations has been investigated in previous works. Here, we focus on spherical perturbations in the spherical clouds. Our goal here is to examine the conditions for the occurrence of TI through the thermally dominated (i.e., gravitationally stable) quasi-static spherical interstellar clouds. First, we obtain the profiles of density, temperature, pressure, and enclosed mass of a symmetric spherical cloud. Then, we use the perturbation method to investigate the linear regime of instability and find its growth rate. Considering spherical perturbations on the quasi-static spherical cloud, instead of a thermal and dynamical equilibrium flat cloud, changes the instability criterion so that we can conclude that sphericalness can increase the occurrence of TI. The results show that in the spherical clouds, perturbations with shorter wavelengths have more chance to grow via TI (i.e., greater growth rates).

Stellar surface heterogeneities, such as spots and faculae, often contaminate exoplanet transit spectra, hindering precise atmospheric characterization. We demonstrate a novel, epoch-based, model-independent method to mitigate stellar contamination, applicable to multi-planet systems with at least one airless planet. We apply this method using quasi-simultaneous transits of TRAPPIST-1 b and TRAPPIST-1 c observed on July 9, 2024, with JWST NIRSpec PRISM. These two planets, with nearly identical radii and impact parameters, are likely either bare rocks or possess thin, low-pressure atmospheres, making them ideal candidates for this technique, as variations in their transit spectra would be primarily attributed to stellar activity. Our observations reveal their transit spectra exhibit consistent features, indicating similar levels of stellar contamination. We use TRAPPIST-1 b to correct the transit spectrum of TRAPPIST-1 c, achieving a 2.5x reduction in stellar contamination at shorter wavelengths. At longer wavelengths, lower SNR prevents clear detection of contamination or full assessment of mitigation. Still, out-of-transit analysis reveals variations across the spectrum, suggesting contamination extends into the longer wavelengths. Based on the success of the correction at shorter wavelengths, we argue that contamination is also reduced at longer wavelengths to a similar extent. This shifts the challenge of detecting atmospheric features to a predominantly white noise issue, which can be addressed by stacking observations. This method enables epoch-specific stellar contamination corrections, allowing co-addition of planetary spectra for reliable searches of secondary atmospheres with signals of 60-250 ppm. Additionally, we identify small-scale cold (2000 K) and warm (2600 K) regions almost uniformly distributed on TRAPPIST-1, with overall covering fractions varying by 0.1% per hour.

In this paper we consider the effect of quantum entanglement and coherence on the radiated intensity from a gas and its absorption capacity at thermal equilibrium or, more generally, under conditions where no population inversion exists. As was shown by Dicke (1954), although entangled states and coherence can lead to superradiance for specific modes of radiation, they can also bring subradiance through significant energy trapping in slow and dark states. While a finite separation between the atoms composing the gas will cause leaking of the trapped energy, we show how the combination of thermal equilibrium and quantum coherence mitigates this effect and leads to significantly reduced radiation intensity from the gas, rendering it dark and collision-less. Furthermore, we show how under the same conditions absorption of a radiation field incident on the gas can lead to higher attenuation levels than those predicted with Beer's law. Beer's law is recovered in the limit of complete decoherence. We apply our analysis to the atomic hydrogen 21 cm line and, considering the gas densities expected in Dark Matter halos, we find that quantum entanglement and coherence can potentially account for some of the Dark Matter known exist in these environments.

We present astrometric results from a Hubble Space Telescope (HST) campaign aimed at determining precise distances for cold Y-type brown dwarfs. Combining observations from a dedicated HST/WFC3 program with archival data, we derive astrometric solutions for 15 nearby Y dwarfs, by linking the high-precision relative astrometry from Hubble to the high-accuracy Gaia DR3 absolute reference system, using stars present in both to anchor the two frames of reference. We reach uncertainties on parallaxes below the 1-mas level for half of the sample, and down to 3 mas for two thirds of the targets, or relative precisions <1% in most cases and 2-5x improvements over previous measurements. For the remaining targets, we achieved slightly lower precisions on parallaxes (5-12 mas, 5-10%), correlated with the lower signal-to-noise of the faintest targets. The precision reached in our derived proper motions is around 0.1-0.4 mas/yr for most targets, and up to 1-2 mas/yr for less precise cases. Our estimated parallaxes and proper motions are generally in good agreement with literature values, and consistent to 1-2 sigma with recent Spitzer-derived parallaxes in most cases. These new astrometric solutions provide important validation of these objects' distances and sky motions, especially given the large disparities seen in previous estimates. Our results demonstrate the power of HST combined with Gaia to measure highly-precise absolute astrometry for faint brown dwarfs, and highlights the limitations reached for the reddest and coldest objects, for which JWST will certainly provide a favourable platform to improve these results.

In the context of kilonova (KN) modeling, the present work focusses on large-scale atomic data and opacity computations for all heavy elements from Ca to Lr, with a special effort on lanthanides and actinides, for a grid of typical KN ejecta conditions between one day and one week after the merger (corresponding to the LTE photosphere phase of the KN ejecta). In order to do so, we used the pseudo-relativistic Hartree-Fock (HFR) method, in which the choice of the interaction configuration model is of crucial importance. In this paper, HFR atomic data and opacities for all elements between Ca (Z = 20) and Lr (Z = 103) are presented, with a special focus on lanthanides and actinides. Besides, we also discuss the contribution of every single element to the total KN ejecta opacity for a given neutron star merger model, depending on their Planck mean opacities and elemental abundances. An important result is that lanthanides are found to not be the dominant sources of opacity, at least on average. The impact on KN light curves of considering such atomic-physics based opacity data instead of typical crude approximation formulae is also evaluated. In addition, the importance of taking the ejecta composition into account directly in the expansion opacity determination (instead of estimating single-element opacities) is highlighted. A database containing all the relevant atomic data and opacity tables has also been created and published online along with this work.

Modeling the orbital dynamics of objects in galactic disks is crucial to understanding the stability and evolution of disk galaxies. While studies of galactic orbits are largely dominated by $N$-body simulations, perturbative analytical models offer a computationally inexpensive and conceptually insightful way of analyzing galactic dynamics. We utilize perturbation theory and the method of matched asymptotics to develop a technique by which the vertical motion of a point mass perpendicular to a thin axisymmetric disk galaxy can be computed analytically to high precision. The objective of this study is to provide an accurate model of non-planar dynamics that could be employed in galactic simulations to bypass computationally expensive integration. We construct a general solution for the dynamics at small $z$ displacements from the plane of the disk using perturbation theory. We solve for the dynamics at very large displacements from the plane of the disk and construct a function that interpolates the solutions at each length scale. To demonstrate the method's utility, we apply it to two commonly used models for the vertical structure of a galactic disk, the exponential model and the isothermal model. Finally, using the principle of adiabatic invariance, we analyze how to relate variations in radial motion to variations in the vertical motion and discuss possible applications of the model in numerical astrophysics.

Glitches are non-Gaussian noise transients originating from environmental and instrumental sources that contaminate data from gravitational wave detectors. Some glitches can even mimic gravitational wave signals from compact object mergers, which are the primary targets of terrestrial observatories. In this study, we present a method to analyze noise transients from the LIGO observatories using Q-transform information combined with t-Distributed Stochastic Neighbor Embedding (t-SNE). We implement classification techniques, examine the influence of parameters on glitch classification, and conduct a week-long daily analysis to track outlier transients over time.

Neutral hydrogen serves as a crucial probe for the Cosmic Dawn and the Epoch of Reionization (EoR). Actual observations of the 21-cm signal often encounter challenges such as thermal noise and various systematic effects. To overcome these challenges, we simulate SKA-Low-depth images and process them with a deep learning method. We utilized foreground residuals acquired by LOFAR during actual North Celestial Pole field observations, thermal and excess variances calculated via Gaussian process regression, and 21-cm signals generated with 21cmFAST for signal extraction tests. Our approach to overcome these foreground, thermal noise, and excess variance components employs a 3D U-Net neural network architecture for image analysis. When considering thermal noise corresponding to 1400 hours of integration, U-Net provides reliable 2D power spectrum predictions, and robustness tests ensure that we get realistic EoR signals. Adding foreground residuals, however, causes inconsistencies below the horizon delay-line. Lastly, evaluating both thermal and excess variance with observations up to 3700 and 14000 hours ensures reliable power spectrum estimations within the EoR window and across nearly all scales, respectively.

In our previous work, we have investigated Galactic cosmic ray (GCR) spectra and anisotropy from 100 GeV to PeV, under anisotropic propagation model with axisymmetric distributed galactic sources. Numerous observational evidence have indicated that the Milky Way is a typical spiral galaxy. In this work, we further utilize anisotropic propagation models with spiral galactic sources to investigate spectra and anisotropy of CRs. During the calculation process, we employ the spatially dependent diffusion (SDP) model with different diffusion coefficients for the inner and outer halo, while the background CR sources is spiral distribution. To better explain the anomalous observations of nuclear spectral hardening at $ {\cal R}\sim$ 200 GV and the complicated energy dependence of anisotropy from GeV to PeV, we introduce the contribution of the Geminga nearby source. Additionally, we incorporate the impact of the local regular magnetic field (LRMF) and the corresponding anisotropic diffusion on large-scale anisotropy within the SDP model. By comparing the spiral and axisymmetric distribution models of background sources, it is found that both of them can well reproduce the CR spectra and anisotropy from 100 GeV-PeV. However, their propagation parameters are different. The diffusion coefficient with spiral distribution of CR sources is larger than that with axisymmetric distribution, and its spectral indices are also harder. Future high-precision measurements of CR anisotropy, such as LHAASO experiment, will be crucial in evaluating the validity of our proposed model.

Gamma-ray bursts (GRBs) have long been proposed as a potential source of high-energy neutrinos. Although no confirmed association between GRBs and neutrinos has been established, meaningful constraints have been placed on GRB prompt emission models. The non-detection of neutrinos, reported by the IceCube collaboration, from both single and stacked GRB events suggests that the radiation zone is likely located at a considerable distance from the central engine, where the photon number density is relatively low. Here, we estimate future GRB-neutrino detection probabilities with more sensitive detectors than IceCube and explore the constraints on models if GRB neutrinos remain undetected despite improved sensitivity. Our findings reveal that if the effective area of a future neutrino detector can be enhanced by a factor of 10 compared to IceCube IC86-II, there is a high likelihood of detecting neutrinos from a GRB 221009A-like event, even in the context of the ICMART model, which exhibits the lowest efficiency in neutrino production. With such an advanced detector (enhanced by a factor of 10) and 5 to 10 years of data accumulation, neutrinos from stacked GRBs should be identifiable, or several popular models for GRB prompt emission (e.g., the dissipative photosphere model and internal shock model) could be effectively ruled out.

Aims. We reported the discovery that a changing-look AGN SDSS J101152.98+544206.4 (J1011+5442 for short) gradually returns to the type 1 state after a short period between 2014 and 2019 in the faint type 1.9 state. Methods. Motivated by the rebrightening in optical and mid-infrared light curves from ZTF and WISE, we obtained the new spectroscopic observations by Xinglong 2.16-m, Lijiang 2.4-m, and MMT 6.5-m optical telescopes in 2024. Results. After changing the optical AGN type from 1 to 1.9 between 2003 and 2015 based on the repeat spectroscopy from the Time Domain Spectroscopic Survey, J1011+5442 returns to its type 1 state in 2024. We detect the significant and very broad Hbeta lines (FWHM > 5000 km/s) based on the new spectra, which suggests that J1011+5442 is in the intermediate state between the dim state in 2015 and the bright state in 2003. The long-term optical and mid-infrared light curves also show a brightening trend between 2019 and 2024 as the broad Hbeta line appears. The time lag of about 100 days between the mid-infrared and optical variability is consistent with the prediction of dust reverberation mapping. Conclusions. The behaviors of the photometric and spectroscopic observations of J1011+5442 are consistent with the argument that the repeating changing-look phenomenon is regulated by the variation of accretion rate.

We compute realistic 3D radiative MHD near-surface models of starspots with substantial penumbrae on cool main-sequence stars using the MURaM simulation code. This work is an improvement on the the previous starspot models in a slab geometry. The umbra, penumbra and the quiet star for all starspots are distinct, not only in intensity and temperature, but also in thermodynamic and velocity structure. These models represent a significant step towards modelling contribution of starspots to stellar lightcurves.

The stellar atmospheric parameters and physical properties of stars in the Kepler Input Catalog (KIC) are of great significance for the study of exoplanets, stellar activity, and asteroseismology. However, despite extensive effort over the past decades, accurate spectroscopic estimates of these parameters are available for only about half of the stars in the full KIC catalog. In our work, by training relationships between photometric colors and spectroscopic stellar parameters from Gaia DR3, the Kepler Issac-Newton Survey, LAMOST DR10, and APOGEE DR17, we have obtained atmospheric-parameter estimates for over 195,000 stars, accounting for 97% of the total sample of KIC stars. We obtain 1{\sigma} uncertainties of 0.1 dex on metallicity [Fe/H], 100 K on effective temperature $T_{\mathrm{eff}}$ , and 0.2 dex on surface gravity log $g$. In addition, based on these atmospheric parameters, we estimated the ages, masses, radii, and surface gravities of these stars using the commonly adopted isochrone-fitting approach. The resulting precisions are 20% for ages, 0.1 $M_{\odot}$ for masses, 0.01 $R_{\odot}$ for radii and 0.1 dex for surface gravities. These accurate parameters are expected to provide valuable insights for future studies on various fields.

We present a new description of the 7.7~$\mu$m region towards the high-mass star-forming region IRAS 23385+6053 taken from open James Webb Space Telescope Mid-Infrared Instrument Medium Resolution Spectrometer (JWST MIRI/MRS) data. This area is commonly attributed to the $\nu_4$ deformation mode of methane ice. For the first time gaseous and solid methane were analyzed simultaneously in IRAS 23385+6053. The band at 7.58--7.8 $\mu$m (1320--1280 cm$^{-1}$) is interpreted as a wide solid absorption methane feature overlapped by the sharp features of the methane emission. We report the detection of gaseous methane and estimate its emitting area radius~$R$, temperature~$T$ and column density~$N$ as $R=2940$~au, $T=103^{+13}_{-11}$~K, and $N=0.78_{+6.18}^{-0.64}\times10^{17}$~cm$^{-2}$, correspondingly. The ice content was analyzed with the laboratory spectra dataset of methane in different molecular environments obtained on the Ice Spectroscopy Experimental Aggregate (ISEAge). We were able to describe the wide feature of solid methane with the following laboratory spectra: CH$_4$~:~CO$_2$~=~1~:~5 (at $27.4^{+6.0}_{-10.8}$~K) and CH$_4$~:~H$_2$O~=~1~:~10 (at $8.4^{+16.4}_{-1.7}$~K) deposited at 6.7~K and warmed up at a rate of 0.5 K per minute. The derived column densities are $N_{\text{CH}_4}$(CO$_2$)~=~$2.97^{+0.37}_{-0.57}\times10^{17}$~cm$^{-2}$ and $N_{\text{CH}_4}$(H$_2$O)~=~$1.02^{+0.46}_{-0.27}\times10^{17}$~cm$^{-2}$. According to the best fit solid methane is mostly surrounded by CO$_2$ rather than H$_2$O. The residuals analysis reveals the unassigned region at 1283--1297~cm$^{-1}$ (7.71--7.79~$\mu$m) which is tentatively assigned to nitrous oxide (N$_2$O) in various environments.

We report a high-redshift ($z=1.404$) tidal disruption event (TDE) candidate in SDSS J000118.70+003314.0 (SDSS J0001), which is a quasar with apparent broad Mg~{\sc ii} emission line. The long-term variability in its nine-year photometric $ugriz$-band light curves, obtained from the SDSS Stripe82 and the PHOTOOBJALL databases, can be described by the conventional TDE model. Our results suggest that the TDE is a main-sequence star with mass of $1.905_{-0.009}^{+0.023}{\rm M_\odot}$ tidally disrupted by a black hole (BH) with mass {$6.5_{-2.6}^{+3.5}\times10^7{\rm M_\odot}$}. The BH mass is about 7.5 times smaller than the virial BH mass derived from the broad Mg~{\sc ii} emission line, which can be explained by non-virial dynamic properties of broad emission lines from TDEs debris. Furthermore, we examine the probability that the event results from intrinsic variability of quasars, which is about $0.009\%$, through applications of the DRW/CAR process. Alternative explanations for the event are also discussed, such as the scenarios of dust obscurations, microlensing and accretion. Our results provide clues to support that TDEs could be detectable in broad line quasars as well as in quiescent galaxies, and to indicate the variability of some active galactic nuclei may be partly attributed to central TDEs.

In the era of exoplanet studies with JWST, the transiting, hot gas giant WASP-52 b provides an excellent target for atmospheric characterization through transit spectroscopy. WASP-52 b orbits an active K-type dwarf recognized for its surface heterogeneities, such as star-spots and faculae, which offers challenges to atmospheric characterization via transmission spectroscopy. Previous transit observations have detected active regions on WASP-52 through crossing events in transit light-curves and via the spectral imprint of unocculted magnetic regions on transmission spectra. Here, we present the first JWST observations of WASP-52 b. Our JWST NIRISS/SOSS transit observation, obtained through the GTO 1201 Program, detects two clear spot-crossing events that deform the 0.6-2.8 $\mu$m transit light-curves of WASP-52 b. We find that these two occulted spots combined cover about 2.4 % of the stellar surface and have temperatures about 400-500 K colder than the stellar photosphere. Our NIRISS/SOSS transmission spectrum is best-fit by an atmosphere with H$_2$O (10.8 $\sigma$), He (7.3 $\sigma$, with evidence of an escaping tail at $\sim$ 2.9 $\sigma$), hints of K (2.5 $\sigma$), and unocculted star-spots and faculae (3.6 $\sigma$). The retrieved H$_2$O abundance ($\log$ H$_2$O $\approx -4 \pm 1$) is consistent with a subsolar or solar atmospheric metallicity for two independent data reductions. Our results underscore the importance of simultaneously modelling planetary atmospheres and unocculted stellar heterogeneities when interpreting transmission spectra of planets orbiting active stars and demonstrate the necessity of considering different stellar contamination models that account for both cold and hot active regions.

Finding Dyson rings around distant pulsars may involve identifying light curve features that have not been previously identified. Previous studies covered the detection of a ring structure uniformly brightened by the central pulsar, mostly in infrared light. Here, more complex light curves are explored, which arise inherently from the pulsar beam spot's commonly predicted superluminal speed. These speeds may cause multiple images of the pulsar's spot on the Dyson ring to appear simultaneously to a distant observer, and so feature bright creation and annihilation events. Therefore, it is possible that even if Dyson ring structures had been observed previously, they might have remained unnoticed. Similar light curve features may appear on naturally occurring dust rings around pulsars that reflect detectable pulsar radiation.

We study the general relativistic transonic accretion flow around the primary black hole, which forms the circumprimary disc (CPD), within a BBH system. The BBH spacetime is characterized by the mass ratio ($q$) and the separation distance ($z_2$) between the two black holes. We numerically solve the radial momentum and energy equations to obtain the accretion solutions. It is observed that the CPD can exhibit shock solutions, which exist for a wide range parameter space spanned by flow angular momentum ($\lambda$) and energy ($E$). We find that the shock parameter space is modified by $q$ and $z_2$. Investigations show that $q$ and $z_2$ also affect various shock properties, such as density compression and temperature compression across the shock fronts. Moreover, we calculate the spectral energy distributions (SEDs) of the CPD and examine how the SEDs are modified by $q$ and $z_2$ for both shock-free and shock-induced accretion solutions. SED is found to be nearly independent of the binary parameters. We depict that the mutual interaction between two black holes in a binary system hardly alters their individual accretion features. This study replicates the existing predictions in a different context of binary accretion.

We investigate the emission of vector radiation by superconducting cosmic string loops, deriving general relations to characterize the vector radiation emission efficiency, and study its impact on the evolution of loops. Building on these results, we compute the stochastic gravitational wave background generated by a chiral superconducting cosmic string network. Our analysis reveals that strong coupling between superconducting cosmic strings and the vector field may lead to a substantial suppression of the gravitational wave signal, while moderate coupling may still produce a detectable signal. We demonstrate that, in this intermediate limit, the presence of superconductivity in cosmic strings may help reconcile their gravitational wave spectrum with pulsar timing array data for large enough values of current.

The direct observation of intermediate-mass black holes (IMBH) populations would not only strengthen the possible evolutionary link between stellar and supermassive black holes, but unveil the details of the pair-instability mechanism and elucidate their influence in galaxy formation. Conclusive observation of IMBHs remained elusive until the detection of gravitational-wave (GW) signal GW190521, which lies with high confidence in the mass gap predicted by the pair-instability mechanism. Despite falling in the sensitivity band of current GW detectors, IMBH searches are challenging due to their similarity to transient bursts of detector noise, known as glitches. In this proof-of-concept work, we combine a matched-filter algorithm with a Machine Learning (ML) method to differentiate IMBH signals from non-transient burst noise, known as glitches. In particular, we build a multi-layer perceptron network to perform a multi-class classification of the output triggers of matched-filter. In this way we are able to distinguish simulated GW IMBH signals from different classes of glitches that occurred during the third observing run (O3) in single detector data. We train, validate, and test our model on O3a data, reaching a true positive rate of over $90\%$ for simulated IMBH signals. In O3b, the true positive rate is over $70\%$. We also combine data from multiple detectors to search for simulated IMBH signals in real detector noise, providing a significance measure for the output of our ML method.

Searches for annihilating dark matter are often designed with a specific dark matter candidate in mind. However, the space of potential dark matter models is vast, which raises the question: how can we search for dark matter without making strong assumptions about unknown physics. We present a model-independent approach for measuring dark matter annihilation ratios and branching fractions with $\gamma$-ray event data. By parameterizing the annihilation ratios for seven different channels, we obviate the need to search for a specific dark matter candidate. To demonstrate our approach, we analyse simulated data using the \texttt{GammaBayes} pipeline. Given a 5$\sigma$ signal, we reconstruct the annihilation ratios for five dominant channels to within 95% credibility. This allows us to reconstruct dark matter annihilation/decay channels without presuming any particular model, thus offering a model-independent approach to indirect dark matter searches in $\gamma$-ray astronomy. This approach shows that for masses between 0.3-5 TeV we can probe values below the thermal relic velocity annihilation weighted cross-section allowing a 2$\sigma$ detection for 525 hours of simulated observation data by the Cherenkov Telescope Array Observatory of the Galactic Centre.

Significant progress has been made over the past decades towards unveiling the sources of the most energetic particles in nature, the ultra-high-energy cosmic rays (UHECRs). Despite these advancements, the exact astrophysical sites capable of accelerating these particles to such extreme energies remain largely unknown. Moreover, the mechanisms by which they achieve these extreme energies are poorly understood. Here, I provide a concise overview of the theory underlying the acceleration and propagation of UHECRs. I then critically discuss three recent results that could help unveil their origins: the reported excess around Centaurus A, the correlation with starburst galaxies, and the efforts to jointly model the energy spectrum, composition, and arrival directions. Finally, I discuss strategies for advancing this field, emphasising the need for refined theoretical models, the challenges in building them, and the potential for new observatories to shed light on the mysteries of UHECRs.

X-ray Astrophysics, which addresses extreme physics in extreme conditions, is particularly well suited for answering questions related to known physics. Reversely tiny effects, but integrated along sidereal distances, allow to probe extensions of known physics or even new physics. The new window into polarimetry in this energy band, opened by the Imaging X-ray Polarimetry Explorer (IXPE) a NASA-ASI Small Explorer mission launched on 9th December 2021 enables an entirely novel approach, whether used alone or in combination with standard observables such as light curves and spectra and with data in other wavelengths. In this paper, we review IXPE's results after nearly three years of successful operation, focusing on their implications for key questions in Fundamental Physics.

Accreting X-ray pulsars, located in X-ray binaries, are neutron stars with magnetic fields as strong as $B\sim10^{12\text{--}13}$ G. This review offers a concise overview of the accretion and radiation processes of X-ray pulsars and summarizes their rich observational features, particularly focusing on complex and variable temporal phenomena, spectral properties, and evolution, the new window for X-ray polarimetry and multi-wavelength advances. We also briefly discuss other related systems, i.e., gamma-ray binaries and pulsating ultraluminous X-ray sources.

V445 Puppis is the only known example of a helium nova, where a layer of helium-rich gas accretes onto the surface of a white dwarf in a cataclysmic variable, with runaway helium burning making for the nova event. Speculatively, helium nova can provide one path to produce a Type Ia supernova (SNIa), within the larger framework of single-degenerate models. Relatively little has been known about V445 Pup, with this work reporting the discovery of the orbital period near 1.87 days. The companion star is 2.65$\pm$0.35 R$_{\odot}$ in radius as an evolved giant star stripped of its outer hydrogen envelope. The orbital period immediately before the 2000 eruption was $P_{\rm pre}$=1.871843$\pm$0.000014 days, with a steady period change of (-0.17$\pm$0.06)$\times$10$^{-8}$ from 1896--1995. The period immediately after the nova eruption was $P_{\rm post}$=1.873593$\pm$0.000034 days, with a $\dot{P}$ of ($-$4.7$\pm$0.5)$\times$10$^{-8}$. The fractional orbital period change ($\Delta P/P$) is $+$935$\pm$27 ppm. This restricts the mass of the gases ejected in the nova eruption to be $\gg$0.001M$_{\odot}$, and much greater than the mass accreted to trigger the nova. So the white dwarf is losing mass over each eruption cycle, and will not become a SNIa. Further, for V445 Pup and helium novae in general, I collect observations from 136 normal SNIa, for which any giant or sub-giant companion star would have been detected, yet zero companions are found. This is an independent proof that V445 Pup and helium novae are not SNIa progenitors.

Aims. This study aims to investigate the deviation of the intensity ratio of the \ion{Si}{IV} 1394 Å and 1403 Å emission lines from the expected value of 2 in the optically thin regime, as observed in many recent studies. Methods. We analyzed the integrated intensity ratio ($R$) and the wavelength-dependent ratio ($r(\Delta\lambda)$) in a small bifurcated eruption event observed by the Interface Region Imaging Spectrograph (IRIS). Results. Despite the relatively complex line profiles, most of the intensity ratio $R$ of \ion{Si}{IV} lines remained greater than 2 in the loops. The ratio $r(\Delta\lambda)$ varied in the line core and wings, changing distinctly from 2.0 to 3.3 along the wavelength. At certain positions, the \ion{Si}{IV} 1394 Å and 1403 Å lines exhibited different Doppler velocities. Conclusions. When diagnosing the spectra of small active region events, not only the impact of opacity but also the influence of resonance scattering should be considered. We propose that the ratio $r(\Delta\lambda)$ can serve as an indicator of the resonance scattering and opacity effect of the \ion{Si}{IV} line.

We investigate the effects of local features in the inflationary potential on the preheating dynamics after inflation. We show that a small feature in the potential can enhance the resonance and bring the radiation-like state equation during preheating despite the inflationary potential being a quadratic one. Such localized features may naturally arise due to various physical effects without altering the large-scale predictions of the original model for cosmic microwave background (CMB) observables. We demonstrate that these features effectively introduce localized higher-power terms in the potential, significantly influencing the preheating dynamics $\unicode{x2013}$ a phenomenon we term potential surge preheating. We outline the resulting modifications in energy distribution among different components. We further show that these small-scale features leave detectable imprints in the form of gravitational wave signals. These signals influence CMB measurements of the effective number of relativistic species, $N_{\mathrm{eff}}$, offering a way to reconstruct the shape of the inflaton potential at small scales. Finally, we argue that these modifications to the scalar potential provide a framework to explore preheating dynamics and the fragmentation of scalar fields using simple scalar potentials.

Red supergiants (RSGs) are essential to understanding the evolution and the contribution to the interstellar medium of massive stars. However, the number of identified RSGs within the Milky Way is still limited mainly due to the difficulty of measuring stellar extinction and distance. The release of approximately one million RVS spectra in Gaia DR3 presents new opportunity for identifying Galactic RSGs, because the equivalent width of the calcium triplet lines (EW(CaT)) in the spectra is an excellent indicator of stellar surface gravity. This work uses the RVS spectra with signal-to-noise ratio (SNR) greater than 100 to search for the Galactic RSGs. The dwarf stars and red giants are removed and the RSG candidates are selected by the location in the EW(CaT) vs. BP-RP diagram. The early-type RSG candidates (K0-M2) are then identified by BP-RP > 1.584 and EW(CaT) > 1.1 nm. To identify late-type RSG candidates (after M2), the criteria of the average equivalent widths of TiO in the XP spectra (EW(TiO)) > 10 nm, the color index K-W3 < 0.5 and the period-amplitude sequence from Gaia DR3 LPV catalog are further applied to reduce the contamination of late-type red giants and asymptotic giant branch stars. This method yields 30 early-type (K0-M2) and 6196 late-type (after M2) RSG candidates, which is a significant increase to the present Galactic RSG sample. The application of this approach to the spectra with SNR > 50 results in 48 early-type and 11,491 late-type RSG candidates. This preliminary analysis paves the way for more extensive research with Gaia DR4 when larger spectral datasets are expected to significantly enhance our understanding of Galactic RSG populations.

Extremely low-mass white dwarfs (ELM WDs) are helium-core white dwarfs with masses less than 0.3 $M_{\odot}$. Short-period ELM WD binaries that exhibit ellipsoidal variations may harbor heavier companions, either massive white dwarfs or millisecond pulsars (MSPs). In this study, we selected $\sim$ 12,000 ELM WDs or their candidates, and searched for ellipsoidal-like lightcurves with orbital periods shorter than one day, by using the public data from Zwicky Transient Facility. Finally 23 such systems were found, with 17 being newly discovered. We selected nine high-priority targets likely to evolve from the Roche-lobe overflow channel and estimated their companion masses from the extracted ellipsoidal variation amplitude. Among them, the four targets have companion masses exceeding 1 $M_{\odot}$. We performed a search for radio pulsations from six of these targets by using Five-hundred-meter Aperture Spherical radio Telescope. However, no convincing radio pulsed signals were found, resulting in upper limits for the radio flux at around 8 $\mu$Jy. Given the non-detection of radio pulsations from a total of 11 similar systems, the fraction of ellipsoidal ELM WDs around MSPs is estimated to be below 15$^{+6}_{-3}$%. We anticipate that multi-wavelength studies of more ellipsoidal-like ELM WDs will further constrain the fraction.

In this study, we investigated the anisotropy of diffusive Ultra-High Energy Cosmic Rays (UHECRs) by employing three cosmological models: two models from the $f(R, T)$ gravity theory and the other is the standard $\Lambda$CDM model. The primary objective of this work was to ascertain the role of the $f(R, T)$ gravity theory in comprehending the anisotropy of UHECRs without implicitly endorsing the conventional cosmology. We parameterized the magnetic field and the source distance in anisotropy calculations to align well with the observational data from the Pierre Auger Observatory for all the cosmological models. An uncertainty band is presented along with the $\chi^2$ test for all cosmological models to demonstrate the goodness of fitting. Our findings revealed that the amplitude of the anisotropy is highly sensitive to these cosmological models. Notably, the $f(R, T)$ models exhibited a lower amplitude of anisotropy (i.e., more isotropy), while the $\Lambda$CDM model predicted a comparatively higher amplitude at most of the energies considered.

Lorentz invariance violation is a feature of several quantum gravity models in which Lorentz symmetry is broken at high energies, leading to potential changes in particle behavior and interactions. In this study, we investigate vacuum Cherenkov radiation, a reaction in which an electron spontaneously emits a photon. This process, forbidden when considering unbroken Lorentz symmetry, is a phenomenological consequence of some quantum gravity models. We derive, for the first time, the spectra for the vacuum Cherenkov reaction, and confirm our results numerically. These results can be used to derive limits on Lorentz invariance violation.

Among the various implications of the spiral arms, it has been observed in the Milky Way disc that the distribution of radial actions from the Gaia Data Release 3 exhibits structures that may be related to the spiral arms. Our goal is to investigate the relationship between regions of low radial action identified in simulated discs and the location of the spiral arms, such as that suggested in Palicio et al. (2023) for the Galaxy. For a sample of 23 simulated spiral galaxies, we modelled the axisymmetric component of their gravitational potential to compute the radial action of their stellar particles using the Stäckel fudge. The spatial distribution of the radial action was then compared to the location of the spiral arms, identified as overdensities in the stellar surface density using a kernel density estimator. Our analysis reveals a strong correlation between the radial action distribution and the spiral arms in 18 of 23 simulated galaxies. However, notable discrepancies are observed in the remaining five, since they are profoundly out-of-equilibrium systems, such as galaxies influenced by external interactions or spiral arms still in the process of winding up. We have confirmed that, in general, there is a tendency of spatial correlation between spiral arms and stellar populations featuring low values of the radial action observed in Gaia DR3 data by Palicio et al. (2023). However, discrepancies between features in the radial action distribution and the spiral structure can be interpreted as signatures of recent disturbances, a scenario applicable to the Milky Way. Furthermore, populations at least as old as 3 Gyr trace the spiral arms with no significant misalignment across age bins. A linear relation between the maximum value of the radial action of the spiral arms and the vertical scale-length is found, which is also satisfied by the Milky Way.

We sought to determine which are the main hydrogenation paths of acetaldehyde (CH3CHO). As a partially unsaturated molecule, CH3CHO can have links with more hydrogenated species, like ethanol (C2H5OH) or with more unsaturated ones, like ketene (H2CCO). We used highly accurate quantum chemical calculations to determine the reaction rate constants for the CH3CHO + H/D reaction. Our theoretical results are confronted against our experiments on the hydrogenation and deuteration of CH3CHO ice. We find that acetaldehyde resists hydrogenation, with only a 10\% of conversion to products different than CH3CHO. This is due to a predominance of H-abstraction at the HCO moiety, with reaction rate constants up to four orders of magnitude higher than the next possible reaction channel, that is hydrogenation at the aldehydic carbon. The formed CH3CO radical experiences barrierless or nearly barrierless reactions in all possible reaction positions, reforming CH3CHO and creating a closed loop that protects the molecule against hydrogenation. We constrain the branching ratios for the second reaction from experiments. Our experiments agree with the calculations and from the combination of both we can explain the presence of H2CCO, CO, CH4, C2H5OH, H2CO or CH3OH as minor products at the end of the reaction. We provide recommendations for future modeling efforts. Our results show limited destruction of acetaldehyde, reinforcing the vision of this molecule as an abundant and resilient COM. From the experiments, we are not able to observe the reactive desorption of this molecule. Our results align with other modeling works, showing that the link between CH3CHO and C2H5OH is not direct. Finally, our results can explain the excess of CH3CDO found in prestellar cores.

The disk instability model suggests that the disk remains in a hot state during standstill. The thermal-tidal instability model posits that superoutbursts require a mass ratio of ( q = M2/M2 < 0.25-0.33). However, this paper presents evidence of superoutbursts and positive superhumps (PSHs) during a standstill of AT Cnc (q > 0.33 ), based on both space- and ground-based photometric data from sky surveys. Notably, the PSHs evolve gradually prior to the onset of a superoutburst, suggesting that an eccentric, prograde-precessing disk forms first, with the superoutburst occurring as the radius of the accretion disk continues to expand. These observations indicate that, during the standstill, the disk radius in AT Cnc not only surpasses the tidal truncation radius and the 3:1 resonance radius, but is still undergoing this http URL analysis further reveals that superoutbursts are triggered by oscillations, which appear to be a characteristic feature of special dwarf nova outbursts. Additionally, we observe that special outbursts develop from normal outbursts, implying that certain stages of the disk may not reach a fully hot state during the standstill. Furthermore, the PSH amplitude was found to be correlated with special outbursts, revealing that the localized thermal instability had a significant impact on the evolution of the eccentric disk. These findings provide the first detailed observational evidence of superoutbursts and PSHs occurring during standstill, offering important new insights into the classification of dwarf novae and the underlying mechanisms of outbursts.

As a major cooling line of interstellar gas, the far-infrared 158 {\mu}m line from singly ionised carbon [CII] is an important tracer of various components of the interstellar medium in galaxies across all spatial and morphological scales. Yet, there is still not a strong constraint on the origins of [CII] emission. In this work, we derive the resolved [CII] star formation rate relation and aim to unravel the complexity of the origin of [CII]. We used the Field-Imaging Far-Infrared Line Spectrometer on board the Stratospheric Observatory for Infrared Astronomy to map [CII] in three nearby star-forming galaxies at sub-kiloparsec scales, namely, NGC 3627, NGC 4321, and NGC 6946, and we compared these [CII] observations to the galactic properties derived from complementary data from the literature. We find that the relationship between the [CII] fine structure line and star formation rate shows variations between the galaxies as well as between different environments within each galaxy. Our results show that the use of [CII] as a tracer for star formation is much more tangled than has previously been suggested within the extragalactic literature, which typically focuses on small regions of galaxies and/or uses large-aperture sampling of many different physical environments. As found within resolved observations of the Milky Way, the picture obtained from [CII] observations is complicated by its local interstellar medium conditions. Future studies will require a larger sample and additional observational tracers, obtained on spatial scales within galaxies, in order to accurately disentangle the origin of [CII] and calibrate its use as a star formation tracer.

Over the past decades, advancements in observational cosmology have introduced us in an era of precision cosmology, dramatically enhancing our understanding of the Universe's history as well as bringing new tensions to light. Observations of the Cosmic Microwave Background, large-scale structure, and distant galaxies have provided unprecedented insights into the processes that shaped our Universe. This PhD thesis contributes to this research by exploring how these cosmic observables can be leveraged to constrain new physics beyond the standard $\Lambda$-Cold Dark Matter model.

With the discovery of gravitational waves (GWs), the disks of Active Galactic Nuclei (AGN) have emerged as an interesting environment for hosting a fraction of their sources. AGN disks are conducive to forming both long and short Gamma-Ray Bursts (GRBs), and their anticipated cosmological occurrence within these disks has potential to serve as an independent tool for probing and calibrating the population of stars and compact objects within them, and their contribution to the GW-detected population. In this study, we employ Monte Carlo methods in conjunction with models for GRB electromagnetic emission in extremely dense media to simulate the cosmological occurrence of both long and short GRBs within AGN disks, while also estimating their detectability across a range of wavelengths, from gamma-rays to radio frequencies. {We investigate two extreme scenarios: ``undiffused", in which the radiation escapes without significant scattering (i.e. if the progenitor has excavated a funnel within the disk), and ``diffused", in which the radiation is propagated through the high-density medium, potentially scattered and absorbed. {In the diffused case,} we find that the majority of detectable GRBs are likely to originate from relatively low redshifts, and from the outermost regions of large supermassive black hole (SMBH) masses, $\gtrsim 10^{7.5} \rm M_{\odot}$. In the undiffused case, we expect a similar trend, but with a considerable contribution from the intermediate regions of lower SMBH masses. Detectable emission is generally expected to be dominant in prompt $\gamma$-rays if diffusion is not dominant, and X-ray afterglow if diffusion is important; however, the nature of the dominant observable signal highly depends on the specific AGN disk model, hence making GRBs in AGN disks also potential probes of the structure of the disks themselves.

As part of the BISTRO survey, we present JCMT 850 $\mu$m polarimetric observations towards the Orion Integral-Shaped Filament (ISF) that covers three portions known as OMC-1, OMC-2, and OMC-3. The magnetic field threading the ISF seen in the JCMT POL-2 map appears as a tale of three: pinched for OMC-1, twisted for OMC-2, and nearly uniform for OMC-3. A multi-scale analysis shows that the magnetic field structure in OMC-3 is very consistent at all the scales, whereas the field structure in OMC-2 shows no correlation across different scales. In OMC-1, the field retains its mean orientation from large to small scales, but shows some deviations at small scales. Histograms of relative orientations between the magnetic field and filaments reveal a bimodal distribution for OMC-1, a relatively random distribution for OMC-2, and a distribution with a predominant peak at 90$^\circ$ for OMC-3. Furthermore, the magnetic fields in OMC-1 and OMC-3 both appear to be aligned perpendicular to the fibers, which are denser structures within the filament, but the field in OMC-2 is aligned along with the fibers. All these suggest that gravity, turbulence, and magnetic field are each playing a leading role in OMC-1, 2, and 3, respectively. While OMC-2 and 3 have almost the same gas mass, density, and non-thermal velocity dispersion, there are on average younger and fewer young stellar objects in OMC-3, providing evidence that a stronger magnetic field will induce slower and less efficient star formation in molecular clouds.

There is growing evidence for star formation inside outflows of active galactic nuclei (AGNs). The formed stars are injected into bulges and give rise to perturbation of bulges. In this paper, we investigate the issues of non-rotating, spherically symmetric bulges under the perturbation of fast, massive outflows with stars formed inside. We show that the potential perturbation of outflows together with injection and dynamical friction of these stars could drive bulge oscillations. Still, we find non-zero radial velocity of bulges will be driven by the episodic outflows of AGNs and after the AGN quenched, the radial velocity will tend to zero within a timescale $\sim\tau_{\rm AGN}$, which is the AGN's lifetime. For some typical values of bulges and AGNs, we find the expansion and contraction velocities are of a few $10\,\rm km\,s^{-1}$ for $10^{10}\,M_\odot$ bulges and mass outflowing rate $500\,M_\odot/\rm yr$, which would give observational signatures.

It has been found that some quasars are undergoing quasi-periodic variations (most of them with damped amplitudes) in optical bands from long-term monitoring campaigns, but how to explain the origin of such light curve variations still remains an open question. In this paper, we use the warped accretion disks model to explain the quasi-periodical variations. This model employs a free-bending wave traveling in an accretion disk which causes the orientation of the central part of the disk to oscillate from the line of sight, resulting in a quasi-periodical variation. We numerically solve the governing equation of warp propagation and calculate the simulated R-band light curves, finding that the periodical light curves generated by this model have damped amplitudes. To compare with observations, we select SDSSJ134820.42+194831.5 as a preliminary example from a sample of periodic quasar candidates by combining CRTS with other public survey data, and fitted its light curve with different observational angles. Our result gives a reduced $\chi^{2}\simeq 2.4$, implying that the model might give insights to future application of warped disk model.

CubeSats require robust OBDH solutions in harsh environments. The Demoiselle OBC, featuring a radiation-tolerant APSoC and layered FSW, supports reuse, in-orbit updates, and secure operations. To be validated through ITASAT2 and SelenITA, it ensures fault tolerance, flexibility, and compatibility with emerging technologies. This architecture establishes a foundation for long-lasting, scalable OBDH systems in future Brazilian CubeSat missions, ensuring long-term reliability and adaptability.

The main goal of the present research is to make a summary of analytical equations that can be found to calculate a swing-by maneuver in the three-dimensional space. Analytical equations based in the patched conics approximation are showed and they allow to calculate the variation in velocity, angular momentum, energy and inclination of the spacecraft that is involved in this maneuver. Based on that, it is possible to obtain expressions for particular cases, like the planar and the polar maneuver. The most important properties of this maneuver can demonstrated using thos equations, like: in the case of planar maneuver the variation in inclination can be only 180, 0, and -180 degrees; the variation in inclination is symmetric with respect to the out of plane angle; a passage by the poles changes only the y-component of the angular momentum, keeping the energy and the inclination of the trajectory unchanged. The results show several maneuvers.

The radio pulsar PSR J0742-2822 is known to exhibit rapid changes between different pulse profile states that correlate with changes in its spin-down rate. However, the connection between these variations and the glitch activity of the pulsar remains unclear. We aim to study the evolution of the pulse profile and spin-down rate of PSR J0742-2822 in the period MJD 58810-60149 (November 2019 to July 2023), which includes the glitch on MJD 59839 (September 2022). In particular, we look for pulse profile or spin-down changes associated with the 2022 glitch. We observed PSR J0742-2822 with high cadence from the Argentine Institute of Radio astronomy (IAR) between November 2019 and July 2023. We used standard timing tools to characterize the times of arrival of the pulses and study the pulsar rotation, and particularly, the oscillations of $\dot \nu$. We also study the evolution of the pulse profile. For both of them, we compare their behavior before and after the 2022 glitch. With respect to $\dot \nu$, we found oscillations diminished in amplitude after the glitch. We found four different components contributing to the pre-glitch $\dot \nu$ oscillations, and only one component after the glitch. About the emission, we found the pulse profile has two main peaks. We detected an increase in the $W_{50}$ of the total pulse profile of $\sim$12% after the glitch and we found the amplitude of the trailing peak increased with respect to the amplitude of the leading one after the glitch. We found significant changes in the pulse profile and the spin-down rate of PSR J0742-2822 after its 2022 glitch. These results suggest that there is a strong coupling between the internal superfluid of the neutron star and its magnetosphere, and that pulse profile changes may be led by this coupling instead of being led purely by magnetospheric effects.

A phase transition in the dark sector (DS) presents a promising explanation for the stochastic gravitational wave (GW) signals detected in recent observations by Pulsar Timing Arrays (PTAs). Instead of focusing solely on fitting data with phenomenological parameters, we systematically delineate simple, underlying dark sector (DS) models at the microscopic, Lagrangian level and uncover the conditions required to yield the GW spectrum observed by PTAs. We also illustrate the possibilities of the DS cosmology, which may include a dark matter candidate with a $\mathcal{O}$(MeV) mass.

The Laser Interferometer Space Antenna (LISA), a spaceborne gravitational wave (GW) detector set to launch in 2035, will observe several stochastic GW backgrounds in the mHz frequency band. At least one of these signals -- arising from the tens of millions of unresolved white dwarf binaries in the Milky Way -- is expected to be highly anisotropic on the sky. We evaluate the angular resolution of LISA and its ability to characterize anisotropic stochastic GW backgrounds (ASGWBs) using the Bayesian Spherical Harmonic formalism in the Bayesian LISA Inference Package (BLIP). We use \blip to simulate and analyze ASGWB signals in LISA across a large grid in total observing time, ASGWB amplitude, and angular size. We consider the ability of the \blip anisotropic search algorithm to both characterize single point sources and to separate two point sources on the sky, using a full-width half-max (FWHM) metric to measure the quality and spread of the recovered spatial distributions. We find that the number of spherical harmonic coefficients used in the anisotropic search model is the primary factor that limits the search's angular resolution. Notably, this trend continues until computational limitations become relevant around $\ell_{\mathrm{max}}=16$; this exceeds the maximum angular resolution achieved by other map-making techniques for LISA ASGWBs.

We have re-examined Mukhanov parametrization for inflationary equation of state $1+\omega=\frac{\beta}{({N}+1)^\alpha}$, in the light of Planck 2018 results and latest bound on tensor-to-scalar ratio employing Hamilton-Jacobi formalism. We have found the current observational values of scalar spectral index and tensor-to-scalar ratio can be used efficiently to constrain the model parameters. The recent findings of $r<0.032$ has been used to put an upper bound on one of the model parameter. Whereas the $1-\sigma$ bound of the scalar spectral index $0.9607\leq n_{_S}\leq 0.9691$ along with the upper bound of scalar-to-tensor ratio provided restriction on the other model parameter $1.50\leq\alpha\leq2.20$. This bound depends on the number of e-foldings still left before the end of inflation and whenever $1.50\leq\alpha\leq2.20$ we can find appropriate values of the other model parameter $\beta$ so that the observational predictions are in tune with the latest available inflationary observables.

Gravitational waveform (GW) models are a core ingredient for the analysis of compact binary mergers observed by current ground-based interferometers. We focus here on a specific class of such models known as PhenomX, which has gained popularity in recent years thanks to its computational efficiency. We introduce a new description of the ``twisting-up'' mapping underpinning the construction of precessing waveforms within this family. The new description is an adaptation to the frequency domain of a technique previously implemented in time-domain models, where the orbit-averaged post-Newtonian spin-precession dynamics is numerically solved on the fly. We also present an improved version of the gravitational-wave strain amplitudes approximating the signal in the co-precessing frame. We demonstrate that the new description yields improved matches against numerical relativity simulations, with only a modest computational overhead. We also show that the new model can be reliably employed in parameter estimation follow-ups of GW events, returning equivalent or more stringent measurements of the source properties compared to its predecessor.

We investigate a class of ion-scale magnetic solitary structures in the solar wind, characterized by distinct magnetic field enhancements and bipolar rotations over spatial scales of several proton inertial lengths. Previously tentatively identified as Alfvénic solitons, these structures are revisited using high-resolution data from the Solar Orbiter and Parker Solar Probe missions. Using a machine learning-based method, we identified nearly a thousand such structures, providing new insights into their evolution and physical properties. Statistical analysis shows that these structures are more abundant closer to the Sun, with occurrence rates peaking around 30-40 solar radii and declining at greater distances, suggesting that they decay. High-cadence measurements reveal that these structures are predominantly found in low-beta environments, with consistent fluctuations in density, velocity, and magnetic field. Magnetic field enhancements are often accompanied by plasma density drops, which, under near pressure balance, limit field increases. This leads to small fractional field enhancements near the Sun (approximately 0.01 at 20 solar radii), making detection challenging. Magnetic field variance analysis indicates that these structures are primarily oblique to the local magnetic field. Alfvénic velocity-magnetic field correlations suggest that most of these structures propagate sunward in the plasma frame, distinguishing them from typical solar wind fluctuations. We compare these findings with previous studies, discussing possible generation mechanisms and their implications for the turbulent cascade in the near-Sun Alfvénic solar wind. Further high-resolution observations and simulations are needed to fully understand their origins and impacts.

We explore the gravitational wave probes of a two-component dark matter framework, consisting of an $SU(2)_L$ triplet scalar and a Standard Model singlet fermion. The triplet scalar dark matter typically remains underabundant in the region below $\sim 1.9$ TeV, due to the strong $SU(2)_L$ gauge mediated interactions. We introduce a second dark matter component, an $SU(2)_L$ singlet vector-like Dirac fermion, to address this deficit in the dark matter relic abundance within a sub-TeV range. A key aspect of the proposed setup is the potential dark matter inter-conversion between the two components, which impacts the dark matter freeze-out dynamics and relic density of individual dark matter components. In such a scenario, we examine the properties of electroweak phase transition and identify the regions of parameter space that exhibit strong first-order phase transition. We estimate the resulting gravitational wave spectrum and its detectability, which could be probed through the conventional power-law-integrated sensitivity limits and the recently proposed peak-integrated sensitivity curves. Our analysis reveals that a novel region of the model's parameter space, compatible with dark matter observables, can generate a detectable gravitational wave spectrum, observable by upcoming space-based gravitational wave detectors such as LISA, BBO, DECIGO, and DECIGOcorr, while also offering complementary detection prospects in the dark matter and collider experiments.

Gravitational wave observations by ground based detectors such as LIGO and Virgo have transformed astrophysics, enabling the study of compact binary systems and their mergers. However, transient noise artifacts, or glitches, pose a significant challenge, often obscuring or mimicking signals and complicating their analysis. In this work, we extend the Attention-boosted Waveform Reconstruction network to address glitch mitigation, demonstrating its robustness in reconstructing waveforms in the presence of real glitches from the third observing run of LIGO. Without requiring explicit training on glitches, AWaRe accurately isolates gravitational wave signals from data contaminated by glitches spanning a wide range of amplitudes and morphologies. We evaluate this capability by investigating the events GW191109 and GW200129, which exhibit strong evidence of anti-aligned spins and spin precession respectively, but may be adversely affected by data quality issues. We find that, regardless of the potential presence of glitches in the data, AWaRe reconstructs both waveforms with high accuracy. Additionally, we perform a systematic study of the performance of AWaRe on a simulated catalog of injected waveforms in real LIGO glitches and obtain reliable reconstructions of the waveforms. By subtracting the AWaRe reconstructions from the data, we show that the resulting residuals closely align with the background noise that the waveforms were injected in. The robustness of AWaRe in mitigating glitches, despite being trained exclusively on GW signals and not explicitly on glitches, highlights its potential as a powerful tool for improving the reliability of searches and characterizing noise artifacts.

We explore the production of gravitational waves (GW) resulting from a first-order phase transition (FOPT) in a non-minimally coupled `Dark Higgs Inflation' model. Utilizing a dark sector scalar field as the inflaton, we demonstrate how inflationary dynamics naturally set the stage for observable FOPT. These transitions, influenced by thermal and quantum effects, generate GW spectra potentially detectable by observatories such as LISA, DECIGO, the Cosmic Explorer and the Einstein Telescope. Our study highlights the inflaton's dual role in cosmic inflation and early Universe phase transitions, presenting a unified framework to probe physics beyond the Standard Model through gravitational wave astronomy.

We analyze direct numerical simulations of large-scale dynamos in inhomogeneous nonhelically-driven rotating turbulence with and without shear. The forcing is modulated so that the turbulent intensity peaks in the middle of the computational box and drops to nearly zero at the two ends above and below the midplane. A large-scale dynamo is driven by an $\alpha$ effect of opposite signs between the two hemispheres. In the presence of shear, the hemispheric magnetic helicity flux from small-scale fields becomes important and can even overcompensate for the magnetic helicity transferred by the $\alpha$ effect between large and small scales. This effect has not previously been observed in non-shearing simulations. Our numerical simulations show that the hemispheric magnetic helicity fluxes are nearly independent of the magnetic Reynolds number, but those between large and small scales, and the consequent dynamo effect, are still found to decrease with increasing Reynolds number -- just like in nonshearing dynamos. However, in contrast to nonshearing dynamos, where the generated mean magnetic field declines with increasing magnetic Reynolds number, it is now found to remain independent of it. This suggests that catastrophic dynamo quenching is alleviated by the shear-induced hemispheric small-scale magnetic helicity fluxes that can even overcompensate the fluxes between large and small scales and thereby cause resistive contributions.

The inert doublet model (IDM), a minimal extension of the Standard Model (SM), provides a scalar dark matter (DM) candidate that belongs to the additional Higgs doublet. The model faces challenges in achieving the correct relic abundance for compressed spectra and DM masses in the high-mass range. In this work we introduce a $Z_2$-odd singlet vector-like quark (VLQ) into the IDM framework that helps us alleviate these issues and provide new channels of contributions to the relic abundance. The VLQ not only enhances the DM relic abundance for masses above $~550$ GeV but also eases constraints from direct detection experiments by enabling smaller couplings between the inert scalars and the SM Higgs. We analyze the impact of the VLQ on DM phenomenology, including relic density, direct and indirect detection constraints. The results demonstrate that the extended IDM framework not only resolves existing limitations in the compressed spectrum but also offers exciting prospects for detection in current and future collider experiments.

We study symmetry-breaking inflation within the framework of metric-affine gravity. By introducing a non-minimal coupling, $\beta(\phi)\tilde{\cal R}$, between the Holst invariant and the inflaton, both small-field and large-field inflationary predictions can be brought into agreement with the latest observational constraints. Remarkably, even for sub-Planckian vacuum expectation values, appropriately chosen values of $\beta(\phi)$ enable viable inflation, a scenario previously considered unattainable.

We study dark gauge-mediated supersymmetry breaking (dark GMSB) in a theory with a new unbroken $U(1)_{D}$ local symmetry and massless dark photon. Messenger fields charged under both Standard Model and dark gauge symmetries produce new soft supersymmetry-breaking terms due to gauge kinetic mixing between $U(1)_Y$ hypercharge and $U(1)_D$. We show that large kinetic mixing induces significant distortions to the superpartner spectra relative to conventional GMSB. Notably, shifts in the Higgs soft masses impact the conditions for electroweak symmetry breaking, lowering the $\mu$ parameter and yielding a relatively light Higgsino that may be accessible at the LHC. Furthermore, for very simple messenger representations, a very light bino-dark photino mixed state is present in the spectrum, which may be probed through exotic Higgs boson decays at future Higgs factories. We also examine the cosmological and phenomenological consequences of the messengers, the lightest of which is absolutely stable and carries fractional electric charge.

The QCD axion is widely studied as a dark matter (DM) candidate and as a solution to the strong CP problem of the Standard Model. In conventional field-theoretic models, a much larger mass scale than the electroweak (EW) scale is typically introduced to spontaneously break Peccei-Quinn (PQ) symmetry with a large enough axion decay constant, $f_a$, thereby avoiding constraints from star cooling. In this paper, I propose an alternative approach to achieving the large decay constant: a PQ scalar field with a large wave function renormalization constant, analogous to a feebly coupled gauge theory. Other dimensionless parameters are ${O}(1)$ in the unit of the EW scale for the naturalness. This framework predicts a light PQ Higgs boson with a mass $\sim (\mathrm{EW~scale})^2 / f_a$. Exotic particles associated with the PQ anomaly are expected to have masses around the EW scale. The proposed model alleviates both the PQ quality and EW scale fine-tuning problems and introduces interesting axion-PQ Higgs cosmologies, encompassing: slim axion DM from a fat string network, heavy axion DM from PQ Higgs condensate fragmentation, PQ Higgs DM, and axion-PQ Higgs co-DM scenarios. Potential experimental signatures are explored, including fifth-force tests, DM detections, accelerator searches, and gravitational wave observations by employing lattice simulation. Possible extensions of the scenario are also discussed.