Papers for Tuesday, Oct 18 2022

The diffuse, unresolved sky provides most of the photons that the Hubble Space Telescope (HST) receives, yet remains poorly understood. HST Archival Legacy program SKYSURF aims to measure the 0.2-1.7 $\mu$m sky surface brightness (sky-SB) from over 200,000 HST images. We describe two sky-SB measurement algorithms designed for SKYSURF (the Percentile-clip and ProFound Median methods), that are able to recover the input sky-SB from simulated images to within 1% uncertainty, and present measurements estimated using each algorithm on the entire SKYSURF database. Comparing our sky-SB spectral energy distribution to measurements from the literature shows general agreements, but highlights that models of Zodiacal Light at HST wavelengths are likely incomplete. Our SKYSURF spectral energy distribution also reveals a dependence on Sun Angle, indicating non-isotropic scattering of Solar photons. Finally, we update Diffuse Light limits for F125W, F140W, and F160W based on the methods from Carleton et al. (2022). The Diffuse Light limits for both sky-SB measurements algorithms are in good agreement, ranging from 0.006 MJy sr$^{-1}$ (14 nW m$^{-2}$ sr$^{-1}$) to 0.015 MJy sr$^{-1}$ (32 nW m$^{-2}$ sr$^{-1}$). These estimates provide the most stringent all-sky constraints to date in this wavelength range. SKYSURF sky-SB measurements are made public on the official SKYSURF website and will be used to constrain Extragalactic Background Light in future papers.

The diversity of dynamical conditions among exoplanets is now well established. Yet, the relevance of orbital dynamical timescales to biological evolutionary timescales is poorly understood. Given that even minor orbital changes may place significant pressure on any organisms living on a planet, dynamical sculpting has important implications for the putative evolution of life. In this manuscript, we employ a Monte Carlo framework to investigate how a range of exoplanetary dynamical sculpting timescales affects timescales for biological evolution. We proceed with minimal assumptions for how dynamical sculpting proceeds and the emergence and persistence of life. We focus our investigation on M dwarf stars, the most common exoplanetary hosts in the Milky Way. We assign dynamical statuses, dependent on stellar age, to a suite of planetary systems, varying the rate of dynamical disruption within limits that are consistent with present-day planet demographics. We then simulate the observed yield of planets according to the completeness of NASA's Kepler and TESS missions, and investigate the properties of these samples. With this simplified approach, we find that systems hosting multiple transiting planets ought to have, on average, shorter dynamically-uninterrupted intervals than single-transiting systems. However, depending upon the rate of dynamical sculpting, planets orbiting older stars will exhibit the opposite trend. Even modest constraints on stellar age would help identify "older" stars for which this holds. The degree of these effects varies, dependent upon both the intrinsic dynamical demographics of exoplanets and whether we consider planets detected by NASA's Kepler or TESS missions.

We present the first high-resolution zoom-in simulation of a Milky-way-like halo extracted from the Aquarius Project in the Fuzzy Dark Matter (FDM) framework. We use the N-body code AX-GADGET, based on a particle oriented solution of the Schr\"{o}dinger-Poisson equations, able to detail the complexity of structure formation while keeping track of the quantum effects in FDM. The halo shows a cored density profile, with a core size of several kpc for a FDM mass of $m_{\chi}=2.5h \times 10^{-22}\ {\rm eV}/c^2$. A flattening is observed also in the velocity profile, representing a distinct feature of FDM dynamics. We provide a quantitative analysis of the impact of fuzziness on satellites in terms of abundance, mass, distance and velocity distribution functions and their evolution with redshift. Very interestingly, we show that all collapsed structures, despite showing a flat density profile at $z=0$, do not reach the solitonic ground-state at the time of formation: on the contrary, they asymptotically converge to it on a timescale that depends on their mass and formation history. This implies that current limits on FDM mass - obtained by applying simple scaling relations to observed galaxies - should be taken with extreme care, since single objects can significantly deviate from the expected asymptotic behavior during their evolution.

Relativistic jets originating from protomagnetar central engines can lead to long duration gamma-ray bursts (GRBs) and are considered potential sources of ultrahigh-energy cosmic rays and secondary neutrinos. We explore the propagation of such jets through a broad range of progenitors, from stars which have shed their envelopes to supergiants which have not. We use a semi-analytical spindown model for the strongly magnetised and rapidly rotating PNS to investigate the role of central engine properties such as the surface dipole field strength, initial rotation period, and jet opening angle on the interactions and dynamical evolution of the jet-cocoon system. With this model, we determine the properties of the relativistic jet, the mildly-relativistic cocoon, and the collimation shock in terms of system parameters such as the time-dependent jet luminosity, injection angle and density profile of the stellar medium. We also analyse the criteria for a successful jet breakout, the maximum energy that can be deposited into the cocoon by the relativistic jet, and structural stability of the magnetised outflow relative to local instabilities. Lastly, we compute the high-energy neutrino emission as these magnetised outflows burrow through their progenitors. Precursor neutrinos from successful GRB jets are unlikely to be detected by IceCube, which is consistent with the results of previous works. On the other hand, we find high-energy neutrinos may be produced for extended progenitors like blue and red supergiants, and we estimate the detectability of neutrinos with next-generation detectors such as IceCube-Gen2.

We present X-ray spectral analysis of XMM and Chandra observations in the 31.3 deg$^2$ Stripe-82X (S82X) field. Of the 6181 unique X-ray sources in this field, we select and analyze a sample of 2937 candidate active galactic nuclei (AGN) with solid redshifts and sufficient counts determined by simulations. Our results show an observed population with median values of spectral index $\Gamma=1.94_{-0.39}^{+0.31}$, column density log $N_H/\rm{cm}^{-2}=20.7_{-0.5}^{+1.2}$ ($21.6_{-0.9}^{+0.7}$ considering upper limits) and intrinsic, de-absorbed, 2-10 keV luminosity log $L_X/\rm{erg\,s}^{-1}=44.0_{-1.0}^{+0.7}$, in the redshift range 0-4. Correcting for observational biases, we derive the intrinsic, model-independent, fraction of AGN that are obscured ($22\leq\rm{log}\,N_H/\rm{cm}^{-2}<24$), finding a significant increase in the obscured AGN fraction with redshift and a decline with increasing luminosity. The average obscured AGN fraction is $57\pm4\%$ for log $L_X/\rm{erg\,s}^{-1}>43$. This work constrains the AGN obscuration and spectral shape of the still uncertain high-luminosity and high-redshift regimes (log $L_X/\rm{erg\,s}^{-1}>45.5$, $z>3$), where the obscured AGN fraction rises to $64\pm12\%$. The total, obscured, and unobscured X-ray luminosity functions (XLFs) are determined out to $z=4$. We report a luminosity and density evolution of the total XLF, with obscured AGN dominating at all luminosities at $z>2$ and unobscured sources prevailing at log $L_X/\rm{erg\,s}^{-1}>45$ at lower redshifts. Our results agree with the evolutionary models in which the bulk of AGN activity is triggered by gas-rich environments and in a downsizing scenario. Also, the black hole accretion density is found to evolve similarly to the star formation rate density, confirming the co-evolution scenario between AGN and host galaxy, but suggesting different time scales in their growing history.

Polar dust has been found to play an important role in the mid-infrared emission of nearby Seyfert nuclei. If and how often polar dust exists among the quasar population is unknown due to the lack of spatially-resolved observations. In this Letter, we report correlations between the prominence of AGN forbidden line emission (commonly associated with the narrow line region) and the dust mid-IR energy output among the archetypal Palomar-Green quasar sample and other bright type-1 AGNs drawn from the SDSS, Spitzer and WISE archives. The AGN mid-IR color differences traced by WISE W2 ($\sim4.6 \mu m$)$-$W3 ($\sim12 \mu m$) and W2 ($\sim4.6 \mu m$)$-$W4 ($\sim22 \mu m$), and near-IR to mid-IR SEDs constrained with 2MASS, WISE and Spitzer data have clear trends with the relative strength of the forbidden line regions traced by the optical \OIII and mid-IR \OIV emission lines. These observations indicate that, where the lines are strong, a large fraction of the AGN emission at $\lambda\gtrsim5 \mu$m comes from dust in the forbidden line regions. We find that the widely quoted universal AGN template is a result of averaging quasar SEDs with different levels of polar dust emission above the torus output and that the typical intrinsic IR SED of compact torus dust emission alone falls with increasing wavelength past 5 $\mu$m (in $\nu F_\nu$). In addition, the association of polar dust with the forbidden lines suggests an alternative to the receding torus hypothesis for the decrease in infrared output with increasing AGN luminosity.

We present estimates of line-of-sight distortion fields derived from the 95 GHz and 150 GHz data taken by BICEP2, BICEP3, and Keck Array up to the 2018 observing season, leading to cosmological constraints and a study of instrumental and astrophysical systematics. Cosmological constraints are derived from three of the distortion fields concerning gravitational lensing from large-scale structure, polarization rotation from magnetic fields or an axion-like field, and the screening effect of patchy reionization. We measure an amplitude of the lensing power spectrum $A_L^{\phi\phi}=0.95 \pm 0.20$. We constrain polarization rotation, expressed as the coupling constant of a Chern-Simons electromagnetic term $g_{a\gamma} \leq 2.6 \times 10^{-2}/H_I$, where $H_I$ is the inflationary Hubble parameter, and an amplitude of primordial magnetic fields smoothed over 1 Mpc $B_{1\text{Mpc}} \leq 6.6 \;\text{nG}$ at 95 GHz. We constrain the root mean square of optical-depth fluctuations in a simple "crinkly surface" model of patchy reionization, finding $A^\tau<0.19$ ($2\sigma$) for the coherence scale of $L_c=100$. We show that all of the distortion fields of the 95 GHz and 150 GHz polarization maps are consistent with simulations including lensed-$\Lambda$CDM, dust, and noise, with no evidence for instrumental systematics. In some cases, the EB and TB quadratic estimators presented here are more sensitive than our previous map-based null tests at identifying and rejecting spurious B-modes that might arise from instrumental effects. Finally, we verify that the standard deprojection filtering in the BICEP/Keck data processing is effective at removing temperature to polarization leakage.

We evolve two high-resolution general relativistic magnetohydrodynamic (GRMHD) simulations of advection-dominated accretion flows around non-spinning black holes (BHs), each over a duration $\sim 3\times 10^5\,GM_{\rm BH}/c^3$. One model captures the evolution of a weakly magnetized (SANE) disk and the other a magnetically arrested disk (MAD). Magnetic flux eruptions in the MAD model push out gas from the disk and launch strong winds with outflow efficiencies at times reaching $10\%$ of the incoming accretion power. Despite the substantial power in these winds, average mass outflow rates remain small out to a radius $\sim100\,GM_{\rm BH}/c^2$, only reaching $\sim 60-80\%$ of the horizon accretion rate. The average outward angular momentum transport is primarily radial in both modes of accretion, but with a clear distinction: magnetic flux eruption-driven disk winds cause a strong vertical flow of angular momentum in the MAD model, while for the SANE model, the magnetorotational instability (MRI) moves angular momentum mostly equatorially through the disk. Further, we find that the MAD state is highly transitory and non-axisymmetric, with the accretion mode often changing to a SANE-like state following an eruption before reattaining magnetic flux saturation with time. The Reynolds stress changes direction during such transitions, with the MAD (SANE) state showing an inward (outward) stress, possibly pointing to intermittent MRI-driven accretion in MADs. Pinning down the nature of flux eruptions using next-generation telescopes will be crucial in understanding the flow of mass, magnetic flux and angular momentum in sub-Eddington accreting BHs like M87$^*$ and Sagittarius A$^*$.

The PhotoDissociation Region Toolbox provides comprehensive, easy-to-use, public software tools and models that enable an understanding of the interaction of the light of young, luminous, massive stars with the gas and dust in the Milky Way and in other galaxies. It consists of an open-source Python toolkit and photodissociation region models for analysis of infrared and millimeter/submillimeter line and continuum observations obtained by ground-based and sub-orbital telescopes, and astrophysics space missions. Photodissociation regions (PDRs) include all of the neutral gas in the ISM where far-ultraviolet photons dominate the chemistry and/or heating. In regions of massive star formation, PDRs are created at the boundaries between the H II regions and neutral molecular cloud, as photons with energies 6 eV $ < h \nu < $ 13.6 eV photodissociate molecules and photoionize metals. The gas is heated by photo-electrons from small grains and large molecules and cools mostly through far-infrared fine-structure lines like [O I] and [C II]. The models are created from state-of-the art PDR codes that includes molecular freeze-out; recent collision, chemical, and photo rates; new chemical pathways, such as for oxygen chemistry; and allow for both clumpy and uniform media. The models predict the emergent intensities of many spectral lines and FIR continuum. The tools find the best-fit models to the observations and provide insights into the physical conditions and chemical makeup of the gas and dust. The PDR Toolbox enables novel analysis of data from telescopes such as ISO, Spitzer, Herschel, STO, SOFIA, SWAS, APEX, ALMA, and JWST.

The nonresonant streaming instability (Bell instability) plays a pivotal role in the acceleration and confinement of cosmic rays (CRs); yet, the exact mechanism responsible for its saturation and the magnitude of the final amplified magnetic field have not been assessed from first-principles. Using a survey of hybrid simulations (with kinetic ions and fluid electrons), we study the evolution of the Bell instability as a function of the parameters of the CR population. We find that, at saturation, the magnetic pressure in the amplified field is comparable with the initial CR momentum flux (i.e., the anisotropic CR pressure). These results provide a predictive prescription for the total magnetic field amplification expected in the many astrophysical environments where the Bell instability is important.

The cross-correlation function and template matching techniques have dominated the world of precision radial velocities for many years. Recently, a new technique, named line-by-line, has been developed as an outlier resistant way to efficiently extract radial velocity content from high resolution spectra. We apply this new method to archival HARPS and CARMENES datasets of the K2-18 system. After reprocessing the HARPS dataset with the line-by-line framework, we are able to replicate the findings of previous studies. Furthermore, by splitting the full wavelength range into sub-domains, we were able to identify a systematic chromatic correlation of the radial velocities in the reprocessed CARMENES dataset. After post-processing the radial velocities to remove this correlation, as well as rejecting some outlier nights, we robustly uncover the signal of both K2-18b and K2-18c, with masses that agree with those found from our analysis of the HARPS dataset. We then combine both the HARPS and CARMENES velocities to refine the parameters of both planets, notably resulting in a revised mass and period for K2-18c of $6.99^{+0.96}_{-0.99}$M$_{Earth}$ and $9.2072\pm0.0065$d, respectively. Our work thoroughly demonstrates the power of the line-by-line technique for the extraction of precision radial velocity information.

The standard perturbation theory (SPT) approach to gravitational clustering is based on a fluid approximation of the underlying Vlasov-Poisson dynamics, taking only the zeroth and first cumulant of the phase-space distribution function into account (density and velocity fields). This assumption breaks down when dark matter particle orbits cross and leads to well-known problems, e.g. an anomalously large backreaction of small-scale modes onto larger scales that compromises predictivity. We extend SPT by incorporating second and higher cumulants generated by orbit crossing. For collisionless matter, their equations of motion are completely fixed by the Vlasov-Poisson system, and thus we refer to this approach as Vlasov Perturbation Theory (VPT). Even cumulants develop a background value, and they enter the hierarchy of coupled equations for the fluctuations. The background values are in turn sourced by power spectra of the fluctuations. The latter can be brought into a form that is formally analogous to SPT, but with an extended set of variables and linear as well as non-linear terms, that we derive explicitly. In this paper, we focus on linear solutions, which are far richer than in SPT, showing that modes that cross the dispersion scale set by the second cumulant are highly suppressed. We derive stability conditions on the background values of even cumulants from the requirement that exponential instabilities be absent. We also compute the expected magnitude of averaged higher cumulants for various halo models and show that they satisfy the stability conditions. Finally, we derive self-consistent solutions of perturbations and background values for a scaling universe and study the convergence of the cumulant expansion. The VPT framework provides a conceptually straightforward and deterministic extension of SPT that accounts for the decoupling of small-scale modes.

We present non-linear solutions of Vlasov Perturbation Theory (VPT), describing gravitational clustering of collisionless dark matter with dispersion and higher cumulants induced by orbit crossing. We show that VPT can be cast into a form that is formally analogous to standard perturbation theory (SPT), but including additional perturbation variables, non-linear interactions, and a more complex propagation. VPT non-linear kernels have a crucial decoupling property: for fixed total momentum, the kernels becomes strongly suppressed when any of the individual momenta cross the dispersion scale into the non-linear regime. This screening of UV modes allows us to compute non-linear corrections to power spectra even for cosmologies with very blue power-law input spectra, for which SPT diverges. We compare predictions for the density and velocity divergence power spectra as well as the bispectrum at one-loop order to N-body results in a scaling universe with spectral indices $-1\leq n_s\leq +2$. We find a good agreement up to the non-linear scale for all cases, with a reach that increases with the spectral index $n_s$. We discuss the generation of vorticity as well as vector and tensor modes of the velocity dispersion, showing that neglecting vorticity when including dispersion would lead to a violation of momentum conservation. We verify momentum conservation when including vorticity, and compute the vorticity power spectrum at two-loop order, necessary to recover the correct large-scale limit with slope $n_w=2$. Comparing to our N-body measurements confirms the cross-over from $k^4$ to $k^2$ scaling on large scales. Our results provide a proof-of-principle that perturbative techniques for dark matter clustering can be systematically improved based on the known underlying collisionless dynamics.

We compare stellar mass surface density, metallicity, age, and line-of-sight velocity dispersion profiles in massive ($M_*\geq10^{10.5}\,\mathrm{M_\odot}$) present-day early-type galaxies (ETGs) from the MaNGA survey with simulated galaxies from the TNG100 simulation of the IllustrisTNG suite. We find an excellent agreement between the stellar mass surface density profiles of MaNGA and TNG100 ETGs, both in shape and normalisation. Moreover, TNG100 reproduces the shapes of the profiles of stellar metallicity and age, as well as the normalisation of velocity dispersion distributions of MaNGA ETGs. We generally also find good agreement when comparing the stellar profiles of central and satellite galaxies between MaNGA and TNG100. An exception is the velocity dispersion profiles of very massive ($M_*\gtrsim10^{11.5}\,\mathrm{M_\odot}$) central galaxies, which, on average, are significantly higher in TNG100 than in MaNGA ($\approx50\,\mathrm{km\,s^{-1}}$). We study the radial profiles of $\mathit{in}$-$\mathit{situ}$ and $\mathit{ex}$-$\mathit{situ}$ stars in TNG100 and discuss the extent to which each population contributes to the observed MaNGA profiles. Our analysis lends significant support to the idea that high-mass ($M_*\gtrsim10^{11}\,\mathrm{M_\odot}$) ETGs in the present-day Universe are the result of a merger-driven evolution marked by major mergers that tend to homogenise the stellar populations of the progenitors in the merger remnant.

Recently, extra motivation has been given to the investigations of an unresolved problem of millisecond pulsars (MSPs) produced by the recycling process, as an apparent role of the accretion-induced collapse (AIC) in white dwarfs (WDs) was suggested to this concern. I have found that the distribution of the orbital periods of binary MSPs in the Galactic disk closely follows an exponential distribution. I have also determined the best-fit mean value of Nobs by fitting our data with an exponential distribution for the MSP population. As a result, it can be stated that reaching the Chandrasekhar limit may cause an explosion of a massive WD as a Type Ia supernova (in the case of a CO WD) or an ignition of a ONeMg WD, and possibly merging in some CO WDs, all resulting in peculiar MSP systems. A possible formation scenario, where the system has a circular orbit during this evolutionary stage, is discussed.

Between 2014 December 31 and 2015 March 17, the OSIRIS cameras on Rosetta documented the growth of a 140m wide and 0.5m deep depression in the Hapi region on Comet 67P/Churyumov-Gerasimenko. This shallow pit is one of several that later formed elsewhere on the comet, all in smooth terrain that primarily is the result of airfall of coma particles. We have compiled observations of this region in Hapi by the microwave instrument MIRO on Rosetta, acquired during October and November 2014. We use thermophysical and radiative transfer models in order to reproduce the MIRO observations. This allows us to place constraints on the thermal inertia, diffusivity, chemical composition, stratification, extinction coefficients, and scattering properties of the surface material, and how they evolved during the months prior to pit formation. The results are placed in context through long-term comet nucleus evolution modelling. We propose that: 1) MIRO observes signatures that are consistent with a solid-state greenhouse effect in airfall material; 2) CO2 ice is sufficiently close to the surface to have a measurable effect on MIRO antenna temperatures, and likely is responsible for the pit formation in Hapi observed by OSIRIS; 3) the pressure at the CO2 sublimation front is sufficiently strong to expel dust and water ice outwards, and to compress comet material inwards, thereby causing the near-surface compaction observed by CONSERT, SESAME, and groundbased radar, manifested as the "consolidated terrain" texture observed by OSIRIS.

We report the results of James Webb Space Telescope/NIRCam observations of 19 (sub)millimeter (submm/mm) sources detected by the Atacama Large Millimeter Array (ALMA). The accurate ALMA positions allowed unambiguous identifications of their NIRCam counterparts. Taking gravitational lensing into account, these represent 16 distinct galaxies in three fields and constitute the largest sample of its kind to date. The counterparts' spectral energy distributions from rest-frame ultraviolet to near infrared provide photometric redshifts ($1<z<4.5$) and stellar masses ($M_*>10^{10.5}$ Msol), which are similar to sub-millimeter galaxy (SMG) hosts studied previously. However, our sample is fainter in submm/mm than the classic SMG samples are, and our sources exhibit a wider range of properties. They have dust-embedded star-formation rates as low as 10 Msol yr$^{-1}$, and the sources populate both the star-forming main sequence and the quiescent categories. The deep NIRCam data allow us to study the rest-frame near-IR morphologies. Excluding two multiply imaged systems and one quasar, the majority of the remaining sources are disk-like and show either little or no disturbance. This suggests that secular growth is a potential route for the assembly of high-mass disk galaxies. While a few hosts have large disks, the majority have small disks (median half-mass radius of 1.6 kpc). At this time, it is unclear whether this is due to the prevalence of small disks at these redshifts or some unknown selection effects of deep ALMA observations. A larger sample of ALMA sources with NIRCam observations will be able to address this question.

We present the resolved properties of the $z=2.82$ Hyper Luminous Infrared Galaxy (HyLIRG) HS170850.1, the brightest 850$\mu$m source found in the SCUBA-2 followup to the Keck Baryonic Structure Survey fields (S$_{\rm 850 \mu m}=$19.5 mJy), and amongst the most luminous starbursts known at any redshift. Using the IRAM-NOEMA interferometer in the highest resolution A-configuration, we resolve the source into two components separated by $\sim$8 kpc, visible as blue shifted and red shifted $^{12}$CO(5-4) lines, exhibiting the expected kinematic properties of a major merger between two gas-rich galaxies. The combined merger system is traced over 2.3$''$ or 18.3 kpc. Each component of the merger shows ordered gas motions suggestive of a massive, turbulent disk. We measure the masses of the blue and red disks as (1.5 $\pm$ 0.2) $\times10^{11}$ M$_\odot$ and (0.71 $\pm$ 0.22) $\times10^{11}$ M$_\odot$ respectively. The more massive disk component shows broad wings in the CO line, offset by $\sim$3 kpc from the disk centroid along the major axis, and extending to velocities $\sim\pm$1000 km $ \rm{s^{-1}}$ from systemic velocity. We interpret this as either a possible bipolar outflowing component, or more likely a warping or tidal structure in the CO disk. Comparing the properties of HS170850.1 to other submillimeter detected galaxies with comparably bright 850$\mu$m luminosities suggests that ongoing gas-rich mergers, or at least a clustered/group environment lead to these most extreme starburst phases.

Aims: Our aim is to provide an observational view of the old Disk structure of the Milky Way galaxy using the distribution of red clump stars. The spiral arms, warp structure, and other asymmetries present in the Disk are re-visited using a systematic study of red clump star counts over the disk of the Galaxy. Methods: We developed a method to systematically extract the red clump stars from 2MASS ($J-K_s, ~J$) colour-magnitude diagram of $1^\circ \times 1^\circ$ bins in $\ell \times b$ covering the range $40^\circ \le \ell \le 320^\circ$ and $-10^\circ \le b \le 10^\circ$. 2MASS data continues to be important since it is able to identify and trace the red clump stars to much farther distances than any optical survey of the Disk. The foreground star contamination in the selected sample is removed by utilising the accurate astrometric data from Gaia EDR3 and astrophysical parameters from Gaia DR3. Results: We have generated a face-on-view (XY-plane) of the Galaxy depicting the density distribution and count ratio above and below the Galactic plane. The resulting overdensity of red clump stars traces the continuous morphology of the Outer arm from the second to the third Galactic quadrant. This is the first study to map the Outer arms across the disk using red clump stars. Through this study, we are able to trace the Outer arm well into the 3rd Galactic quadrant for the first time. Apart from the spiral structures, we also see a wave-like asymmetry above and below the Galactic plane with respect to longitudes, indicating a warp structure. The warp structure is studied systematically by tracing the ratio of red clump stars above and below the Galactic plane. We provide the first direct observational evidence of the asymmetry in the Outer spiral arms confirming that the spiral arms traced by the older population are also warped similar to the Disk.

We implement a multi-group and discrete-ordinate neutrino transport model in spherical symmetry which allows to simulate collective neutrino oscillations by including realistic collisional rates in a self-consistent way. We utilize this innovative model, based on strategic parameter rescaling, to study a recently proposed collisional flavor instability caused by the asymmetry of emission and absorption rates between $\nu_e$ and $\bar\nu_e$ for four different static backgrounds taken from different stages in a core-collapse supernova simulation. Our results confirm that collisional instabilities generally exist around the neutrinosphere during the SN accretion and post-accretion phase, as suggested by [arXiv:2104.11369]. However, the growth and transport of flavor instabilities can only be fully captured by models with global simulations as done in this work. With minimal ingredient to trigger collisional instabilities, we find that the flavor oscillations and transport mainly affect (anti)neutrinos of heavy lepton flavors around their decoupling sphere, which then leave imprints on their energy spectra in the free-streaming regime. For electron (anti)neutrinos, their properties remain nearly intact. We also explore various effects due to the decoherence from neutrino-nucleon scattering, artificially enhanced decoherence from emission and absorption, neutrino vacuum mixing, and inhomogeneous matter profile, and discuss the implication of our work.

Gaia HyperVelocity Star (HVS) kinematic observations favor a local escape velocity of ~700 km/s, nearly forty-percent greater than conventional estimates. Combining HVS and dwarf galaxy satellite data reveal the global escape velocity profile for the Galaxy smoothly traces an unbroken Keplerian decline from the central bar to the most remote satellite galaxy. We reveal a robust upper bound in baryonic mass discrepancy (maximal relative accelerations) linked to the virial theorem and obtain a fundamental and universal mass discrepancy-acceleration relation for virialized compact cosmic objects.

The rich absorption spectra of Sun-like stars are enticing probes for variations in the fine-structure constant, $\alpha$, which gauges the strength of electromagnetism. While individual line wavelengths are sensitive to $\alpha$, they are also sensitive to physical processes in the stellar atmospheres, which has precluded their use so far. Here we demonstrate a new, differential approach using solar twins: velocity separations between close pairs of transitions are compared across stars with very similar physical properties, strongly suppressing astrophysical and instrumental systematic errors. We utilise 423 archival exposures of 18 solar twins from the High-Accuracy Radial velocity Planetary Searcher (HARPS), in which calibration errors can be reduced to $\lesssim$3 m/s. For stars with $\approx$10 high signal-to-noise ratio spectra ($\ge$200 per pixel), velocity separations between pairs are measured with $\approx$10 m/s statistical precision. A companion paper assesses a range of systematic error sources using 130 stars, with a greater range of stellar parameters, providing accurate corrections for astrophysical effects and a residual, intrinsic star-to-star scatter of 0-13 m/s. Within these uncertainties, we find no evidence for velocity separation differences in 17 transition pairs between solar twins. In a second companion paper, this is found to limit local ($\lesssim$50 pc) variations in $\alpha$ to $\approx$50 parts per billion, $\sim$2 orders of magnitude less than other Galactic constraints.

Sun-like stars are a new probe of variations in the fine-structure constant, $\alpha$, via the solar twins approach: velocity separations of close pairs of absorption lines are compared between stars with very similar stellar parameters, i.e. effective temperature, metallicity and surface gravity within 100K, 0.1 dex and 0.2 dex of the Sun's values. Here we assess possible systematic errors in this approach by analysing $\gtrsim$10,000 archival exposures from the High-Accuracy Radial velocity Planetary Searcher (HARPS) of 130 stars covering a much broader range of stellar parameters. We find that each transition pair's separation shows broad, low-order variations with stellar parameters which can be accurately modelled, leaving only a small residual, intrinsic star-to-star scatter of 0-33 m/s (average $\approx$7 m/s, $\approx$10$^{-4}$\r{A} at 5000\r{A}). This limits the precision available from a single pair in one star. We consider potential systematic errors from a range of instrumental and astrophysical sources (e.g. wavelength calibration, charge transfer inefficiency, stellar magnetic activity, line blending) and conclude that variations in elemental abundances, isotope ratios and stellar rotational velocities may explain this star-to-star scatter. Finally, we find that the solar twins approach can be extended to solar analogues - within 300K, 0.3 dex and 0.4 dex of the Sun's parameters - without significant additional systematic errors, allowing a much larger number of stars to be used as probes of variation in $\alpha$, including at much larger distances.

Magnetic reconnection in the relativistic regime has been proposed as an important process for the efficient production of nonthermal particles and high-energy emissions. Using fully kinetic particle-in-cell simulations, we investigate how guide-field strength and domain size affect characteristic spectral features and acceleration processes. We study two stages of acceleration: energization up until the injection energy $\gamma_{\rm inj}$ and further acceleration that generates a power-law spectrum. Stronger guide fields increase the power-law index and $\gamma_{\rm inj}$, which suppresses acceleration efficiency. These quantities seemingly converge with increasing domain size, suggesting that our findings can be extended to large-scale systems. We find that three distinct mechanisms contribute to acceleration during injection: particle streaming along the parallel electric field, Fermi reflection, and the pickup process. Fermi and pickup processes, related to the electric field perpendicular to the magnetic field, govern the injection for weak guide fields and larger domains. Meanwhile, parallel electric fields are important for injection in the strong guide field regime. In the post-injection stage, we find that perpendicular electric fields dominate particle acceleration in the weak guide field regime, whereas parallel electric fields control acceleration for strong guide fields. These findings will help explain the nonthermal acceleration and emissions in high-energy astrophysics, including black hole jets and pulsar wind nebulae.

The use of the Audio Spectrogram Transformer (AST) model for gravitational-wave data analysis is investigated. The AST machine-learning model is a convolution-free classifier that captures long-range global dependencies through a purely attention-based mechanism. In this paper a model is applied to a simulated dataset of inspiral gravitational wave signals from binary neutron star coalescences, built from five distinct, cold equations of state (EOS) of nuclear matter. From the analysis of the mass dependence of the tidal deformability parameter for each EOS class it is shown that the AST model achieves a promising performance in correctly classifying the EOS purely from the gravitational wave signals, especially when the component masses of the binary system are in the range $[1,1.5]M_{\odot}$. Furthermore, the generalization ability of the model is investigated by using gravitational-wave signals from a new EOS not used during the training of the model, achieving fairly satisfactory results. Overall, the results, obtained using the simplified setup of noise-free waveforms, show that the AST model, once trained, might allow for the instantaneous inference of the cold nuclear matter EOS directly from the inspiral gravitational-wave signals produced in binary neutron star coalescences.

We constrain the intrinsic Eddington ratio (\lamEdd ) distribution function for local AGN in bins of low and high obscuration (log NH <= 22 and 22 < log NH < 25), using the Swift-BAT 70-month/BASS DR2 survey. We interpret the fraction of obscured AGN in terms of circum-nuclear geometry and temporal evolution. Specifically, at low Eddington ratios (log lamEdd < -2), obscured AGN outnumber unobscured ones by a factor of ~4, reflecting the covering factor of the circum-nuclear material (0.8, or a torus opening angle of ~ 34 degrees). At high Eddington ratios (\log lamEdd > -1), the trend is reversed, with < 30% of AGN having log NH > 22, which we suggest is mainly due to the small fraction of time spent in a highly obscured state. Considering the Eddington ratio distribution function of narrow-line and broad-line AGN from our prior work, we see a qualitatively similar picture. To disentangle temporal and geometric effects at high lamEdd, we explore plausible clearing scenarios such that the time-weighted covering factors agree with the observed population ratio. We find that the low fraction of obscured AGN at high lamEdd is primarily due to the fact that the covering factor drops very rapidly, with more than half the time is spent with < 10% covering factor. We also find that nearly all obscured AGN at high-lamEdd exhibit some broad-lines. We suggest that this is because the height of the depleted torus falls below the height of the broad-line region, making the latter visible from all lines of sight.

We report ALMA observation of a $z\gtrsim10$ galaxy candidate (GHZ1) discovered from the GLASS-JWST Early Release Science Program. Our ALMA program aims to detect the [OIII] emission line at the rest-frame 3393.0062 GHz ($88.36\mu$m) and far-IR continuum emission with the spectral window setup seamlessly covering 26.125 GHz frequency range ($10.10<z<11.14$). A total of 7 hours of on-source integration was employed, using four frequency settings to cover the full range (1.7 hours per setting), with 0.7" angular resolution. No line or continuum is clearly detected, with a 5$\sigma$ limit to the line emission of 0.4 mJy beam$^{-1}$ at 150 km s$^{-1}$ channel$^{-1}$, and to the continuum emission of 30$\mu$Jy beam$^{-1}$. We report marginal spectral and continuum features ($4.4\sigma$ and $2.6\sigma$ peak signal-to-noise ratio, respectively), within 0.17" of the JWST position of GHZ1. These features with an implication of a spectroscopic redshift of $z=10.38$, need to be verified with further observation. Assuming that the best photometric redshift estimate ($z=10.60^{+0.52}_{-0.60}$) is correct, the broadband galaxy spectral energy distribution model for the $3\sigma$ upper limit of the continuum flux from GHZ1 suggests that GHZ1 might have low dust temperature ($T_d<30$K). The $5\sigma$ upper limit of [OIII] line luminosity and the star formation rate estimate with large error, are consistent with the properties of low metallicity galaxies and those of local starburst galaxies. We also report clear detections of continuum emission from six galaxies in the field with the JWST counterparts.

Ultracompact cataclysmic variables (CVs) of the AM CVn type are deemed to be important verification sources for the future space gravitational wave detectors such as the Laser Interferometer Space Antenna (LISA). We model the present-day Galactic population of AM CVn stars with He-star donors. Such a population has long expected to exist, though only a couple of candidates are known. We applied the hybrid method of binary population synthesis (BPS) which combines a simulation of the population of immediate precursors of AM CVn stars by a fast BPS code with subsequent tracking of their evolution by a full evolutionary code. The model predicts that the present birthrate of He-donor AM CVn stars in the Galaxy is $4.6\times 10^{-4}$ per yr and the Galaxy may harbour about 112000 objects of this class which have orbital periods less than 42-43 min. The foreground confusion limit and instrumental noise of LISA prevent the discovery of longer periods systems in gravitational waves. We find that about 500 He-star AM CVns may be detected by LISA with signal-to-noise ratio (S/N)>5 during a 4 yr mission. Within 1 Kpc from the Sun, there may exist up to 130 He-star AM CVns with the periods in the same range, which may serve as verification binaries, if detected in the electromagnetic spectrum. In the Milky Way, there are also about 14800 immediate precursors of AM CVn stars. They are detached systems with a stripped low-mass He-star and a white dwarf companion, out of which about 75 may potentially be observed by LISA during its mission.

Cosmic string networks form during cosmological phase transitions as a consequence of the Kibble mechanism. The evolution of the simplest networks is accurately described by the canonical Velocity Dependent One-Scale (VOS) model. However, numerical simulations have demonstrated the existence of significant quantities of short-wavelength propagation modes on the strings, known as wiggles, which motivated the recent development of a wiggly string extension of the VOS. Here we summarize recent progress in the physical interpretation of this model through a systematic study of the allowed asymptotic scaling solutions of the model. The modeling mainly relies on three mechanisms: the universe's expansion rate, energy transfer mechanisms (e.g., the production of loops and wiggles), and the choice of the scale in which wiggles are coarse-grained. We consider the various limits in which each mechanism dominates and compare the scaling solutions for each case, in order to gain insight into the role of each mechanism in the overall behavior of the network. Our results show that there are three scaling regimes for the wiggliness, consisting of the well-known Nambu-Goto solution, and non-trivial regimes where the amount of wiggliness can grow as the network evolves or, for specific expansion rates, become a constant. We also demonstrate that full scaling of the network is more likely in the matter era than in the radiation epoch, in agreement with numerical simulations.

The precursors of short and long Gamma Ray Bursts (SGRBs and LGRBs) can serve as probes of their progenitors, as well as shedding light on the physical processes of mergers or core-collapse supernovae. Some models predict the possible existence of Quasi-Periodically Oscillations (QPO) in the precursors of SGRBs. Although many previous studies have performed QPO search in the main emission of SGRBs and LGRBs, so far there was no systematic QPO search in their precursors. In this work, we perform a detailed QPO search in the precursors of SGRBs and LGRBs detected by Fermi/GBM from 2008 to 2019 using the power density spectrum (PDS) in frequency domain and Gaussian processes (GP) in time domain. We do not find any convinced QPO signal with significance above 3 $\sigma$, possibly due to the low fluxes of precursors. Finally, the PDS continuum properties of both the precursors and main emissions are also studied for the first time, and no significant difference is found in the distributions of the PDS slope for precursors and main emissions in both SGRBs and LGRBs.

We report the identification of a set of old super metal-rich dwarf stars with orbits of low eccentricity that reach a maximum height from the Galactic plane between ~0.5-1.5 kpc. We discuss their properties to understand their origins. We use data from the internal data release 6 of the Gaia-ESO Survey. We selected stars observed at high resolution with abundances of 21 species of 18 individual elements. We apply hierarchical clustering to group the stars with similar chemical abundances within the complete chemical abundance space. According to their chemical properties, this set of super metal-rich stars can be arranged into five subgroups. Four seem to follow a chemical enrichment flow, where nearly all abundances increase in lockstep with Fe. The fifth subgroup shows different chemical characteristics. All subgroups have the following features: median ages of the order of 7-9 Gyr, Solar or sub-Solar [Mg/Fe] ratios, maximum height between 0.5-1.5 kpc, low eccentricities, and a detachment from the expected metallicity gradient with guiding radius. The high metallicity of our stars is incompatible with a formation in the Solar neighbourhood. Their dynamic properties agree with theoretical expectations that these stars travelled from the inner Galaxy due to blurring and, most importantly, to churning. We suggest that most of this population's stars originated in the Milky Way's inner regions (inner disc and/or the bulge) and later migrated to the Solar neighbourhood. The region from where the stars originated had a complex chemical enrichment history, with contributions from supernovae types Ia and II and possibly asymptotic giant branch stars.

Transmission spectroscopy is a prime technique to study the chemical composition and structure of exoplanetary atmospheres. Strong excess absorption signals have been detected in the optical Na I D1, 2 Fraunhofer lines during transits of hot Jupiters, which are attributed to the planetary atmospheres and allow us to constrain its structure. We study the atmosphere of WASP-7 b by means of high-resolution transit spectroscopy in the sodium lines. We analyze a spectral transit time-series of 89 high-resolution spectra of the hot Jupiter WASP-7 b that was observed using the Ultraviolet and Visual Echelle Spectrograph (UVES). We use the telluric lines for an accurate alignment of the spectra and carry out a telluric correction with molecfit. Stellar magnetic activity is monitored by investigating chromospheric lines such as the Ca II H and K and hydrogen H$\alpha$ lines. Finally, we obtain transmission spectra and light curves for various lines. The star shows no identifiable flares and, if any, marginal changes in activity during our observing run. The sodium transmission spectra and corresponding light curves clearly show signs of the Rossiter-McLaughlin effect (RM) and the stellar center-to-limb variation (CLV) that we model using synthetic spectra. A statistically significant, narrow absorption feature with a line contrast of 0.50$\pm$0.06% (at $\sim 8.3\sigma$ level) and a full width at half maximum (FWHM) of 0.13$\pm$0.02 A is detected at the location of the Na I D$_2$ line. For the Na I D$_1$ line signal, we derive a line contrast of 0.13$\pm$0.04% (at $\sim 3.2\sigma$ level), which we consider a tentative detection. In addition, we provide upper limits for absorption by the hydrogen Balmer lines (H$\alpha$, H$\beta$, and H$\gamma$), K I $\lambda$7699 A, Ca II H and K, and infra-red triplet (IRT) lines.

In black hole X-ray binaries, a non-thermal high-energy component is sometimes detected at energies above 200 keV. The origin of this component is debated and distinct spectral modelizations can lead to different interpretations. High-energy polarimetry measurements with INTEGRAL allow new diagnostics on the physics responsible for the MeV spectral component. In this work, we aim to investigate the high-energy behavior of three bright sources discovered by the MAXI: MAXI J1535-571, MAXI J1820+070 and MAXI J1348-630. We take advantage of their brightness to investigate their soft gamma-ray (0.1-2 MeV) properties with INTEGRAL. We use both spectral and polarimetric approaches to probe their high-energy emission with the aim to bring new constraints on the ~ MeV emission. We first study the spectral characteristics of the sources in the 3-2000 keV using JEM-X, IBIS and SPI with a semi-phenomenological description of the data. We then use IBIS as a Compton telescope in order to evaluate the polarization properties of the sources above 300 keV. A high-energy component is detected during the HIMS and SIMS of MAXI J1535-571, the LHS of MAXI J1820+070 and the LHS of MAXI J1348-630. The components detected in MAXI J1820+070 and MAXI J1348-630 are polarized with a polarization fraction of 26 +/- 9{\deg} and > 56 % in the 300-1000 keV, respectively. With no polarization information for MAXI J1535-571, the component detected could either come from the jets or the corona. In the case of MAXI J1820+070, the extrapolation of the synchrotron spectrum measured in the infrared indicates that the component is likely due to a non-thermal distribution of electrons from a hybrid corona. For MAXI J1348-630, the high fraction of polarization points towards a jets origin, however, we cannot formally conclude without any infrared data giving information on the optically thin part of the synchrotron spectrum.

This work introduces new low-temperature gas opacities, in the range 3.2 <= log(T/K) <= 4.5, computed with the AESOPUS code under the assumption of thermodynamic equilibrium (Marigo &_Aringer_2009). In comparison to the previous version AESOPUS 1.0, we updated and expanded molecular absorption to include 80 species, mostly using the recommended line lists currently available from the ExoMol and HITRAN databases. Furthermore, in light of a recent study, we revised the H- photodetachment cross section, added the free-free absorption of other negative ions of atoms and molecules, and updated the collision-induced absorption due to H2/H2, H2/H, H2/He, and H/He pairs. Using the new input physics, we computed tables of Rosseland mean opacities for several scaled-solar chemical compositions, including Magg et al. (2022)'s most recent one, as well as alpha-enhanced mixtures. The differences in opacity between the new AESOPUS 2.0 and the original AESOPUS 1.0 versions, as well as other sets of calculations, are discussed. The new opacities are released to the community via a dedicated web-page that includes both pre-computed tables for widely used chemical compositions, and a web-interface for calculating opacities on-the-fly for any abundance distribution.

Built in 2004, the Skynet robotic telescope network originally consisted of six 0.4 m telescopes located at the Cerro-Tololo Inter-American Observatory in the Chilean Andes. The network was designed to carry out simultaneous multi-wavelength observations of gamma-ray bursts (GRBs) when they are only tens of seconds old. To date, the network has been expanded to ~20 telescopes, including a 20 m radio telescope, that span four continents and five countries. The Campaign Manager (CM) is a new observing mode that has been developed for Skynet. Available to all Skynet observers, the CM semi-autonomously and indefinitely scales and schedules exposures on the observer's behalf while allowing for modification to scaling parameters in real time. The CM is useful for follow up to various transient phenomena including gravitational-wave events, GRB localizations, young supernovae, and eventually, sufficiently bright Argus Optical Array and Large Synoptic Survey Telescope events.

The M31 recurrent nova M31N 1926-07c has had five recorded eruptions. Well-sampled light curves of the two most recent outbursts, in January of 2020 (M31N 2020-01b) and September 2022 (M31N 2022-09a), are presented showing that the photometric evolution of the two events were quite similar, with peak magnitudes of $R=17.2\pm0.1$ and $R=17.1\pm0.1$, and $t_2$ times of $9.7\pm0.9$ and $8.1\pm0.5$ days for the 2020 and 2022 eruptions, respectively. After considering the dates of the four most recent eruptions (where the cycle count is believed to be known), a mean recurrence interval of $\langle P_\mathrm{rec}\rangle=2.78\pm0.03$ years is found, establishing that M31N 1926-07c has one of the shortest recurrence times known.

We cross-correlate positions of galaxies measured in data from the first three years of the Dark Energy Survey with Compton-$y$-maps generated using data from the South Pole Telescope (SPT) and the {\it Planck} mission. We model this cross-correlation measurement together with the galaxy auto-correlation to constrain the distribution of gas in the Universe. We measure the hydrostatic mass bias or, equivalently, the mean halo bias-weighted electron pressure $\langle b_{h}P_{e}\rangle$, using large-scale information. We find $\langle b_{h}P_{e}\rangle$ to be $[0.16^{+0.03}_{-0.04},0.28^{+0.04}_{-0.05},0.45^{+0.06}_{-0.10},0.54^{+0.08}_{-0.07},0.61^{+0.08}_{-0.06},0.63^{+0.07}_{-0.08}]$ meV cm$^{-3}$ at redshifts $z \sim [0.30, 0.46, 0.62,0.77, 0.89, 0.97]$. These values are consistent with previous work where measurements exist in the redshift range. We also constrain the mean gas profile using small-scale information, enabled by the high-resolution of the SPT data. We compare our measurements to different parametrized profiles based on the cosmo-OWLS hydrodynamical simulations. We find that our data are consistent with the simulation that assumes an AGN heating temperature of $10^{8.5}$K but are incompatible with the model that assumes an AGN heating temperature of $10^{8.0}$K. These comparisons indicate that the data prefer a higher value of electron pressure than the simulations within $r_{500c}$ of the galaxies' halos.

Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of James Webb Space Telescope (JWST) CEERS NIRCam images to present the first quantitative identification and characterization of stellar bars at $z>1$ based on rest-frame NIR F444W images of high resolution ($\sim$1.3 kpc at $z\sim$1-3). We identify stellar bars in these images using quantitative criteria based on ellipse fits. For this pilot study, we present six examples of robustly identified bars at $z>1$ with spectroscopic redshifts, including the two highest redshift bars at $z\sim$2.136 and 2.312 quantitatively identified and characterized to date. The stellar bars at $z\sim$1.1-2.3 presented in our study have projected semi-major axes of 2.9-4.3 kpc and projected ellipticities of 0.41-0.53 in the rest-frame NIR. The barred host galaxies have stellar masses $\sim 1 \times 10^{10}$ to $2 \times 10^{11}$ $M_{\odot}$, star formation rates of $\sim$ 21-295 $M_{\odot}$ yr$^{-1}$, and several have potential nearby companions. Our finding of bars at $z\sim$1.1-2.3 demonstrates the early onset of such instabilities and supports simulations where bars form early in massive dynamically cold disks. It also suggests that if these bars at lookback times of 8-10 Gyr survive out to present epochs, bar-driven secular processes may operate over a long time and have a significant impact on some galaxies by $z \sim$0.

The origin of high-redshift quasars and their supermassive black hole engines is unclear. One promising solution is the collapse of a primordial supermassive star. Observational confirmation of this scenario may be difficult, but a general relativistic instability supernova provides one avenue for such. Previous studies have found that a general relativistic instability supernova has a potentially decades-long plateau phase visible to JWST at high redshift. In this work, we examine stars with mass just below the general relativistic instability supernova mass range. These stars pulsate, ejecting a portion of their envelopes. They then contract quasi-statically back to an equilibrium temperature, at which point they again become unstable and pulsate once more. Because each pulse consumes a small amount of the available nuclear fuel, there exists the possibility of multiple pulsations. We present simulations of the contracting phase, the pulsation, and the light-curve phase. We find that the lower mass pulsating models are even brighter than the higher mass supernovae because the pulsations occur in the late helium burning phase when the stars have extremely large radii. The fact that the pulsations are more luminous and occur in a wider mass range than the supernovae bodes well for observation.

We conduct a $^{12}$C$^{16}$O($J$=1-0) (hereafter CO) mapping survey of 64 galaxies in the Fornax cluster using the ALMA Morita array in cycle 5. CO emission is detected from 23 out of the 64 galaxies. Our sample includes dwarf, spiral and elliptical galaxies with stellar masses of $M_{\rm star}\sim10^{6.3-11.6}$~M$_\odot$. The achieved beam size and sensitivity are $15''\times8''$ and $\sim12$~mJy~beam$^{-1}$ at the velocity resolution of $\sim10$~km~s$^{-1}$, respectively. We study the cold-gas (molecular- and atomic-gas) properties of 38 subsamples with $M_{\rm star}>10^9$~M$_\odot$ combined with literature HI data. We find that: (1) the low star-formation (SF) activity in the Fornax galaxies is caused by the decrease in the cold-gas mass fraction with respect to stellar mass (hereafter, gas fraction) rather than the decrease of the SF efficiency from the cold gas; (2) the atomic-gas fraction is more heavily reduced than the molecular-gas fraction of such galaxies with low SF activity. A comparison between the cold-gas properties of the Fornax galaxies and their environmental properties suggests that the atomic gas is stripped tidally and by the ram pressure, which leads to the molecular gas depletion with an aid of the strangulation and consequently SF quenching. Pre-processes in the group environment would also play a role in reducing cold-gas reservoirs in some Fornax galaxies.

Nonlinear effects are crucial for the propagation of Fast Radio Bursts (FRBs) near the source. We study the filamentation of FRBs in the relativistic winds of magnetars, which are commonly invoked as the most natural FRB progenitors. As a result of filamentation, the particle number density and the radiation intensity develop strong gradients along the direction of the wind magnetic field. A steady state is reached when the plasma pressure balances the ponderomotive force. In such a steady state, particles are confined into periodically spaced thin sheets, and electromagnetic waves propagate between them as in a waveguide. We show that: (i) The dispersion relation resembles that in the initial homogeneous plasma, but the effective plasma frequency is determined by the separation of the sheets, not directly by the mean particle density. (ii) The contribution of relativistic magnetar winds to the dispersion measure of FRBs could be several orders of magnitude larger than previously thought. The dispersion measure of the wind depends on the properties of individual bursts (e.g. the luminosity), and therefore can change significantly among different bursts from repeating FRBs. (iii) Induced Compton scattering is suppressed because most of the radiation propagates in near vacuum regions.

Hot Jupiters are the most studied and easily detectable exoplanets for transit observations.However, the correlation between the atmospheric flow and the emission spectra of such planets is still not understood. Due to huge day-night temperature contrast in hot Jupiter, the thermal redistribution through atmospheric circulation has a significant impact on the vertical temperature-pressure structure and on the emission spectra. In the present work, we aim to study the variation of the temperature-pressure profiles and the emission spectra of such planets due to different amounts of atmospheric heat redistribution. For this purpose, we first derive an analytical relation between the heat redistribution parameter f and the emitted flux from the uppermost atmospheric layers of hot Jupiter. We adopt the three possible values of f under isotropic approximation as 1/4, 1/2, and 2/3 for full-redistribution, semi-redistribution and no-redistribution cases respectively and calculate the corresponding temperature-pressure profiles and the emission spectra. Next, we model the emission spectra for different values of f by numerically solving the radiative transfer equations using the discrete space theory formalism. We demonstrate that the atmospheric temperature-pressure profiles and the emission spectra both are susceptible to the values of the heat redistribution function. A reduction in the heat redistribution yields a thermal inversion in the temperature-pressure profiles and hence increases the amount of emission flux. Finally, we revisits the hot Jupiter XO-1b temperature-pressure profile degeneracy case and show that a non-inversion temperature-pressure profile best explains this observed planetary dayside emission spectra.

In the realistic model of cosmic inflation the inflaton potential should be flat and stable under quantum corrections. It is natural to imagine that there is some symmetry behind and an idea of the inflaton as a Nambu-Goldstone boson of spontaneous breaking of some symmetry has been examined. We give a general formulation of this idea using the non-linear realization of Nambu-Goldstone boson in low-energy effective theory with some explicit symmetry breaking to generate non-trivial potential. The potential is naturally a simple function, typically the mass term of inflaton, and the scenario of "warm inflation" should necessarily be applied under the present observational constraints. We investigate the generation mechanism of necessary thermal dissipation term in inflaton field equation for "warm inflation" without large thermal corrections to inflaton potential. A simple numerical analysis is given to show that this scenario can be realistic.

We model the physical parameters of Solar Cycles 23 and 24 using a nonlinear dynamical mean-field dynamo model that includes the formation and evolution of bipolar magnetic regions (BMR). The Parker-type dynamo model consists of a complete MHD system in the mean-field formulation: the 3D magnetic induction equation, and 2D momentum and energy equations in the anelastic approximation. The initialization of BMR is modeled in the framework of Parker's magnetic buoyancy instability. It defines the depths of BMR injections, which are typically located at the edge of the global dynamo waves. The distribution with longitude and latitude and the size of the initial BMR perturbations are taken from the NOAA database of active regions. The tilt of the perturbations is modeled by random function, and the mean tilt is modeled as a near-surface helicity (alpha-effect) term. The data-driven models are compared with the models calculated for random longitudinal and latitudinal distributions of the initial perturbation. The modeling results are compared with various observed characteristics of the solar cycles, including the magnetic butterfly diagram, the polar magnetic and basal magnetic fluxes, and the probability distributions of the BMR flux on the surface. Our results show that BMR can play a substantial role in the dynamo processes, and affect the strength of the solar cycle. However, the data-driven model shows that the BMR effect alone cannot explain the weak Cycle 24. This weak cycle and the prolonged preceding minimum of magnetic activity were probably caused by a decrease of the turbulent helicity in the bulk of the convection zone during the decaying phase of Cycle 23.

We present the result of a broadband (0.5-70 keV) X-ray spectral analysis of the late-merger galaxy Mrk 739, which contains a dual active galactic nucleus (AGN), Mrk 739E and Mrk 739W, with a separation of $\sim$3.4 kpc. The spectra obtained with NuSTAR, Chandra, XMM-Newton and Swift/BAT are simultaneously analyzed by separating the contributions from the two AGNs and extended emission with the Chandra data. To evaluate the reflection components from the AGN tori, we consider two models, a phenomenological one (pexrav and zgauss) and a more physically motivated one (XCLUMPY; Tanimoto et al. 2019). On the basis of the results with XCLUMPY, we find that the AGNs in Mrk 739E and Mrk 739W have intrinsic 2-10 keV luminosities of $1.0 \times 10^{43}$ and $7.5 \times 10^{41}\ \rm{erg}\ \rm{s}^{-1}$ absorbed by hydrogen column densities of $N_{\rm{H}} < 6.5 \times 10^{19}\ \rm{cm}^{-2}$ and $N_{\rm{H}} = 6.9^{+3.2}_{-1.7} \times 10^{21}\ \rm{cm}^{-2}$, respectively. The torus covering fraction of the material with $N_{\rm{H}} > 10^{22} \rm{cm}^{-2}$ in Mrk 739E, $C_{\rm{T}}^{(22)} < 0.50$ at a 90% confidence limit, is found to be smaller than those found for late-merger ultra/luminous infrared galaxies, $C_{\rm{T}}^{(22)} = 0.71\pm0.16$ (mean and standard deviation; Yamada et al. 2021). Considering the small star formation rate of Mrk 739E, we suggest that the gas-to-mass ratio of the host galaxy is an important parameter to determine the circumnuclear environment of an AGN in late merger.

Boron nuclei in cosmic rays (CRs) are believed to be mainly produced by the fragmentation of heavier nuclei, such as carbon and oxygen, via collisions with the interstellar matter. Therefore, the boron-to-carbon flux ratio (B/C) and the boron-to-oxygen flux ratio (B/O) are very essential probes of the CR propagation. The energy dependence of the B/C ratio from previous balloon-borne and space-based experiments can be well described by a single power-law up to about 1 TeV/n within uncertainties. This work reports direct measurements of B/C and B/O in the energy range from 10 GeV/n to 5.6 TeV/n with 6 years of data collected by the Dark Matter Particle Explorer, with high statistics and well controlled systematic uncertainties. The energy dependence of both the B/C and B/O ratios can be well fitted by a broken power-law model rather than a single power-law model, suggesting the existence in both flux ratios of a spectral hardening at about 100 GeV/n. The significance of the break is about $5.6\sigma$ and $6.9\sigma$ for the GEANT4 simulation, and $4.4\sigma$ and $6.9\sigma$ for the alternative FLUKA simulation, for B/C and B/O, respectively. These results deviate from the predictions of conventional turbulence theories of the interstellar medium, which point toward a change of turbulence properties of the interstellar medium (ISM) at different scales or novel propagation effects of CRs, and should be properly incorporated in the indirect detection of dark matter via anti-matter particles.

The titanium abundances of late-type stars are important tracers of Galactic formation history. However, abundances inferred from Ti I and Ti II lines can be in stark disagreement in very metal-poor giants. Departures from local thermodynamic equilibrium (LTE) have a large impact on the minority neutral species and thus influences the ionisation imbalance, but satisfactory non-LTE modelling for both dwarfs and giants has not been achieved in previous literature. The reliability of titanium abundances is reassessed in benchmark dwarfs and giants using a new non-LTE model and one-dimensional (1D) model atmospheres. A comprehensive model atom was compiled with a more extended level structure and newly published data for inelastic collisions between Ti I and neutral hydrogen. In 1D LTE, the Ti I and Ti II lines agree to within $0.06$ dex for the Sun, Arcturus, and the very metal-poor stars HD84937 and HD140283. For the very metal-poor giant HD122563, the Ti I lines give an abundance that is $0.47$ dex lower than that from Ti II. The 1D non-LTE corrections can reach $+0.4$ dex for individual Ti I lines and $+0.1$ dex for individual Ti II lines, and reduce the overall ionisation imbalance to $-0.17$ dex for HD122563. However, it also increases the imbalance for the very metal-poor dwarf and sub-giant to around $0.2$ dex. Using 1D non-LTE reduces the ionisation imbalance in very metal-poor giants but breaks the balance of other very metal-poor stars, consistent with the conclusions in earlier literature. To make further progress, consistent 3D non-LTE models are needed.

A key goal of the Solar Orbiter mission is to connect elemental abundance measurements of the solar wind enveloping the spacecraft with EUV spectroscopic observations of their solar sources, but this is not an easy exercise. Observations from previous missions have revealed a highly complex picture of spatial and temporal variations of elemental abundances in the solar corona. We have used coordinated observations from Hinode and Solar Orbiter to attempt new abundance measurements with the SPICE (Spectral Imaging of the Coronal Environment) instrument, and benchmark them against standard analyses from EIS (EUV Imaging Spectrometer). We use observations of several solar features in AR 12781 taken from an Earth-facing view by EIS on 2020 November 10, and SPICE data obtained one week later on 2020 November 17; when the AR had rotated into the Solar Orbiter field-of-view. We identify a range of spectral lines that are useful for determining the transition region and low coronal temperature structure with SPICE, and demonstrate that SPICE measurements are able to differentiate between photospheric and coronal Mg/Ne abundances. The combination of SPICE and EIS is able to establish the atmospheric composition structure of a fan loop/outflow area at the active region edge. We also discuss the problem of resolving the degree of elemental fractionation with SPICE, which is more challenging without further constraints on the temperature structure, and comment on what that can tell us about the sources of the solar wind and solar energetic particles.

The most significant feature in the cosmic-ray (CR) nuclei spectra is the spectral hardening at a few hundred GV. Whether the hardening of the different nuclei species are same or not is important for constructing CR source and propagation models. In this work, we collect the latest released AMS-02 CR nuclei spectra of primary species (proton, helium, carbon, oxygen, neon, magnesium, silicon, and iron), secondary species (lithium, beryllium, boron, and fluorine), and hybrid species (nitrogen, sodium, and aluminum), and study the break positions and the spectral index differences (less and greater than the break rigidity) of the spectral hardening quantitatively. The results show us that the CR nuclei spectral hardening at a few hundred GV has hybrid origins. In detail, the dominating factors of the spectral hardening for primary and secondary CR nuclei species are different: the former comes from the superposition of different kinds of CR sources, while the latter comes from the propagation process. Both of these factors influence all kinds of CR nuclei spectra, just with different weights.

We predict the 21-cm global signal and power spectra during the Epoch of Reionisation using the MERAXES semi-analytic galaxy formation and reionisation model, updated to include X-ray heating and thermal evolution of the intergalactic medium. Studying the formation and evolution of galaxies together with the reionisation of cosmic hydrogen using semi-analytic models (such as MERAXES) requires N-body simulations within large volumes and high mass resolutions. For this, we use a simulation of side-length $210$ $h^{-1}$ Mpc with $4320^3$ particles resolving dark matter haloes to masses of $5\times10^8$ $h^{-1}$ $M_\odot$. To reach the mass resolution of atomically cooled galaxies, thought the dominant population contributing to reionisation, at $z=20$ ($\sim 2\times10^7$ $h^{-1}$ $M_\odot$) we augment this simulation using the DARKFOREST Monte-Carlo merger tree algorithm to achieve an effective particle count of $\sim10^{12}$. Using this augmented simulation we explore the impact of mass resolution on the predicted reionisation history as well as the impact of X-ray heating on the 21-cm global signal and the 21-cm power-spectra. We also explore the cosmic variance of 21-cm statistics within $70^{3}$ $h^{-3}$ Mpc$^3$ sub-volumes. We find that the midpoint of reionisation varies by $\Delta z\sim0.8$ and that the cosmic variance on the power spectrum is underestimated by a factor of $2-4$ at $k\sim 0.1-0.4$ Mpc$^{-1}$ due to the non-Gaussian signal. To our knowledge, this work represents the first model of both reionisation and galaxy formation which resolves low-mass atomically cooled galaxies while simultaneously sampling sufficiently large scales necessary for exploring the effects of X-rays in the early Universe.

The HD108236 system was first announced with the detection of four small planets based on TESS data. Shortly after, the transit of an additional planet with a period of 29.54d was serendipitously detected by CHEOPS. In this way, HD108236 (V=9.2) became one of the brightest stars known to host five small transiting planets (R$_p$<3R$_{\oplus}$). We characterize the planetary system by using all the data available from CHEOPS and TESS space missions. We use the flexible pointing capabilities of CHEOPS to follow up the transits of all the planets in the system, including the fifth transiting body. After updating the host star parameters by using the results from Gaia eDR3, we analyzed 16 and 43 transits observed by CHEOPS and TESS, respectively, to derive the planets physical and orbital parameters. We carried out a timing analysis of the transits of each of the planets of HD108236 to search for the presence of transit timing variations. We derived improved values for the radius and mass of the host star (R$_{\star}$=0.876$\pm$0.007 R$_{\odot}$ and M$_{\star}$=0.867$_{-0.046}^{+0.047}$ M$_{\odot}$). We confirm the presence of the fifth transiting planet f in a 29.54d orbit. Thus, the system consists of five planets of R$_b$=1.587$\pm$0.028, R$_c$=2.122$\pm$0.025, R$_d$=2.629$\pm$0.031, R$_e$=3.008$\pm$0.032, and R$_f$=1.89$\pm$0.04 [R$_{\oplus}$]. We refine the transit ephemeris for each planet and find no significant transit timing variations for planets c, d, and e. For planets b and f, instead, we measure significant deviations on their transit times (up to 22 and 28 min, respectively) with a non-negligible dispersion of 9.6 and 12.6 min in their time residuals. We confirm the presence of planet f and find no significant evidence for a potential transiting planet in a 10.9d orbital period, as previously suggested. Full abstract in the PDF file.

67P/Churyumov-Gerasimenko is a Jupiter-family comet that was the target of the Rosetta mission, the first mission to successfully orbit and land a probe on a comet. This mission was accompanied by a large ground-based observing campaign. We have developed a pipeline to calibrate and measure photometry of comet 67P during its 2016 perihelion passage, making use of all visible wavelength broadband imaging collected across a wide range of facilities. The pipeline calibrates the brightness of the comet to a common photometric system (Pan-STARRS 1) using background stars within the field allowing for compilation and comparison of multiple data sets. Results follow the predictions based on previous apparitions: 67P shows no obvious change in activity levels from orbit-to-orbit and coma colours remain constant throughout the apparition. We detected an outburst on 2015 August 22 of $\sim$0.14 mag. The brightness and estimated mass of this outburst puts it in line with the outbursts directly observed on the nucleus by Rosetta. An in situ outburst was observed at the same time as the one seen from the ground, however linking these two events directly remains challenging.

The current architecture of a given multi-planetary system is a key fingerprint of its past formation and dynamical evolution history. Long-term follow-up observations are key to complete their picture. In this paper we focus on the confirmation and characterization of the components of the TOI-969 planetary system, where TESS detected a Neptune-size planet candidate in a very close-in orbit around a late K-dwarf star. We use a set of precise radial velocity observations from HARPS, PFS and CORALIE instruments covering more than two years in combination with the TESS photometric light curve and other ground-based follow-up observations to confirm and characterize the components of this planetary system. We find that TOI-969 b is a transiting close-in ($P_b\sim 1.82$ days) mini-Neptune planet ($m_b=9.1^{+1.1}_{-1.0}$ M$_{\oplus}$, $R_b=2.765^{+0.088}_{-0.097}$ R$_{\oplus}$), thus placing it on the {lower boundary} of the hot-Neptune desert ($T_{\rm eq,b}=941\pm31$ K). The analysis of its internal structure shows that TOI-969 b is a volatile-rich planet, suggesting it underwent an inward migration. The radial velocity model also favors the presence of a second massive body in the system, TOI-969 c, with a long period of $P_c=1700^{+290}_{-280}$ days and a minimum mass of $m_{c}\sin{i_c}=11.3^{+1.1}_{-0.9}$ M$_{\rm Jup}$, and with a highly-eccentric orbit of $e_c=0.628^{+0.043}_{-0.036}$. The TOI-969 planetary system is one of the few around K-dwarfs known to have this extended configuration going from a very close-in planet to a wide-separation gaseous giant. TOI-969 b has a transmission spectroscopy metric of 93, and it orbits a moderately bright ($G=11.3$ mag) star, thus becoming an excellent target for atmospheric studies. The architecture of this planetary system can also provide valuable information about migration and formation of planetary systems.

Planets can carve gaps in the surface density of protoplanetary discs. The formation of these gaps can reduce the corotation torques acting on the planets. In addition, gaps can halt the accretion of solids onto the planets as dust and pebbles can be trapped at the edge of the gap. This accumulation of dust could explain the origin of the ring-like dust structures observed using high-resolution interferometry. In this work we provide an empirical scaling relation for the depth of the gap cleared by a planet on an eccentric orbit as a function of the planet-to-star mass ratio $q$, the disc aspect ratio $h$, Shakura-Sunyaev viscosity parameter $\alpha$, and planetary eccentricity $e$. We construct the scaling relation using a heuristic approach: we calibrate a toy model based on the impulse approximation with 2D hydrodynamical simulations. The scaling reproduces the gap depth for moderate eccentricities ($e\leq 4h$) and when the surface density contrast outside and inside the gap is $\leq 10^{2}$. Our framework can be used as the basis of more sophisticated models aiming to predict the radial gap profile for eccentric planets.

The Data Release 4 of the Atacama Cosmology Telescope (ACT) shows an agreement with an Harrison-Zel'dovich primordial spectrum ($n_s=1.009 \pm 0.015$), introducing a tension with a significance of $99.3\%$ CL with the results from the Planck satellite. The discrepancy on the value of the scalar spectral index is neither alleviated with the addition of large scale structure information nor with the low multipole polarization data. We discuss possible avenues to alleviate the tension relying on either neglecting polarization measurements from ACT or in extending the inflationary sector of the theory.

SPICA-FT is part of the CHARA/SPICA instrument which combines a visible 6T fibered instrument (SPICAVIS) with a H-band 6T fringe sensor. SPICA-FT is a pairwise ABCD integrated optics combiner. The chip is installed in the MIRC-X instrument. The MIRC-X spectrograph could be fed either by the classical 6T fibered combiner or by the SPICA-FT integrated optics combiner. SPICA-FT also integrates a dedicated fringe tracking software, called the opd-controller communicating with the main delay line through a dedicated channel. We present the design of the integrated optics chip, its implementation in MIRC-X and the software architecture of the group-delay and phase-delay control loops. The final integrated optics chip and the software have been fully characterized in the laboratory. First on-sky tests of the integrated optics combiner began in 2020. We continue the on-sky tests of the whole system (combiner + software) in Spring and Summer 2022. We present the main results, and we deduce the preliminary performance of SPICA-FT.

The Sun shows a wide range of temporal variations, from a few seconds to decades and even centuries, broadly classified into two classes short-term and Long-term. The solar dynamo mechanism is believed to be responsible for these global changes happening in the Sun. Hence, many dynamo models have been proposed to explain the observed behaviour of the Sun. This thesis is primarily focused on studying the \lt\ variation of the Sun and provides various inputs to the solar dynamo models. With a renewed interest in the subject, several automatic techniques have been developed for extensive data analysis as applied to long-term datasets and presented in this thesis. This approach provides better consistency and eliminates human subjectivity, which has been a normal practice in the past. The variation of penumbra to umbra area ratio, q, observed here, will provide constraints in sunspot simulations. In addition, the absence of any difference in the behaviour of small and big spots does not support the idea of the global and local dynamo. Two classes of BMRs observed in the magnetograms further verify this behaviour. The importance of the NSSL is not studied so well in the context of solar dynamo models, but it will be worth waiting to see its significance for understanding solar dynamo. Finally, the indication of tilt quenching presented here needs to be further verified using the more comprehensive data set, including stronger cycles.

At the ultra-high densities existing in the core of neutron stars, it is expected that a phase transition from baryonic to deconfined quark matter may occur. Such a phase transition would affect the underlying equation of state (EoS) as well as the observable astrophysical properties of neutron stars. Comparison of EoS model predictions with astronomical data from multi-messenger signals then provides us an opportunity to probe the behaviour of dense matter. In this work, we restrict the allowed parameter space of EoS models in neutron stars for both nucleonic (relativistic mean field model) and quark matter (bag model) sectors by imposing state-of-the-art constraints from nuclear calculations, multi-messenger astrophysical data and perturbative QCD (pQCD). We systematically investigate the effect of each constraint on the parameter space of uncertainties using a cut-off filter scheme, as well as the correlations among the parameters and with neutron star astrophysical observables. Using the constraints, we obtain limits for maximum NS mass, maximum central density, as well as for NS radii and tidal deformability. Although pQCD constraints are only effective at very high densities, they significantly reduce the parameter space of the quark model. We also conclude that astrophysical data supports high values of the bag parameter B and disfavors the existence of a pure quark matter core in hybrid stars.

We present observations of the historic transient 4U 1730-22 as observed with the Neutron Star Interior Composition Explorer (NICER). After remaining in quiescence since its 1972 discovery, this X-ray binary showed renewed outburst activity in 2021 and 2022. We observed 4U 1730-22 extensively with NICER, detecting a total of 17 thermonuclear X-ray bursts. From a spectroscopic analysis, we find that these X-ray bursts can be divided into a group of bright and weak bursts. All bright bursts showed $1\sim2$ second rise times and a photospheric radius expansion phase, while the weak bursts showed a slower $\sim5$ second rise with a tendency for concave shapes. From the photospheric radius expansion flux, we estimate the source distance at $6.9\pm0.2$ kpc. We consider various interpretations for our observations and suggest that they may be explained if accreted material is burning stably at the stellar equator, and unstable ignition occurs at a range of higher latitudes.

With a possible angular resolution down to 0.1-0.2 millisecond of arc using the 330 m baselines and the access to the 600-900 nm spectral domain, the CHARA Array is ideally configured for focusing on precise and accurate fundamental parameters of stars. CHARA/SPICA (Stellar Parameters and Images with a Cophased Array) aims at performing a large survey of stars all over the Hertzsprung-Russell diagram. This survey will also study the effects of the different kinds of variability and surface structure on the reliability of the extracted fundamental parameters. New surface-brightness-colour relations will be extracted from this survey, for general purposes on distance determination and the characterization of faint stars. SPICA is made of a visible 6T fibered instrument and of a near-infrared fringe sensor. In this paper, we detail the science program and the main characteristics of SPICA-VIS. We present finally the initial performance obtained during the commissioning.

Theory predicts that the temperature of the X-ray-emitting gas ($\sim$10$^{6}$ K) detected from planetary nebulae (PNe) is a consequence of mixing or thermal conduction when in contact with the ionized outer rim ($\sim$10$^{4}$ K). Gas at intermediate temperatures ($\sim$10$^{5}$ K) can be used to study the physics of the production of X-ray-emitting gas, via C IV, N V and O VI ions. Here we model the stellar atmosphere of the CSPN of NGC 1501 to demonstrate that even this hot H-deficient [WO4]-type star cannot produce these emission lines by photoionization. We use the detection of the C IV lines to assess the physical properties of the mixing region in this PNe in comparison with its X-ray-emitting gas, rendering NGC 1501 only the second PNe with such characterization. We extend our predictions to the hottest [WO1] and cooler [WC5] spectral types and demonstrate that most energetic photons are absorbed in the dense winds of [WR] CSPN and highly ionized species can be used to study the physics behind the production of hot bubbles in PNe. We found that the UV observations of NGC 2452, NGC 6751 and NGC 6905 are consistent with the presence mixing layers and hot bubbles, providing excellent candidates for future X-ray observations.

The long-term cyclicity and temporal succession of glacial-periglacial (or deglacial) periods or epochs are keynotes of Quaternary geology on Earth. Relatively recent work has begun to explore the histories of the mid- to higher-latitudinal terrain of Mars, especially in the northern hemisphere, for evidence of similar cyclicity and succession in the Mid to Late Amazonian Epoch. Here, we carry on with this work by focusing on Protonilus Mensae [PM] (43-490 N, 37-590 E). More specifically, we discuss, describe and evaluate an area within PM that straddles a geological contact between two ancient units: [HNt], a Noachian-Hesperian Epoch transition unit; and [eHT] an early Hesperian Epoch transition unit. Dark-toned terrain within the eHt unit (HiRISE image ESP_028457_2255) shows continuous coverage by structures akin to clastically-sorted circles [CSCs]. The latter are observed in permafrost regions on Earth where the freeze-thaw cycling of surface and/or near-surface water is commonplace and cryoturbation is not exceptional. The crater-size frequency distribution of the dark-toned terrain suggests a minimum age of ~100 Ma and a maximum age of ~1 Ga. The age estimates of the candidate CSCs fall within this dispersion. Geochronologically, this places the candidate CSCs amongst the oldest periglacial landforms identified on Mars so far.

The High Resolution Imager (HRI_EUV) telescope of the Extreme Ultraviolet Imager (EUI) instrument onboard Solar Orbiter has observed EUV brightenings, so-called campfires, as fine-scale structures at coronal temperatures. The goal of this paper is to compare the basic geometrical (size, orientation) and physical (intensity, lifetime) properties of the EUV brightenings with regions of energy dissipation in a non-potential coronal magnetic field simulation. In the simulation, HMI line-of-sight magnetograms are used as input to drive the evolution of solar coronal magnetic fields and energy dissipation. We applied an automatic EUV brightening detection method to EUV images obtained on 30 May 2020 by the HRI_EUV telescope. We applied the same detection method to the simulated energy dissipation maps from the non-potential simulation to detect simulated brightenings. We detected EUV brightenings with density of 1.41x10^{-3} brightenings/Mm^2 in the EUI observations and simulated brightenings between 2.76x10^{-2} - 4.14x10^{-2} brightenings/Mm^2 in the simulation, for the same time range. Although significantly more brightenings were produced in the simulations, the results show similar distributions of the key geometrical and physical properties of the observed and simulated brightenings. We conclude that the non-potential simulation can successfully reproduce statistically the characteristic properties of the EUV brightenings (typically with more than 85% similarity); only the duration of the events is significantly different between observations and simulation. Further investigations based on high-cadence and high-resolution magnetograms from Solar Orbiter are under consideration to improve the agreement between observation and simulation.

Dark deposits visible from orbit appear in the Martian south polar region during the springtime. These are thought to form from explosive jets of carbon dioxide gas breaking through the thawing seasonal ice cap, carrying dust and dirt which is then deposited onto the ice as dark 'blotches', or blown by the surface winds into streaks or 'fans'. We investigate machine learning (ML) methods for automatically identifying these seasonal features in High Resolution Imaging Science Experiment (HiRISE) satellite imagery. We designed deep Convolutional Neural Networks (CNNs) that were trained and tested using the catalog generated by Planet Four, an online citizen science project mapping the south polar seasonal deposits. We validated the CNNs by comparing their results with those of ISODATA (Iterative Self-Organizing Data Analysis Technique) clustering and as expected, the CNNs were significantly better at predicting the results found by Planet Four, in both the area of predicted seasonal deposits and in delineating their boundaries. We found neither the CNNs or ISODATA were suited to predicting the source point and directions of seasonal fans, which is a strength of the citizen science approach. The CNNs showed good agreement with Planet Four in cross-validation metrics and detected some seasonal deposits in the HiRISE images missed in the Planet Four catalog; the total area of seasonal deposits predicted by the CNNs was 27% larger than that of the Planet Four catalog, but this aspect varied considerably on a per-image basis.

We use the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation to study how Lyman Continuum (LyC) photons escape from galaxies and the observational signatures of this escape. We define two classes of LyC leaker: Bursty Leakers and Remnant Leakers, based on their star formation rates (SFRs) that are averaged over 10 Myr (SFR$_{10}$) or 100 Myr (SFR$_{100}$). Both have $f_{\rm esc}>20\%$ and experienced an extreme burst of star formation, but Bursty Leakers have ${\rm SFR_{10}>SFR_{100}}$, while Remnant Leakers have ${\rm SFR_{10}<SFR_{100}}$. The maximum SFRs in these bursts were typically $\sim100$ times greater than the SFR of the galaxy prior to the burst, a rare $2\sigma$ outlier among the general high-redshift galaxy population. Bursty Leakers are qualitatively similar to ionization-bounded nebulae with holes, exhibiting high ionization parameters and typical HII region gas densities. Remnant Leakers show properties of density-bounded nebulae, having normal ionization parameters but much lower HII region densities. Both types of leaker exhibit [CII]$_{\rm 158\mu m}$ deficits on the [CII]-SFR$_{100}$ relation, while only Bursty Leakers show deficits when SFR$_{10}$ is used. We predict that [CII] luminosity and SFR indicators such as H$\alpha$ and M$_{\rm 1500\r{A}}$ can be combined to identify both types of LyC leaker and the mode by which photons are escaping. These predictions can be tested with [CII] observations of known $z=3-4$ LyC leakers. Finally, we show that leakers with $f_{\rm esc}>20\%$ dominate the ionizing photon budget at $z\gtrsim7.5$ but the contribution from galaxies with $f_{\rm esc}<5\%$ becomes significant at the tail-end of reionization.

Strong gravitational lensing is a unique observational tool for studying the dark and luminous mass distribution both within and between galaxies. Given the presence of substructures, current strong lensing observations demand more complex mass models than smooth analytical profiles, such as power-law ellipsoids. In this work, we introduce a continuous neural field to predict the lensing potential at any position throughout the image plane, allowing for a nearly model-independent description of the lensing mass. We apply our method on simulated Hubble Space Telescope imaging data containing different types of perturbations to a smooth mass distribution: a localized dark subhalo, a population of subhalos, and an external shear perturbation. Assuming knowledge of the source surface brightness, we use the continuous neural field to model either the perturbations alone or the full lensing potential. In both cases, the resulting model is able to fit the imaging data, and we are able to accurately recover the properties of both the smooth potential and of the perturbations. Unlike many other deep learning methods, ours explicitly retains lensing physics (i.e., the lens equation) and introduces high flexibility in the model only where required, namely, in the lens potential. Moreover, the neural network does not require pre-training on large sets of labelled data and predicts the potential from the single observed lensing image. Our model is implemented in the fully differentiable lens modeling code Herculens.

The Canadian Hydrogen Mapping Experiment (CHIME) is a radio telescope located in British Columbia, Canada. The large field of view (FOV) of $\sim$ 200 square degrees has enabled the CHIME/FRB instrument to produce the largest FRB catalog to date. The large FOV also allows CHIME/FRB to be an exceptional pulsar and Rotating Radio Transient (RRAT) finding machine, despite saving only the metadata information of incoming Galactic events. We have developed a pipeline to search for pulsars/RRATs using DBSCAN, a clustering algorithm. Output clusters are then inspected by a human for pulsar/RRAT candidates and follow-up observations are scheduled with the more sensitive CHIME/Pulsar instrument. The CHIME/Pulsar instrument is capable of a near-daily search mode observation cadence. We have thus developed the CHIME/Pulsar Single Pulse Pipeline to automate the processing of CHIME/Pulsar search mode data. We report the discovery of 21 new Galactic sources, with 14 RRATs, 6 regular slow pulsars and 1 binary system. Owing to CHIME/Pulsar's daily observations we have obtained timing solutions for 8 of the 14 RRATs along with all the regular pulsars. This demonstrates CHIME/Pulsar's ability at finding timing solutions for transient sources.

We present long-term (4-10 years) trends of light pollution observed at 26 locations, covering rural, intermediate and urban sites, including the three major European metropolitan areas of Stockholm, Berlin and Vienna. Our analysis is based on i) night sky brightness (NSB) measurements obtained with Sky Quality Meters (SQMs) and ii) a rich set of atmospheric data products provided by the European Centre for Medium-Range Weather Forecasts. We describe the SQM data reduction routine in which we filter for moon- and clear-sky data and correct for the SQM "aging" effect using an updated version of the twilight method of Puschnig et al. (2021). Our clear-sky, aging-corrected data reveals short- and long-term (seasonal) variations due to atmospheric changes. To assess long-term anthropogenic NSB trends, we establish an empirical atmospheric model via multi-variate penalized linear regression. Our modeling approach allows to quantitatively investigate the importance of different atmospheric parameters, revealing that surface albedo and vegetation have by far the largest impact on zenithal NSB. Additionally, the NSB is sensitive to black carbon and organic matter aerosols at urban and rural sites respectively. Snow depth was found to be important for some sites, while the total column of ozone leaves impact on some rural places. The average increase in light pollution at our 11 rural sites is 1.7 percent per year. At our nine urban sites we measure an increase of 1.8 percent per year and for the remaining six intermediate sites we find an average increase of 3.7 percent per year. These numbers correspond to doubling times of 41, 39 and 19 years. We estimate that our method is capable of detecting trend slopes shallower/steeper than 1.5 percent per year.

The statistical power of weak lensing measurements is principally driven by the number of high redshift galaxies whose shapes are resolved. Conventional wisdom and physical intuition suggest this is optimised by deep imaging at long (red or near IR) wavelengths, to avoid losing redshifted Balmer break and Lyman break galaxies. We use the synthetic Emission Line EL-COSMOS catalogue to simulate lensing observations using different filters, from various altitudes. Here were predict the number of exposures to achieve a target z > 0.3 source density, using off-the-shelf and custom filters. Ground-based observations are easily better at red wavelengths, as (more narrowly) are space-based observations. However, we find that SuperBIT, a diffraction-limited observatory operating in the stratosphere, should instead perform its lensing-quality observations at blue wavelengths.

Population III stars were the first stars to form after the Big Bang, and are believed to have made the earliest contribution to the metal content of the universe beyond the products of the Big Bang Nucleosynthesis. These stars are theorized to have had extremely short lifespans, and therefore would only be observable at high redshifts ($z \geq 3-17$) and faint apparent magnitudes ($m_{AB} \gtrsim 40$). The direct detection of Population III stars therefore remains elusive. However, the recently launched James Webb Space Telescope (JWST) may be capable of detecting stars in the relevant magnitude range in the event of favorable gravitational lensing. Theoretical models are required to interpret these future observations. In this study, new evolutionary models and non-equilibrium model atmospheres were used to characterize the observable properties of zero-age main sequence Population III stars. The calculated models cover a wide range of possible Population III stellar masses, from the minimum mass predicted by star formation studies to the maximum mass capable of maintaining hydrostatic equilibrium. Synthetic photometry and theoretical color-magnitude diagrams were calculated for the bands of the Near-Infrared Camera (NIRCam) on JWST. The final results are compared to the scales of known lensing events and JWST magnitude limits. The purpose of this study is to calculate the observable parameters of Population III stars in the most optimal JWST bands in order to provide a theoretical foundation for anticipated future observations of this stellar population.

Precise measurements of the boron-to-carbon and boron-to-oxygen ratios by DAMPE show clear hardenings around $100$ GeV/n, which provide important implications on the production, propagation, and interaction of Galactic cosmic rays. In this work we investigate a number of models proposed in literature in light of the DAMPE findings. These models can roughly be classified into two classes, the propagation effect and the source effect. Among these models discussed, we find that the re-acceleration of cosmic rays during their propagation by random magnetohydrodynamic waves may not reproduce significant enough hardenings of B/C and B/O, and an additional spectral break of the diffusion coefficient is required. The other models can properly explain the hardenings of the ratios. However, depending on simplifications assumed, the models differ in their quality in reproducing the data in a wide energy range. The models with significant re-acceleration effect will under-predict low-energy antiprotons but over-predict low-energy positrons, and the models with secondary production at sources over-predict high-energy antiprotons. For all models high-energy positron excess exists.

Few spacecraft have studied the dynamics of Venus' deep atmosphere, which is needed to understand the interactions between the surface and atmosphere. Recent global simulations suggest a strong effect of the diurnal cycle of surface winds on the depth of the planetary boundary layer. We propose to use a turbulent-resolving model to characterize the Venus boundary layer and the impact of surface winds for the first time. Simulations were performed in the low plain and high terrain at the Equator and noon and midnight. A strong diurnal cycle is resolved in the high terrain, with a convective layer reaching 7 km above the local surface and vertical wind of 1.3 m/s. The boundary layer depth in the low plain is consistent with the observed wavelength of the dune fields. At noon, the resolved surface wind field for both locations is strong enough to lift dust particles and engender micro-dunes. Convective vortices are resolved for the first time on Venus.

Solar gravity modes are considered as the {\it Rosetta Stone} for probing and subsequently deciphering the physical properties of the solar inner-most layers. Recent claims of positive detection therefore shed some new light on the long-standing issue of estimating solar gravity mode amplitudes. In this article, our objective is to review the theoretical efforts intended to predict solar gravity mode amplitudes. Because most of these studies assumed analogous driving and damping properties to those for the observed acoustic modes, we also provide a short overview of our current knowledge for these modes in the Sun and solar-type stars (which show solar-like oscillations) before diving into the specific problem of solar gravity modes. Finally, taking recent estimates into account, we conclude and confirm that the low-frequency domain (typically between $10\,\mu$Hz and $100\,\mu$Hz) is certainly more suited to focus on for detecting solar gravity modes. More precisely, around $60\,\mu$Hz, the theoretical estimates are slightly lower than the observational detection threshold as provided by the GOLF (Global Oscillations at Low Frequencies) instrument by about a factor of two only. This is typically within the current uncertainties associated with theoretical estimates and should motivate us for improving our knowledge on turbulence in the whole solar convective region, which is key for improving the accuracy of $g$-mode amplitude estimates. The recent detection of solar inertial modes (Gizon et al. 2021) combined with the continuous development of numerical simulations provide interesting prospects for future studies.

We examine in greater detail five events previously identified as being sources of strong transient coronal outflows in a solar polar region in Hinode/EUV Imaging Spectrometer (EIS) Doppler data. Although relatively compact or faint and inconspicuous in Hinode/Soft X-ray Telescope (XRT) soft-X-ray (SXR) images and in Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA) EUV images, we find that all of these events are consistent with being faint coronal X-ray jets. The evidence for this is that the events result from eruption of minifilaments of projected sizes spanning 5000 -- 14,000 km and with erupting velocities spanning 19 -- 46 km/s, which are in the range of values observed in cases of confirmed X-ray polar coronal hole jets. In SXR images, and in some EUV images, all five events show base brightenings, and faint indications of a jet spire that (in four of five cases where determinable) moves away from the brightest base brightening; these properties are common to more obvious X-ray jets. For a comparatively low-latitude event, the minifilament erupts from near (<~few arcsec) a location of near-eruption-time opposite-polarity magnetic-flux-patch convergence, which again is consistent with many observed coronal jets. Thus, although too faint to be identified as jets a priori, otherwise all five events are identical to typical coronal jets. This suggests that jets may be more numerous than recognized in previous studies, and might contribute substantially to solar wind outflow, and to the population of magnetic switchbacks observed in Parker Solar Probe (PSP) data.

Venus clouds host a convective layer between roughly 50 and 60 km that mixes heat, momentum, and chemical species. Observations and numerical modelling have helped to understand the complexity of this region. However, the impact on chemistry is still not known. Here, we use for the first time a three-dimensional convection-resolving model with passive tracers to mimic SO$_2$ and H$_2$O for two latitudinal cases. The tracers are relaxed towards a vertical profile in agreement with measured values, with a timescale varying over several orders of magnitude. The vertical mixing is quantified, it is strong for a relaxation timescale high in front of the convective timescale, around 4 hours. The spatial and temporal variability of the tracer due to the convective activity is estimated, with horizontal structures of several kilometres. At the Equator, the model is resolving a convective layer at the cloud top (70 km) suggested by some observations, the impact of such turbulent activity on chemical species is accounted for the first time. From the resolved convective plumes, a vertical eddy diffusion is estimated, consistent with past estimations from in-situ measurements, but several orders of magnitude higher than values used in 1D chemistry modelling. The results are compared to observations, with some spatial and temporal variability correlation, suggesting an impact of the convective layer on the chemical species.

T Tauri stars are low-mass pre-main sequence stars that are intrinsically variable. Due to the intense magnetic fields they possess, they develop dark spots on their surface that, because of rotation, introduce a periodic variation of brightness.In addition, the presence of surrounding disks could generate flux variations by variable extinction or accretion. Both can lead to a brightness decrease or increase, respectively. Here, we have compiled a catalog of light curves for 379 T Tauri stars in the Lagoon Nebula (M8) region, using VVVX survey data in the Ks-band. All these stars were already classified as pre-MS stars based on other indicators. The data presented here are spread over a period of about eight years, which gives us a unique follow-up time for these sources at this wavelength. The light curves were classified according to their degree of periodicity and asymmetry, to constrain the physical processes responsible for their variation. Periods were compared with the ones found in literature, on a much shorter baseline. This allowed us to prove that for 126 stars, the magnetically active regions remain stable for several years. Besides, our near-IR data were compared with the optical Kepler/K2 light curves, when available, giving us a better understanding of the mechanisms responsible for the brightness variations observed and how they manifest at different bands. We found that the periodicity in both bands is in fairly good agreement, but the asymmetry will depend on the amplitude of the bursts or dips events and the observation cadence.

Recent astrophysical transient Swift J1913.1+1946 is possibly associated with the gamma-ray burst GRB 221009A at the redshift z=0.151. The transient was accompanied by very high-energy gamma rays up to 18 TeV observed by LHAASO and a photon-like air shower of 251 TeV detected by Carpet-2. These energetic gamma rays cannot reach us from the claimed distance of the source because of the pair production on cosmic background radiation. If the identification and redshift measurements are correct, one would require new physics to explain the data. One possibility invokes axion-like particles (ALPs) which mix with photons but do not attenuate on the background radiation. Here we explore the ALP parameter space and find that the ALP-photon mixing in the Milky Way, and not in the intergalactic space, may help to explain the observations. However, given the low Galactic latitude of the event, misidentification with a Galactic transient remains an undiscarded explanation.

We analyse the connection between the star formation quenching of galaxies and their location in theoutskirts of clusters in the redshift range $z=[0,2]$ by estimating the fraction of red galaxies. More specifically, we focus on galaxies that infall isotropically from those that are infalling alongside filaments. We use a sample of galaxies obtained from the semi-analytic model of galaxy formation SAG applied to the MultiDark simulation. {\textsc{mdpl2}}. In agreement with observational results, we find that the infall regions show levels of star formation that are intermediate between those of galaxies in clusters and in the field. Moreover, we show that, in the redshift range [0-0.85], the quenching of the star formation is stronger in the filamentary region than in the isotropic infall region. We also study the fraction of red galaxies as a function of the normalised distance to the cluster centre and find that, for radii $R/R_{200}> 3 $, the fraction of red galaxies in the filamentary region is considerably larger than in the isotropic infall region. From the analysis of properties of the main progenitors of galaxies identified at $z = 0$, we find that they have different evolutionary behaviours depending on the stellar mass and environment. Our results confirm the observational findings that suggest that the infall regions of clusters play an important role in the pre-processing of galaxies along most of the evolutionary history of galaxies.

The recently re-discovered open cluster Stock 2, located roughly 375 pc away and about 400 Myr old, has the potential to be an exciting new testbed for our understanding of stellar evolution. We present results from a spectroscopic campaign to characterize stars near the cluster's main-sequence turnoff; our goal is to identify candidate chemically peculiar stars among the cluster's A stars. We obtained echelle spectra for 64 cluster members with ESPaDOnS on the 3.6-m Canada-France-Hawaii Telescope, Mauna Kea Observatory, USA, and for six stars with SOPHIE on the 1.93-m telescope at the Observatoire de Haute-Provence, France. We complemented these new observations with those of 13 high-mass cluster members from the HARPS-N archive; our overall sample is of 71 stars. We derived the fundamental parameters (Teff, log g, [M/H]) as well as vsini for our sample using the Sliced Inverse Regression (SIR) technique, and then used iSpec to derive individual abundances of 12 chemical species. With these abundance determinations, we identified nine A stars with anomalous levels of Sc, Ca, and other metallic lines. Follow-up observations of these Am candidates with a known age can transform them into benchmarks for evolutionary models that include atomic diffusion and help build a better understanding of the complex interactions between macroscopic and microscopic processes in stellar interiors.

f(R) gravity is one of the simplest viable modifications to General Relativity: it passes local astrophysical tests, predicts both the early-time cosmic inflation and the late-time cosmic acceleration, and also describes dark matter. In this paper, we probe cosmic magnification on large scales in f (R) gravity, using the well-known Hu-Sawicki model as an example. Our results indicate that at redshifts z < 3, values of the model exponent n > 1 lead to inconsistent behaviour in the evolution of scalar perturbations. Moreover, when relativistic effects are taken into account in the large scale analysis, our results show that as z increases, large-scale changes in the cosmic magnification angular power spectrum owing to integral values of n tend to share a similar pattern, while those of decimal values tend to share another. This feature could be searched for in the experimental data, as a potential ``smoking gun" for the given class of gravity models. Furthermore, we found that at z = 1 and lower, relativistic effects lead to a suppression of the cosmic magnification on large scales in f(R) gravity, relative to the concordance model; whereas, at z > 1, relativistic effects lead to a relative boost of the cosmic magnification. In general, relativistic effects enhance the potential of the cosmic magnification as a cosmological probe.

In flat space and at finite temperature, there are two regimes of false vacuum decay in quantum field theory. At low temperature, the decay proceeds through thermally-assisted tunneling described by periodic Euclidean solutions -- bounces -- with non-trivial time dependence. On the other hand, at high temperature the bounces are time-independent and describe thermal jumps of the field over the potential barrier. We argue that only solutions of the second type are relevant for false vacuum decay catalyzed by a black hole in equilibrium with thermal bath. The argument applies to a wide class of spherical black holes, including $d$-dimensional AdS/dS-Schwarzschild black holes and Reissner-Nordstr\"om black holes sufficiently far from criticality. It does not rely on the thin-wall approximation and applies to multi-field scalar theories.

In this work the influence of the post-Newtonian corrections to the equations of stellar structure is analysed. The post-Newtonian Lane-Emden equation follows from the corresponding momentum density balance equation. From a polytropic equation of state the solutions of the Lane-Endem equations in the Newtonian and post-Newtonian theories are determined and the physical quantities for the \textit{Sun}, for the white dwarf \textit{Sirius B} and for neutron stars with masses $M\simeq1.4M_\odot, 1.8M_\odot$ and $2.0M_\odot$ are calculated. It is shown that the post-Newtonian corrections to the fields of mass density, pressure and temperature are negligible for the \textit{Sun} and \textit{Sirius B}, but for stars with strong fields the differences become important. For the neutron stars analysed here the central pressure and the central temperature which follow from the post-Newtonian Lane-Emden equation are about fifty to sixty percent greater than those of the Newtonian theory and the central mass density is about three to four percent smaller.

We discuss production of QCD axions in a novel scenario, which assumes time-varying scale of Peccei-Quinn symmetry breaking. The latter decreases as the Universe's temperature at early times and eventually stabilises at a large constant value. Such behavior is caused by the portal interaction between the complex field carrying Peccei-Quinn charge and a Higgs-like scalar, which is in thermal equilibrium with primordial plasma. In this scenario, axions are efficiently produced during the parametric resonance decay of the complex Peccei-Quinn field, relaxing to the minimum of its potential in the radiation-dominated stage. Notably, this process is not affected by the Universe's expansion rate and allows to generate the required abundance of dark matter independently of an axion mass. This remains the case in the narrow parametric resonance regime corresponding to small oscillation amplitudes. Phenomenological constraints on the model parameter space strongly depend on the number density of radial field fluctuations, which are also generically excited along with axions. Even if produced with a small abundance relative to axions, radial fluctuations can have a sizeable impact on the BBN and the CMB through their decay products contributing to dark radiation. Existing constraints on dark radiation bound Peccei-Quinn symmetry breaking scale below $10^{8}-10^{9}~\mbox{GeV}$ - this paves the way for ruling out our scenario with the near future searches for axions.

The physics potential of detecting $^8$B solar neutrinos is exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the potential low background level, $^8$B solar neutrinos would be observable in the CC and NC interactions on $^{13}$C for the first time. By virtue of optimized event selections and muon veto strategies, backgrounds from the accidental coincidence, muon-induced isotopes, and external backgrounds can be greatly suppressed. Excellent signal-to-background ratios can be achieved in the CC, NC and ES channels to guarantee the $^8$B solar neutrino observation. From the sensitivity studies performed in this work, we show that one can reach the precision levels of 5%, 8% and 20% for the $^8$B neutrino flux, $\sin^2\theta_{12}$, and $\Delta m^2_{21}$, respectively, using ten years of JUNO data. It would be unique and helpful to probe the details of both solar physics and neutrino physics. In addition, when combined with SNO, the world-best precision of 3% is expected for the $^8$B neutrino flux measurement.

Fundamental constants such as masses and coupling constants of elementary particles can have small temporal and spatial variations in the scalar field dark matter model. These variations entail time oscillations of other constants, such as the Bohr and nuclear magnetons, Bohr radius and the hyperfine structure constant. In the presence of an external magnetic field, these oscillations induce hyperfine transitions in atoms and molecules. We determine the probability of magnetic dipole hyperfine transitions, caused by the oscillating fundamental constants, and propose an experiment that could detect the scalar field dark matter through this effect. This experiment may be sensitive to the scalar field dark matter with mass in the range $1\,\mu\text{eV}<m_\phi<100\,\mu\text{eV}$.

Recently, arXiv:2210.05659 shows that a photon-axion like particle (ALP) oscillation can boost the survival rate of the high energy photons associated with the observed GRB221009A event. Here, we show that the proposed ALP is consistent with the electroweak axion with an anomaly free $Z_{10}$ Froggatt-Nielsen symmetry.

We present a general method for estimating the number of particles impinging on a segmented counter or, in general, on a counter with sub-units. We account for unresolved particles, i.e., the effect of two or more particles hitting the same sub-unit almost simultaneously. To achieve full time resolution we account for the dead time that occurs after the first time-bin of a particle signal. This general counting method can be applied to counting muons in existing detectors like the Underground Muon Detector of the Pierre Auger Observatory. We therefore use the latter as a study case to test the performance of our method and to compare it to other methods from literature. Our method proves to perform with little bias, and also provides an estimate of the number of particles as a function of time (as seen by the detector) to a single time-bin resolution. In this context, the new method can be useful for reconstructing parameters sensitive to cosmic ray mass, which are key to unveiling the origin of cosmic rays.

Extrasolar circumbinary planets are so called because they orbit two stars instead of just one; to date, an increasing number of such planets have been discovered with a variety of techniques. If the orbital frequency of the hosting stellar pair is much higher than the planetary one, the tight stellar binary can be considered as a matter ring current generating its own post-Newtonian stationary gravitomagnetic field through its orbital angular momentum. It affects the orbital motion of a relatively distant planet with Lense-Thirring-type precessional effects which, under certain circumstances, may amount to a significant fraction of the static, gravitoelectric ones, analogous to the well known Einstein perihelion precession of Mercury, depending only on the masses of the system's bodies. Instead, when the gravitomagnetic field is due solely to the spin of each of the central star(s), the Lense-Thirring shifts are several orders of magnitude smaller than the gravitoelectric ones. In view of the growing interest in the scientific community about the detection of general relativistic effects in exoplanets, the perspectives of finding new scenarios for testing such a further manifestation of general relativity might be deemed worth of further investigations.

Geoneutrinos are a unique tool that brings to the surface information about our planet, in particular, its radiogenic power, insights formation and chemical composition. To date, only the KamLAND and Borexino experiments observed geoneutrino, with the former characterized by low concentration of heat-producing elements in the Earth in contrast to the latter that sets tight upper limits on the power of a georeactor hypothesized. With respect to the results yielded therefrom, a small discrepancy has been identified. On this account, next generation experiments like JUNO are needed if it is to provide definitive results with respect to the Earth's radiogenic power, and to fully exploit geoneutrinos to better understand deep Earth. An accurate a priori prediction of the crustal contribution plays an important role in enabling the translation of a particle physics measurement into geo-scientific questions. The existing GIGJ model of JUNO only focused on constructing a geophysical model of the local crust, without local geochemical data. Another existing JULOC includes both data, but only able to be achieved for the top layer of the upper crust, not in deep vertical. This paper reports on the development of JUNO's first 3-D integrated model, JULOC-I, which combines seismic, gravity, rock sample and thermal flow data with new building method, solved the problem in vertical depth. JULOC-I results show higher than expected geoneutrino signals are mainly attributable to higher U and Th in southern China than that found elsewhere on Earth. Moreover, the high level of accuracy of the JULOC-I model, complemented by 10 years of experimental data, indicates that JUNO has an opportunity to test different mantle models. Predictions by JULOC-I can be tested after JUNO goes online and higher accuracy local crustal model continue to play an important role to improve mantle measurements precision.

We investigate the collision-induced flavor instability in homogeneous, isotropic, dense neutrino gases in the two-flavor mixing scenario with energy-dependent scattering. We find that the growth rate of such an instability, if it exists, is the negative average of the flavor-decohering collision rates of the neutrinos weighted by the electron lepton number distribution of the neutrinos. This growth rate is independent of the neutrino mass-splitting, the matter density, and the neutrino density, although the initial amplitude of the unstable oscillation mode can be suppressed by a large matter density. Our results suggest that neutrinos are likely to experience collision-induced flavor conversions deep inside a core-collapse supernova even when both the fast and slow collective flavor oscillations are suppressed. However, this phenomenon may not occur in a neutron star merger because the electron antineutrinos have a larger average energy and more abundance than the electron neutrinos in such an environment.

The interaction of binary black hole mergers with their environments can be studied using numerical relativity simulations. These start only a short finite time before merger, at which point appropriate initial conditions must be imposed. A key task is therefore to identify the configuration that is appropriate for the binary and its environment at this stage of the evolution. In this work we study the behaviour of wave dark matter around equal mass black hole binaries, finding that there is a preferred, quasi-stationary profile that persists and grows over multiple orbits, in contrast to heavier mass dark matter where any overdensity tends to be dispersed by the binary motion. Whilst different initial configurations converge to the preferred quasi-stationary one after several orbits, unwanted transient oscillations are generated in the process, which may impact on the signal in short simulation runs. We also point out that naively superimposing the matter onto a circular binary results in artificially eccentric orbits due to the matter backreaction, which is an effect of the initial conditions and not a signature of dark matter. We discuss the further work required so that comparison of waveforms obtained with environments to vacuum cases can be done in a meaningful way.

Gravitational waves provide us with an extraordinary tool to study the matter inside neutron stars. In particular, the postmerger signal probes an extreme temperature and density regime and will help reveal information about the equation of state of supranuclear-dense matter. Although current detectors are most sensitive to the signal emitted by binary neutron stars before the merger, the upgrades of existing detectors and the construction of the next generation of detectors will make postmerger detections feasible. For this purpose, we present a new analytical, frequency-domain model for the inspiral-merger-postmerger signal emitted by binary neutron stars systems. The inspiral and merger part of the signals are modeled with IMRPhenomD$\_$NRTidalv2, and we describe the main emission peak of postmerger with a three-parameter Lorentzian, using two different approaches: one in which the Lorentzian parameters are kept free, and one in which we model them via quasi-universal relations. We test the performance of our new complete waveform model in parameter estimation analyses, both with simulated signals and numerical relativity waveforms. We investigate the performance of different detector networks to determine the improvement that future detectors will bring to our analysis. We consider Advanced LIGO+ and Advanced Virgo+, KAGRA, and LIGO-India. We also study the possible impact of a detector with high sensitivity in the kilohertz band like NEMO, and finally we compare these results to the ones we obtain with third-generation detectors, the Einstein Telescope and the Cosmic Explorer.