Papers for Friday, Oct 21 2022

It is reported that the Large High Altitude Air Shower Observatory (LHAASO) observed thousands of very-high-energy photons up to $\sim$18 TeV from GRB 221009A. We study the survival rate of these photons by considering the fact that they are absorbed by the extragalactic background light. By performing a set of $10^6$ Monte-Carlo simulations, we explore the parameter space allowed by current observations and find that the probability of predicting that LHAASO observes at least one photons of 18 TeV from GRB 221009A within 2000 seconds is 4--5\%. Hence, it is still possible for the standard physics to interpret LHAASO's observation in the energy range of several TeV. Our method can be straightforwardly generalized to study more data sets of LHAASO and other experiments in the future.

Young eclipsing binaries (EBs) are powerful probes of early stellar evolution. Current models are unable to simultaneously reproduce the measured and derived properties that are accessible for EB systems (e.g., mass, radius, temperature, luminosity). In this study we add a benchmark EB to the pre-main sequence population with our characterization of TOI 450 (TIC 77951245). Using \gaia\ astrometry to identify its comoving, coeval companions, we confirm TOI 450 is a member of the $\sim$40 Myr Columba association. This eccentric ($e=0.2969$), equal-mass ($q=1.000$) system provides only one grazing eclipse. Despite this, our analysis achieves the precision of a double-eclipsing system by leveraging information in our high-resolution spectra to place priors on the surface-brightness and radius ratios. We also introduce a framework to include the effect of star spots on the observed eclipse depths. Multi-color eclipse light curves play a critical role in breaking degeneracies between the effects of star spots and limb-darkening. Including star spots reduces the derived radii by $\sim$2\% from a unspotted model ($>2\sigma$) and inflates the formal uncertainty in accordance with our lack of knowledge regarding the star spot orientation. We derive masses of 0.1768($\pm$0.0004) and 0.1767($\pm$0.0003) $M_\odot$, and radii of 0.345($\pm$0.006) and 0.346($\pm$0.006) $R_\odot$ for the primary and secondary, respectively. We compare these measurements to multiple stellar evolution isochones, finding good agreement with the association age. The MESA MIST and SPOTS ($f_{\rm s}=0.17$) isochrones perform the best across our comparisons, but detailed agreement depends heavily on the quantities being compared.

Strong lensing at the galaxy-scale is useful for numerous applications in Astrophysics and Cosmology. Some of the principal applications are studying the mass structure of elliptical galaxies, their formation and evolution, constraining the stellar initial mass function, and measuring cosmological parameters. Since the first discovery of a galaxy-scale strong lens in the eighties, this field has come a long way in terms of data quality and availability, and techniques to model the data. In this review article, we describe the most common methodologies to model lensing observables of galaxy-scale strong lenses, especially the imaging data due to it being the most available and informative source of lensing observables. We review the main results from the literature in astrophysical and cosmological applications of galaxy-scale strong lenses. We also discuss the current limitations of the data and methodologies, and provide a future outlook of the expected development and improvements in both aspects in the near future.

2S 0921$-$630 is an eclipsing low-mass X-ray binary (LMXB) with an orbital period of $\sim$ 9 days. Past X-ray observations have revealed that 2S 0921$-$630 has an extended accretion disk corona (ADC), from which most of the X-rays from the system are emitted. We report the result of our Suzaku archival data analysis of 2S 0921$-$630. The average X-ray spectrum is reproduced with a blackbody emission ($kT_{\rm bb} \sim 0.3$ keV) Comptonized by optically-thick gas (``Compton cloud''; optical depth $\tau \sim 21$) with a temperature of $\sim 2$ keV, combined with thirteen emission lines. We find that most of the emission lines correspond to highly ionized atoms: O, Ne, Mg, Si, S, Ar, and Fe. A K$\alpha$ emission line and an absorption edge of semi-neutral iron (Fe I -- XVII) are also detected. The semi-neutral iron K$\alpha$ line is significantly broad with a width of $0.11 \pm 0.02$ keV in sigma, which corresponds to the Doppler broadening by the Kepler motion at a radius of $\sim 10^9$ cm. We suggest that the observed semi-neutral iron line originates at the inner part of the accretion disk in the immediate outside of the Compton cloud, i.e., the Compton cloud may have a radius of $\sim 10^9$ cm.

The unprecedented number of gravitational lenses expected from new-generation facilities such as the ESA Euclid telescope and the Vera Rubin Observatory makes it crucial to rethink our classical approach to lens-modelling. In this paper, we present LeMoN (Lens Modelling with Neural networks): a new machine-learning algorithm able to analyse hundreds of thousands of gravitational lenses in a reasonable amount of time. The algorithm is based on a Bayesian Neural Network: a new generation of neural networks able to associate a reliable confidence interval to each predicted parameter. We train the algorithm to predict the three main parameters of the Singular Isothermal Ellipsoid model (the Einstein radius and the two components of the ellipticity) by employing two simulated datasets built to resemble the imaging capabilities of the Hubble Space Telescope and the forthcoming Euclid satellite. In this work, we assess the accuracy of the algorithm and the reliability of the estimated uncertainties by applying the network to several simulated datasets of 10.000 images each. We obtain accuracies comparable to previous studies present in the current literature and an average modelling time of just 0.5s per lens. Finally, we apply the LeMoN algorithm to a pilot dataset of real lenses observed with HST during the SLACS program, obtaining unbiased estimates of their SIE parameters. The code is publicly available on GitHub (https://github.com/fab-gentile/LeMoN).

The Hawaii Survey Field SSA22 is the fourth deepest Chandra X-ray field. To allow for the fullest exploration of this field, we present new optical spectroscopy from Keck/DEIMOS and Keck/LRIS, which, in combination with the literature, brings the spectroscopic completeness of the 2--8 keV sample to 62%. We also make optical spectral classifications and estimate photometric redshifts for the sources without spectroscopic redshifts. We then determine hard X-ray luminosity functions (XLFs) for the full sample of Active Galactic Nuclei (AGNs), as well as for the broad-line AGNs (BLAGNs) and the non-BLAGNs separately. Our XLF for the full sample is in good agreement with the literature, showing relatively strong evolution over the redshift range $0.25\le z < 4$. The XLFs for the BLAGNs and the non-BLAGNs imply distinct evolution with redshift, with BLAGNs becoming increasingly dominant at higher redshifts and X-ray luminosities.

Understanding stellar composition is fundamental not only to our comprehension of the galaxy, especially chemical evolution, but it can also shed light on the interior structure and mineralogy of exoplanets, which are formed from the same material as their host stars. Unfortunately, the underlying mathematics describing stellar mass fractions and stellar elemental abundances is difficult to parse, fragmented across the literature, and contains vexing omissions that makes any calculation far from trivial, especially for non-experts. In this treatise, we present clear mathematical formalism and clarification of inherent assumptions and normalizations within stellar composition measurements, which facilitates the conversion from stellar mass fractions to elemental abundances to molar ratios, including error propagation. We also provide an example case study of HIP 544 to further illustrate the provided equations. Given the important chemical association between stars, as well as the interdisciplinary relationship between stars and their planets, it is vital that stellar mass fractions and abundance data be more transparent and accessible to people within different sub-fields and scientific disciplines.

We perform a search for galaxy-galaxy strong lens systems using a convolutional neural network (CNN) applied to imaging data from the first public data release of the DECam Local Volume Exploration Survey (DELVE), which contains $\sim 520$ million astronomical sources covering $\sim 4,000$ $\mathrm{deg}^2$ of the southern sky to a $5\sigma$ point-source depth of $g=24.3$, $r=23.9$, $i=23.3$, and $z=22.8$ mag. Following the methodology of similar searches using DECam data, we select a catalog of $\sim 11$ million extended astronomical sources with $16 < g < 22$, $17.2 < r < 22$, $15 < i < 21$, $0 < g - i < 3$, and $-0.2 < g - r < 1.75$. After scoring with our CNN, the highest scoring 50,000 images were visually inspected and assigned a score on a scale from 0 (definitely not a lens) to 3 (very probable lens). We present a list of 617 strong lens candidates, 599 of which are previously unreported. We additionally highlight 8 potential quadruply lensed quasars from this sample. Due to the location of our search footprint in the northern Galactic cap ($b > 10$ deg) and southern celestial hemisphere (${\rm Dec.}<0$ deg), our candidate list has little overlap with other existing ground-based searches. Where our search footprint does overlap with other searches, we find a significant number of high-quality candidates which were previously unidentified, indicating a degree of orthogonality in our methodology. We categorize our candidates using their human-assigned scores and report properties including apparent magnitude and Einstein radius estimated from the image separation.

Recent observations have shown that the environmental quenching of galaxies at z ~ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been relatively few comparisons at higher redshifts to date. Here we confront three state-of-the-art suites of simulations (BAHAMAS+MACSIS, EAGLE+Hydrangea, IllustrisTNG) with state-of-the-art observations of the field and cluster environments from the COSMOS/UltraVISTA and GOGREEN surveys, respectively, at z ~ 1 to assess the realism of the simulations and gain insight into the evolution of environmental quenching. We show that while the simulations generally reproduce the stellar content and the stellar mass functions of quiescent and star-forming galaxies in the field, all the simulations struggle to capture the observed quenching of satellites in the cluster environment, in that they are overly efficient at quenching low-mass satellites. Furthermore, two of the suites do not sufficiently quench the highest-mass galaxies in clusters, perhaps a result of insufficient feedback from AGN. The origin of the discrepancy at low stellar masses (Mstar <~ 1E10 Msun), which is present in all the simulations in spite of large differences in resolution, feedback implementations, and hydrodynamical solvers, is unclear. The next generation of simulations, which will push to significantly higher resolution and also include explicit modelling of the cold interstellar medium, may help to shed light on the low-mass tension.

The radiosotope $^{44}$Ti is produced through $\alpha$-rich freezeout and explosive helium burning in type Ia supernovae (SNe Ia). In this paper, we discuss how the detection of $^{44}$Ti, either through late-time light curves of SNe Ia, or directly via gamma rays, can uniquely constrain the origin of SNe Ia. In particular, building upon recent advances in the hydrodynamical simulation of helium-ignited double white dwarf binaries, we demonstrate that the detection of $^{44}$Ti in a nearby SN Ia or in a young galactic supernova remnant (SNR) can discriminate between the double-detonation and double-degenerate channels of sub-Chandrasekhar (sub-$M_{\rm Ch}$) and near-Chandrasekhar (near-$M_{\rm Ch}$) SNe Ia. In addition, we predict that the late-time light curves of calcium-rich transients are entirely dominated by $^{44}$Ti.

Searching for Earth-sized planets in data from Kepler's extended mission (K2) is a niche that still remains to be fully exploited. The TFAW survey is an ongoing project that aims to re-analyze all light curves in K2 C1-C8 and C12-C18 campaigns with a wavelet-based detrending and denoising method, and the period search algorithm TLS to search for new transit candidates not detected in previous works. We have analyzed a first subset of 24 candidate planetary systems around relatively faint host stars (10.9 < $K_{p}$ < 15.4) to allow for follow-up speckle imaging observations. Using VESPA and TRICERATOPS, we statistically validate six candidates orbiting four unique host stars by obtaining false-positive probabilities smaller than 1% with both methods. We also present 13 vetted planet candidates that might benefit from other, more precise follow-up observations. All of these planets are sub-Neptune-sized, with two validated planets and three candidates with sub-Earth sizes, and have orbital periods between 0.81 and 23.98 days. Some interesting systems include two ultra-short-period planets, three multi-planetary systems, three sub-Neptunes that appear to be within the small planet Radius Gap, and two validated and one candidate sub-Earths (EPIC 210706310, EPIC 210768568, and EPIC 246078343) orbiting metal-poor stars.

Time-domain analysis of an archival XMM-Newton observation unveiled a very unusual variability pattern in the soft X-ray emission of PSR J1311-3430, a black widow millisecond pulsar in a tight binary (P_B=93.8 min) with a very low-mass (M~0.01 Msun) He companion star, known to show flaring emission in the optical and in the X-rays. A series of six pulses with a regular recurrence time of ~124 min is apparent in the 0.2-10 keV light curve of the system, also featuring an initial, bright flare and a quiescent phase lasting several hours. The X-ray spectrum does not change when the pulses are seen and is consistent with a power law with photon index Gamma~1.6, also describing the quiescent emission. The peak luminosity of the pulses is of several 10^32 erg/s. Simultaneous observations in the U band with the Optical Monitor onboard XMM and in the g' band from the Las Cumbres Observatory do not show any apparent counterpart of the pulses and only display the well-known orbital modulation of the system. We consider different hypotheses to explain the recurrent pulses: we investigate their possible analogy with other phenomena already observed in this pulsar and in similar systems and we also study possible explanations related to the interaction of the energetic pulsar wind with intrabinary material, but we found none of these pictures to be convincing. We identify simultaneous X-ray observations and optical spectroscopy as a possible way to constrain the nature of the phenomenon.

We present an expansion of FLEET, a machine learning algorithm optimized to select transients that are most likely to be tidal disruption events (TDEs). FLEET is based on a random forest algorithm trained on the light curves and host galaxy information of 4,779 spectroscopically classified transients. For transients with a probability of being a TDE, \ptde$>0.5$, we can successfully recover TDEs with a $\approx40$\% completeness and a $\approx30$\% purity when using the first 20 days of photometry, or a similar completeness and $\approx50$\% purity when including 40 days of photometry. We find that the most relevant features for differentiating TDEs from other transients are the normalized host separation, and the light curve $(g-r)$ color during peak. Additionally, we use FLEET to produce a list of the 39 most likely TDE candidates discovered by the Zwicky Transient Facility that remain currently unclassified. We explore the use of FLEET for the Legacy Survey of Space and Time on the Vera C. Rubin Observatory (\textit{Rubin}) and the \textit{Nancy Grace Roman Space Telescope} (\textit{Roman}). We simulate the \textit{Rubin} and \textit{Roman} survey strategies and estimate that $\sim 10^4$ TDEs could be discovered every year by \textit{Rubin}, and $\sim200$ TDEs per year by \textit{Roman}. Finally, we run FLEET on the TDEs in our \textit{Rubin} survey simulation and find that we can recover $\sim 30$\% of those at a redshift $z <0.5$ with \ptde$>0.5$. This translates to $\sim3,000$ TDEs per year that FLEET could uncover from \textit{Rubin}. FLEET is provided as a open source package on GitHub https://github.com/gmzsebastian/FLEET

In November 2019 we began operating FLEET (Finding Luminous and Exotic Extragalactic Transients), a machine learning algorithm designed to photometrically identify Type I superluminous supernovae (SLSNe) in transient alert streams. Using FLEET, we spectroscopically classified 21 of the 50 SLSNe identified worldwide between November 2019 and January 2022. Based on our original algorithm, we anticipated that FLEET would achieve a purity of about 50\% for transients with a probability of being a SLSN, \pslsn$>0.5$; the true on-sky purity we obtained is closer to 80\%. Similarly, we anticipated FLEET could reach a completeness of about 30\%, and we indeed measure an upper limit on the completeness of $\approx 33$\%. Here, we present FLEET 2.0, an updated version of FLEET trained on 4,780 transients (almost 3 times more than in FLEET 1.0). FLEET 2.0 has a similar predicted purity to FLEET 1.0, but outperforms FLEET 1.0 in terms of completeness, which is now closer to $\approx 40$\% for transients with \pslsn$>0.5$. Additionally, we explore possible systematics that might arise from the use of FLEET for target selection. We find that the population of SLSNe recovered by FLEET is mostly indistinguishable from the overall SLSN population, in terms of physical and most observational parameters. We provide FLEET as an open source package on GitHub https://github.com/gmzsebastian/FLEET

Double white dwarfs (DWDs) will be the most numerous gravitational-wave (GW) sources for the Laser Interferometer Space Antenna (LISA). Most of the Galactic DWDs will be unresolved and will superpose to form a confusion noise foreground, the dominant LISA noise source around $\sim 0.5\mathrm{-}3\,\mathrm{mHz}$. A small fraction of these sources will stand out from the background and be individually detectable. Uniquely among GW sources, a handful of these binaries will be known in advance from electromagnetic observations and will be guaranteed sources of detectable GWs in the LISA band; these are known as verification binaries (VBs). High-cadence photometric surveys are continuously discovering new VB systems, and their number will continue to grow ahead of the launch of LISA. We analyse, in a fully Bayesian framework, all the currently known VBs with the latest design requirements for the LISA mission. We explore what can be expected from GW observations, both alone and in combination with electromagnetic observations, and estimate the VB's time to detection in the early months of LISA operations. We also show how VBs can be analysed in the (realistic) case where their GW signals compete with many other unknown binary signals (both resolved and unresolved) from the Galactic population of DWDs.

We have identified 7 objects in Gaia Early DR3 associated with 6 International Pulsar Timing Array (PTA) pulsars. We combine both pulsar timing based parallax distance measurements with Gaia parallaxes to their companions, in an effort to improve distance measurements. We confirm previous cross-match findings for pulsars J0437-4715, J1012+5307, J1024-0719, J1732-5049, and J1843-1113. We no longer find the companion to J1949+3106 even though it was detected with a $> 3\sigma$ confidence in Gaia DR2. It is therefore likely that any definite cross-match must be made with $>3\sigma$ confidence, casting further doubt on companion identified with J1843-1113. We also report a low signal-to-noise detection of two candidates matched with J1747-4036 which we will continue to monitor. Encouragingly, the average signal-to-noise increase between Gaia DR2 and Gaia EDR3 is $\sim 50\%$, and the errors on the binary distances have improved by $10\%$.

We report the detection of three large millimeter flaring events from the nearby Sun-like, $\epsilon$ Eridani, found in archival ALMA 12m and ACA observations at 1.33 mm taken from 2015 January 17-18 and 2016 October 24-November 23, respectively. This is the first time that flares have been detected from a Sun-like star at millimeter wavelengths. The largest flare among our data was detected in the ALMA observations on 2015 January 17 from 20:09:10.4-21:02:49.3 (UTC) with a peak flux density of 28 $\pm$ 7 mJy and a duration of 9 sec. The peak brightness of the largest flare is $ 3.4 \pm 0.9 \times 10^{14}$ erg s$^{-1}$Hz$^{-1}$, a factor of $>50\times$ times brighter than the star's quiescent luminosity and $>10\times$ brighter than solar flares observed at comparable wavelengths. We find changes in the spectral index (F$_\nu\propto\nu^\alpha$) at the flare peak, with $\alpha$ = 1.81 $\pm$ 1.94 and a lower limit on the fractional linear polarization $|Q/I| = $ 0.08 $\pm$ 0.12. This positive spectral index is more similar to millimeter solar flares, differing from M dwarf flares also detected at millimeter wavelengths that exhibit steeply negative spectral indices.

The origin of ultra-high energy cosmic rays (UHECRs) is one of the most mystifying issues in astroparticle physics. It has been suggested that gamma-ray bursts (GRBs) are excellent acceleration sites for cosmic rays. The propagation of UHECRs from the GRB host galaxy to the Earth should generate delayed secondary photons and neutrinos. Here we present a dedicated search for delayed UHECR and neutrino emission centered around the position of nearby GRB 980425/SN 1998bw. Located at a distance of 36.9 Mpc, GRB 980425/SN 1998bw is well within the Greisen-Zatsepin-Kuzmin (GZK) distance horizon. We find no evidence for UHECR or neutrino clustering around the GRB 980425/SN 1998bw position between 2004 and 2020. Under ideal propagation conditions, we propose that it might be possible to detect an excess from delayed UHECRs around GRB 980425/SN 1998bw within the next 100 years if the intergalactic magnetic field (IGMF) strength is $B \leq 3 \times 10^{-13}$ G.

The Frequency Agile Solar Radiotelescope (FASR) has been strongly endorsed as a top community priority by both Astronomy & Astrophysics Decadal Surveys and Solar & Space Physics Decadal Surveys in the past two decades. Although it was developed to a high state of readiness in previous years (it went through a CATE analysis and was declared ``doable now"), the NSF has not had the funding mechanisms in place to fund this mid-scale program. Now it does, and the community must seize this opportunity to modernize the FASR design and build the instrument in this decade. The concept and its science potential have been abundantly proven by the pathfinding Expanded Owens Valley Solar Array (EOVSA), which has demonstrated a small subset of FASR's key capabilities such as dynamically measuring the evolving magnetic field in eruptive flares, the temporal and spatial evolution of the electron energy distribution in flares, and the extensive coupling among dynamic components (flare, flux rope, current sheet). The FASR concept, which is orders of magnitude more powerful than EOVSA, is low-risk and extremely high reward, exploiting a fundamentally new research domain in solar and space weather physics. Utilizing dynamic broadband imaging spectropolarimetry at radio wavelengths, with its unique sensitivity to coronal magnetic fields and to both thermal plasma and nonthermal electrons from large flares to extremely weak transients, the ground-based FASR will make synoptic measurements of the coronal magnetic field and map emissions from the chromosphere to the middle corona in 3D. With its high spatial, spectral, and temporal resolution, as well as its superior imaging fidelity and dynamic range, FASR will be a highly complementary and synergistic component of solar and heliospheric capabilities needed for the next generation of solar science.

Multiply lensed sources experience a relative time delay in the arrival of photons. This effect can be used to measure absolute distances and the Hubble constant ($H_0$) and is known as time-delay cosmography. The methodology is independent of the local distance ladder and early-universe physics and provides a precise and competitive measurement of $H_0$. With upcoming observatories, time-delay cosmography can provide a 1% precision measurement of $H_0$ and can decisively shed light on the current reported 'Hubble tension'. This paper presents the theoretical background and the current techniques applied for time-delay cosmographic studies and the measurement of the Hubble constant. The paper describes the challenges and systematics in the different components of the analysis and strategies to mitigate them. The current measurements are discussed in context and the opportunities with the anticipated data sets in the future are laid out.

We present star formation rates based on cold and ionized gas measurements of Mrk 266 (NGC 5256), a system composed of two colliding gas-rich galaxies, each hosting an active galactic nucleus. Using $^{12}$CO (1-0) observations with the Combined Array for Research in Millimeter-Wave Astronomy (CARMA), we find a total H$_2$ mass in the central region of $1.1\pm0.3\times10^{10}$ $M_\odot$ which leads to a possible future star formation rate of $25\pm10 M_\odot$ yr$^{-1}$. With the Fourier Transform Spectrograph (SITELLE) on the Canada-France-Hawaii Telescope, we measure an integrated H$\alpha$ luminosity and estimate a present-day star formation rate of $15\pm2 M_\odot$ yr$^{-1}$ in the core of the system (avoiding the two active nuclei). These results confirm that Mrk 266 is an intermediate stage merger with a relatively high recent star formation rate and enough molecular gas to sustain it for a few hundred million years. Inflowing gas associated with the merger may have triggered both the starburst episode and two AGN but the two galaxy components differ: the region around the SW nucleus appears to be more active than the NE nucleus, which seems relatively quiet. We speculate that this difference may originate in the properties of the interstellar medium in the two systems.

We numerically model evolution of magnetic fields inside a neutron star under the influence of ambipolar diffusion in the weak-coupling mode in the one-fluid MHD approximation. Our simulations are three-dimensional and performed in spherical coordinates. Our model covers the neutron star core and includes crust where the magnetic field decay is due to Ohmic decay. We discover an instability of poloidal magnetic field under the influence of ambipolar diffusion. This instability develops in the neutron star core and grows on a timescale of 0.2 dimensionless times, reaching saturation by 2 dimensionless times. The instability leads to formation of azimuthal magnetic field with azimuthal wavenumber $m=14$ (at the moment of saturation) which keeps merging and reaches $m=4$ by 16 dimensionless times. Over the course of our simulations (16 dimensionless times) the surface dipolar magnetic field decays, reaching 20 percent of its original value and keeps decaying. The decay timescale for the total magnetic energy is six dimensionless times. The ambipolar diffusion induces electric currents in the crust where these currents dissipate efficiently. Strong electric currents in the crust lead to heating, which could correspond to luminosities of $\approx 10^{29}$ erg s$^{-1}$ during hundreds of Myrs for an initial magnetic field of $10^{14}$ G. Ambipolar diffusion leads to formation of small-scale magnetic fields at the neutron star surface.

We present mosaic images of the Small Magellanic Cloud (SMC) observed with the Spitzer IRAC 3.6 $\mu$m and 4.5 $\mu$m bands over two epochs, 2017 August 25 to 2017 September 13, and 2017 November 24 to 2018 February 12. The survey region comprises $\sim$30 square degrees covering the SMC and the Bridge to the Large Magellanic Cloud. The region is covered by 52 $\sim$1$.\!\!^\circ$1$\times$1$.\!\!^\circ$1 tiles, with each tile including images in each band for both separate and combined epochs. The mosaics are made in individual tangent projections in J2000 coordinates. The angular pixel size is 0$.\!\!^{\prime\prime}$6 with a resolution (FWHM) of $\sim$2$.\!\!^{\prime\prime}$0. We describe processing to correct or mitigate residual artifacts and remove background discontinuities. The mosaic images are publicly available at the Infrared Science Archive (IRSA).

Synchronous binary asteroids can experience libration about their tidally-locked equilibrium, which will result in energy dissipation. This is an important topic to the Asteroid Impact and Deflection Assessment, where excitation caused by the DART kinetic impact in the Didymos binary asteroid system may be reduced through dissipation before Hera arrives to survey the effects of the impact. We develop a numeric model for energy dissipation in binary asteroids to explore how different system configurations affect the rate of energy dissipation. We find tumbling within the synchronous state eliminates a systematic trend in libration damping on short timescales (several years), but not over long times (hundreds of years) depending on the material conditions. Furthermore, damping of libration, eccentricity, and fluctuations in the semimajor axis are primarily dependent on the stiffness of the secondary, whereas the semimajor axis secular expansion rate is dictated by the stiffness of the primary, as expected. Systems experiencing stable planar libration in the secondary can see a noticeable reduction in libration amplitude after only a few years depending on the stiffness of the secondary, and thus dissipation should be considered during Hera's survey of Didymos. For a very dissipative secondary undergoing stable libration, Hera may be able to calculate the rate of libration damping in Dimorphos and therefore constrain its tidal parameters.

Developed right at the beginning of the space age in the 1940s, the proportional counter was the first detector used in X-ray astronomy and stayed its workhorse for almost four decades. Although the principle of such a detector seems to be rather simple, over time it underwent considerable performance improvements and the lifetime under orbital conditions could be extended tremendously. Particularly the invention of position-sensitive proportional counters provided new and sophisticated methods to discriminate background and thus enabled observations of much weaker sources. A leap forward in position resolution was achieved with the advent of microchannel plate (MCP) detectors in the 1970s. In contrary to gas filled detectors, they provide no considerable energy resolution but feature spatial resolutions reaching down to a few ten micrometer, fitting ideally the angular resolution of the novel grazing incidence imaging X-ray telescopes upcoming at that time. Even today, both types of detectors are still relevant in space-based astronomy. However, in case of MCPs new developments focus on the far and extreme ultraviolet wavelength range, while the Chandra X-ray observatory is most likely the last mission applying this technology for X-rays. In contrast, compact detectors with gas electron multiplier (GEM) foils and micropattern readout are currently under heavy development for the soft X-ray range, since they allow for the first time to measure polarization in X-rays over a broad energy range. This chapter presents the principles of proportional counters and MCP detectors, highlights the respective performance characteristics, and summarizes their most important applications in X-ray astronomy.

We continue examining statistical data assimilation (SDA), an inference methodology, to infer solutions to neutrino flavor evolution, for the first time using real - rather than simulated - data. The model represents neutrinos streaming from the Sun's center and undergoing a Mikheyev-Smirnov-Wolfenstein (MSW) resonance in flavor space, due to the radially-varying electron number density. The model neutrino energies are chosen to correspond to experimental bins in the Sudbury Neutrino Observatory (SNO) and Borexino experiments, which measure electron-flavor survival probability at Earth. The procedure successfully finds consistency between the observed fluxes and the model, if the MSW resonance - that is, flavor evolution due to solar electrons - is included in the dynamical equations representing the model.

We present an analysis of Epoch of Reionization data from Phase II of the Murchison Widefield Array using the \texttt{simpleDS} delay spectrum pipeline. This analysis is complementary to that presented in Li et al. (2019), which analyzes the same observations using the FHD/$\varepsilon$ppsilon imaging pipeline. This represents the first time that both principal types of 21 cm cosmology power spectrum estimation approaches have been applied to the same data set. Our limits on the 21 cm power spectrum amplitude span a range in $k$ space of $|k| < 1~h_{100}{\rm Mpc}^{-1}$ with a lowest measurement of $\Delta^2(k) \leq$ $4.58\times10^3$ mK$^2$ at $k = 0.190 h_{100}\rm{Mpc}^{-1}$ and $z = 7.14$. In order to achieve these limits, we need to mitigate a previously unidentified common mode systematic in the data set. If not accounted for, this systematic introduces an overall \emph{negative} bias that can make foreground contaminated measurements appear as stringent, noise-limited constraints on the 21 cm signal amplitude. The identification of this systematic highlights the risk in modeling systematics as positive-definite contributions to the power spectrum and in ``conservatively'' interpreting all measurements as upper limits.

We develop an optimal Bayesian solution for jointly inferring secondary signals in the Cosmic Microwave Background (CMB) originating from gravitational lensing and from patchy screening during the epoch of reionization. This method is able to extract full information content from the data, improving upon previously considered quadratic estimators for lensing and screening. We forecast constraints using the Marginal Unbiased Score Expansion (MUSE) method, and show that they are largely dominated by CMB polarization, and depend on the exact details of reionization. For models consistent with current data which produce the largest screening signals, a detection (3\,$\sigma$) of the cross-correlation between lensing and screening is possible with SPT-3G, and a detection of the auto-correlation is possible with CMB-S4. Models with the lowest screening signals evade the sensitivity of SPT-3G, but are still possible to detect with CMB-S4 via their lensing cross-correlation.

This paper provides analyses of the emission beam structure of 76 ``B''-named pulsars within the Arecibo sky. Most of these objects are included in both the Gould & Lyne and LOFAR High Band surveys and thus complement our other works treating various parts of these populations. These comprise a further group of mostly well studied pulsars within the Arecibo sky that we here treat similarly to those in Olszanski et.al. and extend our overall efforts to study all of the pulsars in both surveys. The analyses are based on observations made with the Arecibo Telescope at 327 MHz and 1.4 GHz. Many have been observed at frequencies down to 100 MHz using either LOFAR or the Pushchino Radio Astronomy Observatory as well as a few with the Long Wavelength Array at lower frequencies. This work uses the Arecibo observations as a foundation for interpreting the low frequency profiles and emission-beam geometries. We attempt to build quantitative geometric emission-beam models using the core/double-cone topology, while reviewing the evidence of previous studies and arguments for previous classifications on these sources. These efforts were successful for all but two pulsars, and interesting new subpulse modulation patterns were identified in a number of the objects. We interpret the Arecibo pulsar population in the context of the entire population of ``B'' pulsars.

The Transiting Exoplanet Survey Satellite mission identified a deep and asymmetric transit-like signal with a periodicity of 4.5 days orbiting the M4 dwarf star TOI-3884. The signal has been confirmed by follow-up observations collected by the ExTrA facility and Las Cumbres Observatory Global Telescope, which reveal that the transit is chromatic. The light curves are well modelled by a host star having a large polar spot transited by a 6-R$_{\oplus}$ planet. We validate the planet with seeing-limited photometry, high-resolution imaging, and radial velocities. TOI-3884 b, with a radius of $6.00 \pm 0.18$ R$_{\oplus}$, is the first sub-Saturn planet transiting a mid-M dwarf. Owing to the host star's brightness and small size, it has one of the largest transmission spectroscopy metrics for this planet size and becomes a top target for atmospheric characterisation with the James Webb Space Telescope and ground-based telescopes.

Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bursts come from two magnetars, SGR 1935$+$2154 and Swift J1818.0$-$1607. We also include three other electromagnetic burst events detected by Fermi GBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper bounds on the root-sum-square of the integrated gravitational-wave strain that reach $2.2 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at 100 Hz for the short-duration search and $8.7 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at $450$ Hz for the long-duration search, given a detection efficiency of 50%. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to $1.8 \times 10^{-22}$ $/\sqrt{\text{Hz}}$. Using the estimated distance to each magnetar, we derive upper bounds on the emitted gravitational-wave energy of $3.2 \times 10^{43}$ erg ($7.3 \times 10^{43}$ erg) for SGR 1935$+$2154 and $8.2 \times 10^{42}$ erg ($2.8 \times 10^{43}$ erg) for Swift J1818.0$-$1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935$+$2154 with available fluence information. The lowest of these ratios is $3 \times 10^3$.

Images obtained with single-conjugate adaptive optics (AO) show spatial variation of the point spread function (PSF) due to both atmospheric anisoplanatism and instrumental aberrations. The poor knowledge of the PSF across the field of view strongly impacts the ability to take full advantage of AO capabilities. The AIROPA project aims to model these PSF variations for the NIRC2 imager at the Keck Observatory. Here, we present the characterization of the instrumental phase aberrations over the entire NIRC2 field of view and we present a new metric for quantifying the quality of the calibration, the fraction of variance unexplained (FVU). We used phase diversity measurements obtained on an artificial light source to characterize the variation of the aberrations across the field of view and their evolution with time. We find that there is a daily variation of the wavefront error (RMS of the residuals is 94~nm) common to the whole detector, but the differential aberrations across the field of view are very stable (RMS of the residuals between different epochs is 59~nm). This means that instrumental calibrations need to be monitored often only at the center of the detector, and the much more time-consuming variations across the field of view can be characterized less frequently (most likely when hardware upgrades happen). Furthermore, we tested AIROPA's instrumental model through real data of the fiber images on the detector. We find that modeling the PSF variations across the field of view improves the FVU metric by 60\% and reduces the detection of fake sources by 70\%.

We study ejection mechanisms for two kinds of steady jets: one observed from black hole binaries in the low/hard state and the other from SS433. The specific energy of the ejected gas is required to be positive for the jets to get to infinity, while that of the accreted gas is naively considered to be negative at the outermost boundary of the accretion flow. To reconcile the opposite sign of the specific energies, we propose a situation where two layers exist in the accretion flow and one layer receives energy from the other sufficiently for the specific energy to be positive. For the steady jets in the low/hard state, the accretion ring at the outermost end of the accretion flow is considered to yield two-layer flow in which a geometrically thick ADAF sandwiches a geometrically thin accretion disk and the thin disk is supposed to turn to another ADAF on the inner side. The energy transfer is expected to occur through turbulent mixing between the two layers and the upper layer is discussed to have the positive specific energy large enough for the terminal velocity to be $\sim$ 0.1 $c$. For the steady jets from SS433, a slim disk is argued to separate into two stratified layers due to the photon diffusion in the direction perpendicular to the equatorial plane under the advection dominated situation. In this case, the specific energy of the upper layer is expected to be positive such that the terminal velocity exceeds 0.2 $c$. The jet ejection process near the black hole is investigated commonly to both the two-layer cases and predicts the jet opening angle becomes as small as 2$^{\circ}$.

Interaction of three-dimensional magnetic fields, turbulence, and self-gravity in the molecular cloud is crucial in understanding star formation but has not been addressed so far. In this work, we target the low-mass star-forming region L1688 and use the spectral emissions of $^{12}$CO, $^{13}$CO, C$^{18}$O, and H I, as well as polarized dust emissions. To obtain the 3D direction of the magnetic field, we employ the polarization fraction analysis developed in Hu & Lazarian (2022). In combining with the plane-of-the-sky (POS) magnetic field strength derived from the Davis-Chandrasekhar-Fermi (DCF) method and the new Differential Measure Analysis (DMA) technique introduced in Lazarian et al. (2022), we present the first measurement of L1688's three-dimensional magnetic field, including its orientation and strength. We find that L1688's magnetic field has two statistically different inclination angles. The low-intensity tail has an inclination angle $\approx55^\circ$ on average, while that of the central dense clump is $\approx30^\circ$. We find the global mean value of total magnetic field strength is $B_{\rm tot}\approx$ 135 uG from DCF and $B_{\rm tot}\approx$ 75 uG from DMA. We use the velocity gradient technique (VGT) to separate the magnetic fields' POS orientations associated with L1688 and its foreground/background. The magnetic fields' orientations are statistically coherent. The probability density function of H$_2$ column density and VGT reveal that L1688 potentially is undergoing gravitational contraction at large scale $\approx1.0$ pc and gravitational collapse at small scale $\approx0.2$ pc. The gravitational contraction mainly along the magnetic field resulting in an approximate power-law relation $B_{\rm tot}\propto n_{\rm H}^{1/2}$ when volume density $n_{\rm H}$ is less than approximately $6.0\times10^3$ cm$^{-3}$.

We undertook a search for new dwarf galaxies in the vicinity of relatively isolated nearby galaxies with distances $D < 12$ Mpc and stellar masses in the $2\times10^{11}-3\times10^8~M_{\odot}$ interval, using the data from the DESI Legacy Imaging Surveys. Around the 46 considered Local Volume galaxies, $67$ new candidates for satellites of these galaxies were found. About half of them are classified as spheroidal dwarfs of low surface brightness. The new galaxies are included in the Local Volume database (LVGDB), which now contains 1421 objects, being 63% more than the Updated Nearby Galaxy Catalog.

Fast radio bursts (FRBs) are intense bursts of radio emission with durations of milliseconds. Although researchers have found them happening frequently all over the sky, they are still in the dark to understand what causes the phenomena because the existing radio observatories have encountered certain challenges during the discovery of FRB progenitors. The construction of Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is being proposed to solve these challenges. We simulate mock Galactic FRB-like events by applying a range of spatial distributions, pulse widths and luminosity functions. The effect of turbulent Interstellar Medium (ISM) on the detectability of FRB-like events within the Milky Way plane is considered to estimate the dispersion measure and pulse scattering of mock events. We evaluate the fraction of FRB-like events in the Milky Way that are detectable by BURSTT and compare the result with those by Survey for Transient Astronomical Radio Emission 2 (STARE2) and Galactic Radio Explorer (GReX). We find that BURSTT could increase the detection rate by more than two orders of magnitude compared with STARE2 and GReX, depending on the slope of luminosity function of the events. We also investigate the influence of the specifications of BURSTT on its detection improvement. This leads to the fact that greatly higher sensitivity and improved coverage of the Milky Way plane have significant effects on the detection improvement of BURSTT. We find that the upgrade version of BURSTT, BURSTT-2048 could increase the detection rate of faint Galactic FRB-like events by a factor of 3.

In the previous paper [arXiv:2210.10435], the nonlinear perturbation theory of cosmological density field is generalized to include the tensor-valued bias of astronomical objects, such as spins and shapes of galaxies and any other tensors of arbitrary ranks which are associated with objects that we can observe. We apply this newly developed method to explicitly calculate nonlinear power spectra and correlation functions both in real space and in redshift space. Multi-dimensional integrals that appeared in loop corrections are reduced to combinations of the one-dimensional Hankel transforms, thanks to the spherical basis of the formalism, and the final expressions are numerically evaluated in a very short time using an algorithm of the fast Fourier transforms such as FFTLog. As an illustrative example, numerical evaluations of loop corrections of the power spectrum and correlation function of the rank-2 tensor field are demonstrated with a simple model of tensor bias.

The abundance of galaxy clusters is a sensitive probe to the amplitude of matter density fluctuations, the total amount of matter in the Universe as well as its expansion history. Inferring correct values and accurate uncertainties of cosmological parameters requires accurate knowledge of cluster abundance statistics, encoded in the likelihood function. In this paper, we test the accuracy of cluster abundance likelihoods used in the literature, namely the Poisson and Gaussian likelihoods as well as the more complete description of the Gauss-Poisson Compound likelihood. This is repeated for a variety of binning choices and analysis setups. In order to evaluate the accuracy of a given likelihood, this work compares individual posterior covariances to the covariance of estimators over the 1000 simulated dark matter halo catalogs obtained from PINOCCHIO algorithm. We find that for Rubin/LSST or Euclid-like surveys the Gaussian likelihood gives robust constraints over a large range of binning choices. The Poisson likelihood, that does not account for sample covariance, always underestimates the errors on the parameters, even when the sample volume is reduced or only high-mass clusters are considered. We find no benefit in using the more complex Gauss-Poisson Compound likelihood as it gives essentially the same results as the Gaussian likelihood, but at a greater computational cost. Finally, in this ideal setup, we note only a small gain on the parameter error bars when using a large number of bins in the mass-redshift plane.

The content and distribution of cool interstellar medium (ISM, <30K) can indicate the evolutionary mechanisms that transform late type to early type galaxies (ETGs). To investigate this, ALMA observations of 12CO[2-1] line emission were obtained for five dusty ETGs from a complete sample in low-density environments. Four of the ETGs have massive (approximately 10^9 Msolar), extended molecular gas reservoirs with effective radii of approximately 3 to 5 kpc. This work provides a kinematic and structural analysis of these observations, to explore possible evolutionary mechanisms. Axisymmetric or bisymmetric kinematic models were fitted to observations of molecular gas discs, to quantify the dominant structures present and highlight additional structures or asymmetries. Integral Field Unit (IFU) observations of these ETGs were also examined where available. Two of the ETGs, GAMA64646 and 622305, appear to have undergone tidal disturbance leading to molecular gas discs and/or star-forming inner rings. GAMA272990 may have undergone a merger, leading to an elliptical galaxy with an embedded star-forming molecular gas disc. GAMA622429 has probably undergone a minor merger, indicated by asymmetry in molecular gas distribution and disturbance in optical images. The remaining ETG, GAMA177186, was affected by source confusion from an offset source which could be a high mass, dust- and gas- rich object at high redshift. Overall, it appears that a high proportion of dusty ETGs in low-density environments have massive, extended molecular gas reservoirs, and have undergone some kind of interaction as part of their recent evolution. Secular evolution can then (re-)transform the ETGs from star-forming to passive galaxies.

Although high energetic radiation from flares is a potential threat to exoplanet atmospheres and may lead to surface sterilization, it might also provide the extra energy for low-mass stars needed to trigger and sustain prebiotic chemistry. We investigate two flares on TRAPPIST-1, an ultra-cool dwarf star that hosts seven exoplanets of which three lie within its habitable zone. The flares are detected in all four passbands of the MuSCAT2 allowing a determination of their temperatures and bolometric energies. We analyzed the light curves of the MuSCAT1 and MuSCAT2 instruments obtained between 2016 and 2021 in $g,r,i,z_\mathrm{s}$-filters. We conducted an automated flare search and visually confirmed possible flare events. We studied the temperature evolution, the global temperature, and the peak temperature of both flares. For the first time we infer effective black body temperatures of flares that occurred on TRAPPIST-1. The black body temperatures for the two TRAPPIST-1 flares derived from the SED are consistent with $T_\mathrm{SED} = 7940_{-390}^{+430}$K and $T_\mathrm{SED} = 6030_{-270}^{+300}$K. The flare black body temperatures at the peak are also calculated from the peak SED yielding $T_\mathrm{SEDp} = 13620_{-1220}^{1520}$K and $T_\mathrm{SEDp} = 8290_{-550}^{+660}$K. We show that for the ultra-cool M-dwarf TRAPPIST-1 the flare black body temperatures associated with the total continuum emission are lower and not consistent with the usually adopted assumption of 9000-10000 K. This could imply different and faster cooling mechanisms. Further multi-color observations are needed to investigate whether or not our observations are a general characteristic of ultra-cool M-dwarfs. This would have significant implications for the habitability of exoplanets around these stars because the UV surface flux is likely to be overestimated by the models with higher flare temperatures.

The analysis of the colour of artificial lights at night has an impact on diverse fields, but current data sources have either limited resolution or scarce availability of images for a specific region. In this work, we propose crowdsourced photos of streetlights as an alternative data source: for this, we designed NightUp Castelldefels, a pilot for a citizen science experiment aimed at collecting data about the colour of streetlights. In particular, we extract the colour from the collected images and compare it to an official database, showing that it is possible to classify streetlights according to their colour from photos taken by untrained citizens with their own smartphones. We also compare our findings to the results obtained from one of the current sources for this kind of study. The comparison highlights how the two approaches give complementary information about artificial lights at night in the area. This work opens a new avenue in the study of the colour of artificial lights at night with the possibility of accurate, massive and cheap data collection.

Polycyclic aromatic hydrocarbons (PAHs) and fullerenes play a major role in the physics and chemistry of the interstellar medium. Based on a number of recent experimental and theoretical investigations we developed a model in which PAHs are subject to photo-dissociation (carbon and hydrogen loss) and hydrogenation. We take into account that dehydrogenated PAHs may fold into closed structures -- fullerenes. Fullerenes, in their turn, can be also hydrogenated, becoming fulleranes, and photo-dissociated, losing carbon and hydrogen atoms. The carbon loss leads to shrinking of fullerene cages to smaller ones. We calculate the abundance of PAHs and fullerenes of different sizes and hydrogenation level depending on external conditions: the gas temperature, intensity of radiation field, number density of hydrogen atoms, carbon atoms, and electrons. We highlight the conditions, which are favourable for fullerene formation from PAHs, and we conclude that this mechanism works not only in H-poor environment but also at modest values of hydrogen density up to 10$^{4}$~cm$^{-3}$. We found that fulleranes can be formed in the ISM, although the fraction of carbon atoms locked in them can be maximum around 10$^{-9}$. We applied our model to two photo-dissociation regions, Orion Bar and NGC 7023. We compare our estimates of the fullerene abundance and synthetic band intensities in these objects with the observations and conclude that our model gives good results for the closest surroundings of ionising stars. We also demonstrate that additional fullerene formation channels should operate along with UV-induced formation to explain abundance of fullerenes far from UV sources.

Both the dynamics and the observational properties of relativistic jets are determined by their interaction with the ambient medium. A crucial role is played by the contact discontinuity at the jet boundary, which in the presence of jet collimation may become subject to Rayleigh-Taylor instability (RTI) and Richtmyer-Meshkov instability (RMI). Here, we study the evolution of these instabilities in non-magnetized relativistic jets using special relativistic three-dimensional hydrodynamic simulations. We show that the growth of initial perturbations is consistent with relativistic RTI operating in the jet collimation region. The contribution of RMI becomes important downstream from the collimation shock in agreement with theoretical expectations. Both instabilities reach non-linear scales above the shock convergence point and trigger strong turbulence, mixing jet with ambient matter. We devise an analytic solution for the mixing rate and show that it is sensitive to the external density gradients. Our results may be applied to different types of astrophysical objects. In particular, different contribution of interface instabilities is a natural explanation for the observed dichotomy between FR-I and FR-II radiogalaxies. The rapid slow-down in the jet of M87 is consistent with baryon entrainment from the circumnuclear matter with the observed density distribution. In microquasars, baryon loading triggered by interface instabilities is a probable reason for the low observed Lorentz factors. We show that the observed variability in gamma-ray bursts cannot come from mixing driven by interface instabilities and likely originates from the engine, suggesting the presence of magnetic fields in the jet.

Stars are born with magnetic fields, but the distribution of their initial field strengths remains uncertain. We combine observations with theoretical models of magnetic field evolution to infer the initial distribution of magnetic fields for AB stars in the mass range of 1.6 - 3.4 M$_{\odot}$. We tested a variety of distributions with different shapes and found that a distribution with a mean of $\sim$800 G and a full width of $\sim$600 G is most consistent with the observed fraction of strongly magnetized stars as a function of mass. Our most-favored distribution is a Gaussian with a mean of $\mu$ = 770 G and standard deviation of $\sigma$ = 146 G. Independent approaches to measure the typical field strength suggest values closer to 2 - 3 kG, a discrepancy that could suggest a mass-dependent and bimodal initial field distribution, or an alternative theoretical picture for the origin of these magnetic fields.

Ground-based observing time is precious in the era of exoplanet follow-up and characterization, especially in high-precision radial velocity instruments. Blind-search radial velocity surveys thus require a dedicated observational strategy in order to optimize the observing time, which is particularly crucial for the detection of small rocky worlds at large orbital periods. We develop an algorithm with the purpose of improving the efficiency of radial velocity observations in the context of exoplanet searches, and we apply it to the K-dwarfs Orbited By habitable Exoplanets (KOBE) experiment. We aim at accelerating exoplanet confirmations or, alternatively, rejecting false signals as early as possible in order to save telescope time and increase the efficiency of both blind-search surveys and follow-up of transiting candidates. Once a minimum initial number of radial velocity datapoints is reached in such a way that a periodicity starts to emerge according to generalized Lomb-Scargle (GLS) periodograms, that period is targeted with the proposed algorithm, named $\texttt{KOBEsim}$. The algorithm selects the next observing date that maximizes the Bayesian evidence for such periodicity in comparison with a model with no Keplerian orbits. By means of simulated data, we prove that the algorithm accelerates the exoplanet detection, needing $29 - 33\,\%$ less observations and $41 - 47\,\%$ less timespan of the full dataset for low-mass planets ($m_{\rm p}\,<\,10\,M_{\oplus}$) in comparison with a conventional monotonic cadence strategy. The enhancement in the number of datapoints for $20\,M_{\oplus}$ planets is also appreciable, $16\,\%$. We also test $\texttt{KOBEsim}$ with real data for a particular KOBE target, and for the confirmed planet $HD~102365\,b$. Both of them demonstrate that the strategy is capable of speeding up the detection up to a factor of $2$.

Forecast of optical turbulence and atmospheric parameters relevant for ground-based astronomy is becoming an important goal for telescope planning and AO instruments optimization in several major telescope. Such detailed and accurate forecast is typically performed with numerical atmospheric models. Recently short-term forecasts (a few hours in advance) are also being provided (ALTA project) using a technique based on an autoregression approach, as part of a strategy that aims to increase the forecast accuracy. It has been proved that such a technique is able to achieve unprecedented performances so far. Such short-term predictions make use of the numerical model forecast and real-time observations. In recent years machine learning (ML) techniques also started to be used to provide an atmospheric and turbulence forecast. Preliminary results indicate however an accuracy not really competitive with respect to the autoregressive method or even prediction by persistence. This technique might be applicable joint to atmospheric model. It is therefore interesting to investigate the main features of their performances and characteristics (also because there is a great number of algorithms potentially accessible) to understand if results achieved so far can be further improved using ML. In this study we focus on a purely machine learning application to short term forecast (1-2 hours) of astroclimatic and other atmospheric parameters above VLT.

Characterizing the PSF of adaptive optics instruments is of paramount importance both for instrument design and observation planning/optimization. Simulation software, such as PASSATA, have been successfully utilized for PSF characterization in instrument design, which make use of standardized atmospheric turbulence profiles to produce PSFs that represent the typical instrument performance. In this contribution we study the feasibility of using such tool for nowcast application (present-time forecast), such as the characterization of an on-sky measured PSF in real observations. Specifically we will analyze the performance of the simulation software in characterizing the real-time PSF of two different state-of-the-art SCAO adaptive optics instruments: SOUL at the LBT, and SAXO at the VLT. The study will make use of on-sky measurements of the atmospheric turbulence and compare the results of the simulations to the measured PSF figures of merit (namely the FHWM and the Strehl Ratio) retrieved from the instrument telemetry in real observations. Our main goal in this phase is to quantify the level of uncertainly of the AO simulations in reproducing real on-sky observed PSFs with an end-to-end code (PASSATA). In a successive phase we intend to use a faster analytical code (TIPTOP). This work is part of a wider study which aims to use simulation tools joint to atmospheric turbulence forecasts performed nightly to forecast in advance the PSF and support science operations of ground-based telescopes facilities. The 'PSF forecast' option might therefore be added to ALTA Center or the operational forecast system that will be implemented soon at ESO.

ALTA project has been active since 2016, providing, at LBT observatory site, forecasts of atmospheric parameters, such as temperature, wind speed and direction, relative humidity and precipitable water vapor, and optical turbulence parameters, such as seeing, wavefront coherence time and isoplanatic angle with the final goal to support nightly the science operation of the LBT. Besides to the forecasts, during the years ALTA has been collecting statistics on the atmospheric conditions which can be used to draw a very accurate characterization of the climatology of the telescope site located on top of Mount Graham, Arizona. Such characterization can be used both for the optimization and calibration of the forecast model and as a reference for a model validation. The climatology of these parameters is conceived to be a further output of ALTA that will be upgraded on the website with time and it will be able to put in evidence trends at short as well as long time scales. In this contribution we present a climatological description of all the atmospheric parameters relevant for ground-based astronomy in order to provide to the scientific community a robust reference of the bserving conditions at LBT. The study is performed using on-site measurements provided by DIMM and atmospheric sensors over several years and made available in the telescope telemetry data.

The preliminary detections of the gamma-ray burst 221009A up to 18 TeV by LHAASO and up to 251 TeV by Carpet 2 have been reported through Astronomer's Telegrams and Gamma-ray Coordination Network circulars. Since this burst is at redshift $z=0.15$, these photons may at first seem to have a low probability to avoid pair production off of background radiation fields and survive to reach detectors on Earth. By extrapolating the reported $0.1-1.0$\ GeV LAT spectrum from this burst to higher energies and using this to limit the intrinsic spectrum of the burst, we show that the survival of the 18 TeV photon detected by LHAASO is not unlikely with many recent extragalactic background light models, although the detection of a 251 TeV event is still very unlikely. This can be resolved if Lorentz invariance is violated at an energy scale $E_{\rm QG}> 49 E_{\rm Planck}$ in the linear ($n=1$) case, and $E_{\rm QG}> 10^{-6}E_{\rm Planck}$ in the quadratic ($n=2$) case (95% confidence limits), where $E_{\rm Planck}$ is the Planck energy. This could potentially be the first evidence for subluminal Lorentz invariance violation.

We describe a new spectrophotometer for the Birmingham Solar Oscillations Network (BiSON), based on a next generation observation platform, BiSON:NG, a significantly miniaturised system making use of inexpensive consumer-grade hardware and off-the-shelf components, where possible. We show through system modelling and simulation, along with a summer observing campaign, that the prototype instrument produces data on the Sun's low-degree acoustic (p-mode) oscillations that are of equal quality and can be seamlessly integrated into the existing network. Refreshing the existing ageing hardware, and the extended observational network potential of BiSON:NG, will secure our ongoing programme of high-quality synoptic observations of the Sun's low-degree oscillations (e.g., for seismic monitoring of the solar cycle at a "whole Sun" level).

The presence of possible correlations between stellar rotation rate $\Omega$ and the frequency of the activity cycle $\omega_\mathrm{cyc}$ is still much debated. We implement a new Bayesian classification algorithm based on a simultaneous regression analysis of multiple scaling laws and we demonstrate the existence of two different scalings in the $\log_{10} \omega_\mathrm{cyc}$ -- $\log_{10} \Omega$ plane for an extended Mt.~Wilson sample of 67 stars. Thanks to metallicity measurements obtained from both ESA Gaia and high-resolution spectroscopy, we argue that the origin of this dichotomy is likely related to the chemical composition: stars whose magnetic cycle frequency increases with rotation rate are less metallic than stars whose magnetic cycle frequency decreases with stellar rotation rates. On the contrary, no clear difference in chromospheric magnetic activity indicators characterizes the two branches.

The introduction of the Rossby number (R$_0$), which incorporates the convective turnover time ($\tau$), in 1984 was a pioneering idea for understanding the correlation between stellar rotation and activity. The convective turnover time, which cannot be measured directly, is often inferred using existing $\tau$-mass or $\tau$-color relations, typically established based on an ensemble of different types of stars by assuming that $\tau$ is a function of mass. In this work, we use {\it Gaia} Early Data Release 3 to demonstrate that the masses used to establish one of the most cited $\tau$-mass relations are overestimated for G type dwarfs and significantly underestimated for late M dwarfs, offsets that affect studies using this $\tau$-mass relation to draw conclusions. We discuss the challenges of creating such relations then and now. In the era of {\it Gaia} and other large datasets, stars used to establish these relations require characterization in a multi-dimensional space, rather than via the single-characteristic relations of the past. We propose that new multi-dimensional relations should be established based on updated theoretical models and all available stellar parameters for different interior structures from a set of carefully vetted single stars, so that the convective turnover time can be estimated more accurately.

We design a new observable, the expansion rate fluctuation $\eta$, to characterize deviations from the linear relation between redshift and distance in the local universe. We also show how to compress the resulting signal into spherical harmonic coefficients in order to better decipher the structure and symmetries of the anisotropies in the local expansion rate. We apply this analysis scheme to several public catalogs of redshift-independent distances, the Cosmicflows-3 and Pantheon data sets, covering the redshift range $0.01<z<0.05$. The leading anisotropic signal is stored in the dipole. Within the standard cosmological model, it is interpreted as a bulk motion ($307 \pm 23$ km/s) of the entire local volume in a direction aligned at better than $4$ degrees with the bulk component of the Local Group velocity with respect to the CMB. This term alone, however, provides an overly simplistic and inaccurate description of the angular anisotropies of the expansion rate. We find that the quadrupole contribution is non-negligible ($\sim 50\%$ of the anisotropic signal), in fact, statistically significant, and signaling a substantial shearing of gravity in the volume covered by the data. In addition, the 3D structure of the quadrupole is axisymmetric, with the expansion axis aligned along the axis of the dipole. Implications for the determination of the $H_0$ parameter are also discussed.

The recent detections of binary stellar mass black hole mergers by the LIGO and Virgo Collaborations suggest that such mergers are common occurrences. Galaxy mergers further indicate that supermassive black holes in centers of galaxies also merge and are typically expected to have had at least one merger in their lifetime, possibly many. In the presence of a jet, these mergers are almost always accompanied by a change of the jet direction and a connected jet precession motion, leading to interactions of the jet with ambient matter and producing very high-energy particles, and consequently high-energy gamma-rays and neutrinos. In this work, we investigate the possibility under which conditions such mergers could be the sources of the diffuse astrophysical neutrino flux measured by the IceCube Neutrino Observatory. The main free parameters in the calculation concern the frequency of the mergers and the fraction of energy that is transferred from the gravitationally released energy to neutrinos. We show that the merger rate for SMBBHs must lie between $\sim 10^{-7}$ and $10^{-5}$ Gpc$^{-3}$ yr$^{-1}$. The ratio of energy going to neutrinos during such mergers lies then between $\sim 10^{-6} - 3\cdot 10^{-4}$. For stellar mass BBH mergers, the rate needs to be $\sim 10-100$ Gpc$^{-3}$ yr$^{-1}$ and the expected ratio of neutrino to gravitational wave energy lies in a comparable range as for SMBBHs, $\sim 2 \cdot 10^{-5} - 10^{-3}$. These values lie in a reasonable parameter range, so that the production of neutrinos at the level of the detected neutrino flux is a realistic possibility.

We report the discovery of a dark companion to 2MASS J15274848+3536572 with an orbital period of 6.14 hours. Combining the radial velocity of LAMOST observation and the modelling of the multi-band light curve, one has a mass function of 0.135 Msun, an inclination of 43.94+0.33-0.21, and a mass ratio of 0.58+0.048-0.018, which demonstrate the binary nature of a dark companion with mass of 1.01+-0.08Msun and a main-sequence K star of 0.59+-0.05 Msun. LAMOST optical spectra at a range of orbital phase reveals extra peaked Halpha emission that suggests the presence of an accretion disk. The dark companion does not seem to be a white dwarf because the lack of any observed dwarf nova outbursts contradicts with the disk instability model in long-term data archive. Alternatively, we propose a scenario that the dark companion is a neutron star, but we have not detected radio pulsation or single pulse from the system with the FAST (Five hundred meter Aperture Spherical radio Telescope), which hints a radio quiet compact object. If the dark companion is identified as a neutron star, it will be nearest (118 pc) and lightest neutron star. Furthermore, a kinematic analysis of the systems orbit in the galaxy may suggest its supernova event is associated with the radionuclide 60Fe signal observed from the deep-sea crusts. This radio-quiet and X-ray dim nearby neutron star may resemble an XDINS (X-ray dim isolated neutron star), but in a binary.

It is postulated in Einstein's relativity that the speed of light in vacuum is a constant for all observers. However, the effect of quantum gravity could bring an energy dependence of light speed, and a series of previous researches on high-energy photon events from gamma-ray bursts (GRBs) and active galactic nuclei (AGNs) suggest a light speed variation $v(E)=c\left(1-E / E_{\mathrm{LV}}\right)$ with $E_{\mathrm{LV}}=3.6 \times 10^{17} ~\mathrm{GeV}$. From the newly detected gamma-ray burst GRB 221009A, we find that a $99.3~$GeV photon detected by Fermi-LAT is coincident with the sharp spike in the light curves detected by Fermi-GBM and HEBS under the above scenario of light speed variation, suggesting that this high energy photon was emitted at the same time with a sharp spike of low energy photon emission at the GRB source. Thus this highest energy photon event detected by Fermi-LAT during the prompt emission of gamma ray bursts supports the linear form modification of light speed in cosmological space.

We report on the discovery of two new molecules, HCCCHCCC and HCCCCS, towards the starless core TMC-1 in the Taurus region from the QUIJOTE line survey in the 31.1-50.2 GHz frequency range. We identify a total of twenty-nine lines of HCCCHCCC and six rotational transitions of HCCCCS. The rotational quantum numbers range from Ju=10 up to 15 and Ka <=2 for HCCCHCCC and Ju=21/2 up to 31/2 for HCCCCS. We derived a column density for HCCCHCCC of N=(1.3+/-0.2)x10^11 cm-2 with a rotational temperature of 6+/-1 K, while for HCCCCS we derived N=(9.5+/-0.8)x10^10 cm-2 and Trot =10+/-1 K. The abundance of HCCCHCCC is higher than that of its recently discovered isomer, l-H2C6. If we compare HCCCCS with its related molecules, HCS and HCCS, we obtain abundance ratios HCS/HCCCCS=58 and HCCS/HCCCCS=7.2. We investigated the formation of these two molecules using chemical modelling calculations. The observed abundances can be accounted for by assuming standard gas-phase formation routes involving neutral-neutral reactions and ion-neutral reactions.

The ages of the oldest and most metal-poor stars in the Milky Way bear important information on the age of the Universe and its standard model. We analyze a sample of 28 extremely metal-poor field stars in the solar vicinity culled from the literature and carefully determine their ages. To this aim, we critically make use of Gaia data to derive their distances and associated uncertainties. Particular attention has been paid to the estimate of the reddening and its effect on the derivation of stellar ages. We employed different reddenings and super-impose isochrones from different sources in the stars color-magnitude diagram built up with different photometric systems. We highlight subtle metallicity effects when using the Johnson photometry for low metallicity stars and finally adopt Gaia photometry. An automatic fitting method is devised to assign ages to each individual star taking into account the uncertainties in the input parameters. The mean age of the sample turns out to be $13.9 \pm 0.5$ Gyr using Padova isochrones, and $13.7 \pm 0.4$ Gyr using BASTI isochrones. We found also a group of very metal-poor stars ($\left[\frac{Fe}{H}\right]$: -2.7 -2.0 dex) with relatively young ages, in the range 8 --10 Gyr.

Most secular effects produced by stellar bars strongly depend on the pattern speed. Unfortunately, it is also the most difficult observational parameter to estimate. In this work, we measured the bar pattern speed of 97 Milky-Way Analogue galaxies from the MaNGA survey using the Tremaine-Weinberg method. The sample was selected by constraining the stellar mass and morphological type. We improve our measurements by weighting three independent estimates of the disc position angle. To recover the disc rotation curve, we fit a kinematic model to the H$_\alpha$ velocity maps correcting for the non-circular motions produced by the bar. The complete sample has a smooth distribution of the bar pattern speed ($\Omega_{Bar}=28.14^{+12.30}_{-9.55}$ km s$^{-1}$ kpc $^{-1}$), corotation radius ($R_{CR} = 7.82^{+3.99}_{-2.96}$ kpc) and the rotation rate ($\mathcal{R} = 1.35^{+0.60}_{-0.40}$). We found two sets of correlations: (i) between the bar pattern speed, the bar length and the logarithmic stellar mass (ii) between the bar pattern speed, the disc circular velocity and the bar rotation rate. If we constrain our sample by inclination within $30 \degree < i < 60 \degree$ and relative orientation $20\degree<|PA_{disc}-PA_{bar} |<70\degree$, the correlations become stronger and the fraction of ultra-fast bars is reduced from 20\% to 10\% of the sample. This suggest that a significant fraction of ultra-fast bars in our sample could be associated to the geometric limitations of the TW-method. By further constraining the bar size and disc circular velocity, we obtain a sub-sample of 25 Milky-Way analogues galaxies with distributions $\Omega_{Bar}=30.48^{+10.94}_{-6.57}$ km s$^{-1}$ kpc$^{-1}$, $R_{CR} = 6.77^{+2.32}_{-1.91}$ kpc and $\mathcal{R} = 1.45^{+0.57}_{-0.43}$, in good agreement with the current estimations for our Galaxy.

Anomalous large-angle cosmic microwave background anisotropies motivate further searches for cosmic topology. We demonstrate that for generic topologies of spatially flat spacetimes, off-diagonal correlations between microwave background harmonic coefficients over a wide range of scales encode significant topological information, even if the topology scale substantially exceeds the diameter of the observable Universe. Observational searches have so far considered only a small subset of testable topologies, and current limits on the topology scale are much weaker than generally understood.

Understanding the formation and evolution of ring galaxies, galaxies with an atypical ring-like structure, will improve understanding of black holes and galaxy dynamics as a whole. Current catalogs of rings are extremely limited: manual analysis takes months to accumulate an appreciable sample of rings and existing computational methods are vastly limited in terms of accuracy and detection rate. Without a sizable sample of rings, further research into their properties is severely restricted. This project investigates the usage of a convolutional neural network (CNN) to identify rings from unclassified samples of galaxies. A CNN was trained on a sample of 100,000 simulated galaxies, transfer learned to a sample of real galaxies and applied to a previously unclassified dataset to generate a catalog of rings which was then manually verified. Data augmentation with a generative adversarial network (GAN) to simulate images of galaxies was also used. A catalog of 1151 rings was extracted with 7.4 times the precision and 15.4 times the detection rate of conventional algorithms. The properties of these galaxies were then estimated from their photometry and compared to the Galaxy Zoo 2 catalog of rings. With upcoming surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time obtaining images of billions of galaxies, similar models could be crucial in classifying large populations of rings to better understand the peculiar mechanisms by which they form and evolve.

The origin of ultra-high-energy cosmic rays (UHECRs, E $> 10^{18}$ eV) is one of the great mysteries of modern astrophysics. It has been suggested that UHECRs could be accelerated in gamma-ray bursts (GRBs) and engine-driven supernovae (SNe). Here we report the discovery of a 1.4 teraelectronvolt (TeV) photon offset 0.97$^{\circ}$ from the site of the nearby (36.9 megaparsecs) GRB 980425/SN 1998bw explosion. The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope detected the TeV emission on 17 November 2018, more than 20 years after the original GRB 980425/SN 1998bw trigger. TeV detections at high Galactic latitudes by the LAT are extremely rare, with an average of 6 events per year. We propose that the delayed TeV emission is consistent with ultra-high-energy cosmic rays and/or electron-positron pairs from GRB 980425/SN 1998bw being deflected by the intergalactic magnetic field (IGMF) and subsequently cascading into secondary gamma rays. Based on the arrival time delay of the TeV emission, we estimate an IGMF strength of order $B \simeq 10^{-12}$-$10^{-13}$ Gauss. This result supports the possibility of UHECR acceleration in GRB 980425/SN 1998bw and suggests that most detected UHECRs are produced in local GRB/SNe within 200 Mpc. In addition, secondary photons from UHECRs out to 0.9-31 Gpc may also offer an explanation for extragalactic background photons with energies $\geq 1$ TeV detected by the Fermi LAT.

We have searched for thermal gyro-synchrotron radio emission from a sample of five radio-loud stars whose X-ray coronae contain a hot ($T>10^7$ K) thermal component. We used the JVLA to measure Stokes I and V/I spectral energy distributions (SEDs) over the frequency range 15--45 GHz, determining the best-fitting model parameters using power-law and thermal gyro-synchrotron emission models. The SEDs of the three chromospherically active binaries (Algol, UX Arietis, HR 1099) were well-fit by a power-law gyro-synchrotron model, with no evidence for a thermal component. However, the SEDs of the two pre-main-sequence (PMS) stars (V410 Tau, HD 283572) had a circularly polarized enhancement above 30 GHz that was inconsistent with a pure power-law distribution. These spectra were well-fit by summing the emission from an extended coronal volume of power-law gyro-synchrotron emission and a smaller region with thermal plasma and a much stronger magnetic field emitting thermal gyro-synchrotron radiation. We used Bayesian inference to estimate the physical plasma parameters of the emission regions (characteristic size, electron density, temperature, power-law index, and magnetic field strength and direction) using the independently VLBI-measured radio sizes, X-ray luminosity, and magnetic field strength as priors, where available. The derived parameters were well-constrained but highly degenerate. The best-fitting temperatures for both PMS stars were $\sim0.5$ dex higher than the X-ray-derived temperatures. We argue that the power-law and thermal volumes in the PMS stars are probably not co-spatial and speculate that they may arise from two distinct regions, the stellar corona and the inner edge of their accretion disc, respectively.

We show that gravitational waves can act as waveguides for electromagnetic radiation, that is if the latter is initially aligned with the gravitational waves, then the alignment will survive during the propagation. The analysis is performed using the Hamiltonian formalism and the Jacobi equation for null geodesics and conditions for certain cases of polarization of the waves are obtained. The effect of waveguiding by the gravitational waves can influence the interpretation of associated gravitational and electromagnetic wave events, since the latter cannot necessarily obey the inverse square decay law for intensity.

Positivity bounds represent nontrivial limitations on effective field theories (EFTs) if those EFTs are to be completed into a Lorentz-invariant, causal, local, and unitary framework. While such positivity bounds have been applied in a wide array of physical contexts to obtain useful constraints, their application to inflationary EFTs is subtle since Lorentz invariance is spontaneously broken during cosmic inflation. One path forward is to employ a $\textit{Breit parameterization}$ to ensure a crossing-symmetric and analytic S-matrix in theories with broken boosts. We extend this approach to a theory with multiple fields, and uncover a fundamental obstruction that arises unless all fields obey a dispersion relation that is approximately lightlike. We then apply the formalism to various classes of inflationary EFTs, with and without isocurvature perturbations, and employ this parameterization to derive new positivity bounds on such EFTs. For multifield inflation, we also consider bounds originating from the generalized optical theorem and demonstrate how these can give rise to stronger constraints on EFTs compared to constraints from traditional elastic positivity bounds alone. We compute various shapes of non-Gaussianity (NG), involving both adiabatic and isocurvature perturbations, and show how the observational parameter space controlling the strength of NG can be constrained by our bounds.

The right-handed (RH) Higgs-induced neutrino mixing (RHINO) model explains neutrino masses and origin of matter in the universe within a unified picture. The mixing, effectively described by a dimension five operator, is responsible both for the production of dark neutrinos, converting a small fraction of seesaw neutrinos acting as source, and for their decays. We show that including the production of source neutrinos from Higgs portal interactions, their abundance can thermalise prior to the onset of source-dark neutrino oscillations, resulting into an enhanced production of dark neutrinos that thus can play the role of decaying dark matter (DM) for a much higher seesaw scale. This can be above the sphaleron freeze-out temperature and as high as $\sim 100\,{\rm TeV}$, so that strong thermal resonant leptogenesis for the generation of the matter-antimatter asymmetry is viable. We obtain a $\sim 1\,{\rm TeV}$--$1\,{\rm PeV}$ allowed dark neutrino mass range. Intriguingly, their decays can also explain a neutrino flux excess at ${\cal O}(100\,{\rm TeV})$ energies recently confirmed by the IceCube collaboration analysing 7.5yr HESE data. Our results also point to an effective scale for Higgs portal interactions nicely identifiable with the grandunified scale and many orders of magnitude below the effective scale for the mixing. We explain this hierarchy in a UV-complete model with a very heavy fermion as mediator: the first scale corresponds to the fundamental scale of new physics, while the second is much higher because of a very small coupling identifiable with a symmetry breaking parameter. Therefore, RHINO realises a simple unified model of neutrino masses and origin of matter in the universe currently under scrutiny at neutrino telescopes and potentially embeddable within a grandunified model.

Nucleon-nucleon short-range correlations (SRCs) induce a high momentum tail (HMT) in the single-nucleon momentum distribution function $n_{\v{k}}^J(\rho,\delta)$ in cold neutron-rich matter. While there are clear experimental evidences that the SRC/HMT effects are different for neutrons and protons and their strengths depend strongly on the isospin asymmetry of finite nuclei mostly based on electron-nucleus scattering experiments, much less is known experimentally about the SRC/HMT effects in the dense neutron-rich matter. To facilitate further explorations of SRC/HMT effects in dense neutron-rich matter especially with heavy-ion reactions involving high-energy radioactive beams as well as multimessenger observations of neutron stars and their mergers, by incorporating the SRC-induced HMT in $n_{\v{k}}^J(\rho,\delta)$ into a Gogny-like energy density functional we study SRC/HMT effects on the equation of state (EOS) especially its symmetry energy term and single-nucleon potential in the dense asymmetric nucleonic matter (ANM). Using a parametrization as a surrogate for the momentum-dependent kernel in the Gogny-like energy density functional (EDF) we derive analytical expressions for all components of the ANM EOS and their characteristics (e.g., magnitude, slope and curvature as well as nucleon effective mass) at saturation density $\rho_0$ as well as the momentum-dependent single-nucleon optical potential in neutron-rich matter using parameters characterizing nuclear interactions as well as the size, shape and isospin dependence of the HMT at $\rho_0$. Some consequences of the SRC/HMT effects on properties of neutron stars are also studied.

We derive the weak field limit of scalar-Gauss-Bonnet theory and place novel bounds on the parameter space using terrestrial and space-based experiments. In order to analyze the theory in the context of a wide range of experiments, we compute the deviations from Einstein gravity around source masses with planar, cylindrical, and spherical symmetry. We find a correction to the Newtonian potential around spherical and cylindrical sources that can be larger than PPN corrections sufficiently close to the source. We use this to improve on laboratory constraints on the scalar-Gauss-Bonnet coupling parameter $\Lambda$ by two orders of magnitude. Present laboratory and Solar System bounds reported here are superseded by tests deriving from black holes.

A white light illuminated etalon is a valuable resource for spectrograph calibration in radial velocity exoplanet detection, and other astronomical applications. These etalons benefit from low drift (${\sim} 10^{-11}$/day) and well-characterized stability of their mode structure. However, measurements of several etalon systems across bandwidths greater than $500$ nm indicate that the modes exhibit complex, wavelength-dependent drift. Surprisingly, modes in different regions of the spectrum were found to drift in different directions, implying that the optical length of the etalon is getting both longer and shorter at the same time. In this paper, we provide a solution to this puzzling observation. With Fresnel analysis and the transfer matrix method, we model the reflective phase of the multi-layer dielectric mirrors in the etalon used as a calibrator for the Habitable Zone Planet Finder (HPF). We use this phase to calculate the etalon mode positions and are able to reproduce the observed oscillatory chromatic drift of the etalon's mode spectrum across 800-1300 nm. Despite the complexity of the mirror structure, our modeling indicates that the gradual relaxation of the outermost layers of the etalon mirrors is the dominant source of the observed behavior. We also model the effect of temperature, incident angle alignment variations, and manufacturing tolerances, and show that they are likely not causes of the chromatic mode frequency shifts. Our work highlights techniques that can be employed in the design of broad bandwidth mirrors for future etalons to make them more useful for the highest precision astronomical spectroscopy.

We theoretically and numerically investigate the ion acoustic feature of collective Thomson scattering (CTS) in two-stream plasmas. When the electron distribution functions of two components overlap each other, the theoretical spectrum shows asymmetry that is not seen in the spectra from Maxwell distributions because of the electron Landau damping. Numerical simulations support the theoretical spectra. We also present the effect of two-stream type instability in the ion acoustic feature showing the opposite trend to the overlapped case resulting from the asymmetry in electron distribution function caused by the instability. Our results show that the two-stream plasmas have significant effects on CTS spectra and the waves resulting from instability can be observed in the ion acoustic feature.

Future gravitational wave observations are potentially sensitive to new physics corrections to the Higgs potential once the first-order electroweak phase transition arises. We study the SMEFT dimension-six operator effects on the Higgs potential, where three types of effects are taken into account: (i) SMEFT tree level effect on $\varphi^6$ operator, (ii) SMEFT tree level effect on the wave function renormalization of the Higgs field, and (iii) SMEFT top-quark one-loop level effect. The sensitivity of future gravitational wave observations to these effects is numerically calculated by performing a Fisher matrix analysis. We find that the future gravitational wave observations can be sensitive to (ii) and (iii) once the first-order electroweak phase transition arises from (i). The sensitivities of the future gravitational wave observations are also compared with those of future collider experiments.

In this work we shall study a new class of attractor models which we shall call generalized $R^p$-attractor models. This class of models is based on a generalization of the Einstein frame potential of $R^p$ $f(R)$ gravity models in the Jordan frame. We present the attractor properties of the corresponding non-minimally coupled Jordan frame theory, and we calculate the observational indices of inflation in the Einstein frame. As we show, there is a large class of non-minimally coupled scalar theories, with an arbitrary non-minimal coupling which satisfies certain conditions, that yield the same Einstein frame potential, this is why these models are characterized attractors. As we demonstrate, the generalized $R^p$-attractor models are viable and well fitted within the Planck constraints. This includes the subclass of the generalized $R^p$-attractor models, namely the Einstein frame potential of $R^p$ inflation in the Jordan frame, a feature also known in the literature. We also highlight an important issue related to the $R^p$ inflation in the Jordan frame, which is known to be non-viable. By conformal invariance, the $R^p$ inflation model should also be viable in the Jordan frame, which is not. We pinpoint the source of the problem using two different approaches in the $f(R)$ gravity Jordan frame, and as we demonstrate, the problem arises in the literature due to some standard simplifications made for the sake of analyticity. We demonstrate the correct way to analyze $R^p$ inflation in the Jordan frame, using solely the slow-roll conditions.

Several different factors influence the seasonal cycle of a planet. This study uses a general circulation model and an energy balance model (EBM) to assess the parameters that govern the seasonal cycle. We define two metrics to describe the seasonal cycle, $\phi_s$, the latitudinal shift of the maximum temperature, and $\Delta T$, the maximum seasonal temperature variation amplitude. We show that alongside the expected dependence on the obliquity and orbital period, where seasonality generally strengthens with an increase in these parameters, the seasonality depends in a nontrivial way on the rotation rate. While the seasonal amplitude decreases as the rotation rate slows down, the latitudinal shift, $\phi_s$, shifts poleward. A similar result occurs in a diffusive EBM with increasing diffusivity. These results suggest that the influence of the rotation rate on the seasonal cycle stems from the effect of the rotation rate on the atmospheric heat transport.

We investigate collective Thomson scattering (CTS) in two-stream non-equilibrium plasmas analytically, numerically and experimentally. In laboratory astrophysics, CTS is a unique tool to obtain local plasma diagnostics. While the standard CTS theory assumes plasmas to be linear, stationary, isotropic and equilibrium, it is often nonlinear, non-stationary, anisotropic, and non-equilibrium in high energy phenomena relevant to laboratory astrophysics. We theoretically calculate and numerically simulate the CTS spectra in two-stream plasmas as a typical example of non-equilibrium system in space and astrophysical plasmas. The simulation results show the feasibility to diagnose two-stream instability directly via CTS measurements. In order to confirm the non-equilibrium CTS analysis, we have been developing experimental system with high repetition rate table top laser for laboratory astrophysics.

Pulsar Timing Array (PTA) experiments are expected to be sensitive to gravitational waves (GWs) emitted by individual supermassive black hole binaries (SMBHBs) inspiralling along eccentric orbits. We compare the computational cost of different methods of computing the PTA signals induced by relativistic eccentric SMBHBs, namely approximate analytic expressions, Fourier series expansion, post-circular expansion, and numerical integration. We show that the fastest method for evaluating PTA signals is by using the approximate analytic expressions, providing up to a $\sim$50 times improvement in computational performance over the alternative methods. We investigate the accuracy of the approximate analytic expressions by employing a mismatch metric valid for PTA signals. We show that this method is accurate within the region of the binary parameter space that is of interest to PTA experiments. We introduce a spline-based method for further accelerating the PTA signal evaluations for narrow-band PTA datasets. These results are crucial for searching for eccentric SMBHBs in large PTA datasets. We have implemented these results in the GWecc package and can be readily accessed from the popular ENTERPRISE package to search for such signals in PTA datasets.