Papers for Thursday, Feb 09 2023

Short-period comet 108P/Ciffreo is known for its peculiar double morphology, in which the nucleus is accompanied by a co-moving, detached, diffuse 'blob'. We report new observations of 108P/Ciffreo taken with the Hubble Space Telescope and the Nordic Optical Telescope and use them to determine the cause of this unusual morphology. The separation and the longevity of the blob across several orbits together rule out the possibility of a single, slow-moving secondary object near the primary nucleus. We use a model of coma particle dynamics under the action of solar gravity and radiation pressure to show that the blob is an artifact of the turn-around of particles ejected sunward and repelled by sunlight. Numerical experiments limit the range of directions which can reproduce the morphology and explain why the co-moving blob appearance is rare.

Star formation and stellar feedback are interlinked processes that redistribute energy and matter throughout galaxies. When young, massive stars form in spatially clustered environments, they create pockets of expanding gas termed superbubbles. As these processes play a critical role in shaping galaxy discs and regulating the baryon cycle, measuring the properties of superbubbles provides important input for galaxy evolution models. With wide coverage and high angular resolution ($\sim$50-150 pc) of the PHANGS-ALMA $^{12}$CO (2-1) survey, we can now resolve and identify a statistically representative number of superbubbles with molecular gas in nearby galaxies. We identify superbubbles by requiring spatial correspondence between shells in CO with stellar populations identified in PHANGS-HST, and combine the properties of the stellar populations with CO to constrain feedback models and quantify their energetics. We visually identify 325 cavities across 18 PHANGS-ALMA galaxies, 88 of which have clear superbubble signatures (unbroken shells, central clusters, kinematic signatures of expansion). We measure their radii and expansion velocities using CO to dynamically derive their ages and the mechanical power driving the bubbles, which we use to compute the expected properties of the parent stellar populations driving the bubbles. We find consistency between the predicted and derived stellar ages and masses of the stellar populations if we use a supernova (SN) model that injects energy with a coupling efficiency of 5-12%, whereas wind models fail to explain stellar ages we measure. Not only does this confirm molecular gas accurately traces superbubble properties, but it also provides key observational constraints for superbubble models. We also find evidence that the bubbles are sweeping up gas as they expand and speculate that these sites have the potential to host new generations of stars.

Scientific Complementary Metal-Oxide-Semiconductor (CMOS) detectors have developed quickly in recent years, owing to their low cost, high availability and some advantages over CCDs, such as high frame rate or typically lower readout noise. With the development of the first back-illuminated models, these sensors started to be used in astronomy, so it is worth studying their characteristics, advantages and weaknesses. In this paper, we present the results of the laboratory characterization of the IMX455M and IMX411M sensors, integrated into the QHY600 and QHY411 cameras respectively. These are large (36$\times$24 and 54$\times$40 mm) native 16-bit sensors with 3.76 $\mu$m pixels sensitive in the optical range. Their quantum efficiency has been found to peak at 80% at 475 nm, 40% at 700 nm and 10% at 900 nm. Their linearity and photon transfer performance have been evaluated, as well as their dark behaviour. They showed a low dark current, but also the presence of warm pixels of about 0.024% in the QHY600 and 0.005% in the QHY411, which were proved to be stable and linear with exposure time. We have analysed in detail the effect of random telegraph noise, also called Salt & Pepper noise, one of the most important issues to be addressed with these two sensors, since it affects around 2% of the pixels in each exposure. Sky tests are also presented and the effect of this noise on astronomical images is discussed.

The Type I active galactic nucleus (AGN) ESO 511-G030, a formerly bright and soft-excess dominated source, has been observed in 2019 in the context of a multi-wavelength monitoring campaign. However, in these novel exposures, the source was found in a $\sim$10 times lower flux state, without any trace of the soft-excess. Interestingly, the X-ray weakening corresponds to a comparable fading of the UV suggesting a strong link between these components. The UV/X-ray spectral energy distribution (SED) of ESO 511-G030 shows remarkable variability. We tested both phenomenological and physically motivated models on the data finding that the overall emission spectrum of ESO 511-G030 in this extremely low flux state is the superposition of a power-law-like continuum ($\Gamma\sim$1.7) and two reflection components emerging from hot and cold matter. has Both the primary continuum and relativistic reflection are produced in the inner regions. The prominent variability of ESO 511-G030 and the lack of a soft-excess can be explained by the dramatic change in the observed accretion rate, which dropped from an L/L$_{\rm Edd}$ of 2\% in 2007 to 0.2\% in 2019. The X-ray photon index also became harder during the low flux 2019 observations, perhaps as a result of a photon starved X-ray corona.

In the past decades, several neutron stars (NSs), particularly pulsars, with mass $M>2M_\odot$ have been observed. On the other hand, the existence of massive white dwarfs (WDs), even violating Chandrasekhar mass-limit, was inferred from the peak luminosities of type Ia supernovae. Hence, there is a generic question of the origin of massive compact objects. Here we explore the existence of massive, magnetized, rotating NSs with soft and steep equation of states (EoSs) by solving axisymmetric stationary stellar equilibria in general relativity. For our purpose, we consider the Einstein equation solver for stellar structure XNS code. Such rotating NSs with magnetic field and rotation axes misaligned, hence with non-zero obliquity angle, can emit continuous gravitational waves (GW), which can be detected by upcoming detectors, e.g., Einstein Telescope, etc. We discuss the decays of magnetic field, angular velocity and obliquity angle with time, due to angular momentum extraction by GW and dipole radiation, which determine the timescales related to the GW emission. Further, in the Alfv\'en timescale, a differentially rotating, massive proto-NS rapidly settles into an uniformly rotating, less massive NS due to magnetic braking and viscosity. These explorations suggest that detecting massive NSs is challenging and sets a timescale for detection. We calculate the signal-to-noise ratio of GW emission, which confirms that any detector cannot detect them immediately, but detectable by Einstein Telescope, Cosmic Explorer over months of integration time, leading to direct detection of NSs.

In this paper we present a chemical and kinematic analysis of two ultra-faint dwarf galaxies (UFDs), Aquarius II (Aqu~II) and \text{Bo\"{o}tes II} (Boo~II), using Magellan/IMACS spectroscopy. We present the largest sample of member stars for Boo~II (12), and the largest sample of red-giant-branch members with metallicity measurements for Aqu~II (8). In both UFDs, over 80\% of targets selected based on $Gaia$ proper motions turned out to be spectroscopic members. In order to maximize the accuracy of stellar kinematic measurements, we remove the identified binary stars and RR Lyrae variables. For Aqu~II we measure a systemic velocity of $-65.3 \pm 1.8$ km s$^{-1}$ and a metallicity of [Fe/H] = $-2.57^{+0.17}_{-0.17}$. When compared with previous measurements, these values display a $\sim 6$ km s$^{-1}$ difference in radial velocity and a decrease of 0.27 dex in metallicity. Similarly for Boo~II, we measure a systemic velocity of $-130.4^{+1.4}_{-1.1}$ km s$^{-1}$, more than 10 km s$^{-1}$ different from the literature, a metallicity almost 1 dex smaller at [Fe/H] = $-2.71^{+0.11}_{-0.10}$, and a velocity dispersion 3 times smaller at $\sigma_{v_{\rm hel}} = 2.9^{+1.6}_{-1.2}$ km s$^{-1}$. Additionally, we derive systemic proper motion parameters and model the orbits of both UFDs. Finally, we highlight the extremely dark matter dominated nature of Aqu~II and compute the J-factor for both galaxies to aid searches of dark matter annihilation. Despite the small size and close proximity of Boo~II, it is an intermediate target for the indirect detection of dark matter annihilation due to its low velocity dispersion and corresponding low dark matter density.

Multifrequency studies of galaxy clusters are crucial for inferring their dynamical states and physics. Moreover, these studies allow us to investigate cluster-embedded sources, whose evolution is affected by the physical and dynamical condition of the cluster itself. So far, these kinds of studies have been preferentially conducted on clusters visible from the northern hemisphere due to the high-fidelity imaging capabilities of ground-based radio interferometers located there. In this paper, we conducted a multifrequency study of the poorly known galaxy cluster Abell 3718. We investigated the unknown origin of an extended radio source with a length of $\sim$612 kpc at 943 MHz detected in images from the Evolutionary Map of the Universe (EMU) and POlarisation Sky Survey of the Universe's Magnetism (POSSUM) surveys. We analyzed optical and X-ray data to infer the dynamical state of the cluster and, in particular, the merger activity. We conducted a radio spectral index study from 943 MHz up to 9 GHz. We also evaluated the polarization properties of the brightest cluster-embedded sources to understand if they are related to the radio emission observed on larger scales. [Abstract truncated due to arxiv limit! Please see the pdf version]

We present initial results of an ongoing survey with the Karl G. Jansky Very Large Array targeting the CO($J$ = 1-0) transition in a sample of 30 submillimeter-selected, dusty star-forming galaxies at $z =$ 2-5 with existing mid--$J$ CO detections from ALMA and NOEMA, of which 17 have been fully observed. We detect CO(1-0) emission in 11 targets, along with three tentative ($\sim$1.5-2$\sigma$) detections; three galaxies are undetected. Our results yield total molecular gas masses of 6-23$\times$10$^{10}$ ($\alpha_\mathrm{CO}$/1) M$_\odot$, with gas mass fractions, $f_\mathrm{gas}$=$M_\mathrm{mol}$/($M_*$+$M_\mathrm{mol}$), of 0.1-0.8 and a median depletion time of (140$\pm$70) Myr. We find median CO excitation ratios of $r_{31}$ = 0.75$\pm$0.39 and $r_{41}$ = 0.63$\pm$0.44, with a significant scatter. We find no significant correlation between the excitation ratio and a number of key parameters such as redshift, CO(1-0) line width or $\Sigma_\mathrm{SFR}$. We only find a tentative positive correlation between $r_{41}$ and the star-forming efficiency, but we are limited by our small sample size. Finally, we compare our results to predictions from the SHARK semi-analytical model, finding a good agreement between the molecular gas masses, depletion times and gas fractions of our sources and their SHARK counterparts. Our results highlight the heterogeneous nature of the most massive star-forming galaxies at high-redshift, and the importance of CO(1--0) observations to robustly constrain their total molecular gas content and ISM properties.

Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two environmental processes. We compare simulations with and without extinction of photoevaporation-driving radiation. We find that the majority of discs in our simulations are considerably shielded from photoevaporation-driving radiation for at least 0.5 Myr after the formation of the first massive stars. Radiation shielding increases disc lifetimes by an order of magnitude and can let a disc retain more solid material for planet formation. The reduction in external photoevaporation leaves discs larger and more easily dynamically truncated, although external photoevaporation remains the dominant mass loss process. Finally, we find that the correlation between disc mass and projected distance to the most massive nearby star (often interpreted as a sign of external photoevaporation) can be erased by the presence of less massive stars that dominate their local radiation field. Overall, we find that the presence and dynamics of gas in embedded clusters with massive stars is important for the evolution of protoplanetary discs.

The radial transport, or drift, of dust has taken a critical role in giant planet formation theory. However, it has been challenging to identify dust drift pile ups in the hard-to-observe inner disk. We find that the IM Lup disk shows evidence that it has been shaped by an episode of dust drift. Using radiative transfer and dust dynamical modeling we study the radial and vertical dust distribution. We find that high dust drift rates exceeding 110 M_earth/Myr are necessary to explain both the dust and CO observations. Furthermore, the bulk of the large dust present in the inner 20 au needs to be vertically extended, implying high turbulence alpha_z > 10^{-3} and small grains (0.2-1 mm). We suggest that this increased level of particle stirring is consistent with the inner dust-rich disk undergoing turbulence triggered by the vertical shear instability. The conditions in the IM Lup disk imply that giant planet formation through pebble accretion is only effective outside 20 au. If such an early, high turbulence inner region is a natural consequence of high dust drift rates, then this has major implications for understanding the formation regions of giant planets including Jupiter and Saturn.

Modeling the dynamics of the formation and evolution of protostellar disks as well as the history of stellar mass accretion typically involve the numerical solution of complex systems of coupled differential equations. The resulting mass accretion history of protostars is known to be highly episodic due to recurrent instabilities and also exhibits short timescale flickering. By leveraging the strong predictive abilities of neural networks, we extract some of the critical temporal dynamics experienced during the mass accretion including periods of instability. Particularly, we utilize a novel form of the Echo-State Neural Network (ESN), which has been shown to efficiently deal with data having inherent nonlinearity. We introduce the use of Optimized-ESN (Opt-ESN) to make model-independent time series forecasting of mass accretion rate in the evolution of protostellar disks. We apply the network to multiple hydrodynamic simulations with different initial conditions and exhibiting a variety of temporal dynamics to demonstrate the predictability of the Opt-ESN model. The model is trained on simulation data of $\sim 1-2$ Myr, and achieves predictions with a low normalized mean square error ($\sim 10^{-5}$ to $10^{-3}$) for forecasts ranging between 100 and 3800 yr. This result shows the promise of the application of machine learning based models to time-domain astronomy.

Context: Most stars in the galactic stellar population are low-mass stars. Very low-mass (VLM) stars are a subset of the low-mass stars typically defined in terms of the stellar masses ranging from 0.6 M_sun to the hydrogen-burning limit of about 0.075 M_sun. Aim: The observational studies of VLM binaries can provide effective diagnostics for testing the VLM formation scenarios. The small size of VLMs makes them suitable candidates to detect planets around them in the habitable zone. Methods: In this work, using the high-resolution near-infrared adaptive optics imaging from the NaCo instrument installed on the Very Large Telescope, we report the discovery of a new binary companion to the M-dwarf LP 1033-31 and also confirm the binarity of LP 877-72. We have characterized both stellar systems and estimated the properties of their individual components. Results and Conclusions: We have found that LP 1033-31 AB with the spectral type of M4.5+M4.5 has a projected separation of 6.7+/-1.3 AU. On the other hand, with the spectral type of M1+M4, the projected separation of LP 877-72 AB is estimated to be 45.8+/-0.3 AU. We further investigated the masses, surface gravity, radii, and effective temperature of the detected components. The orbital period of LP 1033-31 and LP 877-72 systems are estimated to be ~28 and ~349 yr, respectively. Our analysis suggests that there is a possibility of finding up to `two' exoplanets around LP 877-72 B. In contrast, the maximum probabilities of hosting exoplanets around LP 877-72 A, LP 1033-31 A, and LP 1033-31 B are estimated to be only ~50%.

The distribution of galaxies and clusters of galaxies on the mega-parsec scale of the Universe follows an intricate pattern now famously known as the Large-Scale Structure or the Cosmic Web. To study the environments of this network, several techniques have been developed that are able to describe its properties and the properties of groups of galaxies as a function of their environment. In this work we analyze the previously introduced framework: 1-Dimensional Recovery, Extraction, and Analysis of Manifolds (1-DREAM) on N-body cosmological simulation data of the Cosmic Web. The 1-DREAM toolbox consists of five Machine Learning methods, whose aim is the extraction and modelling of 1-dimensional structures in astronomical big data settings. We show that 1-DREAM can be used to extract structures of different density ranges within the Cosmic Web and to create probabilistic models of them. For demonstration, we construct a probabilistic model of an extracted filament and move through the structure to measure properties such as local density and velocity. We also compare our toolbox with a collection of methodologies which trace the Cosmic Web. We show that 1-DREAM is able to split the network into its various environments with results comparable to the state-of-the-art methodologies. A detailed comparison is then made with the public code DisPerSE, in which we find that 1-DREAM is robust against changes in sample size making it suitable for analyzing sparse observational data, and finding faint and diffuse manifolds in low density regions.

We present an empirical measurement of the dark count rate seen in a large-format MKID array identical to those currently in use at observatories such as Subaru on Maunakea. This work provides compelling evidence for their utility in future experiments that require low-count rate, quiet environments such as dark matter direct detection. Across the bandpass from 0.946-1.534 eV (1310-808 nm) an average count rate of $(1.847\pm0.003)\times10^{-3}$ photons/pixel/s is measured. Breaking this bandpass into 5 equal-energy bins based on the resolving power of the detectors we find the average dark count rate seen in an MKID is $(6.26\pm0.04)\times10^{-4}$ photons/pixel/s from 0.946-1.063 eV and $(2.73\pm0.02)\times10^{-4}$ photons/pixel/s at 1.416-1.534eV. Using lower-noise readout electronics to read out a single MKID pixel we demonstrate that the events measured while the detector is not illuminated largely appear to be a combination of real photons, possible fluorescence caused by cosmic rays, and phonon events in the array substrate. We also find that using lower-noise readout electronics on a single MKID pixel we measure a dark count rate of $(9.3\pm0.9)\times10^{-4}$ photons/pixel/s over the same bandpass (0.946-1.534 eV) With the single-pixel readout we also characterize the events when the detectors are not illuminated and show that these responses in the MKID are distinct from photons from known light sources such as a laser, likely coming from cosmic ray excitations.

Ground-based 21-cm experiments targeting the global signal from the periods of Cosmic Dawn (CD) and Epoch of Reionization (EoR) are susceptible to adverse effects presented by i) the ionosphere ii) antenna chromaticity induced by objects in its vicinity iii) terrestrial radio frequency interference (RFI). Terrestrial RFI is particularly challenging as the FM radio band spanning over 88-108 MHz lies entirely within the frequency range of the CD/EoR experiments ($\sim 40-200$ MHz). Multiple space-based experiments have been proposed to operate in the radio-quiet zone on the lunar farside. An intermediate option in cost and complexity is an experiment operating in space in an orbit around Earth, which readily alleviates the first two challenges. However, the effect of RFI in Earth's orbit on the detection of global signal needs to be quantitatively evaluated. We present STARFIRE -- Simulation of TerrestriAl Radio Frequency Interference in oRbits around Earth -- an algorithm that provides an expectation of FM seeded RFI at different altitudes over Earth. Using a limited set of publicly available FM transmitter databases, which can be extended by the user community, we demonstrate the use of the STARFIRE framework to generate a three-dimensional spatio-spectral cube of RFI as would be measured in Earth orbit. Applications of STARFIRE include identifying minimum RFI orbits around Earth, producing RFI spectra over a particular location, and generating RFI heatmaps at specific frequencies for a range of altitudes. STARFIRE can be easily adapted for different frequencies, altitudes, antenna properties, RFI databases, and combined with astrophysical sky-models. This can be used to estimate the effect of RFI on the detection of global 21-cm signal from Earth-orbit, and hence for sensitivity estimates and experiment design of an Earth orbiting CD/EoR detection experiment.

The composition, including the ionic charge states and elemental abundances of heavy elements, within interplanetary coronal mass ejections (ICMEs) has tight correlations with their source regions and eruption processes. This can help analyze the eruption mechanisms and plasma origins of CMEs, and deepen our understanding of energetic solar activities. The active regions and quiet-Sun regions have different physical properties, thus from a statistical point of view, ICMEs originating from the two types of regions should exhibit different compositional characteristics. To demonstrate the differences comprehensively, we conduct survey studies on the ionic charge states of five elements (Mg, Fe, Si, C, and O) and the relative abundances of six elements (Mg/O, Fe/O, Si/O, C/O, Ne/O, and He/O) within ICMEs from 1998 February to 2011 August through the data of advanced composition explorer. The results show that ICMEs from active regions have higher ionic charge states and relative abundances than those from quiet-Sun regions. For the active-region ICMEs, we further analyze the relations between their composition and flare class, and find a positive relationship between them, i.e., the higher classes of the associated flares, the larger means of ionic charge states and relative abundances (except the C/O) within ICMEs. As more (less) fractions of ICMEs originate from active regions around solar maximum (minimum), and active-region ICMEs usually are associated with higher-class flares, our studies might answer why ICME composition measured near 1 au exhibits the solar cycle dependence.

In the present work we derive an analytical expression for the mass density of an object with spherical symmetry, whose corresponding potential allows obtaining a circular velocity around it that is in agreement with the observed rotation curve of galaxies. The rotation curve of our galaxy is analyzed, determining the properties of the central object, its radius and mass whose value obtained is very close to that reported in the literature.

We present photometric and spectroscopic observations of the extraordinary gamma-ray burst (GRB) 221009A in search of an associated supernova. Some past GRBs have shown bumps in the optical light curve that coincide with the emergence of supernova spectral features, but we do not detect any significant light curve features in GRB~221009A, nor do we detect any clear sign of supernova spectral features. Using two well-studied GRB-associated supernovae (SN~2013dx, $M_{r,max} = -19.54$; SN~2016jca, $M_{r,max} = -19.04$) at a similar redshift as GRB~221009A ($z=0.151$), we modeled how the emergence of a supernova would affect the light curve. If we assume the GRB afterglow to decay at the same rate as the X-ray data, the combination of afterglow and a supernova component is fainter than the observed GRB brightness. For the case where we assume the best-fit power law to the optical data as the GRB afterglow component, a supernova contribution should have created a clear bump in the light curve, assuming only extinction from the Milky Way. If we assume a higher extinction of $E(B-V)$=$1.74$ mag (as has been suggested elsewhere), the supernova contribution would have been hard to detect, with a limit on the associated supernova of $M_{r,max} \approx-$19.54. We do not observe any clear supernova features in our spectra, which were taken around the time of expected maximum light. The lack of a bright supernova associated with GRB~221009A may indicate that the energy from the explosion is mostly concentrated in the jet, leaving a lower energy budget available for the supernova.

Purpose: The High-Energy X-ray telescope (HE), one of the three main payloads of the \textit{Insight}-HXMT mission, is composed of eighteen NaI(Tl)/CsI(Na) phoswich detectors, where NaI(Tl) serves as the primary detector covering 20--250\,keV, and CsI(Na) is used as an active shield detector to suppress the background of NaI(Tl) and also serves as an all-sky gamma-ray burst monitor covering 0.2--3\,MeV. In this paper, we review the in-orbit performance of HE in the first 5 years since \textit{Insight}-HXMT was launched on June 15, 2017. Methods: The major performances we concern include the gain and energy resolution of NaI(Tl) and CsI(Na) detectors, the performance of pulse-shape-discriminator (PSD) and system dead-time. In this work, we investigate these performances mainly using the data of blank-sky observations and the data when the telescope in earth occultation. Results: The overall performance of HE/NaI(Tl) is very stable in the first 5 years, whereas the gain of HE/CsI(Na) shows a continuously increasing trend and should be calibrated regularly. Conclusion: In general, HE is still in good health and well-calibrated status after five-year's operation. The in-orbit performance of HE has no significant deviation from expectation. HE is expected to be in operation healthily for another several years of extended mission life.

We present a new version of our analytical model of the spatial interstellar extinction variations within the nearest kiloparsec. This model treats the 3D dust distribution as a superposition of three overlapping layers: (1) the layer along the Galactic midplane, (2) the layer in the Gould Belt, and (3) the layer passing through the Cepheus and Chamaeleon dust cloud complexes. In each layer the dust density decreases exponentially with increasing distance from the midplane of the layer. In addition, there are sinusoidal longitudinal extinction variations along the midplane of each layer. We have found the most probable values of 29 parameters of our model using four data sets: the 3D stellar reddening maps by Gontcharov (2017), Lallement et al. (2019), and Green et al. (2019) and the extinctions inferred by Anders et al. (2022) for 993291 giants from the Gaia EDR3. All of the data give similar estimates of the model parameters. The extinction for a star or a point in space is predicted by our model with an accuracy from $\sigma(A_\mathrm{V})=0.07$ to 0.37 for high and low Galactic latitudes, respectively. The natural fluctuations of the dust medium dominate in these values. When ignoring the fluctuations of the medium, the average extinction for an extended object (a galaxy, a star cluster, a dust cloud) or a small region of space is predicted by our model with an accuracy from $\sigma(A_\mathrm{V})=0.04$ to 0.15 for high and low Galactic latitudes, respectively. Green et al. (2019) and Anders et al. (2022) give in unison an estimate of $A_\mathrm{V}=0.12^m$ for the extinction at high latitudes across the whole Galactic dust half-layer above or below the Sun with the natural fluctuations of the medium $\sigma(A_\mathrm{V})=0.06^m$. Our model is a step to explain how a substantial amount of dust ended up far from the Galactic midplane.

The detection of planets around other stars by the measurement of the stellar Radial Velocity (RV) variations benefits from improvements of dedicated spectrographs, allowing to achieve a precision of 1 ms$^{-1}$ or better. Spectral intervals within which stellar lines are contaminated by telluric lines are classically excluded from the RV processing. We aim at estimating the potential improvement of telluric absorption removal and subsequent extension of the useful spectral domain on the precision of radial velocity measurements. We developed a correction method based on the on-line web service TAPAS, allowing to determine a synthetic atmospheric transmission spectrum for the time and location of observations. This method was applied to the telluric H$_{2}$O and O$_2$ absorption removal from a series of 200 ESPRESSO consecutive exposures of the K2.5V star HD40307, available in ESO archives. We calculated the radial velocity using the standard Cross-Correlation Function (CCF) method and Gaussian fit of the CCF, with uncorrected spectra and the ESPRESSO standard stellar binary mask on one hand, and telluric-corrected spectra and an augmented binary mask with 696 additional lines on the other hand. We find that the precision of radial velocity measurements is improved in the second case, with a reduction of the average formal error from 1.04 ms$^{-1}$ down to 0.78 ms$^{-1}$ in the case of these ESPRESSO data and this stellar type for the red arm. Using an estimator of the minimal error based on photon noise limit applied to the full CCF, the error is reduced from 0.89 ms$^{-1}$ down to 0.78 ms$^{-1}$. This corresponds to a significant decrease of about 35\% of observing time to reach the same precision in the red part.

The continuous wavelet transform (CWT) is very useful for processing signals with intricate and irregular structures in astrophysics and cosmology. It is crucial to propose precise and fast algorithms for the CWT. In this work, we review and compare four different fast CWT algorithms for the 1D signals, including the FFTCWT, the V97CWT, the M02CWT, and the A19CWT. The FFTCWT algorithm implements the CWT using the Fast Fourier Transform (FFT) with a computational complexity of $\mathcal{O}(N\log_2N)$ per scale. The rest algorithms achieve the complexity of $\mathcal{O}(N)$ per scale by simplifying the CWT into some smaller convolutions. We illustrate explicitly how to set the parameters as well as the boundary conditions for them. To examine the actual performance of these algorithms, we use them to perform the CWT of signals with different wavelets. From the aspect of accuracy, we find that the FFTCWT is the most accurate algorithm, though its accuracy degrades a lot when processing the non-periodic signal with zero boundaries. The accuracy of $\mathcal{O}(N)$ algorithms is robust to signals with different boundaries, and the M02CWT is more accurate than the V97CWT and A19CWT. From the aspect of speed, the $\mathcal{O}(N)$ algorithms do not show an overall speed superiority over the FFTCWT at sampling numbers of $N\lesssim10^6$, which is due to their large leading constants. Only the speed of the V97CWT with real wavelets is comparable to that of the FFTCWT. However, both the FFTCWT and V97CWT are substantially less efficient in processing the non-periodic signal because of zero padding. Finally, we conduct wavelet analysis of the 1D density fields, which demonstrate the convenience and power of techniques based on the CWT. We publicly release our CWT codes as resources for the community.

A generalization of inflationary $\alpha$-attractor models was recently proposed by Kallosh and Linde, in which the potential involves logarithmic functions of the inflaton so that the derivative of the potential but not potential itself has a singularity. We find that the models can lead to viable inflationary observables even without the pole in the kinetic term. Also, the generalization with a pole order other than two does not significantly change the functional form of the potential. This allows a systematic analysis of the predictions of this class of models. The models typically predict the scalar spectral index $n_s$ around 0.97 and values of the tensor-to-scalar ratio $r$ observable by LiteBIRD. Taking advantage of the relatively large $n_s$, we discuss the modification of the potential to produce primordial black holes as the whole dark matter and gravitational waves induced by curvature perturbations detectable by LISA and BBO/DECIGO, while keeping $n_s$ in agreement with the Planck/BICEP/Keck data.

The X-shooter Spectral Library (XSL) is a large empirical stellar library used as a benchmark for the development of stellar population models. The inclusion of $\alpha$-elements abundances is crucial to disentangling the chemical evolution of any stellar system. The aim of this paper is to provide a catalogue of high-precision and accurate magnesium and calcium abundances from a wide variety of stars well distributed in the Hertzsprung-Russell (HR) diagram. We originally performed an analysis of the derived Mg and Ca abundances for medium-resolution spectra of 611 stars from the XSL Data Release 2. For this purpose, we used the GAUGUIN automated abundance estimation code to fit the ultraviolet-blue (UVB) and visible (VIS) spectra. We tested the consistency of the atmospheric parameters and chemical abundances with the Gaia DR3 and the AMBRE Project datasets. We have finally obtained precise [Mg/Fe] and [Ca/Fe] abundances for 192 and 217 stars respectively, from which 174 stars have measurements in both elements. The stars cover a broad range of effective temperature 4000 < T$_{\rm eff}$ < 6500K, surface gravity 0.3 < log(g) < 4.8~cm s$^{\rm -2}$, and metallicity -2.5 < [Fe/H] < +0.4 dex. We find an excellent agreement with the abundance estimates from the AMBRE:HARPS and the Gaia/RVS analysis. Moreover, the resulting abundances reproduce a plateau in the metal-poor regime followed by a decreasing trend even at supersolar metallicities, as predicted by Galactic chemical evolution models. This catalogue is suitable for improving the modelling of evolutionary stellar population models with empirical $\alpha$-enhancements, which could significantly contribute to the analysis of external galaxies abundances in the near future.

We investigate the sensitivity of a universe's nuclear entropy after Big Bang nucleosynthesis (BBN) to variations in both the baryon-to-photon ratio and the temporal evolution of cosmological expansion. Specifically, we construct counterfactual cosmologies to quantify the degree by which these two parameters must vary from those in our Universe before we observe a substantial change in the degree of fusion, and thus nuclear entropy, during BBN. We find that, while the post-BBN nuclear entropy is indeed linked to baryogenesis and the Universe's expansion history, the requirement of leftover light elements does not place strong constraints on the properties of these two cosmological processes.

MAXI J1535-571 is a black-hole X-ray binary that in 2017 exhibited a very bright outburst which reached a peak flux of up to 5 Crab in the 2-20 keV band. Given the high flux, several X-ray space observatories obtained unprecedented high signal-to-noise data of key parts of the outburst. In our previous paper we studied the corona of MAXI J1535-571 in the hard-intermediate state (HIMS) with Insight-HXMT. In this paper we focus on the study of the corona in the soft-intermediate state (SIMS) through the spectral-timing analysis of 26 NICER detections of the type-B quasi-periodic oscillations (QPOs). From simultaneous fits of the energy, rms and lag spectra of these QPOs with our time-dependent Comptonization model, we find that in the SIMS the corona size is ~ 6500 km and vertically extended. We detect a narrow iron line in the energy spectra, which we interpret to be due to the illumination of the outer part of the accretion disk by this large corona. We follow the evolution of the corona and the radio jet during the HIMS-SIMS transition, and find that the jet flux peaks after the time when the corona extends to its maximum vertical size. The jet flux starts to decay after the corona contracts vertically towards the black hole. This behavior points to a connection between the X-ray corona and the radio jet similar to that seen in other sources.

Serving as the progenitors of electromagnetic and gravitational wave transients, massive stars have received renewed interest in recent years. However, many aspects of their birth and evolution remain opaque, particularly in the context of binary interactions. The centre of our galaxy hosts a rich cohort of very massive stars, which appear to play a prominent role in the ecology of the region. In this paper we investigate the binary properties of the Arches cluster, which is thought to host a large number of very massive stars. A combination of multi-epoch near-IR spectroscopy and photometry was utilised to identify binaries. 13 from 36 cluster members meet our criteria to be classed as RV variable. Combining the spectroscopic data with archival radio and X-ray observations - to detect colliding wind systems - provides a lower limit to the binary fraction of ~43%; increasing to >50% for the O-type hypergiants and WNLha. Dynamical and evolutionary masses reveal the primaries to be uniformly massive (>50M$_{\odot}$). Where available, orbital analysis reveals a number of short period, highly eccentric binaries, which appear to be pre-interaction systems. Such systems are X-ray luminous, with 80% above an empirical bound of $(L_{\rm x}/L_{\rm bol})\sim10^{-7}$ and their orbital configurations suggest formation and evolution via a single star channel; however, we cannot exclude a binary formation channel for a subset. Qualitative comparison to surveys of lower mass OB-type stars confirms that the trend to an extreme binary fraction (>60%) extends to the most massive stars currently forming in the local Universe.

This paper presents a new method to estimate systematic errors in the maximum-likelihood regression of count data. The method is applicable in particular to X-ray spectra in situations where the Poisson log-likelihood, or the Cash goodness-of-fit statistic, indicate a poor fit that is attributable to overdispersion of the data. Overdispersion in Poisson data is treated as an intrinsic model variance that can be estimated from the best-fit model, using the maximum-likelihood Cmin statistic. The paper also studies the effects of such systematic errors on the Delta C likelihood-ratio statistic, which can be used to test for the presence of a nested model component in the regression of Poisson count data. The paper introduces an overdispersed chi-square distribution that results from the convolution of a chi-square distribution that models the usual Delta C statistic, and a zero-mean Gaussian that models the overdispersion in the data. This is proposed as the distribution of choice for the Delta C statistic in the presence of systematic errors. The methods presented in this paper are applied to XMM-Newton data of the quasar 1ES 1553+113 that were used to detect absorption lines from an intervening warm-hot intergalactic medium (WHIM). This case study illustrates how systematic errors can be estimated from the data, and their effect on the detection of a nested component, such as an absorption line, with the Delta C statistic.

The mass and distance of a binary black hole (BBH) are fundamental parameters to measure in gravitational-wave (GW) astronomy. It is well-known that the measurement is affected by cosmological redshift, and recent works also showed that Doppler and gravitational redshifts could further affect the result if the BBH coalesces close to a supermassive black hole (SMBH). Here we consider the additional lensing effect induced by the nearby SMBH on the measurement. We compute the null geodesics originating within $10$ gravitational radii of a Kerr SMBH to determine the redshift and magnification of the GWs emitted by the BBH. We find a positive correlation between redshift and demagnification, which results in a positive correlation between the mass and distance of the BBH in the detector frame. More importantly, we find a higher probability for the signal to appear redshifted and demagnified to a distant observer, rather than blueshifted and magnified. Based on these results, we show that a binary at a cosmological redshift of $z_{\rm cos}=(10^{-2}-10^{-1})$ and composed of BHs of $(10-20)M_\odot$ could masquerade as a BBH at a redshift of $z_{\rm cos}\sim1$ and containing BHs as large as $(44-110)M_\odot$. In the case of extreme demagnification, we also find that the same BBH could appear to be at $z_{\rm cos}>10$ and contain subsolar-mass BHs. Such an effect, if not accounted for, could bias our understanding of the origin of the BHs detected via GWs.

We present a comprehensive simulation-based study of the BayesEoR code for 21 cm power spectrum recovery when analytically marginalizing over foreground parameters. To account for covariance between the 21 cm signal and contaminating foreground emission, BayesEoR jointly constructs models for both signals within a Bayesian framework. Due to computational constraints, the forward model is constructed using a restricted field-of-view (FoV) in the image domain. When the only EoR contaminants are noise and foregrounds, we demonstrate that BayesEoR can accurately recover the 21 cm power spectrum when the component of sky emission outside this forward-modelled region is downweighted by the beam at the level of the dynamic range between the foreground and 21 cm signals. However, when all-sky foreground emission is included along with a realistic instrument primary beam with sidelobes above this threshold extending to the horizon, the recovered power spectrum is contaminated by unmodelled sky emission outside the restricted FoV model. Expanding the combined cosmological and foreground model to cover the whole sky is computationally prohibitive. To address this, we present a modified version of BayesEoR that allows for an all-sky foreground model, while the modelled 21 cm signal remains only within the primary FoV of the telescope. With this modification, it will be feasible to run an all-sky BayesEoR analysis on a sizeable compute cluster. We also discuss several future directions for further reducing the need to model all-sky foregrounds, including wide-field foreground subtraction, an image-domain likelihood utilizing a tapering function, and instrument primary beam design.

Based on published data, we have assembled a sample of 126 radio stars with the trigonometric parallaxes and proper motions measured by VLBI and available in the Gaia DR3 catalogue (in fact, Gaia EDR3). Our analysis of the Gaia--VLBI proper motion differences for 84 radio stars based on the model of solid-body mutual rotation has revealed no rotation components differing significantly from zero, $(\omega_x,\omega_y,\omega_z)=(0.06,0.08,-0.10)\pm(0.06,0.07,0.08)$ mas yr$^{-1}.$ Based on the trigonometric parallax differences for 90 stars, we have obtained a new estimate of the systematic offset between the optical and radio frames, $\Delta\pi=-0.022\pm0.017$~mas, and showed that the parallax scale factor is close to unity, $b=1.001\pm0.002$.

Ultra-massive H-rich (DA spectral type) white dwarf stars ($M_{\star} > 1.05 M_{\odot}$) are expected to be substantially crystallized by the time they reach the ZZ Ceti instability strip ($T_{\rm eff} \sim 12\,000$ K). Crystallization leads to a separation of $^{16}$O and $^{20}$Ne (or $^{12}$C and $^{16}$O) in the core of ultra-massive WDs, which strongly impacts their pulsational properties. An additional factor to take into account when modeling the evolution and pulsations of WDs in this range of masses are the relativistic effects, which induce changes in the cooling times and the stellar masses derived from the effective temperature and surface gravity. Given the arrival of large amounts of photometric data from space missions such as {\it Kepler}/{\it K2} and {\it TESS}, it is important to assess the impact of General Relativity in the context of pulsations of ultra-massive ZZ Ceti stars. In this work, we present results of Newtonian gravity($g$)-mode pulsation calculations in evolutionary ultra-massive WD models computed in the frame of the General Relativity theory.

AGILE is a space mission launched in 2007 devoted to high-energy astrophysics. The AGILE Team is involved in the multi-messenger campaigns to send and receive science alerts about transient events in the shortest time possible. For this reason, the AGILE Team developed several real-time analysis pipelines to analyse data and follow-up external science alerts. However, the results obtained by these pipelines are preliminary and must be validated with manual analyses that are the bottleneck of the workflow. To speed up the scientific analysis performed by scientists, the AGILE Team developed the AGILEScience mobile application (for iOS and Android devices) that offers to the AGILE Team a password-protected section used to visualise the results of automated pipelines. We present in this contribution an improved functionality of the AGILEScience application that aims to enable the AGILE Team to execute a full scientific analysis using their mobile devices. When the analysis is completed, the system sends an email to notify the user that can visualise the results (e.g. plots, tables, and HTML pages) through the application. The possibility to perform scientific analysis from a mobile device enables the AGILE researchers to perform fast scientific analyses remotely to validate the preliminary results obtained with the automated pipelines. This workflow reduces the overall reaction time of the AGILE Team for the follow-up of transient phenomena.

We investigate prospects for the detection of high-energy neutrinos produced in the prolonged jets of short gamma-ray bursts (sGRBs). The X-ray lightcurves of sGRBs show extended emission components lasting for 100-1000 seconds, which are considered to be produced by prolonged engine activity. Jets by prolonged engine activity should interact with photons in the cocoon formed by the jet propagation inside the ejecta of neutron star mergers. We calculate neutrino emission from jets by prolonged engine activity, taking account of the interaction between photons provided from the cocoon and cosmic rays accelerated in the jets. We find that IceCube-Gen2, a future neutrino telescope, with the second-generation gravitational wave detectors will probably be able to observe neutrino signals associated with gravitational waves with around 10 years of operation, regardless of the assumed value of the Lorentz factor of the jets. Neutrino observations may enable us to constrain the dissipation region of the jets. We apply this model to GRB 211211A, a peculiar long GRB whose origin may be a binary neutron-star merger. Our model predicts that IceCube is unlikely to detect any associated neutrino, but a few similar events will be able to put a meaningful constraint on the physical quantities of the prolonged engine activities.

Earth-sized exoplanets that transit nearby, late spectral type red dwarfs will be prime targets for atmospheric characterization in the coming decade. Such systems, however, are difficult to find via wide-field transit surveys like Kepler or TESS. Consequently, the presence of such transiting planets is unexplored and the occurrence rates of short-period Earth-sized planets around late M dwarfs remain poorly constrained. Here, we present the deepest photometric monitoring campaign of 22 nearby late M dwarf stars, using data from over 500 nights on seven 1-2 meter class telescopes. Our survey includes all known single quiescent northern late M dwarfs within 15 pc. We use transit-injection-and-recovery tests to quantify the completeness of our survey, successfully identify most ($>80\%$) transiting short-period (0.5-1 d) super-Earths ($R > 1.9 R_\oplus$), and are sensitive ($\sim50\%$) to transiting Earth-sized planets ($1.0-1.2 R_\oplus$). Our high sensitivity to transits with a near-zero false positive rate demonstrates an efficient survey strategy. Our survey does not yield a transiting planet detection, yet it provides the most sensitive upper limits on transiting planets orbiting our target stars. Finally, we explore multiple hypotheses about the occurrence rates of short-period planets (from Earth-sized planets to giant planets) around late M dwarfs. We show, for example, that giant planets at short periods ($<1$ day) are uncommon around our target stars. Our dataset provides some insight into occurrence rates of short-period planets around TRAPPIST-1-like stars, and our results can help test planetary formation and system evolution models, as well as guide future observations of nearby late M dwarfs.

In view of the next LIGO-Virgo-KAGRA Observing period O4 (to start in Spring 2023), we address the question of the ability of the interferometers network to discriminate among different neutron stars equation of states better than what was possible with the observation of the binary neutron stars merger GW170817. We show that the observation of an event similar to GW170817 during O4 would allow to resolve the dimensionless effective tidal deformability $\tilde{\Lambda}$ within an uncertainty 7 times better than the one obtained in O2. Thanks to the expected increase in sensitivities, we show that any GW170817-like single-event within a distance of 100 Mpc would imply significantly improved constraints of the neutron stars equations of state. We also illustrate the important impact of the noise in the analysis of the signal, showing how it can impact the effective tidal deformability probability density function for large signal-to-noise ratio.

The galaxy distribution in dark matter-dominated halos is expected to approximately trace the details of the underlying dark matter substructure. In this paper we introduce halo `core-tracking' as a way to efficiently follow the small-scale substructure in cosmological simulations and apply the technique to model the galaxy distribution in observed clusters. The method relies on explicitly tracking the set of particles identified as belonging to a halo's central density core, once a halo has attained a certain threshold mass. The halo cores are then followed throughout the entire evolution of the simulation. The aim of core-tracking is to simplify substructure analysis tasks by avoiding the use of subhalos and, at the same time, to more easily account for the so-called ``orphan'' galaxies, which have lost substantial dark mass due to tidal stripping. We show that simple models based on halo cores can reproduce the number and spatial distribution of galaxies found in optically-selected clusters in the Sloan Digital Sky Survey. We also discuss future applications of the core-tracking methodology in studying the galaxy-halo connection.

The Atacama Large Millimeter/submillimeter Array (ALMA) serendipitously detected H$_2$O $J_{Ka, Kc} = 10_{2,9} - 9_{3,6}$ emission at 321 GHz in NGC 1052. This is the first submillimeter maser detection in a radio galaxy and the most luminous 321-GHz H$_2$O maser known to date with the isotropic luminosity of 1090 $L_{\odot}$. The line profile consists of a broad velocity component with FWHM $= 208 \pm 12$ km s$^{-1}$ straddling the systemic velocity and a narrow component with FWHM $= 44 \pm 3$ km s$^{-1}$ blueshifted by 160 km s$^{-1}$. The profile is significantly different from the known 22-GHz $6_{1,6} - 5_{2,3}$ maser which shows a broad profile redshifted by 193 km s$^{-1}$. The submillimeter maser is spatially unresolved with a synthesized beam of $0^{\prime \prime}.68 \times 0^{\prime \prime}.56$ and coincides with the continuum core position within 12 pc. These results indicate amplification of the continuum emission through high-temperature ($>1000$ K) and dense ($n({\rm H}_2{\rm O}) > 10^4$ cm$^{-3}$) molecular gas in front of the core

We measure the galaxy-ellipticity (GI) correlations for the Slogan Digital Sky Survey DR12 LOWZ and CMASS samples with the shape measurements from the DESI Legacy Imaging Surveys. We model the GI correlations in an N-body simulation with our recent accurate stellar-halo mass relation from the Photometric object Around Cosmic webs (PAC) method. The large data set and our accurate modeling turns out an accurate measurement of the alignment angle between central galaxies and their host halos. We find that the alignment of central elliptical galaxies with their host halos increases monotonically with galaxy stellar mass or host halo mass, which can be well described by a power law for the massive galaxies. We also find that central elliptical galaxies are more aligned with their host halos when they evolve to a lower redshift. In contrast, central disk galaxies are aligned with their host halos about 10 times more weakly in the GI correlation. These results have important implications for intrinsic alignment (IA) correction in weak lensing studies, IA cosmology, and theory of massive galaxy formation.

Methanol is the most complex molecule securely identified in interstellar ices and is a key chemical species for understanding chemical complexity in astrophysical environments. Important aspects of the methanol ice photochemistry are still unclear such as the branching ratios and photo-dissociation cross-sections at different temperatures and irradiation fluxes. This work aims at a quantitative agreement between laboratory experiments and astrochemical modelling of the CH3OH ice UV photolysis. This work also allows us to better understand which processes govern the methanol ice photochemistry present in laboratory experiments. We use ProDiMo to simulate the conditions of laboratory measurements. The simulations start with simple chemistry consisting only of methanol ice and helium to mimic the residual gas in the experimental chamber. A surface chemical network enlarged by photo-dissociation reactions is used to study the chemical reactions within the ice. Additionally, different surface chemistry parameters (surface competition, tunnelling, thermal diffusion and reactive desorption) are adopted to check those that reproduce the experimental results. The chemical models with ProDiMo can reproduce the methanol ice destruction via UV photodissociation at temperatures of 20, 30, 50 and 70 K as observed in the experiments. We note that the results are sensitive to different branching ratios after photolysis and to the mechanisms of reactive desorption. In the simulations of a molecular cloud at 20 K, we observed an increase in the methanol gas abundance of one order of magnitude, with a similar decrease in the solid-phase abundance. Comprehensive astrochemical models provide new insights into laboratory experiments as the quantitative understanding of the processes that govern the reactions within the ice. Ultimately, those insights can help to better interpret astronomical observations.

Present-day Mars is cold and dry, but mineralogical and morphological evidence shows that liquid-water existed on the surface of ancient Mars. In order to explain this evidence and assess ancient Mars's habitability, one must understand the size and composition of the ancient atmosphere. Here we place constraints on the ancient Martian atmosphere by modeling the coupled, self-consistent evolution of atmospheric CO2, N2, and Ar on Mars from 3.8 billion years ago (Ga) to the present. Our model traces the evolution of these species' abundances and isotopic composition caused by atmospheric escape, volcanic outgassing, and crustal interaction. Using a Markov-Chain Monte Carlo method to explore a plausible range of parameters, we find hundreds of thousands of model solutions that recreate the modern Martian atmosphere. These solutions indicate that Mars's atmosphere contained 0.3-1.5 bar CO2 and 0.1-0.5 bar N2 at 3.8 Ga. The global volume of deposited carbonates critically determines the ancient atmospheric composition. For example, a ~1 bar CO2 ancient atmosphere with 0.2-0.4 bar N2 requires ~0.9 bar CO2 deposited in carbonates primarily in open-water systems. With the joint analysis of C, N, and Ar isotopes, we refine the constraints on the relative strengths of outgassing and sputtering, leading to an indication of a reduced early mantle from which the outgassing is sourced. Our results indicate that a CO2-N2 atmosphere with a potential H2 component on ancient Mars is consistent with Mars's geochemical evolution and may explain the evidence for its past warm and wet climate.

We present an analysis of a volume-complete sample of 363 mid-to-late M dwarfs within 15 pc of the Sun with masses between 0.1 and 0.3 M$_\odot$ observed by TESS within Observation Sectors 1 to 42. The median mass of the stars in this sample is 0.17 M$_\odot$. We search the TESS 2-minute cadence light curves for transiting planets with orbital periods below 7 days using a modified Box-Least Squares (BLS) algorithm and recover all 6 known planets within the sample as well as a likely planet candidate orbiting LHS 475 (TESS Object of Interest 910.01). Each of these planets is consistent with a terrestrial composition, with planet radii ranging from 0.91 R$_\oplus$ to 1.31 R$_\oplus$. In addition, we perform a transit injection and recovery analysis for each of the 363 stars to characterize the transit detection sensitivity as a function of planet radius, insolation, and orbital period. We obtain a cumulative occurrence rate of $0.61^{+0.24}_{-0.19}$ terrestrial planets per M dwarf with radii above 0.5 R$_\oplus$ and orbital periods between 0.4-7 days. We find that for comparable insolations, planets larger than 1.5 R$_\oplus$ (sub-Neptunes and water worlds) are significantly less abundant around mid-to-late M dwarfs compared to earlier-type stars, while the occurrence rate of terrestrial planets is comparable to that of more massive M dwarfs. We estimate that overall, terrestrials outnumber sub-Neptunes around mid-to-late M dwarfs at a ratio of 14 to 1, in contrast to GK dwarfs where they are roughly equinumerous. We place a $1\sigma$ upper limit of 0.07 planets larger than 1.5 R$_\oplus$ per mid-to-late M dwarf, within the orbital period range of 0.5-7 days. We find evidence for a downturn in occurrence rates for planet radii below 0.9 R$_\oplus$, suggesting that Earth-sized and larger terrestrials may be more common around mid-to-late M dwarfs.

The amount of baryons hosted in the disks of galaxies is lower than expected based on the mass of their dark-matter halos and the fraction of baryon-to-total matter in the universe, giving rise to the so called galaxy missing-baryon problem. The presence of cool circum-galactic matter gravitationally bound to its galaxy's halo up to distances of at least ten times the size of the galaxy's disk, mitigates the problem but is far from being sufficient for its solution. It has instead been suggested, that the galaxy missing baryons may hide in a much hotter gaseous phase of the circum-galactic medium, possibly near the halo virial temperature and co-existing with the cool phase. Here we exploit the best available X-ray spectra of known cool circum-galactic absorbers of L$^*$ galaxies to report the first direct high-statistical-significance ($5.3-6.8\sigma$) detection of associated OVII and NVI absorption in the stacked XMM-Newton and Chandra spectra of three quasars. We show that these absorbers trace hot medium in the X-ray halo of these systems, at logT(in k)$\simeq 5.88-6.1$ K (comprising the halo virial temperature T$_{vir} = 10^6$ K). We estimate a mass of the X-ray halo M$_{hot-CGM}\simeq (1.4-1.6)\times 10^{11}$ M$_{\odot}$, corresponding, for these systems, to a galaxy missing baryon fraction $\xi = M_{hot-CGM}/M_{missing}\simeq 0.99-1.13$ and thus closing the galaxy baryon census in typical L$^*$ galaxies. Our measurements contribute significantly to the solution of the long-standing galaxy missing baryon problem and to the understanding of the continuous cycle of baryons in-and-out of galaxies throughout the life of the universe.

Acoustic wave heating is believed to contribute significantly to the missing energy input required to maintain the solar chromosphere in its observed state. We studied the propagation of waves above the acoustic cutoff in the upper photosphere into the chromosphere with ultraviolet and optical spectral observations interpreted through comparison with three dimensional radiative magnetohydrodynamic (rMHD) \emph{Bifrost} models to constrain the heating contribution from acoustic waves in the solar atmosphere. Sit-and-stare observations taken with the Interface Region Imaging Spectrograph (IRIS) and data from the Interferometric BIdimensional Spectrograph (IBIS) were used to provide the observational basis of this work. We compared the observations with synthetic observables derived from the Bifrost solar atmospheric model. Our analysis of the \emph{Bifrost} simulations show that internetwork and enhanced network regions exhibit significantly different wave propagation properties, which are important for the accurate wave flux estimates. The inferred wave energy fluxes based on our observations are not sufficient to maintain the solar chromosphere. We point out that the systematics of the modeling approaches in the literature lead to differences which could determine the conclusions of this type of studies, based on the same observations.

Details of the magnetic field in the quiet Sun chromosphere are key to our understanding of essential aspects of the solar atmosphere. We aim to determine the longitudinal magnetic field component (B_lon) of quiet Sun regions depending on their size. We estimated B_lon by applying the weak-field approximation (WFA) to high-spatial-resolution Ca II 854.2 nm data taken with the Swedish 1m Solar Telescope. Specifically, we analyzed the estimates inferred for different spectral ranges using the data at the original cadence and temporally integrated signals. The longitudinal magnetic field in each considered plasma structure correlates with its size. Using a spectral range restricted to the line core leads to chromospheric longitudinal fields varying from 50 G at the edges to 150-500 G at the center of the structure. These values increase as the spectral range widens due to the photospheric contribution. However, the difference between this contribution and the chromospheric one is not uniform for all structures. Small and medium-sized concentrations show a steeper height gradient in B_lon compared to their chromospheric values, so estimates for wider ranges are less trustworthy. Signal addition does not alleviate this situation as the height gradients in B_lon are consistent with time. Finally, despite the amplified noise levels that deconvolving processes may cause, data restored with the destretching technique show similar results, though are affected by smearing. We obtained B_lon estimates similar to those previously found, except for large concentrations and wide spectral ranges. In addition, we report a correlation between the height variation of B_lon compared to the chromospheric estimates and the concentration size. This correlation affects the difference between the photospheric and chromospheric magnetic flux values and the reliability of the estimates for wider spectral ranges.

We present an early analysis to search for the high redshift galaxies using the deepest public \emph{JWST} imaging to date; the NGDEEP field. This data consists of 6-band NIRCam imaging on the Hubble Ultra Deep Field-Par2, covering a total area of 6.8 arcmin$^{2}$. Based on our initial reduction of the first half of this survey, we reach 5 $\sigma$ depths up to mag = 29.5--29.9 between $1-5$~\textmu m. Such depths present an unprecedented opportunity to begin exploring the early Universe with \emph{JWST}. As such, we find high redshift galaxies in this field by examining the spectral energy distribution of these systems and present 18 new $z > 8$ systems identified using two different photometric redshift codes: \lephare\ and \eazy\, combined with other significant criteria. The highest redshift object in our sample is at $z=15.57^{+0.39}_{-0.38}$, which has a blue beta slope of $\beta=-3.25^{+0.41}_{-0.46}$ and a very low inferred stellar mass of M$_{*} = 10^{7.39}$~\solm\,. We also discover a series of faint, low-mass dwarf galaxies with M$_{*} < 10^{8.5}$~\solm at $z \sim 9$ that have blue colors and UV slopes. The structure of these galaxies is such that they all have very flat surface brightness profiles and small sizes of $< 1 \mathrm{kpc}$. We also compare our results to theory, finding no significant disagreement with some CDM based models. The discovery of these objects, most of which are low luminosity and inferred stellar mass, demonstrate the power of probing continuously deeper into the Universe. These observations will point the way to future observations of even deeper or similar deep but wider area surveys and critical need for \emph{JWST} deep fields to explore this aspect of the early Universe.

Gravitational waves from binary neutron star (BNS) mergers can constrain nuclear matter models predicting the neutron star's equation of state (EOS). Matter effects on the inspiral-merger signal are encoded in the multipolar tidal polarizability parameters, whose leading order combination is sufficient to capture to high accuracy the key features of the merger waveform (e.g.~the merger frequency). Similar EOS-insensitive relations exist for the post-merger signal and can be used to model the emission from the remnant. Several works suggested that the appearance of new degrees of freedom or phase transitions in high-density post-merger matter can be inferred by observing a violation of these EOS-insensitive relations. Here, we demonstrate a Bayesian method to test such an EOS-insensitive relation between the tidal polarizability parameters (or any other equivalent parameter) and the dominant post-merger frequency, using information either up to merger or from the post-merger signal. Technically, the method is similar to tests of General Relativity with binary black holes that verify the inspiral-merger-ringdown consistency. However, differently from the latter, BNS pre/post-merger consistency tests are conceptually less informative and they only address the consistency (or the breaking) of the assumed EOS-insensitive relation. Specifically, we discuss how such tests cannot conclusively discriminate between an EOS not respecting such relation and the appearance of new degrees of freedom (or phase transitions) in high-density matter.

We present closed-form solutions for plunging geodesics in the extended Kerr spacetime using Boyer-Lindquist coordinates. Our solutions directly solve for the dynamics of generic timelike plunges, we also specialise to the case of test particles plunging from a precessing innermost stable circular orbit (ISSO). We find these solutions in the form of elementary and Jacobi elliptic functions parameterized by Mino time. In particular, we demonstrate that solutions for the ISSO case can be determined almost entirely in terms of elementary functions, depending only on the spin parameter of the black hole and the radius of the ISSO. This extends recent work on the case of equatorial plunges from the innermost stable circular orbit. Furthermore, we introduce a new equation that characterizes the radial inflow from the ISSO to the horizon, taking into account the inclination. For ease of application, our results have been implemented in the KerrGeodesics package in the Black Hole Perturbation Toolkit.

We explore the dynamics of FLRW cosmologies which consist of dark matter, radiation and dark energy with a quadratic equation of state. Standard cosmological singularities arise due to energy conditions which are violated by dark energy, therefore we focus our analysis on non-singular bouncing and cyclic cosmologies, in particular focusing on the possibility of closed models always having a bounce for any initial conditions. We analyse the range of dynamical behaviour admitted by the system, and find a class of closed models that admit a non-singular bounce, with early- and late-time accelerated expansion connected by a decelerating phase. In all cases, we find the bouncing models are only relevant when dark matter and radiation appear at a certain energy scale, and so require a period such as reheating. We then investigate imposing an upper bound on the dark matter and radiation, such that their energy densities cannot become infinite. We find that bounces are always the general closed model, and a class of models exist with early- and late-time acceleration, connected by a decelerating phase. We also consider parameter values for the dark energy component, such that the discrepancy between the observed value of $\Lambda$ and the theoretical estimates of the contributions to the effective cosmological constant expected from quantum field theory would be explained. However, we find that the class of models left does not allow for an early- and late-time accelerated expansion, connected by a decelerating period where large-scale structure could form. Nonetheless, our qualitative analysis serves as a basis for the construction of more realistic models with realistic quantitative behaviour.

It is shown how self-reproduction can be easily avoided in the inflationary universe, even when inflation starts at Planck scales. This is achieved by a simple coupling of the inflaton potential with a mimetic field. In this case, the problem of fine-tuning of the initial conditions does not arise, while eternal inflation and the multiverse with all their widely discussed problems are avoided.

We perform a detailed comparison of the dynamics of cosmic string loops obtained in cosmological field theory simulations with their expected motion according to the Nambu-Goto action. We demonstrate that these loops follow the trajectories predicted within the NG effective theory except in regions of high curvature where energy is emitted from the loop in the form of massive radiation. This energy loss continues for all the loops studied in this simulation until they self-intersect or become small enough that they annihilate and disappear well before they complete a single oscillation. We comment on the relevance of this investigation to the interpretation of the results from cosmological field theory simulations as well as their extrapolation to a cosmological context.

Triple-$\alpha$ reaction rates have been determined well with the sequential process via the narrow resonances. However, the direct triple-$\alpha$ process at off-resonant energies still remains in unsolved problems. In the present report, the direct triple-$\alpha$ contribution is estimated with a non-adiabatic method, and it is found to be 10$^{-15}$--10$^{-3}$ pb order in photodisintegration cross sections of $^{12}$C(2$^+_1 \rightarrow$ 0$^+$) for $0.15 < E < 0.35$ MeV. This is far below the values predicted by the recent adiabatic models. In spite of the large difference, the derived rates are found to be concordant with NACRE at the helium burning temperatures.

String cosmology models predict a relic background of gravitational-wave (GW) radiation in the early universe. The GW energy spectrum of radiated power increases rapidly with the frequency, and therefore it becomes a potential and meaningful observation object for high-frequency GW detector. We focus on the stochastic background generated by superinflation in string theory and search for such signal in the observing data of Advanced LIGO and Virgo O1$\sim$O3 runs in a Bayesian framework. We do not find the existence of the signal, and thus put constraints on the GW energy density. Our results indicate that at $f=100\,\text{Hz}$, the fractional energy density of GW background is less than $1.7\times10^{-8}$ and $2.1\times10^{-8}$ for dilaton-string and dilaton only cases respectively, and further rule out the parameter space restricted by the model itself due to the non-decreasing dilaton and stable cosmology background ($\beta$ bound).

The equation of state (EOS) of nuclear dense matter plays a crucial role in many astrophysical phenomena associated with neutron stars (NSs). Fluid oscillations are one of the most fundamental properties therein. NSs support a family of gravity $g$-modes, which are related to buoyancy. We study the gravity $g$-modes caused by composition gradient and density discontinuity in the framework of pseudo-Newtonian gravity. The mode frequencies are calculated in detail and compared with Newtonian and general-relativistic (GR) solutions. We find that the $g$-mode frequencies in one of the pseudo-Newtonian treatments can approximate remarkably well the GR solutions, with relative errors in the order of $1\%$. Our findings suggest that, with much less computational cost, pseudo-Newtonian gravity can be utilized to accurately analyze oscillation of NSs constructed from an EOS with a first-order phase transition between nuclear and quark matter, as well as to provide an excellent approximation of GR effects in core-collapse supernova (CCSN) simulations.

A portable monoenergetic 24 keV neutron source based on the $^{124}$Sb-$^9$Be photoneutron reaction and an iron filter has been constructed and characterized. The coincidence of the neutron energy from SbBe and the low interaction cross-section with iron (mean free path up to 29 cm) makes pure iron specially suited to shield against gamma rays from $^{124}$Sb decays while letting through the neutrons. To increase the $^{124}$Sb activity and thus the neutron flux, a $>$1 GBq $^{124}$Sb source was produced by irradiating a natural Sb metal pellet with a high flux of thermal neutrons in a nuclear reactor. The design of the source shielding structure makes for easy transportation and deployment. A hydrogen gas proportional counter is used to characterize the neutrons emitted by the source and a NaI detector is used for gamma background characterization. At the exit opening of the neutron beam, the characterization determined the neutron flux in the energy range 20-25 keV to be 5.36$\pm$0.20 neutrons per cm$^2$ per second and the total gamma flux to be 213$\pm$6 gammas per cm$^2$ per second (numbers scaled to 1 GBq activity of the $^{124}$Sb source). A liquid scintillator detector is demonstrated to be sensitive to neutrons with incident kinetic energies from 8 to 17 keV, so it can be paired with the source as a backing detector for neutron scattering calibration experiments. This photoneutron source provides a good tool for in-situ low energy nuclear recoil calibration for dark matter experiments and coherent elastic neutrino-nucleus scattering experiments.

Ultralight Axion Like Particle (ALP) can mediate a long range monopole-dipole macroscopic force between Earth and Sun if Earth is treated as a polarized source. There are about $10^{42}$ number of polarized electrons in Earth due to the presence of the geomagnetic field. The monopole-dipole interactions between electrons in Earth and nucleons in Sun can influence the perihelion precession of Earth, gravitational light bending and Shapiro time delay. The contribution of monopole-dipole potential is limited to be no larger than the measurement uncertainty. We obtain the first bound on monopole-dipole strength from single astrophysical observations. The perihelion precession of Earth puts the stronger bound on monopole-dipole coupling strength as $g_Sg_P\lesssim 3.62\times 10^{-11}$ for the ALP of mass $m_a\lesssim 1.35\times 10^{-18}~\rm{eV}$

The impacts of various symmetry energy parameters on the properties of neutron stars (NSs) have been recently investigated, and the outcomes are at variance, as summarized in Table III of Phys. Rev. D 106, 063005 (2022). We have systematically analyzed the correlations of slope and curvature parameters of symmetry energy at the saturation density ($\rho_0=0.16 \text{fm}^{-3}$) with the tidal deformability and stellar radius of non-spinning neutron stars in the mass range of $1.2 - 1.6 M_\odot$ using a large set of minimally constrained equations of state (EoSs). The EoSs at low densities correspond to the nucleonic matter and are constrained by empirical ranges of a few low-order nuclear matter parameters from the finite nuclei data and the pure neutron matter EoS from chiral effective field theory. The EoSs at high densities ($\rho > 1.5 - 2\rho_0$) are obtained by a parametric form for the speed of sound that satisfies the causality condition. Several factors affecting the correlations between the NS properties and the individual symmetry energy parameters usually encountered in the literature are considered. These correlations are quite sensitive to the choice of the distributions of symmetry energy parameters and their interdependence. But, variations of NS properties with the pressure of $\beta -$ equilibrated matter at twice the saturation density remain quite robust which maybe due to the fact that the pressure depends on the combination of multiple nuclear matter parameters that describe the symmetric nuclear matter as well as the density dependence of the symmetry energy. Our results are practically insensitive to the behavior of EoS at high densities.

In this work, we consider a hilltop version of the (supersymmetric) sneutrino hybrid inflation where the right-handed sneutrino field plays the role of the inflaton field. This model is a type III hilltop inflation that can produce a spectral index $n_s=0.96$ which fits perfectly to experimental observations without fine-tuning of parameters. We also briefly consider nonthermal leptogenesis via the decay of the right-handed sneutrino inflaton field after inflation.

The emerging field of quantum sensors and electronics for fundamental physics is introduced, emphasising the role of thin-film superconducting devices. Although the next generation of ground-based and space-based experiments requires the development of advanced technology across the whole of the electromagnetic spectrum, this article focuses on ultra-low-noise techniques for radio to far-infrared wavelengths, where existing devices fall short of theoretical limits. Passive circuits, detectors and amplifiers are described from classical and quantum perspectives, and the sensitivities of detector-based and amplifier-based instruments discussed. Advances will be achieved through refinements in existing technology, but innovation is essential. The needed developments go beyond engineering and relate to theoretical studies that bring together concepts from quantum information theory, quantum field theory, classical circuit theory, and device physics. This article has been written to introduce graduate-level scientists to quantum sensor physics, rather than as a formal review.

We describe the earliest measurements of the DT fusion cross section commissioned by the Manhattan Project, first at Purdue University in 1943 and then at Los Alamos 1945-6 and later, in 1951-2. The Los Alamos measurements led to the realization that a 3/2$^+$ resonance in the DT system enhances the fusion cross section by a factor of one hundred at energies relevant to applications. This was a transformational discovery, making the quest for terrestrial fusion energy possible. The earliest measurements were reasonably accurate given the technology of the time and the scarcity of tritium, and were quickly improved to provide cross section data accurate to just a few percent. We provide a previously-unappreciated insight: that DT fusion was first reported in Ruhlig's 1938 University of Michigan experiment and likely influenced Konopinski in 1942 to suggest its usefulness for thermonuclear technologies. We report on preliminary work to repeat the 1938 measurement, and our simulations of that experiment. We also present some work by Fermi, from his 1945 Los Alamos lectures, showing that he used the S-factor concept about a decade before it was introduced by nuclear astrophysicists.