Papers for Monday, Feb 27 2023

Galaxy morphology is a key parameter in galaxy evolution studies. The enormous number of galaxies which current and future surveys will observe demand of automated methods for morphological classification. Supervised learning techniques have been successfully used for the morphological classification of galaxies from different datasets, including Sloan Digital Sky Survey (SDSS), Mapping Galaxies with Apache Point Observatory (MaNGA) or Dark Energy Survey (DES). With these proceedings, we release the morphological catalogue for a sample of 670,000 SDSS galaxies based on the deep learning models trained on SDSS RGB images with morphological labels from human-based classification catalogues. The released catalogue includes binary classifications (early-type versus late-type, elliptical versus lenticular, identification of edge-on and barred galaxies) plus a T-Type. The classifications also include k-fold based uncertainties. This is, as of today, the largest catalogue including a T-Type classification. As an example of the scientific potential of this classification, we show how the location of the galaxies in the star formation - stellar mass plane (SFR-M$^{*}$) depends on morphology. This is the first time the SFR-M$^{*}$ relation is combined with T-Type information for such a large sample of galaxies.

We combine an unprecedented MaNGA sample of over 3,000 passive galaxies in the stellar mass range 10^{9}-10^{12} Msun with the Sloan Digital Sky Survey group catalog by Tinker to quantify how central and satellite formation, quantified by radial profiles in stellar age, [Fe/H], and [Mg/Fe], depends on the stellar mass of the galaxy (M*) and the mass of the host halo (Mh). After controlling for M* and Mh, the stacked spectra of centrals and satellites beyond the effective radius (r_e) show small, yet significant differences in multiple spectral features at the 1% level. According to spectral fitting with the code alf, a primary driver of these differences appears to be [Mg/Fe] variations, suggesting that stellar populations in the outskirts of satellites formed more rapidly than the outer populations of centrals. To probe the physical mechanisms that may be responsible for this signal, we examined how satellite stellar populations depend on Mh. We find that satellites in high-Mh halos show older stellar ages, lower [Fe/H], and higher [Mg/Fe] compared to satellites in low-Mh halos, especially for M*=10^{9.5}-10^{10.5} Msun. These signals lend support to environmentally driven processes that quench satellite galaxies, although variations in the merger histories of central and satellite galaxies also emerge as a viable explanation.

Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model - supported and benchmarked with large-scale three-dimensional particle-in-cell simulations - that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factor $\gamma\gtrsim 3\sigma$ (here, $\sigma$ is the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is $t_{\rm acc}\sim \gamma \,\eta_{\rm rec}^{-1}\omega_{\rm c}^{-1}\simeq 20\,\gamma\,\omega_{\rm c}^{-1}$, where $\eta_{\rm rec}\simeq0.06$ is the inflow speed in units of the speed of light and $\omega_{\rm c}=eB_0/mc$ is the gyrofrequency in the upstream magnetic field. They leave the region of active energization after $t_{\rm esc}$, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measure $t_{\rm esc}$ in our simulations and find that $t_{\rm esc}\sim t_{\rm acc}$ for $\sigma\gtrsim {\rm few}$. This leads to a universal (i.e., $\sigma$-independent) power-law spectrum $dN_{\rm free}/d\gamma\propto \gamma^{-1}$ for the particles undergoing active acceleration, and $dN/d\gamma\propto \gamma^{-2}$ for the overall particle population. Our results help shedding light on the ubiquitous presence of power-law particle and photon spectra in astrophysical non-thermal sources.

We take advantage of the NIRCam photometric observations available as part of the Cosmic Evolution Early Release Science survey (CEERS) to identify and analyse very red sources in an effort to discover very dusty star forming galaxies. We select red galaxies as objects with a S/N>3 at 4.4 $\mu$m and a S/N<2 in all JWST and HST filters at $\lambda\leq2\mu$m, which corresponds to [F200W]-[F444W]>1.2 considering CEERS depths. This selection is ideal to identify very dusty (Av>1 mag) galaxies with stellar masses between $10^6$ to $10^{10}\, \rm M_{\odot}$ at z<5, more massive dusty galaxies at z=5-18 and galaxies at z>18 due to the Lyman absorption, independently of their dust extinction. Our sample of F200W-dropouts contains no strong candidates at z>6.5, instead it consists almost completely (~82%) of z<2 low-mass galaxies, with a median stellar mass of $10^{7.3} \rm M_{\odot}$. These galaxies show an exceptional dust extinction with median value of Av=4.9 mag, completely unexpected given their low stellar mass. The remaining galaxies, which are at z<6.5, show similar large dust extinction (Av>1), but they are generally more massive $>10^{7.5}\rm M_{\odot}$.

Cosmological analyses of second-order weak lensing statistics require precise and accurate covariance estimates. One significant but sometimes misunderstood component of the covariance is the so-called super-sample covariance (SSC). The SSC is regularly defined as the covariance part capturing all impact of modes outside of a survey area. However, we show here that this intuition is incorrect for real-space statistics. We derive the covariance of second-order shear statistics from first principles. For this, we use an estimator in real space without relying on an estimator for the power spectrum. All parts of the covariance, not only the SSC, depend on the power- and trispectrum at all modes, including those larger than the survey. The defining feature of the SSC is not its dependence on ``super-survey modes'' but its behaviour under the `large-field approximation', i.e., the limiting case of a very large survey area. While the non-SSC parts of the covariance scale with the inverse survey area in this limit, the SSC completely vanishes. We also show that it is generally impossible to transform an estimate for the power spectrum covariance to the covariance of a real-space statistic. Such a transformation is only possible in the limiting case of the `large-field approximation'. Additionally, we find that the total covariance of a real-space statistic can be completely estimated using correlation functions measured only inside the survey area. Consequently, estimating covariances of real-space statistics, in principle, does not require information outside the survey boundaries.

AM CVn systems are binary star systems with orbital periods less than 70 minutes in which a white dwarf accretes matter from a companion star, which must be either a stripped helium burning star, or a white dwarf of lower mass than the accretor. Here, we present the discoveries of two of these systems in which there is mass transfer from the lighter white dwarf or helium star onto a strongly magnetized heavier white dwarf. These represent the first clear example of magnetized accretion in ultracompact binaries. These systems, along with similar systems that are slightly more widely separated, and that have not started to transfer mass yet, are expected to be the primary source of gravitational waves to be detected by space-based gravitational wave observatories. The presence of strong magnetic fields can substantially affect both the evolution of the binaries, and also the particular wave forms of the gravitational waves themselves, and understanding these magnetic effects is vital for understanding what to expect from the Laser Interferometer Space Antenna.

The anticorrelated distributions of temperature and density of protons are a well-known property of the solar wind. Nevertheless, it is unclear till now if they are formed by some kind of the universal physical mechanism? Unfortunately, a straightforward comparison of the characteristic relaxation times for the temperature and density, on the one hand, and pressure, on the other hand, encounters the problem of inapplicability of the hydrodynamical approach in the situation when the free-path length of the protons is considerably greater than the spatial scale of the structures under consideration. To resolve this problem, some kinds of the MHD turbulence - reducing the effective free paths - are usually assumed. In the present paper, we use an alternative approach based on the electrostatic (Langmuir) turbulence, described by the mathematical formalism of the spin-type Hamiltonians, which was actively discussed in the recent time in the literature on statistical physics. As follows from the corresponding calculations, formation of the anticorrelated distributions of temperature and density is a universal property of the strongly-nonequilibrium plasmas governed by the spin-type Hamiltonians when they gradually approach the thermodynamic equilibrium. So, just this phenomenon could be responsible for the anticorrelations observed in the solar wind.

Context. Magnetic fields play a fundamental role in the dynamical evolution of protoplanetary disks, in particular via magnetically induced disk winds. The magnetic field structure at the disk surface is crucial for driving the disk winds; however, it is still poorly understood observationally. Aims. We explore a new method to probe the magnetic field structure at the disk surface using near-infrared (NIR) circular polarization. Near-infrared circular polarization arises when unpolarized stellar light is scattered by magnetically aligned grains at the disk surface. In this study, we aim to clarify to what extent the observed circular polarization pattern can be used to diagnose the magnetic field structure. Methods. We first calculated light scattering properties of aligned spheroids, and the results were then used to create expected observational images of the degree of circular polarization at a NIR wavelength. Results. Magnetically aligned grains can produce circular polarization, particularly when the field configuration deviates from a purely toroidal field. We find that disk azimuthal dependence of the degree of circular polarization tends to exhibit a double peaked profile when the field structure is favorable for driving disk winds by centrifugal force. We also find that even if the disk is spatially unresolved, a net circular polarization can possibly be nonzero. We also show that the amplitude of circular polarization is strongly dependent on grain composition and axis ratio. Conclusions. Our results suggest that circular polarization observations would be useful to study the magnetic field structure and dust properties at the disk surface.

Recent detections of potentially habitable exoplanets around sunlike stars demand increased exploration of the physical conditions that can sustain life, by whatever methods available. Insight into these conditions can be gained by considering the multiverse hypothesis; in a multiverse setting, the probability of living in our universe depends on assumptions made about the factors affecting habitability. Various proposed habitability criteria can be systematically considered to rate each on the basis of their compatibility with the multiverse, generating predictions which can both guide expectations for life's occurrence and test the multiverse hypothesis. Here, we evaluate several aspects of planetary habitability, and show that the multiverse does indeed induce strong preferences among them. We find that the notion that a large moon is necessary for habitability is untenable in the multiverse scenario, as in the majority of parameter space, moons are not necessary to maintain stable obliquity. Further, we consider various proposed mechanisms for water delivery to the early Earth, including delivery from asteroids, both during giant planet formation and a grand tack, delivery from comets, and oxidation of a primary atmosphere by a magma ocean. We find that, depending on assumptions for how habitability depends on water content, some of these proposed mechanisms are disfavored in the multiverse scenario by Bayes factors of up to several hundred.

We construct accurate emulators for the projected and redshift space galaxy correlation functions and excess surface density as measured by galaxy-galaxy lensing, based on Halo Occupation Distribution (HOD) modeling. Using the complete Mira-Titan suite of 111 $N$-body simulations, our emulators vary over eight cosmological parameters and include the effects of neutrino mass and dynamical dark energy. We demonstrate that our emulators are sufficiently accurate for the analysis of the BOSS DR12 CMASS galaxy sample over the range 0.5 < r < 50 Mpc/h. Furthermore, we show that our emulators are capable of recovering unbiased cosmological constraints from realistic mock catalogs over the same range. Our mock catalog tests show the efficacy of combining small scale galaxy-galaxy lensing with redshift space clustering and that we can constrain the growth rate and \sigma_8 to 7% and 4.5% respectively for a CMASS-like sample using only the measurements covered by our emulator. With the inclusion of a CMB prior on H_0, this reduces to a 2% measurement on the growth rate.

Globular clusters are among the oldest stellar populations in the Milky Way; consequently, they also host some of the oldest known stellar-mass black holes, providing insight into black hole formation and evolution in the early ($z\gtrsim 2$) Universe. Recent observations of supermassive black holes in elliptical galaxies have been invoked to suggest the possibility of a cosmological coupling between astrophysical black holes and the surrounding expanding Universe, offering a mechanism for black holes to grow over cosmic time, and potentially explaining the origin of dark energy. In this paper, I show that the mass functions of the two radial velocity black hole candidates in NGC 3201 place strong constraints on the cosmologically-coupled growth of black holes. In particular, the amount of coupling required to explain the origin of dark energy would either require both NGC 3201 black holes to be nearly face on (a configuration with probability of at most $10^{-4}$) or one of the BHs would need to have formed with a mass below that of the most massive neutron stars ($2.2M_{\odot}$). This emphasizes that these and other detached black hole-star binaries can serve not only as laboratories for compact object and binary astrophysics, but as constraints on the long-term evolution of astrophysical black holes.

Over-luminous type Ia supernovae (SNe Ia) show peculiar observational features, for which an explosion of a super-massive white dwarf (WD) beyond the classical Chandrasekhar-limiting mass has been suggested, largely based on their high luminosities and slow light-curve evolution. However, their observational features are diverse, with a few extremely peculiar features whose origins have not been clarified; strong and persisting C II lines, late-time accelerated luminosity decline and red spectra, and a sub-day time-scale initial flash clearly identified so far at least for three over-luminous SNe Ia. In the present work, we suggest a scenario that provides a unified solution to these peculiarities, through hydrodynamic and radiation transfer simulations together with analytical considerations; a C+O-rich envelope (~0.01 - 0.1 Msun) attached to an exploding WD. Strong C II lines are created within the shocked envelope. Dust formation is possible in the late phase, providing a sufficient optical depth thereafter. The range of the envelope mass considered here predicts an initial flash with time-scale of ~0.5 - 3 days. The scenario thus can explain some of the key diverse observational properties by a different amount of the envelope, but additional factors are also required; we argue that the envelope is distributed in a disc-like structure, and also the ejecta properties, e.g., the mass of the WD, plays a key role. Within the context of the hypothesized super-Chandrasekhar-mass WD scenario, we speculatively suggest a progenitor WD evolution including a spin-up accretion phase followed by a spin-down mass-ejection phase.

In the mixed-mode asteroseismology of subgiants and red giants, the coupling between the p- and g-mode cavities must be understood well in order to derive localised estimates of interior rotation from measurements of mode multiplet rotational splittings. There exist now two different descriptions of this coupling: one based on an asymptotic quantisation condition, and the other arising from coupling matrices associated with "acoustic molecular orbitals". We examine the analytic properties of both, and derive closed-form expressions for various quantities -- such as the period-stretching function $\tau$ -- which previously had to be solved for numerically. Using these, we reconcile both formulations for the first time, deriving relations by which quantities in each formulation may be translated to and interpreted within the other. This yields an information criterion for whether a given configuration of mixed modes meaningfully constrains the parameters of the asymptotic construction, which is likely not satisfied by the majority of first-ascent red giant stars in our observational sample. Since this construction has been already used to make rotational measurements of such red giants, we examine -- through a hare-and-hounds exercise -- whether, and how, such limitations affect existing measurements. While averaged estimates of core rotation seem fairly robust, template-matching using the asymptotic construction has difficulty reliably assigning rotational splittings to individual multiplets, or estimating mixing fractions $\zeta$ of the most p-dominated mixed modes, where such estimates are most needed. We finally discuss implications for extending the two-zone model of radial differential rotation, e.g. via rotational inversions, with these methods.

No existence of convincing signals has become a ``new normal'' in the Dark Matter (DM) direct detection community for decades. Among other possibilities, the ``new normal'' might indicate that in addition to the traditional NR (Nuclear Recoil) events, ER (Electron Recoil) ones could also result from DM interactions, which have been considered as backgrounds historically. Further, we argue that ER and NR-like DM signals could co-exist in a detector's same dataset. So in total, there would be three scenarios we can search for DM candidate signals: (i) NR excess only, (ii) ER excess only, and (iii) both ER and NR excess. To effectively identify possible DM signals under the three scenarios, a DM detector should (a) have the minimum ER and NR backgrounds and (b) be capable of discriminating ER events from NR ones. Accordingly, we introduce the newly established project, ALETHEIA (A Liquid hElium Time projection cHambEr In dArk matter), which implements liquid helium TPCs (Time Projection Chamber) to hunt for low-mass ($\sim $100s MeV/c$^2$ - 10 GeV/c$^2$) DM. Presumably, the LHe TPC technology would have the minimum intrinsic backgrounds and strong ER/NR discrimination, therefore, be capable of identifying any kind of DM induced excess if it exists at all.

We propose a new method to constrain the size of the dusty torus in broad-line active galactic nuclei (AGNs) using optical and mid-infrared (MIR) ensemble structure functions (SFs). Because of the geometric dilution of the torus, the mid-infrared response to optical continuum variations has suppressed variability with respect to the optical that depends on the geometry (e.g., size, orientation and opening angle) of the torus. More extended tori have steeper MIR SFs with respect to the optical SF. We demonstrate the feasibility of this SF approach using simulated AGN light curves and a geometric torus model. While it is difficult to use SFs to constrain the orientation and opening angle of the torus due to model degeneracies and insensitivity of the SF on these parameters, the size of the torus can be well determined using this method. Applying this method to the ensemble SFs measured for 587 SDSS quasars with both optical and MIR light curves, we measure a best-fit torus $R-L$ relation of ${\rm log}\,R_{\rm eff}\, {(\rm pc)} = 0.49_{-0.03}^{+0.04} \times {\rm log}\,(L_{\rm bol}/10^{46}\,\rm erg\,s^{-1}) -0.40_{-0.01}^{+0.01}$, which is in good agreement with dust reverberation mapping measurements in AGNs. Compared with the dust reverberation mapping technique, the SF method is much less demanding in data quality and can be applied to any optical+MIR light curves for which a lag measurement may not be possible, as long as the variability process and torus structure are stationary. While this SF method does not extract and utilize all information contained in the light curves (i.e., the transfer function), it provides an intuitive interpretation for the observed trends of AGN MIR SFs compared with optical SFs.

The isotopic signature of atmospheric carbon offers a unique tracer for the history of the Martian atmosphere and the origin of organic matter on Mars. Photolysis of CO$_{2}$ is known to induce strong isotopic fractionation of carbon between CO$_{2}$ and CO. However, its effect on the carbon isotopic compositions in the Martian atmosphere remains uncertain. Here we develop a 1-D photochemical model considering isotopic fractionation via photolysis of CO$_{2}$ to estimate the vertical profiles of the carbon isotopic compositions of CO and CO$_{2}$ in the Martian atmosphere. We find that CO is depleted in $^{13}$C compared with CO$_{2}$ at each altitude due to the fractionation via CO$_{2}$ photolysis: the minimum value of $\delta ^{13}$C in CO is about $-170$ per mil under the standard eddy diffusion setting. This result supports the hypothesis that fractionated atmospheric CO is responsible for the production of the $^{13}$C-depleted organic carbon in Martian sediments detected by Curiosity Rover through the conversion of CO into organic materials and their deposition on the surface. The photolysis and transport-induced fractionation of CO we report here leads to a $\sim 15$ % decrease in the amount of inferred atmospheric loss when combined with the present-day fractionation of the atmosphere and previous studies of carbon escape to space. The fractionated isotopic composition of CO in the Martian atmosphere may be observed by ExoMars Trace Gas Orbiter (TGO) and ground-based telescopes, and escaping ion species produced by the fractionated carbon-bearing species may be detected by Martian Moons eXploration (MMX) in the future.

Solar flares fall into two types with eruptive ones associated with coronal mass ejection (CME) and confined ones without CME. To explore whether there are pre-flare conditions in terms of magnetic energy and helicity that can effectively determine the types of flares, here we analyzed a suite of related parameters of the reconstructed pre-flare coronal magnetic field of major solar flares, either eruptive or confined, from 2011 to 2017 near the solar disk center. The investigated parameters include the extensive-type quantities such as the total magnetic energy $E_T$, the potential energy $E_P$, the free energy $E_F$, the relative helicity $H_R$, and the non-potential helicity $H_J$, as well as the intensive-type indices $E_F/E_P$, $|H_J/H_R|$, $|H_R/\phi^{\prime2}|$ and $|H_J/\phi^{\prime2}|$, where $\phi^{\prime}$ is half of the total unsigned magnetic flux. We have the following key findings: (1) None of the extensive parameters can effectively distinguish the eruptive and confined potential of the pre-flare coronal fields, though the confined events have averagely larger values; (2) All the intensive parameters have significantly larger average and median values for eruptive flares than the confined events, which indicates that the field for eruptive flares have overall higher degree of non-potentiality and complexity than that of the confined flares; (3) The energy ratio $E_F/E_P$ and the normalized non-potential helicity $|H_J/\phi^{\prime2}|$, which are strongly correlated with each other, have among the highest capability of distinguishing the fields that possibly produce a major eruptive or confined flare, as over 75\% of all the events are successfully discriminated between eruptive and confined flares by using critical values of $E_F/E_P\ge0.27$ and $|H_J/\phi^{\prime2}|\ge0.009$.

We have collected 10025 foreground-background quasar pairs with projected distances $d_p<500$ kpc from the large quasar catalog of the SDSS DR16Q. We investigate the properties of the Mg II absorption lines with $W_r>0.15$ \AA\ around foreground quasars, including both the LOS (line-of-sights of foreground quasars) and transverse (TRA, perpendicular to the LOS) absorptions. Both the equivalent width (the correlation coefficient $\rho=-0.915$ and the probability $P < 10^{-4}$ of no correlation) and incident rate ($\rho=-0.964$ and $P < 10^{-6}$) of TRA \Mgii\ absorption lines are obviously anti-correlated with projected distance. The incident rate of TRA \Mgii\ absorption lines is obviously ($>4\sigma$) greater than that of LOS \Mgii\ absorption lines at projected distances $d_p<200$ kpc, while the TRA and LOS \Mgii\ both have similar ($<3\sigma$) incident rates at scales $d_p>200$ kpc. The anisotropic radiation from quasars would be the most possible interpretation for the anisotropic absorption around quasars. This could also indicate that the quasar radiation is not obviously impacting the gas halos of quasars at scales $d_p>200$ kpc.

Hot subdwarfs in close binaries with either M dwarf, brown dwarf or white dwarf companions show unique light variations. In hot subdwarf binaries with M dwarf or brown dwarf companions we can observe the so-called reflection effect and in hot subdwarfs with close white dwarf companions ellipsoidal modulation and/or Doppler beaming. Aims. The analysis of these light variations can be used to derive the mass and radius of the companion and hence determine its nature. Thereby we assume the most probable sdB mass and the radius of the sdB derived by the fit of the spectral energy distribution and the Gaia parallax. In the high signal-to-noise space-based light curves from the Transiting Exoplanet Survey Satellite and the K2 space mission, several reflection effect binaries and ellipsoidal modulation binaries have been observed with much better quality than possible for ground-based observations. The high quality of the light curves allowed us to analyse a large sample of sdB binaries with M dwarf or white dwarf companions using lcurve. For the first time we can constrain the absolute parameters of 19 companions of reflection effect systems covering periods from 2.5 to 19 hours and companion masses from the hydrogen burning limit to early M dwarfs. Moreover, we could determine the mass of eight white dwarf companion to hot subdwarf binaries showing ellipsoidal modulations, covering a so far unexplored period range from 7 to 19 hours. The derived masses of the white dwarf companions show that all but two of the white dwarf companions are most likely helium-core white dwarfs. Combining our results with previously measured rotation velocities allowed us to derive the rotation period of seven sdBs in short-period binaries. In four of those systems the rotation period of the sdB agrees with a tidally locked orbit, in the other three systems the sdB rotates significantly slower.

Methane is thought to have been an important greenhouse gas during the Archean, although its potential warming has been found to be limited at high concentrations due to its high shortwave absorption. We use the Met Office Unified Model, a general circulation model, to further explore the climatic effect of different Archean methane concentrations. Surface warming peaks at a pressure ratio CH$_4$:CO$_2$ of approximately 0.1, reaching a maximum of up to 7 K before significant cooling above this ratio. Equator-to-pole temperature differences also tend to increase up to pCH$_4$ $\leq$300 Pa, which is driven by a difference in radiative forcing at the equator and poles by methane and a reduction in the latitudinal extend of the Hadley circulation. 3D models are important to fully capture the cooling effect of methane, due to these impacts of the circulation.

The link between the fate of most massive stars and the resulting supernova (SN) explosion is still debated, also because of the ambiguity of the light curve powering mechanisms. When they explode as SNe, the light curve luminosity is typically sustained by a central engine (radioactive decay, magnetar spin-down, or fallback accretion). However, since massive stars eject considerable amounts of material during their evolution, there can be significant contribution from interaction with the previously-ejected circumstellar medium (CSM). Reconstructing the progenitor configuration at the time of explosion requires a detailed analysis of the long-term photometric and spectroscopic evolution of the related transient. In this paper we present the results of our follow-up campaign of SN 2020faa that, because of the high luminosity and peculiar slow light curve, is candidate for having a massive progenitor. We present the spectro-photometric dataset and investigate different options to explain the unusual observed properties. We compute the bolometric luminosity of the supernova and the evolution of its temperature, radius, and expansion velocity. We also fit the observed light curve with a multi-component model to infer information on the progenitor and the explosion mechanism. Reasonable parameters are inferred for SN 2020faa with a magnetar of energy Ep = 1.5(+0.5,-0.2) x 10^50 erg and spin-down time t(spin)=15+/-1 d, a shell mass M(shell) = 2.4(+0.5,-0.4) Msun and kinetic energy Ekin(shell) = 0.9(+0.5,-0.3) x 10^51 erg, and a core with M(core) = 21.5(+1.4,-0.7) Msun and Ekin(core) = 3.9(+0.1,-0.4) x 10^51 erg. In addition, we need an extra source to power the luminosity of the second peak. We find that hidden interaction with either a CSM disc or with delayed, choked jets is a viable mechanism to supply the required energy.

Active Galactic Nuclei (AGN) are remarkable astronomical sources emitting over the whole electromagnetic spectrum, with different bands providing unique windows on distinct sub-structures and their related physics. AGN come in a large number of types only partially related to intrinsic differences. I highlight here the most important AGN classes, namely jetted and non-jetted, radiatively efficient and inefficient, and face-on and edge-on, the source types selected by different bands together with the most important selection effects and biases, and the underlying emission processes, emphasising the gamma-ray band. I then conclude with a look at some open issues in AGN research and at the main new astronomical facilities, which will provide us with new data to tackle them.

Phobos is the target of the return sample mission Martian Moons eXploration by JAXA that will analyze in great details the physical and compositional properties of the satellite from orbit, from the surface and in terrestrial laboratories, giving clues about its formation. Some models propose that Phobos and Deimos were formed after a giant impact giving rise to an extended debris disk. Assuming that Phobos formed from a cascade of disruptions and re-accretions of several parent bodies in this disk, and that they are all characterized by a low material cohesion, Hesselbrock & Milton (2017) have showed that a recycling process may happen during the assembling of Phobos, by which Phobos' parents are destroyed into a Roche-interior ring and reaccreted several times. In the current paper we explore in details the recycling model, and pay particular attention to the characteristics of the disk using 1D models of disk/satellite interactions. In agreement with previous studies we confirm that, if Phobos' parents bodies are gravitational aggregates (rubble piles), then the recycling process does occur. However, Phobos should be accompanied today by a Roche-interior ring. Furthermore, the characteristics of the ring are not reconcilable with today`s observations of Mars' environment, which put stringent constraints on the existence of a ring around Mars. The recycling mechanism may or may not have occurred at the Roche limit for an old moon population, depending on their internal cohesion. However, the Phobos we see today cannot be the outcome of such a recycling process.

We present an adapted version of the code HII-CHI-Mistry-UV (P\'erez-Montero & Amor\'in 2017) to derive chemical abundances from emission lines in the ultraviolet, for use in narrow line regions (NLR) of Active Galactic Nuclei (AGN). We evaluate different ultraviolet emission line ratios and how different assumptions about the models, including the presence of dust grains, the shape of the incident spectral energy distribution, or the thickness of the gas envelope around the central source, may affect the final estimates as a function of the set of emission lines used. We compare our results with other published recipes for deriving abundances using the same emission lines and show that deriving the carbon-to-oxygen abundance ratio using CIII] $\lambda$ 1909 \r{A} and OIII] $\lambda$ 1665 \r{A} emission lines is a robust indicator of the metal content in AGN that is nearly independent of the model assumptions, similar to the case of star-forming regions. Moreover, we show that a prior determination of C/O allows for a much more precise determination of the total oxygen abundance using carbon UV lines, as opposed to assuming an arbitrary relationship between O/H and C/O, which can lead to non-negligible discrepancies.

M dwarfs show the highest rocky planet occurrence among all spectral types, in some instances within the Habitable Zone. Because some of them are very active stars, they are often subject to frequent and powerful flaring, which can be a double-edged sword in regard of exoplanet habitability. On one hand, the increased flux during flare events can trigger the chemical reactions that are necessary to build the basis of prebiotic chemistry. On the other hand, sufficiently strong flares may erode exoplanets' atmospheres and reduce their UV protection. Recent observations of flares have shown that the flaring flux can be x100 times stronger in UV than in the optical. UV is also preferable to constrain more accurately both the prebiotic abiogenesis and the atmospheric erosion. For these reasons, we are developing a CubeSat payload concept to complement current flare surveys operating in the optical. This CubeSat will observe a high number of flaring M dwarfs, following an all-sky scanning law coverage, both in the UV and the optical to better understand the different effective temperatures as wavelengths and flaring status go. This will complement the bright optical flares data acquired from the current ground-based, high-cadence, wide FoV surveys. Another scientific planned goal is to conduct few-minute after-the-flare follow-up optical ground-based time-resolved spectroscopy, that will be triggered by the detection of UV flares in space on board of the proposed CubeSat. Finally, the study of M dwarfs stellar activity in the UV band will provide useful data for larger forthcoming missions that will survey exoplanets, such as PLATO, ARIEL, HabEx and LUVOIR.

The recent extension of the Hubble diagram of Supernovae and quasars to redshifts much higher than 1 prompted a revived interest in non-parametric approaches to test cosmological models and to measure the expansion rate of the Universe. In particular, it is of great interest to infer model-independent constraints on the possible evolution of the dark energy component. Here we present a new method, based on a Neural Network Regression, to analyze the Hubble Diagram in a completely non-parametric, model-independent fashion. We first validate the method through simulated samples with the same redshift distribution as the real ones, and discuss the limitations related to the "inversion problem" for the distance-redshift relation. We then apply this new technique to the analysis of the Hubble diagram of Supernovae and quasars. We confirm that the data up to $z \sim 1-1.5$ are in agreement with a flat ${\Lambda}CDM$ model with ${\Omega}_M \sim 0.3$, while $\sim 5$-sigma deviations emerge at higher redshifts. A flat ${\Lambda}CDM$ model would still be compatible with the data with ${\Omega}_M > 0.4$. Allowing for a generic evolution of the dark energy component, we find solutions suggesting an increasing value of ${\Omega}_M$ with the redshift, as predicted by interacting dark sector models.

The Shower Database (SD) of the Meteor Data Center (MDC) has been operating for 15 years and is used by the entire community of meteor astronomers. It contains meteor showers categorised in individual lists on the basis of their status. Since the inception of the SD, no objective rules for moving showers between individual lists have been established. The content of the SD has not yet been checked for the correctness of the meteor data contained therein. Our aims are (1) to formulate criteria for nominating meteor showers for established status, (2) to improve the rules for the removal of showers, (3) to verify and enhance the content of the SD, and (4) to improve the user area of the MDC SD. The criteria for moving showers from the Working list to the Lists of established or removed Showers were generated using an empirical evaluation of their impact on the registered showers. The correctness of the parameters of each stream included in the SD was checked by comparing them with the values given in the source publications. We developed a set of criteria for nominating showers to be established. We objectified rules for the temporary and permanent removal of meteor showers from the Working list. Both of our proposed new procedures were approved by a vote of the commission F1 of the IAU in July 2022. We verified more than $1350$ data records of the MDC SD and introduced $\sim$1700 corrections. We included new parameters for shower characterisation. As a result of our verification procedure, 117 showers have been moved to the List of removed showers. As of October 2022, the SD contains 923 showers, 110 of which are in the List of established Showers and 813 are in the Working list. We also improved the user area of the SD and added a simple tool to allow a quick check of the similarity of a new shower to those in the database.

The maximum of a solar cycle contain two or more peaks, known as Gnevyshev peaks. Studies of this property of solar cycles may help for better understanding the solar dynamo mechanism. We analysed the 13-month smoothed monthly mean Version-2 international sunspot number (SN) during the period 1874-2017 and found that there exists a good correlation between the amplitude (value of the main and highest peak) and the value of the second maximum (value of the second highest peak) during the maximum of a solar cycle.Using this relationship and the earlier predicted value 86 (~92) of the amplitude of Solar Cycle 25, here we predict a value ~73 (~79) for the second maximum of Solar Cycle 25. The ratio of the predicted second maximum to the amplitude is found to be 0.85, almost the same as that of Solar Cycle 24.The cosine fits to the values of the peaks that occurred first and second during the maxima of Solar Cycles 12-24 suggests that in Solar Cycle 25 the second maximum would occur before the main maximum, the same as in Solar Cycle 24. However, these fits suggest ~106 and ~119 for the second maximum and the amplitude of Solar Cycle 25, respectively. We analysed the polar-fields data measured in Wilcox Observatory during Solar Cycles 20-24 and obtained a value ~125 for the amplitude of Solar Cycle 25. This is slightly larger-whereas the value ~86 (~92) predicted from the aforementioned relationship is slightly smaller-than the observed amplitude of Solar Cycle~24. This difference is discussed briefly.

Accretion disc coronae (ADC) sources are very high inclination neutron star or black hole binaries, where the outer accretion flow completely blocks a direct view of the central X-ray source. Instead, their observed weak X-rays are produced by scattering, line and recombination emission from highly ionised gas (the ADC) surrounding the source. The origin of this scattering material is still under debate. We use the ADC source 2S 0921-630 (V395 Car) to test whether it is consistent with a thermal-radiative wind produced by the central X-ray source illuminating the outer disc. This wind is clearly visible in blueshifted absorption lines in less highly inclined systems, where the source is seen directly through this material. Using the phenomenological photoionised plasma model, we first characterise the parameter that creates emission lines observed in 2S0921 in XMM-Newton and Chandra data. Then we run the Monte Carlos radiation transfer simulation to get scattered/reprocessed emissions in the wind, with density and velocity structure obtained from the previous work. Our model agrees with all the wind emission lines in the Chandra high and medium energy grating spectra for an intrinsic source luminosity of L > 0.2 LEdd. This result strongly favours thermal-radiative winds as the origin of the ADC. We also show how high-resolution spectra via microcalorimeters can provide a definitive test by detecting blueshifted absorption lines.

The local radio-loud AGN population is dominated by compact sources named FR0s. These sources show features, for example the host type, the mass of the supermassive black hole (SMBH), and the multi-band nuclear characteristics, that are similar to those of FRI radio galaxies. However, in the radio band, while FR0 and FRI share the same nuclear properties, the kiloparsec-scale diffuse component dominant in FRI is missing in FR0s. With this project we would like to study the parsec-scale structure in FR0s in comparison with that of FRI sources. To this end we observed 18 FR0 galaxies with the VLBA at 1.5 and 5 GHz and/or with the EVN at 1.7 GHz and produced detailed images at milliarcsec resolution of their nuclear emission to study the jet and core structure. All sources have been detected but one. Four sources are unresolved, even in these high-resolution images; jets have been detected in all other sources. We derived the distribution of the jet-to-counter-jet ratio of FR0s and found that it is significantly different from that of FRIs, suggesting different jet bulk speed velocities. Combining the present data with published data of FR0 with VLBI observations, we derive that the radio structure of FR0 galaxies shows strong evidence that parsec-scale jets in FR0 sources are mildly relativistic with a bulk velocity on the order of 0.5c or less. A jet structure with a thin inner relativistic spine surrounded by a low-velocity sheath could be in agreement with the SMBH and jet launch region properties.

Rotating radio transients (RRATs) are neutron stars that emit detectable radio bursts sporadically. They are outliers in the neutron star population, in many observable properties, but by their nature are practically difficult to study in depth. In this paper, we present the results from 1408 h of observations of RRAT candidates using the Irish Low Frequency Array (LOFAR) station at 150 MHz. As of October 2022, this census involved observing 113 sources, leading to 29 detections which were then followed up systematically. Single-pulse emission was detected from 25 sources, and periodic emission from 14 sources. 18 sources were found to have emission behaviour that is not discussed in prior works using LOFAR instruments. Four novel or modified source periods have been determined, ranging from 1.507-3.921 s, and 8 new or updated phase-coherent pulsar timing ephemerides have been produced using detected bursts. One unexpected single-pulse with a clearly-Galactic dispersion measure was detected as a part of this work, but has not been re-detected in follow-up observations. Observations are ongoing to expand the number of observed sources and further characterise and improve ephemerides for the detected sources. This census has demonstrated the capability for international LOFAR stations to detect, monitor and characterise a significant fraction of these unique sources.

We present analysis of 63 nearby ($<$ 44 Mpc) early-type galaxies hosting nuclear star clusters using the recently discovered parameter Central Intensity Ratio (CIR$_I$) determined from near-infra-red (3.6 $\mu$m) observations with the Infra-red-array-camera of \emph{Spitzer} space telescope. The CIR$_I$, when combined with filters involving age and $B-K$ colour of host galaxies, helps identify two distinct classes of galaxies hosting nuclear star clusters. This is independently verified using Gaussian Mixture Model. CIR shows a positive trend with faint, low mass, and blue galaxies in the sample, while the opposite is true for bright, high mass, and red galaxies, albeit with large scatter. The variation of CIR$_I$ with central velocity dispersion, absolute B band magnitude, dynamical mass, and stellar mass of host galaxies suggests that the mass of nuclear star clusters increases with that of host galaxies, for faint, low mass, young and blue galaxies in the sample. In bright, high-mass, old and red galaxies, on the other hand, the evolution of nuclear star clusters appears complex, with no apparent trends.The analysis also reveals that redder galaxies ($B-K > 3.76$) are more likely to be dominated by the central black-hole than the nuclear star clusters, while for bluer galaxies ($B-K < 3.76$) in the sample the situation is quite opposite.

When large overdensities gravitationally collapse in the early universe, they lead to primordial black holes (PBH). Depending on the exact model of inflation leading to necessary large perturbations at scales much smaller than scales probed at the Cosmic Microwave Background (CMB) surveys, PBHs of masses $\lesssim$$10^3 M_{\odot}$ can be formed sometime between the end of inflation and nucleosynthesis. However, the lack of a direct probe for the exact expansion history of the universe in this duration introduces uncertainties in the PBH formation process. The presence of alternate cosmological evolution for some duration after inflation affects the relation between (i) PBH mass and the scale of the collapsing overdensity; and (ii) PBH abundance and amplitude of the overdensities. In this review, the non-standard cosmological epochs relevant for a difference in PBH production are motivated and discussed. The importance of developing the framework of PBH formation in non-standard epochs is discussed from a phenomenological point of view, with particular emphasis on the advances in gravitational wave (GW) phenomenology, since abundant PBHs are always accompanied by large induced GWs. PBH formation in general non-standard epochs is also reviewed including the mathematical formalism. Specific examples, such as PBH formation in a kinetic energy dominated epoch and an early matter dominated epoch, are discussed with figures showing higher PBH abundances as compared to the production in standard radiation domination.

Sun-like stars shed angular momentum due to the presence of magnetised stellar winds. Magnetohydrodynamic models have been successful in exploring the dependence of this "wind-braking torque" on various stellar properties, however the influence of surface differential rotation is largely unexplored. As the wind-braking torque depends on the rotation rate of the escaping wind, the inclusion of differential rotation should effectively modulate the angular momentum-loss rate based on the latitudinal variation of wind source regions. In order to quantify the influence of surface differential rotation on the angular momentum-loss rate of the Sun, we exploit the dependence of the wind-braking torque on the effective rotation rate of the coronal magnetic field. This quantity is evaluated by tracing field lines through a Potential Field Source Surface (PFSS) model, driven by ADAPT-GONG magnetograms. The surface rotation rates of the open magnetic field lines are then used to construct an open-flux weighted rotation rate, from which the influence on the wind-braking torque can be estimated. During solar minima, the rotation rate of the corona decreases with respect to the typical solid-body rate (the Carrington rotation period is 25.4 days), as the sources of the solar wind shift towards the slowly-rotating poles. With increasing activity, more solar wind emerges from the Sun's active latitudes which enforces a Carrington-like rotation. The effect of differential rotation on the Sun's current wind-braking torque is found to be small. The wind-braking torque is ~10-15% lower during solar minimum, than assuming solid body rotation, and a few percent larger during solar maximum. For more rapidly-rotating Sun-like stars, differential rotation may play a more significant role, depending on the configuration of the large-scale magnetic field.

We here study the multi-band properties of a kiloparsec-size superbubble in the late-type spiral galaxy NGC628. The superbubble is the largest of many holes seen in the early release images using JWST/MIRI filters that trace the Polycyclic Aromatic Hydrocarbon (PAH) emissions. The superbubble is located in the interarm region ~3 kpc from the galactic center in the south-east direction. The shell surrounding the superbubble is detected in HI, CO, and Halpha with an expansion velocity of 12 km/s, and contains as much as 2x10^7 Msun of mass in gas that is mostly in molecular form. We find a clear excess of blue, bright stars inside the bubble as compared to the surrounding disk on the HST/ACS images. These excess blue, bright stars are part of a stellar population of 10^5 Msun mass that is formed over the last 50 Myr in different star formation episodes, as determined from an analysis of color-magnitude diagrams using a Bayesian technique. The mechanical power injected by the massive stars of these populations is sufficient to provide the energy necessary for the expansion of the shell gas. Slow and steady, rather than violent, injection of energy is probably the reason for the maintenance of the shell structure over the kiloparsec scale. The expanding shell is currently the site for triggered star formation as inferred from the JWST 21 micron (F2100W filter) and the Halpha images.

Galactic compact binaries with orbital periods shorter than a few hours emit detectable gravitational waves at low frequencies. Their gravitational wave signals can be detected with the future Laser Interferometer Space Antenna (LISA). Crucially, they may be useful in the early months of the mission operation in helping to validate LISA's performance in comparison to pre-launch expectations. We present an updated list of 48 candidate LISA binaries with measured properties, for which we derive distances based on Gaia Data release 3 astrometry. Based on the known properties from electromagnetic observations, we predict the LISA detectability after 1, 3, 6, and 48 months with state-of-the-art Bayesian analysis methods. We distinguish between verification and detectable binaries as being detectable after 3 and 48 months respectively. We find 16 verification binaries and 21 detectable sources, which triples the number of known LISA binaries over the last few years. These include detached double white dwarfs, AM CVn binaries, one ultracompact X-ray binary and two hot subdwarf binaries. We find that across this sample the gravitational wave amplitude is expected to be measured to $\approx10\%$ on average, while the inclination is expected to be determined with $\approx15^\circ$ precision. For detectable binaries these average errors increase to $\approx50\%$ and to $\approx40^\circ$ respectively.

We announce the discovery of activity, in the form of a distinct cometary tail, emerging from main-belt asteroid 2015 VA108. Activity was first identified by volunteers of the Citizen Science project Active Asteroids (a NASA Partner). We uncovered one additional image from the same observing run which also unambiguously shows 2015 VA108 with a tail oriented between the anti-solar and anti-motion vectors that are often correlated with activity orientation on sky. Both publicly available archival images were originally acquired UT 2015 October 11 with the Dark Energy Camera (DECam) on the Blanco 4 m telescope at the Cerro Tololo Inter-American Observatory (Chile) as part of the Dark Energy Camera Legacy Survey. Activity occurred near perihelion and, combined with its residence in the main asteroid belt, 2015 VA108 is a candidate main-belt comet, an active asteroid subset known for volatile sublimation.

The Y-dwarf WISE 1828+2650 is one of the coldest known Brown Dwarfs with an effective temperature of $\sim$300 K. Located at a distance of just 10 pc, previous model-based estimates suggest WISE1828+2650 has a mass of $\sim$5-10 Mj, making it a valuable laboratory for understanding the formation, evolution and physical characteristics of gas giant planets. However, previous photometry and spectroscopy have presented a puzzle with the near-impossibility of simultaneously fitting both the short (0.9-2.0 microns) and long wavelength (3-5 microns) data. A potential solution to this problem has been the suggestion that WISE 1828+2650 is a binary system whose composite spectrum might provide a better match to the data. Alternatively, new models being developed to fit JWST/NIRSpec and MIRI spectroscopy might provide new insights. This article describes JWST/NIRCam observations of WISE 1828+2650 in 6 filters to address the binarity question and to provide new photometry to be used in model fitting. We also report Adaptive Optics imaging with the Keck 10 m telescope. We find no evidence for multiplicity for a companion beyond 0.5 AU with either JWST or Keck. Companion articles will present low and high resolution spectra of WISE 1828+2650 obtained with both NIRSpec and MIRI.

The development and implementation of GEAR-RT, a radiative transfer solver using the M1 closure in the open source code SWIFT, is presented, and validated using standard tests for radiative transfer. GEAR-RT is modeled after RAMSES-RT (Rosdahl et al. 2013) with some key differences. Firstly, while RAMSES-RT uses Finite Volume methods and an Adaptive Mesh Refinement (AMR) strategy, GEAR-RT employs particles as discretization elements and solves the equations using a Finite Volume Particle Method (FVPM). Secondly, GEAR-RT makes use of the task-based parallelization strategy of SWIFT, which allows for optimized load balancing, increased cache efficiency, asynchronous communications, and a domain decomposition based on work rather than on data. GEAR-RT is able to perform sub-cycles of radiative transfer steps w.r.t. a single hydrodynamics step. Radiation requires much smaller time step sizes than hydrodynamics, and sub-cycling permits calculations which are not strictly necessary to be skipped. Indeed, in a test case with gravity, hydrodynamics, and radiative transfer, the sub-cycling is able to reduce the runtime of a simulation by over 90%. Allowing only a part of the involved physics to be sub-cycled is a contrived matter when task-based parallelism is involved, and is an entirely novel feature in SWIFT. Since GEAR-RT uses a FVPM, a detailed introduction into Finite Volume methods and Finite Volume Particle Methods is presented. In astrophysical literature, two FVPM methods are written about: Hopkins (2015) have implemented one in their GIZMO code, while the one mentioned in Ivanova et al. (2013) isn't used to date. In this work, I test an implementation of the Ivanova et al. (2013) version, and conclude that in its current form, it is not suitable for use with particles which are co-moving with the fluid, which in turn is an essential feature for cosmological simulations.

The past decade has seen an outstanding development of nonthermal particle acceleration in magnetic reconnection in magnetically-dominated systems, with clear signatures of power-law energy distributions as a common outcome of first-principles kinetic simulations. Here we propose a semi-analytical model for systematically investigating nonthermal particle acceleration in reconnection. We show particle energy distributions are well determined by particle injection, acceleration, and escape processes. Using a series of kinetic simulations, we accurately evaluate the energy- and time-dependent model coefficients. The resulting spectral characteristics, including the spectral index and lower and upper bounds of the power-law distribution, agree well with the simulation results. Finally, we apply the model to predict the power-law indices and break energies in astrophysical reconnection systems.

The history of accretion and differentiation processes in the planetesimals is provided by various groups of meteorites. Sampling different parent body layers, they reveal the circumstances of the metal-silicate segregation and the internal structures of the protoplanets. The ungrouped achondrite Erg Chech 002 (EC 002) added to the suite of samples from primitive igneous crusts. Here we present models that utilize thermo-chronological data for EC 002 and fit the accretion time and size of its parent body to these data. The U-corrected Pb-Pb, Al-Mg, and Ar-Ar ages used imply a best-fit planetesimal with a radius of 20-30 km that formed at 0.1 Ma after CAIs. Its interior melted early and differentiated by 0.5 Ma, allowing core and mantle formation with a transient lower mantle magma ocean, and a melt fraction of <25 % at the meteorite layering depth. EC 002 formed from this melt at a depth of 0.8 km in a partially differentiated region covered by an undifferentiated crust. By simulating collisions with impactors of different sizes and velocities, we analyzed the minimum ejection conditions of EC 002 from its original parent body and the surface composition of the impact site. The magma ocean region distinct from the layering depth of EC 002 implyes that it was not involved in the EC 002 genesis. Our models estimate closure temperatures for the Al-Mg ages as 1060 K to 1200 K. A fast parent body cooling attributes the late Ar-Ar age to a local reheating by another, late impact.

From the TNG50 cosmological simulation we build a sample of 191 well-resolved barred galaxies with a stellar mass $\logMstar > 10$ at $z=0$. We search for box/peanut bulges (BPs) in this sample, finding them in 55\% of cases. We compute $\fbp$, the BP probability for barred galaxies as a function of $\Mstar$, and find that this rises to a plateau, as found in observations of nearby galaxies. The transition mass where $\fbp$ reaches half the plateau value is $\logMstar = 10.14$, consistent with the observational value within measurement errors. We show that this transition in $\fbp$ can be attributed to the youth of the bars at low $\Mstar$, which is a consequence of downsizing of galaxies. Young bars, being generally shorter and weaker, have not yet had time to form BPs. At high mass, while we find a plateau, the value is at $\sim 60\%$ whereas observations saturate at $100\%$. We attribute this difference to excessive heating in TNG50, due to merger activity and to numerical resolution effects. BPs in TNG50 tend to occur in galaxies with more quiescent merger histories. As a result, the main driver of whether a bar hosts a BP in TNG50 is not the galaxy mass, but how long and strong the bar is. Separating the BP sample into those that have visibly buckled and those that have not, we find that fully half of BP galaxies show clear signs of buckling, despite the excessive heating and limited vertical resolution of TNG50.

Primordial non-Gaussianities (PNG) leave unique signatures in the bispectrum of the large-scale structure. With upcoming galaxy surveys set to improve PNG constraints by at least one order of magnitude, it is important to account for any potential contamination. In our work we show how to include wide-angle effects into the 3-dimensional observed galaxy bispectrum. We compute the leading wide-angle corrections to the monopole, finding that they could mimic local PNG with an amplitude of $f_{\rm NL} = \mathcal{O}\left( 0.1 \right)$. We also compute the dipole induced by wide-angle effects, whose amplitude is a few-percent of the flat-sky monopole. Our formalism can be readily adapted to realistic survey geometries and to include relativistic effects, which may become relevant at high redshifts.

We report the discovery of a remarkable Ly$\alpha$ emitting galaxy at $z=7.278$, JADES-GS+53.16746-27.7720 (shortened to JADES-GS-z7-LA), with EW$_0$(Ly$\alpha$) $\approx400 \pm 90$A and UV magnitude $-16.7$. The spectroscopic redshift is confirmed via rest-frame optical lines [O II], H$\beta$ and [O III] in its JWST/NIRSpec Micro-Shutter Assembly (MSA) spectrum. The Ly$\alpha$ line is detected in both lower resolution ($R\sim100$) PRISM as well as medium resolution ($R\sim1000$) G140M grating spectra. The Ly$\alpha$ FWHM in the grating is $\approx360$ km s$^{-1}$ and the line peaks within $120$ km s$^{-1}$ of the systemic redshift, indicative of very little neutral gas or dust within the galaxy. We estimate the Ly$\alpha$ escape fraction to be $\sim100\%$. JADES-GS-z7-LA has a [O III]/[O II] ratio (O32) of $8.8 \pm 1.1$ and ([O III]+[O II])/H$\beta$ ratio (R23) of $9.6\pm2.2$, consistent with low metallicity and high ionization parameters. Deep NIRCam imaging also revealed a close companion source (separated by $0.23''$), which exhibits similar photometry to that of JADES-GS-z7-LA, with a photometric excess in the F410M NIRCam image consistent with [O III]+H$\beta$ emission at the same redshift. The spectral energy distribution of JADES-GS-z7-LA indicates a `bursty' star-formation history, with a low stellar mass of $10^{7.15}$ M$_\odot$. The only explanation of the high EW Ly$\alpha$ emission seen in JADES-GS-z7-LA is if it resides in an ionized bubble with radius $\gtrsim 3$ pMpc. Owing to the faint nature of JADES-GS-z7-LA, we show that it is incapable of single-handedly ionizing a bubble large enough. Therefore, we suggest that JADES-GS-z7-LA (and possibly the companion source) may be a part of a larger overdensity, presenting direct evidence of overlapping ionized bubbles at $z>7$, allowing us to study the process of reionization across both small and large scales.

Purely collisionless Dark Matter Only (DMO) structure formation simulations predict that Dark Matter (DM) haloes are typically prolate in their centres and spheroidal towards their outskirts. The addition of gas cooling transforms the central DM shape to be rounder and more oblate. It is not clear, however, whether such shape transformations occur in `ultra-faint' dwarfs, which have extremely low baryon fractions. We present the first study of the shape and velocity anisotropy of ultra-faint dwarf galaxies that have gas mass fractions of $f_{\rm gas}(r<R_{\rm half}) < 0.06$. These dwarfs are drawn from the Engineering Dwarfs at Galaxy formation's Edge (EDGE) project, using high resolution simulations (spatial and mass resolution of 3 pc and $120$ M$_\odot$, respectively) that allow us to resolve DM halo shapes within the half light radius ($\sim 100$ pc). We show that gas-poor ultra-faints ($M_{\rm 200c} \leqslant 1.5\times10^9$ M$_\odot$; $f_{\rm gas} < 10^{-5}$) retain their pristine prolate DM halo shape even when gas, star formation and feedback are included. This could provide a new and robust test of DM models. By contrast, gas-rich ultra-faints ($M_{\rm 200c} > 3\times10^9$ M$_\odot$; $f_{\rm gas} > 10^{-4}$) become rounder and more oblate within $\sim 10$ half light radii. Finally, we find that most of our simulated dwarfs have significant radial velocity anisotropy that rises to $\tilde{\beta} > 0.5$ at $R \gtrsim 3 R_{\rm half}$. The one exception is a dwarf that forms a rotating gas/stellar disc because of a planar, major merger. Such strong anisotropy should be taken into account when building mass models of gas-poor ultra-faints.

Circumstellar disks do not evolve in isolation, as about half of solar-type stars were born in binary or multiple systems. Resolving disks in binary systems provides the opportunity to examine the influence of stellar companions on the outcomes of planet formation. We aim to investigate and compare disks in stellar multiple systems with near-infrared scattered-light imaging as part of the Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS) program. We used polarimetric differential imaging with SPHERE/IRDIS at the VLT to search for scattered light from the circumstellar disks in three multiple systems, CHX 22, S CrA, and HP Cha. We performed astrometric and orbit analyses for the stellar companions using archival HST, VLT/NACO, and SPHERE data. Combined with the age and orbital constraints, the observed disk structures provide insights into the evolutionary history and the impact of the stellar companions. The small grains in CHX 22 form a tail-like structure surrounding the close binary, which likely results from a close encounter and capture of a cloudlet. S CrA shows intricate structures (tentative ringed and spiral features) in the circumprimary disk as a possible consequence of perturbations by companions. The circumsecondary disk is truncated and connected to the primary disk via a streamer, suggesting tidal interactions. In HP Cha, the primary disk is less disturbed and features a tenuous streamer, through which the material flows towards the companions. The comparison of the three systems spans a wide range of binary separation (50 - 500 au) and illustrates the decreasing influence on disk structures with the distance of companions. This agrees with the statistical analysis of exoplanet population in binaries, that planet formation is likely obstructed around close binary systems, while it is not suppressed in wide binaries.

Asteroid restructuring uses robotics, self replication, and mechanical automatons to autonomously restructure an asteroid into a large rotating space station. The restructuring process makes structures from asteroid oxide materials; uses productive self-replication to make replicators, helpers, and products; and creates a multiple floor station to support a large population. In an example simulation, it takes 12 years to autonomously restructure a large asteroid into the space station. This is accomplished with a single rocket launch. The single payload contains a base station, 4 robots (spiders), and a modest set of supplies. Our simulation creates 3000 spiders and over 23,500 other pieces of equipment. Only the base station and spiders (replicators) have advanced microprocessors and algorithms. These represent 21st century technologies created and trans-ported from Earth. The equipment and tools are built using in-situ materials and represent 18th or 19th century technologies. The equipment and tools (helpers) have simple mechanical programs to perform repetitive tasks. The resulting example station would be a rotating framework almost 5 kilometers in diameter. Once completed, it could support a population of over 700,000 people. Many researchers identify the high launch costs, the harsh space environment, and the lack of gravity as the key obstacles hindering the development of space stations. The single probe addresses the high launch cost. The autonomous construction eliminates the harsh space environment for construction crews. The completed rotating station provides radiation protection and centripetal gravity for the first work crews and colonists.

The cosmic microwave background (CMB) is a significant source of knowledge about the origin and evolution of our universe. However, observations of the CMB are contaminated by foreground emissions, obscuring the CMB signal and reducing its efficacy in constraining cosmological parameters. We employ deep learning as a data-driven approach to CMB cleaning from multi-frequency full-sky maps. In particular, we develop a graph-based Bayesian convolutional neural network based on the U-Net architecture that predicts cleaned CMB with pixel-wise uncertainty estimates. We demonstrate the potential of this technique on realistic simulated data based on the Planck mission. We show that our model accurately recovers the cleaned CMB sky map and resulting angular power spectrum while identifying regions of uncertainty. Finally, we discuss the current challenges and the path forward for deploying our model for CMB recovery on real observations.

Unresolved sources of gravitational waves can create a stochastic gravitational wave background (SGWB) which may have intrinsic or extrinsic anisotropies. The angular power spectrum is a well-suited estimator for characterizing diffuse anisotropic distributions in the sky. Here we estimate the first model-independent all-sky all-frequency (ASAF) SGWB angular power spectra in the 20-1726 Hz frequency range from the third observing run (O3) of the Advanced LIGO and Advanced Virgo detectors. We develop a method to use the spectrum's signal-to-noise ratio (SNR) as the detection statistic and show that the distribution of the statistic obtained from the data agrees with the analytical model. Since we find the data to be consistent with noise, $95\%$ confidence Bayesian upper limits are set on the angular power spectra, ranging from $C_\ell^{1/2}\leq(3.1\times10^{-9}-0.76) \text{sr}^{-1}$. We also introduce a method to combine the narrowband angular power spectra to obtain estimators for broadband SGWB. These results can directly constrain theoretical models which predict the SGWB angular power spectra and for estimating or constraining the corresponding parameters. In addition, the results and the techniques introduced in this work can be useful for performing correlation-based searches, for instance, with electromagnetic observations.

We generalize the Rindler-Ishak (2007) result for the lensing deflection angle in an SdS spacetime, to the case of a general spherically symmetric fluid beyond the cosmological constant. We thus derive an analytic expression to first post-Newtonian order for the lensing deflection angle in a general asymptotically flat static spherically symmetric metric of the form $ ds^2 = f(r)dt^{2} -\frac{dr^{2}}{f(r)}-r^{2}(d\theta ^2 +\sin ^2 \theta d\phi ^2)$ with $f(r) = 1 - 2m \frac{r_0}{r}-\sum_{i} b_i\; r_0^{-q_i}\; \left( \frac{r_0}{r}\right)^{q_i}$ where $r_0$ is the lensing impact parameter, $b_i\ll r_0^{q_i}$, $m$ is the mass of the lens and $q_i$ are real arbitrary constants related to the properties of the fluid that surrounds the lens or to modified gravity. This is a generalization of the well known Kiselev black hole metric which involves a single fluid. The approximate analytic expression of the deflection angle is verified by an exact numerical derivation and in special cases it reduces to results of previous studies. The density and pressure of the spherically symmetric fluid that induces this metric is derived in terms of the constants $b_i$. The Kiselev case of a Schwarzschild metric perturbed by a general spherically symmetric dark fluid (eg vacuum energy) is studied in some detail and consistency with the special case of Rindler Ishak result is found for the case of a cosmological constant background. Observational data of the Einstein radii from distant clusters of galaxies lead to observational constraints on the constants $b_i$ and through them on the density and pressure of dark fluids, field theories or modified gravity theories that could induce this metric.

Detecting stochastic background radiation of cosmological origin is an exciting possibility for current and future gravitational-wave (GW) detectors. However, distinguishing it from other stochastic processes, such as instrumental noise and astrophysical backgrounds, is challenging. It is even more delicate for the space-based GW observatory LISA since it cannot correlate its observations with other detectors, unlike today's terrestrial network. Nonetheless, with multiple measurements across the constellation and high accuracy in the noise level, detection is still possible. Previous studies have always assumed that instrumental noise has a fixed and known spectral shape. In this study, we challenge this crucial assumption and assume that the single-link interferometric noises have an arbitrary and unknown spectrum. We investigate possible ways of separating instrumental and GW contributions by using realistic LISA data simulations with time-varying arms and second-generation time-delay interferometry. By fitting a generic spline model to the interferometer noise and a power-law template to the signal, we can detect GW stochastic backgrounds up to energy density levels comparable with fixed-shape models. We also demonstrate that we can probe a region of the GW background parameter space that today's detectors cannot access.