Papers for Thursday, Sep 21 2023

We present observational evidence of the correlation between the high-mass slope of the stellar initial mass function (IMF) in young star clusters and their stellar surface density, $\sigma_{*}$. When the high-mass end of the IMF is described by a power law of the form $dN/d{\rm log}{M_{*}}\propto M_{*}^{-\Gamma}$, the value of $\Gamma$ is seen to weakly decrease with increasing $\sigma_{*}$, following a $\Gamma=1.31~\sigma_{*}^{-0.095}$ relation. We also present a model that can explain these observations. The model is based on the idea that the coalescence of protostellar cores in a protocluster forming clump is more efficient in high density environments where cores are more closely packed. The efficiency of the coalescence process is calculated as a function of the parental clump properties and in particular the relation between its mass and radius as well as its core formation efficiency. The main result of this model is that the increased efficiency of the coalescence process leads to shallower slopes of the IMF in agreement with the observations of young clusters, and the observations are best reproduced with compact protocluster forming clumps. These results have significant implications for the shape of the IMF in different Galactic and extragalactic environments and have very important consequences for galactic evolution.

Recent work on the characterization of small exoplanets has allowed us to accumulate growing evidence that the sub-Neptunes with radii greater than $\sim2.5\,R_\oplus$ often host H$_2$/He-dominated atmospheres both from measurements of their low bulk densities and direct detections of their low mean-molecular-mass atmospheres. However, the smaller sub-Neptunes in the 1.5-2.2 R$_\oplus$ size regime are much less understood, and often have bulk densities that can be explained either by the H$_2$/He-rich scenario, or by a volatile-dominated composition known as the "water world" scenario. Here, we report the detection of water vapor in the transmission spectrum of the $1.96\pm0.08$ R$_\oplus$ sub-Neptune GJ9827d obtained with the Hubble Space Telescope. We observed 11 HST/WFC3 transits of GJ9827d and find an absorption feature at 1.4$\mu$m in its transit spectrum, which is best explained (at 3.39$\sigma$) by the presence of water in GJ9827d's atmosphere. We further show that this feature cannot be caused by unnoculted star spots during the transits by combining an analysis of the K2 photometry and transit light-source effect retrievals. We reveal that the water absorption feature can be similarly well explained by a small amount of water vapor in a cloudy H$_2$/He atmosphere, or by a water vapor envelope on GJ9827d. Given that recent studies have inferred an important mass-loss rate ($>0.5\,$M$_\oplus$/Gyr) for GJ9827d making it unlikely to retain a H-dominated envelope, our findings highlight GJ9827d as a promising water world candidate that could host a volatile-dominated atmosphere. This water detection also makes GJ9827d the smallest exoplanet with an atmospheric molecular detection to date.

Galaxy cluster mergers are natural consequences of the structure formation in the Universe. Such events involve a large amount of energy ($\sim 10^{63}$ erg) dissipated during the process. Part of this energy can be channelled in particle acceleration and magnetic field amplification, enhancing non-thermal emission of the intra- and inter-cluster environment. Recently, low-frequency observations have detected a bridge of diffuse synchrotron emission connecting two merging galaxy clusters, Abell 399 and Abell 401. Such a result provides clear observational evidence of relativistic particles and magnetic fields in-between clusters. In this work, we have used LOw Frequency ARray (LOFAR) observations at 144 MHz to study for the first time the polarized emission in the A399-A401 bridge region. No polarized emission was detected from the bridge region. Assuming a model where polarization is generated by multiple shocks, depolarization can be due to Faraday dispersion in the foreground medium with respect to the shocks. We constrained its Faraday dispersion to be greater than 0.10 rad m$^{-2}$ at 95% confidence level, which corresponds to an average magnetic field of the bridge region greater than 0.46 nG (or 0.41 nG if we include regions of the Faraday spectrum that are contaminated by Galactic emission). This result is largely consistent with the predictions from numerical simulations for Mpc regions where the gas density is $\sim 300$ times larger than the mean gas density.

The strength and morphology of the Sun's magnetic field evolves significantly during the solar cycle, with the overall polarity of the Sun's magnetic field reversing during the maximum of solar activity. Long-term changes are also observed in sunspot and geomagnetic records, however systematic magnetic field observations are limited to the last four cycles. We investigate the long-term evolution of the Sun's magnetic field, and the influence this has on the topology and rotation of the solar corona. The Sun's photospheric magnetic field was decomposed into spherical harmonics using synoptic Carrington magnetograms from 1) WSO, 2) MDI onboard the SOHO, and 3) HMI onboard the SDO. The time-evolution of the spherical harmonic coefficients was used to explore the variation of the Sun's magnetic field, focusing on the large-scale modes. PFSS extrapolations of the photospheric field were computed to follow topological changes in the corona. The footpoints of the Sun's open magnetic field vary between the polar coronal holes and activity driven features such as active regions, and equatorial coronal holes. Consequently, the mean rotation rate of the solar wind is modulated during each cycle by the latitudinal variation of open field footpoints, with slower rotation during minima and faster (Carrington-like) rotation during maxima. Thisc variation is sensitive to cycle to cycle differences in the polar field strengths and hemispherical flux emergence rates, with the ratio of quadrupole to dipole energy following a similar variation. Cycle 23 maintained a larger fraction of quadrupolar energy in the declining phase, which kept the sources of open magnetic flux closer to the equator, extending the period of faster equator-ward connectivity. The ratio of quadrupole to dipole energy could be a useful proxy when examining the impact of differential rotation on the coronae of other Sun-like stars.

Using an updated and significantly augmented sample of Cepheid and TRGB distances to 28 nearby spiral and irregular galaxies, covering a wide range of metallicities, we have searched for evidence of a correlation of the zero-point of the Cepheid Period-Luminosity relation with HII region (gas-phase) metallicities. Our analysis, for the 21 galaxies closer than 12.5 Mpc, results in the following conclusions: (1) The zero points of the Cepheid and TRGB distance scales are in remarkably good agreement, with the mean offset in the zero points of the most nearby distance-selected sample being close to zero, Delta mod_o(Cepheid - TRGB) = -0.026 +\- 0.015 mag (for an I-band TRGB zero point of M_I = -4.05 mag); however, for the more distant sample, there is a larger offset between the two distance scales, amounting to -0.073 +/- 0.057 mag. (2) The individual differences, about that mean, have a measured scatter of +/- 0.068~mag. (3) We find no statistically significant evidence for a metallicity dependence in the Cepheid distance scale using the reddening-free W(V,VI) period-luminosity relation: Delta mod_o (Cepheid - TRGB) = -0.022 (+/- 0.015) \times ([O/H]-8.50) - 0.003 (+/- 0.007)

Context. Open clusters are groups of coeval stars sharing properties such as distance and metallicity, and they are key to understanding stellar evolution. Aims. Our main goal is to study the evolution of open clusters with a special focus on the universality of the luminosity function. Methods. We applied an upgraded version of the convergent point technique on about 50 open clusters. The selection of cluster members was based purely on the exquisite astrometry of the Gaia DR3 and Hipparcos catalogues in the five-dimensional or full six-dimensional space. Results. We present updated lists of bona fide members of ~50 open clusters within 500 pc and younger than 1 Gyr, exploiting the full depth of the third Gaia data release complemented by Hipparcos at the bright end, excluding regions in the Galactic plane. Our catalogues also are complemented by optical and infrared photometry from the major large-scale public surveys. All the data will be made available on a dedicated webpage with interactive plots and a direct link to Aladin and Vizier hosted at the Centre de Donn\'ees de Strasbourg. We derived luminosity functions for all bound clusters and compared them in three age groups of ~50 Myr, ~150 Myr, and ~600 Myr, discussing similarities and differences to constrain their dynamical evolution. Conclusions. Luminosity functions of clusters at 50 Myr are more likely similar to each other and show a greater degree of similarity than older clusters. We explain this observation with the universal luminosity function within the volume of our sample (500 pc). Luminosity functions of clusters with ages similar to the Pleiades or Hyades are more diverse, perhaps due to internal dynamical evolution, but more work is needed to provide additional evidence.

The cosmic large-scale structure (LSS) provides a unique testing ground for connecting fundamental physics to astronomical observations. Modelling the LSS requires numerical $N$-body simulations or perturbative techniques that both come with distinct shortcomings. Here we present the first unified numerical approach, enabled by new time integration and discreteness reduction schemes, and demonstrate its convergence at the field level. In particular, we show that our simulations (1) can be initialised directly at time zero, and (2) can be made to agree with high-order Lagrangian perturbation theory in the fluid limit. This allows fast, self-consistent, and UV-complete forward modelling of LSS observables.

de Grijs and Bono compiled 211 independent measurements of the distance to galaxy M87 in the Virgo cluster from 15 different tracers and reported the arithmetic mean of a subset of this compilation as the best estimate of the distance. We compute three different central estimates -- the arithmetic mean, weighted mean, and the median -- and corresponding statistical uncertainty for the full data set as well as two sub-compilations. We find that for all three central estimates the error distributions show that the data sets are significantly non-Gaussian. As a result, we conclude that that the median is the most reliable of the three central estimates, as median statistics does not assume Gaussianity. We use median statistics to determine the systematic error on the distance by analyzing the scatter in the 15 tracer subgroup distances. From the 211 distance measurements, we recommend a summary M87 distance modulus of $31.08^{+0.04}_{-0.05}$ (statistical) ${}^{+0.04}_{-0.06}$ (systematic) mag, or combining the two errors in quadrature $31.08^{+0.06}_{-0.08}$ mag, rounded to $16.4^{+0.5}_{-0.6}$ Mpc, all at $68.27\%$ significance.

We investigate sharply outlined features recorded in solar magnetic field tracers. It is shown that the magnetic boundaries of a sunspot do not coincide with the photometric ones. Moreover, there is no clear magnetic boundary around sunspots. Thus, the widely accepted concept of a magnetic tube with clearly pronounced borders is not always correct and should be used with caution. It is also shown that even in the periods of complete absence of visible spots on the Sun, there are magnetic fields over 800 Gauss. The nature of these strong magnetic fields remains unclear; they may originate at relatively small depths under the photosphere.

Giant planet migration appears widespread among planetary systems in our Galaxy. However, the timescales of this process, which reflect the underlying dynamical mechanisms, are not well constrained, even within the solar system. Since planetary migration scatters smaller bodies onto intersecting orbits, it would have resulted in an epoch of enhanced bombardment in the solar system's asteroid belt. To accurately and precisely quantify the timescales of migration, we interrogate thermochronologic data from asteroidal meteorites, which record the thermal imprint of energetic collisions. We present a database of 40K-40Ar system ages from chondrite meteorites and evaluate it with an asteroid-scale thermal code coupled to a Markov chain Monte Carlo inversion. Simulations require bombardment in order to reproduce the observed age distribution and identify a bombardment event beginning ~11 million years after the Sun formed. Our results associate a giant planet instability in our solar system with the dissipation of the gaseous protoplanetary disk.

The evolution of a planet's atmosphere depends strongly on its host star's properties. When their host stars are younger, planets can experience stronger winds and EUV emissions. This is particularly true for planets orbiting M-dwarfs due to their close proximity to the host star. To determine if these planets retain an atmosphere, we consider the impacts from stellar wind and EUV fluxes in driving atmospheric escape throughout the planet's lifetime. For this, we determined the atmospheric mass loss due to stellar wind and photoevaporation on 4 planets in close orbit and 34 in their star's habitable zone (HZ). The M-dwarf host stars' wind velocity, density, and EUV flux were calculated through rotation period and X-ray flux scaling over time. The mass loss rate due to stellar wind and photoevaporation was then computed as a function of time and accumulated throughout the planet's age to determine the total atmospheric mass loss of the planet's initial H/He envelope. We find that for HZ planets at orbits $<$ 0.1 AU, stellar wind can only remove $\leq 1\%$ of the H/He envelope, while photoevaporation is essential for completely removing the H/He envelope of most targets. Moreover, due to either mechanism, most planets orbiting at $>$ 0.1 AU do not have their primordial envelope stripped. Overall, out of the 38 planets studied, 13 were predicted to have lost the primordial envelope due to photoevaporation, while 2 planets lost the envelope due to both stellar wind and photoevaporation.

We recently constructed the G4Jy-3CRE, a catalog of extragalactic radio sources based on the GLEAM 4-Jy (G4Jy) sample, with the aim of increasing the number of powerful radio galaxies and quasars with similar selection criteria to those of the revised release of the Third Cambridge catalog (3CR). The G4Jy-3CRE consists of a total of 264 radio sources mainly visible from the Southern Hemisphere. Here, we present an initial X-ray analysis of 89 G4Jy-3CRE radio sources with archival X- ray observations from the Neil Gehrels Swift Observatory. We reduced a total of 615 Swift observations, for about 0.89 Msec of integrated exposure time, we found X-ray counterparts for 61 radio sources belonging to the G4Jy-3CRE, 11 of them showing extended X-ray emission. The remaining 28 sources do not show any X-ray emission associated with their radio cores. Our analysis demonstrates that X-ray snapshot observations, even if lacking uniform exposure times, as those carried out with Swift, allow us to (i) verify and/or re ne the host galaxy identi cation; (ii) discover the extended X-ray emission around radio galaxies of the intracluster medium when harbored in galaxy clusters, as the case of G4Jy 1518 and G4Jy 1664, and (iii) detect X-ray radiation arising from their radio lobes, as for G4Jy 1863.

VERITAS observed the bright blazar 1ES 2344+514 during two flaring periods, one from Dec. 17 to Dec. 18, 2015 (MJD 57373-57374) with a peak flux of ~60% of the Crab and another from Nov. 28 to Dec. 3, 2021 (MJD 59546-59551) with a peak flux of ~20% of the Crab. This blazar, located at a redshift of z = 0.044, is classified as an extreme high-frequency-peaked BL Lacertae object (HBL). It is known to be variable, including several previous day-scale flares: Whipple on Dec. 20, 1995, VERITAS on Dec. 7, 2007, and MAGIC on Aug. 11, 2016. The VERITAS near-nightly monitoring of 1ES 2344+514 during the 2015-2016 and 2021-2022 seasons provides good coverage of the pre- and post-flare flux as well as the rise/fall time of the flares. We present the multiwavelength light curves of each flare as well as the very high-energy spectra in the two flare states and the two pre-flare states.

We combine archival ALMA data targeting the Hubble Ultra Deep Field (HUDF) to produce the deepest currently attainable 1-mm maps of this key, extragalactic survey field. Combining all existing data in Band 6, our deepest map covers 4.2arcmin^2, with a beamsize of 1.49"x1.07" at an effective frequency of 243GHz (1.23mm). It reaches an rms of 4.6uJy/beam, with 1.5arcmin^2 below 9.0uJy/beam, an improvement of >5% over the best previously published map and 50% improvement in some regions. We also make a wider, but shallower map, covering 25.4arcmin^2. We detect 45 galaxies in the deep map down to 3.6sigma, including 10 more 1-mm sources than previously detected. 38 of these galaxies have a JWST ID from the JADES NIRCam imaging and the new sources are typically faint and red. A stacking analysis on the positions of ALMA-undetected JADES galaxies yields detections for z<4 and stellar masses from 10^(8.4) to 10^(10.4)Msun, extracting 10% of additional stacked signal from our map compared to previous analyses. Detected sources and stacking contribute (10.0+/-0.5)Jy/deg^2 of the cosmic infrared background (CIB) at 1.23mm. Although this is short of the (uncertain) background level of about 20Jy/deg^2, after taking into account intrinsic fluctuations in the CIB, our measurement is consistent with the background if the HUDF is a mild (~2sigma) negative fluctuation. This suggests that within the HUDF, JWST may have detected essentially all of the galaxies that contribute to the CIB. Our stacking analysis predicts that the field contains around 60 additional galaxies with 1.23mm flux densities averaging around 15uJy, and over 300 galaxies at the few uJy level. However, the contribution of these fainter more modestly-obscured objects to the background is small, and converging, as anticipated from the now well-established strong correlation between galaxy stellar mass and obscured star formation.

The exploration of the low acceleration $a<a_{0}$ regime, where $a_{0}=1.2 \times 10^{-10}$m s$^{-1}$ is the acceleration scale of MOND around which gravitational anomalies at galactic scale appear, has recently been extended to the much smaller mass and length scales of local wide binaries thanks to the availability of the {\it Gaia} catalogue. Statistical methods to test the underlying structure of gravity using large samples of such binary stars and dealing with the necessary presence of kinematic contaminants in such samples have also been presented. However, an alternative approach using binary samples carefully selected to avoid any such contaminants, and consequently much smaller samples, has been lacking a formal statistical development. In the interest of having independent high quality checks on the results of wide binary gravity tests, we here develop a formal statistical framework for treating small, clean, wide binary samples in the context of testing modifications to gravity of the form $G \to \gamma G$. The method is validated through extensive tests with synthetic data samples, and applied to recent {\it Gaia} DR3 binary star observational samples of relative velocities and internal separations on the plane of the sky, $v_{2D}$ and $r_{2D}$, respectively. Our final results for a high acceleration $r_{2D}<0.01$pc region are of $\gamma=1.000 \pm 0.096$, in full accordance with Newtonian expectations. For a low acceleration $r_{2D}>0.01$pc region however, we obtain $\gamma=1.512 \pm 0.199$, inconsistent with the Newtonian value of $\gamma=1$ at a $2.6 \sigma$ level, and much more indicative of MOND AQUAL predictions of close to $\gamma=1.4$.

Negative superhumps (NSHs) are signals a few percent shorter than the orbital period of a binary star and are considered to originate from the reverse precession of the tilted disk. Based on TESS photometry, we find nine new cataclysmic variable stars (CVs) with NSHs. Three (ASAS J1420, TZ Per, and V392 Hya) of these stars similar to AH Her still have NSHs during dwarf nova outbursts, and the NSH amplitude varies with the outburst. The variation in the radius of the accretion disk partially explains this phenomenon. However, it does not explain the rebound of the NSH amplitude after the peak of the outburst and the fact that the NSH amplitude of the quiescence is sometimes not the largest, and it is necessary to combine the disk instability model (DIM) and add other ingredients. Therefore, we suggest that the variation of NSH amplitude with outburst can be an essential basis for studying the origin of NSHs and improving the DIM. The six ( ASASSN-V J1137, ASASSN-V J0611, 2MASS J0715, LAMOST J0925, ASASSN-17qj, and ZTF18acakuxo) remaining stars have been poorly studied, and for the first time we determine their orbital periods, NSHs and Superorbital signal (SOR) periods. The NSH periods and amplitudes of ASASSN-V J1137 and ASASSN-17qj vary with the SOR, and based on the comparison of the observations with the theory, we suggest that a single change in tilted disk angle does not explain the observations of the SOR and that other ingredients need to be considered as well.

HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measurements using precision radial velocity monitoring calibrated with a suite of state-of-the-art absolute and relative wavelength references. MODHIS is the counterpart to HISPEC for the Thirty Meter Telescope and is being developed in parallel with similar scientific goals. In this proceeding, we provide a brief overview of the current design of both instruments, and the requirements for the two spectrographs as guided by the scientific goals for each. We then outline the current science case for HISPEC and MODHIS, with focuses on the science enabled for exoplanet discovery and characterization. We also provide updated sensitivity curves for both instruments, in terms of both signal-to-noise ratio and predicted radial velocity precision.

A phase-resolved analysis on the X-ray spectrum of Ultra-Luminous X-ray Pulsar (ULXP) NGC 300 ULX-1 is performed with data taken with XMM-Newton and NuSTAR on 2016 December 16th. In addition to the classical phase-restricting analysis, a method developed in active galactic nuclei studies is newly employed for ULXP. It has revealed that the pulsation cycle of the source can be divided into two intervals in terms of X-ray variability. This suggests the rotating flow consists of at least two representative emission regions. Furthermore, the new method successfully decomposed the spectrum into an independent pair in each interval. One is an unchanging-component spectrum that can be reproduced by a standard disk model with a $720^{+220}_{-120}$ km inner radius and a $0.25\pm0.03$ keV peak temperature. The other is the spectrum of the component that coincides with the pulsation. This was explained with a Comptonization of a $0.22^{+0.2}_{-0.1}$ keV blackbody and exhibited a harder photon index in the brighter phase interval of two. The results are consistent with a picture that the pulsating emission originates from a funnel-like flow formed within the magnetosphere, and the inner flow exhibiting a harder continuum is observed exclusively when the opening cone points to the observer.

Recently, it has been shown via the application of the Efron-Petrosian technique that the gamma-ray fluxes of pulsars scale with distance according to $F \propto D^{-3/2}$, thereby pointing to a violation of the inverse-square law. We apply the same procedure as in these works to the radio fluxes of pulsars detected in the Parkes multi-beam survey as well as for a synthetic population of pulsars whose fluxes scale with distance based on the inverse-square law. We find that both the observed data and the synthetic pulsar population show the same trends for the Efron-Petrosian statistics. Furthermore, lower values of the distance exponent are favored compared to the distance exponent of two (corresponding to inverse-square law) even for the synthetic pulsar population. Therefore, we conclude that the Efron-Petrosian method cannot be used to infer a violation of the inverse-square law for radio pulsar fluxes.

In this work, we investigate the optimal map-making technique for the linear system $\bm{d}=\bm{A}\bm{x}+\bm{n}$ while carefully taking into account singularities that may come from either the covariance matrix $\bm{C} = \<\bm{n}\bm{n}^t\>$ or the main matrix $\bm{A}$. We first describe the general optimal solution, which is quite complex, and then use the modified pseudo inverse to create a near-optimal solution, which is simple, robust, and can significantly alleviate the unwanted noise amplification during map-making. The effectiveness of the nearly optimal solution is then compared to that of the naive co-adding solution and the standard pseudo inverse solution, showing noticeable improvements. Interestingly, all one needs to get the near-optimal solution with singularity is just a tiny change to the traditional optimal solution that is designed for the case without singularity.

The understanding of the solar magnetic coronal structure is tightly linked to the shape of open field regions, specifically coronal holes. A dynamically evolving coronal hole coincides with the local restructuring of open to closed magnetic field, which leads to changes in the interplanetary solar wind structure. By investigating the dynamic evolution of a fast-tilting coronal hole, we strive to uncover clues about what processes may drive its morphological changes, which are clearly visible in EUV filtergrams. Using combined 193A and 195A EUV observations by AIA/SDO and EUVI/STEREO_A, in conjunction with line-of-sight magnetograms taken by HMI/SDO, we track and analyze a coronal hole over 12 days to derive changes in morphology, area and magnetic field. We complement this analysis by potential field source surface modeling to compute the open field structure of the coronal hole. We find that the coronal hole exhibits an apparent tilting motion over time that cannot solely be explained by solar differential rotation. It tilts at a mean rate of ~3.2{\deg}/day that accelerates up to ~5.4{\deg}/day. At the beginning of May, the area of the coronal hole decreases by more than a factor of three over four days (from ~13 * 10^9 km^2 to ~4 * 10^9 km^2), but its open flux remains constant (~2 * 10^20 Mx). Further, the observed evolution is not reproduced by modeling that assumes the coronal magnetic field to be potential. In this study, we present a solar coronal hole that tilts at a rate that has yet to be reported in literature. The rate exceeds the effect of the coronal hole being advected by either photospheric or coronal differential rotation. Based on the analysis we find it likely that this is due to morphological changes in the coronal hole boundary caused by ongoing interchange reconnection and the interaction with a newly emerging ephemeral region in its vicinity.

We report a first detection, at very high significance ($25\sigma$), of the cross-correlation between cosmic shear and the diffuse X-ray background, using data from the Dark Energy Survey and the ROSAT satellite. The X-ray cross-correlation signal is sensitive to the distribution of the surrounding gas in dark matter haloes. This allows us to use our measurements to place constraints on key physical parameters that determine the impact of baryonic effects in the matter power spectrum. In particular, we determine the mass of haloes in which feedback has expelled half of their gas content on average to be $\log_{10}(M_c/M_\odot)=13.643^{+0.081}_{-0.12}$, and the polytropic index of the gas to be $\Gamma = 1.231^{+0.015}_{-0.011}$. This represents a first step in the direct use of X-ray cross-correlations to obtain improved constraints on cosmology and the physics of the intergalactic gas.

We present a comprehensive spectro-temporal analysis of five ultraluminous X-ray sources (ULXs) with central object likely being a black hole, using archival {\it XMM-Newton} observations. These sources, namely NGC 1313 X-1, NGC 5408 X-1, NGC 6946 X-1, M82 X-1 and IC 342 X-1, reveal short-term variability with fractional variance of $1.42-27.28\%$ and exhibit QPOs with frequency $\nu_{\rm QPO} \sim 8-667$ mHz. Long-term evolution of ULXs energy spectra ($0.3 - 10$ keV; excluding M82 X$-$1) are described satisfactorily with a model combination that comprises a thermal Comptonization component (\texttt{nthComp}, yielding $\Gamma_{\rm nth} \sim 1.48 - 2.65$, $kT_{\rm e} \sim 1.62 - 3.76$ keV, $\tau \sim 8 - 20$, y-par $\sim 1.16 - 6.24$) along with a standard disc component (\texttt{diskbb}, $kT_{\rm in} \sim 0.16 - 0.54$ keV). We find that these ULXs generally demonstrate anti-correlation between disc luminosity and inner disc temperature as $L_{\rm disc} \propto T_{\rm in}^\alpha$, where $\alpha = - 3.58 \pm 0.04$ for NGC 1313 X-1 and IC 342 X-1, $\alpha = - 8.93 \pm 0.11$ for NGC 6946 X-1, and $\alpha = - 10.31 \pm 0.10$ for NGC 5408 X-1. We also obtain a linear correlation between bolometric luminosity $L_{\rm bol}$ and $\Gamma_{\rm nth}$ that indicates spectral softening of the sources when $L_{\rm bol}$ increases. We observe that in presence of QPO, Comptonized seed photon fraction varies in between $\sim 5 - 20 \%$, while the Comptonized flux contribution ($50 - 90\%$) dominates over disc flux. Utilizing $\nu_{\rm QPO}$ and $L_{\rm bol}$, we constrain ULXs mass by varying their spin ($a_{\rm k}$) and accretion rate ($\dot m$). We find that NGC 6946 X-1 and NGC 5408 X-1 seem to accrete at sub-Eddington accretion rate provided their central sources are rapidly rotating, whereas IC 342 X-1 and NGC 1313 X-1 can accrete in sub/super-Eddington limit irrespective to their spin values.

The sunspot penumbra is usually observed in the photosphere and it is of particular interest for its magneto-convection which seems to transport the heat from the top of the convection zone into the solar atmosphere. It is well known that the penumbra magnetic field extends in the upper layers of the solar atmosphere forming the so called super-penumbra. Thanks to the application of the Self Organizing Map technique to a spectral dataset containing monochromatic images acquired along the Ca II 854.2 nm and H$\alpha$ 656.28 nm lines, we were able to segment the penumbra and to measure the plasma velocity along the chromospheric portions of penumbral filaments. We found that the head, body and tail of penumbral filaments show vertical flows compatible with the persistence of the Evershed flow. Instead, the inverse Evershed flow has been observed only in the outer portion of the super-penumbra. We found that two opposite Evershed regimes work next to each other, without overlapping and both contribute to the downflow around sunspots. These results confirm the uncombed model of the sunspot penumbra and provide some hints that the downflow around sunspots may be ascribed to the magnetic field dragging the plasma down.

Hard X-/soft Gamma-ray astronomy (> 100 keV) is a crucial field for the study of important astrophysical phenomena such as the 511 keV positron annihilation line in the Galactic center region and its origin, gamma-ray bursts, soft gamma-ray repeaters, nuclear lines from SN explosions and more. However, several key questions in this field require sensitivity and angular resolution that are hardly achievable with present technology. A new generation of instruments suitable to focus hard X-/soft Gamma-rays is necessary to overcome the technological limitations of current direct-viewing telescopes. One solution is using Laue lenses based on Bragg's diffraction in a transmission configuration. To date, this technology is in an advanced stage of development and further efforts are being made in order to significantly increase its technology readiness level (TRL). To this end, massive production of suitable crystals is required, as well as an improvement of the capability of their alignment. Such a technological improvement could be exploited in stratospheric balloon experiments and, ultimately, in space missions with a telescope of about 20 m focal length, capable of focusing over a broad energy pass-band. We present the latest technological developments of the TRILL (Technological Readiness Increase for Laue Lenses) project, supported by ASI, devoted to the advancement of the technological readiness of Laue lenses. We show the method we developed for preparing suitable bent Germanium and Silicon crystals and the latest advancements in crystals alignment technology.

A new detection system for X-/Gamma-ray broad energy passband detectors for astronomy has been developed. This system is based on Silicon Drift Detectors (SDDs) coupled with scintillator bars; the SDDs act as a direct detector of soft (<30 keV) X-ray photons, while hard X-/Gamma-rays are stopped by the scintillator bars and the scintillation light is collected by the SDDs. With this configuration, it is possible to build compact, position sensitive detectors with unprecedented energy passband (2 keV - 10/20 MeV). The X and Gamma-ray Imaging Spectrometer (XGIS) on board the THESEUS mission, selected for Phase 0 study for M7, exploits this innovative detection system. The Wide Field Monitor - Imager and Spectrometer (WFM-IS) of the ASTENA (Advanced Surveyor of Transient Events and Nuclear Astrophysics) mission concept consists of 12 independent detection units, also based on this new technology. For the WFM-IS, a coded mask provides imaging capabilities up to 150 keV, while above this limit the instrument will act as a full sky spectrometer. However, it is possible to extend imaging capabilities above this limit by alternatively exploiting the Compton kinematics reconstruction or by using the information from the relative fluxes measured by the different cameras. In this work, we present the instrument design and results from MEGAlib simulations aimed at evaluating the effective area and the imaging performances of the WFM-IS above 150 keV.

Shocks are ubiquitous in astrophysical sources, many of which involve relativistic bulk motions, leading to the formation of relativistic shocks. Such relativistic shocks have so far been studied mainly in one dimension, for simplicity, but the complex nature of the relevant astrophysical flows often requires higher dimensional studies. Here we study the two-dimensional problem of the reflection of a planer shock off of a wall for a general incidence angle and a cold unshocked medium. We use primarily relativistic hydrodynamic numerical simulations, and elaborately compare the results to an analytic treatment. The simulations are performed both in the rest frame S of the unshocked fluid, where the dimensionless proper speed of the singly shocked fluid is $u_1=\Gamma_1\beta_1$ and the shock incidence angle is $\alpha_1$, and in the rest frame S$^\prime$ of the point P of intersection of the incident shock and the wall for regular reflection (RR). Good agreement is obtained between the simulations in these two frames and with the analytic solution. The establishment of a steady flow in frame S$^\prime$ is explored, along with the transition between the strong and weak shock RR solutions. The transition line between RR and Mach reflection (MR) is studied numerically in the $u_1$-$\alpha_1$ plane and found to coincide with the analytic detachment/sonic line. The flow properties along the sonic line are investigated in detail focusing on how they vary between the Newtonian and relativistic limits.

We employ the corrected Gaia Early Data Release 3 (EDR3) photometric data and spectroscopic data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR7 to assemble a sample of approximately 0.25 million FGK dwarf photometric standard stars for the 12 J-PLUS filters using the Stellar Color Regression (SCR) method. We then independently validated the J-PLUS DR3 photometry, and uncovered significant systematic errors: up to 15 mmag in the results of Stellar Locus (SL) method, and up to 10 mmag mainly caused by magnitude-, color-, and extinction-dependent errors of the Gaia XP spectra with the Gaia BP/RP (XP) Synthetic Photometry (XPSP) method. We have also further developed the XPSP method using the corrected Gaia XP spectra by Huang et al. (2023) and applied it to the J-PLUS DR3 photometry. This resulted in an agreement of 1-5 mmag with the SCR method, and a two-fold improvement in the J-PLUS zero-point precision. Finally, the zero-point calibration for around 91% of the tiles within the LAMOST observation footprint is determined through the SCR method, with the remaining approximately 9% of tiles outside this footprint relying on the improved XPSP method. The re-calibrated J-PLUS DR3 photometric data establishes a solid data foundation for conducting research that depends on high-precision photometric calibration.

Small metal-containing molecules have been detected and recognized as one of the hybrid species efficiently formed in space; especially in the circumstellar envelopes of evolved stars. It has been predicted also that more complex hybrid species like those formed by metals and fullerenes (metallofullerenes) could be present in such circumstellar environments. Recently, quantum-chemical simulations of metallofullerenes have shown that they are potential emitters contributing to the observed mid-IR spectra in the fullerene-rich circumstellar environments of different types of evolved stars. Here we present the individual simulated mid-IR (~5-50 um) spectra of twenty-eight metallofullerene species; both neutral and charged endo- and exohedral metallofullerenes for seven different metals (Li, Na, K, Ca, Mg, Ti, and Fe) have been considered. The changes induced by the metal-C60 interaction on the intensity and position of the spectral features are highlighted using charge density difference maps and electron density partitioning. Our calculations identify the fundamental IR spectral regions where, depending on the metal binding nature, there should be a major spectral contribution from each of the metallofullerenes. The metallofullerenes IR spectra are made publicly available to the astronomical community, especially James Webb Space Telescope users, for comparisons that could eventually lead to the detection of these species in space.

We present the results of our search for highly variable active galactic nuclei (AGNs) the X-ray flux from which changed by more than an order of magnitude during the SRG/eROSITA all-sky survey. Using the eROSITA data obtained in the period from December 2019 to February 2022, we have found 1325 sources the X-ray flux from which in the 0.3-2.3 keV energy band changed by more than a factor of 10 at a confidence level of at least 99.73 %. Of them, 635 objects have been classified as AGNs or AGN candidates. We describe the procedure of searching for highly variable sources and the selection of extragalactic objects among them and describe the statistical properties of the produced catalog. We provide a catalog of 49 sources for which a statistically significant flux in their low state was detected. For the latter we provide their light curves and X-ray spectra and discuss in detail the most interesting of them.

We investigate cosmological consequences of a generalized early dark energy (EDE) model where a scalar field behaves as dark energy at various cosmological epochs for a broad range of parameters such as the energy scale and the initial field value. We consider power-law and axion-type potentials for such an EDE field and study how it affects the cosmological evolution. We show that gravitational wave background can be significantly enhanced to be detected in future observations such as LISA and DECIGO in some parameter space. Implications of the EDE model are also discussed for a scenario where a blue-tilted inflationary tensor power spectrum can explain the recent NANOGrav 15-year signal. We argue that the bounds on the reheating temperature can be relaxed compared to the case of the standard thermal history.

Our study focuses on analysing the coronal, transition and chromospheric activity of four rapidly rotating stars located within 50 pc in the solar neighbourhood. We have used the multi-wavelength capabilities of AstroSat, to investigate the outer atmospheres of AB Dor, BO Mic, DG CVn and GJ 3331. These stars, classified as M and K type active stars, are known for their short rotation periods, leading to increased surface magnetic activity. Our soft X-ray observations provide the coronal properties such as emission measures, temperatures and elemental coronal abundances. We report the detection of X-ray flares from AB Dor, BO Mic, and DG CVn, while UV light curves reveal flares in both BO Mic and DG CVn.

A long-standing question related to nova eruptions is how these eruptions can lead to the formation of dust despite the ostensibly inhospitable environment for dust within the hot, irradiated ejecta. Novae in systems such as the symbiotic binary RS Oph offers a articularly clear view of some nova shocks and any associated dust production. Here we use spectropolarimetric monitoring of the RS Oph starting two days after its eruption in 2021 Aug. to show that: dust was present in the RS Oph system as early as two days into the 2021 eruption; the spatial distribution of this early dust was asymmetric, with components both aligned with and perpendicular to the orbital plane of the binary; between two and nine days after the start of the eruption, this early dust was gradually destroyed; and dust was again created, aligned roughly with the orbital plane of the binary, more than 80 days after the start of the outburst, most likely as a result of shocks that arose as the ejecta interacted with circumbinary material concentrated in the orbital plane. Modelling of X-rays and very-high energy GeV and TeV emission from RS Oph days to months into the 2021 eruption suggests that collisions between the ejecta and the circumbinary material may have led to shock formation in two regions: the polar - perpendicular to the orbital plane where collimated outflows have been observed after prior eruptions, and a circumbinary torus in the orbital plane. The observations described here indicate that dust formed in approximately the same two regions, supporting the connection between shocks and dust in novae and revealing a very early onset of asymmetry. The spectropolarimetric signatures of RS Oph in the first week into the 2021 outburst indicate: polarized flux across the H{\alpha} emission line and position angle orientation relative to the radio axis are similar to the spectropolarimetric signatures of AGNs.

The reflection of shock waves from a wall or other obstacles has been extensively studied in the Newtonian limit. In this limit the flow can always be analysed in a frame $S'$ where it is in steady state, which greatly simplifies the problem. Consider a planar incident shock front moving at a velocity $v_{s1}$ along its normal and making an angle $\alpha_1$ relative to a planar reflecting wall, in the lab frame $S$ where the unshocked fluid is at rest. The point $P$ of its intersection with the wall moves along the wall at a velocity $v_p=v_{s1}/\sin\alpha_1$, and a boost along the wall at $v_p$ transforms to the steady-state rest frame $S'$. However, there is a "super-luminal" regime where $v_p>c$ and therefore such a steady-state rest frame does not exist (as well as in the limiting case $v_p=c$). In the Newtonian regime it corresponds to very small incidence angles, $\alpha_1\approx\sin\alpha_1<\beta_{s1}=v_{s1}/c\ll1$, and is therefore usually ignored. However, in the relativistic regime (where $\Gamma_{s1}=(1-\beta_{s1}^2)^{-1/2}\gg1$), which is of interest in astrophysics, it corresponds to almost all incidence angles, $\tan\alpha_1<\Gamma_{s1}\beta_{s1}\approx\Gamma_{s1}$. Therefore, a new approach is needed to solve this problem. We formulate the integral conservation laws in the lab frame $S$ for the case of regular reflection (RR), either a. for a fixed volume, or b. for a fixed fluid, showing the equivalence between them and with the steady-state oblique shock jump conditions in frame $S'$ in the sub-luminal regime ($v_p<c$). Using this formalism we find that both the weak and strong shock RR solutions are bounded by the detachment line on the higher incidence angles side, while the strong shock solution is also bounded by the luminal line on the lower incidence angles side, and exists only between these two critical lines, in the sub-luminal attachment region.

We present a sample of SRG/eROSITA X-ray sources located in the eastern Galactic hemisphere (0<l<180 deg), with significant proper motions according to GAIA eDR3 measurements and whose extragalactic nature has been confirmed. The catalog consists of 248 extragalactic sources with spectroscopically measured redshifts. It includes all objects available in the Simbad database and matched to the identified optical component within a radius of 0.5 arcsec. Additionally, the catalog includes 18 sources with the spectral redshift measurements based on observations at the Russian-Turkish 1.5-m telescope RTT-150. The sources of the catalog are AGNs of various types (Sy1, Sy2, LINER), quasars, radio galaxies, and star-forming galaxies. The imitation of significant proper motions can be explained (previously known in astrometry as the VIM effect) by the presence of transient events on the line of sight in the field of view of AGN nuclei and quasars (within the GAIA resolution element). Such astrophysical phenomena may be the supernovae outbursts, tidal destruction events in AGNs with double nuclei, variability of large-mass supergiants, the presence of O-B associations in field of view of variable brightness AGN, etc. A model of flares with a fast rise and exponential decay profile allows to describe the variable positional parameters of most similar sources observed in GAIA. This cross-matching approach of the X-ray source catalogs of the SRG/eROSITA observatory and the optical catalog of the GAIA observatory can be used as an independent technique for detecting transient events in the neighborhood of AGN core (on scales of several hundred parsecs in the picture plane).

The next-generation ground-based gamma-ray Cherenkov Telescope Array Observatory (CTAO) will consist of imaging atmospheric Cherenkov telescopes (IACTs) of three different sizes distributed in two sites. The Large-Sized Telescopes will cover the low-energy end of the CTA energy range, starting at about 20 GeV. After its first years of operation at the CTA northern site, the Large-Sized Telescope prototype (LST-1) is in the final stage of its commissioning phase, having collected a significant amount of scientific data to date. In this contribution, we present the physics performance of the telescope using low-zenith Crab Nebula observations and Monte Carlo simulations fine-tuned accordingly. We show performance figures of merit such as the energy threshold, effective area, energy and angular resolution, and sensitivity based on the standard Hillas-parameters approach and following the source-independent and dependent analysis methods. The analysis threshold is estimated at 30 GeV. The energy resolution is around 30%, and the angular resolution is 0.3 degrees at 100 GeV. The best integral sensitivity of LST-1 is about 1.1% of the Crab Nebula flux above 250 GeV for 50 hours of observations. We also show the spectral energy distribution and light curve from Crab Nebula observations, which agree with results from other IACTs and link smoothly with Fermi-LAT when considering statistical and systematic uncertainties near the energy threshold.

We investigate the possible anisotropy of the universe using the most up-to-date type Ia supernovae, i.e. the Pantheon+ compilation. We fit the full Pantheon+ data with the dipole-modulated $\Lambda$CDM model, and find that it is well consistent with a null dipole. We further divide the full sample into several subsamples with different high-redshift cutoff $z_c$. It is shown that the dipole appears at $2\sigma$ confidence level only if $z_c\leq 0.1$, and in this redshift region the dipole is very stable, almost independent of the specific value of $z_c$. For $z_c=0.1$, the dipole amplitude is $D=1.0_{-0.4}^{+0.4}\times 10^{-3}$, pointing towards $(l,b)=(334.5_{\ -21.6^{\circ}}^{\circ +25.7^{\circ}},16.0_{\ -16.8^{\circ}}^{\circ +27.1^{\circ}})$, which is about $65^{\circ}$ away from the CMB dipole. This implies that the full Pantheon+ is consistent with a large-scale isotropic universe, but the low-redshift anisotropy couldn't be purely explained by the peculiar motion of the local universe.

Because early black holes (BHs) grew to $\sim10^{9} ~M_\odot$ in less than 1 Gyr of cosmic time, BH seeding models face stringent constraints. To efficiently constrain the parameter space of possible seeding criteria, we combine the advantages of the cosmological IllustrisTNG (TNG) simulations with the flexibility of semi-analytic modeling. We identify TNG galaxies as BH seeding sites based on various criteria including a minimum gas mass of $10^7$-$10^9~M_\odot$, total host mass of $10^{8.5}$-$10^{10.5}~M_\odot$, and a maximum gas metallicity of $0.01 - 0.1 ~Z_\odot$. Each potential host is assigned a BH seed with a probability of $0.01 - 1$; these BHs are then traced through the TNG galaxy merger tree. This approach improves upon the predictive power of the simple TNG BH seeding prescription, especially in the low-mass regime at high redshift, and it is readily adaptable to other cosmological simulations. Most of our seed models predict $z\lesssim4$ BH mass densities that are consistent with empirical data as well as the TNG BHs. However, high-redshift BH number densities can differ by factors of $\sim$ 10 - 100 between models. In most models, $\lesssim10^5~M_\odot$ BHs substantially outnumber heavier BHs at high redshifts. Mergers between such BHs are prime targets for gravitational-wave detection with LISA. The $z=0$ BH mass densities in most models agree well with observations, but our strictest seeding criteria fail at high redshift. Our findings strongly motivate the need for better empirical constraints on high-$z$ BHs, and they underscore the significance of recent AGN discoveries with JWST.

We use the distance sum rule (DSR) method to constrain the spatial curvature of the Universe with a large sample of 161 strong gravitational lensing (SGL) systems, whose distances are calibrated from the Pantheon compilation of type Ia supernovae (SNe Ia) using deep learning. To investigate the possible influence of mass model of the lens galaxy on constraining the curvature parameter $\Omega_k$, we consider three different lens models. Results show that a flat Universe is supported in the singular isothermal sphere (SIS) model with the parameter $\Omega_k=0.049^{+0.147}_{-0.125}$. While in the power-law (PL) model, a closed Universe is preferred at $\sim 3\sigma$ confidence level, with the parameter $\Omega_k=-0.245^{+0.075}_{-0.071}$. In extended power-law (EPL) model, the 95$\%$ confidence level upper limit of $\Omega_k$ is $<0.011$. As for the parameters of the lens models, constrains on the three models indicate that the mass profile of the lens galaxy could not be simply described by the standard SIS model.

GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverage of the skymap beginning within 1 hour and 3 hours respectively of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 190425 has been suggested, given their spatial and temporal coincidence. The smaller sky localisation area of FRB 190425 and its dispersion measure have led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 +/- 10 Mpc. Our optical imaging covered the galaxy 6.0 hrs after GW190425 was detected and 3.5 hrs after the FRB 190425. No optical emission was detected and further imaging at +1.2 and +13.2 days also revealed no emission. If the FRB 190425 and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 190425 from a BNS merger involves a supramassive magnetised neutron star spinning down by dipole emission on the timescale of hours. We show that magnetar enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB-GW link for this source.

The Cherenkov Telescope Array (CTA) will be the next-generation observatory in the field of very-high-energy (20 GeV to 300 TeV) gamma-ray astroparticle physics. The traditional approach to data analysis in this field is to apply quality cuts, optimized using Monte Carlo simulations, on the data acquired to maximize sensitivity. Subsequent steps of the analysis typically use the surviving events to calculate one set of instrument response functions (IRFs) to physically interpret the results. However, an alternative approach is the use of event types, as implemented in experiments such as the Fermi-LAT. This approach divides events into sub-samples based on their reconstruction quality, and a set of IRFs is calculated for each sub-sample. The sub-samples are then combined in a joint analysis, treating them as independent observations. In previous works we demonstrated that event types, classified using Machine Learning methods according to their expected angular reconstruction quality, have the potential to significantly improve the CTA angular and energy resolution of a point-like source analysis. Now, we validated the production of event-type wise full-enclosure IRFs, ready to be used with science tools (such as Gammapy and ctools). We will report on the impact of using such an event-type classification on CTA high-level performance, compared to the traditional procedure.

In the last years there has been a substantial increase in the number of the reported massive and luminous star-forming regions with related explosive outflows thanks to the superb sensitivity and angular resolution provided by the new radio, infrared, and optical facilities. Here, we report one more explosive outflow related with the massive and bright star-forming region IRAS 12326$-$6245 using Band 6 sensitive and high angular resolution ($\sim$0.2$"$) Atacama Large Millimeter/Submillimeter Array (ALMA) observations. We find over 10 molecular and collimated well-defined streamers, with Hubble-Lemaitre like expansion motions, and pointing right to the center of a dusty and molecular shell (reported for the first time here) localized in the northern part of the UCHII region known as G301.1A. The estimated kinematic age, and energy for the explosion are $\sim$700 yrs, and 10$^{48}$ erg, respectively. Taking into account the recently reported explosive outflows together with IRAS 12326$-$6245, we estimate an event rate of once every 90 yr in our Galaxy, similar to the formation rate of massive stars.

We report the discovery of a hot Jupiter on a 3.28-day orbit around a 1.08 M$_{Sun}$ G0 star that is the secondary component in a loose binary system. Based on follow-up radial velocity observations of TOI-858 B with CORALIE on the Swiss 1.2 m telescope and CHIRON on the 1.5 m telescope at the Cerro Tololo Inter-American Observatory (CTIO), we measured the planet mass to be $1.10\pm 0.08$ M$_{J}$ . Two transits were further observed with CORALIE to determine the alignment of TOI-858 B b with respect to its host star. Analysis of the Rossiter-McLaughlin signal from the planet shows that the sky-projected obliquity is $\lambda = 99.3\pm 3.8$. Numerical simulations show that the neighbour star TOI-858 A is too distant to have trapped the planet in a Kozai-Lidov resonance, suggesting a different dynamical evolution or a primordial origin to explain this misalignment. The 1.15 Msun primary F9 star of the system (TYC 8501-01597-1, at $\rho$ ~11") was also observed with CORALIE in order to provide upper limits for the presence of a planetary companion orbiting that star.

The stochastic gravitational wave background (SGWB) recently detected by the PTA collaborations could be the gravitational waves (GWs) induced by curvature perturbations. However, primordial black holes (PBHs) might be overproduced if the SGWB is explained by the GWs induced by the curvature perturbations that follow the Gaussian distribution. This motivates models associated with the non-Gaussianity of the curvature perturbations that suppress the PBH production rate. In this work, we show that the axion curvaton model can produce the curvature perturbations that induce GWs for the detected SGWB while preventing the PBH overproduction with the non-Gaussianity.

Most stars form in highly clustered environments within molecular clouds, but eventually disperse into the distributed stellar field population. Exactly how the stellar distribution evolves from the embedded stage into gas-free associations and (bound) clusters is poorly understood. We investigate the long-term evolution of stars formed in the STARFORGE simulation suite -- a set of radiation-magnetohydrodynamic simulations of star-forming turbulent clouds that include all key stellar feedback processes inherent to star formation. We use Nbody6++GPU to follow the evolution of the young stellar systems after gas removal. We use HDBSCAN to define stellar groups and analyze the stellar kinematics to identify the true bound star clusters. The conditions modeled by the simulations, i.e., global cloud surface densities below 0.15 g cm$^{-2}$,, star formation efficiencies below 15%, and gas expulsion timescales shorter than a free fall time, primarily produce expanding stellar associations and small clusters. The largest star clusters, which have $\sim$1000 bound members, form in the densest and lowest velocity dispersion clouds, representing $\sim$32 and 39% of the stars in the simulations, respectively. The cloud's early dynamical state plays a significant role in setting the classical star formation efficiency versus bound fraction relation. All stellar groups follow a narrow mass-velocity dispersion power law relation at 10 Myr with a power law index of 0.21. This correlation result in a distinct mass-size relationship for bound clusters. We also provide valuable constraints on the gas dispersal timescale during the star formation process and analyze the implications for the formation of bound systems.

We present detailed morphology measurements for 8.67 million galaxies in the DESI Legacy Imaging Surveys (DECaLS, MzLS, and BASS, plus DES). These are automated measurements made by deep learning models trained on Galaxy Zoo volunteer votes. Our models typically predict the fraction of volunteers selecting each answer to within 5-10\% for every answer to every GZ question. The models are trained on newly-collected votes for DESI-LS DR8 images as well as historical votes from GZ DECaLS. We also release the newly-collected votes. Extending our morphology measurements outside of the previously-released DECaLS/SDSS intersection increases our sky coverage by a factor of 4 (5,000 to 19,000 deg$^2$) and allows for full overlap with complementary surveys including ALFALFA and MaNGA.

After the initial suggestion by Zahnle and Walker (1987) that the torque accelerating the spin rate of the Earth and produced by the heating of the atmosphere by the Sun could counteract the braking lunir-solar gravitational torque in the Precambrian, several authors have recently revisited this hypothesis. In these studies, it is argued that the geological evidences of the past spin state of the Earth play in favor of this atmospheric tidal locking of the length of the day (LOD). In the present review of the recent literature, we show that the drawn conclusions depend crucially on the consideration of the stromatolite geological LOD estimates obtained by Pannella at 1.88 and 2.0 Ga, which are subject to large uncertainties. When only the most robust cyclostatigraphic estimates of the LOD are retained, the LOD locking hypothesis is not supported. Moreover, the consideration of the published General Circulation Model numerical simulations and of new analytical models for the thermal atmospheric tides suggest that the atmospheric tidal resonance, which is the crucial ingredient for the LOD locking in the Precambrian, was never of sufficiently large amplitude to allow for this tidal LOD lock.

The non-detection of a coma surrounding 1I/`Oumuamua, the first discovered interstellar object (ISO), has prompted a variety of hypotheses to explain its nongravitational acceleration. Given that forthcoming surveys are poised to identify analogues of this enigmatic object, it is prudent to devise alternative approaches to characterization. In this study, we posit X-ray spectroscopy as a surprisingly effective probe of volatile ISO compositions. Heavily ionized metals in the solar wind interact with outgassed neutrals and emit high-energy photons in a process known as charge exchange, and charge exchange induced X-rays from comets and planetary bodies have been observed extensively in our Solar System. We develop a model to predict the X-ray flux of an ISO based on its chemical inventory and ephemeris. We find that while standard cometary constituents, such as H$_2$O, CO$_2$, CO, and dust are best probed via optical or infrared observations, we predict strong X-ray emission generated by charge exchange with extended comae of H$_2$ and N$_2$ -- species which lack strong infrared fluorescence transitions. We find that XMM-Newton would have been sensitive to charge exchange emission from 1I/`Oumuamua during the object's close approach to Earth, and that constraints on composition may have been feasible. We argue for follow-up X-ray observations of newly discovered ISOs with close-in perihelia. Compositional constraints on the general ISO population could reconcile the apparently self-conflicting nature of 1I/`Oumuamua, and provide insight into the earliest stages of planet formation in extrasolar systems.

Centaurs are minor solar system bodies with orbits transitioning between those of Trans-Neptunian Scattered Disk objects and Jupiter Family comets. 39P/Oterma is a frequently active Centaur that has recently held both Centaur and JFC classifications and was observed with the JWST NIRSpec instrument on 2022 July 27 UTC while it was 5.82 au from the Sun. For the first time, CO$_2$ gas emission was detected in a Centaur, with a production rate of Q$_{CO_2}$ = (5.96 $\pm$ 0.80) $\times$ 10$^{23}$ molecules s$^{-1}$. This is the lowest detection of CO$_2$ of any Centaur or comet. CO and H$_2$O were not detected down to constraining upper limits. Derived mixing ratios of Q$_{CO}$/Q$_{CO_2}$ $\leq$2.03 and Q$_{CO_2}$/Q$_{H_2O}$ $\geq$0.60 are consistent with CO$_2$ and/or CO outgassing playing large roles in driving the activity, but not water, and show a significant difference between the coma abundances of 29P/Schwassmann-Wachmann 1, another Centaur at a similar heliocentric distance, which may be explained by thermal processing of 39P's surface during its previous Jupiter-family comet orbit. To help contextualize the JWST data we also acquired visible CCD imaging data on two dates in July (Gemini North) and September (Lowell Discovery Telescope) 2022. Image analysis and photometry based on these data are consistent with a point source detection and an estimated effective nucleus radius of 39P in the range of $R_{nuc}= $2.21 to 2.49~km.

Cosmology offers opportunities to test Dark Matter independently of its interactions with the Standard Model. We study the imprints of long-range forces acting solely in the dark sector on the distribution of galaxies, the so-called Large Scale Structure (LSS). We derive the strongest constraint on such forces from a combination of Planck and BOSS data. Along the way we consistently develop, for the first time, the Effective Field Theory of LSS in the presence of new dynamics in the dark sector. We forecast that future surveys will improve the current bound by an order of magnitude.

Expansions in the oscillation modes of tidally perturbed bodies provide a useful framework for representing tidally induced flows. However, recent work has demonstrated that such expansions produce inaccurate predictions for secular orbital evolution when mode damping rates are computed independently. We explore the coupling of collectively driven modes by frictional and viscous dissipation, in tidally perturbed bodies that are both non-rotating and rigidly rotating. This exploration leads us to propose an alternative approach to treating the damping of tidally driven oscillations that accounts for dissipative mode coupling, but which does not require any information beyond the eigenfunctions and eigenfrequencies of adiabatic modes.

Spontaneous symmetry breaking in grand unified theories is thought to have produced an exceedingly large number of magnetic monopoles in the early Universe. In the absence of suppression or annihilation, these very massive particles should be dominating the cosmic energy budget today, but none has ever been found. Inflation was invented in part to dilute their number, thereby rendering their density undetectable by current instruments. Should the inflationary paradigm not survive, however, the ensuing disagreement between theory and observation would constitute a cosmological `monopole problem' and create further tension for any extension to the standard model of particle physics. But as is also true for all horizon problems, a monopole overabundance emerges only in cosmologies with an initial period of deceleration. We show that the alternative Friedmann-Lemaitre-Robertson-Walker cosmology known as the R_h=ct universe completely eliminates all such anomalies rather trivially and naturally, without the need for an inflated expansion. We find that the monopole energy density today would be completely undetectable in R_h=ct. Evidence continues to grow that the zero active mass condition from general relativity ought to be an essential ingredient in LCDM.

We study the cosmological collider signatures in the Higgs-$R^2$ inflation model. We consider two distinct types of signals: one originating from the inflaton coupling to Standard Model fermions and gauge bosons, and another arising from the isocurvature mode interaction with the inflaton. In the former case, we determine that the signal magnitude is likely too small for detection by upcoming probes, primarily due to suppression by both the Planck scale and slow-roll parameters. However, we provide a detailed computation of the signal which could be potentially applicable to various Higgs inflation variants. For the isocurvature mode signals, we observe that the associated couplings remain unsuppressed when the isocurvature mode is relatively light or comparable to the inflationary scale. In this case, we study the Higgs-$R^2$ inflation parameter space that corresponds to the quasi-single-field inflation regime and find that the signal strength could be as large as $|f_{\rm NL}| > 1$, making Higgs-$R^2$ inflation a viable candidate for observation by future 21-cm surveys.

We show that Milky Way white dwarfs are excellent targets for dark matter (DM) detection. Using Fermi and H.E.S.S. Galactic center gamma-ray data, we investigate sensitivity to DM annihilating within white dwarfs into long-lived or boosted mediators and producing detectable gamma rays. Depending on the Galactic DM distribution, we set new constraints on the spin-independent scattering cross section down to $10^{-45}-10^{-41}$ cm$^2$ in the sub-GeV DM mass range, which is multiple orders of magnitude stronger than existing limits. For a generalized NFW DM profile, we find that our white dwarf constraints exceed spin-independent direct detection limits across most of the sub-GeV to multi-TeV DM mass range, achieving sensitivities as low as about $10^{-46}$ cm$^2$. In addition, we improve earlier versions of the DM capture calculation in white dwarfs, by including the low-temperature distribution of nuclei when the white dwarf approaches crystallization. This yields smaller capture rates than previously calculated by a factor of a few up to two orders of magnitude, depending on white dwarf size and the astrophysical system.

The most general tree-level boundary correlation functions of quantum fields in inflationary spacetime involve multiple exchanges of massive states in the bulk, which are technically difficult to compute due to the multi-layer nested time integrals in the Schwinger-Keldysh formalism. On the other hand, correlators with multiple massive exchanges are well motivated in cosmological collider physics, with the original quasi-single-field inflation model as a notable example. In this work, with the partial Mellin-Barnes representation, we derive a simple rule, called family-tree decomposition, for directly writing down analytical answers for arbitrary nested time integrals in terms of multi-variable hypergeometric series. We present the derivation of this rule together with many explicit examples. This result allows us to obtain analytical expressions for general tree-level inflation correlators with multiple massive exchanges. As an example, we present the full analytical results for a range of tree correlators with two massive exchanges.

Several modified gravity theories predict a possible time variation of the propagation speed of gravitational waves (GW) which could be tested with multimessenger astronomy. For this purpose we derive the relation between the redshift dependence of the propagation speed of GWs and the time delay between the detection of GWs and electromagnetic waves (EMWs) emitted by the same source. For theories with Einstein frame minimal matter-gravity coupling (EMC) the propagation speed of GWs can be jointly constrained by the time delay between GWs and EWs and the GW-EMW luminosity distance ratio, allowing to derive a consistency relation between these two observables. A violation of this consistency relation could support a non-minimal coupling of matter to the Einstein frame metric, as in Chameleon models, allowing to test such theories with multimessenger astronomy.

In this study, my main goal is to examine the nuclear matter properties across a wide range of conditions, such as temperature, density, asymmetry, pressure, and magnetic field. Understanding the effect of these factors on nuclear matter is essential, given their relevance in various phenomena such as heavy-ion collisions, neutron stars, and supernovae. However, due to the absence of a fundamental nuclear theory, we must rely on models to describe the nuclear matter. Predicted properties like neutron star mass-radius relationships, tidal deformability, structure, critical points in the nuclear matter phase diagram etc. depend on the chosen model. This dependence arises because key parameters characterizing any nuclear model or equation of state (EoS) are not precisely known. Therefore, it is crucial to investigate how nuclear matter properties behave under various conditions and in relation to different EoS parameters. To accomplish this, I have examined three distinct forms of nuclear matter: infinite nuclear matter, finite nuclei, and neutron stars, using the effective relativistic mean field model (E-RMF).

Cooling electronic devices to cryogenic temperatures (< 77 K) is crucial in various scientific and engineering domains. Efficient cooling involves the removal of heat generated from these devices through thermal contact with either a liquid cryogen or a dry cryostat cold stage. However, as these devices cool, thermal boundary resistance, also known as Kapitza resistance, hinders the heat flow across thermal interfaces, resulting in elevated device temperatures. In transistors, the presence of passivation layers like Silicon Nitride (SiN) introduces additional interfaces that further impede heat dissipation. This paper investigates the impact of passivation layer thickness on Kapitza resistance at the interface between a solid device and liquid nitrogen. The Kapitza resistance is measured using a capacitance thermometer that has been passivated with SiN layers ranging from 0 to 240 nm. We observe that Kapitza resistance increases with increasing passivation thickness.

We study a class of early universe cosmological models based on Einstein-Cartan gravity and including a higher derivative term corresponding to a power of the Holst scalar curvature. The resulting effective action is basically given by General Relativity and an additional neutral pseudoscalar field (the pseudoscalaron), unequivocally related to the corresponding components of the torsion, that necessarily acquire a dynamics. The induced pseudoscalaron potential provides a realistic inflationary phase together with a very rich postinflationary epoch, resulting by the coupling of the pseudoscalaron to ordinary matter.

In this work, we explicitly construct the vacuum solution of Einstein's equations with prescribed multipole moments. By observing the behavior of the multipole spacetime metric at small distances, we conjecture that for a sufficiently large multipole moment, there is a minimal size below which no object in nature can support such a moment. The examples we have investigated suggest that such minimal size scales as $(M_n)^{1/(n+1)}$ (instead of $(M_n/M)^{1/n}$), where $M$ is the mass and $M_n$ is the $n$th order multipole moment. With the metric of the "multipole spacetime", we analyze the shape of black hole shadow for various multipole moments and discuss the prospects of constraining the moments from shadow observations. In addition, we discuss the shift of gravitational wave phase with respect to those of the Kerr spacetime, for a test particle moving around an object with this set of multipole moments. These phase shifts are required for the program of mapping out the spacetime multipole moments based on gravitational wave observations of extreme mass-ratio inspirals.

The expansion history and thermal physical process that happened in the early Universe before big bang nucleosynthesis (BBN) remains relatively unconstrained by observations. Low reheating temperature universes with normalcy temperatures of $T_\mathrm{RH}\sim 2\,\mathrm{MeV}$ remain consistent with all observations, and accommodate several new physics scenarios that would normally be constrained by high-temperature reheating models, including massive sterile neutrinos. We explore such scenarios' production of keV scale sterile neutrinos and their resulting constraints from cosmological observations. The parameter space for massive sterile neutrinos is much less constrained than in high-$T_\mathrm{RH}$ thermal histories, though several cosmological constraints remain. Such parameter space is the target of several current and upcoming laboratory experiments such as TRISTAN (KATRIN), HUNTER, MAGNETO-$\nu$, and PTOLEMY. Cosmological constraints remain stringent for stable keV-scale sterile neutrinos. However, we show that sterile neutrinos with a dark decay to radiation through a $Z^\prime$ or a new scalar are largely unconstrained by cosmology. In addition, this mechanism of sterile neutrinos with large mixing may provide a solution to the Hubble tension. We find that keV-scale sterile neutrinos are therefore one of the best probes of the untested pre-BBN era in the early Universe and could be seen in upcoming laboratory experiments.