Papers for Monday, May 08 2023

Gravitational wave recoil kicks from merging supermassive black hole binaries can have a profound effect on the surrounding stellar population. In this work, we study the dynamic and kinematic properties of nuclear star clusters following a recoil kick. We show that these post-kick structures present unique signatures that can provide key insight to observational searches for recoiling supermassive black holes. In Akiba & Madigan (2021), we showed that an in-plane recoil kick turns a circular disk into a lopsided, eccentric disk such as the one we observe in the Andromeda nucleus. Building on this work, here we explore many recoil kick angles as well as initial stellar configurations. For a circular disk of stars, an in-plane kick causes strong apsidal alignment with a significant fraction of the disk becoming retrograde at large radii. If initial orbits are highly eccentric, an in-plane kick forms a bar-like structure made up of two anti-aligned lopsided disks. An out-of-plane kick causes clustering in the argument of periapsis, $\omega$, regardless of the initial eccentricity distribution. Initially isotropic configurations form anisotropies in the form of a torus of eccentric orbits oriented perpendicular to the recoil kick. Post-kick surface density and velocity maps are presented in each case to highlight the distinct, observable structures of these systems.

We study the structure and total energy of a strange quark star (SQS) endowed with a strong magnetic field with different rotational frequencies. The MIT bag model is used, with the density-dependent bag constant for the equation of state (EOS). The EOS is computed considering the Landau quantization effect regarding the strong magnetic fields (up to $5\times10^{17}$ G) in the interior of the strange quark star. Using the LORENE library, we calculate the structural parameters of SQS for different setups of magnetic field strengths and rotational frequencies. In each setup, we perform calculations for $51$ stellar configurations, with specified central enthalpy values. We investigate the configurations with the maximum gravitational mass of SQS in each setup. Our models of SQSs are compared in the maximum gravitational mass, binding energy, compactness, and deformation of the star. We show that the gravitational mass might exceed $2.3 M_\odot$ in some models, which is comparable with the mass of the recently detected ``black widow'' pulsar \emph{PSR J0952-0607} and the mass of \emph{GW190814} detected by the LIGO/Virgo collaboration. The deformation and maximum gravitational mass of SQS can be characterized by simple functions that have been fitted to account for variations in both magnetic field strength and frequency. Rapidly rotating strange stars have a minimum gravitational mass given by the equatorial mass-shedding limit.

In preceding papers, Project Lyra has covered many possible trajectory options available to a spacecraft bound for 1I/Oumuamua, including Solar Oberth manoeuvres, Passive Jupiter encounters, Jupiter Oberths, Double Jupiter Gravitational Assists, etc. Because feasibility was the key driver for this analysis, the important question of which launcher to exploit was largely skirted in favour of adopting the most powerful options as being sufficient, though these launchers are clearly not necessary, there being alternative less capable candidates which could be utilised instead. In this paper the various launch options available to Project Lyra are addressed to allow a general overview of their capabilities. It is found that the SpaceX Super-Heavy Starship would be a game-changer for Project Lyra, especially in the context of refuelling in LEO, and furthermore a SpaceX Falcon Heavy Expendable could also be utilised. Other launchers are considered, including Ariane 6 and the future Chinese Long March 9. The importance of the V infinity Leveraging Manoeuvre (VILM) in permitting less capable launchers to nevertheless deliver a payload to Oumuamua is elaborated

We implement a novel formalism to constrain primordial non-Gaussianity of the local type from the large-scale modulation of the small-scale power spectrum. Our approach combines information about primordial non-Gaussianity contained in the squeezed bispectrum and the collapsed trispectrum of large-scale structure together in a computationally amenable and consistent way, while avoiding the need to model complicated covariances of higher $N$-point functions. This work generalizes our recent work, which used a neural network estimate of local power, to the more conventional local power spectrum statistics, and explores using both matter field and halo catalogues from the Quijote simulations. We find that higher $N$-point functions of the matter field can provide strong constraints on $f_{NL}$, but higher $N$-point functions of the halo field, at the halo density of Quijote, only marginally improve constraints from the two-point function.

We present very early photometric and spectroscopic observations of the Type Ia supernova (SN) 2023bee, starting about 8 hours after the explosion, which reveal a strong excess in the optical and nearest UV (U and UVW1) bands during the first several days of explosion. This data set allows us to probe the nature of the binary companion of the exploding white dwarf and the conditions leading to its ignition. We find a good match to the Kasen model in which a main-sequence companion star stings the ejecta with a shock as they buzz past. Models of double detonations, shells of radioactive nickel near the surface, interaction with circumstellar material, and pulsational-delayed detonations do not provide good matches to our light curves. We also observe signatures of unburned material, in the form of carbon absorption, in our earliest spectra. Our radio non-detections place a limit on the mass-loss rate from the putative companion that rules out a red giant but allows a main-sequence star. We discuss our results in the context of other similar Type Ia SNe in the literature.

The study of dynamically cold stellar streams reveals information about the gravitational potential where they reside and provides important constraints on dark matter properties. However, their intrinsic faintness makes detection beyond Local environments highly challenging. Here we report the detection of an extremely faint stellar stream (mu_g,max = 29.5 mag arcsec-2) with an extraordinarily coherent and thin morphology in the Coma Galaxy Cluster. This Giant Coma Stream spans 510 kpc in length and appears as a free-floating structure located at a projected distance of 0.8 Mpc from the center of Coma. We do not identify any potential galaxy remnant or core, and the stream structure appears featureless in our data. We interpret the Giant Coma Stream as being a recently accreted, tidally disrupting dwarf of M* ~ 10^8 Msun, and report a case with similar characteristics within the Illustris-TNG50 simulation. Our work shows the presence of free-floating, extremely faint and thin stellar streams in galaxy clusters, widening the environmental context for their promising future applications in the study of dark matter properties.

Observations of the cosmic 21-cm power spectrum (PS) are starting to enable precision Bayesian inference of galaxy properties and physical cosmology, during the first billion years of our Universe. Here we investigate the impact of common approximations about the likelihood used in such inferences, including: (i) assuming a Gaussian functional form; (ii) estimating the mean from a single realization; and (iii) estimating the (co)variance at a single point in parameter space. We compare "classical" inference that uses an explicit likelihood with simulation based inference (SBI) that estimates the likelihood from a training set. Our forward-models include: (i) realizations of the cosmic 21-cm signal computed with 21cmFAST by varying UV and X-ray galaxy parameters together with the initial conditions; (ii) realizations of the telescope noise corresponding to a 1000 h integration with SKA1-Low; (iii) the excision of Fourier modes corresponding to a foreground-dominated, horizon "wedge". We find that the 1D PS likelihood is well described by a Gaussian accounting for covariances between wavemodes and redshift bins (higher order correlations are small). However, common approaches of estimating the forward-modeled mean and (co)variance from a random realization or at a single point in parameter space result in biased and over-constrained posteriors. Our best results come from using SBI to fit a non-Gaussian likelihood with a Gaussian mixture neural density estimator. Such SBI can be performed with up to an order of magnitude fewer simulations than classical, explicit likelihood inference. Thus SBI provides accurate posteriors at a comparably low computational cost.

We use explainable neural networks to connect the evolutionary history of dark matter halos with their density profiles. The network captures independent factors of variation in the density profiles within a low-dimensional representation, which we physically interpret using mutual information. Without any prior knowledge of the halos' evolution, the network recovers the known relation between the early time assembly and the inner profile, and discovers that the profile beyond the virial radius is described by a single parameter capturing the most recent mass accretion rate. The results illustrate the potential for machine-assisted scientific discovery in complicated astrophysical datasets.

With its cored surface brightness profile, the elliptical galaxy NGC 5419 appears as a typical high-mass early-type galaxy (ETG). However, the galaxy hosts two distinct nuclei in its center. We use high-signal MUSE (Multi-Unit Spectroscopic Explorer) spectral observations and novel triaxial dynamical orbit models to reveal a surprisingly isotropic central orbit distribution in NGC 5419. Recent collisionless simulations of merging massive ETGs suggest a two-phase core formation model, in which the low-density stellar core forms rapidly by supermassive black holes (SMBHs) sinking into the center due to dynamical friction. Only afterwards the SMBHs form a hard binary and the black hole scouring process slowly changes the central orbit distribution from isotropic to tangential. The observed cored density profile, the double nucleus and the isotropic center of NGC 5419 together thus point to an intermediate evolutionary state where the first phase of the core formation has taken place, yet the scouring process is only beginning. This implies that the double nucleus is a SMBH binary. Our triaxial dynamical models indicate a total mass of the two SMBHs in the center of NGC 5419 of MBH = (1.0 +/- 0.08) 10^10 Msol. Moreover, we find that NGC 5419's complex kinematically distinct core (KDC) can be explained by a coherent flip of the orbital rotation direction of stars on tube orbits at ~3kpc distance from the galaxy center together with projection effects. This is also in agreement with merger simulations hosting SMBHs in the same mass regime.

We present the results from our ongoing spectroscopic survey targeting low mass white dwarf binaries, focusing on the southern sky. We used a Gaia DR2 and eDR3 based selection and identified 28 new binaries, including 19 new extremely low mass white dwarfs, one short period, likely eclipsing, DABZ, and two potential LISA binaries. We present orbital and atmospheric parameters for each new binary based on our spectroscopic follow-up. Four of our new binaries show periodic photometric variability in the TESS 2-minute cadence data, including one new eclipsing double-lined spectroscopic binary. Three others show periodic photometric variability in ZTF, including one new eclipsing binary. We provide estimates for the inclinations and scaled component radii for these ZTF variables, based on light curve modeling to our high-speed photometric follow-up observations. Our observations have increased the sample of ELM Survey binaries identified in the southern sky to 41, an increase of 64%. Future time domain surveys, such as BlackGEM and the Vera C. Rubin Observatory Legacy Survey of Space and Time, will efficiently identify the photometric variables in the southern sky and significantly increase the population of southern sky low mass white dwarf binaries, leading to a more complete all-sky population of these systems.

Z Canis Majoris is a fascinating early-type binary with a Herbig Be primary and a FU Orionis-type secondary. Both of the stars exhibit sub-arcsecond jet-like ejecta. In addition, the primary is associated with the extended jet as well as with the large-scale outflow. In this study, we investigate further the nature of the large-scale outflow, which has not been studied since its discovery almost three and a half decades ago. We present proper motion measurements of individual features of the large-scale outflow and determine their kinematical ages. Furthermore, with our newly acquired deep images, we have discovered additional faint arc-shaped features that can be associated with the central binary.

We investigate the statistical properties and the origin of the scatter within the spatially resolved surface brightness profiles of the CHEX-MATE sample, formed by 118 galaxy clusters selected via the SZ effect. These objects have been drawn from the Planck SZ catalogue and cover a wide range of masses, M$_{500}=[2-15] \times 10^{14} $M$_{\odot}$, and redshift, z=[0.05,0.6]. We derived the surface brightness and emission measure profiles and determined the statistical properties of the full sample. We found that there is a critical scale, R$\sim 0.4 R_{500}$, within which morphologically relaxed and disturbed object profiles diverge. The median of each sub-sample differs by a factor of $\sim 10$ at $0.05\,R_{500}$. There are no significant differences between mass- and redshift-selected sub-samples once proper scaling is applied. We compare CHEX-MATE with a sample of 115 clusters drawn from the The Three Hundred suite of cosmological simulations. We found that simulated emission measure profiles are systematically steeper than those of observations. For the first time, the simulations were used to break down the components causing the scatter between the profiles. We investigated the behaviour of the scatter due to object-by-object variation. We found that the high scatter, approximately 110%, at $R<0.4R_{500}$ is due to a genuine difference between the distribution of the gas in the core. The intermediate scale, $R_{500} =[0.4-0.8]$, is characterised by the minimum value of the scatter on the order of 0.56, indicating a region where cluster profiles are the closest to the self-similar regime. Larger scales are characterised by increasing scatter due to the complex spatial distribution of the gas. Also for the first time, we verify that the scatter due to projection effects is smaller than the scatter due to genuine object-by-object variation in all the considered scales. [abridged]

We perform extensive spectroscopy of the supernova remnant N63A in the Large Magellanic Cloud, using $\sim 43$ ks {\it Chandra} archival data. By analysing the spectra of the entire remnant, we determine the abundance distributions for O, Ne, Mg, Si, and Fe. We detect evidence of enhanced O and possibly Ne and Mg in some of the central regions which might indicate an asymmetric distribution of the ejecta. The average O/Ne, O/Mg, and Ne/Mg abundance ratios of the ejecta are in plausible agreement with the nucleosynthesis products from the explosion of a $\sim40$ $M_{\odot}$ progenitor. We estimate an upper limit on the Sedov age of $\sim 5,400\pm200$ yr and explosion energy of $\sim 8.9\pm 1.6\times 10^{51}$ erg for N63A. We discuss the implications of our results for the morphological structure of the remnant, its circumstellar medium and the nature of the progenitor star.

Short and ultra-short planets are a peculiar type of exoplanets with periods as short as a few days or less. Although it is challenging to detect them, already several are observed with many additional candidates. If these planets have formation pathways to their longer period counterparts, they are predicted to reside in multi-planet systems. Thus, gravitational perturbation from potential planetary neighbors may affect their orbital configuration. However, due to their close proximity to their host star, they are also subjects to general relativity precession and torques from the stellar spin quadrupole moment ($J_2$). Here we show that an evolving $J_2$ due to magnetic braking, affects the magnitude and location of secular resonances of the short period planet in a multi-planet system. Thus, driving the short period planet into and out of a secular resonance, exciting the planet's eccentricity and inclination. The high inclination can hinder transit observation, and, in some cases, the high eccentricity may result in an unstable configuration. We propose that evolving $J_2$ in a multi-planet system can be critical in understanding the detectability and stability of short-period planets.

Water and the hydroxyl radical are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are highly excited have been observed in both H$_2$O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution with the ALMA Array. Spectra and maps of H$_2$O and OH observed with an angular resolution of 20 to $\sim$200 mas were obtained at two epochs with the ALMA array. Observations with the Compact Array were also used to check for time variability of water transitions. Radiative transfer models of water were revisited to characterize masing conditions and up-to-date chemical models were used for comparison with our observations. Ten rotational transitions of H$_2$O with energies up to 9000 K were observed in various vibrational states. All but one are new detections in space, and from these we have derived accurate rest frequencies. Hyperfine split $\Lambda$-doubling transitions in v = 0, J = 27/2 and 29/2 levels of the $^2\Pi_{3/2}$ state and, $J = 33/2$ and 35/2 of the $^2\Pi_{1/2}$ state of OH with excitation energies up to 8900 K were also observed. Four of these transitions are new detections in space. Combining our measurements with earlier observations of OH, the v = 0 and v = 1 $\Lambda$-doubling frequencies have been improved. Our H$_2$O maps show compact emission and extensions up to twelve stellar radii or more. The 268.149 GHz emission line of water in the v$_2$ = 2 state is time variable, tends to be masing with dominant radiative pumping, and is widely excited. The widespread but weaker 262.898 GHz water line in v$_2$ = 1 also shows signs of maser emission. Emission and absorption of both H$_2$O and OH reveal an infall of matter and complex kinematics influenced by binarity. From our observed column densities, we derived OH/H$_2$O abundance ratios in a few stars.

Solar coronal mass ejections are the most energetic events in the Solar System. In their standard formation model, a magnetic flux rope builds up into a coronal mass ejection through magnetic reconnection that continually converts overlying, untwisted magnetic flux into twisted flux enveloping the pre-existing rope. However, only a minority of coronal mass ejections carry a coherent magnetic flux rope as their core structure, which casts doubt on the universality of this orderly wrapping process. Here we provide observational evidence of a different formation and eruption mechanism of a magnetic flux rope from an S-shaped thread, where its magnetic flux is fully replaced via flare reconnections. One of the footpoints of the sigmoidal feature slipped and expanded during the formation, and then moved to a completely new place, associated with the highly dynamical evolution of flare ribbons and a twofold increase in magnetic flux through the footpoint, during the eruption. Such a configuration is not predicted by standard formation models or numerical simulations and highlights the three-dimensional nature of magnetic reconnections between the flux rope and the surrounding magnetic field.

The asteroseismic scaling relation, dnu~rho^{0.5}, linking a star's large frequency separation, dnu, and its mean density, rho, is not exact. Yet, it provides a very useful way to obtain fundamental stellar properties. Common ways to make the relation more accurate is to apply correction factors to it. Because the corrections depend on stellar properties, such as mass, Teff, and metallicity, it is customary to interpolate these properties over stellar model grids that include both dnu, measured from adiabatic frequencies of the models, and the models' stellar density; hence linking both sides of the scaling relation. A grid and interpolation tool widely used for this purpose, known as Asfgrid, was published by Sharma & Stello 2016. Here, we present a significant extension of Asfgrid to cover higher- and lower-mass stars and to increase the density of grid points, especially in the low-metallicity regime.

One of the most protracted problems in astronomy has been understanding the evolution of galaxy morphology. Much discussion has surrounded how lenticular galaxies may form a bridging population between elliptical and spiral galaxies. However, with recourse to a galaxy's central black hole mass, accretion-built spiral galaxies have emerged as the bridging population between low-mass lenticular galaxies and the dusty merger-built lenticular galaxies contiguous with elliptical galaxies and `brightest cluster galaxies' in the black hole/galaxy mass diagram. Spiral galaxies, including the Milky Way, appear built from gas accretion and minor mergers onto what were initially lenticular galaxies. These connections are expressed as a new morphology sequence, dubbed the `Triangal', which subsumes elements of the Hubble sequence and the van den Bergh trident and reveals the bridging nature of the often overlooked ellicular galaxies. Furthermore, a quadratic black hole/galaxy mass relation is found to describe ordinary elliptical galaxies. The relation is roughly parallel to the quadratic-like relations observed for the central spheroidal component of spiral galaxies, dust-rich lenticular galaxies, and old dust-poor lenticular galaxies. The brightest cluster galaxies are offset according to expectations from an additional major merger. The findings have implications for feedback from active galactic nuclei, mapping morphology into simulations, and predicting gravitational wave signals from colliding supermassive black holes. A new galaxy speciation model is presented. It disfavours the `monolithic collapse' scenario for spiral, dusty lenticular, and elliptical galaxies. It reveals substantial orbital angular momentum in the Universe's first galaxies and unites dwarf and ordinary `early-type' galaxies.

High synchrotron energy peaked blazar 1ES 1959+650 is studied with Swift and XMM-Newton satellite in total 127 observations during the period June 2018-December 2020. We extensively studied its flux and spectral variability on intra-day and long-term timescales. Discrete correlation function analysis between soft and hard X-ray bands indicates soft as well as hard lags. The results are used to constrain the magnetic field of the emitting region which is found to be 0.64 (0.05) Gauss. On long-term timescales, distribution of fluxes shows lognormality behaviour which could be attributed to minijets-in-a-jet model or might be due to the propagation of relativistic shocks down the jet. The spectral energy distribution around the synchrotron peak is well described by the log parabola model. Spectral parameters like peak energy E$_{p}$, curvature $\beta$ and the peak luminosity L$_{p}$ are derived from spectral analysis. Their correlations are studied to constrain the acceleration processes of the emitting particles. E$_{p}$ shows strong correlation with L$_{p}$ during the high state of the source which indicates spectral changes might be caused by the variations of the average electron energy. Low values of curvature parameter $\beta$ and a weak correlation between E$_{p}$ and ${\beta}$ indicates co-existence of stochastic/statistical acceleration of electrons in the emitting region. Implications of other results are also discussed.

Young planets (< 1 Gyr) are helpful for studying the physical processes occurring at the early stage of planet evolution. TOI-251 b is a recently discovered sub-Neptune orbiting a young G dwarf, which has an imprecise age estimation of 40-320 Myr. We select TOI-251 sibling candidates based on kinematics and spatial proximity to TOI-251, and further use the color-magnitude diagram (CMD) to refine the list and to compare to multiple open clusters. We report stellar rotational period for 321 sibling candidates in a 50 pc radius around TOI-251 by analyzing their stellar light curves, and find a color - rotational period sequence that lie in between the Group X (300 Myr) and Pleiades (120 Myr) members, suggesting an age ~ 200 Myr. A quantitative age analysis by using gyrochronology relations give 204 $\pm$ 45 Myr, consistent with the average Li-age of selected siblings (238 $\pm$ 38 Myr) and the Gaia variability age (193$^{102}_{-54}$ Myr). The detection fraction of comoving candidates that have short rotational period is 68.1%, much higher than the typical value in the field (14% - 16% from Kepler). The overdensity of young stars and consistency in age of stellar siblings suggest a potential young association candidate in the Pheonix-Grus constellation. Though TOI-251 b has a radius larger than most of its field-age counterparts, we are uncertain whether TOI-251 is inflated due to a lack of knowledge on the planet's mass.

Andromeda (M 31) is the nearest giant spiral galaxy to our Milky Way, and over the past few decades, has been dubbed the most massive member of the Local Group. I explore the evolution of the measured mass of M 31 over the past ~80 years, reviewing the different observational and modelling techniques that have developed over time to measure its mass. I discuss the best present-day constraints of the mass of M 31 and the consistency of different techniques.

We present a catalog of distances for 19544 K giants drawn from LAMOST DR8. Most of them are located in the halo of the Milky Way up to ~120~kpc. There are 15% K giants without SDSS photometry, for which we supplements with Pan-STARRS1 (PS1) photometry calibrated to SDSS photometric system. The possible contamination of the red clumps/horizontal branch are removed according to metallicities and colors before the distance determination. Combining the LAMOST spectroscopic metallicities with the SDSS/PS1 photometry, we estimate the absolute magnitudes in SDSS $r-$band, the distance moduli, and the corresponding uncertainties through an Bayesian approach devised by Xue et al. (2014) for the SEGUE halo K-giants. The typical distance precision is about 11%. The stars in the catalog lie in a region of 4-126 kpc from the Galactic center, of which with 6, 320 stars beyond 20 kpc and 273 stars beyond 50 kpc, forming the largest spectroscopic sample of distant tracers in the Milky Way halo so far.

We model the Parker instability in vertically stratified isothermal gas using non-ideal MHD three-dimensional simulations. Rotation, especially differential, more strongly and diversely affects the nonlinear state than the linear stage (where we confirm the most important conclusions of analytical models), and stronger than any linear analyses predict. Steady state magnetic fields are stronger and cosmic ray energy density higher than in comparable nonrotating systems. Transient gas outflows induced by the nonlinear instability persist longer, of order 2 Gyr, with rotation. Stratification combined with (differential) rotation drives helical flows, leading to mean-field dynamo. Consequently, the nonlinear state becomes oscillatory (while both the linear instability and the dynamo are non-oscillatory). The horizontal magnetic field near the midplane reverses its direction propagating to higher altitudes as the reversed field spreads buoyantly. The spatial pattern of the large-scale magnetic field may explain the alternating magnetic field directions in the halo of the edge-on galaxy NGC 4631. Our model is unique in producing a large-scale magnetic structure similar to such observation. Furthermore, our simulations show that the mean kinetic helicity of the magnetically driven flows has the sign opposite to that in the conventional non-magnetic flows. This has profound consequences for the nature of the dynamo action and large-scale magnetic field structure in the coronae of spiral galaxies which remain to be systematically explored and understood. We show that the energy density of cosmic rays and magnetic field strength are not correlated at scales of order a kiloparsec.

We summarize in this paper the spectro-polarimetric methods used at the Pic du Midi Turret Dome in spectroscopic or imagery mode. The polarimeters and spectrograph allow the cartography of solar magnetic fields at high spatial resolution through the Zeeman effect or measurements of the unresolved turbulent magnetic fields in the quiet Sun through the Hanle effect. We describe in this paper the optical capabilities of the successive versions of the polarimeters operating since 2003, and we present new results of magnetic field analysis with the CaII K 3933.7 {\AA} spectral line.

The effects of the interaction between Type Ia supernova ejecta and their circumstellar wind on the photometric properties of Type Ia supernovae are investigated. We assume that a hydrogen-rich, dense, and extended circumstellar matter (CSM) is formed by the steady mass loss of their progenitor systems. The CSM density is assumed to be proportional to r^{-2}. When the mass-loss rate is above 1e-4 Msun/yr with a wind velocity of 100 km/s, CSM interaction results in an early flux excess in optical light-curves within 4 days of explosion. In these cases, the optical colour quickly evolves to the blue. The ultraviolet flux below 3000 A is found to have a persistent flux excess compared to Type Ia supernovae as long as CSM interaction continues. Type Ia supernovae with progenitor mass-loss rates between 1e-4 and 1e-3 Msun/yr may not have a CSM that is dense enough to affect spectra to make them Type Ia-CSM, but they may still result in Type Ia supernovae with an early optical flux excess. Because they have a persistent ultraviolet flux excess, ultraviolet light curves around the luminosity peak would be significantly different from those with a low-density CSM.

Disc accretion rate onto low mass protostar FU Ori suddenly increased hundreds of times 85 years ago and remains elevated to this day. We show that the sum of historic and recent observations challenges existing FU Ori models. We build a theory of a new process, Extreme Evaporation (EE) of young gas giant planets in discs with midplane temperatures exceeding 30, 000 K. Such temperatures are reached in the inner 0.1 AU during thermal instability bursts. In our 1D time-dependent code the disc and an embedded planet interact through gravity, heat, and mass exchange. We use disc viscosity constrained by simulations and observations of dwarf novae instabilities, and we constrain planet properties with a stellar evolution code. We show that dusty gas giants born in the outer self-gravitating disc reach the innermost disc in a $\sim$ 10,000 years with radius of $\sim 10 R_J$. We show that their EE rates are $\sim 10^{-5}$ Msun/yr; if this exceeds the background disc accretion activity then the system enters a planet-sourced mode. Like a stellar secondary in mass-transferring binaries, the planet becomes the dominant source of matter for the star, albeit for $\sim$ O(100) years. We find that a $\sim$ 6 Jupiter mass planet evaporating in a disc fed at a time-averaged rate of $\sim 10^{-6}$ Msun/yr appears to explain all that we currently know about FU Ori accretion outburst. More massive planets and/or planets in older less massive discs do not experience EE process. Future FUOR modelling may constrain planet internal structure and evolution of the earliest discs.

It has been suggested that there is evidence for cosmological coupling of black holes (BHs) with an index of $k\approx 3$ and hence the BHs serve as the astrophysical source of the dark energy. The data sample however is limited for the redshifts $\leq 2.5$. Recently, the James Webb Space Telescope (JWST) has detected more than 180 high-redshift Active Galactic Nuclei (AGNs) and quasars. Among the JWST NIRSpec/NIRCam resolved AGNs, three are identified in early-type host galaxies with a redshift $z\sim 4.5-7$. Their $M_{\star}$ and $M_{\rm BH}$, however, are in tension with the prediction of the cosmological coupling of black holes with $k=3$ at a confidence level of $\sim 3\sigma$, which is not in support of the hypothesis that BHs serve as the origin of dark energy. The future observations of high-redshift AGNs by JWST will further test such a hypothesis by identifying more early-type host galaxies in the higher mass range.

Exomoons are expected to orbit gas giant exoplanets just as moons orbit solar system planets. Tidal heating is present in solar system satellites and it can heat up their interior depending on their orbital and interior properties. We aim to identify a Tidally Heated Exomoon's (THEM) orbital parameter space that would make it observable in infrared wavelengths with MIRI/JWST around $\epsilon$ Eridani b. We study the possible constraints on orbital eccentricity and interior properties that a successful THEM detection in infrared wavelengths can bring. We also investigate what exomoon properties need to be independently known in order to place these constraints. We use a coupled thermal-tidal model to find stable equilibrium points between the tidally produced heat and heat transported within a moon. For the latter, we consider a spherical and radially symmetric satellite with heat being transported via magma advection in a sub-layer of melt (asthenosphere) and convection in the lower mantle. We incorporate uncertainties in the interior and tidal model parameters to assess the fraction of simulated moons that would be observable with MIRI. We find that a $2 R_{Io}$ THEM orbiting $\epsilon$ Eridani b with an eccentricity of 0.02, would need to have a semi-major axis of 4 planetary Roche-radii for 100% of the simulations to produce an observable moon. These values are comparable with the orbital properties of gas giant solar system satellites. We place similar constraints for eccentricities up to 0.1. We conclude that if the semi-major axis and radius of the moon are known (eg. with exomoon transits), tidal dissipation can constrain the orbital eccentricity and interior properties of the satellite, such as the presence of melt and the thickness of the melt containing sub-layer.

Cometary dust particles are best preserved remnants of the matter present at the onset of the formation of the Solar System. Space missions, telescopic observations and laboratory analyses advanced the knowledge on the properties of cometary dust. Cometary samples were returned from comet 81P/Wild2 by the Stardust mission. The chondritic (porous) anhydrous interplanetary dust particles and chondritic porous micrometeorites, and the ultracarbonaceous Antarctic micrometeorites (UCAMMs) also show strong evidence for a cometary origin. The composition of cometary dust is generally chondritic, but with high C and N compared with CI. The cometary organic matter is mixed with minor amounts of crystalline and amorphous minerals. The most abundant crystalline minerals are ferromagnesian silicates, refractory minerals and low Ni Fe sulfides are also present. The presence of carbonates in cometary dust is still debated, but a phyllosilicate-like phase was observed in a UCAMM. GEMS phases are usually abundant. Some of the organic matter present in cometary dust particle resembles the insoluble organic matter present in primitive meteorites, but amorphous carbon and exotic (e.g. N-rich) organic phases are also present. The H isotopic composition of the organic matter traces a formation at very low temperatures, in the protosolar cloud or in the outer regions of the protoplanetary disk. The presolar dust concentration in cometary dust can reach about 1%, which is the most elevated value observed in extraterrestrial samples. The differential size distribution of cometary dust in comet trails is well represented by a power-law distribution with a mean power index N typically ranging from -3 to -4. Polarimetric and light scattering studies suggest mixtures of porous agglomerates of sub-micrometer minerals with organic matter. Cometary dust particles have low tensile strength, and low density.

Modern radio telescopes will daily generate data sets on the scale of exabytes for systems like the Square Kilometre Array (SKA). Massive data sets are a source of unknown and rare astrophysical phenomena that lead to discoveries. Nonetheless, this is only plausible with the exploitation of intensive machine intelligence to complement human-aided and traditional statistical techniques. Recently, there has been a surge in scientific publications focusing on the use of artificial intelligence in radio astronomy, addressing challenges such as source extraction, morphological classification, and anomaly detection. This study presents a succinct, but comprehensive review of the application of machine intelligence techniques on radio images with emphasis on the morphological classification of radio galaxies. It aims to present a detailed synthesis of the relevant papers summarizing the literature based on data complexity, data pre-processing, and methodological novelty in radio astronomy. The rapid advancement and application of computer intelligence in radio astronomy has resulted in a revolution and a new paradigm shift in the automation of daunting data processes. However, the optimal exploitation of artificial intelligence in radio astronomy, calls for continued collaborative efforts in the creation of annotated data sets. Additionally, in order to quickly locate radio galaxies with similar or dissimilar physical characteristics, it is necessary to index the identified radio sources. Nonetheless, this issue has not been adequately addressed in the literature, making it an open area for further study.

Measuring the proper motion of the emission component in radio-quiet quasars (RQQs) could help to distinguish between the origins of the radio emission and to understand whether the jet production mechanism is the same in radio-loud quasars (RLQs) and RQQs. PG 1351+640 is one of the few RQQs suitable for proper motion studies: it has two compact components on milli-arcsecond scales, a flat-spectrum core and a steep-spectrum jet; both components are >2 mJy at 5 GHz and are well suited for Very Long Baseline Array (VLBA) observations. We compare recent VLBA observations with that made seventeen years ago and find no significant change in the core-jet separation between 2005 and 2015 (a proper motion of 0.003 mas yr-1). However, the core-jet separation increased significantly between 2015 and 2022, inferring a jet proper motion velocity of 0.063 mas yr-1, which corresponds to an apparent transverse velocity of 0.37c. The result suggests that the jet of the RQQ PG 1351+640 is mildly relativistic and oriented at a relatively small viewing angle.

We report the discovery of the first hot subdwarf B (sdB) star with a massive compact companion in a wide ($P=892.5\pm60.2\,{\rm d}$) binary system. It was discovered based on an astrometric binary solution provided by the Gaia mission Data Release 3. We performed detailed analyses of the spectral energy distribution (SED) as well as spectroscopic follow-up observations and confirm the nature of the visible component as a sdB star. The companion is invisible despite of its high mass of $M_{\rm comp}=1.50_{-0.45}^{+0.37}\,M_{\rm \odot}$. A main sequence star of this mass would significantly contribute to the SED and can be excluded. The companion must be a compact object, either a massive white dwarf or a neutron star. Stable Roche lobe overflow to the companion likely led to the stripping of a red giant and the formation of the sdB, the hot and exposed helium core of the giant. Based on very preliminary data, we estimate that $\sim9\%$ of the sdBs might be formed through this new channel. This binary might also be the prototype for a new progenitor class of supernovae type Ia, which has been predicted by theory.

We discuss the issue of perturbativity in single-field inflationary models with a phase of ultra slow-roll (USR) tailor suited to generate an order-one abundance of primordial black holes (PBHs). More in detail, we impose the condition that loop corrections made up of short-wavelength modes enhanced by the USR dynamics do not alter the tree-level power spectrum of curvature perturbations. In our analysis, the USR phase is preceded and followed by two stages of ordinary slow-roll (SR), and we model the resulting SR/USR/SR dynamics using both instantaneous and smooth transitions. Focusing on scales relevant for CMB observations, we find that it is not possible, with these arguments, to rule out the scenario of PBH formation via USR, not even in the limit of instantaneous transition. However, we also find that loop corrections of short modes on the power spectrum of long modes, even though not large enough to violate perturbativity requirements, remain appreciable and, most importantly, are not tamed in realistic realisations of smooth SR/USR/SR transitions. This makes perturbativity a powerful theoretical tool to constrain USR dynamics. We extend the analysis at any scale beyond those relevant for CMB observations. We find that loop corrections of short modes remain within the few percent if compared to the tree-level power spectrum. However, we also find one notable exception of phenomenological relevance: we show that the so-called dip in the power spectrum of curvature perturbation is an artifact of the tree-level computation.

The photospheric unsigned magnetic flux has been shown to be highly correlated with radial velocity (RV) variations caused by solar surface activity. This activity indicator is therefore a prime candidate to unlock the potential of RV surveys to discover Earth twins orbiting Sun-like stars. We show for the first time how a precise proxy of the unsigned magnetic flux ($\Delta\alpha B^2$) can be obtained from Sun-as-a-star intensity spectra by harnessing the magnetic information contained in over 4000 absorption lines in the wavelength range from 380 to 690 nm. This novel activity proxy can thus be obtained from the same spectra from which RVs are routinely extracted. We derived $\Delta\alpha B^2$ from 500 randomly selected spectra from the HARPS-N public solar data set, which spans from 2015 to 2018. We compared our estimates with the unsigned magnetic flux values from the Solar Dynamics Observatory (SDO) finding excellent agreement (median absolute deviation: 4.9 per cent). The extracted indicator $\Delta\alpha B^2$ correlates with SDO's unsigned magnetic flux estimates on the solar rotational timescale (Pearson correlation coefficient 0.67) and on the three-year timescale of our data set (correlation coefficient 0.91). We find correlations of $\Delta\alpha B^2$ with the HARPS-N solar RV variations of 0.49 on the rotational timescale and 0.78 on the three-year timescale. The Pearson correlation of $\Delta\alpha B^2$ with the RVs is found to be greater than the correlation of the classical activity indicators with the RVs. For solar-type stars, $\Delta\alpha B^2$ therefore represents the best simultaneous activity proxy known to date.

In our work, we examine, for the first time, the possibility of fast and efficient source localization directly from the uvobservations, omitting the recovering of the dirty or clean images. We propose a deep neural network-based framework that takes as its input a low-dimensional vector of sampled uvdata and outputs source positions on the sky. We investigated a representation of the complex-valued input uv-data via the real and imaginary and the magnitude and phase components. We provided a comparison of the efficiency of the proposed framework with the traditional source localization pipeline based on the state-of-the-art Python Blob Detection and Source Finder (PyBDSF) method. The investigation was performed on a data set of 9164 sky models simulated using the Common Astronomy Software Applications (CASA) tool for the Atacama Large Millimeter Array (ALMA) Cycle 5.3 antenna configuration. We investigated two scenarios: (i) noise-free as an ideal case and (ii) sky simulations including noise representative of typical extra-galactic millimeter observations. In the noise-free case, the proposed localization framework demonstrates the same high performance as the state-of-the-art PyBDSF method. For noisy data, however, our new method demonstrates significantly better performance, achieving a completeness level that is three times higher for sources with uniform signal-to-noise (S/N) ratios between 1 and 10, and a high increase in completeness in the low S/N regime. Furthermore, the execution time of the proposed framework is significantly reduced (by factors about 30) as compared to traditional methods that include image reconstructions from the uv-plane and subsequent source detections.

V618 Sgr was previously classified as an R CrB-type variable and later as a possible symbiotic star. Our study aims to analyse the nature of this target, which is currently undergoing significant brightening in properties similar to those of known symbiotic novae. We analyse literature information, photometric observations, and 35 new optical spectra. Our findings strongly suggest that V618 Sgr is an eclipsing symbiotic nova currently in outburst. Additionally, since the star has demonstrated at least two similar brightenings in the past, we propose that V618 Sgr could be the first known galactic symbiotic nova observed in repeated outbursts of this type and may host a relatively massive white dwarf.

I have used high-precision photometry and astrometry from the third data release of Gaia to perform a survey for members of the TW Hya association (TWA). I have identified candidate members that appear to share similar kinematics and ages with bona fide members compiled by Gagne et al. (2017) and I have assessed their membership using radial velocities and spectroscopic diagnostics of age from various sources. My new catalog of adopted members contains 67 Gaia sources in 55 systems. The histogram of spectral types for TWA peaks near M5 (~0.15 Msun), resembling the distributions measured for other nearby young associations. The UVW velocities of its members indicate that the association is expanding. The rate of expansion corresponds to an age of 9.6+0.9/-0.8 Myr. In a Gaia color-magnitude diagram, the members of TWA exhibit well-defined sequences of single stars and unresolved binary stars. The combined sequence of low-mass stars in TWA is indicative of an age of 11.4+1.3/-1.2 Myr when compared to the sequence for Upper Centaurus-Lupus/Lower Centaurus-Crux, for which an age of 20 Myr is assumed. Based on these expansion and isochronal ages, I have adopted an age of 10+/-2 Myr for TWA. Finally, I have used mid-infrared photometry from the Wide-field Infrared Survey Explorer to check for excess emission from circumstellar disks among the TWA members. Fourteen members have detected disks, all of which have been reported in previous studies. The fraction of members at <=M6 (>=0.1 Msun) that have full, transitional, or evolved disks is 10/52=0.19+0.08/-0.06. That value is similar to the fraction previously measured for the Upper Sco association, which is roughly coeval with TWA.

Although neutron star-black hole binaries have been identified through mergers detected in gravitational waves, a pulsar-black hole binary has yet to be detected. While short-period binaries are detectable due to a clear signal in the pulsar's timing residuals, effects from a long-period binary could be masked by other timing effects, allowing them to go undetected. In particular, a long-period binary measured over a small subset of its orbital period could manifest via time derivatives of the spin-frequency incompatible with isolated pulsar properties. We assess the possibility of pulsars having unknown companions in long-period binaries and put constraints on the range of binary properties that may remain undetected in current data, but that may be detectable with further observations. We find that for 35% of canonical pulsars with published higher order derivatives, the precision of measurements is not enough to confidently reject binarity (period greater than ~2 kyr), and that a black-hole binary companion could not be ruled out for a sample of pulsars without published constraints if the period is greater than 1 kyr. While we find no convincing cases in the literature, we put more stringent limits on orbital period and longitude of periastron for the few pulsars with published higher-order frequency derivatives (n greater than 3). We discuss the detectability of candidates and find that a sample pulsar in a 100 yr orbit could be detectable within 5-10 yr.

Planetesimals are born fragile and are subject to destruction by wind erosion as they move through the gas of a protoplanetary disk. In microgravity experiments, we determined the shear stress necessary for erosion of a surface consisting of 1 mm dust pebbles down to 1 Pa ambient pressure. This is directly applicable to protoplanetary disks. Even pebble pile planetesimals with low eccentricities of 0.1 cannot survive inside of 1 au in a minimum-mass solar nebula, and safe zones for planetesimals with higher eccentricities are located even farther out.

Many questions remain about the compressibility of solar wind turbulence with respect to its origins and properties. Low plasma beta (ratio of thermal to magnetic pressure) environments allow for the easier generation of compressible turbulence, enabling study of the relationship between density fluctuations and turbulent Mach number. Utilizing Parker Solar Probe plasma data, we examine the normalized proton density fluctuations $\langle \delta n_p^2 \rangle ^{1/2}/\langle n_p\rangle = \delta {n_p}_{rms}/\langle n_p\rangle$ as a function of turbulent Mach number $M_t$ conditioned on plasma beta and cross helicity. With consideration of statistical error in the parameters computed from in-situ data, we find a general result that $\delta {n_p}_{rms}/\langle n_p\rangle \sim M_t^{1.18 \pm 0.04}$, consistent with both linear-wave theory, and nearly-incompressible turbulence in an inhomogeneous background field. We compare observational results conditioned on plasma beta and cross helicity with 3D magnetohydrodynamic simulations, and observe rather significant similarities with respect to how those parameters affect the proportionality between density fluctuations and turbulent Mach number. This study further investigates the complexity of compressible turbulence as viewed by the density scaling relationship, and may help better understand the compressible environment of the near-Sun solar wind.

Neutron stars are the densest objects in the Universe. They have microscopically homogeneous core and heterogeneous crust. In particular, there may be a specific layer inside neutron stars, the mantle, which consists of substantially non-spherical nuclei immersed in a background of relativistic degenerate electrons and quasi-free neutrons. In this paper we reconsider transverse shear modulus for cylindrical phases of the mantle within the framework of compressible liquid drop model. We demonstrate that transverse shear affects the shape of nuclear clusters: their cross-section becomes elliptical. This effect reduces respective elastic constant. Using a simple model we perform all derivations analytically and obtain the expression for the transverse shear modulus, which can be useful for astrophysical applications.

We report new total intensity visible light high contrast imaging of the TW Hya disk taken with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST). This represents the first published images of the disk with STIS since 2016, when a moving shadow on the disk surface was reported. We continue to see the shadow moving in a counter-clockwise fashion, but in these new images the shadow has evolved into two separate shadows, implying a change in behavior for the occulting structure. Based on radiative transfer models of optically thick disk structures casting shadows, we infer that a plausible explanation for the change is that there are now two misaligned components of the inner disk. The first of these disks is located between 5-6au with an inclination of 5.5\arcdeg and PA of 170\arcdeg, the second between 6-7au with and inclination of 7\arcdeg and PA of 50\arcdeg. Finally, we speculate on the implications of the new shadow structure and determine that additional observations are needed to disentangle the nature of TW Hya's inner disk architecture.

We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249 clusters to 251 GMCs measured from a CARMA $^{12}$CO(1-0) survey with 20 pc resolution. Through two-point correlation analysis, we find that young clusters have a high probability of being near other young clusters, but correlation between GMCs is suppressed by the cloud identification algorithm. By comparing the positions, we find that younger clusters are closer to GMCs than older clusters. Through cross-correlation analysis of the M33 cluster data, we find that clusters are statistically associated when they are $\leq$10 Myr old. Utilizing the high precision ages of the clusters, we find that clusters older than $\approx 18$ Myr are uncorrelated with the molecular ISM. Using the spatial coincidence of the youngest clusters and GMCs in M33, we estimate that clusters spend $\approx$4-6 Myr inside their parent GMC. Through similar analysis, we find that the GMCs in M33 have a total lifetime of $\approx 11$-15 Myr. We also develop a drift model and show that the above correlations can be explained if the clusters in M33 have a 5-10 km s$^{-1}$ velocity dispersion relative to the molecular ISM.

Quasi-Periodic Eruptions (QPEs) are a rare phenomenon in which the X-ray emission from the nuclei of galaxies shows a series of large amplitude flares. Only a handful of QPEs have been observed but the possibility remains that there are as yet undetected sources in archival data. Given the volume of data available a manual search is not feasible, and so we consider an application of machine learning to archival data to determine whether a set of time-domain features can be used to identify further lightcurves containing eruptions. Using a neural network and 14 variability measures we are able to classify lightcurves with accuracies of greater than 94% with simulated data and greater than 98% with observational data on a sample consisting of 12 lightcurves with QPEs and 52 lightcurves without QPEs. An analysis of 83,531 X-ray detections from the XMM Serendipitous Source Catalogue allowed us to recover lightcurves of known QPE sources and examples of several categories of variable stellar objects.

We identify giant molecular clouds (GMCs) associated with HII regions for a sample of 19 nearby galaxies using catalogs of GMCs and H regions released by the PHANGS-ALMA and PHANGS-MUSE surveys, using the overlap of the CO and H{\alpha} emission as the key criterion for physical association. We compare the distributions of GMC and HII region properties for paired and non-paired objects. We investigate correlations between GMC and HII region properties among galaxies and across different galactic environments to determine whether GMCs that are associated with HII regions have significantly distinct physical properties to the parent GMC population. We identify trends between the H{\alpha} luminosity of an HII region and the CO peak brightness and the molecular mass of GMCs that we tentatively attribute to a direct physical connection between the matched objects, and which arise independently of underlying environmental variations of GMC and HII region properties within galaxies. The study of the full sample nevertheless hides a large variability galaxy by galaxy. Our results suggests that at the ~100 pc scales accessed by the PHANGS-ALMA and PHANGS-MUSE data, pre-supernova feedback mechanisms in HII regions have a subtle but measurable impact on the properties of the surrounding molecular gas, as inferred from CO observations.

The locking of lasers to optical cavities is ubiquitously required in the field of precision interferometry such as Advanced LIGO to yield optimal sensitivity. Using higher-order Hermite-Gauss (HG) modes for the main interferometer beam has been a topic of recent study, due to their potential for reducing thermal noise of the test masses. It has been shown however that higher-order HG modes are more susceptible to coupling losses into optical cavities: the misalignment and mode mismatch induced power losses scale as $2n+1$ and $n^{2}+n+1$ respectively with $n$ being the mode index. In this paper we calculate analytically for the first time the alignment and mode mismatch sensing signals for arbitrary higher-order HG modes with both the traditional sensing schemes (using Gouy phase telescopes and quadrant photodetectors) and the more recently proposed radio-frequency jitter-based sensing schemes (using only single element photodiodes). We show that the sensing signals and also the signal-to-shot noise ratios for higher-order HG modes are larger than for the fundamental mode. In particular, the alignment and mode mismatch sensing signals in the traditional sensing schemes scale approximately as $\sqrt{n}$ and $n$ respectively, whereas in the jitter-based sensing schemes they scale exactly as $2n+1$ and $n^{2}+n+1$, respectively, which exactly matches the decrease in their respective tolerances. This potentially mitigates the downside of higher-order HG modes for their suffering from excessive misalignment and mode-mismatch induced power losses.

Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed particle hydrodynamics simulations of two-fluid discs, in the regime where the Stokes number of the particles St>1, to study how the formation of solid clumps depends on the disc-to-star mass ratio, the strength of gravitational instability, and the Stokes number. Gravitational instability of the simulated discs is sustained by local cooling. We find that the ability of the spiral structures to concentrate solids increases with the cooling time, and decreases with the Stokes number, while the relative dynamical temperature between gas and dust of the particles decreases with the cooling time and the disc-to-star mass ratio, and increases with the Stokes number. Dust collapse occurs in a subset of high disc mass simulations, yielding clumps whose mass is close to linear theory estimates, namely 1-10 Earth masses. Our results suggest that if planet formation occurs via this mechanism, the best conditions correspond to near the end of the self-gravitating phase, when the cooling time is long and the Stokes number close to unity.

We present the analysis of archival XMM-Newton and Chandra observations of CH Cyg, one of the most studied symbiotic stars (SySts). The combination of the high-resolution XMM-Newton RGS and Chandra HETG X-ray spectra allowed us to obtain reliable estimates of the chemical abundances and to corroborate the presence of multi-temperature X-ray-emitting gas. Spectral fitting of the medium-resolution XMM-Newton MOS (MOS1+MOS2) spectrum required the use of an additional component not seen in previous studies in order to fit the 2.0-4.0 keV energy range. Detailed spectral modelling of the XMM-Newton MOS data suggests the presence of a reflection component, very similar to that found in active galactic nuclei. The reflection component is very likely produced by an ionised disk (the accretion disk around the white dwarf) and naturally explains the presence of the fluorescent Fe emission line at 6.4 keV while also contributing to the soft and medium energy ranges. The variability of the global X-ray properties of CH Cyg are discussed as well as the variation of the three Fe lines around the 6-7 keV energy range. We conclude that reflection components are needed to model the hard X-ray emission and may be present in most $\beta/\delta$-type SySt.

Odd-indexed higher-order Hermite-Gauss (HG) modes are compatible with 4-quadrant segmented mirrors due to their intensity nulls along the principal axes, which guarantees minimum beam intensity illuminating the bond lines between the segments thus leading to low power loss. However, a misplaced HG beam can cause extra power loss due to the bright intensity spots probing the bond lines. This paper analytically and numerically studies the beam displacement tolerances on a segmented mirror for the $\mathrm{HG_{3,3}}$ mode. We conclude that for "effective" bond lines with 6 $\mu$m width, and the $\mathrm{HG_{3,3}}$ beam size chosen to guarantee 1 ppm clipping loss when centered, the beam can be rotated by roughly 1 degree or laterally displaced by 4% of its beam size while keeping the total power on the bond lines under 1 ppm. We also demonstrate that the constrained beam displacement parameter region that guarantees a given power loss limit, or the beam displacement tolerance, is inversely proportional to the bond line thickness.

Feedback from AGN is known to affect the host galaxy's evolution. In radio AGN, one manifestation of feedback is seen in gas outflows. However, it is still not well understood whether the effect of feedback evolves with the radio AGN life cycle. In this study, we investigate this link using the radio spectral shape as a proxy for the evolutionary stage of the AGN. We used [OIII] emission line spectra to trace the presence of outflows on the ionised gas. Using a sample of uniformly selected 129 radio AGN with $L_\textrm{1.4GHz}\approx10^{23}-10^{26}$ W Hz$^{-1}$, and a mean stacking analysis of the [OIII] profile, we conclude that the ionised gas outflow is linked to the radio spectral shape, and it evolves with the evolution of the radio source. We find that sources with a peak in their radio spectra (optically thick), on average, drive a broad outflow ($FWHM\approx1330\pm418$ km s$^{-1}$) with a velocity $v_\textrm{out}\approx 240$ km s$^{-1}$. However, we detect no outflow in the stacked [OIII] profile of sources without a peak in their radio spectrum. In addition, we find that individual outflow detections are kinematically more extreme in peaked than non-peaked sources. We conclude that radio jets are most effective at driving gas outflows when young, and the outflow is typically short lived. Our stacking analysis shows no significant dependence of the presence of ionised gas outflows on the radio morphology, 1.4 GHz luminosity, optical luminosity and Eddington ratio of these sources. We also identify candidate restarted AGN in our sample, whose [OIII] profiles suggest that they have more disturbed gas kinematics than their evolved counterparts, although the evidence for this is tentative. Our findings support the picture where the impact of AGN feedback changes as the source evolves, and young radio jets interact with the ambient medium, clearing a channel of gas as they expand.

Gravitational freeze-in is a mechanism to explain the observed dark matter relic density if dark matter neither couples to inflation nor to standard model sector. In this work, we study gravitational freeze-in dark matter production during Higgs preheating based on non-perturbative resonance. Using reliable lattice method to handle this non-perturbative process, we show that tachyonic resonance is prohibited by strong back reaction due to Higgs self interaction needed to keep the positivity of potential during preheating, and parameter resonance is viable by tuning the Higgs self-interaction coupling to be small enough in ultraviolet energy scale. We then derive the dark matter relic density under the context of Higgs preheating, and uncover a new dark matter parameter space with dark matter mass larger than inflaton mass, which arises from out-of-equilium Higgs annihilation. Finally, we briefly remark the open question of testing gravitational dark matter.

It has been hypothesized that dark matter is comprised of ultra-light bosons whose collective phenomena can be described as a scalar field undergoing coherent oscillations. Examples include axion and fuzzy dark matter models. In this ultra-light dark matter scenario, the harmonic variation in the field's energy-momentum tensor sources an oscillating component of the gravitational potential that we show can resonantly-excite stellar oscillations. A mathematical framework for predicting the amplitude of these oscillations is developed, which reveals that ultra-light dark matter predominantly excites p-modes of degree $l=1$. An investigation of resonantly-excited solar oscillations is presented, from which we conclude that dark matter-induced oscillations of the Sun are likely undetectable. We discuss prospects for constraining ultra-light dark matter using other stellar objects.

The astrophysical stochastic gravitational wave background (SGWB) originates from numerous faint sub-threshold gravitational wave (GW) signals arising from the coalescing binary compact objects. This background is expected to be discovered from the current (or next-generation) network of GW detectors by cross-correlating the signal between multiple pairs of GW detectors. However, detecting this signal is challenging and the correlation is only detectable at low frequencies due to the arrival time delay between different detectors. In this work, we propose a novel technique, \texttt{Spectrogram Correlated Stacking} (or \texttt{SpeCs}), which goes beyond the usual cross-correlation (and to higher frequencies) by exploiting the higher-order statistics in the time-frequency domain which accounts for the \textit{chirping} nature of the individual events that comprise SGWB. We show that \texttt{SpeCs} improves the signal-to-noise for the detection of SGWB by a factor close to $8$, compared to standard optimal cross-correlation methods which are tuned to measure only the power spectrum of the SGWB signal.\texttt{SpeCs} can probe beyond the power spectrum and its application to the GW data available from the current and next-generation GW detectors would speed up the SGWB discovery.

We present a novel method to constrain the axion-electron coupling constant using the observed calibration of the tip of the red giant branch (TRGB) I band magnitude $M_I$ that fully accounts for uncertainties and degeneracies with stellar input physics.~We simulate a grid of 116,250 models varying initial mass, helium abundance, and metallicity and train a machine learning emulator to predict $M_I$ as a function of these parameters.~Our emulator enables the use of Markov Chain Monte Carlo simulations where the axion-electron coupling $\alpha_{26}$ is varied simultaneously with the stellar parameters. We find that, once stellar uncertainties and degeneracies are accounted for, the region $\alpha_{26} < 2$ is not excluded by empirical TRGB calibrations.~Our work opens up a large region of parameter space currently believed to be excluded.~$\alpha_{26} = 2$ is the upper limit of the parameter space considered by this study, and it is likely that larger values of $\alpha_{26}$ are also unconstrained.~We discuss potential applications of our work to reevaluate other astrophysical probes of new physics.

Multidimensional modification of gravity with a smaller mass scale of the gravitational interaction is considered. Stable by assumption dark matter particles could decay via interactions with virtual black holes. The decay rates of such processes are estimated. It is shown that with the proper fixation of the parameters the decays of these ultra-massive particles can give noticeable contribution to the flux of high energy cosmic rays in particular, near the Greisen-Zatsepin-Kuzmin limit. Such particles can also create neutrinos of very high energies observed in the existing huge underwater or ice-cube detectors.

It is of great interest to understand the equation of state (EOS) of the neutron star (NS), whose core includes highly dense matter. However, there are large uncertainties in the theoretical predictions for the EOS of NS. It is useful to develop a new framework, which is flexible enough to consider the systematic error in theoretical predictions and to use them as a best guess at the same time. We employ a deep neural network to perform a non-parametric fit of the EOS of NS using currently available data. In this framework, the Gaussian process is applied to represent the EOSs and the training set data required to close physical solutions. Our model is constructed under the assumption that the true EOS of NS is a perturbation of the relativistic mean-field model prediction. We fit the EOSs of NS using two different example datasets, which can satisfy the latest constraints from the massive neutron stars, NICER, and the gravitational wave of the binary neutron stars. Given our assumptions, we find that a maximum neutron star mass is $2.38^{+0.15}_{-0.13} M_\odot$ or $2.41^{+0.15}_{-0.14}$ at $95\%$ confidence level from two different example datasets. It implies that the $1.4 M_\odot$ radius is $12.31^{+0.29}_{-0.31}$ km or $12.30^{+0.35}_{-0.37}$ km. These results are consistent with results from previous studies using similar priors. It has demonstrated the recovery of the EOS of NS using a nonparametric model.

Axions and axion-like-particles (ALPs) are characterised by their two-photon coupling, which entails so-called photon-ALP oscillations as photons propagate through a magnetic field. These oscillations lead to distinctive signatures in the energy spectrum of high-energy photons from astrophysical sources, allowing one to probe the existence of ALPs. In particular, photon-ALP oscillations will induce energy dependent oscillatory features, or ``ALP wiggles'', in the photon spectra. We propose to use the discrete power spectrum to search for ALP wiggles and present a model-independent statistical test. By using PKS 2155-304 as an example, we show that the method has the potential to significantly improve the experimental sensitivities for ALP wiggles. Moreover, we discuss how these sensitivities depend on the modelling of the magnetic field. We find that the use of realistic magnetic field models, due to their larger cosmic variance, substantially enhances detection prospects compared to the use of simplified models.

Tilt-to-length coupling was the limiting noise source in LISA Pathfinder between 20 and 200 mHz before subtraction in post-processing. To prevent the adding of sensing noise to the data by the subtraction process, the success of this strategy depended on a previous direct noise reduction by test mass alignment. The exact dependency of the level of tilt-to-length coupling on the set-points of LISA Pathfinder's test masses was not understood until the end of the mission. Here, we present, for the first time, an analytical tilt-to-length coupling model that describes the coupling noise changes due to the realignments. We report on the different mechanisms, namely the lever arm and piston effect as well as the coupling due to transmissive components, and how they contribute to the full coupling. Further, we show that a pure geometric model would not have been sufficient to describe the coupling in LISA Pathfinder. Therefore, we model also the non-geometric tilt-to-length noise contributions. For the resulting coupling coefficients of the full model, we compute the expected error bars based on the known individual error sources. Also, we validated the analytical model against numerical simulations. A detailed study and thorough understanding of this noise are the basis for a successful analysis of the LISA Pathfinder data with respect to tilt-to-length coupling.

The discovery of new, flavor-dependent neutrino interactions would provide compelling evidence of physics beyond the Standard Model. We focus on interactions generated by the anomaly-free, gauged, abelian lepton-number symmetries, specifically $L_e-L_\mu$, $L_e-L_\tau$, and $L_\mu-L_\tau$, that introduce a new matter potential sourced by electrons and neutrons, potentially impacting neutrino flavor oscillations. We revisit, revamp, and improve the constraints on these interactions that can be placed via the flavor composition of the diffuse flux of high-energy astrophysical neutrinos, with TeV-PeV energies, i.e., the proportion of $\nu_e$, $\nu_\mu$, and $\nu_\tau$ in the flux. Because we consider mediators of these new interactions to be ultra-light, lighter than $10^{-10}$ eV, the interaction range is ultra-long, from km to Gpc, allowing vast numbers of electrons and neutrons in celestial bodies and the cosmological matter distribution to contribute to this new potential. We leverage the present-day and future sensitivity of high-energy neutrino telescopes and of oscillation experiments to estimate the constraints that could be placed on the coupling strength of these interactions. We find that, already today, the IceCube neutrino telescope demonstrates potential to constrain flavor-dependent long-range interactions significantly better than existing constraints, motivating further analysis. We also estimate the improvement in the sensitivity due to the next-generation neutrino telescopes such as IceCube-Gen2, Baikal-GVD, KM3NeT, P-ONE, and TAMBO.