Papers for Wednesday, Jul 27 2022

The entropies of the known entities in the universe add to a total which is some twenty orders of magnitude below the holographic limit. Based on an assumption that the entropies should saturate the limit, we suggest that there exists dark matter, in the form of extremely massive primordial black holes, in addition to the dark matter known to exist inside galaxies and clusters of galaxies.

The two-body gravitational lens equation underlying planetary microlensing is usually transformed into a quintic polynomial that can only be solved numerically. Here, I present methods to acquire approximate analytic and exact semi-analytic solutions. First, I propose the pure-shear approximation, which allows one to acquire closed-form magnification solutions that are accurate apart from a small region near the primary star. While previous works on the perturbative picture suggest that the uniform-shear Chang-Refsdal lens only describes the vicinity of planetary caustics and breaks down in the resonant regime, the pure-shear lens formalism allows for all three caustic topologies. I show that the recently proposed offset degeneracy is a direct consequence of the pure-shear approximation. Second, the sole recognition that there always exists one image that is largely unaffected by the presence of the planet allows one to easily factor out the corresponding root from the quintic polynomial, reducing it to an analytically solvable quartic polynomial. This allows one to acquire semi-analytic solutions that are exact. The two analytic simplifications proposed here not only can allow for orders-of-magnitude faster forward models, but also enables the use of gradient-based inference algorithms that provide additional factors of acceleration for the analysis of observed events.

Low-mass compact stellar systems (CSSs; $M_{\star}$ $<$ 10$^{10}$ M$_{\odot}$) are thought to be a mixed bag of objects with various formation mechanisms. Previous surveys of CSSs were biased to relatively high-density environments and cannot provide a complete view of the environmental dependence of the formation of CSSs. We conduct the first-ever unbiased flux-limited census of nearby quiescent CSSs over a total sky area of $\sim$ 200 deg$^{2}$ observed by the GAMA spectroscopic survey. The complete sample includes 82 quiescent CSSs, of which 85\% fall within the stellar mass range of classical compact ellipticals (cEs).\ By quantifying the local environment with the normalized projected distance $D/R_{\rm vir}$ to the nearest luminous neighboring galaxy, we find that these CSSs have a bimodal $D/R_{\rm vir}$ distribution, with one group peaking near $\sim$ 0.1$\times$$R_{\rm vir}$ (satellite) and the other peaking near $\sim$ 10$\times$$R_{\rm vir}$ (field). In contrast to the CSSs, ordinary quiescent galaxies of similar masses have unimodal $D/R_{\rm vir}$ distribution.\ Satellite CSSs are older and more metal-rich than field CSSs on average. The bimodal $D/R_{\rm vir}$ distribution of quiescent CSSs reinforces the existence of two distinct formation channels (tidal stripping and born-to-be) for cEs, and may be understood in two mutually inclusive perspectives, i.e., substantial tidal stripping happens only when satellite galaxies travel sufficiently close to their massive hosts, and there exists an excess of high-density cE-bearing subhalos close to massive halos.

In the first paper of this series, we provide a method, Photometric objects Around Cosmic webs (PAC), that can make full use of the large-area photometric and cosmological spectroscopic surveys. With PAC we can measure the excess surface density $\bar{n}_2w_{\rm{p}}$ of the photometric objects with certain physical properties, such as stellar mass, around the spectroscopic objects. Combining with $w_{\rm{p}}$ measured from cosmological surveys, we can determine $\bar{n}_2$, the galaxy stellar mass function (GSMF), which is completely model independent. We measure the GSMFs in the redshift ranges of $z_s<0.2$, $0.2<z_s<0.4$ and $0.5<z_s<0.7$ down to the stellar mass $M_*=10^{8.2}$, $10^{10.6}$ and $10^{10.6}M_{\odot}$, using the data from the DESI Legacy Imaging Surveys and the spectroscopic samples of Slogan Digital Sky Survey (i.e. Main, LOWZ and CMASS samples). Our results show that there is no evolution of GSMF from $z_s=0.6$ to $z_s=0.1 $ for $M_*>10^{10.6} M_{\odot}$, and that there is a clear up-turn at $M_*\approx 10^{9.5} M_{\odot}$ towards smaller galaxies in the local GMSF at $z_s=0.1$. We provide an accurate double Schechter fit to the local GSMF for the entire range of $M_*$ and a table of our measurements at the three redshifts. Our method can achieve an accurate measurement of GSMF to the stellar mass limit where the spectroscopic sample is already highly incomplete (e.g. $\sim 10^{-3}$) for its target selection.

We present the BAT AGN Spectroscopic Survey (BASS) Near-infrared Data Release 2 (DR2), a study of 168 nearby ($\bar z$ = 0.04, $z$ < 0.6) active galactic nuclei (AGN) from the all-sky Swift Burst Array Telescope X-ray survey observed with Very Large Telescope (VLT)/X-shooter in the near-infrared (NIR; 0.8 - 2.4 $\mu$m). We find that 49/109 (45%) Seyfert 2 and 35/58 (60%) Seyfert 1 galaxies observed with VLT/X-shooter show at least one NIR high-ionization coronal line (CL, ionization potential $\chi$ > 100 eV). Comparing the emission of the [Si vi] $\lambda$1.9640 CL with the X-ray emission for the DR2 AGN, we find a significantly tighter correlation, with a lower scatter (0.37 dex) than for the optical [O iii] $\lambda$5007 line (0.71 dex). We do not find any correlation between CL emission and the X-ray photon index $\Gamma$. We find a clear trend of line blueshifts with increasing ionization potential in several CLs, such as [Si vi] $\lambda$1.9640, [Si x] $\lambda$1.4300, [S viii] $\lambda$0.9915, and [S ix] $\lambda$1.2520, indicating the radial structure of the CL region. Finally, we find a strong underestimation bias in black hole mass measurements of Sy 1.9 using broad H$\alpha$ due to the presence of significant dust obscuration. In contrast, the broad Pa$\alpha$ and Pa$\beta$ emission lines are in agreement with the $M$-$\sigma$ relation. Based on the combined DR1 and DR2 X-shooter sample, the NIR BASS sample now comprises 266 AGN with rest-frame NIR spectroscopic observations, the largest set assembled to date.

The BAT AGN Spectroscopic Survey (BASS) is designed to provide a highly complete census of the key physical parameters of supermassive black holes (SMBHs) that power local active galactic nuclei (AGN) (z<0.3), including their bolometric luminosity, black hole mass, accretion rates, and line-of-sight gas obscuration, and the distinctive properties of their host galaxies (e.g., star formation rates, masses, and gas fractions). We present an overview of the BASS data release 2 (DR2), an unprecedented spectroscopic survey in spectral range, resolution, and sensitivity, including 1449 optical (3200-10000 A) and 233 NIR (1-2.5 um) spectra for the brightest 858 ultra-hard X-ray (14-195 keV) selected AGN across the entire sky and essentially all levels of obscuration. This release provides a highly complete set of key measurements (emission line measurements and central velocity dispersions), with 99.9% measured redshifts and 98% black hole masses estimated (for unbeamed AGN outside the Galactic plane). The BASS DR2 AGN sample represents a unique census of nearby powerful AGN, spanning over 5 orders of magnitude in AGN bolometric luminosity, black hole mass, Eddington ratio, and obscuration. The public BASS DR2 sample and measurements can thus be used to answer fundamental questions about SMBH growth and its links to host galaxy evolution and feedback in the local universe, as well as open questions concerning SMBH physics. Here we provide a brief overview of the survey strategy, the key BASS DR2 measurements, data sets and catalogs, and scientific highlights from a series of DR2-based works.

We present the first cosmological study of a sample of $eROSITA$ clusters, which were identified in the $eROSITA$ Final Equatorial Depth Survey (eFEDS). In a joint selection on X-ray and optical observables, the sample contains $455$ clusters within a redshift range of $0.1<z<1.2$, of which $177$ systems are covered by the public data from the Hyper Suprime-Cam (HSC) survey that enables uniform weak-lensing cluster mass constraints. With minimal assumptions, at each cluster redshift $z$ we empirically model (1) the scaling relations between the cluster halo mass $M$ and the observables, which include the X-ray count rate $\eta$, the optical richness $\lambda$, and the weak-lensing mass $M_{\mathrm{WL}}$, and (2) the X-ray selection in terms of the completeness function $\mathcal{C}$. Using the richness distribution of the clusters, we directly measure the X-ray completeness and adopt those measurements as informative priors for the parameters of $\mathcal{C}$. We validate our cosmology analysis framework using mock data, and then in a blinded analysis of the real data we further validate the joint modeling of the cluster abundance and the weak-lensing mass calibration with a set of consistency tests. After unblinding, we obtain the cosmological constraints $\Omega_{\mathrm{m}} = 0.245^{+0.048}_{-0.058}$, $\sigma_{8} = 0.833^{+0.075}_{-0.063}$ and $S_{8} \equiv \sigma_{8}\left(\Omega_{\mathrm{m}}/0.3\right)^{0.3}= 0.791^{+0.028}_{-0.031}$ in a flat $\Lambda$CDM cosmology. Extending to a flat $w$CDM cosmology leads to the constraint on the equation of state parameter of the dark energy of $w = -1.25\pm 0.47$. The eFEDS constraints are in good agreement with the results from the $Planck$ mission, the galaxy-galaxy lensing and clustering analysis of the Dark Energy Survey, and the cluster abundance analysis of the SPT-SZ survey at a level of $\lesssim1\sigma$. (abridged)

We present ultra-deep Keck/MOSFIRE rest-optical spectra of two star-forming galaxies at z=2.18 in the COSMOS field with bright emission lines, representing more than 20~hours of total integration. The fidelity of these spectra enabled the detection of more than 20 unique emission lines for each galaxy, including the first detection of the auroral [O II]$\lambda\lambda$7322,7332 lines at high redshift. We use these measurements to calculate the electron temperature in the low-ionization O$^+$ zone of the ionized ISM and derive abundance ratios of O/H, N/H, and N/O using the direct method. The N/O and $\alpha$/Fe abundance patterns of these galaxies are consistent with rapid formation timescales and ongoing strong starbursts, in accord with their high specific star-formation rates. These results demonstrate the feasibility of using auroral [O II] measurements for accurate metallicity studies at high redshift in a higher metallicity regime previously unexplored with the direct method in distant galaxies. These results also highlight the difficulty in obtaining the measurements required for direct-method metallicities from the ground. We emphasize the advantages that the JWST/NIRSpec instrument will bring to high-redshift metallicity studies, where the combination of increased sensitivity and uninterrupted wavelength coverage will yield more than an order of magnitude increase in efficiency for multiplexed auroral-line surveys relative to current ground-based facilities. Consequently, the advent of JWST promises to be the beginning of a new era of precision chemical abundance studies of the early universe at a level of detail rivaling that of local galaxy studies.

Studying the physical and chemical properties of cold and dense molecular clouds is crucial for the understanding of how stars form. Under the typical conditions of infrared dark clouds, CO is removed from the gas phase and trapped on to the surface of dust grains by the so-called depletion process. This suggests that the CO depletion factor ($f_{\rm D}$) can be a useful chemical indicator for identifying cold and dense regions (i.e., prestellar cores). We have used the 1.3 mm continuum and C$^{18}$O(2-1) data observed at the resolution of $\sim$5000 au in the ALMA Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES) to construct averaged maps of $f_{\rm D}$ in twelve clumps to characterise the earliest stages of the high-mass star formation process. The average $f_{\rm D}$ determined for 277 of the 294 ASHES cores follows an unexpected increase from the prestellar to the protostellar stage. If we exclude the temperature effect due to the slight variations in the NH$_3$ kinetic temperature among different cores, we explain this result as a dependence primarily on the average gas density, which increases in cores where protostellar conditions prevail. This shows that $f_{\rm D}$ determined in high-mass star-forming regions at the core scale is insufficient to distinguish among prestellar and protostellar conditions for the individual cores, and should be complemented by information provided by additional tracers. However, we confirm that the clump-averaged $f_{\rm D}$ values correlates with the luminosity-to-mass ratio of each source, which is known to trace the evolution of the star formation process.

We present the AGN catalog and optical spectroscopy for the second data release of the Swift BAT AGN Spectroscopic Survey (BASS DR2). With this DR2 release we provide 1425 optical spectra, of which 1181 are released for the first time, for the 858 hard X-ray selected AGN in the Swift BAT 70-month sample. The majority of the spectra (813/1425, 57%) are newly obtained from VLT/Xshooter or Palomar/Doublespec. Many of the spectra have both higher resolution (R>2500, N~450) and/or very wide wavelength coverage (3200-10000 A, N~600) that are important for a variety of AGN and host galaxy studies. We include newly revised AGN counterparts for the full sample and review important issues for population studies, with 44 AGN redshifts determined for the first time and 780 black hole mass and accretion rate estimates. This release is spectroscopically complete for all AGN (100%, 858/858) with 99.8% having redshift measurements (857/858) and 96% completion in black hole mass estimates of unbeamed AGN (outside the Galactic plane). This AGN sample represents a unique census of the brightest hard X-ray selected AGN in the sky, spanning many orders of magnitude in Eddington ratio (Ledd=10^-5-100), black hole mass (MBH=10^5-10^10 Msun), and AGN bolometric luminosity (Lbol=10^40-10^47 ergs/s).

We present new central stellar velocity dispersions for 484 Sy 1.9 and Sy 2 from the second data release of the Swift/BAT AGN Spectroscopic Survey (BASS DR2). This constitutes the largest study of velocity dispersion measurements in X-ray selected, obscured AGN with 956 independent measurements of the Ca H+K and Mg b region (3880-5550A) and the Ca triplet region (8350-8730A) from 642 spectra mainly from VLT/Xshooter or Palomar/DoubleSpec. Our sample spans velocity dispersions of 40-360 km/s, corresponding to 4-5 orders of magnitude in black holes mass (MBH=10^5.5-9.6 Msun), bolometric luminosity (LBol~10^{42-46 ergs/s), and Eddington ratio (L/Ledd~10^{-5}-2). For 281 AGN, our data provide the first published central velocity dispersions, including 6 AGN with low mass black holes (MBH=10^5.5-6.5 Msun), discovered thanks to our high spectral resolution observations (sigma~25 km/s). The survey represents a significant advance with a nearly complete census of hard-X-ray selected obscured AGN with measurements for 99% of nearby AGN (z<0.1) outside the Galactic plane. The BASS AGN have higher velocity dispersions than the more numerous optically selected narrow line AGN (i.e., ~150 vs. ~100 km/s), but are not biased towards the highest velocity dispersions of massive ellipticals (i.e., >250 km/s). Despite sufficient spectral resolution to resolve the velocity dispersions associated with the bulges of small black holes (~10^4-5 Msun), we do not find a significant population of super-Eddington AGN. Using estimates of the black hole sphere of influence, direct stellar and gas black hole mass measurements could be obtained with existing facilities for more than ~100 BASS AGN.

Direct observational measurements of the magnetic field strength in prestellar cores typically find supercritical mass-to-flux ratios, suggesting that the magnetic field is insufficient to prevent gravitational collapse. These measurements suffer from significant uncertainties; an alternative approach is to utilise the sensitivity of prestellar chemistry to the evolutionary history, and indirectly constrain the degree of magnetic support. We combine non-ideal magnetohydrodynamic simulations of prestellar cores with time-dependent chemistry and radiative transfer modelling, producing synthetic observations of the model cores in several commonly-observed molecular lines. We find that molecules strongly affected by freeze-out, such as CS and HCN, typically have much lower line intensities in magnetically subcritical models compared to supercritical ones, due to the longer collapse timescales. Subcritical models also produce much narrower lines for all species investigated. Accounting for a range of core properties, ages and viewing angles, we find that supercritical models are unable to reproduce the distribution of CS and N$_2$H$^+$ line strengths and widths seen in an observational sample, whereas subcritical models are in good agreement with the available data. This suggests that despite presently having supercritical mass-to-flux ratios, prestellar cores form as magnetically subcritical objects.

Galaxies with stellar masses as high as $M_*\sim 10^{11}M_\odot$ have been identified out to redshifts $z\sim 6$, approximately one billion years after the Big Bang. It has been difficult to find massive galaxies at even earlier times, as the Balmer break region, which is needed for accurate mass estimates, is redshifted to wavelengths $>2.5\,\mu$m. Here we make use of the excellent long-wavelength coverage of the JWST early release observations to search for massive galaxies in the first $\approx 750$ million years of cosmic history. We find seven galaxies with $M_*>10^{10}M_\odot$ and $7<z<11$ in the survey area, including two galaxies with $M_*\sim 10^{11}M_\odot$. The stellar mass density in massive galaxies is much higher than anticipated from previous studies based on rest-frame UV-selected samples: a factor of 10-30 at $z\sim 8$ and more than three orders of magnitude at $z\sim 10$. From these first JWST images we infer that the central regions of at least some massive galaxies were already largely in place 500 Myr after the Big Bang, and that massive galaxy formation began extremely early in the history of the Universe. The presence of these galaxies at $z\sim 10$ suggests that galaxies with masses $M_* \sim 5\times 10^9 M_{\odot}$ may be found out to redshifts as high as $z\sim 18$.

Our basic theoretical understanding of the sources of cosmic rays and their propagation through the interstellar medium is hindered by the Sun, that through the solar wind affects the observed cosmic-ray spectra. This effect is known as solar modulation. Recently released cosmic-ray observations from the Alpha Magnetic Spectrometer (AMS-02) and publicly available measurements of the solar wind properties from the Advanced Composition Explorer and the Wilcox observatory allow us to test the analytical modeling of the time-, charge- and rigidity-dependence of solar modulation. We rely on associating measurements on the local heliospheric magnetic field and the heliospheric current sheet's tilt angle, to model the time-dependence and amplitude of cosmic-ray solar modulation. We find evidence for the solar modulation's charge- and rigidity-dependence during the era of solar cycle 24. Our analytic prescription to model solar modulation can explain well the large-scale time-evolution of positively charged cosmic-ray fluxes in the range of rigidities from 1 to 10 GV. We also find that cosmic-ray electron fluxes measured during the first years of cycle 24 are less trivial to explain, due to the complex and rapidly evolving structure of the Heliosphere's magnetic field that they experienced as they propagated inwards.

The SPectrogram Analysis and Cataloguing Environment (SPACE) tool is an interactive python tool designed to label radio emission features of interest in a time-frequency map (called 'dynamic spectrum'). The program uses Matplotlib's Polygon Selector widget to allow a user to select and edit an undefined number of vertices on top of the dynamic spectrum before closing the shape (polygon). Multiple polygons may be drawn on any spectrum, and the feature name along with the coordinates for each polygon vertex are saved into a '.json' file as per the 'Time-Frequency Catalogue' (TFCat) format along with other data such as the feature id, observer name, and data units. This paper describes the first official stable release (version 2.0) of the tool.

Cosmic rays escaping the Milky-Way disk interact with circumgalactic gas which fills the virial volume of our Galaxy. These interactions should produce guaranteed fluxes of energetic diffuse neutrinos and photons observable at the Earth. This neutrino flux would be a plausible contribution to the spectrum measured by the IceCube neutrino observatory: the energy emitted in this way is weakly constrained from cascade gamma rays, since the cascades have no time to develop, but the arrival directions of the neutrinos do not point to the Galactic disk, in agreement with observations. However, previous studies reported very different estimates of the corresponding neutrino flux, so it was unclear if this contribution to the observed spectrum is essential. Here we readdress the calculation of this diffuse neutrino flux component under various assumptions about the cosmic-ray spectrum and propagation in the circumgalactic medium. We find that this contribution to the observed neutrino flux is always subleading.

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) aims to improve constraints on the dark energy equation of state through measurements of large-scale structure at high redshift ($0.8<z<2.5$), while serving as a state-of-the-art fast radio burst detector. Bright galactic foregrounds contaminate the 400--800~MHz HIRAX frequency band, so meeting the science goals will require precise instrument characterization. In this paper we describe characterization of the HIRAX antenna, focusing on measurements of the antenna beam and antenna noise temperature. Beam measurements of the current HIRAX antenna design were performed in an anechoic chamber and compared to simulations. We report measurement techniques and results, which find a broad and symmetric antenna beam for $\nu <$650MHz, and elevated cross-polarization levels and beam asymmetries for $\nu >$700MHz. Noise temperature measurements of the HIRAX feeds were performed in a custom apparatus built at Yale. In this system, identical loads, one cryogenic and the other at room temperature, are used to take a differential (Y-factor) measurement from which the noise of the system is inferred. Several measurement sets have been conducted using the system, involving CHIME feeds as well as four of the HIRAX active feeds. These measurements give the first noise temperature measurements of the HIRAX feed, revealing a $\sim$60K noise temperature (relative to 30K target) with 40K peak- to-peak frequency-dependent features, and provide the first demonstration of feed repeatability. Both findings inform current and future feed designs.

We study the half-light radius versus black hole mass as well as the luminosity versus black hole mass relations in active galactic nuclei (AGN) when the disc is illuminated by the X-ray corona. We used KYNSED, a recently developed spectral model for studying broadband spectral energy distribution in AGN. We considered non-illuminated Novikov-Thorne discs and X-ray illuminated discs based on a Novikov-Thorne temperature radial profile. We also considered the case where the temperature profile is modified by a colour-correction factor. We assumed that the X-ray luminosity is equal to the accretion power that is dissipated to the disc below a transition radius. The half-light radius of X-ray illuminated radii can be up to some three to four times greater than the radius of a standard disc, even for a non-spinning black hole, due to the fact that the absorbed X-rays act as a secondary source of energy, increasing the disc temperature. Non-illuminated discs are consistent with observations, but only at the 2.5sigma level. On the other hand, X-ray illuminated discs can explain both the half-light radius-black hole mass as well as the luminosity-black hole mass relation in AGN, for a wide range of physical parameters. In addition, we show that the observed X-ray luminosity of the gravitationally lensed quasars is fully consistent with the X-ray luminosity that is necessary for heating the disc. X-ray disc illumination was proposed many years ago to explain various features that are commonly observed in the X-ray spectra of AGN. Recently, we showed that X-ray illumination of accretion disc can also explain the observed UV/optical time-lags in AGN, while in this work, we show that the same model can also account for the quasar micro-lensing disc size problem. These results support the hypothesis of the disc X-ray illumination in AGN.

We report the discovery of a candidate galaxy with a photo-z of z~14 in the first epoch of JWST NIRCam imaging from the Cosmic Evolution Early Release Science (CEERS) Survey. We searched for z>12 sources using photo-z's based primarily on the expected Lyman-alpha breaks. Following conservative vetting criteria, we identify a robust source at z_phot=14.3 (+0.4, -1.1, 1-sig uncertainty) with F277W=27.8, and detections in five (two) filters at >5sig (>10sig) significance. This object (dubbed Maisie's Galaxy) exhibits a strong F150W-F200W Lyman-alpha break color of >2.5 mag (1sig), resulting in 99.99% (87%) of the photo-z PDF favoring z>10 (13). The source may be marginally detected by HST in F160W, which would widen the lower-redshift bound to z~12.5. This slightly lower-redshift interpretation would require very strong Ly-alpha emission (>300 A rest equivalent width), requiring an extremely early ionized bubble. The colors of this object are inconsistent with Galactic stars, and it is resolved (r_h=0.1+/-0.01"; 330 pc). All data quality images show no artifacts at the candidate's position. Independent methods are in robust agreement, and we consistently find a strong preference for z>13 with our fiducial photometry. Maisie's Galaxy appears to be fairly high mass (log M*/Msol~8.5) and highly star-forming (log sSFR ~ -7.9/yr) for this early epoch, with a blue rest-UV color (beta~-2.3) indicating little dust, though not extremely low metallicities. While the presence of this source is in tension with most cosmological simulation predictions and may seriously challenge dark matter models with suppressed small-scale power, it is in agreement with empirical extrapolations from lower redshift assuming a smoothly declining cosmic SFR density. Should followup spectroscopy validate this redshift, our Universe was already aglow with fairly massive galaxies less than 300 Myr after the Big Bang.

The environments where galaxies reside play a key role in shaping their star formation histories over cosmic time, yet such environmental effects remain elusive at high redshifts. We update this environmental process by adopting an advanced measure for cluster in-fall time using kinematics and reveal that galaxies experience a gradual decline of star formation after they fall into cluster environments up to $z \sim 1$. This conclusion is drawn from a uniform analysis of a remarkably large sample of 105 clusters and 1626 spectroscopically-confirmed member galaxies from the SPT and ACT Sunyaev-Zel'dovich surveys at 0.26 $< z <$ 1.13. Intriguingly, we find clear evidence for a gradual increase in the mean age ($\sim$ 0.71 $\pm$ 0.4 Gyr based on a 4000 Angstrom break, $\rm D_{\rm n}4000$) of the galaxy's stellar populations with the time spent in the cluster environment. This environmental quenching effect is found regardless of galaxy luminosity (faint or bright) and redshift (low-$z$ or high-$z$), although the exact stellar age of galaxies depends on both parameters at fixed environmental effect. Such a systematic increase of $\rm D_{\rm n}4000$ with in-fall proxy would suggest that galaxies that were accreted into hosts earlier were quenched earlier, due to longer exposure to environmental effects such as ram pressure stripping and strangulation. Thus, our results provide new insights into environmental quenching effects spanning a large range in cosmic time ($\sim 5.2$ Gyr, $z=0.26-1.13$) and demonstrate the power of using a kinematically-derived in-fall time proxy.

In recent years the number of exoplanets has grown considerably. The most successful techniques in these detections are the radial velocity (RV) and planetary transits techniques, the latter significantly advanced by the Kepler, K2 and, more recently, the TESS missions. The detection of exoplanets both by means of transit and by RVs is of importance, because this would allows characterizing their bulk densities, and internal compositions. The Transiting Exoplanet Survey Satellite (TESS) survey offers a unique possibility to search for transits of extrasolar planets detected by RV. In this work, we present the results of the search for transits of planets detected with the radial velocity technique, using the photometry of the TESS space mission. We focus on systems with super-Earths and Neptunes planets on orbits with periods shorter than 30 days. This cut is intended to keep objects with a relatively high transit probability, and is also consistent with duration of TESS observations on a single sector. Given the summed geometric transit probabilities, the expected number of transiting planets is $3.4 \pm 1.8$. The sample contains two known transiting planets. We report null results for the remaining 66 out of 68 planets studied, and we exclude in all cases planets larger than 2.4 R$_{\oplus}$, under the assumption of central transits. The remaining two planets orbit HD~136352 and have been recently been announced.

The increasing sensitivity and resolution of ground-based telescopes have enabled the detection of gas-phase complex organic molecules (COMs) across a variety of environments. Many of the detected species are expected to form on the icy surface of interstellar grains and transfer later into the gas phase. Therefore, icy material is regarded as a primordial source of COMs in the ISM. Upcoming JWST observations of interstellar ices in star-forming regions will reveal IR features of frozen molecules with unprecedented resolution and sensitivity. To identify COM features in the JWST data, lab IR spectra of ices for conditions that simulate interstellar environments are needed. This work presents FTIR spectra (500-4000 cm$^{-1}$/20-2.5$\mu$m, with a resolution of 1 cm$^{-1}$) of methyl cyanide (CH$_3$CN, aka acetonitrile) mixed with H$_2$O, CO, CO$_2$, CH$_4$, and NH$_3$, at temperatures from 15-150 K. The refractive index of pure amorphous CH$_3$CN ice at 15K and the band strength, peak position, and FWHM of selected IR bands are also measured. These bands are: the CH$_3$ sym stretching at 2940.9 cm$^{-1}$, the CN stretching at 2252.2 cm$^{-1}$, a combination of modes at 1448.3 cm$^{-1}$ , the CH$_3$ antisym def. at 1410 cm$^{-1}$ , the CH$_3$ sym def. at 1374.5 cm$^{-1}$, and the CH$_3$ rock at 1041.6 cm$^{-1}$. The lab spectra of CH$_3$CN are compared to observations of ices toward W33A and three low-mass YSOs. Since an unambiguous identification of CH$_3$CN is not possible, upper limits for the CH$_3$CN column density are determined as $\leq 2.4\times 10^{17}$ molecules cm$^{-2}$ for W33A and $5.2 \times 10^{16}$, $1.9\times 10^{17}$, and $3.8\times 10^{16}$ molecules cm$^{-2}$ for EC92, IRAS 03235, and L1455 IRS3, respectively. W.r.t solid H$_2$O, these values correspond to relative abundances of 1.9, 3.1, 1.3, and 4.1\%, for W33A, EC92, IRAS 03235, and L1455 IRS3, respectively.

We present a method to reconstruct the initial linear-regime matter density field from the late-time non-linearly evolved density field in which we channel the output of standard first-order reconstruction to a convolutional neural network (CNN). Our method shows dramatic improvement over the reconstruction of either component alone. We show why CNNs are not well-suited for reconstructing the initial density directly from the late-time density: CNNs are local models, but the relationship between initial and late-time density is not local. Our method leverages standard reconstruction as a preprocessing step, which inverts bulk gravitational flows sourced over very large scales, transforming the residual reconstruction problem from long-range to local and making it ideally suited for a CNN. We develop additional techniques to account for redshift distortions, which warp the density fields measured by galaxy surveys. Our method improves the range of scales of high-fidelity reconstruction by a factor of 2 in wavenumber above standard reconstruction, corresponding to a factor of 8 increase in the number of well-reconstructed modes. In addition, our method almost completely eliminates the anisotropy caused by redshift distortions. As galaxy surveys continue to map the Universe in increasingly greater detail, our results demonstrate the opportunity offered by CNNs to untangle the non-linear clustering at intermediate scales more accurately than ever before.

Convective overshoot mixing is a critical ingredient of stellar structure models, but is treated in most cases by ad hoc extensions of the mixing-length theory for convection. Advanced theories which are both more physical and numerically treatable are needed. Convective flows in stellar interiors are highly turbulent. This poses a number of numerical challenges for the modelling of convection in stellar interiors. We include an effective turbulence model into a 1D stellar evolution code in order to treat non-local effects within the same theory. We use a turbulent convection model which relies on the solution of second order moment equations. We implement this into a state of the art 1D stellar evolution code. To overcome a deficit in the original form of the model, we take the dissipation due to buoyancy waves in the overshooting zone into account. We compute stellar models of intermediate mass main-sequence stars between 1.5 and 8 $M_\odot$. Overshoot mixing from the convective core and modified temperature gradients within and above it emerge naturally as a solution of the turbulent convection model equations. For a given set of model parameters the overshooting extent determined from the turbulent convection model is comparable to other overshooting descriptions, the free parameters of which had been adjusted to match observations. The relative size of the mixed cores decreases with decreasing stellar mass without additional adjustments. We find that the dissipation by buoyancy waves constitutes a necessary and relevant extension of the turbulent convection model in use.

We present estimates for the occurrence rates of hot Jupiters around dwarf stars based on data from the Transiting Exoplanet Survey Satellite (TESS) Prime Mission. We take 97 hot Jupiters orbiting 198,721 AFG dwarf stars (ranging in mass from $0.8M_{\odot}$ to $2.3M_{\odot}$) from an independent search for hot Jupiters using TESS Prime Mission data. We estimate our planet sample's false positive rates as $14\pm7\%$ for A stars, $16\pm6\%$ for F stars, and $0\%$ for G stars. We find hot Jupiter occurrence rates of $0.29 \pm 0.05\%$ for A stars, $0.36 \pm 0.06\%$ for F stars and $0.55 \pm 0.14\%$ for G stars, with a weighted average across AFG stars of $0.33\pm0.04\%$. Our results show a correlation between higher hot Jupiter abundance and lower stellar mass, and are in good agreement with occurrence rates found by Kepler. After correcting for the presence of binaries in the TESS stellar sample, we estimate a single-star hot Jupiter occurrence rate of $0.98\pm0.36\%$ for G stars. This is in agreement with results from radial velocity (RV) surveys, indicating that stellar multiplicity correction is able to resolve the discrepancy between hot Jupiter occurrence rates based on transits and RVs.

Scattering transforms have been successfully used to describe dust polarization for flat-sky images. This paper expands this framework to noisy observations on the sphere with the aim of obtaining denoised Stokes Q and U all-sky maps at 353GHz, as well as a non-Gaussian model of dust polarization, from the Planck data. To achieve this goal, we extend the computation of scattering coefficients to the Healpix pixelation and introduce cross-statistics that allow us to make use of half-mission maps as well as the correlation between dust temperature and polarization. Introducing a general framework, we develop an algorithm that uses the scattering statistics to separate dust polarization from data noise. The separation is validated on mock data, before being applied to the SRoll2 Planck maps at N_side = 256. The validation shows that the statistics of the dust emission, including its non-Gaussian properties, are recovered until l < 700, where, at high Galactic latitudes, the dust power is smaller than that of the dust by two orders of magnitude. On scales where the dust power is lower than one tenth of that of the noise, structures in the output maps have comparable statistics but are not spatially coincident with those of the input maps. Our results on Planck data are significant milestones opening new perspectives for statistical studies of dust polarization and for the simulation of Galactic polarized foregrounds. The Planck denoised maps will be made available together with results from our validation on mock data, which may be used to quantify uncertainties.

High-precision astrometry well beyond the capacities of Gaia will provide a unique way to achieve astrophysical breakthroughs, in particular on the nature of dark matter, and a complete survey of nearby habitable exoplanets. In this contribution, we present the scientific cases that require a flexibly-pointing instrument capable of high astrometric accuracy and we review the best mission profiles that can achieve such observations with the current space technology as well as within the boundary conditions defined by space agencies. We also describe the way the differential astrometric measurement is made using reference stars within the field. We show that the ultimate accuracy can be met without drastic constrains on the telescope stability.

We present a new optical transmission spectrum of the hot Jupiter HAT-P-1b based on two transits observed with the Double Spectrograph (DBSP) on the Palomar 200-inch (P200) telescope. The DBSP transmission spectrum, covering a wavelength range from 325 to 1001 nm, is consistent with that observed with the Hubble Space Telescope (HST), but the former has a finer spectral resolution. The DBSP spectrum alone reveals the presence of a pressure broadened line wing for Na, the line cores for both Na and K, and tentative evidence for H$_2$O. We obtain consistent results from the spectral retrieval analyses performed on the DBSP-only dataset and the DBSP, HST, and Spitzer combined dataset. Our retrievals suggest a mostly clear atmosphere for HAT-P-1b, with a cloud coverage of $22^{+5}_{-3}$% that is dominated by enhanced haze. We derive subsolar abundances for Na, K, and C, and subsolar-to-solar for O. Future observations with James Webb Space Telescope and ground-based high-resolution spectrographs should be able to not only confirm the presence of these species but also stringently constrain the formation and migration pathways for HAT-P-1b.

Recent interferometers (e.g. ALMA and NOEMA) allow us to obtain the detailed brightness distribution of the astronomical sources in 3 dimension (R.A., Dec., frequency). However, the interpixel correlation of the noise due to the limited uv coverage makes it difficult to evaluate the statistical uncertainty of the measured quantities and the statistical significance of the obtained results. The noise correlation properties are characterized by the noise autocorrelation function (ACF). We will present the method for (1) estimating the statistical uncertainty due to the correlated noise in the spatially integrated flux and spectra directly from the noise ACF and (2) simulating the correlated noise to perform a Monte Carlo simulation in image analyses. Our method has potential applications to a range of astronomical images of not only interferometers but also single dish mapping observation and interpolated and resampled optical images.

Measurement of magnetic fields in dense molecular clouds is essential for understanding the fragmentation process prior to star formation. Radio interferometric observations of CCS 22.3 GHz emission, from the starless core TMC-1C, have been carried out with the Karl G. Jansky Very Large Array to search for Zeeman splitting of the line in order to constrain the magnetic field strength. Toward a region offset from the dust peak, we report a detection of the Zeeman splitting of the CCS 2_1 - 1_0 transition, with an inferred magnetic field of ~2 mG. If we interpret the dust peak to be the core of TMC-1C, and the region where we have made a detection of the magnetic field to be the envelope, then our observed value for the magnetic field is consistent with a subcritical mass-to-flux ratio envelope around a core with supercritical mass-to-flux ratio. The ambipolar diffusion timescale for the formation of the core is consistent with the relevant timescale based on chemical modeling of the TMC-1C core. This work demonstrates the potential of deep CCS observation to carry out future measurements of magnetic field strengths in dense molecular clouds and, in turn, understand the role of the magnetic field in star formation.

We present observations of recurrent active region coronal jets and derive their thermal and non-thermal properties, by studying the physical properties of the plasma simultaneously at the base footpoint, and along the outflow of jets. The sample of analyzed solar jets were observed by SDO-AIA in Extreme Ultraviolet and by RHESSI in the X-Ray domain. The main thermal plasma physical parameters: temperature, density, energy flux contributions, etc. are calculated using multiple inversion techniques to obtain the differential emission measure from extreme-ultraviolet filtergrams. The underlying models are assessed, and their limitations and applicability are scrutinized. Complementarily, we perform source reconstruction and spectral analysis of higher energy X-Ray observations to further assess the thermal structure and identify non-thermal plasma emission properties. We discuss a peculiar penumbral magnetic reconnection site, which we previously identified as a ``Coronal Geyser''. Evidence supporting cool and hot thermal emission, and non-thermal emission, is presented for a subset of geyser jets. These active region jets are found to be energetically stronger than their polar counterparts, but we find their potential influence on heliospheric energetics and dynamics to be limited. We scrutinize whether the geyser does fit the non-thermal erupting microflare picture, finding that our observations at peak flaring times can only be explained by a combination of thermal and non-thermal emission models. This analysis of geysers provides new information and observational constraints applicable to theoretical modeling of solar jets.

Incremental particle growth in turbulent protoplanetary nebulae is limited by a combination of barriers that can slow or stall growth. Moreover, particles that grow massive enough to decouple from the gas are subject to inward radial drift which could lead to the depletion of most disk solids before planetesimals can form. Compact particle growth is probably not realistic. Rather, it is more likely that grains grow as fractal aggregates which may overcome this so-called radial drift barrier because they remain more coupled to the gas than compact particles of equal mass. We model fractal aggregate growth and compaction in a viscously evolving solar-like nebula for a range of turbulent intensities $\alpha_{\rm{t}} = 10^{-5}-10^{-2}$. We do find that radial drift is less influential for porous aggregates over much of their growth phase; however, outside the water snowline fractal aggregates can grow to much larger masses with larger Stokes numbers more quickly than compact particles, leading to rapid inward radial drift. As a result, disk solids outside the snowline out to $\sim 10-20$ AU are depleted earlier than in compact growth models, but outside $\sim 20$ AU material is retained much longer because aggregate Stokes numbers there remain lower initially. Nevertheless, we conclude even fractal models will lose most disk solids without the intervention of some leap-frog planetesimal forming mechanism such as the Streaming Instability (SI), though conditions for the SI are generally never satisfied, except for a brief period %for a brief stage around $\sim 0.2$ Myr at the snowline for $\alpha_{\rm{t}}=10^{-5}$.

Formation of the first planetesimals remains an unsolved problem. Growth by sticking must initiate the process, but multiple studies have revealed a series of barriers that can slow or stall growth, most of them due to nebula turbulence. In a companion paper, we study the influence of these barriers on models of fractal aggregate and solid, compact particle growth in a viscously evolving solar-like nebula for a range of turbulent intensities $\alpha_{\rm{t}} = 10^{-5}-10^{-2}$. Here, we examine how disk composition in these same models changes with time. We find that advection and diffusion of small grains and vapor, and radial inward drift for larger compact particles and fractal aggregates, naturally lead to diverse outcomes for planetesimal composition. Larger particles can undergo substantial inward radial migration due to gas drag before being collisionally fragmented or partially evaporating at various temperatures. This leads to enhancement of the associated volatile in both vapor inside, and solids outside, their respective evaporation fronts, or ``snowlines''. In cases of lower $\alpha_{\rm{t}}$, we see narrow belts of volatile or supervolatile material develop in the outer nebula, which could be connected to the bands of ``pebbles" seen by ALMA. Volatile bands, which migrate inwards as the disk cools, can persist over long timescales as their gas phase continues to advect or diffuse outward across its evaporation front. These belts could be sites where supervolatile-rich planetesimals form, such as the rare CO-rich and water-poor comets; giant planets formed just outside the H$_2$O snowline may be enhanced in water.

The onset of the {\it JWST}-era provides a much-improved opportunity to characterize the resolved structure of early star forming systems. Previous observations of $z\gtrsim 6$ galaxies with {\it Spitzer} revealed the presence of old stars and luminous HII regions (via [OIII]+H$\beta$ emission), but the poor resolution stunted our ability to map the location of these components with respect to the star forming regions identified in the rest-UV. As such, it has long been unclear whether the old stars are situated in a separate nuclear component, surrounded by the regions dominating the light in the rest-UV. In this paper, we investigate the internal structure of 12 of the most luminous and massive $z\simeq 6-8$ galaxies in the portion of the EGS field observed with recent {\it JWST}/NIRCam imaging. The systems appear clumpy in the rest-UV, with more than half of the light coming from $\simeq 10^7$ to 10$^{9}$ M$_\odot$ star forming complexes that are $\simeq 150$ - 480 pc in size. Multiple clumps are found in individual galaxies with separations of 0.3 to 4.3 kpc. The clumps tend to be dominated by young stars (median = 23 Myr), but we find large variations in the age of clumps within individual galaxies. The [OIII]+H$\beta$ EW varies significantly across different clumps in galaxies (reflecting differences in stars and gas properties), but the HII regions largely track the UV-bright complexes. We do not find older (and redder) nuclear stellar components that were previously undetected or faint in the UV. Perhaps surprisingly, the rest-optical continuum is just as clumpy as the UV, with emission mostly dominated by the bright star forming complexes. The majority of the stellar mass in bright $6<z<8$ galaxies appears to be contained in the $\gtrsim 150$ pc-scale clumpy star forming complexes, reflecting the very active phase of assembly that is common in reionization-era galaxies.

We examine the stability of hierarchical triple systems using direct $N$-body simulations without adopting a secular perturbation approximation. We estimate their disruption timescales in addition to the mere stable/unstable criterion, with particular attention to the mutual inclination between the inner and outer orbits. First, we improve the fit to the dynamical stability criterion by \citet{Mardling1999,Mardling2001} widely adopted in the previous literature. Especially, we find that that the stability boundary is very sensitive to the mutual inclination; coplanar retrograde triples and orthogonal triples are much more stable and unstable, respectively, than coplanar prograde triples. Next, we estimate the disruption timescales of triples satisfying the stability condition up to $10^9$ times the inner orbital period. The timescales follow the scaling predicted by \citet{Mushkin2020}, especially at high $e_\mathrm{out}$ where their random walk model is most valid. We obtain an improved empirical fit to the disruption timescales, which indicates that the coplanar retrograde triples are significantly more stable than the previous prediction. We furthermore find that the dependence on the mutual inclination can be explained by the energy transfer model based on a parabolic encounter approximation. We also show that the disruption timescales of triples are highly sensitive to the tiny change of the initial parameters, reflecting the genuine chaotic nature of the dynamics of those systems.

Questions regarding how primordial or pristine the comets of the solar system are have been an ongoing controversy. In this review, we describe comets' physical evolution from dust and ice grains in the solar nebula to the contemporary small bodies in the outer solar system. This includes the phases of dust agglomeration, the formation of planetesimals, their thermal evolution and the outcomes of collisional processes. We use empirical evidence about comets, in particular from the Rosetta Mission to comet 67P/Churyumov--Gerasimenko, to draw conclusions about the possible thermal and collisional evolution of comets.

We consider the chaotic motion of low-mass bodies in two-body high-order mean-motion resonances with planets in model planetary systems, and analytically estimate the Lyapunov and diffusion timescales of the motion in multiplets of interacting subresonances corresponding to the mean-motion resonances. We show that the densely distributed (though not overlapping) high-order mean-motion resonances, when certain conditions on the planetary system parameters are satisfied, may produce extended planetary chaotic zones -- "zones of weak chaotization," -- much broader than the well-known planetary connected chaotic zone, the Wisdom gap. This extended planetary chaotic zone covers the orbital range between the 2/1 and 1/1 resonances with the planet. On the other hand, the orbital space inner (closer to the host star) with respect to the 2/1 resonance location is essentially long-term stable. This difference arises because the adiabaticity parameter of subresonance multiplets specifically depends on the particle's orbit size. The revealed effect may control the structure of planetesimal disks in planetary systems: the orbital zone between the 2/1 and 1/1 resonances with a planet should be normally free from low-mass material (only that occasionally captured in the first-order 3/2 or 4/3 resonances may survive); whereas any low-mass population inner to the 2/1 resonance location should be normally long-lived (if not perturbed by secular resonances, which we do not consider in this study).

Multi-planet systems orbiting M dwarfs provide valuable tests of theories of small planet formation and evolution. K2-3 is an early M dwarf hosting three small exoplanets (1.5-2.0 Earth radii) at distances of 0.07-0.20 AU. We measure the high-energy spectrum of K2-3 with HST/COS and XMM-Newton, and use empirically-driven estimates of Ly-alpha and extreme ultraviolet flux. We use EXOFASTv2 to jointly fit radial velocity, transit, and SED data. This constrains the K2-3 planet radii to 4% uncertainty and the masses of K2-3b and c to 13% and 30%, respectively; K2-3d is not detected in RV measurements. K2-3b and c are consistent with rocky cores surrounded by solar composition envelopes (mass fractions of 0.36% and 0.07%), H2O envelopes (55% and 16%), or a mixture of both. However, based on the high-energy output and estimated age of K2-3, it is unlikely that K2-3b and c retain solar composition atmospheres. We pass the planet parameters and high-energy stellar spectrum to atmospheric models. Dialing the high-energy spectrum up and down by a factor of 10 produces significant changes in trace molecule abundances, but not at a level detectable with transmission spectroscopy. Though the K2-3 planets span the small planet radius valley, the observed system architecture cannot be readily explained by photoevaporation or core-powered mass loss. We instead propose 1) the K2-3 planets are all volatile-rich, with K2-3d having a lower density than typical of super-Earths, and/or 2) the K2-3 planet architecture results from more stochastic processes such as planet formation, planet migration, and impact erosion.

The Cadmium Zinc Telluride Imager (CZTI) onboard AstroSat, an open detector above $\sim$100 keV, is a promising tool for the investigation of hard X-ray characteristics of $\gamma$-ray pulsars. A custom algorithm has been developed to detect pulsars from long integration ($\sim$years) of archival data, as reported by us earlier. Here we extend this method to include in the analysis an additional $\sim$20% of the CZTI pixels that were earlier ignored due to their lower gain values. Recent efforts have provided better and more secure calibration of these pixels, demonstrating their higher thresholds and extended energy range up to $\sim$1 MeV. Here we use the additional information provided by these pixels, enabling the construction of pulse profiles over a larger energy range. We compare the profiles of the Crab pulsar at different sub-bands and show that the behaviour is consistent with the extended energy coverage. As detailed spectroscopy over this full band remains difficult due to the limited count rate, we construct hardness ratios which, together with AstroSat Mass Model simulations, are able to constrain the power-law index of the radiation spectrum. We present our results for the phase-resolved spectrum of PSR J0534+2200 and for the total pulsed emission of PSR J1513-5908. The recovered photon indices are found to be accurate to within $\sim 20$%.

Extreme BL Lacs (EHBL) form a subclass of blazars which challenge standard emission scenarios. In a recent study it has been argued that their peculiar properties can be explained if emitting electrons are accelerated in a series of oblique shocks induced by the recollimation of the relativistic jet. However, new 3D simulations of recollimated, weakly magnetized jets reveal that, in correspondence with the first recollimation shock, the flow develops a rapidly growing instability, becomes highly turbulent and decelerates, effectively hampering the formation of the multiple shock structure routinely observed in 2D simulations. Building on these new findings, we propose here a revised scenario for EHBL, in which the emission is produced by electrons accelerated at the recollimation shock and subsequently further energized through stochastic acceleration in the turbulent downstream flow. We apply a simple version of this scenario to the prototypical EHBL 1ES 0229-200, showing that the SED can be satisfactorily reproduced with standard values of the main physical parameters.

Stellar rotation is of key importance for the formation process, evolution, and final fate of massive stars. In this paper we review results from the study of the spin rate properties of a sample of more than 400 Galactic O-type stars surveyed by the IACOB and OWN projects. By combining vsini, Teff, and logg estimates (resulting from a detailed quantitative spectroscopic analysis) with information about the spectroscopic binarity status for an important fraction of the stars in the sample, we provide a renewed overview about how the empirical distribution of projected rotational velocities in the O-star domain depends on mass, evolutionary and binary status. The obtained distributions are then compared with predictions of several state-of-the-art evolutionary models for single stars, as well as from population synthesis simulations including binary interaction, and used to provide hints about the initial velocity distribution of stars with masses in the range ~15-80 Msol.

In close exoplanetary systems, tidal interactions drive orbital and spin evolution of planets and stars over long timescales. Tidally-forced inertial waves (restored by the Coriolis acceleration) in the convective envelopes of low-mass stars and giant gaseous planets contribute greatly to the tidal dissipation when they are excited and subsequently damped (e.g. through viscous friction), especially early in the life of a system. These waves are known to be subject to nonlinear effects, including triggering differential rotation in the form of zonal flows. In this study, we use a realistic tidal body forcing to excite inertial waves through the residual action of the equilibrium tide in the momentum equation for the waves. By performing 3D nonlinear hydrodynamical simulations in adiabatic and incompressible convective shells, we investigate how the addition of nonlinear terms affects the tidal flow properties, and the energy and angular momentum redistribution. In particular, we identify and justify the removal of terms responsible for unphysical angular momentum evolution observed in a previous numerical study. Within our new set-up, we observe the establishment of strong cylindrically-sheared zonal flows, which modify the tidal dissipation rates from prior linear theoretical predictions. We demonstrate that the effects of this differential rotation on the waves neatly explains the discrepancies between linear and nonlinear dissipation rates in many of our simulations. We also highlight the major role of both corotation resonances and parametric instabilities of inertial waves, which are observed for sufficiently high tidal forcing amplitudes or low viscosities, in affecting the tidal flow response.

Severe geomagnetic storms appear to be ordered by the solar cycle in a number of ways. They occur more frequently close to solar maximum and declining phase, are more common in larger solar cycles and show different patterns of occurrence in odd- and even-numbered solar cycles. Our knowledge of the most extreme space weather events, however, comes from the spikes in cosmogenic-isotope ($^{14}$C, $^{10}$Be and $^{36}$Cl) records that are attributed to significantly larger solar energetic particle (SEP) events than have been observed during the space age. Despite both storms and SEPs being driven by solar eruptive phenomena, the event-by-event correspondence between extreme storms and extreme SEPs is low. Thus it should not be assumed a priori that the solar cycle patterns found for storms also hold for SEPs and the cosmogenic-isotope events. In this study we investigate the solar cycle trends in the timing and magnitude of the 67 SEP ground-level enhancements (GLEs) recorded by neutron monitors since the mid 1950s. Using a number of models of GLE occurrence probability, we show that GLEs are around a factor four more likely around solar maximum than around solar minimum, and that they preferentially occur earlier in even-numbered solar cycles than in odd-numbered cycles. There are insufficient data to conclusively determine whether larger solar cycles produce more GLEs. Implications for putative space-weather events in the cosmogenic-isotope records are discussed. We find that GLEs tend to cluster within a few tens of days, likely due to particularly productive individual active regions, and with approximately 11-year separations, owing to the solar cycle ordering. But these timescales do not explain cosmogenic-isotope spikes which require multiple extreme SEP events over consecutive years.

We present a spectroscopic analysis of the GIRAFFE and UVES data collected by the Gaia-ESO survey for the young open cluster NGC 3293. Archive spectra from the same instruments obtained in the framework of the `VLT-FLAMES survey of massive stars' are also analysed. Atmospheric parameters, non-LTE chemical abundances for six elements, or variability information are reported for a total of about 160 B stars spanning a wide range in terms of spectral types (B1 to B9.5) and rotation rate (up to 350 km/s). We take advantage of the multi-epoch observations to detect several binary systems or intrinsically line-profile variables. A deconvolution algorithm is used to infer the current, true (deprojected) rotational velocity distribution. We find a broad, Gaussian-like distribution peaking around 200-250 km/s. Although some stars populate the high-velocity tail, most stars in the cluster appear to rotate far from critical. We discuss the chemical properties of the cluster, including the low occurrence of abundance peculiarities in the late B stars and the paucity of objects showing CN-cycle burning products at their surface. We argue that the former result can largely be explained by the inhibition of diffusion effects because of fast rotation, while the latter is generally in accord with the predictions of single-star evolutionary models under the assumption of a wide range of initial spin rates at the onset of main-sequence evolution. However, we find some evidence for a less efficient mixing in two quite rapidly rotating stars that are among the most massive objects in our sample. Finally, we obtain a cluster age of ~20 Myrs through a detailed, star-to-star correction of our results for the effect of stellar rotation. This is significantly older than previous estimates from turn-off fitting that fully relied on classical, non-rotating isochrones. [abridged]

Differential settling and growth of dust grains impact the structure of the radiative envelopes of gaseous planets during formation. Sufficiently rapid dust growth can result in envelopes with substantially reduced opacities for radiation transport, thereby facilitating planet formation. We revisit the problem and establish that dust settling and grain growth also lead to outer planetary envelopes that are prone to compositional instabilities, by virtue of their inverted mean-molecular weight gradients. Under a variety of conditions, we find that the radiative envelopes of forming planets experience compositional turbulence driven by a semi-transparent version of the thermohaline instability ('fingering convection'). The compositional turbulence seems efficient at mixing dust in the radiative envelopes of planets forming at super-AU distances (say $5$ AU) from a Sun-like star, but not so at sub-AU distances (say $0.2$ AU). Furthermore, compositional layering is favoured only at large (super-AU) distances. Such distinct turbulent regimes for planetary envelopes growing at sub-AU vs. super-AU distances could leave an imprint on the final planets formed.

This work uses multiscale environments and structures of galaxies in the Sloan Digital Sky Survey as consistency checks of the evolution from starburst to quiescence at redshift $z < 0.2$. The environmental indicators include fixed aperture mass overdensities ($\delta_{x\mathrm{Mpc}}$, $x \in \{0.5, 1, 2, 4, 8\}\,h^{-1}$Mpc), $k$-nearest neighbor distances, the tidal parameter, halo mass ($M_h$), and satellite/central classification. The residuals of specific star formation rates ($\Delta\,\mathrm{SSFR}$) is used to select starbursts ($\Delta\,\mathrm{SSFR} > 0.6\,$dex, $N \approx 8,\,600$). Quenched post-starbursts (QPSBs) are selected using H$\alpha < 3\,$angstrom in emission and H$\delta_A > 4\,$ angstrom in absorption ($N \approx 750$). The environments of starbursts and QPSBs are compared with those of active galactic nuclei (AGNs) and inactive galaxies of varying $\Delta\,\mathrm{SSFR}$. The environments of starbursts, AGNs, and QPSBs are unlike the environments of most quiescent galaxies (QGs). About $70\%-90\%$ of starbursts, AGNs with H$\delta_A > 4$, and QPSBs are centrals, $\sim 80\%-90\%$ have $M_h < 10^{13}\,M_\odot$, and only $\sim 2\%-4\%$ have $M_h > 10^{14}\,M_\odot$ or live in clusters. Their $M_h$ and satellite fractions are also different from those of QGs. All QPSBs are matched to some SFGs, starbursts, AGNs, and QGs of similar $M_\star$, environments, concentration indices, and velocity dispersions. A significant fraction ($\sim 20\%-30\%$) of starbursts cannot be matched to QPSBs or QGs. The implications are: (1) some starbursts do not quench rapidly. (2) Satellite-quenching mechanisms operating in high density environments cannot account for most QPSBs. (3) The evolution from starbursts to QPSBs to QGs is not the dominant path at $z < 0.2$. (4) Starbursts are not mainly triggered by tidal interactions.

We present the evolutions of the C-BLUE One family of cameras (formerly introduced as C-MORE), a laser guide star oriented wavefront sensor camera family. Within the Opticon WP2 european funded project, which has been set to develop LGS cameras, fast path solutions based on existing sensors had to be explored to provide working-proven cameras to ELT projects ready for the first light schedule. Result of this study, C-BLUE One is a CMOS based camera with 1600x1100 pixels (9um pitch) and 481 FPS refresh rate. It has been developed to answer most of the needs of future laser based adaptive optics systems (LGS) to be deployed on 20-40m-class telescopes as well as on smaller ones. We present the main features of the camera and measured performance in terms of noise, dark current, quantum efficiency and image quality which are the key parameters for the application. The camera has been declined also in fast smaller format (800x600x1500FPS) and large format (3200X2200x250FPS) to cover most of the AO applications.

We introduce the radiative transfer code Sweep for the cosmological simulation suite Arepo. Sweep is a discrete ordinates method in which radiation transport is carried out by a transport sweep across the entire computational grid. Since Arepo is based on an adaptive, unstructured grid, the dependency graph induced by the sweep dependencies of the grid cells is non-trivial. In order to solve the topological sorting problem in a distributed manner, we employ a task-based-parallelism approach. The main advantage of the sweep method is that the computational cost scales only with the size of the grid, and is independent of the number of sources or the distribution of sources in the computational domain, which is an advantage for radiative transfer in cosmological simulations, where there are large numbers of sparsely distributed sources. We successfully apply the code to a number of physical tests such as the expansion of HII regions, the formation of shadows behind dense objects, the scattering of light, as well as its behavior in the presence of periodic boundary conditions. In addition, we measure its computational performance with a focus on highly parallel, large-scale simulations.

Despite two decades of studies, it is still not clear whether ULX spectral transitions are due to stochastic variability in the wind or variations in the accretion rate or in the source geometry. The compact object is also unknown for most ULXs. In order to place constraints onto such scenarios and on the structure of the accretion disc, we studied the temporal evolution of the spectral components of the variable source NGC 55 ULX-1. Using recent and archival data obtained with the XMM-Newton satellite, we modelled the spectra with two blackbody components which we interpret as thermal emission from the inner accretion flow and the regions around or beyond the spherization radius. The luminosity-temperature (L-T) relation of each spectral component agrees with the L proportional T^4 relationship expected from a thin disc model, which suggests that the accretion rate is close to the Eddington limit. However, there are some small deviations at the highest luminosities, possibly due to an expansion of the disc and a contribution from the wind at higher accretion rates. Assuming that such deviations are due to the crossing of the Eddington or supercritical accretion rate, we estimate a compact object mass of 6-14 Msun, favouring a stellar-mass black hole as the accretor.

We present a spectropolarimetric analysis of the hot star V352Peg. We have acquired 18 spectropolarimetric observations of the star with ESPaDOnS at the CFHT between 2018 and 2019 and completed our dataset with one archival ESPaDOnS measurement obtained in 2011. Our analysis of the spectra shows that the star is on the main sequence and chemically peculiar, i.e. it is a Bp star, with overabundances of iron peak elements (Ti, Cr and Fe) and underabundance of He and O. Through a Least-Square Deconvolution of each spectrum, we extracted the mean Zeeman signature and mean line profile of thousands of spectral lines and detected a magnetic field in V352Peg. By modelling the Stokes I and V profiles and using the Oblique Rotator Model, we determined the geometrical configuration of V352Peg. We also performed Zeeman-Doppler Imaging (ZDI) to provide a more detailed characterization of the magnetic field of V352Peg and its surface chemical distributions. We find a magnetic field that is mainly dipolar, dominantly poloidal, and largely non-axisymmetric with a dipole field strength of $\sim$9 kG and a magnetic axis almost perpendicular to the rotation axis. The strong variability of Stokes I profiles also suggests the presence of chemical spots at the stellar surface.

While previous work explored the fractality and self-organized criticality (SOC) of flares and nanoflares in wavelengths emitted in the solar corona (such as in hard X-rays, soft X-rays, and EUV wavelenghts), we focus here on impulsive phenomena in the photosphere and transition region, as observed with the {\sl Interface Region Imaging Spectrograph (IRIS)} in the temperature range of $T_e \approx 10^4-10^6$ K. We find the following fractal dimensions (in increasing order): $D_A=1.21 \pm 0.07$ for photospheric granulation, $D_A=1.29 \pm 0.15$ for plages in the transition region, $D_A=1.54 \pm 0.16$ for sunspots in the transition region, $D_A=1.59 \pm 0.08$ for magnetograms in active regions, $D_A=1.56 \pm 0.08$ for EUV nanoflares, $D_A=1.76 \pm 0.14$ for large solar flares, and up to $D_A=1.89 \pm 0.05$ for the largest X-class flares. We interpret low values of the fractal dimension ($1.0 \lapprox D_A \lapprox 1.5$) in terms of sparse curvi-linear flow patterns, while high values of the fractal dimension ($1.5 \lapprox D_A \lapprox 2.0$) indicate near space-filling transport processes, such as chromospheric evaporation. Phenomena in the solar transition region appear to be consistent with SOC models, based on their size distributions of fractal areas $A$ and (radiative) energies $E$, which show power law slopes of $\alpha_A^{obs}=2.51 \pm 0.21$ (with $\alpha_A^{theo}=2.33$ predicted), and $\alpha_E^{obs}=2.03 \pm 0.18$ (with $\alpha_E^{theo}=1.80$ predicted).

The aim of this study is to examine dust dynamics on a large scale and investigate the coupling of dust with gas fluid in the star formation process. We propose a method for calculating the dust trajectory in a gravitationally collapsing cloud, where the dust grains are treated as Lagrangian particles and are assumed to be neutral. We perform the dust trajectory calculations in combination with non-ideal magnetohydrodynamics simulation. Our simulation shows that dust particles with a size of $\le 10\,{\rm \mu m}$ are coupled with gas in a star-forming cloud core. We investigate the time evolution of the dust-to-gas mass ratio and the Stokes number, which is defined as the stopping time normalized by the freefall time-scale, and show that large dust grains ($\gtrsim 100\,{\rm \mu m}$) have a large Stokes number (close to unity) and tend to concentrate in the central region (i.e., protostar and rotationally supported disk) faster than do small grains ($\lesssim 10\,{\rm \mu m}$). Thus, large grains significantly increase the dust-to-gas mass ratio around and inside the disk. We also confirm that the dust trajectory calculations, which trace the physical quantities of each dust particle, reproduce previously reported results obtained using the Eulerian approach.

We have investigated the evolution of a system of sparking discharges in the inner acceleration region (IAR) above the pulsar polar cap. The surface of the polar cap is heated to temperatures around $10^6$ K and forms a partially screened gap (PSG) due to thermionic emission of positively charged ions from the stellar surface. The sparks lag behind the co-rotation speed during their lifetimes due to variable $E$x$B$ drift. In a PSG the sparking discharges arise in locations where the surface temperatures go below the critical level ($T_i$) for ions to freely flow from the surface. The sparking commences due to the large potential drop developing along the magnetic field lines in these lower temperature regions and subsequently the back streaming particles heat the surface to $T_i$. The temperature regulation requires the polar cap to be tightly filled with sparks and a continuous presence of sparks is required around its boundary since no heating is possible from the closed field line region. We have estimated the time evolution of the sparking system in the IAR which shows a gradual shift in the spark formation along two distinct directions resembling clockwise and anti-clockwise motion in two halves of the polar cap. Due to the differential shift of the sparking pattern in the two halves, a central spark develops representing the core emission. The temporal evolution of the sparking process was simulated for different orientations of the non-dipolar polar cap and reproduced the diverse observational features associated with subpulse drifting.

We search for detectable signatures of f(R) gravity and its chameleon screening mechanism in the baryonic and dark matter (DM) properties of simulated void galaxies. The enhancement of the gravitational acceleration can have a meaningful impact on the scaling relations as well as on the halo morphology. The galaxy rotational velocity field (calculated with the velocity of the gas disc and the acceleration fields) deviates from the typical values of the Tully-Fisher Relation (TFR) in GR. For a given stellar mass, f(R) gravity tends to produce greater maximum velocities. On the other hand, the mass in haloes in f(R) gravity is more concentrated than their counterparts in GR. This trend changes when the concentration is calculated with the dynamical density profile, which takes into account the unscreened outer regions of the halo. Stellar discs interact with the overall potential well in the central regions, modifying the morphology of the screening regions and reshaping them. We find a trend for galaxies with a more dominant stellar disc to deviate further from round screening regions. We find that small haloes are less triaxial and more round in f(R) than their GR counterparts. The difference between halo morphology becomes smaller in f(R) haloes whose inner regions are screened. These results suggest possible observables that could unveil modified gravity effects on galaxies in voids in future cosmological tests of gravity.

The goal of our study is to provide a reliable classification of variability of A-F stars brighter than 11 mag located in the northern TESS continuous viewing zone. We also aim at thorough discussion about issues in the classification related to the data characteristics and the issues arising from the similar light curve shape generated by different physical mechanisms. We used TESS long- and short-cadence photometric data and corresponding Fourier transform to classify the variability type of the stars. We present a clear and concise classification system that is demonstrated on many examples. We found clear signs of variability in 3025 of 5923 studied stars (51 %). For 1813 of these 3025 stars, we provide a classification. From the classified stars, 64.5 % are pulsating stars of GDOR and DSCT types and their hybrids. We realised that the long- and short-cadence PDCSAP data can differ significantly not only in amplitudes but also in the content of instrumental/data reduction artifacts making the long-cadence data less reliable. We identified a new group of stars showing stable light curves and characteristic frequency spectra pattern (8.5 % of the classified stars). According to the position in the Hertzsprung-Russell diagram, these stars are likely GDOR stars but are about 200 K cooler than GDORs on average and have smaller amplitudes and longer periods in average. We show that without spectroscopic observations, it can be impossible to unambiguously distinguish between ellipsoidal variability and rotational variability. We also apply our methodology to three previous studies and found significant discrepancies in the classification. We demonstrate how difficult the classification of variable A-F stars can be when using only photometric data.

Scientific writing is an important skill for a career as a professional astrophysicist. Very few researchers, however, receive any formal training in how to write scientific research papers of high quality in an efficient manner. This paper (Paper I) is the first of a two-part self-help guide in scientific writing to address this skills gap. Paper I focuses on planning your academic research paper in astronomy. We discuss how to crystallise the ideas that underlie a research project, analyse how the paper can be constructed considering the audience and the chosen journal, and give an overview of the publishing process. Paper II is a detailed description of the different sections that make up a research paper in astronomy and shares the best practice in how to write in English. Whether you are a student writing your first paper or an experienced author, you may find the ideas presented here useful.

With the characterisations of potentially habitable planetary atmospheres on the horizon, the search for biosignatures is set to become a major area of research in the coming decades. To understand the atmospheric characteristics that might indicate alien life we must understand the abiotic characteristics of a planet and how life interacts with its environment. In the field of biogeochemistry, sophisticated models of life-environment coupled systems demonstrate that many assumptions specific to Earth-based life, e.g. specific ATP maintenance costs, are unnecessary to accurately model a biosphere. We explore a simple model of a single-species microbial biosphere that produces CH4 as a byproduct of the microbes' energy extraction - known as a type I biosignature. We demonstrate that although significantly changing the biological parameters has a large impact on the biosphere's total population, such changes have only a minimal impact on the strength of the resulting biosignature, while the biosphere is limited by H2 availability. We extend the model to include more accurate microbial energy harvesting and show that adjusting microbe parameters can lead to a regime change where the biosphere becomes limited by energy availability and no longer fully exploits the available H2, impacting the strength of the resulting biosignature. We demonstrate that, for a nutrient limited biosphere, identifying the limiting nutrient, understanding the abiotic processes that control its abundance, and determining the biosphere's ability to exploit it, are more fundamental for making type I biosignature predictions than the details of the population dynamics of the biosphere.

Polstar is a proposed NASA MIDEX space telescope that will provide high-resolution, simultaneous full-Stokes spectropolarimetry in the far ultraviolet, together with low-resolution linear polarimetry in the near ultraviolet. This observatory offers unprecedented capabilities to obtain unique information on the magnetic and plasma properties of the magnetospheres of hot stars. We describe an observing program making use of the known population of magnetic hot stars to test the fundamental hypothesis that magnetospheres should act to rapidly drain angular momentum, thereby spinning the star down, whilst simultaneously reducing the net mass-loss rate. Both effects are expected to lead to dramatic differences in the evolution of magnetic vs. non-magnetic stars.

The mean parallax of the Pleiades open cluster from the Hipparcos catalog is larger than the true value by approximately 1 mas. The origin of this error, as well as a possible algorithm of correcting it, was proposed by Makarov (2002). The problem is reassessed using the more accurate Gaia data with a focus on the predicted correction to the Pleiades proper motions. The accurately determined differences Gaia - Hipparcos for 52 common stars are close to these estimates within the formal uncertainties for all three parameters, which strongly suggests that the proposed interpretation was correct. With adjustments for the systematic vector field fitted with 126 vector spherical harmonics to degree 7, these differences amount to $(+0.39,\,-0.74)$ mas yr$^{-1}$. The implications of small-scale proper motion and position errors in Hipparcos for present day astrometry are briefly discussed.

Modern astrometric and spectroscopic surveys have revealed a wealth of structure in the phase space of stars in the Milky Way, with evidence of resonance features and non-equilibrium processes. Using Gaia DR3, we present evidence of a new resonance-like feature in the outer disc of the Milky Way. The feature is most evident in the angular momentum distribution of the young Classical Cepheids, a population for which we can derive accurate distances over much of the Galactic disc. We then search for similar features in the outer disc using a much larger sample of red giant stars, as well as a compiled list of over 31 million stars with spectroscopic line-of-sight velocity measurements. While much less evident in these two older samples, the distribution of stars in action-configuration space suggests that resonance features are present here as well. The position of the feature in action-configuration space suggests that the new feature may be related to the Galactic bar, but other possibilities are discussed.

Large constellations of orbiting communication satellites will become an important source of noise for present and future astronomical observatories. Mitigation measures rely on high quality predictive models of the position and expected brightness of these objects. Optical linear imaging polarimetry holds promise as a quantitative tool to improve our understanding of the physics of reflection of sunlight off satellite components and through which models of expected brightness can be improved. We present the first simultaneous short-timescale linear polarimetry and optical photometry observations of a geostationary satellite, using the new MOPTOP imaging polarimeter on the 2m Liverpool Telescope. Our target, telecommunication satellite Thor-6, shows prominent short timescale glint-like features in the lightcurve, some as short as seconds. Our polarimetric observations overlap with several of these micro-glints, and have the cadence required to resolve them. We find that the polarisation lightcurve is remarkably smooth, the short time scale glints are not seen to produce strong polarimetric features in our observation. We show how short timescale polarimetry can further constrain the properties of the components responsible for these micro-glints.

In the past decade, ALMA observations have revealed that a large fraction of protoplanetary discs contains rings in the dust continuum. These rings are the locations where pebbles accumulate, which is beneficial for planetesimal formation and subsequent planet assembly. We investigate the viability of planet formation inside ALMA rings in which pebbles are trapped by either a Gaussian-shape pressure bump or by the strong dust backreaction. Planetesimals form at the midplane of the ring via streaming instability. By conducting N-body simulations, we study the growth of these planetesimals by collisional mergers and pebble accretion. Thanks to the high concentration of pebbles in the ring, the growth of planetesimals by pebble accretion becomes efficient as soon as they are born. We find that type-I migration plays a decisive role in the evolution of rings and planets. For discs where planets can migrate inward from the ring, a steady state is reached where the ring spawns ${\sim}20 M_\oplus$ planetary cores as long as rings are fed with materials from the outer disc. The ring acts as a long-lived planet factory and it can explain the 'fine-tuned' optical depths of the observed dust rings in the DSHARP large program. In contrast, in the absence of a planet removal mechanism (migration), a single massive planet will form and destroy the ring. A wide and massive planetesimals belt will be left at the location of the planet-forming ring. Planet formation in rings may explain the mature planetary systems observed inside debris discs.

Interacting dark energy models have been suggested as alternatives to the standard cosmological model, $\Lambda$CDM. We focus on a phenomenologically interesting class of dark scattering models that is characterised by pure momentum exchange between dark energy and dark matter. This model extends the parameter space with respect to $\Lambda$CDM by two parameters, $w$ and $A$, which define the dark energy equation of state and the strength of the coupling between dark energy and dark matter, respectively. In order to test non-standard cosmologies with Stage-IV galaxy clustering surveys, it is crucial to model mildly nonlinear scales and perform precision vs accuracy tests. We use the Effective Field Theory of Large-Scale Structure, and we perform validation tests by means of an MCMC analysis using a large set of N-body simulations. We find that adding the bispectrum monopole to the power spectrum multipoles improves the constraints on the dark energy parameters by $\sim 30 \%$ for $k_{\mathrm{max}, B}^{l=0} = 0.11$ $h$ Mpc$^{-1}$ without introducing biases in the parameter estimation. We also find that the same improvement can be achieved with more moderate scale cuts and the use of bias relations, or with the addition of the bispectrum quadrupole. Finally, we forecast constraints assuming Stage-IV surveys specifications at $z=1$ in the interacting dark energy scenario: $\sigma_w = 0.08$ and $\sigma_A = 2.51$ b GeV$^{-1}$, as well as for $w$CDM cosmology: $\sigma_w = 0.1$.

The Early Release Observations (EROs) of JWST beautifully demonstrate the promise of JWST in characterizing the universe at cosmic dawn. We analyze the ERO spectra of three $z \sim 8$ galaxies to determine their metallicities, gas temperatures and ionization. These galaxies offer the first opportunity to understand the physical properties of epoch-of-reionization galaxies through detailed rest-optical emission line spectroscopy. We show that these objects have metal abundances $12+\log[O/H] \approx 6.9 - 8.2$, based on both the $T_e$ method and on a recent calibration of the $R_{23}$ metallicity indicator. Since the spectra are some of the earliest science data from JWST, we compare several line ratios with values expected from robust physics, to validate our measurement procedures. We compare the abundances and emission line ratios to a nearby sample of Green Pea galaxies -- a population of nearby emission line galaxies whose UV properties resemble epoch-of-reionization galaxies, and which often have large Lyman continuum escape fractions. The JWST data show striking further similarities between these high redshift galaxies and nearby Green Peas. The $z\sim 8$ galaxies span the metallicity range covered by Green Peas. They also show the compact morphology that is typical of emission line dominated galaxies at all redshifts. Based on these similarities with Green Peas, it is likely that these are the first rest-optical spectra of galaxies that are actively driving cosmological reionization

We present a low-cost ultraviolet to infrared absolute quantum efficiency detector characterization system developed using commercial off-the-shelf components. The key components of the experiment include a light source,a regulated power supply, a monochromator, an integrating sphere, and a calibrated photodiode. We provide a step-by-step procedure to construct the photon and quantum efficiency transfer curves of imaging sensors. We present results for the GSENSE 2020 BSI CMOS sensor and the Sony IMX 455 BSI CMOS sensor. As a reference for similar characterizations, we provide a list of parts and associated costs along with images of our setup.

The location of $\gamma$-ray emitting region in blazars has been an open issue for several decades and is still being debated. We use the Paliya et al. sample of 619 $\gamma$-ray-loud flat-spectrum radio quasars with the available spectral energy distributions, and employ a seed photon factor approach, to locate the $\gamma$-rays production region. This method efficiently set up a relation between the peak frequencies and luminosities for the synchrotron emission and inverse Compton scattering, together with a combination of the energy density and characteristic energy for the external seed photon field, namely, $\sqrt{U_0}/\epsilon_0$, an indicative factor of seed photons (SF) in units of Gauss. By means of comparing it with canonical values of broad-line region and molecular dusty torus, we principally ascertain that the GeV emission is originated far beyond the BLR and close to the DT -- farther out at pc scales from the central black hole, which supports a {\it far-site} scenario for $\gamma$-ray blazars. We probe the idea that inverse Compton scattering of infrared seed photons is happening in the Thomson regime. This approach and our findings are based on the validity of the External Compton model, which is applicable to understand the GeV emission mechanism in FSRQs. However, the completeness of this framework has been challenged by reports of neutrino emission from blazars. Thus we also shed new light on the neutrino production region by using our derived results since blazars are promising neutrino emitters.

The James Webb Space Telescope (JWST) Early Release Observations (EROs) is a set of public outreach products created to mark the end of commissioning and the beginning of science operations for JWST. Colloquially known as the ``Webb First Images and Spectra", these products were intended to demonstrate to the worldwide public that JWST is ready for science, and is capable of producing spectacular results. The package was released on July 12, 2022, and included images and spectra of the galaxy cluster SMACS~J0723.3-7327 and distant lensed galaxies, the interacting galaxy group Stephan's Quintet, NGC 3324 in the Carina star-forming complex, the Southern Ring planetary nebula NGC 3132, and the transiting hot Jupiter WASP 96b. This paper describes the ERO technical design, observations, and scientific processing of data underlying the colorful outreach products.

The sharpest optical images of the R136 cluster in the Large Magellanic Cloud are presented, allowing for the first time to resolve members of the central core, including R136a1, the most massive star known. These data were taken using the Gemini speckle imager Zorro in medium-band filters with effective wavelengths similar to BVRI achieving angular resolutions between 30-40 mas. All stars previously known in the literature, having $V<$16 mag within the central $2''\times2''$ were recovered. Visual companions ($\geq$40 mas; 2000 au) were detected for the WN5h stars R136 a1, and a3. Photometry of the visual companion of a1 suggests it is of mid O spectral type. Based on new photometric luminosities using the resolved Zorro imaging, the masses of the individual WN5h stars are estimated to be between 150-200 $M_{\odot}$, lowering significantly the present-day masses of some of the most massive stars known. These mass estimates are critical anchor points for establishing the stellar upper-mass function.

We present a kinetically mixed dark sector (KMIX) model to address the Hubble and $S_8$ tensions. Inspired from string theory, our model includes two fields: an axion, which plays a role similar to the scalar field in early dark energy models, and a dilaton. This theory differs from other axio-dilaton models aimed at the Hubble tension in that there is necessarily kinetic mixing between the two fields which allows for efficient energy transfer from the axion into the dilaton which has $w\approx1$. As a direct consequence of these dynamics, we find the model does not need to resort to a fine-tuned potential to solve the Hubble tension and naturally accommodates a standard axion potential. Furthermore, the axion will necessarily makeup a small (fuzzy) fraction of $\Omega_{\rm cdm}$ once it begins to oscillate at the bottom of its potential and will suppress the growth of perturbations on scales sensitive to $S_8$. Interestingly, the scale of the potential for the dilaton has to be small, $\lesssim \mathcal{O}(10~{\rm meV})^4$, suggesting the possibility for a connection to dark energy. Implementing the dynamics for the background and perturbations in a modified Boltzmann code we calculate the CMB and matter power spectra for our theory. Exploring the parameter space of our model, we find regions which can accommodate a $\sim 10\%$ increase in $H_0$ from the Planck inferred value and $S_8$ values that are consistent with large-scale structure constraints.

Conventional planet formation theory suggests that chondritic materials have delivered crucial atmospheric and hydrospheric elements such as carbon (C), nitrogen (N), and hydrogen (H) onto primitive Earth. However, recent measurements highlight the significant elemental ratio discrepancies between terrestrial parent bodies and the supposed planet building blocks. Here we present a volatile evolution model during the assembly of Earth and Earth-like planets. Our model includes impact losses, atmosphere-mantle exchange, and time dependent effects of accretion and outgassing calculated from dynamical modeling outcomes. Exploring a wide range of planetesimal properties (i.e., size and composition) as well as impact history informed by N-body accretion simulations, we find that while the degree of CNH fractionation has inherent stochasticity, the evolution of C/N and C/H ratios can be traced back to certain properties of the protoplanet and projectiles. Interestingly, the majority of our Earth-like planets acquire superchondritic final C/N ratios, implying that the volatile elemental ratios on terrestrial planets are driven by the complex interplay between delivery, atmospheric ablation, and mantle degassing.

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they undergo core-collapse supernovae (CCSNe). The resulting dense circumstellar medium (CSM) may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomena were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2-4 km/s convective motions seed the intrinsic convective instability associated with the high entropy material deep in the envelope, enabling gas from deep within the envelope to escape, increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by 10x along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

We try to find conditions, the fulfillment of which allows a universe born in a metastable false vacuum state to survive and not to collapse. The conditions found are in the form of inequalities linking the depending on time $t$ instantaneous decay rate ${\it\Gamma}(t)$ of the false vacuum state and the Hubble parameter $H(t)$. Properties of the decay rate of a quantum metastable states are discussed and then the possible solutions of the conditions found are analyzed and discussed. Within the model considered it is shown that a universe born in the metastable vacuum state has a very high chance of surviving until very late times if the lifetime, $\tau_{0}^{F}$, of the metastable false vacuum state is comparable to (or even shorter), than time when the inflation process begins. Our analysis shows that the instability of the electroweak vacuum does not have to result in the tragic fate of our Universe leading to its death.

The Scintillating Bubble Chamber (SBC) Collaboration is developing liquid-noble bubble chambers for the quasi-background-free detection of low-mass (GeV-scale) dark matter and coherent scattering (CE$\nu$NS) of low-energy (MeV-scale) neutrinos. The first physics-scale demonstrator of this technique, a 10-kg liquid argon bubble chamber dubbed SBC-LAr10, is now being commissioned at Fermilab. This device will calibrate the background discrimination power and sensitivity of superheated argon to nuclear recoils at energies down to 100 eV. A second functionally-identical detector with a focus on radiopure construction is being built for SBC's first dark matter search at SNOLAB. The projected spin-independent sensitivity of this search is approximately $10^{-43}$ cm$^2$ at 1 GeV$/c^2$ dark matter particle mass. The scalability and background discrimination power of the liquid-noble bubble chamber make this technique a compelling candidate for future dark matter searches to the solar neutrino fog at 1 GeV$/c^2$ particle mass (requiring a $\sim$ton-year exposure with non-neutrino backgrounds sub-dominant to the solar CE$\nu$NS signal) and for high-statistics CE$\nu$NS studies at nuclear reactors.

A developing supercritical collisionless shock propagating in a homogeneously magnetized plasma of ambient gas origin having higher uniformity than the previous experiments is formed by using high-power laser experiment. The ambient plasma is not contaminated by the plasma produced in the early time after the laser shot. While the observed developing shock does not have stationary downstream structure, it possesses some characteristics of a magnetized supercritical shock, which are supported by a one-dimensional full particle-in-cell simulation taking the effect of finite time of laser-target interaction into account.

Inverse cascades of kinetic energy and thermal variance in the subset of vertically homogeneous modes in spectral space are found to cause a slow aggregation to a pair of convective supergranules that eventually fill the whole horizontally extended, three-dimensional, turbulent Rayleigh-B\'{e}nard convection layer when a heat flux is prescribed at the top and bottom. An additional weak rotation of the layer around the vertical axis stops this aggregation at a scale that is smaller than the lateral domain extension and ceases the inverse cascade. The resulting characteristic length of the aggregated convection patterns depends on the thermal driving and linearly on the strength of rotation. Our study demonstrates the importance of inverse energy cascades beyond the two-dimensional turbulence case, in a three-dimensional convection flow that is subject to a multi-scale energy injection by thermal plumes and driven by boundary heat fluxes as typically present in natural geo- and astrophysical systems.

Relativistic magnetohydrodynamics (RMHD) provides an extremely useful description of the low-energy long-wavelength phenomena in a variety of physical systems from quark-gluon plasma in heavy-ion collisions to matters in supernovas, compact stars, and early universe. We review the recent theoretical progresses of RMHD, such as a formulation of RMHD from the perspective of magnetic flux conservation using the entropy-current analysis, the nonequilibrium statistical operator approach applied to quantum electrodynamics, and the relativistic kinetic theory. We discuss how the transport coefficients in RMHD are computed in kinetic theory and perturbative quantum field theories. We also explore the collective modes and instabilities in RMHD with a special emphasis on the role of chirality in a parity-odd plasma. We also give some future prospects of RMHD, including the interaction with spin hydrodynamics and the new kinetic framework with magnetic flux conservation.

Recent advances in nuclear theory combined with new astrophysical observations have led to the need for specific theoretical models that actually apply to phenomena on dense-matter physics. At the same time, quantum chromodynamics (QCD) predicts the existence of non-nucleonic degrees of freedom at high densities in neutron-star matter, such as quark matter. Within a confining quark matter model, which consists of homogeneous, neutral 3-flavor interacting quark matter with $\mathcal{O}(m_s^4)$ corrections, we study the structure of compact stars made of a charged perfect fluid in the context of $f(R,T)$ gravity. The system of differential equations that describe the structure of charged compact stars have been derived and solved numerically for a gravity model with $f(R,T)= R+ 2\beta T$. For simplicity, we assume that the charge density is proportional to the energy density, namely, $\rho_{\rm ch} = \alpha \rho$. It is demonstrated that matter-geometry coupling constant $\beta$ and the charge parameter $\alpha$ affect the total gravitational mass and the radius of the star.

The ultra high energy (UHE) cosmic neutrinos are expected to play a pivotal role in the disquisition of physics beyond the standard model of particle physics as well as serve as an ideal cosmic messengers. This epitomizes the selling point of several currently running or planned neutrino telescopes. The UHE cosmic neutrinos usually perambulate gargantuan scales in the extragalactic universe having a magnetic field. If neutrinos have a finite magnetic moment ($\mu_{\nu}$) owing to quantum loop corrections, this may result in spin-flavor oscillations, which can affect the cosmic neutrino flux. Using the current limit, we show that the flux of cosmic neutrinos will reduce by half if they traverse few Mpcs through the intergalactic magnetic field, in the range of $\rm \mu G$ to nG. The condition of averaging out probabilities would not be possible for the neutrino sources within the Milky Way due to insufficient intersteller distances. Moreover, one can safely neglect the effect of $\mu_{\nu}$ if the current upper bound is improved by a few orders of magnitude even if the neutrinos travels through the size of the visible universe.

Generating an ensemble of equations of state that fulfill multimessenger constraints, we statistically determine the properties of dense matter found inside neutron stars (NSs). We calculate the speed of sound and trace anomaly and demonstrate that they are driven towards their conformal values at the center of maximally massive NSs. The local peak of the speed of sound is shown to be located at values of the energy and particle densities which are consistent with deconfinement and percolation conditions in QCD matter. We also analyze fluctuations of the net-baryon number density in the context of possible remnants of critical behavior. We find that the global maxima of the variance of these fluctuations emerge at densities beyond those found in the interiors of NSs.

Given the lack of empirical evidence of weakly interacting dark matter, it is reasonable to look to other candidates such as a confining dark sector with a similar number of particles as the standard model. Twin Higgs mirror matter is one such model that is a twin of the standard model with particles masses 3--6 times heavier than the standard model that solves the hierarchy problem. This generically predicts mirror neutron stars, degenerate objects made entirely of mirror nuclear matter. We find their structure using a realistic equation of state from crust (nuclei) to core (relativistic mean-field model) and scale the particle masses using lattice QCD results. We find that mirror neutron stars have unique signatures that are detectable via gravitational waves and binary pulsars, that provides an intriguing possibility for probing dark matter.