Papers for Monday, Sep 12 2022

The dispersion in chemical abundances provides a very strong constraint on the processes that drive the chemical enrichment of galaxies. Due to its proximity, the spiral galaxy M33 has been the focus of numerous chemical abundance surveys to study the chemical enrichment and dispersion in abundances over large spatial scales. The CHemical Abundances Of Spirals (CHAOS) project has observed $\sim$100 H II regions in M33 with the Large Binocular Telescope (LBT), producing the largest homogeneous sample of electron temperatures (T$_e$) and direct abundances in this galaxy. Our LBT observations produce a robust oxygen abundance gradient of $-$0.037 $\pm$ 0.007 dex/kpc and indicate a relatively small (0.043 $\pm$ 0.015 dex) intrinsic dispersion in oxygen abundance relative to this gradient. The dispersions in N/H and N/O are similarly small and the abundances of Ne, S, Cl, and Ar relative to O are consistent with the solar ratio as expected for $\alpha$-process or $\alpha$-process-dependent elements. Taken together, the ISM in M33 is chemically well-mixed and homogeneously enriched from inside-out with no evidence of significant abundance variations at a given radius in the galaxy. Our results are compared to those of the numerous studies in the literature, and we discuss possible contaminating sources that can inflate abundance dispersion measurements. Importantly, if abundances are derived from a single T$_e$ measurement and T$_e$-T$_e$ relationships are relied on for inferring the temperature in the unmeasured ionization zone, this can lead to systematic biases which increase the measured dispersion up to 0.11 dex.

We consider the effects of an attractive, long-range Yukawa interaction between baryons and dark matter (DM), focusing in particular on temperature and pulsar timing observations of neutron stars (NSs). We show that such a fifth force, with strength modestly stronger than gravity at ranges greater than tens of kilometers (corresponding to mediator masses less than $10^{-11} \text{eV}$), can dramatically enhance dark matter kinetic heating, capture, and pulsar timing Doppler shifts relative to gravity plus short range interactions alone. Using the coldest observed NS and pulsar timing array (PTA) data, we derive limits on fifth force strength over a DM mass range spanning light dark matter up to order solar mass composite DM objects. We also consider an indirect limit by combining bullet cluster limits on the DM self-interaction with weak equivalence principle test limits on baryonic self-interactions. We find the combined indirect limits are moderately stronger than kinetic heating and PTA limits, except when considering a DM subcomponent.

Calibrating the redshift distributions of photometric galaxy samples is essential in weak lensing studies. The self-calibration method combines angular auto- and cross-correlations between galaxies in multiple photometric redshift (photo-$z$) bins to reconstruct the scattering rates matrix between redshift bins. In this paper, we test a recently proposed self-calibration algorithm using the DECaLS Data Release 9 and investigate to what extent the scattering rates are determined. We first mitigate the spurious angular correlations due to imaging systematics by a machine learning based method. We then improve the algorithm for $\chi^2$ minimization and error estimation. Finally, we solve for the scattering matrices, carry out a series of consistency tests and find reasonable agreements: (1) finer photo-$z$ bins return a high-resolution scattering matrix, and it is broadly consistent with the low-resolution matrix from wider bins; (2) the scattering matrix from the Northern Galactic Cap is almost identical to that from Southern Galactic Cap; (3) the scattering matrices are in reasonable agreement with those constructed from the power spectrum and the weighted spectroscopic subsample. We also evaluate the impact of cosmic magnification. Although it changes little the diagonal elements of the scattering matrix, it affects the off-diagonals significantly. The scattering matrix also shows some dependence on scale cut of input correlations, which may be related to a known numerical degeneracy between certain scattering pairs. This work demonstrates the feasibility of the self-calibration method in real data and provides a practical alternative to calibrate the redshift distributions of photometric samples.

Observations of the Milky Way at TeV-PeV energies reveal a bright diffuse flux of hadronic cosmic rays and also bright point sources of gamma rays. If the gamma-ray sources are hadronic cosmic-ray accelerators, then they must also be neutrino sources. However, no neutrino sources have been detected. Where are they? We introduce a new population-based approach to probe Milky Way hadronic PeVatrons, demanding consistency between diffuse and point-source PeV-range data on cosmic rays, gamma rays, and neutrinos. For the PeVatrons, two extreme scenarios are allowed: (1) the hadronic cosmic-ray accelerators and the gamma-ray sources are the same objects, so that bright neutrino sources exist and improved telescopes can detect them, versus (2) the hadronic cosmic-ray accelerators and the gamma-ray sources are distinct, so that there are no detectable neutrino sources. The latter case is possible if hadronic accelerators have sufficiently thin column densities. We quantify present constraints and future prospects, showing how to reveal the nature of the hadronic PeVatrons.

In the third paper from the COOL-LAMPS Collaboration, we report the discovery of COOL J0542-2125, a gravitationally lensed quasar at $z=1.84$, observed as three images due to an intervening massive galaxy cluster at $z=0.61$. The lensed quasar images were identified in a search for lens systems in recent public optical imaging data and have separations on the sky up to 25".9, wider than any previously known lensed quasar. The galaxy cluster acting as a strong lens appears to be in the process of merging, with two sub-clusters separated by $\sim 1$ Mpc in the plane of the sky, and their central galaxies showing a radial velocity difference of $\sim 1000$ km/s. Both cluster cores show strongly lensed images of an assortment of background sources, as does the region between them. A preliminary strong lens model implies masses of $M(<250\ \rm{kpc}) = 1.79^{+0.16} _{-0.01} \times 10^{14} M_{\odot}$ and $M(<250\ \rm{kpc}) = 1.48^{+0.04}_{-0.10} \times 10^{14} M_{\odot}$ for the East and West sub-clusters, respectively. This line of sight is also coincident with a ROSAT ALL-sky Survey source, centered between the two confirmed cluster halos reminiscent of other major cluster-scale mergers.

The kinetic Sunyaev-Zeldovich (kSZ) effect will be an important source of cosmological and astrophysical information in upcoming surveys of the cosmic microwave background (CMB). However, the kSZ effect will also act as the dominant source of noise for several other measurements that use small angular scales in CMB temperature maps, since its blackbody nature implies that standard component separation techniques cannot be used to remove it from observed maps. In this paper, we explore the idea of "de-kSZing": constructing a template for the late-time kSZ effect using external surveys of large-scale structure, and then subtracting this template from CMB temperature maps in order to remove some portion of the kSZ signal. After building intuition for general aspects of the de-kSZing procedure, we perform forecasts for the de-kSZing efficiency of several large-scale structure surveys, including BOSS, DESI, Roman, MegaMapper, and PUMA. We also highlight potential applications of de-kSZing to cosmological constraints from the CMB temperature power spectrum, CMB lensing reconstruction, and the moving-lens effect. While our forecasts predict achievable de-kSZing efficiencies of 10-20% at best, these results are specific to the de-kSZing formalism adopted in this work, and we expect that higher efficiencies are possible using improved versions of this formalism.

In this work we derive constraints on interacting dark matter-dark radiation models from a full-shape analysis of BOSS-DR12 galaxy clustering data, combined with Planck legacy cosmic microwave background (CMB) and baryon acoustic oscillation (BAO) measurements. We consider a set of models parameterized within the effective theory of structure formation (ETHOS), quantifying the lifting of the $S_8$ tension in view of KiDS weak-lensing results. The most favorable scenarios point to a fraction $f\sim 10-100\%$ of interacting dark matter as well as a dark radiation temperature that is smaller by a factor $\xi\sim 0.1-0.15$ compared to the CMB, leading to a reduction of the tension to the $\sim 1\sigma$ level. The temperature dependence of the interaction rate favored by relaxing the $S_8$ tension is realized for a weakly coupled unbroken non-Abelian $SU(N)$ gauge interaction in the dark sector. To map our results onto this $SU(N)$ model, we compute higher-order corrections due to Debye screening. We find a lower bound $\alpha_d\equiv g_d^2/(4\pi)\gtrsim 10^{-8} (10^{-9})$ for dark matter mass $1000 (1)$ GeV for relaxing the $S_8$ tension, consistent with upper bounds from galaxy ellipticities and compatible with self-interactions relevant for small-scale structure formation.

A star destroyed by a supermassive black hole (SMBH) in a tidal disruption event (TDE) enables the study of SMBHs. We propose that the distance within which a star is completely destroyed by a SMBH, defined $r_{\rm t, c}$, is accurately estimated by equating the SMBH tidal field (including numerical factors) to the maximum gravitational field in the star. We demonstrate that this definition accurately reproduces the critical $\beta_{\rm c} = r_{\rm t}/r_{\rm t, c}$, where $r_{\rm t} = R_{\star}\left(M_{\bullet}/M_{\star}\right)^{1/3}$ is the standard tidal radius with $R_{\star}$ and $M_{\star}$ the stellar radius and mass and $M_{\bullet}$ the SMBH mass, for multiple stellar progenitors at various ages, and can be reasonably approximated by $\beta_{\rm c} \simeq \left[\rho_{\rm c}/(4\rho_{\star})\right]^{1/3}$, where $\rho_{\rm c}$ ($\rho_{\star}$) is the central (average) stellar density. We also calculate the peak fallback rate and time at which the fallback rate peaks, finding excellent agreement with hydrodynamical simulations, and also suggest that the partial disruption radius -- the distance at which any mass is successfully liberated from the star -- is $\beta_{\rm partial} \simeq 4^{-1/3} \simeq 0.6$. For given stellar and SMBH populations, this model yields, e.g., the fraction of partial TDEs, the peak luminosity distribution of TDEs, and the number of directly captured stars.

The correspondence between galaxy major mergers and starburst activity is well-established observationally and in simulations of low redshift galaxies. However, the evolution of the properties of interactions and of the galaxies involved suggests that the starburst response of galaxies to merger events could vary across cosmic time. Using the VINTERGATAN cosmological zoom-in simulation of a Milky Way-like galaxy, we show here that starbursts, i.e. episodes of fast star formation, are connected with the onset of tidal compression, itself induced by mergers. However, this compression becomes strong enough to trigger starbursts only after the formation of the galactic disc. As a consequence, starburst episodes are only found during a precise phase of galaxy evolution, after the formation of the disc and until the last major merger. As the depletion time quantifies the instantaneous star formation activity, while the specific star formation rate involves the integrated result of the past activity (via the stellar mass), starburst episodes do not necessarily coincide with elevated specific star formation rate. This suggests that not all starburst galaxies are outliers above the main sequence of galaxy formation.

The diagram comparing the flux ratio of the [OIII] and H$\beta$ emission lines with the total stellar mass of galaxies (also known as the mass-excitation diagram, MEx) has been widely used to classify the ionization mechanism in high redshift galaxies between star formation and active galactic nuclear ones. This diagram was mainly derived using single-fiber spectroscopy from the SDSS-DR7 survey. In this study, we revise this diagram using the central and integrated spectral measurement from the entire Integral Field Spectroscopic MaNGA sample. Our results suggest that along with the physical parameters of this diagram, the equivalent width of the H$\alpha$ emission line is also required to constrain the ionization mechanism of a high-redshifted galaxy. Furthermore, the location of a galaxy in the excitation-mass diagram varies depending on the use of central or integrated properties.

We present a catalog of continuum and emission line properties for 750,414 broad-line quasars included in the Sloan Digital Sky Survey Data Release 16 quasar catalog (DR16Q), measured from optical spectroscopy. These quasars cover broad ranges in redshift ($0.1\lesssim z\lesssim 6$) and luminosity ($44\lesssim \log (L_{\rm bol}/{\rm erg\,s^{-1}})\lesssim 48$), and probe lower luminosities than an earlier compilation of SDSS DR7 quasars. Derived physical quantities such as single-epoch virial black hole masses and bolometric luminosities are also included in this catalog. We present improved systemic redshifts and realistic redshift uncertainties for DR16Q quasars using the measured line peaks and correcting for velocity shifts of various lines with respect to the systemic velocity. About 1%, 1.4%, and 11% of the original DR16Q redshifts deviate from the systemic redshifts by $|\Delta V|>1500\,{\rm km\,s^{-1}}$, $|\Delta V|\in [1000,1500]\,{\rm km\,s^{-1}}$, and $|\Delta V|\in [500,1000]\,{\rm km\,s^{-1}}$, respectively; about $1900$ DR16Q redshifts were catastrophically wrong ($|\Delta V|>10,000\,{\rm km\,s^{-1}}$). We demonstrate the utility of this data product in quantifying the spectral diversity and correlations among physical properties of quasars with large statistical samples.

We present the analysis of $\sim 100$pc-scale compact radio continuum sources detected in 63 local (Ultra) Luminous Infrared Galaxies (U/LIRGs; $L_{\rm IR} \ge 10^{11} L_\odot$), using FWHM $\lesssim 0''.1 - 0''.2$ resolution 15 and 33 GHz observations with the Karl G. Jansky Very Large Array. We identify a total of 133 compact radio sources with effective radii of 8 - 170pc, which are classified into four main categories -- "AGN" (AGN), "AGN/SBnuc" (AGN-starburst composite nucleus), "SBnuc" (starburst nucleus) and "SF" (star-forming clumps) -- based on ancillary datasets and the literature. We find that "AGN" and "AGN/SBnuc" more frequently occur in late-stage mergers and have up to 3 dex higher 33 GHz luminosities and surface densities compared with "SBnuc" and "SF", which may be attributed to extreme nuclear starburst and/or AGN activity in the former. Star formation rates (SFRs) and surface densities ($\Sigma_{\rm SFR}$) are measured for "SF" and "SBnuc" using both the total 33 GHz continuum emission (SFR $\sim 0.14 - 13$ M$_\odot$ yr$^{-1}$, $\Sigma_{\rm SFR} \sim 13 - 1600$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$) and the thermal free-free emission from HII regions (median SFR$_{\rm th} \sim 0.4$ M$_\odot$ yr$^{-1}$, $\Sigma_{\rm SFR_{th}} \sim 44$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$). These values are 1 - 2 dex higher than those measured for similar-sized clumps in nearby normal (non-U/LIRGs). The latter also have much flatter median 15 - 33 GHz spectral index ($\sim -0.08$) compared with "SBnuc" and "SF" ($\sim -0.46$), which may reflect higher non-thermal contribution from supernovae and/or ISM densities in local U/LIRGs that directly result from and/or lead to their extreme star-forming activities on 100\,pc scales.

The Exoplanet Modeling and Analysis Center (EMAC) at NASA Goddard Space Flight Center is a web-based catalog, repository, and integration platform for modeling and analysis resources focused on the study of exoplanet characteristics and environments. EMAC hosts user-submitted resources ranging in category from planetary interior models to data visualization tools. Other features of EMAC include integrated web tools developed by the EMAC team in collaboration with the tools' original author(s) and video demonstrations of a growing number of hosted tools. EMAC aims to be a comprehensive repository for researchers to access a variety of exoplanet resources that can assist them in their work, and currently hosts a growing number of code bases, models, and tools. EMAC is a key project of the NASA GSFC Sellers Exoplanet Environments Collaboration (SEEC) and can be accessed at https://emac.gsfc.nasa.gov.

We present a method to reconstruct the longitudinal profile of electrons in showers using Cherenkov telescopes. We show how the Cherenkov light collected by an array of telescopes can be transformed into the number of electrons as a function of atmospheric depth. This method is validated using air shower and simplified telescope simulations. The reconstruction of the depth in which the shower has the maximum number of electrons ($\mathrm{X_{max}}$) opens the possibility of cosmic ray composition studies with Cherenkov telescopes in the energy range from 10 to 100 TeV. A resolution of less than 16 $\mathrm{g/cm^{2}}$ in the $\mathrm{X_{max}}$ reconstruction is obtained.

Understanding the occurrence of Earth-sized planets in the habitable zone of Sun-like stars is essential to the search for Earth analogues. Yet a lack of reliable Kepler detections for such planets has forced many estimates to be derived from the close-in ($2<P_{\mathrm{orb}}<100$ days) population, whose radii may have evolved differently under the effect of atmospheric mass loss mechanisms. In this work, we compute the intrinsic occurrence rates of close-in super-Earths ($\sim1-2\,R_\oplus$) and sub-Neptunes ($\sim2-3.5\,R_\oplus$) for FGK stars ($0.56-1.63\,M_\odot$) as a function of orbital period and find evidence of two regimes: where super-Earths are more abundant at short orbital periods, and where sub-Neptunes are more abundant at longer orbital periods. We fit a parametric model in five equally populated stellar mass bins and find that the orbital period of transition between these two regimes scales with stellar mass, like $P_\mathrm{trans} \propto M_*^{1.7\pm0.2}$. These results suggest a population of former sub-Neptunes contaminating the population of Gyr-old close-in super-Earths, indicative of a population shaped by atmospheric loss. Using our model to constrain the long-period population of intrinsically rocky planets, we estimate an occurrence rate of $\Gamma_\oplus = 15^{+6}_{-4}\%$ for Earth-sized habitable zone planets, and predict that sub-Neptunes may be $\sim$twice as common as super-Earths in the habitable zone (when normalized over the natural log orbital period and radius range used). Finally, we discuss our results in the context of future missions searching for habitable zone planets.

The passage of sungrazing comets in the solar corona can be a powerful tool to probe the local plasma properties. Here, we carry out a study of the striae pattern appearing in the tail of sungrazing Comet Lovejoy, as observed by the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO) during the inbound and outbound phases of the comet orbit. We consider the images in EUV in the 171 {\AA} bandpass, where emission from oxygen ions O$^{4+}$ and O$^{5+}$ is found. The striae are described as due to a beam of ions injected along the local magnetic field, with the initial beam velocity decaying because of collisions. Also, ion collisional diffusion contributes to ion propagation. Both the collision time for velocity decay and the diffusion coefficient for spatial spreading depend on the ambient plasma density. A probabilistic description of the ion beam density along the magnetic field is developed, where the beam position is given by the velocity decay and the spreading of diffusing ions is described by a Gaussian probability distribution. Profiles of emission intensity along the magnetic field are computed and compared with the profiles along the striae observed by AIA, showing a good agreement for most considered striae. The inferred coronal densities are then compared with a hydrostatic model of the solar corona. The results confirm that the coronal density is strongly spatially structured.

We extend the results obtained in \cite{Piattella_2016, mcvittie_2015} for gravitational lensing in the McVittie metric by including the effect of the transition from the matter-dominated epoch of the Universe to the $\Lambda$-dominated era. We derive a formula that agrees with the previous results for the McVittie metric at lowest order, and compare the lensing angle predictions obtained from the Schwarzschild approximation, the McVittie model and higher order corrections to the McVittie model. In doing this, we test if, beyond the correction from the accelerated expansion of the Universe, there is a need for including the matter content of the Universe in modeling lens systems at the redshifts observerd in lens systems. We investigate if there is a need for a modification of the lens equation from these corrections, and if so, to which order and whether it is measurable. We find that while the effect is of the same order as the one calculated previously, there is no significant contribution to the bending angle, as the 1st order effect is already of order $\mathcal{O}(\theta_O^4)$ in the observed angle.

Some Cataclysmic Variables (CVs) exhibits a very long photometric period (VLPP). We calculate the properties of a hypothetical third body, initially assumed on circular--planar orbit, by matching the modelled VLPP to the observed one of four CVs studied here: {\sl LU Camelopardalis} (LU Cam), QZ Serpentis (QZ Ser), V1007 Herculis (V1007 Her) and BK Lyncis (BK Lyn). The eccentric and low inclination orbits for a third body are considered using analytical results. The chosen parameters of the binary components are based on the orbital period of each CV. The smallest corresponding semi-major axis permitted before the third body's orbit becomes unstable is also calculated. A first-order analytical post-Newtonian correction is applied, and the rate of precession of the pericentre is found, but it can not explain any of the observed VLPP. For the first time, we also estimate the effect of secular perturbations by this hypothetical third body on the mass transfer rate of such CVs. We made sure that the observed and calculated amplitude of variability was comparable too. The mass of the third body satisfying all constrains range from 0.63 to 97 Jupiter masses. Our results show further evidence supporting the hypothesis of a third body in three of these CVs, but only marginally in V1007 Her.

We present and analyze 58 transit light curves of TrES-3b and 98 transit light curves of Qatar-1b observed by Transiting Exoplanet Survey Satellite (TESS), plus two transit light curves of Qatar-1b observed by us using a ground-based 1.23\,m telescope. These light curves are combined with the best-quality light curves taken from the Exoplanet Transit Database (ETD) and literature. The precisely determined mid-transit times from these light curves enable us to obtain the refined orbital ephemerides with improved precision for both hot Jupiters. From the timing analysis, we find an indication for the presence of transit timing variations (TTVs) in both systems. Since the observed TTVs are unlikely to be short-term and periodic, the possibility of additional planets in the orbits close to TrES-3b and Qatar-1b are ruled out. Possible causes of long-term TTVs such as orbital decay, apsidal precession, the Applegate mechanism and line-of-sight acceleration are also examined. However, none of these possibilities are found to explain the observed TTV of TrES-3b. In contrast to this, the line-of-sight acceleration appears to be a plausible explanation for the observed TTV of Qatar-1b. In order to confirm these findings, further high-precision transit and RV observations of both systems would be worthwhile.

The multi-band photometry of the VOICE imaging data enables both shape measurement and photometric redshift estimation which are the two essential quantities for weak lensing analysis. We estimate the Excess Surface Density (ESD; $\Delta\Sigma$) based on galaxy-galaxy measurements around galaxies at lower redshift (0.10<$z_l$<0.35 ) while we select the background sources to be at higher redshift ranging from 0.3 to 1.5. The foreground galaxies are further divided into two major categories according to their color (blue/red), each of which has been further divided into high/low stellar mass bins. Then the halo masses of the samples are estimated by modeling the signals, and the posterior of the parameters are samples via Mote Carlo Markov Chain (MCMC) process. We compare our results with the existing Stellar-to-Halo Mass Relation (SHMR) and find that the blue low stellar mass bin deviates from the SHMR relation whereas all other three samples agrees well with empirical curves. We interpret this discrepancy as the effect of a low star formation efficiency of the low-mass blue dwarf galaxy population dominated in the VOICE-CDFS area.

Using the first epoch of four-band NIRCam observations obtained by the James Webb Space Telescope Prime Extragalactic Areas for Reionization and Lensing Science Program in the Spitzer IRAC Dark Field, we search for F150W and F200W dropouts following the conventional dropout method to select candidate objects at z>11. In 14.2 arcmin^2, we have found 13 F150W dropouts and 8 F200W dropouts, all brighter than 27.5 mag in the band to the red side of the break signature. Most notably, some of these objects are as bright as ~24 mag, which corresponds to M_{UV}<-23 mag at z>11. As they are detected in multiple bands, these must be real objects. If the observed color decrements are due to the expected Lyman break, these objects should be at z>11.7 and z>15.4, respectively. The color diagnostics show that at least 11 F150W dropouts are far away from the usual contaminators encountered in dropout searches (red galaxies at much lower redshifts or brown dwarf stars), and therefore they are legitimate candidates. While the diagnostics of the F200W dropouts are less certain due to the limited number of passbands, at least one of them is not likely a known type of contaminant and the rest are consistent with either z>11 galaxies with evolved stellar populations or old galaxies at z ~ 3 to 8. If a significant fraction of our dropouts are indeed at z>11, we have to face the severe problem of explaining their high luminosities and number densities. Spectroscopic identifications of such objects are urgently needed.

We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST "Prime Extragalactic Areas for Reionization and Lensing Science" ("PEARLS") project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift proto-clusters; and the iconic backlit VV~191 galaxy system to map its dust attenuation. PEARLS also includes NIRISS spectra for one of the NEP fields and NIRSpec spectra of two high-redshift quasars. The main goal of PEARLS is to study the epoch of galaxy assembly, AGN growth, and First Light. Five fields, the JWST NEP Time-Domain Field (TDF), IRAC Dark Field (IDF), and three lensing clusters, will be observed in up to four epochs over a year. The cadence and sensitivity of the imaging data are ideally suited to find faint variable objects such as weak AGN, high-redshift supernovae, and cluster caustic transits. Both NEP fields have sightlines through our Galaxy, providing significant numbers of very faint brown dwarfs whose proper motions can be studied. Observations from the first spoke in the NEP TDF are public. This paper presents our first PEARLS observations, their NIRCam data reduction and analysis, our first object catalogs, the 0.9-4.5 $\mu$m galaxy counts and Integrated Galaxy Light. We assess the JWST sky brightness in 13 NIRCam filters, yielding our first constraints to diffuse light at 0.9-4.5 $\mu$m. PEARLS is designed to be of lasting benefit to the community.

This study investigates elastic deformation driven by the Hall drift in a magnetized neutron-star crust. Although the dynamic equilibrium initially holds without elastic displacement, the magnetic-field evolution changes the Lorentz force over a secular timescale, which inevitably causes the elastic deformation to settle in a new force balance. Accordingly, elastic energy is accumulated, and the crust is eventually fractured beyond a particular threshold. We assume that the magnetic field is axially symmetric, and we explicitly calculate the breakup time, maximum elastic energy stored in the crust, and spatial shear-stress distribution. For the barotropic equilibrium of a poloidal dipole field expelled from the interior core without a toroidal field, the breakup time corresponds to a few years for the magnetars with a magnetic field strength of $\sim 10^{15}$G; however, it exceeds 1 Myr for normal radio pulsars. The elastic energy stored in the crust before the fracture ranges from $10^{41}$ to $10^{45}$ erg, depending on the spatial-energy distribution. Generally, a large amount of energy is deposited in a deep crust. The energy released at fracture is typically $\sim 10^{41}$ erg when the rearrangement of elastic displacements occurs only in the fragile shallow crust. The amount of energy is comparable to the outburst energy on the magnetars.

We imaged the infrared dark clouds (IRDCs) G1.75-0.08 and G11.36+0.80 at 350 $\mu$m and 450 $\mu$m using the ArT\'eMiS bolometer. These data were used in conjunction with our previous 870 $\mu$m observations with the Large APEX BOlometer CAmera (LABOCA). The clumps in G11.36+0.80 were also observed in the N$_2$H$^+(1-0)$ transition with the IRAM 30-metre telescope. G1.75-0.08 was found to be composed of two cold ($\sim14.5$ K), massive (several $\sim10^3$ M$_{\odot}$) clumps that are projectively separated by $\sim3.7$ pc. Both clumps are 70 $\mu$m dark, but they do not appear to be bounded by self-gravity. The G1.75-0.08 filament was found to be subcritical by a factor of $\sim14$ with respect to its critical line mass. G11.36+0.80 was found to be moderately (by a factor of $\sim2$) supercritical and composed of four clumps. The dust temperatures of the clumps are $\sim13-15$ K, and their masses are in the range $\sim 232-633$ M$_{\odot}$. All the clumps are gravitationally bound. The projected, average separation of the clumps is $\sim1$ pc. A configuration that is observed in G1.75-0.08, namely two clumps at the ends of the filament, could be the result of gravitational focussing acting along the cloud. The two clumps fulfil the mass-radius threshold for high-mass star formation. Owing to the location of G1.75-0.08 near the Galactic centre ($\sim270$ pc), environmental effects such as a high level of turbulence, tidal forces, and shearing motions could affect the cloud dynamics. The observed clump separation in G11.36+0.80 can be understood in terms of a sausage instability. The G11.36+0.80 clumps do not lie above the mass-radius threshold for high-mass star formation. The substructure observed in one of the clumps in G11.36+0.80 suggests that the IRDC has fragmented in a hierarchical fashion. This conforms to the filamentary paradigm for Galactic star formation.

The X-IFU (X-ray Integral Field Unit) onboard the large ESA mission Athena (Advanced Telescope for High ENergy Astrophysics), planned to be launched in the mid 2030s, will be a cryogenic X-ray imaging spectrometer operating at 55 mK. It will provide unprecedented spatially resolved high-resolution spectroscopy (2.5 eV FWHM up to 7 keV) in the 0.2-12 keV energy range thanks to its array of TES (Transition Edge Sensors) microcalorimeters of more than 2k pixel. The detection chain of the instrument is developed by an international collaboration: the detector array by NASA/GSFC, the cold electronics by NIST, the cold amplifier by VTT, the WFEE (Warm Front-End Electronics) by APC, the DRE (Digital Readout Electronics) by IRAP and a focal plane assembly by SRON. To assess the operation of the complete readout chain of the X-IFU, a 50 mK test bench based on a kilo-pixel array of microcalorimeters from NASA/GSFC has been developed at IRAP in collaboration with CNES. Validation of the test bench has been performed with an intermediate detection chain entirely from NIST and Goddard. Next planned activities include the integration of DRE and WFEE prototypes in order to perform an end-to-end demonstration of a complete X-IFU detection chain.

In discussing open question in the field of massive stars, I consider their evolution from birth to death. After touching upon massive star formation, which may be bi-modal and not lead to a zero-age main sequence at the highest masses, I consider the consequences of massive stars being close to their Eddington limit. Then, when discussing the effects of a binary companion, I highlight the importance of massive Algols and contact binaries for understanding the consequences of mass transfer, and the role of binaries in forming Wolf-Rayet stars. Finally, a discussion on pair instability supernovae and of superluminous supernovae is provided.

We introduce an improved method for constraining the growth rate of structure with the galaxy overdensity and peculiar velocity power spectrum. This method reduces the modelling systematic error by accounting for the wide-angle effect and the zero-point calibration uncertainty during the modelling process. We also speed up the posterior sampling by around 30 times by first calculating the likelihood at a small number of fiducial points and then interpolating the likelihood values during MCMC sampling. We test the new method on mocks and we find it is able to recover the fiducial growth rate of structure. We applied our new method to the SDSS PV catalogue, which is the largest single peculiar velocity catalogue to date. Our constraint on the growth rate of structure is $f\sigma_8= 0.405_{-0.071}^{+0.076}$ (stat) $\pm$ 0.009 (sys) at the effective redshift of 0.073. Our constraint is consistent with a Planck 2018 cosmological model, $f\sigma_8$ = 0.448, within one standard deviation. Our improved methodology will enable similar analysis on future data, with even larger sample sizes and covering larger angular areas on the sky.

Orbital motion in binary and planetary systems is the main source of precise stellar and planetary mass measurements, and joint analysis of data from multiple observational methods can both lift degeneracies and improve precision. We set out to measure the masses of individual stars in binary systems using all the information brought by the Hipparcos and Gaia absolute astrometric missions. We present BINARYS, a tool which uses the Hipparcos and Gaia absolute astrometric data and combines it with relative astrometry and/or radial velocity measurements to determine the orbit of a binary system. It rigorously combines the Hipparcos and Gaia data (here EDR3), and it can use the Hipparcos Transit Data as needed for binaries where Hipparcos detect significant flux from the secondary component. It also support the case where Gaia resolved the system, giving an astrometric solution for both components. We determine model-independent individual masses for the first time for three systems: the two mature binaries Gl~494 ($M_1=0.584 \pm 0.003 M_{\odot}$ and $M_2=87 \pm 1 M_{\textrm{Jup}}$) and HIP~88745 ($M_1=0.96 \pm 0.02 M_{\odot}$ and $M_2= 0.60^{+ 0.02 }_{- 0.01 } M_{\odot}$), and the younger AB Dor member GJ~2060 ($M_1=0.60 ^{+ 0.06}_{- 0.05} M_{\odot}$ and $M_2=0.45 ^{+ 0.06}_{- 0.05}M_{\odot}$). The latter provides a rare test of evolutionary model predictions at young ages in the low stellar-mass range and sets a lower age limit of 100~Myr for the moving group.

We present UVIT/Astrosat UV photometry of the RSG population of the Small Cloud galaxy (SMC). As RSGs are extremely faint in the far-UV, these observations directly probe potential companion stars. From a sample of 861 SMC RSGs, we find 88 have detections at far-UV wavelengths: a clear signature of binarity. Stellar parameters are determined for both components, which allows us to study - for the first time - the mass-ratio (q) distribution of RSG binary systems. We find a flat mass-ratio distribution best describes the observations up to M{RSG}~15 Msun. We account for our main observing bias (i.e. the limiting magnitude of the UVIT survey) to determine the intrinsic RSG binary fraction of 18.8+/-1.5 %, for mass-ratios in the range 0.3 < q < 1.0 and orbital periods approximately in the range 3 < log P[ days] < 8.

We report the discovery of giant (50-100 kpc) [O II] emitting nebulae with the Multi-Unit Spectroscopic Explorer (MUSE) in the field of TXS 0206-048, a luminous quasar at z=1.13. An archival, down-the-barrel UV spectrum of the quasar shows absorption at velocities coincident with those of the extended nebulae, enabling new insights into inflows and outflows around the quasar host. One nebula exhibits a filamentary morphology extending over 120 kpc from the halo toward the quasar and intersecting with another nebula surrounding the quasar host with a radius of 50 kpc. The filamentary nebula has line-of-sight velocities >300 km/s from nearby galaxies but matches that of the nebula surrounding the quasar host where they intersect, consistent with filamentary accretion of cool inter- or circum-galactic medium or cooling hot halo gas. The kinematics of the nebulae surrounding the quasar host are unusual and complex, with one redshifted and one blue-shifted spiral-like structure. The nebular emission velocities at 5-10 kpc from the quasar match those of inflowing absorbing gas observed in a UV spectrum of the quasar. Together, the extended nebulae and associated redshifted absorption represent a compelling case of cool, filamentary gas accretion from halo scales into the extended interstellar medium and toward the nucleus of a massive quasar host galaxy at intermediate redshift. The inflow rate implied by the combination of emission and absorption constraints is orders-of-magnitude below levels required to sustain the quasar's radiative luminosity, indicating highly anisotropic or highly variable accretion.

We study here the phenomena of decayless kink oscillations in a system of active region (AR) coronal loops. Using high resolution observations from two different instruments, namely the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory, we follow these AR loops for an hour each on three consecutive days. Our results show significantly more resolved decayless waves in the higher-resolution EUI data compared with the AIA data. Furthermore, the same system of loops exhibits many of these decayless oscillations on Day-2, while on Day-3, we detect very few oscillations and on Day-1, we find none at all. Analysis of photospheric magnetic field data reveals that at most times, these loops were rooted in sunspots, where supergranular flows are generally absent. This suggests that supergranular flows, which are often invoked as drivers of decayless waves, are not necessarily driving such oscillations in our observations. Similarly, our findings also cast doubt on other possible drivers of these waves, such as a transient driver or mode conversion of longitudinal waves near the loop footpoints. In conclusion, through our analysis we find that none of the commonly suspected sources proposed to drive decayless oscillations in active region loops seems to be operating in this event and hence, the search for that elusive wave driver needs to continue.

The Sun exhibits a depletion in $^{17,18}$O relative to $^{16}$O by 6 % compared to the Earth and Moon$^{1}$. The origin of such a non-mass-dependent isotope fractionation has been extensively debated since the three-isotope-analysis$^{2}$ became available in 1970's. Self-shielding$^{3,4}$ of CO molecules against UV photons in the solar system's parent molecular cloud has been suggested as a source of the non-mass-dependent effect, in which a $^{17,18}$O-enriched oxygen was trapped by ice and selectively incorporated as water into planet-forming materials$^{5}$. The truth is that the Earth-Moon and other planetary objects deviate positively from the Sun by ~6 % in their isotopic compositions. A stunning exception is the magnetite/sulfide symplectite found in Acfer 094 meteorite, which shows 24 % enrichment in $^{17,18}$O relative to the Sun$^{6}$. Water does not explain the enrichment this high. Here we show that the SO and SO$_2$ molecules in the molecular cloud, ~106 % enriched in $^{17,18}$O relative to the Sun, evolved through the protoplanetary disk and planetesimal stages to become a sulfuric acid, 24 % enriched in $^{17,18}$O. The sulfuric acid provided a cryofluid environment in the planetesimal and by itself reacted with ferric iron to form an amorphous ferric-hydroxysulfate-hydrate, which eventually decomposed into the symplectite by shock. We indicate that the Acfer-094 symplectite and its progenitor, sulfuric acid, is strongly coupled with the material evolution in the solar system since the days of our molecular cloud.

Context. Carbonaceous nano-grains play a fundamental role in the physico-chemistry of the interstellar medium (ISM) and especially of photon-dominated regions (PDRs). Their properties vary with the local physical conditions and affect the local chemistry and dynamics. Aims. We aim to highlight the evolution of carbonaceous nano-grains in three different PDRs and propose a scenario of dust evolution as a response to the physical conditions. Methods. We used Spitzer/IRAC (3.6, 4.5, 5.8, and 8 $\mu$m) and Spitzer/MIPS (24 $\mu$m) together with Herschel/PACS (70 $\mu$m) to map dust emission in IC63 and the Orion Bar. To assess the dust properties, we modelled the dust emission in these regions using the radiative transfer code SOC together with the THEMIS dust model. Results. Regardless of the PDR, we find that nano-grains are depleted and that their minimum size is larger than in the diffuse ISM (DISM), which suggests that the mechanisms that lead nano-grains to be photo-destroyed are very efficient below a given critical size limit. The evolution of the nano-grain dust-to-gas mass ratio with both G0 and the effective temperature of the illuminating star indicates a competition between the nano-grain formation through the fragmentation of larger grains and nano-grain photo-destruction. We modelled dust collisions driven by radiative pressure with a classical 1D approach to show that this is a viable scenario for explaining nano-grain formation through fragmentation and, thus, the variations observed in nano-grain dust-to-gas mass ratios from one PDR to another. Conclusions. We find a broad variation in the nano-grain dust properties from one PDR to another, along with a general trend of nano-grain depletion in these regions. We propose a viable scenario of nano-grain formation through fragmentation of large grains due to radiative pressure-induced collisions.

Light curves of young star systems show photometric variability due to different kinematic, and physical processes. One of the main contributors to the photometric variability is the changing mass accretion rate, which regulates the interplay between the forming young star and the protoplanetary disk. We collected high-resolution spectroscopy in eight different epochs, as well as ground-based and space-borne multi-epoch optical and infrared photometry of WX Cha, an M0 binary system, with an almost edge-on disk (i = 87degrees) in the Chamaeleon I star-forming region. Spectroscopic observations cover 72 days, the ground-based optical monitoring covers 42 days while space-borne TESS photometry extends for 56 days. The multi-wavelength light curves exhibit quasi-periodic variability of 0.35 - 0.53 mag in the near-infrared, and of 1.3 mag in g band. We studied the variability of selected emission lines that trace the accretion, computed the accretion luminosity and the mass accretion rate using empirical relations and obtained values of the accretion luminosity between 1.6 and 3.2 Lsun and mass accretion rate between 3.31x10{-7} Msun/yr and 7.76x10^{-7} Msun/yr. Our results show that WX Cha is accreting at a rate larger than what is typical for T Tauri stars in the same star-forming region with the same stellar parameters. We theorize that this is due to the higher disk mass of WX Cha than what is usual for stars with similar stellar mass, and to the binary nature of the system. Daily changes in the accretion luminosity and in the extinction can explain the photometric variability.

Double detonations of sub-Chandrasekhar mass white dwarfs are a promising explosion scenario for Type Ia supernovae, whereby a detonation in a surface helium shell triggers a secondary detonation in a carbon-oxygen core. Recent work has shown that low mass helium shell models reproduce observations of normal SNe Ia. We present 3D radiative transfer simulations for a suite of 3D simulations of the double detonation explosion scenario for a range of shell and core masses. We find light curves broadly able to reproduce the faint end of the width-luminosity relation shown by SNe Ia, however, we find that all of our models show extremely red colours, not observed in normal SNe Ia. This includes our lowest mass helium shell model. We find clear Ti II absorption features in the model spectra, which would lead to classification as peculiar SNe Ia, as well as line blanketing in some lines of sight by singly ionised Cr and Fe-peak elements. Our radiative transfer simulations show that these explosion models remain promising to explain peculiar SNe Ia. Future full non-LTE simulations may improve the agreement of these explosion models with observations of normal SNe Ia.

The dependence of galaxy clustering on local density provides an effective method for extracting non-Gaussian information from galaxy surveys. The two-point correlation function (2PCF) provides a complete statistical description of a Gaussian density field. However, the late-time density field becomes non-Gaussian due to non-linear gravitational evolution and higher-order summary statistics are required to capture all of its cosmological information. Using a Fisher formalism based on halo catalogues from the Quijote simulations, we explore the possibility of retrieving this information using the density-split clustering (DS) method, which combines clustering statistics from regions of different environmental density. We show that DS provides more precise constraints on the parameters of the $\nu \Lambda$CDM model compared to the 2PCF, and we provide suggestions for where the extra information may come from. DS improves the constraints on the sum of neutrino masses by a factor of $8$ and by factors of 5, 3, 4, 6, and 6 for $\Omega_m$, $\Omega_b$, $h$, $n_s$, and $\sigma_8$, respectively. We compare DS statistics when the local density environment is estimated from the real or redshift-space positions of haloes. The inclusion of DS autocorrelation functions, in addition to the cross-correlation functions between DS environments and haloes, recovers most of the information that is lost when using the redshift-space halo positions to estimate the environment. We discuss the possibility of constructing simulation-based methods to model DS clustering statistics in different scenarios.

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

Triton is the largest satellite of Neptune and probably a Kuiper Belt Object that was captured by the planet. It has a tenuous nitrogen atmosphere similar to the one of Pluto and may be an ocean world. The Neptunian system has only been visited by Voyager 2 in 1989. Over the last few years, the demand for a new mission to the Ice Giants and their systems has increased so that a theoretical basis to prepare for such a mission is important. We aim to develop a photochemical model of Triton's atmosphere with an up-to-date chemical scheme, as previous photochemical models date back to the post-flyby years. This is done to better understand the mechanisms governing Triton's atmospheric chemistry and highlight the critical parameters having a significant impact on the atmospheric composition. We also study model uncertainties to find what chemical studies are necessary to improve the modeling of Triton's atmosphere. We adapted a model of Titan's atmosphere to Triton's conditions. We first used Titan's chemical scheme before updating it to better model Triton's atmosphere. Once the nominal results were obtained, we studied model uncertainties with a Monte-Carlo procedure. Then, we performed global sensitivity analyzes to identify the reactions responsible for model uncertainties. With the nominal results, we determined the composition of Triton's atmosphere and studied the main chemical processes. We highlighted key chemical reactions that are the most important for the overall chemistry. We also identified some key parameters having a significant impact on the results. Uncertainties are large for most of the main atmospheric species as the atmospheric temperature is very low. We identified key uncertainty reactions that have the largest impact on the results uncertainties. These reactions must be studied in priority in order to improve the significance of our results.

The evidence of a relation between the mass accretion rate and the disk mass is established for young, Class II pre-main sequence stars. This observational result opened an avenue to test theoretical models and constrain the initial conditions of the disk formation, fundamental in the understanding of the emergence of planetary systems. However, it is becoming clear that the planet formation starts even before the Class II stage, in disks around Class 0 and I protostars. We show for the first time evidence for a correlation between the mass accretion rate and the disk mass for a large sample of Class I young stars located in nearby (< 500 pc) star-forming regions. We fit our sample, finding that the Class I objects relation has a slope flatter than Class II stars, and have higher mass accretion rates and disk masses. The results are put in context of the disk evolution models.

We present the first empirical constraints on the turbulent velocity field of the diffuse circumgalactic medium around four luminous QSOs at $z\!\approx\!0.5$--1.1. Spatially extended nebulae of $\approx\!50$--100 physical kpc in diameter centered on the QSOs are revealed in [OII]$\lambda\lambda\,3727,3729$ and/or [OIII]$\lambda\,5008$ emission lines in integral field spectroscopic observations obtained using MUSE on the VLT. We measure the second- and third-order velocity structure functions (VSFs) over a range of scales, from $\lesssim\!5$ kpc to $\approx\!20$--50 kpc, to quantify the turbulent energy transfer between different scales in these nebulae. While no constraints on the energy injection and dissipation scales can be obtained from the current data, we show that robust constraints on the power-law slope of the VSFs can be determined after accounting for the effects of atmospheric seeing, spatial smoothing, and large-scale bulk flows. Out of the four QSO nebulae studied, one exhibits VSFs in spectacular agreement with the Kolmogorov law, expected for isotropic, homogeneous, and incompressible turbulent flows. The other three fields exhibit a shallower decline in the VSFs from large to small scales but with loose constraints, in part due to a limited dynamic range in the spatial scales in seeing-limited data. For the QSO nebula consistent with the Kolmogorov law, we determine a turbulence energy cascade rate of $\approx\!0.2$ cm$^{2}$ s$^{-3}$. We discuss the implication of the observed VSFs in the context of QSO feeding and feedback in the circumgalactic medium.

Magnetically confined mountains on accreting neutron stars are candidates for producing continuous gravitational waves. We formulate a magnetically confined mountain on a neutron star with strong multipole magnetic fields and obtain some sequences of numerical solutions. We find that the mass ellipticity of the mountain increases by one order of magnitude if the neutron star has strong multipole magnetic fields. As matter accretes on to the magnetic pole, the size of the mountain increases and the magnetic fields are buried. If the neutron star has a dipole magnetic field, the dipole magnetic field is buried and transformed into multipole components. By contrast, if the neutron star has both dipole and strong multipole magnetic fields, the multipole magnetic fields are buried and transformed into a negative dipole component. We also calculate magnetically confined mountains with toroidal magnetic fields and find that the ellipticity becomes slightly smaller when the mountain has toroidal magnetic fields. If the multipole magnetic fields are buried, they sustain the intense toroidal magnetic field near the stellar surface, and the ratio of the toroidal magnetic field to the poloidal magnetic field is close to 100. The hidden strong toroidal magnetic fields are sustained by the buried multipole magnetic fields.

First-order cosmological phase transitions in the early Universe source sound waves and, subsequently, a background of stochastic gravitational waves. Currently, predictions of these gravitational waves rely heavily on simulations of a Higgs field coupled to the plasma of the early Universe, the former providing the latent heat of the phase transition. Numerically, this is a rather demanding task since several length scales enter the dynamics. From smallest to largest, these are the thickness of the Higgs interface separating the different phases, the shell thickness of the sound waves, and the average bubble size. In this work, we present an approach to perform Higgsless simulations in three dimensions, producing fully nonlinear results, while at the same time removing the hierarchically smallest scale from the lattice. This significantly reduces the complexity of the problem and contributes to making our approach highly efficient. We provide spectra for the produced gravitational waves for various choices of wall velocity and strength of the phase transition, as well as introduce a fitting function for the spectral shape.

Using VIPERS, we estimated a set of SFR based on photometric and spectroscopic data. We used, as estimators, photometric bands from ultraviolet to mid-infrared, and spectral lines. Assuming a reference SFR obtained from the spectral energy distribution reconstructed with Code Investigating GALaxy Emission, we estimated the reliability of each band as an SFR tracer. We used GSWLC to trace the dependence of these SFR calibrators with redshift. The far and near UV, u-band and 24-$\mu$m bands, as well as $L_{TIR}$, are found to be good SFR tracers up to $z\sim0.9$ with a strong dependence on the attenuation prescription used for the bluest bands (scatter of SFR of 0.26, 0.14, 0.15, 0.23, and 0.24dex for VIPERS, and 0.25, 0.24, 0.09, 0.12, and 0.12dex for GSWLC). The 8-$\mu$m band provides only a rough estimate of the SFR as it depends on metallicity and polycyclic aromatic hydrocarbon properties (scatter of 0.23dex for VIPERS). We estimated the scatter of rest-frame luminosity estimations from CIGALE to be 0.26, 0.14, 0.12, 0.15, and 0.20dex for FUV, NUV, ugriz, K$_{\mathrm{s}}$, and 8-24$\mu$m-$L_{\mathrm{TIR}}$). At intermediate redshift, the H$\beta$ line is a reliable SFR tracer (scatter of 0.19dex) and the [OII] line gives an equally good estimation when the metallicity from the $R_{23}$ parameter is taken into account (0.17 for VIPERS and 0.20dex for GSWLC). A calibration based on [OIII] retrieves the SFR only when additional information such as the metallicity or the ionization parameter of galaxies are used (0.26 for VIPERS and 0.20dex for GSWLC), diminishing its usability as a direct SFR tracer. Based on rest-frame luminosities estimated with CIGALE, we propose our own set of calibrations from FUV, NUV, u-band, 8, 24$\mu$m, $L_{TIR}$, H$\beta$, [OII], and [OIII].

Astrophysical data, in the domains of time, involve a wide range of stellar variability phenomena, among them the magnetic activity of the order of a few hours until the signature of an extra-solar planet which can cover a scale of time of a few days until tens of years. Numerous instruments are being developed to detect Earth-sized exoplanets. Exoplanets with this dimension challenge scientific instrumentation and the field of research in the data processing. In this context, our study offers a powerful framework to explain dynamical properties as a function of timescale in light curves with the planetary signal. For that, we selected the stellar target Kepler-30 to test our methods and procedures. In this sense, we investigate the multifractal behavior of the Kepler-30 system composed of a sun-like star with a rotation period of ~16 days and three planets with masses between 2 Earth and 2.5 Jupiter masses. Furthermore, this system has an orbital period varying from 29 to 143 days and orbits almost coplanar. This system is highly interesting because starspots dynamics are strongly affected by the passing of a planet in front of the star. We used about 1600 days of high-precision photometry collected by the Kepler mission to investigate the quasi-periodic variation caused by the rotation of the star and the effect of spot evolution as a function of timescale. We applied indexes extract from multifractal analysis to model the flux rotational modulation induced by active regions. Our results that stellar flux variations in Kepler-30 star caused by rotational modulation can be replicated in detail with just four recent-known multifractal indexes. These indexes will greatly simplify spot modelling of current TESS and future PLATO data.

The SKA pulsar search pipeline will be used for real time detection of pulsars. Modern radio telescopes such as SKA will be generating petabytes of data in their full scale of operation. Hence experience-based and data-driven algorithms become indispensable for applications such as candidate detection. Here we describe our findings from testing a state of the art object detection algorithm called Mask R-CNN to detect candidate signatures in the SKA pulsar search pipeline. We have trained the Mask R-CNN model to detect candidate images. A custom annotation tool was developed to mark the regions of interest in large datasets efficiently. We have successfully demonstrated this algorithm by detecting candidate signatures on a simulation dataset. The paper presents details of this work with a highlight on the future prospects.

We compute the effect of rigid rotation on the non-relativistic bound states. The energy levels of the bound states increase with the angular velocity of rotation until at certain value of the angular velocity they are completely pushed out into the continuum which corresponds to dissociation of the bound states. When the angular velocity exceeds the critical value at which the ground state disappears into the continuum, no bound state is possible. This effect should have important consequences for the phenomenology of the quark-gluon plasma. One of the ways to study it experimentally is to observe the electromagnetic radiation emitted by a rotating bound state. We compute the corresponding intensity of electromagnetic radiation and show that it strongly depends on the angular velocity of rotation.

Recently, as a generalization of general relativity, a gravity theory has been proposed in which gravitational field equations are described by the Cotton tensor. That theory allows an additional contribution to the gravitational potential of a point mass that rises linearly with radius as $\Phi = -GM/r + \gamma r/2$, where $G$ is the Newton constant. The coefficients $M$ and $\gamma$ are the constants of integration and should be determined individually for each physical system. When applied to galaxies, the coefficient $\gamma$, which has the dimension of acceleration, should be determined for each galaxy. This is the same as having to determine the mass $M$ for each galaxy. If $\gamma$ is small enough, the linear potential term is negligible at short distances, but can become significant at large distances. In fact, it may contribute to the extragalactic systems. In this paper, we derive the effective field equation for Cotton gravity applicable to extragalactic systems. We then use the effective field equation to numerically compute the gravitational potential of a sample of 84 rotating galaxies. The 84 galaxies span a wide range, from stellar disk-dominated spirals to gas-dominated dwarf galaxies. We do not assume the radial density profile of the stellar disk, bulge, or gas; we use only the observed data. We find that the rotation curves of 84 galaxies can be explained by the observed distribution of baryons. This is due to the flexibility of Cotton gravity to allow the integration constant $\gamma$ for each galaxy. In the context of Cotton gravity, "dark matter" is in some sense automatically included as a curvature of spacetime. Consequently, even galaxies that have been assumed to be dominated by dark matter do not need dark matter.

This is the final report from the Snowmass 2021 Neutrino Frontier Topical Group on Neutrinos from Natural Sources. It covers a broad range of neutrino sources, from low-energy neutrinos from the early universe to ultra high-energy sources. We divide this report by source, and discuss the motivations for pursuing searches in each case, the current state of the field, and the prospects for future theoretical and experimental developments. We consider neutrinos produced in the early universe; solar neutrinos; geoneutrinos; supernova neutrinos, including the diffuse supernova neutrino background (DSNB); neutrinos produced in the atmosphere; and high-energy astrophysical neutrinos.

The Laser Interferometer Space Antenna (LISA) mission, scheduled for launch in the mid-2030s, is a gravitational wave observatory in space designed to detect sources emitting in the millihertz band. LISA is an ESA flagship mission, currently entering the Phase B development phase. It is expected to help us improve our understanding about our Universe by measuring gravitational wave sources of different types, with some of the sources being at very high redshifts $z\sim 20$. On the 23rd of February 2022 we organized the 1$^\mathrm{st}$ {\it LISA in Greece Workshop}. This workshop aimed to inform the Greek scientific and tech industry community about the possibilities of participating in LISA science and LISA mission, with the support of the Hellenic Space Center (HSC). In this white paper, we summarize the outcome of the workshop, the most important aspect of it being the inclusion of $15$ Greek researchers to the LISA Consortium, raising our total number to $22$. At the same time, we present a road-map with the future steps and actions of the Greek Gravitational Wave community with respect to the future LISA mission.

Dark matter that interacts strongly with baryons can avoid the stringent dark matter direct detection constraints, because, like baryons, they are likely to be absorbed when traversing the rocks, leading to a suppressed flux in deep underground labs. Such strongly interacting dark matter, however, can be probed by dark matter experiments or other experiments operated on the ground level or in the atmosphere. In this paper, we carry out systematic analysis of two of these experiments, XQC and CSR, to compute the experimental constraints on the strongly interacting dark matter, where three scenarios are considered: (1) spin-independent and spin-dependent interactions; (2) different velocity dependent cross sections; (3) different dark matter mass fractions. Some of the scenarios are first analyzed in the literature. We find that the XQC exclusion region has some non-trivial dependencies on the various parameters and the limits in the spin-dependent case is quite different from the spin-independent case. We also identify a peculiar region in the parameter space where the XQC constraint disappears, in the case where the interaction cross section is proportional to the square of the velocity. We further compare our XQC and CSR limits to other experimental constraints, and find that a large parameter space is allowed by various experiments if the dark matter mass fraction is sufficiently small, $f_\chi\lesssim 10^{-4}$.

Fluctuation terms and higher moments of a quantum state imply corrections to the classical equations of motion that may have implications in early-universe cosmology, for instance in the state-dependent form of effective potentials. In addition, space-time properties are relevant in cosmology, in particular when combined with quantum corrections required to maintain general covariance in a consistent way. Here, an extension of previous investigations of static quasiclassical space-time models to dynamical ones is presented, describing the evolution of 1-dimensional space as in the classical Lemaitre--Tolman--Bondi models. The corresponding spatial metric has two independent components, both of which are in general subject to quantum fluctuations. The main result is that individual moments from both components are indeed required for general covariance to be maintained at a semiclassical level, while quantum correlations between the components are less relevant.