Papers for Wednesday, Aug 02 2023

Understanding atmospheric escape in close-in exoplanets is critical to interpreting their evolution. We map out the parameter space over which photoevaporation and core-powered mass loss dominate atmospheric escape. Generally, the transition between the two regimes is determined by the location of the Bondi radius (i.e. the sonic point of core-powered outflow) relative to the penetration depth of XUV photons. Photoevaporation dominates the loss when the XUV penetration depth lies inside the Bondi radius ($R_{XUV}<R_B$) and core-powered mass-loss when XUV radiation is absorbed higher up in the flow ($R_B<R_{XUV}$). The transition between the two regimes occurs at a roughly constant ratio of the planet's radius to its Bondi radius, with the exact value depending logarithmically on planetary and stellar properties. In general, core-powered mass-loss dominates for lower-gravity planets with higher equilibrium temperatures, and photoevaporation dominates for higher-gravity planets with lower equilibrium temperatures. However, planets can transition between these two mass-loss regimes during their evolution, and core-powered mass loss can ``enhance'' photo-evaporation over a significant region of parameter space. Interestingly, a planet that is ultimately stripped by core-powered mass-loss has likely only ever experienced core-powered mass-loss. In contrast a planet that is ultimately stripped by photoevaporation could have experienced an early phase of core-powered mass-loss. Applying our results to the observed super-Earth population suggests that it contains significant fractions of planets where each mechanism controlled the final removal of the H/He envelope, although photoevaporation appears to be responsible for the final carving of the exoplanet radius-valley.

Our current knowledge of star-forming metallicity relies primarily on gas-phase oxygen abundance measurements. This may not allow one to accurately describe differences in stellar evolution and feedback driven by variations in iron abundance. $\alpha$-elements (such as oxygen) and iron are produced by sources that operate on different timescales and the link between them is not straightforward. We explore the origin of the [O/Fe] - specific SFR (sSFR) relation, linking chemical abundances to galaxy formation timescales. This relation is followed by star-forming galaxies across redshifts according to cosmological simulations and basic theoretical expectations. Its apparent universality makes it suitable for trading the readily available oxygen for iron abundance. The relation is determined by the relative iron production efficiency of core-collapse and type Ia supernovae and the delay time distribution of the latter -- uncertain factors that could be constrained empirically with the [O/Fe]-sSFR relation. We compile and homogenise a literature sample of star-forming galaxies with observational iron abundance determinations to place first constraints on the [O/Fe]-sSFR relation over a wide range of sSFR. The relation shows a clear evolution towards lower [O/Fe] with decreasing sSFR and a flattening above log(sSFR/yr)>-9. The result is broadly consistent with expectations, but better constraints are needed to inform the models. We independently derive the relation from old Milky Way stars and find a remarkable agreement between the two, as long as the recombination-line absolute oxygen abundance scale is used in conjunction with stellar metallicity measurements.

The recent detection of TeV neutrino emission from the nearby active galaxy NGC 1068 by IceCube suggests that AGN could make a sizable contribution to the total high-energy cosmic neutrino flux. The absence of TeV gamma rays from NGC 1068, indicates neutrino production originates in the innermost region of the AGN. Disk-corona models predict a correlation between neutrinos and keV X-rays in Seyfert galaxies, a subclass of AGN to which NGC 1068 belongs. Using 10 years of IceCube through-going track events, we report results from searches for neutrino signals from 27 additional sources in the Northern Sky by studying both the generic single power-law spectral assumption and spectra predicted by the disk-corona model. Our results show excesses of neutrinos associated with two sources, NGC 4151 and CGCG 420-015, at 2.7$\sigma$ significance, and at the same time constrain the collective neutrino emission from our source list.

We present the discovery of the eclipsing double white dwarf (WD) binary WDJ 022558.21-692025.38 that has an orbital period of 47.19 min. Following identification with the Transiting Exoplanet Survey Satellite, we obtained time-series ground based spectroscopy and high-speed multi-band ULTRACAM photometry which indicate a primary DA WD of mass 0.40 +- 0.04 Msol and a 0.28 +- 0.02 Msol mass secondary WD, which is likely of type DA as well. The system becomes the third-closest eclipsing double WD binary discovered with a distance of approximately 400 pc and will be a detectable source for upcoming gravitational wave detectors in the mHz frequency range. Its orbital decay will be measurable photometrically within 10 yrs to a precision of better than 1%. The fate of the binary is to merge in approximately 41 Myr, likely forming a single, more massive WD.

The lower energy peak of the spectral energy distribution of blazars has commonly been ascribed to synchrotron radiation from relativistic particles in the jets. Despite the consensus regarding jet emission processes, the particle acceleration mechanism is still debated. Here, we present the first X-ray polarization observations of PG 1553+113, a high-synchrotron-peak blazar observed by the Imaging X-ray Polarimetry Explorer (IXPE). We detect an X-ray polarization degree of $(10\pm2)\%$ along an electric-vector position angle of $\psi_X=86^{\circ}\pm8^{\circ}$. At the same time, the radio and optical polarization degrees are lower by a factor of $\sim$3. During our IXPE pointing, we observed the first orphan optical polarization swing of the IXPE era, as the optical angle of PG 1553+113 underwent a smooth monotonic rotation by about 125$^\circ$, with a rate of $\sim$17 degrees per day. We do not find evidence of a similar rotation in either radio or X-rays, which suggests that the X-ray and optically emitting regions are separate or, at most, partially co-spatial. Our spectro-polarimetric results provide further evidence that the steady-state X-ray emission in blazars originates in a shock-accelerated and energy-stratified electron population.

Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, basic demographics of clumps and their host galaxies remain uncertain. Using the HST/WFC3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (UVCANDELS), we detect and analyze giant star-forming clumps in galaxies at 0.5 $\leq$ z $\leq$ 1, connecting two epochs when clumps are common (at cosmic high-noon, z $\sim$ 2) and rare (in the local universe). We construct a clump sample whose rest-frame 1600 {\AA} luminosity is 3 times higher than the most luminous local HII regions (M$_{UV} \leq -$16 AB). In our sample, 35 $\pm$ 3$\%$ of low-mass galaxies (log[M$_{*}$/M$_{\odot}$] $<$ 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22 $\pm$ 3$\%$ and 22 $\pm$ 4$\%$ for intermediate (10 $\leq$ log[M$_{*}$/M$_{\odot}$] $\leq$ 10.5) and high-mass (log[M$_{*}$/M$_{\odot}$] $>$ 10.5) galaxies in agreement with previous studies. When compared to similar-mass non-clumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher SFRs and bluer rest-frame U-V colors, while high-mass clumpy SFGs tend to be larger than non-clumpy SFGs. However, clumpy and non-clumpy SFGs have similar S\'ersic index, indicating a similar underlying density profile. Furthermore, we investigate how UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher-specific SFR galaxies and are found closer to their hosts' galactic center.

We present JWST/NIRSpec observations of a highly magnified star candidate at a photometric redshift of $z_{\mathrm{phot}}\simeq4.8$, previously detected in JWST/NIRCam imaging of the strong lensing (SL) cluster MACS J0647+7015 ($z=0.591$). The spectroscopic observation allows us to precisely measure the redshift of the host arc at $z_{\mathrm{spec}}=4.758\pm0.004$, and the star's spectrum displays clear Lyman- and Balmer-breaks commensurate with this redshift. A fit to the spectrum suggests a B-type super-giant star of surface temperature $T_{\mathrm{eff,B}}\simeq15000$ K with either a redder F-type companion ($T_{\mathrm{eff,F}}\simeq6250$K) or significant dust attenuation ($A_V\simeq0.82$) along the line of sight. We also investigate the possibility that this object is a magnified young globular cluster rather than a single star. We show that the spectrum is in principle consistent with a star cluster, which could also accommodate the lack of flux variability between the two epochs. However, the lack of a counter image and the strong upper limit on the size of the object from lensing symmetry, $r\lesssim0.5$ pc, could indicate that this scenario is somewhat less likely -- albeit not completely ruled out by the current data. The presented spectrum seen at a time when the Universe was only $\sim1.2$ Gyr old showcases the ability of JWST to study early stars through extreme lensing.

We characterize the stellar mass of J2239+0207, a z~6.25 sub-Eddington quasar (M_1450=-24.6), using dedicated JWST/NIRCam medium-band observations of a nearby PSF star to remove the central point source and reveal the underlying galaxy emission. We detect the host galaxy in two bands longward of the Balmer break, obtaining a stellar mass of ~10^10 M_sun, more than an order of magnitude less than this quasar's existing measured [C II] dynamical mass. We additionally calculate the mass of J2239+0207's central supermassive black hole using JWST/NIRSpec IFU observations, and determine that the black hole is ~15 times more massive than predicted by the local M_BH-M* relation, similar to many high-redshift quasars with dynamical masses determined via millimeter-wave line widths. We carefully consider potential selection effects at play, and find that even when z~6 quasars are compared to a local sample with similarly determined dynamical masses, many of the high-redshift quasars appear to possess overmassive black holes. We conclude z~6 quasars are likely to have a larger spread about the M_BH-M* relation than observed in the local Universe.

Techniques to retrieve the atmospheric properties of exoplanets via direct observation of their reflected light have often been limited in scope due to computational constraints imposed by the forward-model calculations. We have developed a new set of techniques which significantly decreases the time required to perform a retrieval while maintaining accurate results. We constructed a grid of 1.4 million pre-computed geometric albedo spectra valued at discrete sets of parameter points. Spectra from this grid are used to produce models for a fast and efficient nested sampling routine called PSGnest. Beyond the upfront time to construct a spectral grid, the amount of time to complete a full retrieval using PSGnest is on the order of seconds to minutes using a personal computer. An extensive evaluation of the error induced from interpolating intermediate spectra from the grid indicates that this bias is insignificant compared to other retrieval error sources, with an average coefficient of determination between interpolated and true spectra of 0.998. We apply these new retrieval techniques to help constrain the optimal bandpass centers for retrieving various atmospheric and bulk parameters from a LuvEx-type mission observing several planetary archetypes. We show that spectral observations made using a 20\% bandpass centered at 0.73 microns can be used alongside our new techniques to make detections of $H_2O$ and $O_2$ without the need to increase observing time beyond what is necessary for a signal-to-noise ratio of 10. The methods introduced here will enable robust studies of the capabilities of future observatories to characterize exoplanets.

We describe a novel method to compute the components of dynamo tensors from direct magnetohydrodynamic (MHD) simulations. Our method relies upon an extension and generalisation of the standard H\"ogbom CLEAN algorithm widely used in radio astronomy to systematically remove the impact of the strongest beams onto the corresponding image. This generalisation, called the Iterative Removal of Sources (IROS) method, has been adopted here to model the turbulent electromotive force (EMF) in terms of the mean magnetic fields and currents. Analogous to the CLEAN algorithm, IROS treats the time series of the mean magnetic field and current as beams that convolve with the dynamo coefficients which are treated as (clean) images to produce the EMF time series (the dirty image). We apply this method to MHD simulations of galactic dynamos, to which we have previously employed other methods of computing dynamo coefficients such as the test-field method, the regression method, as well as local and non-local versions of the singular value decomposition (SVD) method. We show that our new method reliably recovers the dynamo coefficients from the MHD simulations. It also allows priors on the dynamo coefficients to be incorporated easily during the inversion, unlike in earlier methods. Moreover, using synthetic data, we demonstrate that it may serve as a viable post-processing tool in determining the dynamo coefficients, even when the power of additive noise to the EMF is twice as much the actual EMF.

Spatiotemporal techniques for signal coordination with actively transmitting extraterrestrial civilizations, without the need for prior communication, can constrain technosignature searches to a significantly smaller coordinate space. With the variable star catalog from Gaia Data Release 3, we explore two related signaling strategies: the SETI Ellipsoid, and that proposed by Seto, which are both based on the synchronization of transmissions with a conspicuous astrophysical event. This dataset contains more than 10 million variable star candidates with light curves from the first three years of Gaia's operational phase, between 2014 and 2017. Using four different historical supernovae as source events, we find that less than 0.01% of stars in the sample have crossing times, the times at which we would expect to receive synchronized signals on Earth, within the date range of available Gaia observations. For these stars, we present a framework for technosignature analysis that searches for modulations in the variability parameters by splitting the stellar light curve at the crossing time.

We present the detection of 661 known pulsars observed with the Australian SKA Pathfinder (ASKAP) telescope at 888 MHz as a part of the Rapid ASKAP Continuum Survey (RACS). Detections were made through astrometric coincidence and we estimate the false alarm rate of our sample to be ~0.5%. Using archival data at 400 and 1400 MHz, we estimate the power law spectral indices for the pulsars in our sample and find that the mean spectral index is -1.78 +/- 0.6. However, we also find that a single power law is inadequate to model all the observed spectra. With the addition of the flux densities between 150 MHz and 3 GHz from various imaging surveys, we find that up to 40% of our sample shows deviations from a simple power law model. Using Stokes V measurements from the RACS data, we measured the circular polarization fraction for 9% of our sample and find that the mean polarization fraction is ~10% (consistent between detections and upper limits). Using the dispersion measure (DM) derived distance we estimate the pseudo luminosity of the pulsars and do not find any strong evidence for a correlation with the pulsars' intrinsic properties.

We present Kepler exoplanet occurrence rates for planets between $0.5-16$ R$_\oplus$ and between $1-400$ days. To measure occurrence, we use a non-parametric method via a kernel density estimator and use bootstrap random sampling for uncertainty estimation. We use a full characterization of completeness and reliability measurements from the Kepler DR25 catalog, including detection efficiency, vetting completeness, astrophysical- and false alarm reliability. We also include more accurate and homogeneous stellar radii from Gaia DR2. In order to see the impact of these final Kepler properties, we revisit benchmark exoplanet occurrence rate measurements from the literature. We compare our measurements with previous studies to both validate our method and observe the dependence of these benchmarks on updated stellar and planet properties. For FGK stars, between $0.5-16$ R$_\oplus$ and between $1-400$ days, we find an occurrence of $1.52\pm0.08$ planets per star. We investigate the dependence of occurrence as a function of radius, orbital period, and stellar type and compare with previous studies with excellent agreement. We measure the minimum of the radius valley to be $1.78^{+0.14}_{-0.16}$ R$_\oplus$ for FGK stars and find it to move to smaller radii for cooler stars. We also present new measurements of the slope of the occurrence cliff at $3-4$ R$_\oplus$, and find that the cliff becomes less steep at long orbital period. Our methodology will enable us to constrain theoretical models of planet formation and evolution in the future.

Neutron stars are identified as pulsars, X-ray binary components, central objects of supernovae remnants, or isolated thermally emitting sources, and at distances beyond 120 pc. A population extrapolation suggests 10$^3$ objects within that boundary. Potentially, neutron stars could continuously emit gravitational waves at sensitivity reach of present instrumentation. As part of our Search for the Nearest Neutron Stars ``Five Seasons'' project, we search for nearby resolved neutron stars. Based on expected fluxes and magnitudes of thermally cooling neutron stars and pulsars, we selected sources in Gaia DR3. The sources have $G$-band absolute magnitudes $M_G>16$ mag, parallax signal-to-noise ratios greater than two, and colours $G_{BP}-G<0.78$ and $G-G_{RP}<0.91$ mag for power-law emitters of flux $F_{\nu} \propto \nu^{-\alpha_{\nu}}$ with spectral indices $\alpha_{\nu}<3$. The photometric region overlaps with that of white dwarfs, in confluence with most known pulsars in binaries having white dwarf companions. We looked for counterparts in gamma-ray, X-ray, ultraviolet, radio, optical, and infrared catalogues. We find about two X-ray-, 15 ultraviolet-, one radio probable counterparts, and at least four sources with power-law profiles at the ultraviolet-optical(-infrared). Because the sources have $G\gtrapprox20$ mag, we rely on Gaia DR3 single-source parameters. We identify possible binaries based on photoastrometric parameters, visual companions, and flux excesses. Some emission components suggest small thermal radii. Source types, neutron star content, and properties require further inquiry.

Time-domain observations are increasingly important in astronomy, and are often the only way to study certain objects. The volume of time-series data is increasing dramatically as new surveys come online - for example, the Vera Rubin Observatory will produce 15 terabytes of data per night, and its Legacy Survey of Space and Time (LSST) is expected to produce five-year lightcurves for $>10^7$ sources, each consisting of 5 photometric bands. Historically, astronomers have worked with Fourier-based techniques such as the Lomb-Scargle periodogram or information-theoretic approaches; however, in recent years Bayesian and data-driven approaches such as Gaussian Process Regression (GPR) have gained traction. However, the computational complexity and steep learning curve of GPR has limited its adoption. `pgmuvi` makes GPR of multi-band timeseries accessible to astronomers by building on cutting-edge open-source machine-learning libraries, and hence `pgmuvi` retains the speed and flexibility of GPR while being easy to use. It provides easy access to GPU acceleration and Bayesian inference of the hyperparameters (e.g. the periods), and is able to scale to large datasets.

This paper presents a comprehensive X-ray and multiwavelength study of three ultraluminous X-ray sources (ULXs) in NGC 1566 using available archival Chandra, Swift-XRT, James Webb Space Telescope, and Hubble Space Telescope observations. The X-ray count rates of the ULX-1 and ULX-2 and ULX-3 were seen to vary by over three orders of magnitude while count rates of ULX-3 varied almost an order of magnitude in long-term light curves. X-ray luminosities of ULX-1 and ULX-2 exceed 10^40 erg/s in 0.3 10 keV energy band as seen in hyperluminous X-ray sources. Obvious evidence for the existence of bi-modality features in the count rate distribution of ULX-1 which may indicate the presence of the propeller regime. For ULX-1, a 201-day super-orbital or orbital periodicity candidate was detected. The transition track of ULX-1 and ULX-2 looks like for atoll sources from the hardness intensity diagrams while the transition track of ULX-3 is most likely spectral transition seen in black holes. Resulting from precise astrometric calculations, a unique both optical and NIR (Near Infrared) counterpart, and two both optical and NIR counterparts were identified for ULX-1 and ULX-2 while a unique optical counterpart was determined for ULX-3. Scenarios for the possible donor of the ULX-1 indicate that it could be either a red supergiant or a red giant. Findings for the optical counterpart of ULX-3 may show that it is a typical OB-type or blue supergiant. The emission counterparts of ULX-2 appear to be dominated by the outer region of the accretion disk.

Solar flares are driven by release of the free magnetic energy and its conversion to other forms of energy -- kinetic, thermal, and nonthermal. Quantification of partitions between these energy components and their evolution is needed to understand the solar flare phenomenon including nonthermal particle acceleration, transport, and escape and the thermal plasma heating and cooling. The challenge of remote sensing diagnostics is that the data are taken with finite spatial resolution and suffer from line-of-sight (LOS) ambiguity including cases when different flaring loops overlap and project one over the other. Here we address this challenge by devising a data-constrained evolving 3D model of a multi-loop SOL2014-02-16T064620 solar flare of GOES class C1.5. Specifically, we employed a 3D magnetic model validated earlier for a single time frame and extended it to cover the entire flare evolution. For each time frame we adjusted the distributions of the thermal plasma and nonthermal electrons in the model, such as the observables synthesized from the model matched the observations. Once the evolving model has been validated this way, we computed and investigated the evolving energy components and other relevant parameters by integrating over the model volume. This approach removes the LOS ambiguity and permits to disentangle contributions from the overlapping loops. It reveals new facets of electron acceleration and transport, as well as heating and cooling the flare plasma in 3D. We find signatures of substantial direct heating of the flare plasma not associated with the energy loss of nonthermal electrons.

Crystalline materials are promising candidates as substrates or high-reflective coatings of mirrors to reduce thermal noises in future laser interferometric gravitational wave detectors. However, birefringence of such materials could degrade the sensitivity of gravitational wave detectors, not only because it can introduce optical losses, but also because its fluctuations create extra phase noise in the arm cavity reflected beam. In this paper, we analytically estimate the effects of birefringence and its fluctuations in the mirror substrate and coating for gravitational wave detectors. Our calculations show that the requirements for the birefringence fluctuations in silicon substrate and AlGaAs coating will be in the order of $10^{-8}$ rad/$\sqrt{\rm Hz}$ and $10^{-10}$ rad/$\sqrt{\rm Hz}$ at 100~Hz, respectively, for future gravitational wave detectors. We also point out that optical cavity response needs to be carefully taken into account to estimate optical losses from depolarization.

In a previous work, we have identified for the first time the spin of the dominant black hole of a binary from its jet properties. Analysing Very Long Baseline Array (VLBA) observations of the quasar S5 1928+738, taken at $15$ GHz during $43$ epochs between $1995.96$ and $2013.06$, we showed that the inclination angle variation of the inner ($<2$ mas) jet symmetry axis naturally decomposes into a periodic and a monotonic contribution. The former emerges due to the Keplerian orbital evolution, while the latter is interpreted as the signature of the spin-orbit precession of the jet emitting black hole. In this paper, we revisit the analysis of the quasar S5 1928+738 by including new $15$ GHz VLBA observations extending over $29$ additional epochs, between $2013.34$ and $2020.89$. The extended data set confirms our previous findings which are further supported by the flux density variation of the jet. By applying an enhanced jet precession model that can handle arbitrary spin orientations $\kappa$ with respect to the orbital angular momentum of a binary supermassive black hole system, we estimate the binary mass ratio as $\nu=0.21\pm0.04$ for $\kappa=0$ (i.e. when the spin direction is perpendicular to the orbital plane) and as $\nu=0.32\pm0.07$ for $\kappa=\pi/2$ (i.e. when the spin lies in the orbital plane). We estimate more precisely the spin precession velocity, cutting half of its uncertainty from $(-0.05\pm0.02)^\mathrm{o} \mathrm{yr}^{-1}$ to $(-0.04\pm0.01)^\mathrm{o} \mathrm{yr}^{-1}$.

The planned in-ice radio array of IceCube-Gen2 at the South Pole will provide unprecedented sensitivity to ultra-high-energy (UHE) neutrinos in the EeV range. The ability of the detector to measure the neutrino's energy and direction is of crucial importance. This contribution presents an end-to-end reconstruction of both of these quantities for both detector components of the hybrid radio array ('shallow' and 'deep') using deep neural networks (DNNs). We are able to predict the neutrino's direction and energy precisely for all event topologies, including the electron neutrino charged-current interactions, which are more complex due to the LPM effect. This highlights the advantages of DNNs for modeling the complex correlations in radio detector data, thereby enabling a measurement of the neutrino energy and direction. We discuss how we can use normalizing flows to predict the PDF for each individual event which allows modeling the complex non-Gaussian uncertainty contours of the reconstructed neutrino direction. Finally, we discuss how this work can be used to further optimize the detector layout to improve its reconstruction performance.

The IceCube Neutrino Observatory has discovered a diffuse neutrino flux of astrophysical origin and measures its properties in various detection channels. With more than 10 years of data, we use multiple data samples from different detection channels for a combined fit of the diffuse astrophysical neutrino spectrum. This leverages the complementary information of different neutrino event signatures. For the first time, we use a coherent modelling of the signal and background, as well as the detector response and corresponding systematic uncertainties. The detector response is continuously varied during the simulation in order to generate a general purpose Monte Carlo set, which is central to our approach. We present a combined fit yielding a measurement of the diffuse astrophysical neutrino flux properties with unprecedented precision.

We present new JWST-NIRSpec IFS data for the luminous infrared galaxy NGC7469: a nearby (70.6Mpc) active galaxy with a Sy 1.5 nucleus that drives a highly ionized gas outflow and a prominent nuclear star-forming ring. Using the superb sensitivity and high spatial resolution of the JWST instrument NIRSpec-IFS, we investigate the role of the Seyfert nucleus in the excitation and dynamics of the circumnuclear gas. Our analysis focuses on the [Fe ii], H2, and hydrogen recombination lines that trace the radiation/shocked-excited molecular and ionized ISM around the AGN. We investigate the gas excitation through H2/Br{\gamma} and [Fe ii]/Pa\b{eta} emission line ratios and find that photoionization by the AGN dominates within the central 300 pc of the galaxy and together with a small region show ing signatures of shock-heated gas; these shock-heated regions are likely associated with a compact radio jet. In addition, the velocity field and velocity dispersion maps reveal complex gas kinematics. Rotation is the dominant feature, but we also identify non-circular motions consistent with gas inflows as traced by the velocity residuals and the spiral pattern in the Pa{\alpha} velocity dispersion map. The inflow is consistent with the mass outflow rate and two orders of magnitude higher than the AGN accretion rate. The compact nuclear radio jet has enough power to drive the highly ionized outflow. This scenario suggests that the inflow and outflow are in a self-regulating feeding-feedback process, with a contribution from the radio jet helping to drive the outflow.

Stars exhibit a bewildering variety of rapidly variable behaviors ranging from explosive magnetic flares to stochastically changing accretion to periodic pulsations or rotation. The principal Rubin Observatory Legacy Survey of Space and Time (LSST) surveys will have cadences too sparse and irregular to capture many of these phenomena. We propose here a LSST micro-survey to observe a single Galactic field, rich in unobscured stars, in a continuous sequence of 30 second exposures for one long winter night in a single photometric band. The result will be a unique dataset of $\sim 1$ million regularly spaced stellar light curves (LCs). The LCs will constitute a comprehensive collection of late-type stellar flaring, but also other classes like short-period binary systems and cataclysmic variables, young stellar objects and ultra-short period exoplanets. An unknown variety of anomalous Solar System, Galactic and extragalactic variables and transients may also be present. A powerful array of statistical procedures can be applied to individual LCs from the long-standing fields of time series analysis, signal processing and econometrics. Dozens of `features' describing the variability can be extracted and the ensemble of light curves can be subject to advanced machine learning clustering procedures. This will give a unique, authoritative, objective taxonomy of the rapidly variable sky derived from identically cadenced LCs. This micro-survey is best performed early in the Rubin Observatory program, and the results can inform the wider community on the best approaches to variable star identification and classification from the sparse, irregular cadences that dominate the planned surveys.

Popular wisdom suggests that measuring the tensor to scalar ratio $r$ on CMB scales is a "proof of inflation" since one generic prediction is a scale-invariant tensor spectrum while alternatives predict $r$ that is many orders of magnitude below the sensitivity of future experiments. A bouncing Universe with sourced fluctuations allows for nearly scale-invariant spectra of both scalar and tensor perturbations challenging this point of view. Past works have analyzed the model until the bounce, under the assumption that the bounce will not change the final predictions. In this work, we discard this assumption. We explicitly follow the evolution of the Universe and fluctuations across the bounce until reheating. The evolution is stable, and the existence of the sourced fluctuations does not destroy the bounce. The bounce enhances the scalar spectrum while leaving the tensor spectrum unchanged. The enhancement depends on the duration of the bounce - a shorter bounce implies a larger enhancement. The model matches current observations and predicts any viable tensor-to-scalar ratio $r\lesssim 10^{-2}$, which may be observed in upcoming CMB experiments. Hence, a measurement of $r$ will no longer be a "proof of inflation'', and a Sourced Bounce is a viable paradigm with distinct predictions.

The absorption signals of metastable He in HD 209458b and several other exoplanets can be explained via escaping atmosphere model with a subsolar He/H ratio. The low abundance of helium can be a result of planet formation if there is a small amount of helium in their primordial atmosphere. However, another possibility is that the low He/H ratio is caused by the process of mass fractionation of helium in the atmosphere. In order to investigate the effect of the fractionation in the hydrogen-helium atmosphere, we developed a self-consistent multi-fluid 1D hydrodynamic model based on the well-known open-source MHD code PLUTO. Our simulations show that a lower He/H ratio can be produced spontaneously in the multi-fluid model. We further modeled the transmission spectra of He 10830 lines for HD 209458b in a broad parameter space. The transmission spectrum of the observation can be fitted in the condition of 1.80 times the X-ray and extreme-ultraviolet flux of the quiet Sun. Meanwhile, the ratio of the escaping flux of helium to hydrogen, $F_{He}/F_{H}$, is 0.039. Our results indicate that the mass fractionation of helium to hydrogen can naturally interpret the low He/H ratio required by the observation. Thus, in the escaping atmosphere of HD 209458b, decreasing the abundance of helium in the atmosphere is not needed even if its He abundance is similar to that of the Sun. The simulation presented in this work hints that in the escaping atmosphere, mass fractionation can also occur on other exoplanets, which needs to be explored further.

Compact elliptical (cE) galaxies remain an elusively difficult galaxy class to study. Recent observations have suggested that isolated and host-associated cEs have different formation pathways, while simulation studies have also shown different pathways can lead to a cE galaxy. However a solid link has not been established, and the relative contributions of each pathway in a cosmological context remains unknown. Here we combine a spatially-resolved observational sample of cEs taken from the SAMI galaxy survey with a matched sample of galaxies within the IllustrisTNG cosmological simulation to establish an overall picture of how these galaxies form. The observed cEs located near a host galaxy appear redder, smaller and older than isolated cEs, supporting previous evidence for multiple formation pathways. Tracing the simulated cEs back through time, we find two main formation pathways; 32 $\pm$ 5 percent formed via the stripping of a spiral galaxy by a larger host galaxy, while 68 $\pm$ 4 percent formed through a gradual build-up of stellar mass in isolated environments. We confirm that cEs in different environments do indeed form via different pathways, with all isolated cEs in our sample having formed via in-situ formation (i.e. none were ejected from a previous host), and 77 $\pm$ 6 percent of host-associated cEs having formed via tidal stripping. Separating them by their formation pathway, we are able to reproduce the observed differences between isolated and host-associated cEs, showing that these differences can be fully explained by the different formation pathways dominating in each environment.

When subhalos infall into galaxy clusters, their gas content is ram pressure stripped by the intracluster medium (ICM) and may turn into cometary tails. We report the discovery of two spectacular X-ray double tails in a single galaxy cluster, Z8338, revealed by 70 ks Chandra observations. The brighter one, with an X-ray bolometric luminosity of $3.9 \times 10^{42}{\rm\ erg\ s}^{-1}$, is a detached tail stripped from the host halo and extended at least 250 kpc in projection. The head of the detached tail is a cool core with the front tip of the cold front $\sim$ 30 kpc away from the nucleus of its former host galaxy. The cooling time of the detached cool core is $\sim 0.3$ Gyr. For the detached gas, the gravity of the once-associated dark matter halo further enhances the Rayleigh-Taylor (RT) instability. From its survival, we find that a magnetic field of a few $\mu$G is required to suppress the hydrodynamic instability. The X-ray temperature in the tail increases from 0.9 keV at the front tip to 1.6 keV in the wake region, which suggests the turbulent mixing with the hotter ICM. The fainter double X-ray tail, with a total X-ray luminosity of $2.7 \times 10^{42}{\rm\ erg\ s}^{-1}$, appears to stem from the cool core of a subcluster in Z8338, and likely was formed during the ongoing merger. This example suggests that X-ray cool cores can be displaced and eventually destroyed by mergers, while the displaced cool cores can survive for some extended period of time.

The volume of research on fast radio bursts (FRBs) observation have been seeing a dramatic growth. To facilitate the systematic analysis of the FRB population, we established a database platform, Blinkverse (https://blinkverse.alkaidos.cn), as a central inventory of FRBs from various observatories and with published properties, particularly dynamic spectra from FAST, CHIME, GBT, Arecibo, etc. Blinkverse thus not only forms a superset of FRBCAT, TNS, and CHIME/FRB, but also provides convenient access to thousands of FRB dynamic spectra from FAST, some of which were not available before. Blinkverse is regularly maintained and will be updated by external users in the future. Data entries of FRBs can be retrieved through parameter searches through FRB location, fluence, etc., and their logical combinations. Interactive visualization was built into the platform. We analyzed the energy distribution, period analysis, and classification of FRBs based on data downloaded from Blinkverse. The energy distributions of repeaters and non-repeaters are found to be distinct from one another.

The puzzling complexity of terrestrial biomolecules is driving the search for complex organic molecules in the Interstellar Medium (ISM) and serves as a motivation for many in situ studies of reservoirs of extraterrestrial organics from meteorites and interplanetary dust particles (IDPs) to comets and asteroids. Comet 67P/Churyumov-Gerasimenko (67P) -- the best-studied comet to date -- has been visited and accompanied for two years by the European Space Agency's Rosetta spacecraft. Around 67P's perihelion and under dusty conditions, the high-resolution mass spectrometer on board provided a spectacular glimpse into this comet's chemical complexity. For this work, we analyzed in unprecedented detail the O-bearing organic volatiles. In a comparison of 67P's inventory to molecules detected in the ISM, in other comets, and in Soluble Organic Matter (SOM) extracted from the Murchison meteorite, we also highlight the (pre)biotic relevance of different chemical groups of species. We report first evidence for abundant extraterrestrial O-bearing heterocycles (with abundances relative to methanol often on the order of 10% with a relative error margin of 30-50%) and various representatives of other molecule classes such as carboxylic acids and esters, aldehydes, ketones, and alcohols. Like with the pure hydrocarbons, some hydrogenated forms seem to be dominant over their dehydrogenated counterparts. An interesting example is tetrahydrofuran (THF) as it might be a more promising candidate for searches in the ISM than the long-sought furan itself. Our findings not only support and guide future efforts to investigate the origins of chemical complexity in space, but also they strongly encourage studies of, e.g., the ratios of unbranched vs. branched and hydrogenated vs. dehydrogenated species in astrophysical ice analogs in the laboratory as well as by modeling.

J1044+0353 is considered a local analog of the young galaxies that ionized the intergalactic medium at high-redshift due to its low mass, low metallicity, high specific star formation rate, and strong high-ionization emission lines. We use integral field spectroscopy to trace the propagation of the starburst across this small galaxy using Balmer emission- and absorption-line equivalent widths and find a post-starburst population (~ 15 - 20 Myr) roughly one kpc east of the much younger, compact starburst (~ 3 - 4 Myr). Using the direct electron temperature method to map the O/H abundance ratio, we find similar metallicity (1 to 3 sigma) between the starburst and post-starburst regions but with a significant dispersion of about 0.3 dex within the latter. We also map the Doppler shift and width of the strong emission lines. Over scales several times the size of the galaxy, we discover a velocity gradient parallel to the galaxy's minor axis. The steepest gradients (~ 30 $\mathrm{km \ s^{-1} \ kpc^{-1}}$) appear to emanate from the oldest stellar association. We identify the velocity gradient as an outflow viewed edge-on based on the increased line width and skew in a biconical region. We discuss how this outflow and the gas inflow necessary to trigger the starburst affect the chemical evolution of J1044+0353. We conclude that the stellar associations driving the galactic outflow are spatially offset from the youngest association, and a chemical evolution model with a metal-enriched wind requires a more realistic inflow rate than a homogeneous chemical evolution model.

Fully analytical dynamical models usually have an infinite extent, while real star clusters, galaxies, and dark matter haloes have a finite extent. The standard method for generating dynamical models with a finite extent consists of taking a model with an infinite extent and applying a truncation in binding energy. This method, however, cannot be used to generate models with a pre-set analytical mass density profile. We investigate the self-consistency and dynamical properties of a family of power-law spheres with a general tangential Cuddeford (TC) orbital structure. By varying the density power-law slope $\gamma$ and the central anisotropy $\beta_0$, these models cover a wide parameter space in density and anisotropy profiles. We explicitly calculate the phase-space distribution function for various parameter combinations, and interpret our results in terms of the energy distribution of bound orbits. We find that truncated power-law spheres can be supported by a TC orbital structure if and only if $\gamma \geqslant 2\beta_0$, which means that the central density slope-anisotropy inequality is both a sufficient and a necessary condition for this family. We provide closed expressions for structural and dynamical properties such as the radial and tangential velocity dispersion profiles, which can be compared against more complex numerical modelling results. This work significantly adds to the available suite of self-consistent dynamical models with a finite extent and an analytical description.

Black hole X-ray binaries are ideal environments to study the accretion phenomena in strong gravitational potentials. These systems undergo dramatic accretion state transitions and analysis of the X-ray spectra is used to probe the properties of the accretion disc and its evolution. In this work, we present a systematic investigation of $\sim$1800 spectra obtained by RXTE PCA observations of GRO J1655-40 and LMC X-3 to explore the nature of the accretion disc via non-relativistic and relativistic disc models describing the thermal emission in black-hole X-ray binaries. We demonstrate that the non-relativistic modelling throughout an outburst with the phenomenological multi-colour disc model DISKBB yields significantly lower and often unphysical inner disc radii and correspondingly higher ($\sim$50-60\%) disc temperatures compared to its relativistic counterparts KYNBB and KERRBB. We obtained the dimensionless spin parameters of $a_{*}=0.774 \pm 0.069 $ and $a_{*}=0.752 \pm 0.061 $ for GRO J1655-40 with KERRBB and KYNBB, respectively. We report a spin value of $a_{*}=0.098 \pm 0.063$ for LMC X-3 using the updated black hole mass of 6.98 ${M_{\odot}}$. Both measurements are consistent with the previous studies. Using our results, we highlight the importance of self-consistent modelling of the thermal emission, especially when estimating the spin with the continuum-fitting method which assumes the disc terminates at the innermost stable circular orbit at all times.

Internal shocks are a leading candidate for the dissipation mechanism that powers the prompt $\gamma$-ray emission in gamma-ray bursts (GRBs). In this scenario a compact central source produces an ultra-relativistic outflow with varying speeds, causing faster parts or shells to collide with slower ones. Each collision produces a pair of shocks -- a forward shock (FS) propagating into the slower leading shell and a reverse shock (RS) propagating into the faster trailing shell. The RS's lab-frame speed is always smaller, while the RS is typically stronger than the FS, leading to different conditions in the two shocked regions that both contribute to the observed emission. We show that optically-thin synchrotron emission from both (weaker FS + stronger RS) can naturally explain key features of prompt GRB emission such as the pulse shapes, time-evolution of the $\nu{}F_\nu$ peak flux and photon-energy, and the spectrum. Particularly, it can account for two features commonly observed in GRB spectra: (i) a sub-dominant low-energy spectral component (often interpreted as ``photospheric''-like), or (ii) a doubly-broken power-law spectrum with the low-energy spectral slope approaching the slow cooling limit. Both features can be obtained while maintaining high overall radiative efficiency without any fine-tuning of the physical conditions.

In this study, we perform a halo-finder code comparison between Rockstar and CompaSO. Based on our previous analysis aiming at quantifying resolution of $N$-body simulations by exploiting large (up to $N=4096^3$) simulations of scale-free cosmologies run using Abacus, we focus on convergence of the HMF, 2PCF and mean radial pairwise velocities of halo centres selected with the aforementioned two algorithms. We establish convergence, for both Rockstar and CompaSO, of mass functions at the $1\%$ precision level and of the mean pairwise velocities (and also 2PCF) at the $2\%$ level. At small scales and small masses, we find that Rockstar exhibits greater self-similarity, and we also highlight the role played by the merger-tree post-processing of CompaSO halos on their convergence. Finally, we give resolution limits expressed as a minimum particle number per halo in a form that can be directly extrapolated to LCDM.

Non-thermal emission from relativistic electrons gives insight into the strength and morphology of intra-cluster magnetic fields, as well as providing powerful tracers of structure formation shocks. Emission caused by Cosmic Ray (CR) protons on the other hand still challenges current observations and is therefore testing models of proton acceleration at intra-cluster shocks. Large-scale simulations including the effects of CRs have been difficult to achieve and have been mainly reduced to simulating an overall energy budget, or tracing CR populations in post-processing of simulation output and has often been done for either protons or electrons. We use an efficient on-the-fly Fokker-Planck solver to evolve distributions of CR protons and electrons within every resolution element of our simulation. The solver accounts for CR acceleration at intra-cluster shocks, based on results of recent PIC simulations, re-acceleration due to shocks and MHD turbulence, adiabatic changes and radiative losses of electrons. We apply this model to zoom simulations of galaxy clusters, recently used to show the evolution of the small-scale turbulent dynamo on cluster scales. For these simulations we use a spectral resolution of 48 bins over 6 orders of magnitude in momentum for electrons and 12 bins over 6 orders of magnitude in momentum for protons. We present preliminary results about a possible formation mechanism for Wrong Way Radio Relics in our simulation.

Lunar impact flash (LIF) observations typically occur in R, I, or unfiltered light, and are only possible during night, targeting the night side of a 10-60% illumination Moon, while >10{\deg} above the observers horizon. This severely limits the potential to observe, and therefore the number of lower occurrence, high energy impacts observed is reduced. By shifting from the typically used wavelengths to the J-Band Short-Wave Infrared, the greater spectral radiance for the most common temperature (2750 K) of LIFs and darker skies at these wavelengths enables LIF monitoring to occur during the daytime, and at greater lunar illumination phases than currently possible. Using a 40.0 cm f/4.5 Newtonian reflector with Ninox 640SU camera and J-band filter, we observed several stars and lunar nightside at various times to assess the theoretical limits of the system. We then performed LIF observations during both day and night to maximise the chances of observing a confirmed LIF to verify the methods. We detected 61 >5{\sigma} events, from which 33 candidate LIF events could not be discounted as false positives. One event was confirmed by multi-frame detection, and by independent observers observing in visible light. While this LIF was observed during the night, the observed signal can be used to calculate the equivalent Signal-to-Noise ratio for a similar daytime event. The threshold for daylight LIF detection was found to be between Jmag=+3.4+-0.18 and Jmag=+5.6+-0.18 (Vmag=+4.5 and Vmag=+6.7 respectively at 2750 K). This represents an increase in opportunity to observe LIFs by almost 500%.

The spatial fluctuations in the tomographic maps of the redshifted 21-cm signal from the Cosmic Dawn (CD) crucially depend on the size and distribution of the regions with gas temperatures larger than the radio background temperature. In this article, we study the morphological characteristics of such emission regions and their absorption counterparts using the shape diagnostic tool {\sc surfgen2}. Using simulated CD brightness temperature cubes of the 21-cm signal, we find that the emission regions percolate at stages with the filling factor of the emission regions $FF_{\rm emi}\gtrsim 0.15$. Percolation of the absorption regions occurs for $FF_{\rm abs}\gtrsim 0.05$. The largest emission and absorption regions are topologically complex and highly filamentary for most parts of the CD. The number density of these regions as a function of the volume shows the power-law nature with the power-law indexes $\approx -2$ and $-1.6$ for the emission and absorption regions, respectively. Overall, the planarity, filamentarity and genus increase with the increase of the volume of both emission and absorption regions.

We present the results of observations of complex powerful type II burst associated with narrow Earth-directed CME, which was ejected at around 11 UT on 31 May 2013. The observations were performed by radio telescope UTR-2, which operated as local interferometer, providing the possibility of detection of the spatial parameters of the radio emission source. There are also polarization data from URAN-2 radio telescope. The CME was detected by two space-born coronagraphs SOHO/LASCO/C2 and STEREO/COR1-BEHIND, and was absolutely invisible for STEREO-AHEAD spacecraft. The associated type II burst consisted of two successive parts of quite different appearance on the dynamic spectrum. The first burst was narrow in frequency, had cloudy structure and was completely unpolarized while the second one represented rich herring-bone structure and exposed high degree of circular polarization. Both parts of the whole event reveal band splitting and well distinguished harmonic structure. The positions and sizes of the sources of the type II burst were found using cross-correlation functions of interferometer bases. The sources of the type II bursts elements were found to be of about 15 arcmin in size in average, with the smallest ones reaching as low as 10 arcmin. Corresponding brightness temperatures were estimated. In most cases these temperatures were between $10^{11}$ and $10^{12}$ K with maximum value as high as $10^{14}$ K. The spatial displacement of the source was measured and model independent velocities of the type II burst sources were determined.

We propose a new method for finding galaxy protoclusters that is motivated by structure formation theory, and is also directly applicable to observations. Protoclusters are defined as the galaxy groups whose virial mass $M_{\rm vir} < 10^{14}\,M_{\odot}$ at their epochs but would exceed that limit by $z=0$. They are distinguished from clusters, groups of galaxies whose virial mass currently exceeds $10^{14}\,M_{\odot}$. According to these definitions there can be a mixture of clusters and protoclusters at a given epoch. The future mass that a protocluster would acquire at $z=0$ is estimated using the spherical collapse model. The centers of protoclusters are identified using the critical overdensity for collapse by $z=0$ that is predicted by the spherical collapse model, and the physical size of protoclusters is defined by the overdensity corresponding to the turnaround radius. We use the cosmological hydrodynamical simulation Horizon Run 5 (HR5) to calibrate this prescription and demonstrate its performance. We find that the protocluster identification method suggested in this study is quite successful. Its application to the high redshift HR5 galaxies shows a tight correlation between the mass within the protocluster regions identified according to the spherical collapse model and the final mass to be found within the cluster at $z=0$, meaning that the regions can be regarded as the bona-fide protoclusters with high reliability.

Earth-like planets orbiting M-dwarf stars, M-Earths, are currently the best targets to search for signatures of life. Life as we know it requires water. The habitability of M-Earths is jeopardized by water loss to space: high flux from young M-dwarf stars can drive the loss of 3--20 Earth oceans from otherwise habitable planets. We develop a 0-D box model for Earth-mass terrestrial exoplanets, orbiting within the habitable zone, which tracks water loss to space and exchange between reservoirs during an early surface magma ocean phase and the longer deep-water cycling phase. A key feature is the duration of the surface magma ocean, assumed concurrent with the runaway greenhouse. This timescale can discriminate between desiccated planets, planets with desiccated mantles but substantial surface water, and planets with significant water sequestered in the mantle. A longer-lived surface magma ocean helps M-Earths retain water: dissolution of water in the magma provides a barrier against significant loss to space during the earliest, most active stage of the host M-dwarf, depending on the water saturation limit of the magma. Although a short-lived basal magma ocean can be beneficial to surface habitability, a long-lived basal magma ocean may sequester significant water in the mantle at the detriment of surface habitability. We find that magma oceans and deep-water cycling can maintain or recover habitable surface conditions on Earth-like planets at the inner edge of the habitable zone around late M-dwarf stars -- these planets would otherwise be desiccated if they form with less than ${\sim}$10 terrestrial oceans of water.

Galaxy clusters are promising targets for indirect detection of dark matter thanks to the large dark matter content. Using 14 years of Fermi-LAT data from seven nearby galaxy clusters, we obtain constraints on the lifetime of decaying very heavy dark matter particles with masses ranging from $10^3$ GeV to $10^{16}$ GeV. We consider a variety of decaying channels and calculate prompt gamma rays and electrons/positrons from the dark matter. Furthermore, we take into account electromagnetic cascades induced by the primary gamma rays and electrons/positrons, and search for the resulting gamma-ray signals from the directions of the galaxy clusters. We adopt a Navarro-Frenk-White profile of the dark matter halos, and use the profile likelihood method to set lower limits on the dark matter lifetime at a 95% confidence level. Our results are competitive with those obtained through other gamma-ray observations of galaxy clusters and provide complementary constraints to existing indirect searches for decaying very heavy dark matter.

The ultra-short-period super-Earth 55 Cancri e has a measured radius of 1.8 Earth radii. Previous thermal phase curve observations suggest a strong temperature contrast between the dayside and nightside of around 1000 K with the hottest point shifted $41\pm12$ degrees east from the substellar point, indicating some degree of heat circulation. The dayside (and potentially even the nightside) is hot enough to harbour a magma ocean. We use results from general circulation models (GCMs) of atmospheres to constrain the surface temperature contrasts. There is still a large uncertainty on the vigour and style of mantle convection in super-Earths, especially those that experience stellar irradiation large enough to harbour a magma ocean. In this work, we aim to constrain the mantle dynamics of the tidally locked lava world 55 Cancri e. Using the surface temperature contrasts as boundary condition, we model the mantle flow of 55 Cancri e using 2D mantle convection simulations and investigate how the convection regimes are affected by the different climate models. We find that large super-plumes form on the dayside if that hemisphere is covered by a magma ocean and the nightside remains solid or only partially molten. Cold material descends into the deep interior on the nightside, but no strong downwellings form. In some cases, the super-plume also moves several tens of degrees towards the terminator. A convective regime where the upwelling is preferentially on the dayside might lead to preferential outgassing on that hemisphere which could lead to the build-up of atmospheric species that could be chemically distinct from the nightside.

Kilonovae are a novel class of astrophysical transients, and the only observationally-confirmed site of rapid neutron capture nucleosynthesis (the r-process) in the Universe. To date, only a handful of kilonovae have been detected, with just a single spectroscopically-observed event (AT 2017gfo). Spectra of AT 2017gfo provided evidence for the formation of elements heavier than iron; however, these spectra were collected during the first ~ 10 days, when emission from light r-process elements dominates the observations. Heavier elements, if synthesised, are expected to shape the late-time evolution of the kilonova, beyond the phases for which we have spectral observations. Here we present spectroscopic observations of a rapidly-reddening thermal transient, following the gamma-ray burst, GRB 230307A. Early (2.4 day) optical spectroscopy identifies the presence of a hot (T ~ 6700 K) thermal continuum. By 29 days, this component has expanded and cooled significantly (T ~ 640 K), yet it remains optically thick, indicating the presence of high-opacity ejecta. We show that these properties can only be explained by the merger of compact objects, and further, leads us to infer the production of the heavy lanthanide elements. We identify several spectral features (in both absorption and emission), whose cause can be explained by newly-synthesised heavy elements. This event marks only the second recorded spectroscopic evidence for the synthesis of r-process elements, and the first to be observed at such late times.

Kilonovae are a rare class of astrophysical transients powered by the radioactive decay of nuclei heavier than iron, synthesized in the merger of two compact objects. Over the first few days, the kilonova evolution is dominated by a large number of radioactive isotopes contributing to the heating rate. On timescales of weeks to months, its behavior is predicted to differ depending on the ejecta composition and merger remnant. However, late-time observations of known kilonovae are either missing or limited. Here we report observations of a luminous red transient with a quasi-thermal spectrum, following an unusual gamma-ray burst of long duration. We classify this thermal emission as a kilonova and track its evolution up to two months after the burst. At these late times, the recession of the photospheric radius and the rapidly-decaying bolometric luminosity ($L_{\rm bol}\propto t^{-2.7\pm 0.4}$) support the recombination of lanthanide-rich ejecta as they cool.

Detecting H2O in exoplanet atmospheres is the first step on the path to determining planet habitability. Coronagraphic design currently limits the observing strategy used to detect H2O, requiring the choice of specific bandpasses to optimize abundance constraints. In order to examing the optimal observing strategy for initial characterization of habitable planets using coronagraph-based direct imaging, we quantify the detectability of H2O as a function of signal-to-noise ratio (SNR) and molecular abundance across 25 bandpasses in the visible wavelength range (0.5-1 micron). We use a pre-constructed grid consisting of 1.4 million geometric albedo spectra across a range of abundance and pressure, and interpolate the produce forward models for an efficient nested sampling routine, PSGnest. We first test the detectability of H2O in atmospheres that mimix a modern-Earth twin, and then expand to examine a wider range of H2O abundances; for each abundance value, we constrain the optimal 20% bandpasses based on the effective signal-to-noise ratio (SNR) of the data. We present our findings of H2O detectability as functions of SNR, wavelength, and abundance, and discuss how to use these results for optimizing future coronographic instrument design. We find that there are specific points in wavelength where H2o can be detected down to 0.74 micron with moderate-SNR data for abundances at the upper end of Earth's presumed historical values, while at 0.9 micron, detectability is possible with low-SNR data at modern Earth abundances of H2O.

The chaotic nature of outer space and the limitation of visual displays command for much more than visual display of information and for the integration of other sensorial modalities during data exploration. Haptic real time devices may enrich the detection of astronomical events that otherwise would escape the eyes. Departing from the Harvard Astronomy Lab Orchestar (color Arduino) we present the work in progress of the Proof of Concept (PoC) of a sensitive yet simple device to Bluetooth transfer real time color into haptic motion built by Adafruit components. We assemble 2 Adafruit nRF52840 feather express with the RGB color sensor and haptic driver respectively to trigger vibrations according to the color variation from external light sources. In addition, a transparent hexagon cover will be mounted on the color sensor to maximise absorbed light from the telescope. The device aims to be a translator for people to see hear and feel the hidden information from the original data set. We also present its application in the calculation of complex astrophysics quantities such as the masses of solar coronal mass ejections.

The presence of blueshifted absorption lines in the X-ray spectra of Black Hole Low Mass X-ray Binaries is the telltale of massive outflows called winds. These signatures are found almost exclusively in soft states of high-inclined systems, hinting at equatorial ejections originating from the accretion disk and deeply intertwined with the evolution of the outburst patterns displayed by these systems. In the wake of the launch of the new generation of X-ray spectrometers, studies of wind signatures remain mostly restricted to single sources and outbursts, with some of the recent detections departing from the commonly expected behaviors. We thus give an update to the current state of iron band absorption lines detections, through the analysis of all publicly available XMM-$Newton$-PN and $Chandra$-HETG exposures of known Black Hole Low-Mass X-ray Binary candidates. Our results agree with previous studies, with wind detections exclusively found in dipping, high-inclined sources, and almost exclusively in bright ($L_{X}>0.01L_{Edd}$) soft ($HR<0.8$) states, with blueshift values generally restricted to few 100 km s$^{-1}$. The line parameters indicate similar properties between objects and outbursts of single sources, and despite more than 20 years of data, very few sources have the HID sampling necessary to properly study the evolution of the wind during single outbursts. We provide an online tool with details of the wind signatures and outburst evolution data for all sources in the sample.

Exoplanets can be detected with various observational techniques. Among them, radial velocity (RV) has the key advantages of revealing the architecture of planetary systems and measuring planetary mass and orbital eccentricities. RV observations are poised to play a key role in the detection and characterization of Earth twins. However, the detection of such small planets is not yet possible due to very complex, temporally correlated instrumental and astrophysical stochastic signals. Furthermore, exploring the large parameter space of RV models exhaustively and efficiently presents difficulties. In this review, we frame RV data analysis as a problem of detection and parameter estimation in unevenly sampled, multivariate time series. The objective of this review is two-fold: to introduce the motivation, methodological challenges, and numerical challenges of RV data analysis to nonspecialists, and to unify the existing advanced approaches in order to identify areas for improvement.

We study the $\Delta$ resonance inside neutron stars by considering relativistic mean fields containing strange mesons ($\sigma^*$ and $\phi$). $\sigma^*$ and $\phi$ shift the critical density of hyperons to a lower density region, while shifting the threshold of $\Delta$ resonance to a higher density region, in this respect an early appearance of $\Delta$ resonances is crucial to guarantee the stability of the branch of hyperonized star with the difference of the coupling parameter $x_{\sigma \Delta}$ constrained based on the QCD rules in nuclear matter. The $\Delta$ resonance produces a softer equation of state in the low density region, which makes the tidal deformability and radius consistent with the observation of GW170817. As the addition of new degrees of freedom will lead to a softening of the equation of state, the ${\sigma}$-cut scheme, which states the softening of EOS can be slowed down if one assumes a limited decrease of the ${\sigma}$-meson strength at ${\rho_B}$($\rho_B > \rho_0$), finally we get a maximum mass neutron star with $\Delta$ resonance heavier than 2$M_{\odot}$.

The Eisenhart lift allows to formulate the dynamics of a scalar field in a potential as pure geodesic motion in a curved field-space manifold involving an additional fictitious vector field. Making use of the formalism in the context of inflation, we show that the main inflationary observables can be expressed in terms of the geometrical properties of a two-dimensional uplifted field-space manifold spanned by the time derivatives of the scalar and the temporal component of the vector. This allows to abstract from specific potentials and models and describe inflation solely in terms of the flow of geometric quantities. Our findings are illustrated through several inflationary examples previously considered in the literature.

We sketch the main features of the Noether Symmetry Approach, a method to reduce and solve dynamics of physical systems by selecting Noether symmetries, which correspond to conserved quantities. Specifically, we take into account the vanishing Lie derivative condition for general canonical Lagrangians to select symmetries. Furthermore, we extend the prescription to the first prolongation of the Noether vector. It is possible to show that the latter application provides a general constraint on the infinitesimal generator $\xi$, related to the spacetime translations. This approach can be used for several applications. In the second part of the work, we consider a gravity theory including the coupling between a scalar field $\phi$ and the Gauss-Bonnet topological term $\mathcal{G}$. In particular, we study a gravitational action containing the function $F(\mathcal{G}, \phi)$ and select viable models by the existence of symmetries. Finally, we evaluate the selected models in a spatially-flat cosmological background and use symmetries to find exact solutions.

Cosmological gravitational wave backgrounds (CGWBs) are the conglomeration of unresolved gravitational wave signals from early Universe sources, which make them a promising tool for cosmologists. Because gravitons decouple from the cosmic plasma early on, one can consider interactions between gravitons and any particle species that were present in the very early Universe. Analogous to the cosmic microwave background (CMB), elastic scattering on any cosmological background will induce small distortions in its energy density spectrum. We then quantify the magnitude of these spin-dependent spectral distortions when attributed to the dark matter (DM) in the early Universe. Lastly, we give estimates for potentially measurable distortions on CGWBs due to gravitational scattering by primordial black holes.

The macroscopic properties of compact stars in modified gravity theories can be significantly different from the general relativistic (GR) predictions. Within the gravitational context of scalar-tensor theories, with a scalar field $\phi$ and coupling function $\Phi(\phi)= \exp[2\phi/\sqrt{3}]$, we investigate the hydrostatic equilibrium structure of neutron stars for the simple potential $V(\phi)= \omega\phi^2/2$ defined in the Einstein frame (EF). From the scalar field in the EF, we also interpret such theories as $f(R)$ gravity in the corresponding Jordan frame (JF). The mass-radius relations, proper mass, and binding energy are obtained for a polytropic equation of state (EoS) in the JF. Our results reveal that the maximum-mass values increase substantially as $\omega$ gets smaller, while the radius and mass decrease in the low-central-density region as we move further away from the pure GR scenario. Furthermore, a cusp is formed when the binding energy is plotted as a function of the proper mass, which indicates the appearance of instability. Specifically, we find that the central-density value where the binding energy is a minimum corresponds precisely to $dM/d\rho_c^J = 0$ on the $M(\rho_c^J)$-curve.

The correlation between the charge radii differences in mirror nuclei pairs and the neutron skin thickness has been studied with the so-called finite range simple effective interaction over a wide mass region. The so far precisely measured charge radii difference data within their experimental uncertainty ranges in the 34Ar-34S, 36Ca-36S, 38Ca-38Ar, and 54Ni-54Fe mirror pairs are used to ascertain an upper limit for the slope parameter of the nuclear symmetry energy L $\approx$ 100 MeV. This limiting value of L is found to be consistent with the upper bound of the NICER PSR J0740+6620 constraint at 1$\sigma$ level for the radius R$_{1.4}$ of 1.4 M$_\odot$ neutron stars. The lower bound of the NICER R$_{1.4}$ data constrains the lower limit of L to $\approx$ 70 MeV. Within the range for L = 70-100 MeV the tidal deformability $\Lambda^{1.4}$ constraint, which is extracted from the GW170817 event at 2$\sigma$ level, and the recent PREX-2 and CREX data on the neutron skin thickness are discussed.

The gravitational constant variation means the breakdown of the strong equivalence principle. As the cornerstone of general relativity, the validity of general relativity can be examined by studying the gravitational constant variation. Such variations have the potential to affect both the generation and propagation of gravitational waves. In this paper, our focus lies on the effect of gravitational constant variation specifically on the propagation of gravitational waves. We employ two analytical methods, namely based on the Fierz-Pauli action and the perturbation of Einstein-Hilbert action around Minkowski spacetime, both leading to the the same gravitational wave equation. By solving this equation, we find the effects of gravitational constant variation on gravitational wave propagation. The result is consistent with previous investigations based on Maxwell-like equations for gravitational waves. Notably, we find that small variations in the gravitational constant result in an amplitude correction at the leading order and a phase correction at the sub-leading order for gravitational waves. These results provide valuable insights for probing gravitational constant variation and can be directly applied to gravitational wave data analysis.

Approximate all-terrain spacetimes for astrophysical applications are presented. The metrics possess five relativistic multipole moments, namely mass, rotation, mass quadrupole, charge, and magnetic dipole moment. All these spacetimes approximately satisfy the Einstein-Maxwell field equations. The first metric is generated by means of the Hoenselaers-Perjes method from the given relativistic multipoles. The second metric is a perturbation of the Kerr-Newman metric, which makes it a relevant approximation for astrophysical calculations. The last metric is an extension of the Hartle-Thorne metric that is important for obtaining internal models of compact objects perturbatively. These spacetimes are relevant to infer properties of compact objects from astrophysical observations. Furthermore, the numerical implementations of these metrics are straightforward, making them versatile for simulating these potential astrophysical applications.

We show that warm inflation can be successfully realized in the high temperature regime through dissipative interactions between the inflaton and a single fermionic degree of freedom, provided that the latter's mass is an oscillatory function of the inflaton field value. We demonstrate, in particular, that despite the consequent large amplitude oscillations of the eta slow-roll parameter, their effect is, on average, sufficiently suppressed to allow for a slow-roll trajectory. In addition, we demonstrate that, even though this also induces a parametric resonance that amplifies inflaton perturbations, this has a negligible effect on CMB scales in the relevant parametric range. Hence, the "Warm Little Inflaton" scenario can be realized with one less fermionic degree of freedom and no need of imposing an additional discrete interchange symmetry.

We study the cosmological implications of gravity models which break diffeomorphisms (Diff) invariance down to transverse diffeomorphisms (TDiff). We start from the most general gravitational action involving up to quadratic terms in derivatives of the metric tensor and identify TDiff models as the only stable theories consistent with local gravity tests. These models propagate an additional scalar graviton and although they are indistinguishable from GR at the post-Newtonian level, their cosmological dynamics exhibits a rich phenomenology. Thus we show that the model includes standard $\Lambda$CDM as a solution when the extra scalar mode is not excited, but different cosmological evolutions driven by the new term are possible. In particular, we show that for a soft Diff breaking, the new contribution always behaves as a cosmological constant at late times. When the extra contribution is not negligible, generically its evolution either behaves as dark energy or tracks the dominant background component. Depending on the initial conditions, solutions in which the universe evolves from an expanding to a contracting phase, eventually recollapsing are also possible.

The current observational properties of neutron stars have not definitively ruled out the possibility of dark matter. In this study, we primarily focus on exploring correlations between the dark matter model parameters and different neutron star properties using a rich set of EOSs. We adopt a two-fluid approach to calculate the properties of neutron stars. For the nuclear matter EOS, we employ several realistic EOS derived from the relativistic mean field model (RMF), each exhibiting varying stiffness and composition. In parallel, we look into the dark matter EOS, considering fermionic matter with repulsive interaction described by a relativistic mean field Lagrangian. A reasonable range of parameters is sampled meticulously. Interestingly, our results reveal a promising correlation between the dark matter model parameters and stellar properties, particularly when we ignore the uncertainties in the nuclear matter EOS. However, when introducing uncertainties in the nuclear sector, the correlation weakens, suggesting that the task of conclusively constraining any particular dark matter model might be challenging using global properties alone, such as mass, radius, and tidal deformability. Notably, we find that dark-matter admixed stars tend to have higher central baryonic density, potentially allowing for non-nucleonic degrees of freedom or direct Urca processes in stars with lower masses. There is also a tantalizing hint regarding the detection of stars with the same mass but different surface temperatures, which may indicate the presence of dark matter. With our robust and extensive dataset, we delve deeper and demonstrate that even in the presence of dark matter, the semi-universal C-Love relation remains intact.

Conversion of electromagnetic energy into magnetohydrodynamic energy occurs when the electric conductivity changes from negligible to finite values. This process is relevant during the epoch of reheating of the early Universe at the end of inflation and before the emergence of the radiation-dominated era. We find that conversion into kinetic and thermal energies is primarily the result of electric energy dissipation and that the magnetic energy plays only a secondary role in this process. This means that, since electric energy dominates over magnetic energy during inflation and reheating, significant amounts of electric energy can be converted into magnetohydrodynamic energy when conductivity emerges early enough, before the relevant length scales become stable.