Papers for Thursday, Dec 08 2022

The vision of 2030STEM is to address systemic barriers in institutional structures and funding mechanisms required to achieve full inclusion in Science, Technology, Engineering, and Mathematics (STEM) and provide leadership opportunities for individuals from underrepresented populations across STEM sectors. 2030STEM takes a systems-level approach to create a community of practice that affirms diverse cultural identities in STEM. This is the first in a series of white papers based on 2030STEM Salons - discussions that bring together visionary stakeholders in STEM to think about innovative ways to infuse justice, equity, diversity, and inclusion into the STEM ecosystem. Our salons identify solutions that come from those who have been most affected by systemic barriers in STEM. Our first salon focused on the power of social media to accelerate inclusion and diversity efforts in STEM. Social media campaigns, such as the #XinSTEM initiatives, are powerful new strategies for accelerating change towards inclusion and leadership by underrepresented communities in STEM. This white paper highlights how #XinSTEM campaigns are redefining community, and provides recommendations for how scientific and funding institutions can improve the STEM ecosystem by supporting the #XinSTEM movement.

Ground-based synoptic solar observations provide critical contextual data used to model the large-scale state of the heliosphere. The next decade will see a combination of ground-based telescopes and space missions that will study our Sun's atmosphere microscopic processes with unprecedented detail. This white paper describes contextual observations from a ground-based network needed to fully exploit this new knowledge of the underlying physics that leads to the magnetic linkages between the heliosphere and the Sun. This combination of a better understanding of small-scale processes and the appropriate global context will enable a physics-based approach to Space Weather comparable to Terrestrial Weather forecasting.

The Schwabe (~11 yr) value for the annual sunspot number is sometimes uncritically applied to other measures of solar activity, direct and indirect, including the 10.7 cm radio flux, the inflow of galactic cosmic rays, solar flare frequency, terrestrial weather, and components of space climate, with the risk of a resulting loss of information. The ruling (Babcock) hypothesis and its derivatives link the sunspot cycle to dynamo processes mediated by differential solar rotation, but despite 60 years of observation and analysis the ~11 yr periodicity remains difficult to model; the possible contribution of planetary dynamics is undergoing a revival. The various solar sequences that genuinely display an ~11 yr cycle stand to benefit from an understanding of its periodicity that goes beyond statistical kinship. The outcome could ironically prompt the demotion of sunspots from their dominant historical role in favour of other possible indicators of solar cyclicity, such as the solar wind flux and its isotopic signatures, even if they are less accessible.

With the IRAM 30-m telescope, we performed observations of the $J$ = 2-1 transitions of CS, C$^{33}$S, C$^{34}$S, C$^{36}$S, $^{13}$CS, $^{13}$C$^{33}$S, and $^{13}$C$^{34}$S as well as the $J$ = 3-2 transitions of C$^{33}$S, C$^{34}$S, C$^{36}$S, and $^{13}$CS toward a large sample of 110 HMSFRs. We measured the $^{12}$C/$^{13}$C, $^{32}$S/$^{34}$S, $^{32}$S/$^{33}$S, $^{32}$S/$^{36}$S, $^{34}$S/$^{33}$S, $^{34}$S/$^{36}$S, and $^{33}$S/$^{36}$S abundance ratios with rare isotopologues of CS, thus avoiding significant saturation effects. With accurate distances obtained from parallax data, we confirm previously identified $^{12}$C/$^{13}$C and $^{32}$S/$^{34}$S gradients as a function of galactocentric distance (RGC). In the CMZ $^{12}$C/$^{13}$C ratios are higher than suggested by a linear fit to the disk values as a function of RGC. While $^{32}$S/$^{34}$S ratios near the Galactic center and in the inner disk are similar, this is not the case for $^{12}$C/$^{13}$C, when comparing central values with those near RGC of 5 kpc. As was already known, there is no $^{34}$S/$^{33}$S gradient but the average ratio of 4.35~$\pm$~0.44, derived from the $J$ = 2-1 transition lines of C$^{34}$S and C$^{33}$S, is well below previously reported values. A comparison between solar and local interstellar $^{32}$S/$^{34}$S and $^{34}$S/$^{33}$S ratios suggests that the solar system may have been formed from gas with a particularly high $^{34}$S abundance. For the first time, we report positive gradients of $^{32}$S/$^{33}$S, $^{34}$S/$^{36}$S, $^{33}$S/$^{36}$S and $^{32}{\rm S}/^{36}{\rm S}$ in our Galaxy. The predicted $^{12}$C/$^{13}$C ratios from the latest GCE models are in good agreement with our results, while $^{32}$S/$^{34}$S and $^{32}$S/$^{36}$S ratios show larger differences at larger RGC, $^{32}$S/$^{33}$S ratios show an offset across the entire inner 12 kpc of the Milky Way.

Early galaxies were the radiation source for reionization, with the photoheating feedback from the reionization process expected to reduce the efficiency of star formation in low mass haloes. Hence, to fully understand reionization and galaxy formation, we must study their impact on each other. The THESAN project has so far aimed to study the impact of galaxy formation physics on reionization, but here we present the new THESAN simulations with a factor 50 higher resolution ($m_{\rm b} \approx 10^4$~M$_\odot$) that aim to self-consistently study the back-reaction of reionization on galaxies. By resolving haloes with virial temperatures $T_{\rm vir} < 10^4$~K, we are able to demonstrate that simplistic, spatially-uniform, reionization models are not sufficient to study early galaxy evolution. Comparing the self-consistent THESAN model (employing fully coupled radiation hydrodynamics) to a uniform UV background, we are able to show that galaxies in THESAN are predicted to be larger in physical extent (by a factor $\sim 2$), less metal enriched (by $\sim 0.2$~dex), and less abundant (by a factor $\sim 10$ at $M_{\rm 1500}~=~-10$) by $z=5$. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation $z \gtrsim 3$ must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties.

In preparation to make the most of our own planned James Webb Space Telescope investigations, we take advantage of publicly available calibration and early-science observations to independently derive and test a geometric-distortion solution for NIRCam detectors. Our solution is able to correct the distortion to better than ~0.2 mas. Current data indicate that the solution is stable and constant over the investigated filters, temporal coverage, and even over the available filter combinations. We successfully tested our geometric-distortion solution matching the JWST and archive HST catalogues. We considered three different applications: (i) cluster-field separation for the stars in the globular cluster M92; (ii) measuring the internal proper motions for M92's stars; (iii) measuring the internal proper motions for the stars in the Large Magellanic Cloud system. While we were not able to detect significant variations of the geometric distortion solution over 22 days, it is clear that more data are still necessary to have a better understanding of the instrument and to characterise the solution to a higher level of accuracy. To our knowledge, the here-derived geometric-distortion solution for NIRCam is the best available and we publicly release it, as many other investigations could potentially benefit from it. Along with our geometric-distortion solution, we also release a Python tool to convert the raw-pixels coordinates of each detector into distortion-free positions, and also to put all the ten detectors of NIRCam into a common reference system.

K-corrections, conversions between flux in observed bands to flux in rest-frame bands, are critical for comparing galaxies at various redshifts. These corrections often rely on fits to empirical or theoretical spectral energy distribution (SED) templates of galaxies. However, the templates limit reliable K-corrections to regimes where SED models are robust. For instance, the templates are not well-constrained in some bands (e.g., WISE W4), which results in ill-determined K-corrections for these bands. We address this shortcoming by developing an empirically-driven approach to K-corrections as a means to mitigate dependence on SED templates. We perform a polynomial fit for the K-correction as a function of a galaxy's rest-frame color determined in well-constrained bands (e.g., rest-frame (g-r)) and redshift, exploiting the fact that galaxy SEDs can be described as a one parameter family at low redshift (0.01 < z < 0.09). For bands well-constrained by SED templates, our empirically-driven K-corrections are comparable to the SED fitting method of Kcorrect and SED template fitting employed in the GSWLC-M2 catalogue (the updated medium-deep GALEX-SDSS-WISE Legacy Catalogue). However, our method dramatically outperforms the available SED fitting K-corrections for WISE W4. Our method also mitigates incorrect template assumptions and enforces the K-correction to be 0 at z = 0. Our K-corrected photometry and code are publicly available.

We present an analysis of multi-epoch spectroscopic and photometric observations of BW Scl, which is believed to be one of the best period-bouncer candidates. We detected multiple irradiation-induced emission lines from the donor star allowing the radial velocity variations to be measured with high accuracy. Also, using the absorption Mgii 4481 line originated in the photosphere of the accreting white dwarf (WD), we measured the radial velocity semi-amplitude of the WD and its gravitational redshift. We find that the WD has a mass of 0.92$\pm$0.04 M$_\odot$, while the donor is a low-mass object with a mass of 0.054$\pm$0.008 M$_\odot$, well below the hydrogen-burning limit. Using NIR data, we put an upper limit on the effective temperature of the donor to be $\leq$1600 K, corresponding to a brown dwarf of T spectral type. The optically thin accretion disc in BW Scl has a very low luminosity $\lesssim$4 $\times 10^{30}$ erg s$^{-1}$ which corresponds to a very low mass accretion rate of $\lesssim$6 $\times 10^{-13}$ M$_\odot$ year$^{-1}$. The outer parts of the disc have a low density allowing the stream to flow down to the inner disc regions. The brightest part of the hotspot is located close to the circularization radius of the disc. The hotspot is optically thick and has a complex, elongated structure. Despite the relatively high temperature of the WD (14750-15000 K), we suggest that BW Scl has already passed the minimum period and is now evolving back towards longer periods. Thus, BW Scl is a period bouncer.

Large-scale outflows are believed to be an important mechanism in the evolution of galaxies. We can determine the impact of these outflows by studying either current galaxy outflows and their effect in the galaxy or by studying the effect of past outflows on the gas surrounding the galaxy. In this work, we examine the CO(7-6), [CI]\,($^{3} \rm P_{1} \rightarrow {\rm ^3 P}_{0}$), H$_2$O 2$_{11}$--2$_{02}$ and dust continuum emission of 15 extremely red quasars (ERQs) at z$\sim$2.3 using ALMA. By investigating the radial surface brightness profiles of both the individual sources and the stacked emission, we detect extended cold gas and dust emission on scales of $\sim$14 kpc in CO(7-6), [CI](2-1), and dust continuum. This is the first time that the presence of a large amount of molecular gas was detected on large, circum-galactic medium scales around quasar host galaxies using [CI] extended emission. We estimate the dust and molecular gas mass of these halos to be 10$^{7.6}$ and 10$^{10.6}$ M$_\odot$, indicating significant dust and molecular gas reservoirs around these extreme quasars. By estimating the timescale at which this gas can reach these distances by molecular gas outflows (7-32 Myr), we conclude that these halos are a relic of past AGN or starburst activity, rather than an effect of the current episode of extreme quasar activity.

The central engine of GRB170817A post-merger to GW170817 is probed by GW-calorimetry and event timing, applied to a post-merger descending chirp which can potentially break the degeneracy in spin-down of a neutron star or black hole remnant by the relatively large energy reservoir in the angular momentum, $E_J$, of the latter according to the Kerr metric. This analysis derives from model-agnostic spectrograms with equal sensitivity to ascending and descending chirps generated by time-symmetric butterfly matched filtering. The sensitivity was calibrated by response curves generated by software injection experiments. The statistical significance for candidate emission from the central engine of GRB170817A is expressed by probabilities of false alarm (PFA; type I errors) derived from an event-timing analysis. PDFs were derived for start-time $t_s$, identified via high-resolution image analyses of the available spectrograms. For merged (H1,L1)-spectrograms of the LIGO detectors, a PFA $p_1$ derives from causality in $t_s$ given GW170817-GRB17081A. A statistically independent confirmation is presented in individual H1 and L1 analyses, in a second PFA $p_2$ of consistency in their respective observations of $t_s$. A combined PFA derives from their product since mean and (respectively) difference in timing are statistically independent. Applied to GW170817-GRB170817A, PFAs of event timing in $t_s$ produce $p_1=8.3\times 10^{-4}$ and $p_2=4.9\times 10^{-5}$ of a post-merger output ${\cal E}_{GW}\simeq 3.5\%M_\odot c^2$ ($p_1p_2=4.1\times 10^{-8}$, equivalent $Z$-score 5.48). ${\cal E}_{GW}$ exceeds $E_J$ of the hyper-massive neutron star in the immediate aftermath of GW170817, yet it is consistent with $E_J$ rejuvenated in delayed gravitational collapse to a Kerr black hole. Similar emission may be expected from energetic core-collapse supernovae producing black holes. (Abbr.)

Spectroscopic detection of transient narrow emission lines (flash-ionisation features) traces the presence of circumstellar material (CSM) around massive stars exploding as core-collapse supernovae. Transient emission lines disappearing shortly after the SN explosion suggest that the spatial extent of this material is compact; hence implying that the progenitor star experienced episodes of enhanced mass loss shortly prior to explosion. The early light curves of Type II supernovae (SNe II) are assumed to be initially powered by shock-cooling emission. Additional luminosity may arise from interaction via shocks with the CSM, if it is present. We performed a systematic survey of SNe II discovered within less than two days from explosion during the first phase of the Zwicky Transient Facility (ZTF) survey. We gathered the early light curves and spectra of thirty SNe II. The measured fraction of events showing emission line evidence for CSM (>30% at 95% confidence level) indicates that elevated mass loss is a common process occurring in massive stars. We also measure the rise time and peak magnitude of each event. We find that SNe II showing spectroscopic evidence for CSM interaction at early time are not significantly brighter, nor bluer, nor more slowly rising than those who do not. This implies that the CSM in these events is likely optically thin, and therefore that the CSM interaction does not contribute significantly to their early continuum emission. We also introduce for the first time a measurement of the timescale of appearance of flash ionisation features. Most SN show flash features for ~5 days. A rarer population of events with timescales longer than 10 days, seem to be brighter and rise longer, thus making them a potential bridging population between regular SNe II and strongly-interacting SNe IIn.

We present radial gas-phase metallicity profiles, gradients, and break radii at redshift $z = 0 - 3$ from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, $r_{\rm Break}$, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of $\sim$2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that $r_{\rm Break}$ occurs where the gas mixing time-scale is $\sim$10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic disks can be thought of in two-parts: a star-forming inner disk with a steep gradient and a mixing-dominated outer disk with a flat gradient, with the break radius marking the region of transition between them.

Primordial black hole (PBHs) is interesting to people for its ability of explaining dark matter as well as supermassive astrophysical objects. In normal inflation scenario, the generation of PBHs usually requires enhanced power spectrum of scalar perturbation at the end of inflation era, which is expected when the dispersion relation of the inflaton field gets modified. In this work, we study a kind of inflation model called "{\it DBI-inspired non-minimal kinetic coupling" (DINKIC)} model, where the dispersion relation is modified by a square root existing in the field Lagrangian. We discuss the enhancement of scalar power spectrum due to the modified dispersion relation, as well as the abundance of PBHs produced by the Press-Schechter collapse mechanism. We also discuss the formation of Scalar-Induced Gravitational Waves (SIGWs) by linear scalar perturbations.

Type Ibn supernovae (SNe) are a rare class of stellar explosions whose progenitor systems are not yet well determined. We present and analyze observations of the Type Ibn SN 2019kbj, and model its light curve in order to constrain its progenitor and explosion parameters. SN 2019kbj shows roughly constant temperature during the first month after peak, indicating a power source (likely CSM interaction) that keeps the continuum emission hot at ~15000K. Indeed, we find that the radioactive decay of Ni56 is disfavored as the sole power source of the bolometric light curve. A radioactive decay + circumstellar-material (CSM) interaction model, on the other hand, does reproduce the bolometric emission well. The fits prefer a uniform-density CSM shell rather than CSM due to a steady mass-loss wind, similar to what is seen in other Type Ibn SNe. The uniform-density CSM shell model requires ~0.1 solar masses of Ni56 and ~1 solar mass of total ejecta to reproduce the light curve. SN 2019kbj differs in this manner from another Type Ibn SN with derived physical parameters, SN 2019uo, for which an order of magnitude lower Ni56 mass and larger ejecta mass were derived. This points towards a possible diversity in SN Ibn progenitor systems and explosions.

Solar macrospicules are beam-like cool plasma ejections of size in-between spicules and coronal jets, which can elucidate potential connections between plasma jetting activity at different scales. With high-resolution observations from the {\em New Vacuum Solar Telescope} and Solar Dynamic Observatory, we investigate the origin of five groups of recurrent active-region macrospicules. Before the launch of each macrospicule, we detect a compact bright patch (BP) at its base where a newly emerging dipole contacts and cancel with the pre-existing ambient field. The spectral diagnosis from the {\em Interface Region Imaging Spectrograph} at one of BPs reveals signatures of reconnection at the lower atmosphere. Multiwavelength imaging of these BPs show that they mainly occur at the rising phase of the flux emergence and slowly ascend from the lower to the upper chromosphere. Remarkable macrospicules occur and fade out once the BPs appear and decay from the AIA 304 A images, respectively. We suggest that these macrospicules and related BPs form in a common reconnection process, in which the increasing reconnection height between the emerging dipole and the ambient field results in the observed variations from BPs to macrospicules. Interestingly, most macrospicules show similar characteristics to larger-scale coronal jets and/or smaller-scale spicules, i.e., the rotating motions, the presence of minifilaments and BPs before the eruptions, and magnetic flux emergence and cancellation. We conclude that the formation mechanism of macrospicules should be the same as spicules and coronal jets, i.e., solar jetting phenomena at different scales share the same physical mechanism in association with magnetic reconnection.

We report the discovery of a rare early-type dwarf galaxy (dE), SDSS J125651.47+163024.2 (hereafter dE1256), possessing a tidal feature that was likely built up by accretion of an even smaller dwarf galaxy. dE1256 is located in a nearly isolated environment, at the outskirt of the Virgo cluster. A detailed morphological examination reveals that the accreted stellar population is mainly deposited in the outer part of dE1256, where the tidal tail is most prominent. The inner part of dE1256 is perfectly modeled with a simple S\'ersic function of index n = 0.63 and half-light radius R$_{h}$ = 0.6 kpc, but in contrast, the entire galaxy has a size of R$_{h}$ = 1.2 kpc. The mass ratio between the host and the putative accreted dwarf galaxy is calculated to be 5:1, assuming that the observed two components, inner S\'ersic, and outer tidal tail residual, represent the host's and accreted galaxy's stellar populations, respectively. We suggest that while the accretion contributes only 20% of the overall stellar population, the size of dE1256 grew by a factor of two via the accretion event. Our results provide, for the first time, strong observational evidence that a dE is undergoing a two-phase growth, a common phenomenon for massive galaxies.

Frequent, low-latency measurements of the Earth's rotation phase, UT1$-$UTC, critically support the current estimate and short-term prediction of this highly variable Earth Orientation Parameter (EOP). Very Long Baseline Interferometry (VLBI) Intensive sessions provide the required data. However, the Intensive UT1$-$UTC measurement accuracy depends on the accuracy of numerous models, including the VLBI station position. Intensives observed with the Maunakea (Mk) and Pie Town (Pt) stations of the Very Long Baseline Array (VLBA) illustrate how a geologic event (i.e., the $M_w$ 6.9 Hawai`i Earthquake of May 4th, 2018) can cause a station displacement and an associated offset in the values of UT1$-$UTC measured by that baseline, rendering the data from the series useless until it is corrected. Using the non-parametric Nadaraya-Watson estimator to smooth the measured UT1$-$UTC values before and after the earthquake, we calculate the offset in the measurement to be 75.7 $\pm$ 4.6 $\mu$s. Analysis of the sensitivity of the Mk-Pt baseline's UT1$-$UTC measurement to station position changes shows that the measured offset is consistent with the 67.2 $\pm$ 5.9 $\mu$s expected offset based on the 12.4 $\pm$ 0.6 mm total coseismic displacement of the Maunakea VLBA station determined from the displacement of the co-located global navigation satellite system (GNSS) station. GNSS station position information is known with a latency on the order of tens of hours, and thus can be used to correct the a priori position model of a co-located VLBI station such that it can continue to provide accurate measurements of the critical EOP UT1$-$UTC as part of Intensive sessions. The VLBI station position model would likely not be updated for several months. This contrast highlights the benefit of co-located GNSS and VLBI stations in support of the monitoring of UT1$-$UTC with single baseline Intensives. Abridged.

The dwarf irregular galaxy HIPASS J1131-31 was discovered as a source of HI emission at low redshift in such close proximity of a bright star that we call it Peekaboo. The galaxy resolves into stars in images with Hubble Space Telescope, leading to a distance estimate of 6.8+-0.7 Mpc. Spectral optical observations with the Southern African Large Telescope reveal HIPASS J1131-31 to be one of the most extremely metal-poor galaxies known with the gas-phase oxygen abundance 12+log(O/H) = 6.99+-0.16 dex via the direct [OIII] 4363 line method and 6.87+-0.07 dex from the two strong line empirical methods. The red giant branch of the system is tenuous compared with the prominence of the features of young populations in the color-magnitude diagram, inviting speculation that star formation in the galaxy only began in the last few Gyr.

Sources of high-energy cosmic rays are presently unknown, but can be constrained in various ways. Some of these constraints can be graphically presented on the so-called Hillas diagram. Previous versions of this diagram determined the range of geometrical sizes and magnetic fields of potential astrophysical accelerators, taking into account geometrical criteria and radiation losses. In this work, we update the Hillas diagram for protons, taking into account the losses associated with the $p\gamma$ interaction, relating the allowed regions to the source electromagnetic luminosity. The strongest constraints are obtained for bright compact sources such as the central regions of active galactic nuclei.

It is clearly established that the Sun has an 11-year cycle that is caused by its internal magnetic field. This cycle is also observed in a sample of M dwarfs. In the framework of exoplanet detection or atmospheric characterisation of exoplanets, the activity status of the host star plays a crucial role, and inactive states are preferable for such studies. This means that it is important to know the activity cycles of these stars. We study systematic long-term variability in a sample of 211 M dwarfs observed with CARMENES, the high-resolution optical and near-infrared spectrograph at Calar Alto Observatory. In an automatic search using time series of different activity indicators, we identified 26 stars with linear or quadratic trends or with potentially cyclic behaviour. Additionally, we performed an independent search in archival R$^{\prime}_{\rm HK}$ data collected from different instruments whose time baselines were usually much longer. These data are available for a subset of 186 of our sample stars. Our search revealed 22 cycle candidates in the data. We found that the percentage of stars showing long-term variations drops dramatically to the latest M dwarfs. Moreover, we found that the pseudo-equivalent width (pEW) of the H$\alpha$ and Ca ii infrared triplet more often triggers automatic detections of long-term variations than the TiO index, differential line width, chromatic index, or radial velocity. This is in line with our comparison of the median relative amplitudes of the different indicators. For stars that trigger our automatic detection, this leads to the highest amplitude variation in R$^{\prime}_{\rm HK}$, followed by pEW(H$\alpha$), pEW(Ca ii IRT), and the TiO index.

Based on the rotating turbulent thermal convection model and using the rotating equivalent temperature assumption and new convection criterion, this paper analyzed the repression of the vorticity gradient on the heat transport and explained that the formation of sunspots originated from the variation of the vorticity gradient in the solar troposphere.

We present the first results from the Southern Hemisphere Parallax Interferometric Radio Astrometry Legacy Survey (\spirals): $10\mu$as-accurate parallaxes and proper motions for two southern hemisphere 6.7 GHz methanol masers obtained using the inverse MultiView calibration method. Using an array of radio telescopes in Australia and New Zealand, we measured the trigonometric parallax and proper motions for the masers associated with the star formation regions G232.62+00.99 of $\pi = 0.610\pm0.011$~mas, $\mu_x=-2.266\pm0.021$~mas/yr and $\mu_y=2.249\pm0.049$~mas/yr, which implies its distance to be $d=1.637\pm0.029$~kpc. These measurements represent an improvement in accuracy by more than a factor of 3 over the previous measurements obtained through Very Long Baseline Array observations of the 12~GHz methanol masers associated with this region. We also measure the trigonometric parallax and proper motion for G323.74--00.26 as $\pi = 0.364\pm0.009$~mas, $\mu_x=-3.239\pm0.025$~mas/yr and $\mu_y=-3.976\pm0.039$~mas/yr, which implies a distance of $d=2.747\pm0.068$~kpc. These are the most accurate measurements of trigonometric parallax obtained for 6.7~GHz class II methanol masers to date. We confirm that G232.62+00.99 is in the Local arm and find that G323.74--00.26 is in the Scutum-Centaurus arm. We also investigate the structure and internal dynamics of G323.74--00.26

Diffuse interstellar bands comprise hundreds of absorption features in the ISM. Most DIBs are observed in the optical, but some are in the IR. We observed 76 early-type stars at R=50,000 and S/N ratios of several hundreds using CRIRES. We measure DIBs around 1318, 1527, 1561, 1565, 1567, 1574 and/or 1624 nm. We detect a total of 6 DIB features and 17 likely stellar features assisted by a CMFGEN model. We also measured the DIBs at 1318 and 1527 nm using X-shooter towards ten Ceph. variables with 3.2 < E(B-V) < 6.5 and 4 stars at low values of water vapour. Correlation coeffs. of 0.73-0.96 are found comparing NIRDIB eq. width vs. E(B-V) and with r > 0.8 when comparing the NIR and optical DIBs 5705, 5780, 6203, 6283 and 6269 A. The 5797 A DIB is less well correlated with the NIDIBs. The "C60+" DIB at 9632 A is not well correlated with the 1318 nm DIB. Partial correlation coefficients using E(B-V) as the covariate were also determined. For stars earlier than B2, the 1318 nm DIB is affected by an emission line on its blue wing, likely stellar in nature, although we cannot rule out interstellar/circumstellar origin for example caused by by a DIB in emission. The 1318 nm DIB has a red wing and is reasonably well fitted by two gaussians. Neither the component ratios nor separation are correlated with 5780/5797 or E(B-V). EW(1318 nm) correlates with HI with EW(1318 nm)/E(B-V) decreasing with f(H2). Five pairs of stars within 1 am show similar 1318 nm DIB profiles. Variation in 1318 nm is seen in HD 145501/145502 and HD 168607/168625 pairs. CRIRES data for 17 stars separated by 6-14 months and 2 X-shooter sightlines separated by 9.9 yr were analysed. No time-variability is detected in the 5780, 5797 A, 1318 nm or 1527 nm DIBs. Tentative time variation is observed in the C60+ DIBs at 9577 and 9632 A towards HD 183143 although very close to the noise level with confirmation required.

In this paper, we propose a model in which a spectator field non-minimally couples to an inflaton field and the power spectrum of the perturbation of the spectator field at small scales is dramatically enhanced by the sharp feature in the form of non-minimal coupling. At or after the end of inflation, the perturbation of the spectator field is converted into curvature perturbation and leads to the formation of primordial black holes (PBHs). Furthermore, for example, we consider three phenomenological models for generating PBHs with mass function peaked at $\sim10^{-12}M_\odot$ and representing all the cold dark matter in our Universe and find that the scalar induced gravitational waves generated by the curvature perturbation can be detected by the future space-borne gravitational-wave detectors such as Taiji, TianQin and LISA.

We present a new gamma ray energy reconstruction method based on Random Forest to be commonly used for the data analysis of the MAGIC Telescopes, a system of two Imaging Atmospheric Cherenkov Telescopes. The energy resolution with the new energy reconstruction improves compared to the one obtained with the LUTs method. For standard observations i.e. dark conditions with pointing zenith (Zd) less than 35 deg for a point-like source, the energy resolution goes from $\sim 20\%$ at 100 GeV to $\sim 10\%$ at a few TeV. In addition, the new method suppresses the outlier population in the energy error distribution, which is thus better described by a Gaussian distribution. The new energy reconstruction method enhances the reliability especially for the sources with steep spectra, in higher energies and/or in observations at higher Zd pointings. We validate the new method in different ways and demonstrate some cases of its remarkable benefit in spectral analysis with simulated observation data.

We aim to analyse the abundance pattern of 169 Barium (Ba) stars, using machine learning techniques and the AGB final surface abundances predicted by Fruity and Monash stellar models. We developed machine learning algorithms that use the abundance pattern of Ba stars as input to classify the initial mass and metallicity of its companion star using stellar model predictions. We use two algorithms: the first exploits neural networks to recognise patterns and the second is a nearest-neighbour algorithm, which focuses on finding the AGB model that predicts final surface abundances closest to the observed Ba star values. In the second algorithm we include the error bars and observational uncertainties to find the best fit model. The classification process is based on the abundances of Fe, Rb, Sr, Zr, Ru, Nd, Ce, Sm, and Eu. We selected these elements by systematically removing s-process elements from our AGB model abundance distributions, and identifying those whose removal has the biggest positive effect on the classification. We excluded Nb, Y, Mo, and La. Our final classification combines the output of both algorithms to identify for each Ba star companion an initial mass and metallicity range. With our analysis tools we identify the main properties for 166 of the 169 Ba stars in the stellar sample. The classifications based on both stellar sets of AGB final abundances show similar distributions, with an average initial mass of M = 2.23 MSun and 2.34 MSun and an average [Fe/H] = -0.21 and -0.11, respectively. We investigated why the removal of Nb, Y, Mo, and La improves our classification and identified 43 stars for which the exclusion had the biggest effect. We show that these stars have statistically significant different abundances for these elements compared to the other Ba stars in our sample. We discuss the possible reasons for these differences in the abundance patterns.

The MRS is one of the four observing modes of JWST/MIRI. Using JWST in-flight data of unresolved (point) sources, we can derive the MRS absolute spectral response function (ASRF) starting from raw data. Spectral fringing plays a critical role in the derivation and interpretation of the MRS ASRF. In this paper, we present an alternative way to calibrate the data. Firstly, we aim to derive a fringe correction that accounts for the dependence of the fringe properties on the MIRI pupil illumination and detector pixel sampling of the point spread function. Secondly, we aim to derive the MRS ASRF using an absolute flux calibrator observed across the full 5 to 28 $\mu$m wavelength range of the MRS. Thirdly, we aim to apply the new ASRF to the spectrum of a G dwarf and compare with the output of the JWST/MIRI default data reduction pipeline. Finally, we examine the impact of the different fringe corrections on the detectability of molecular features in the G dwarf and K giant. The absolute flux calibrator HD 163466 (A-star) is used to derive tailored point source fringe flats at each of the default dither locations of the MRS. The fringe-corrected point source integrated spectrum of HD 163466 is used to derive the MRS ASRF using a theoretical model for the stellar continuum. A cross-correlation is run to quantify the uncertainty on the detection of CO, SiO, and OH in the K giant and CO in the G dwarf for different fringe corrections. The point-source-tailored fringe correction and ASRF are found to perform at the same level as the current corrections, beating down the fringe contrast to the sub-percent level, whilst mitigating the alteration of real molecular features. The same tailored solutions can be applied to other MRS unresolved targets. A pointing repeatability issue in the MRS limits the effectiveness of the tailored fringe flats is at short wavelengths.

It is known that gap opening depends on the disc's viscosity; however, eccentricity damping formulas have only been derived at high viscosities, ignoring partial gap opening. We aim at obtaining a simple formula to model $e$-damping of the type-I regime in low viscosity discs, where even small planets may start opening partial. We perform high resolution 2D locally isothermal hydrodynamical simulations of planets with varying masses on fixed orbits in discs with varying aspect ratios and viscosities. We determine the torque and power felt by the planet to derive migration and eccentricity damping timescales. We first find a lower limit to the gap depths below which vortices appear; this happens roughly at the transition between type-I and type-II regimes. For the simulations that remain stable, we obtain a fit to the observed gap depth in the limit of vanishing eccentricities that is similar to the one currently used in the literature but is accurate down to $\alpha=3.16\times 10^{-5}$. We record the $e$-damping efficiency as a function of the observed gap depth and $e$: when the planet has opened a deep enough gap, a linear trend is observed independently of $e$; at shallower gaps this linear trend is preserved at low $e$, while it deviates to more efficient damping when $e$ is comparable to the disc's scale height. Both trends can be understood on theoretical grounds and are reproduced by a simple fitting formula. Our combined fits yield a simple recipe to implement type-I $e$-damping in $N$-body for partial gap opening planets that is consistent with high-resolution 2D hydro-simulations. The typical error of the fit is of the order of a few percent, and lower than the error of type-I torque formulas widely used in the literature. This will allow a more self-consistent treatment of planet-disc interactions of the type-I regime for population synthesis models at low viscosities.

High-resolution Hubble Space Telescope optical observations have been used to analyze the stellar population and the structure of the poorly investigated bulge globular cluster NGC 6316. We constructed the first high-resolution reddening map in the cluster direction, which allowed us to correct the evolutionary sequences in the color magnitude diagram (CMD) for the effects of differential reddening. A comparison between the CMDs of NGC 6316 and 47 Tucanae revealed strikingly similar stellar populations, with the two systems basically sharing the same turn-off, sub-giant branch, and horizontal branch morphologies, indicating comparable ages. The red giant branch in NGC 6316 appears slightly bluer than in 47 Tucanae, suggesting a lower metal content. This has been confirmed by the isochrone fitting of the observed CMD, which provided us with updated values of the cluster age, distance, average color excess, and metallicity. We estimated an absolute age of 13.1 $\pm$ 0.5 Gyr, consistent with the age of 47 Tucanae, an average color excess E(B-V) = 0.64 $\pm$ 0.01, and a true distance modulus (m-M)0 = 15.27 $\pm$ 0.03 that sets the cluster distance at 11.3 kpc from the Sun. In addition, the photometric estimate of the cluster metallicity suggests [Fe/H]$\approx$ -0.9, which is $\sim$ 0.2 dex smaller than that of 47 Tucanae. We also determined the gravitational center and the density profile of the system from resolved stars. The latter is well reproduced by a King model. Our results confirm that NGC 6316 is another extremely old relic of the assembly history of the Galaxy.

We present multi object optical spectroscopy of the galaxies in the environment of 12 low-redshift (z < 0.5) quasars and of 11 inactive massive galaxies chosen to match the properties of the quasar host galaxies to probe physical association and possible events of recent star formation. The quasars are selected from a sample of QSOs in the SDSS Stripe82 region for which both the host galaxy and the large scale environments were previously investigated. The new observations complement those reported in our previous works on close companion galaxies of nearby quasars. For the whole dataset we find that for about half (19 out of 44 ) of the observed QSOs there is at least one associated companion galaxy. In addition to the new spectroscopic observations, we add data from the SDSS database for the full sample of objects. We find that the incidence of companion galaxies in the fields of QSO (17%) is not significantly different from that of inactive galaxies (19%) similar to quasar hosts in redshift and mass. Nevertheless, the companions of quasars exhibit more frequently emission lines than those of inactive galaxies, suggesting a moderate link between the nuclear activity and recent star formation in their environments.

The Atmospheric Remote-sensing Infrared Exoplanet Large-survey (Ariel) is the first space mission dedicated to measuring the chemical composition and thermal structures of thousands of transiting exoplanets. Ariel was adopted in 2020 as the M4 mission in ESA "Cosmic Vision" program, with launch expected in 2029. The mission will operate from the Sun-Earth Lagrange Point L2. The scientific payload consists of two instruments: a high resolution spectrometer in the waveband 1.95-7.8 microns, and a fine guidance system / visible photometer / low resolution near-infrared spectrometer. The instruments are fed a collimated beam from an unobscured, off-axis Cassegrain telescope. Instruments and telescope will operate at a temperature below 50 K. Telescope mirrors and supporting structures will be realized in aerospace-grade aluminum. Given the large aperture of the primary mirror (0.6 m$^2$), it is a choice of material that requires careful optical and opto-mechanical design, and technological advances in the three areas of mirror substrate thermal stabilization, optical surface polishing and optical coating. This thesis presents the work done by the author in these areas, as member of the team responsible for designing and manufacturing the telescope and mirrors, starting with a systematic review of the optical and opto-mechanical requirements and design choices of the Ariel telescope, in the context of previous development work and scientific goals and requirements of the mission. The review then progresses with opto-mechanical design, examining the most important choices in terms of structural and thermal design, and with a statistical analysis of the deformations of the optical surface of the telescope mirrors and of their alignment in terms of rigid body motions. The details of the qualification work on thermal stabilization, polishing and coating are then presented.

Liquid water is a critical component of habitability. However, the production and stability of surficial liquid water can be challenging on planets outside the Habitable Zone and devoid of adequate greenhouse warming. On such cold, icy exoEarths, basal melting of regional, global ice sheets by geothermal heat provides an alternative means of forming liquid water. Here, we model the thermophysical evolution of ice sheets to ascertain the geophysical conditions that allow liquid water to be produced and maintained at temperatures above the pressure controlled freezing point of water ice on exoEarths. We show that even with a modest, Moon like geothermal heat flow, subglacial oceans of liquid water can form at the base of and within the ice sheets on exoEarths. Furthermore, subglacial oceans may persist on exoEarths for a prolonged period due to the billion year half lives of heat producing elements responsible for geothermal heat. These subglacial oceans, often in contact with the planets crust and shielded from the high energy radiation of their parent star by thick ice layers, may provide habitable conditions for an extended period.

Recently, we proposed a varying speed of light (VSL) model to solve various late-time cosmological problems \cite{Lee:2020zts}. This model has one free parameter, $b$, to characterize the time variation of the speed of light as a function of a scale factor, $c = c_0a^{b/4}$. Time variations of various physical constants and quantities have different powers of scale factor as a function of $b$ to satisfy all known local physics laws, including special relativity, thermodynamics, and electromagnetic force. This model is based on the Robertson-Walker metric and satisfies the isotropic and homogeneous 3-space required by the cosmological principle. Adiabaticity is a necessary condition to keep homogeneity and isotropy because a net energy flux would falsify the isotropy if there is a preferential energy flow direction. It also might forge homogeneity if the outward (inward) flux is isotropic. Thus, any VSL model should also preserve an adiabatic expansion condition to be a viable model. Also, it provides an additional condition for constraining physical constants.

Recent cosmological analyses with large-scale structure and weak lensing measurements, usually referred to as 3$\times$2pt, had to discard a lot of signal-to-noise from small scales due to our inability to precisely model non-linearities and baryonic effects. Galaxy-galaxy lensing, or the position-shear correlation between lens and source galaxies, is one of the three two-point correlation functions that are included in such analyses, usually estimated with the mean tangential shear. However, tangential shear measurements at a given angular scale $\theta$ or physical scale $R$ carry information from all scales below that, forcing the scale cuts applied in real data to be significantly larger than the scale at which theoretical uncertainties become problematic. Recently there have been a few independent efforts that aim to mitigate the non-locality of the galaxy-galaxy lensing signal. Here we perform a comparison of the different methods, including the Y transformation described in Park et al. (2021), the point-mass marginalization methodology presented in MacCrann et al. (2020) and the Annular Differential Surface Density statistic described in Baldauf et al. (2010). We do the comparison at the cosmological constraints level in a noiseless simulated combined galaxy clustering and galaxy-galaxy lensing analysis. We find that all the estimators perform equivalently using a Rubin Observatory Legacy Survey of Space and Time (LSST) Year 1 like setup. This is because all the estimators project out the mode responsible for the non-local nature of the galaxy-galaxy lensing measurements, which we have identified as $1/R^2$. We finally apply all the estimators to DES Y3 data and confirm that they all give consistent results.

Star-formation is one of the main processes that shape galaxies, defining its stellar population and metallicity production and enrichment. It is nowadays known that this process is ruled by a set of relations that connect three parameters: the molecular gas mass, the stellar mass and the star-formation rate itself. These relations are fulfilled at a wide range of scales in galaxies, from galaxy wide to kpc-scales. At which scales they are broken, and how universal they are (i.e., if they change at different scales or for different galaxy types) it is still an open question. We explore here how those relations compare at different scales using as proxy the new analysis done using Integral Field Spectroscopy data and CO observations data from the EDGE-CALIFA survey and the AMUSSING++ compilation.

During its 53-day polar orbit around Jupiter, Juno often crosses the boundaries of the Jovian magnetosphere (namely the magnetopause and bow shock). From the boundary locations, the upstream solar wind dynamic pressure can be inferred, which in turn illustrates the state of compression or relaxation of the system. The aim of this study is to examine Jovian radio emissions during magnetospheric compressions, in order to determine the relationship between the solar wind and Jovian radio emissions. In this paper, we give a complete list of bow shock and magnetopause crossings (from June 2016 to August 2022), along with some extra informations (e.g. solar wind dynamic pressure and position of the standoff distances inferred from Joy et al. (2002)). We then select two compression events that occur in succession (inferred from magnetopause crossings) and we present a case study of the response of the Jovian radio emissions. We demonstrate that magnetospheric compressions lead to the activation of new radio sources. Newly activated broadband kilometric emissions are observed almost simultaneously to compression of the magnetosphere, with sources covering a large range of longitudes. Decametric emission sources are seen to be activated more than one rotation later only at specific longitudes and dusk local times. Finally, the activation of narrowband kilometric radiation is not observed during the compression phase, but when the magnetosphere is in its expansion phase.

Gas phase Elemental abundances in molecular CloudS (GEMS) is an IRAM 30m large program aimed at determining the elemental abundances of carbon (C), oxygen (O), nitrogen (N), and sulfur (S) in a selected set of prototypical star-forming filaments. In particular, the elemental abundance of S remains uncertain by several orders of magnitude and its determination is one of the most challenging goals of this program. We have carried out an extensive chemical modeling of the fractional abundances of CO, HCO$^+$, HCN, HNC, CS, SO, H$_2$S, OCS, and HCS$^+$ to determine the sulfur depletion toward the 244 positions in the GEMS database. These positions sample visual extinctions from A$_V$ $\sim$ 3 mag to $>$50 mag, molecular hydrogen densities ranging from a few 10$^3$~cm$^{-3}$ to 3$\times$10$^6$~cm$^{-3}$, and T$_k$ $\sim$ 10$-$35 K. Most of the positions in Taurus and Perseus are best fitted assuming early-time chemistry, t=0.1 Myr, $\zeta_{H_2}$$\sim$ (0.5$-$1)$\times$10$^{-16}$ s$^{-1}$, and [S/H]$\sim$1.5$\times$10$^{-6}$. On the contrary, most of the positions in Orion are fitted with t=1~Myr and $\zeta_{H_2}$$\sim$ 10$^{-17}$ s$^{-1}$. Moreover, $\sim$40% of the positions in Orion are best fitted assuming the undepleted sulfur abundance, [S/H]$\sim$1.5$\times$10$^{-5}$. Our results suggest that sulfur depletion depends on the environment. While the abundances of sulfur-bearing species are consistent with undepleted sulfur in Orion, a depletion factor of $\sim$20 is required to explain those observed in Taurus and Perseus. We propose that differences in the grain charge distribution in the envelopes of the studied clouds might explain these variations. The shocks associated with past and ongoing star formation could also contribute to enhance [S/H] in Orion.

Models of the protostellar source, B335, are developed using axisymmetric three-dimensional models to resolve conflicts found in one-dimensional models. The models are constrained by a large number of observations, including ALMA, Herschel, and Spitzer data. Observations of the protostellar source B335 with ALMA show red-shifted absorption against a central continuum source indicative of infall in the HCO$^+$ and HCN $J = 4\rightarrow 3$ transitions. The data are combined with a new estimate of the distance to provide strong constraints to three-dimensional radiative transfer models including a rotating, infalling envelope, outflow cavities, and a very small disk. The models favor ages since the initiation of collapse between $3 \times 10^4$ and $4 \times 10^4$ yr for both the continuum and the lines, resolving a conflict found in one-dimensional models. The models under-predict the continuum emission seen by ALMA, suggesting an additional component such as a pseudo-disk. The best-fitting model is used to convert variations in the 4.5 $\mu m$ flux in recent years into a model for a variation of a factor of 5-7 in luminosity over the last 8 years.

Aims. The project aims to understand better the role of wide brown dwarf companions on planetary systems. Methods. We obtained high-resolution spectra of six bright stars with co-moving wide substellar companions with the SONG, CARMENES, and STELLA high-resolution spectrographs. We used these spectra to derive radial velocities together with a complete set of stellar physical parameters. We then investigated radial velocities signals and discussed the fraction of planets in such systems. We also re-analyzed the ages of our targets, which were used to derive the physical parameters of wide brown dwarf companions. Finally, a compilation of systems with known planets from the literature is considered along with our sample to search for possible peculiarities in their parameter distributions. Results. Based on the derived ages of six observed systems, we re-computed the masses of the wide companions, confirming their substellar nature. We confirmed planets in the HD 3651 and HIP 70849 systems and found a new planetary candidate in the HD 46588 system. In our survey, which is sensitive mostly to Neptune-mass planets at short periods of a few days and Saturn-mass planets at longer periods of hundreds of days, we derived a frequency of planets orbiting stars with wide brown dwarf companions below 70% with the uncertainties included. Comparing the parameter distributions of our sample with single stars, we observe the enhancement of planets with short periods below six days in systems with a wide stellar companion. Finally, planets in systems with wide BD companions follow their own eccentricity distribution with a maximum at $\sim0.65$ and have periods larger than 40 days, masses larger than $0.1\,M_J$, and eccentricities larger than 0.4.

21-cm HI4PI survey data are used to study the anomalous-velocity hydrogen gas associated with high-velocity cloud Complex M. These high-sensitivity, high-resolution, high-dynamic-range data show that many of the individual features, including MI, MIIa, and MIIb, are components of a long, arched filament that extends from about (l, b) = (105{\deg}, 53{\deg}) to (l, b) = (196{\deg}, 55{\deg}). Maps at different velocities, results from Gaussian analysis, and observations of associated high-energy emission make a compelling case that the MI cloud and the arched filament are physically interacting. If this is the case, we can use the distance to MI, 150 pc as reported by Schmelz & Verschuur (2022), to set the distance to Complex M. The estimated mass of Complex M is then about 120 solar masses and the energy implied using the observed line-of-sight velocity, -85 km/s, is 8.4 x 10^48 ergs. Integrating over 4{\pi} steradians, the total energy for a spherically symmetrical explosion is estimated to be 1.9 x 10^50 ergs, well within the energy budget of a typical supernova.

In 2016-2017, TUS, the world's first experiment for testing the possibility of registering ultra-high energy cosmic rays (UHECRs) by their fluorescent radiation in the night atmosphere of Earth was carried out. Since 2019, the Russian-Italian fluorescence telescope (FT) Mini-EUSO ("UV Atmosphere") has been operating on the ISS. The stratospheric experiment EUSO-SPB2, which will employ an FT for registering UHECRs, is planned for 2023. We show how a simple convolutional neural network can be effectively used to find track-like events in the variety of data obtained with such instruments.

Many hot and ultra-hot Jupiters have inflated radii, implying that their interiors retain significant entropy from formation. These hot interiors lead to an enhanced internal heat flux that impinges upon the atmosphere from below. In this work, we study the effect of this hot interior on the atmospheric circulation and thermal structure of hot and ultra-hot Jupiters. To do so, we incorporate the population-level predictions from evolutionary models of hot and ultra-hot Jupiters as input for a suite of General Circulation Models (GCMs) of their atmospheric circulation with varying semi-major axis and surface gravity. We conduct simulations with and without a hot interior, and find that there are significant local differences in temperature of up to hundreds of Kelvin and in wind speeds of hundreds of m s$^{-1}$ or more across the observable atmosphere. These differences persist throughout the parameter regime studied, and are dependent on surface gravity through the impact on photosphere pressure. These results imply that the internal evolution and atmospheric thermal structure and dynamics of hot and ultra-hot Jupiters are coupled. As a result, a joint approach including both evolutionary models and GCMs may be required to make robust predictions for the atmospheric circulation of hot and ultra-hot Jupiters.

For a satellite with an irregular shape, which is the common shape among asteroids, the well-known spin-orbit resonance problem could be changed to a spin-orbit coupling problem since a decoupled model does not accurately capture the dynamics of the system. In this paper, having provided a definition for close binary asteroid systems, we explore the structure of the phase space in a classical Hamiltonian model for spin-orbit coupling in a binary system. To map out the geography of resonances analytically and the cartography of resonances numerically, we reformulate a fourth-order gravitational potential function, in Poincare variables, via Stokes coefficients. For a binary system with a near-circular orbit, isolating the Hamiltonian near each resonance yields the pendulum model. Analysis of the results shows the geographical information, including the location and width of resonances, is modified due to the prominent role of the semi-major axis in the spin-orbit coupling model but not structurally altered. However, this resulted in modified Chirikov criterion to predict onset of large-scale chaos. For a binary system with arbitrary closed orbit, we thoroughly surf in the phase space via cartography of resonances created by fast Lyapunov indicator (FLI) maps. The numerical study confirms the analytical results, provides insight into the spin-orbit coupling, and shows some bifurcations in the secondary resonances which can occur due to material transfer. Also, we take the (65803) Didymos binary asteroid as a case to show analytical and numerical results.

We present the photometric analysis of three ultra-short period total eclipsing binaries in contact configuration, CRTS_J172718.0+431624, OGLE-BLG-ECL-000104, and OGLE-BLG-ECL-000012, mined from massive astronomical surveys. Using the available archival light curves (LCs) from Vista Variables in the Via Lactea (VVV), Optical Gravitational Lensing Experiment (OGLE), Zwicky Transient Facility (ZTF) and Catalina Sky Survey (CSS) in different passbands and new multi-band photometric observations with the 2.3 m Aristarchos telescope at Helmos Observatory, their relative physical parameters were derived. We explored the parameter space by using PIKAIA genetic algorithm optimizer. The best photometric solution and error budget estimation was adopted for each system through MCMC sampling of the global optimum. The approximate absolute parameters were derived for each contact system adopting an empirical mass-luminosity relation. All three systems have a mass ratio lower than 0.5. The exchange between primary and secondary depths of CRTS_J172718.0+431624 during the years 2016-2022 may be due to spot activity. In addition, we present a detailed analysis of the first well characterized shortest period contact eclipsing binary with total eclipses known so far (OGLE-BLG-ECL-000104). Thanks to VVV and OGLE LCs, new distances were derived for OGLE-BLG-ECL-000104 and OGLE-BLG-ECL-000012 using empirical period-luminosity relations. The origin and evolutionary status of all three ultra-short period contact binaries are thoroughly discussed in the context of the detached-binary formation channel.

Conditions in the outer protoplanetary disk during Solar System formation were thought to be favorable for the formation of amorphous water ice (AWI),a glassy phase of water ice. However, subsequent collisional processing could have shock crystallized any AWI present. Here we use the iSALE shock physics hydrocode to simulate impacts between large icy bodies at impact velocities relevant to these collisional environments, and then feed these results into a custom-built AWI crystallization script, to compute how much AWI crystallizes/survives these impact events. We find that impact speeds between icy bodies post-planet migration (i.e., between trans-Neptunian Objects or TNOs) are too slow to crystallize any meaningful fraction of AWI. During planet migration, however, the amount of AWI that crystallizes is highly stochastic: relatively little AWI crystallizes at lower impact velocities (less than ~2 km/s), yet most AWI present in the bodies (if equal sized) or impactor and impact site (if different sizes) crystallizes at higher impact velocities (greater than ~4 km/s). Given that suspected impact speeds during planet migration were ~2-4 km/s, this suggests that primordial AWI's ability to survive planet migration is highly stochastic. However, if proto-EKB objects and their fragments experienced multiple impact events, nearly all primordial AWI could have crystallized; such a highly collisional proto-EKB during planet migration is consistent with the lack of any unambiguous direct detection of AWI on any icy body. Ultimately, primordial AWI's survival to the present day depends sensitively on the proto-EKB's size-frequency distribution, which is currently poorly understood.

We study the detectability of gravitational-wave signals from sub-solar mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as $78.4 \%$ for an equal mass binary system of chirp mass $0.17 \, \textrm{M}_{\odot}$, in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for sub-solar mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of sub-solar mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for sub-solar mass binary neutron star signals.

Arguments are provided for the reality of the quantum vacuum fields. A polarization correlation experiment with two maximally entangled photons created by spontaneous parametric down-conversion is studied in the Weyl-Wigner formalism, that reproduces the quantum predictions. An interpretation is proposed in terms of stochastic processes assuming that the quantum vacuum fields are real. This proves that local realism is compatible with a violation of Bell inequalities, thus rebutting the claim that it has been refuted by experiments. Entanglement appears as a correlation between fluctuations of a signal field and vacuum fields. Key words; local realism, Bell inequalities, entangled photons, Weyl-Wigner, loopholes,vacuum fields

Current generation of detectors using noble gases in liquid phase for direct dark matter search and neutrino physics need large area photosensors. Silicon based photo-detectors are innovative light collecting devices and represent a successful technology in these research fields. %of direct dark matter search detectors based on liquified noble gases. The DarkSide collaboration started a dedicated development and customization of SiPM technology for its specific needs resulting in the design, production and assembly of large surface modules of 20$\times$20 cm$^{2}$ named Photo Detection Unit for the DarkSide-20k experiment. Production of a large number of such devices, as needed to cover about 20 m$^{2}$ of active surface inside the DarkSide-20k detector, requires a robust testing and validation process. In order to match this requirement a dedicated test facility for the photosensor test was designed and commissioned at INFN-Naples laboratory. The first commissioning test was successfully performed in 2021. Since then a number of testing campaigns were performed. Detailed description of the facility is reported as well as results of some tests.

Plasma current filamentation of an ultrarelativistic electron beam impinging on an overdense plasma is investigated, with emphasis on radiation-induced electron polarization. Particle-in-cell simulations provide the classification and in-depth analysis of three different regimes of the current filaments, namely, the normal filament, abnormal filament, and quenching regimes. We show that electron radiative polarization emerges during the instability along the azimuthal direction in the momentum space, which significantly varies across the regimes. We put forward an intuitive Hamiltonian model to trace the origin of the electron polarization dynamics. In particular, we discern the role of nonlinear transverse motion of plasma filaments, which induces asymmetry in radiative spin flips, yielding an accumulation of electron polarization. Our results break the conventional perception that quasi-symmetric fields are inefficient for generating radiative spin-polarized beams, suggesting the potential of electron polarization as a source of new information on laboratory and astrophysical plasma instabilities.

Motivated by events in which black holes can lose their environment due to tidal interactions in a binary system, we develop a waveform model in which the tidal deformability interpolates between a finite value (dressed black hole) at relatively low frequency and a zero value (naked black hole) at high frequency. We then apply this model to the example case of a black hole dressed with an ultralight scalar field and investigate the detectability of the tidal Love number with the Einstein Telescope. We show that the parameters of the tidal deformability model could be measured with high accuracy, providing a useful tool to understand dynamical environmental effects taking place during the inspiral of a binary system.

A crossover QCD phase transition in the early Universe, involving a formation scenario of stable strangeon nuggets is studied. The Polyakov-Nambu-Jona-Lasinio model is applied to calculate the thermodynamics of the QCD phase with u, d, s quarks, and the relativistic mean-field model describes the hadronic matter. The crossover phase transition from quarks to hadrons occurred at cosmic temperature of T~170 MeV, and those two phases are connected in a three-window model. Due to quark's non-perturbative coupling, quark clusters with net strangeness (i.e., strangeons) and then strangeon nuggets could form during the transition process. A distribution function of the nugget baryon number, A, is introduced to describe the nuggets' number density. All the strangeon nuggets with A>A_c are considered to be stable, where the critical number, A_c, is determined by both the weak and strong interactions. A non-relativistic equation of state is applied to calculate the thermodynamics of stable nuggets. The calculation shows that the thermodynamical contributions (pressure, entropy, etc.) of the stable strangeon nuggets are negligible. The resultant mass density of the strangeon nuggets survival from the early Universe is comparable to the dark matter, that indicates a possible explanation of the cold dark matter without introducing any exotic particles beyond the standard model.

Recently, the LHAASO collaboration has observed the gamma rays of energies up to ten TeV from the gamma-ray burst GRB221009A, which has stimulated the community of astronomy, particle physics and astrophysics to propose various possible interpretations. In this paper, we put forward a viable scenario that neutrinos are produced together with TeV photons in the gamma-ray burst and gradually decay into the axion-like particles, which are then converted into gamma rays in the galactic magnetic fields. In such a scenario, the tension between previous axion-like particle interpretations and the existing observational constraints on the relevant coupling constant and mass can be relaxed.

We study spin effects in the neutrino gravitational scattering by a supermassive black hole with a magnetized accretion disk having a finite thickness. We exactly describe the propagation of ultrarelativistic neutrinos on null geodesics and solve the spin precession equation along each neutrino trajectory. The interaction of neutrinos with the magnetic field is owing to the nonzero diagonal magnetic moment. Additionally, neutrinos interact with plasma of the accretion disk electroweakly within the Fermi approximation. These interactions are obtained to change the polarization of incoming neutrinos, which are left particles. The fluxes of scattered neutrinos, proportional to the survival probability of spin oscillations, are derived for various parameters of the system. In particular, we are focused on the matter influence on the outgoing neutrinos flux. The possibility to observe the predicted effects for astrophysical neutrinos is briefly discussed.

The standard model for galaxy bias is built in a Newtonian framework, and several attempts have been made in the past to put it in a relativistic framework. The focus of past works was, however, to use the same Newtonian formulation, but to provide its interpretation in a relativistic framework by either fixing a gauge condition or transforming to a local coordinate system. Here we demonstrate that these reverse-engineered approaches fail, because they do not respect the diffeomorphism symmetry in general relativity. Hence we need to take a different approach to the problem and develop a covariant formulation of galaxy bias model that is diffeomorphism compatible. We consider a simple toy model for galaxy bias and discuss the implication for measuring the primordial non-Gaussianity by using the flawed standard model for galaxy bias.

Globular clusters contain multiple stellar populations, with some previous generation of stars polluting the current stars with heavier elements. Understanding the history of globular clusters is helpful in understanding how galaxies merged and evolved and therefore constraining the site or sites of this historic pollution is a priority. The acceptable temperature and density conditions of these polluting sites depend on critical reaction rates. In this paper, three experimental studies helping to constrain astrophysically important reaction rates are briefly discussed.

The total energy and other bound state properties of the ground (bound) $1^{1}S$-state in the Ps$^{-}$ ion are determined to very high accuracy. Our best variational energy for the ground state in this ion equals $E$ = -0.262005070232980107770402018838 $a.u.$ For this three-body ion we have evaluated (to very high accuracy) the rates of two-, three-, four- and five-photon annihilation. We also discuss some problems which currently exist in accurate computations of the rate of one-photon annihilation $\Gamma_{1 \gamma}$. Highly accurate computations of a number of singular and quasi-singular bound state properties in the Ps$^{-}$ ion are also performed and discussed. By investigating the sources of annihilation $\gamma-$quanta in the universe we have arrived to the conclusion about the high-temperature limit in optics. This can be formulated by the following statement: due to the electromagnetic instability of the vacuum, it is impossible to see (directly) any object heated to temperatures above 350 - 400 $keV$. In reality, instead of such an object an observer will see only an intense flow of annihilation $\gamma-$quanta, electrons and positrons. This phenomenon can be called the annihilation shielding of overheated matter and it is of great interest in Galactic astrophysics.

We present a simple procedure to obtain universal bounds for quantities of cosmological interest, such as the number of $e$-folds during inflation, reheating, and radiation, as well as the reheating temperature. The main assumption is to represent each of the various epochs of evolution of the universe as being due to a single substance changing instantaneously into the next, describing a new era of evolution of the universe. This assumption, commonly used to obtain solutions of the Friedmann equations for simple cosmological models, is implemented here to find model-independent bounds on cosmological quantities of interest. In particular, we find that the bound $N_k\approx 56$ for $-\frac{1}{3} < \omega_{re} < \frac{1}{3}$ is very robust as an upper bound on the number of $e$-folds during inflation and also as a lower bound when $\omega_{re} > \frac{1}{3}$, where $\omega_{re}$ is the effective equation of state parameter during reheating. These are model-independent results that any single-field model of inflation should satisfy. As an example, we illustrate the two approaches with the basic $\alpha$ attractor model and show how they complement each other.

Charged-current neutrino processes such as $\nu_e + n \rightleftharpoons p + e^-$ and $\bar\nu_e + p \rightleftharpoons n + e^+$ destroy the flavor coherence among the weak-interaction states of a single neutrino and thus damp its flavor oscillation. In a dense neutrino gas such as that inside a core-collapse supernova or the black hole accretion disk formed in a compact binary merger, however, these "collision" processes can trigger large flavor conversion in cooperation with the strong neutrino-neutrino refraction. We show that there exist two types of collisional flavor instability in a homogeneous and isotropic neutrino gas which are identified by the dependence of their real frequencies on the neutrino density $n_\nu$. The instability transitions from one type to the other and exhibits a resonance-like behavior in the region where the net electron lepton number of the neutrino gas is negligible. In the transition region, the flavor instability grows exponentially at a rate $\propto n_\nu^{1/2}$. We find that the neutrino gas in the black hole accretion disk is susceptible to the collision-induced flavor conversion where the neutrino densities are the highest. As a result, large amounts of heavy-lepton flavor neutrinos can be produced through flavor conversion, which can have important ramifications in the subsequent evolution of the remnant.

We construct new classes of modified theories in which the matter sector couples with the Einstein tensor, namely we consider direct couplings of the latter to the energy-momentum tensor, and to the derivatives of its trace. We extract the general field equations and we apply them in a cosmological framework, obtaining the Friedmann equations, whose extra terms give rise to an effective dark energy sector. At the background level we show that we can successfully describe the usual thermal history of the universe, with the sequence of matter and dark-energy epochs, while the dark-energy equation-of-state parameter can lie in the phantom regime, tending progressively to -1 at present and future times. Furthermore, we confront the theory with Cosmic Chronometer data, showing that the agreement is very good. Finally, we perform a detailed investigation of scalar and tensor perturbations, and extracting an approximate evolution equation for the matter overdensity we show that the predicted behavior is in agreement with observations.

A large portion of astronomy's carbon footprint stems from fossil fuels supplying the power demand of astronomical observatories. Here, we explore various isolated low-carbon power system setups for the newly planned Atacama Large Aperture Submillimeter Telescope, and compare them to a business-as-usual diesel power generated system. Technologies included in the designed systems are photovoltaics, concentrated solar power, diesel generators, batteries, and hydrogen storage. We adapt the electricity system optimization model highRES to this case study and feed it with the telescope's projected energy demand, cost assumptions for the year 2030 and site-specific capacity factors. Our results show that the lowest-cost system with LCOEs of $116/MWh majorly uses photovoltaics paired with batteries and fuel cells running on imported green hydrogen. Some diesel generators run for backup. This solution would reduce the telescope's power-side carbon footprint by 95% compared to the business-as-usual case.