Papers for Tuesday, May 02 2023

Big data has become the norm in astronomy, making it an ideal domain for computer science research. Astronomers typically classify galaxies based on their morphologies, a practice that dates back to Hubble (1936). With small datasets, classification could be performed by individuals or small teams, but the exponential growth of data from modern telescopes necessitates automated classification methods. In December 2013, Winton Capital, Galaxy Zoo, and the Kaggle team created the Galaxy Challenge, which tasked participants with developing models to classify galaxies. The Kaggle Galaxy Zoo dataset has since been widely used by researchers. This study investigates the impact of colour space transformation on classification accuracy and explores the effect of CNN architecture on this relationship. Multiple colour spaces (RGB, XYZ, LAB, etc.) and CNN architectures (VGG, ResNet, DenseNet, Xception, etc.) are considered, utilizing pre-trained models and weights. However, as most pre-trained models are designed for natural RGB images, we examine their performance with transformed, non-natural astronomical images. We test our hypothesis by evaluating individual networks with RGB and transformed colour spaces and examining various ensemble configurations. A minimal hyperparameter search ensures optimal results. Our findings indicate that using transformed colour spaces in individual networks yields higher validation accuracy, and ensembles of networks and colour spaces further improve accuracy. This research aims to validate the utility of colour space transformation for astronomical image classification and serve as a benchmark for future studies.

We develop a neural network based pipeline to estimate masses of galaxy clusters with a known redshift directly from photon information in X-rays. Our neural networks are trained using supervised learning on simulations of eROSITA observations, focusing in this paper on the Final Equatorial Depth Survey (eFEDS). We use convolutional neural networks which are modified to include additional information of the cluster, in particular its redshift. In contrast to existing work, we utilize simulations including background and point sources to develop a tool which is usable directly on observational eROSITA data for an extended mass range from group size halos to massive clusters with masses in between $10^{13}M_\odot<M<10^{15}M_\odot.$ Using this method, we are able to provide for the first time neural network mass estimation for the observed eFEDS cluster sample from Spectrum-Roentgen-Gamma/eROSITA observations and we find consistent performance with weak lensing calibrated masses. In this measurement, we do not use weak lensing information and we only use previous cluster mass information which was used to calibrate the cluster properties in the simulations. When compared to simulated data, we observe a reduced scatter with respect to luminosity and count-rate based scaling relations. We comment on the application for other upcoming eROSITA All-Sky Survey observations.

Radio jets and the lobes they inflate are common in cool-core clusters and are expected to play a critical role in regulating the heating and cooling of the intracluster medium (ICM). This is an inherently multi-scale problem, and much effort has been made to understand the processes governing the inflation of lobes and their impact on the cluster, as well as the impact of the environment on the jet-ICM interaction, on both macro- and microphysical scales. Developments of new numerical techniques and improving computational resources have seen simulations of jet feedback in galaxy clusters become ever more sophisticated. This ranges from modelling ICM plasma physics processes such as the effects of magnetic fields, cosmic rays and viscosity to including jet feedback in cosmologically evolved cluster environments in which the ICM thermal and dynamic properties are shaped by large-scale structure formation. In this review, we discuss the progress made over the last ~decade in capturing both the macro- and microphysical processes in numerical simulations, highlighting both the current state of the field as well as open questions and potential ways in which these questions can be addressed in the future.

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars.

Current theories of dynamical friction on galactic bars are based either on linear perturbation theory, which is valid only in the fast limit where the bar changes its pattern speed rapidly, or on adiabatic theory, which is applicable only in the slow limit where the bar's pattern speed is near-constant. In this paper, we study dynamical friction on galactic bars spinning down at an arbitrary speed, seamlessly connecting the fast and slow limits. We treat the bar-halo interaction as a restricted $N$-body problem and solve the collisionless Boltzmann equation using the angle-averaged Hamiltonian. The phase-space distribution and density wakes predicted by our averaged model are in excellent agreement with full 3D simulations. In the slow regime where resonant trapping occurs, we show that, in addition to the frictional torque, angular momentum is transferred directly due to the migration of the trapped phase-space: trapped orbits comoving with the resonance typically gain angular momentum, while untrapped orbits leaping over the trapped island lose angular momentum. Due to the negative gradient in the distribution function, gainers typically outnumber the losers, resulting in a net negative torque on the perturber. The torque due to the untrapped orbits was identified by Tremaine & Weinberg, who named the phenomenon dynamical feedback. Here we derive the complete formula for dynamical feedback, accounting for both trapped and untrapped orbits. Using our revised formula, we show that dynamical feedback can comprise up to $30\%$ of the total torque on the Milky Way's bar.

We assemble the largest CIV absorption line catalogue to date, leveraging machine learning, specifically Gaussian processes, to remove the need for visual inspection for detecting CIV absorbers. The catalogue contains probabilities classifying the reliability of the absorption system within a quasar spectrum. Our training set was a sub-sample of DR7 spectra that had no detectable CIV absorption in a large visually inspected catalogue. We used Bayesian model selection to decide between our continuum model and our absorption-line models. Using a random hold-out sample of 1301 spectra from all of the 26,030 investigated spectra in DR7 CIV catalogue, we validated our pipeline and obtained an 87% classification performance score. We found good purity and completeness values, both ~80%, when a probability of ~95% is used as the threshold. Our pipeline obtained similar CIV redshifts and rest equivalent widths to our training set. Applying our algorithm to 185,425 selected quasar spectra from SDSS DR12, we produce a catalogue of 113,775 CIV doublets with at least 95% confidence. Our catalogue provides maximum a posteriori values and credible intervals for CIV redshift, column density, and Doppler velocity dispersion. We detect CIV absorption systems with a redshift range of 1.37 $\!-\!$ 5.1, including 33 systems with a redshift larger than 5 and 549 absorbers systems with a rest equivalent width greater than 2 A at more than 95% confidence. Our catalogue can be used to investigate the physical properties of the circumgalactic and intergalactic media.

We present new CO ($J=5-4$ and $7-6$) and [CI] ($^3P_2\,-\, ^3P_1$ and $^3P_1\,-\, ^3P_0$) emission line observations of the star-forming galaxy D49 at the massive end of the Main Sequence at $z=3$. We incorporate previous CO ($J=3-2$) and optical-to-millimetre continuum observations to fit its spectral energy distribution (SED). Our results hint at high-$J$ CO luminosities exceeding the expected location on the empirical correlations with the infrared luminosity. [CI] emission fully consistent with the literature trends is found. We do not retrieve any signatures of a bright active galactic nucleus that could boost the $J=5-4,\,7-6$ lines in either the infrared or X-ray bands, but warm photon-dominated regions, shocks or turbulence could in principle do so. We suggest that mechanical heating could be a favourable mechanism able to enhance the gas emission at fixed infrared luminosity in D49 and other main-sequence star-forming galaxies at high redshift, but further investigation is necessary to confirm this explanation. We derive molecular gas masses from dust, CO, and [CI] that all agree within the uncertainties. Given its large star formation rate (SFR) $\sim 500~M_\odot~{\rm yr}^{-1}$ and stellar mass $>10^{11.5}~M_\odot$, the short depletion time scale of $<0.3$ Gyr might indicate that D49 is experiencing its last growth spurt and will soon transit to quiescence.

While the dipole magnetic field axis of neutron stars is usually postulated to cross the star's centre, it may be displaced from this location, as it has been recently indicated in the millisecond pulsar J0030+0451. Under these conditions, the electromagnetic rocket effect may be activated, where the magnetic field exerts a net force, accelerating the star. This post-natal kick mechanism relies on asymmetric electromagnetic radiation from an off-centre dipole may be relevant to the high spatial velocities of pulsars $\sim 10^{3}$ $\rm{km/s}$. Here, we explore its impact in young pulsars associated with supernova remnants and we compare the observational data on characteristic quantities, such as the braking index and proper motion, with results obtained from the rocket effect. Using a Markov Chain Monte Carlo analysis, we explore the required conditions, for the initial spin periods and the distance between the magnetic axis and the star's center, so that the velocity kick due to the rocket effect approaches the present velocity. We find that the electromagnetic rocket effect can account for typical pulsar transverse velocities assuming an initial spin period of 3.8 $\rm{ms}$ and a dipole field whose distance from the centre of the star is approximately 7 $\rm{km}$. We also explore the influence of the rocket effect on the braking index of a neutron star, and we find that for the sample studied this impact is minimal. Finally, we apply the rocket effect model on the pulsars J0030+0451 and J0538+2817, which are likely candidates for this mechanism.

Gamma-ray bursts (GRBs), both long and short, are explosive events whose inner engine is generally expected to be a black hole or a highly magnetic neutron star (magnetar) accreting high density matter. Recognizing the nature of GRB central engines, and in particular the formation of neutron stars (NSs), is of high astrophysical significance. A possible signature of NSs in GRBs is the presence of a plateau in the early X-ray afterglow. Here we carefully select a subset of long and short GRBs with a clear plateau, and look for an additional NS signature in their prompt emission, namely a transition between accretion and propeller in analogy with accreting, magnetic compact objects in other astrophysical sources. We estimate from the prompt emission the minimum accretion luminosity below which the propeller mechanism sets in, and the NS magnetic field and spin period from the plateau. We demonstrate that these three quantities obey the same universal relation in GRBs as in other accreting compact objects switching from accretion to propeller. This relation provides also an estimate of the radiative efficiency of GRBs, which we find to be several times lower than radiatively efficient accretion in X-ray binaries and in agreement with theoretical expectations. These results provide additional support to the idea that at least some GRBs are powered by magnetars surrounded by an accretion disc.

We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed $c_{\rm s}$, Alfv\'en speed $v_{\rm A}$, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity onto the fluid-frame wavevector is larger than the projection of the Alfv\'en speed onto the same direction, i.e., shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed $v$, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{re} \equiv v/v_{\rm f\perp} \sqrt{(1-v_{\rm f\perp}^2)/(1-v^2)} \cos\theta=\sqrt{2}$, where $v_{rm f\perp}=[v_A^2+c_{\rm s}^2(1-v_A^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and $\theta$ is the laboratory-frame angle between the flow velocity and the wavevector projection onto the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes -- both for single objects and for merging binaries -- where the Alfv\'en speed may approach the speed of light.

Context. Giant planets play a major role in multiple planetary systems. Knowing their demographics is important to test their overall impact on planetary systems formation. It is also important to test their formation processes. Recently, three radial velocity surveys have established radial distributions of giant planets. All show a steep increase up to 1-3 au, and two suggest a decrease beyond. Aims. We aim at understanding the limitations associated with the characterization of long-period giant radial velocity planets, and to estimate their impact on the radial distribution of these planets. Methods. We revisit the results obtained by two major surveys that derived such radial distributions, using the RV data available at the time of the surveys as well as, whenever possible, new data. Results. We show that the radial distributions published beyond (5-8 au) are not secure. More precisely, the decrease of the radial distribution beyond the peak at 1-3 au is not confirmed.

The astrophysical environments capable of triggering heavy-element synthesis via rapid neutron capture (the r-process) remain uncertain. While binary neutron star mergers (NSMs) are known to forge r-process elements, certain rare supernovae (SNe) have been theorized to supplement, or even dominate, r-production by NSMs. However, the most direct evidence for such SNe, unusual reddening of the emission caused by the high opacities of r-process elements, has not been observed. Recent work identified the distribution of r-process material within the SN ejecta as a key predictor of the ease with which signals associated with r-process enrichment could be discerned. Though this distribution results from hydrodynamic processes at play during the SN explosion, thus far it has been treated only in a parameterized way. We use hydrodynamic simulations to model how disk winds, the alleged locus of r-production in rare SNe, mix with initially r-process-free ejecta. We study mixing as a function of the wind mass and duration and of the initial SN explosion energy, and find that it increases with the first two of these and decreases with the third. This suggests that SNe accompanying the longest long-duration gamma-ray bursts are promising places to search for signs of r-process enrichment. We use semianalytic radiation transport to connect hydrodynamics to electromagnetic observables, allowing us to assess the mixing level at which the presence of r-process material can be diagnosed from SN light curves. Analytic arguments constructed atop this foundation imply that a wind-driven r-process-enriched SN model is unlikely to explain standard energetic SNe.

Recent gravitational-wave transient catalogs have used $p_\mathrm{astro}$, the probability that a gravitational-wave candidate is astrophysical, to select interesting candidates for further analysis. Unlike false alarm rates, which exclusively capture the statistics of the instrumental noise triggers, $p_\mathrm{astro}$ incorporates the rate at which triggers are generated by both astrophysical signals and instrumental noise in estimating the probability that a candidate is astrophysical. Multiple search pipelines can independently calculate $p_\mathrm{astro}$, each employing a specific data reduction. While the range of $p_\mathrm{astro}$ results can help indicate the range of uncertainties in its calculation, it complicates interpretation and subsequent analyses. We develop a statistical formalism to calculate a $\textit{unified } p_\mathrm{astro}$ for gravitational-wave candidates, consistently accounting for triggers from all pipelines, thereby incorporating extra information about a signal that is not available with any one single pipeline. We demonstrate the properties of this method using a toy model and by application to the publicly available list of gravitational-wave candidates from the first half of the third LIGO--Virgo--KAGRA observing run. Adopting a unified $p_\mathrm{astro}$ for future catalogs would provide a simple and easy-to-interpret selection criterion that incorporates a more complete understanding of the strengths of the different search pipelines

SkyPortal is an open-source platform designed to efficiently discover interesting transients, manage follow-up, perform characterization, and visualize the results, all in one application. By enabling fast access to archival and catalog data, cross-matching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, SkyPortal has been operating at scale for > 2 yr for the Zwicky Transient Facility Phase II community, with hundreds of users, containing tens of millions of time-domain sources, interacting with dozens of telescopes, and enabling community reporting. While SkyPortal emphasizes rich user experiences (UX) across common frontend workflows, recognizing that scientific inquiry is increasingly performed programmatically, SkyPortal also surfaces an extensive and well-documented API system. From backend and frontend software to data science analysis tools and visualization frameworks, the SkyPortal design emphasizes the re-use and leveraging of best-in-class approaches, with a strong extensibility ethos. For instance, SkyPortal now leverages ChatGPT large-language models (LLMs) to automatically generate and surface source-level human-readable summaries. With the imminent re-start of the next-generation of gravitational wave detectors, SkyPortal now also includes dedicated multi-messenger features addressing the requirements of rapid multi-messenger follow-up: multi-telescope management, team/group organizing interfaces, and cross-matching of multi-messenger data streams with time-domain optical surveys, with interfaces sufficiently intuitive for the newcomers to the field. (abridged)

We present a study of a directional search for Dark Matter boosted forward when scattered by cosmic-ray nuclei, using a module of the NEWSdm experiment. The boosted Dark Matter flux at the edge of the Earth's atmosphere is expected to be pointing to the Galactic Center, with a flux 15 to 20 times larger than in the transverse direction. The module of the NEWSdm experiment consists of a 10 kg stack of Nano Imaging Trackers, i.e.~newly developed nuclear emulsions with AgBr crystal sizes down to a few tens of nanometers. The module is installed on an equatorial telescope. The relatively long recoil tracks induced by boosted Dark Matter, combined with the nanometric granularity of the emulsion, result in an extremely low background. This makes an installation at the INFN Gran Sasso laboratory, both on the surface and underground, viable. A comparison between the two locations is made. The angular distribution of nuclear recoils induced by boosted Dark Matter in the emulsion films at the surface laboratory is expected to show an excess with a factor of 3.5 in the direction of the Galactic Center. This excess allows for a Dark Matter search with directional sensitivity. The surface laboratory configuration prevents the deterioration of the signal in the rock overburden and it emerges as the most powerful approach for a directional observation of boosted Dark Matter with high sensitivity. We show that, with this approach, a 10 kg module of the NEWSdm experiment exposed for one year at the Gran Sasso surface laboratory can probe Dark Matter masses between 1 keV/c$^2$ and 1 GeV/c$^2$ and cross-section values down to $10^{-30}$~cm$^2$ with a directional sensitive search.

We report the discovery of 2M0056-08 as an equal-mass eclipsing binary (EB), comprising two red straggler stars (RSSs) with an orbital period of 33.9 d. Both stars have masses of 1.419 Msun, identical to within 0.2%. Both stars appear to be in the early red-giant phase of evolution; however, they are far displaced to cooler temperatures and lower luminosities compared to standard stellar models. The broadband spectral energy distribution shows NUV excess and X-ray emission, consistent with chromospheric and coronal emission from magnetically active stars; indeed, the stars rotate more rapidly than typical red giants and they evince light curve modulations due to spots. These modulations also reveal the stars to be rotating synchronously with one another. There is evidence for excess FUV emission and long-term modulations in radial-velocities; it is not clear whether these are also attributable to magnetic activity or if they reveal a tertiary companion. Stellar evolution models modified to account for the effects of spots can reproduce the observed radii and temperatures of the RSSs. If the system possesses a white dwarf tertiary, then mass-transfer scenarios could explain the manner by which the stars came to possess such remarkably identical masses and by which they came to be sychronized. However, if the stars are presumed to have been formed as identical twins, and they managed to become tidally synchronized as they evolved toward the red giant branch, then all of the features of the system can be explained via activity effects, without requiring a complex dynamical history.

The apparent shape of galaxy clustering depends on the adopted cosmology used to convert observed redshift to comoving distance, the $r(z)$ relation, as it changes the line elements along and across the line of sight differently. The Alcock-Paczy\'nski (AP) test exploits this property to constrain the expansion history of the universe. We present an extensive review of past studies on the AP test. We adopt an extended AP test method introduced by Park et al. (2019), which uses the full shape of redshift-space two-point correlation function (CF) as the standard shape, and apply it to the SDSS DR7, BOSS, and eBOSS LRG samples covering the redshift range up to $z=0.8$.We calibrate the test against the nonlinear cosmology-dependent systematic evolution of the CF shape using the Multiverse simulations. We focus on examining whether or not the flat $\Lambda$CDM `concordance' model is consistent with observation. We constrain the flat $w$CDM model to have $w=-0.892_{-0.050}^{+0.045}$ and $\Omega_m=0.282_{-0.023}^{+0.024}$ from our AP test alone, which is significantly tighter than the constraints from the BAO or SNe I$a$ methods by a factor of 3 - 6. When the AP test result is combined with the recent BAO and SNe I$a$ results, we obtain $w=-0.903_{-0.023}^{+0.023}$ and $\Omega_m=0.285_{-0.009}^{+0.014}$. This puts a strong tension with the flat $\Lambda$CDM model with $w=-1$ at $4.2\sigma$ level. Consistency with $w=-1$ is obtained only when the Planck CMB observation is combined. It remains to see if this tension between observations of galaxy distribution at low redshifts and CMB anisotropy at the decoupling epoch becomes greater in the future studies and leads us to a new paradigm of cosmology.

Recently the Galactic and Extra-galactic All-sky Murchison Widefield Array survey has published 27 new candidate radio supernova remnants (SNRs) which are located within the longitude ranges of 345{\deg} < l < 60{\deg} and 180{\deg} < l < 240{\deg}. To search for the gamma-ray counterparts of these candidate radio SNRs, we analyzed 14 years of {\it Fermi}-LAT data in the energy range of 1 - 300 GeV. There are three promising SNRs; G18.9$-$1.2, G23.1$+$0.1, and G28.3$+$0.2, which we detected at a significance level of $\sim$9$\sigma$, $\sim$13$\sigma$, and $\sim$12$\sigma$, respectively. Here we report the results of our morphological and spectral analyses of G18.9$-$1.2, G23.1$+$0.1, and G28.3$+$0.2. No extended gamma-ray emission is detected for any of these SNRs. Our analysis of the 3 SNRs' {\it Fermi}-LAT gamma-ray emission showed that their best-fit positions (if assumed point-like) overlap with the locations of the corresponding GLEAM counterparts.

The Sun constantly releases radiation and plasma into the heliosphere. Sporadically, the Sun launches solar eruptions such as flares and coronal mass ejections (CMEs). CMEs carry away a huge amount of mass and magnetic flux with them. An Earth-directed CME can cause serious consequences to the human system. It can destroy power grids/pipelines, satellites, and communications. Therefore, accurately monitoring and predicting CMEs is important to minimize damages to the human system. In this study we propose an ensemble learning approach, named CMETNet, for predicting the arrival time of CMEs from the Sun to the Earth. We collect and integrate eruptive events from two solar cycles, #23 and #24, from 1996 to 2021 with a total of 363 geoeffective CMEs. The data used for making predictions include CME features, solar wind parameters and CME images obtained from the SOHO/LASCO C2 coronagraph. Our ensemble learning framework comprises regression algorithms for numerical data analysis and a convolutional neural network for image processing. Experimental results show that CMETNet performs better than existing machine learning methods reported in the literature, with a Pearson product-moment correlation coefficient of 0.83 and a mean absolute error of 9.75 hours.

In Luparello et al. 2023, a new and hitherto unknown CMB foreground was detected. A systematic decrease in Cosmic Microwave Background (CMB) temperatures around nearby large spiral galaxies points to an unknown interaction with CMB photons in a sphere up to several projected Mpc around these galaxies. We investigate to which extent this foreground may impact the CMB fluctuations map and create the so-called CMB anomalies. Using the observed temperature decrements around the galaxies, and making some general assumptions about the unknown interaction, we propose a common radial temperature profile. By assigning this profile to nearby galaxies in the redshift range $z=[0.004,0.02]$ we create a foreground map model. We find a remarkable resemblance between this temperature model map based on nearby galaxies and the Planck CMB map. Out of 1000 simulated maps, none of them show such a strong correlation with the foreground map over both large and small angular scales. In particular, the quadrupole, octopole, as well as $\ell=4$ and $\ell=5$ modes correlate with the foreground map to high significance. Furthermore, one of the most prominent temperature decrements in the foreground map coincides with the position of the CMB cold spot. The largest scales of the CMB and thereby the cosmological parameters, may have important changes after proper corrections of this foreground component. However, reliable CMB corrected maps can only be derived when suitable physical mechanisms are proposed and tested.

Blue supergiants (BSGs) are important objects to study the intermediate phases of massive star evolution, helping to constrain evolutionary models. However, the lack of a holistic study of a statistically significant and unbiased sample of these objects makes several long-standing questions about their nature to remain unsolved. The present and other upcoming papers of the IACOB series are focused in studying - from a pure empirical point of view - a sample of 500 Galactic O9 - B9 stars with luminosity classes I and II (plus 250 late O- and early B-type stars with luminosity classes III, IV and V) and covering distances up to 4 kpc from the Sun. We compile an initial set of 11000 high-resolution spectra of 1600 Galactic late O- and B-type stars. We use a new novel spectroscopic strategy based on a simple fitting of the Hbeta line to select stars in a specific region of the spectroscopic HR diagram. We evaluate the completeness of our sample using the Alma Luminous Star catalog (ALS III) and Gaia-DR3 data. We show the benefits of the proposed strategy for identifying BSGs descending from stellar objects born as O-type stars, in the context of single star evolution. The resulting sample reaches a high level of completeness with respect to the ALS III catalog, gathering the 80% for all-sky targets brighter than Bmag < 9 located within 2 kpc. However, we identify the need for new observations in specific regions of the Southern hemisphere. In conclusion, we have explored a very fast and robust method to select BSGs, providing a valuable tool for large spectroscopic surveys like WEAVE-SCIP or 4MIDABLE-LR, and highlighting the risk of using spectral classifications from the literature. Upcoming works will make use of this large and homogeneous spectroscopic sample to study specific properties of these stars in detail. We initially provide first results about their rotational properties.

Solar flares and coronal mass ejections (CMEs) are thought to be the most powerful events on the Sun. They can release energy as high as 10^32 erg in tens of minutes,and could produce solar energetic particles (SEPs) in the interplanetary space. We explore global energy budgets of solar major eruptions on 6 September 2017, including the energy partition of a powerful solar flare, the energy budget of the accompanied CME and SEPs. In the wavelength range shortward of 222 nm, a major contribution of the flare radiated energy is in the soft X-ray (SXR) 0.1-7 nm domain. The flare energy radiated at wavelengths of Ly-alpha and middle ultraviolet is larger than that radiated in the extreme ultraviolet wavelength, but it is much less than that radiated in the SXR waveband. The total flare radiated energy could be comparable to the thermal and nonthermal energies. The energies carried by the major flare and its accompanied CME are roughly equal, and they are both powered by the magnetic free energy in the AR NOAA 12673. Moreover, the CME is efficient in accelerating SEPs, and that the prompt component (whether it comes from the solar flare or the CME) contributes only a negligible fraction.

We present an in-depth analysis of a newly proposed correlation function in visibility space, between the $E$ and $B$ modes of the linear polarization, hereafter the $EB$-correlation, for a set of time-averaged GRMHD simulations compared with the phase map from different semi-analytic models as well as the Event Horizon Telescope (EHT) 2017 data for M87* source. We demonstrate that the phase map of the time-averaged $EB$-correlation contains novel information that might be linked to the BH spin, accretion state and the electron temperature. A detailed comparison with a semi-analytic approach with different azimuthal expansion modes shows that to recover the morphology of the real/imaginary part of the correlation function and its phase, we require higher orders of these azimuthal modes. To extract the phase features, we propose to use the Zernike polynomial reconstruction developing an empirical metric to break degeneracies between models with different BH spins that are qualitatively similar. We use a set of different geometrical ring models with various magnetic and velocity field morphologies and show that both the image space and visibility based $EB$-correlation morphologies in MAD simulations can be explained with simple fluid and magnetic field geometries as used in ring models. SANEs by contrast are harder to model, demonstrating that the simple fluid and magnetic field geometries of ring models are not sufficient to describe them owing to higher Faraday Rotation depths. A qualitative comparison with the EHT data demonstrates that some of the features in the phase of $EB$-correlation might be well explained by the current models for BH spins as well as electron temperatures, while others may require a larger theoretical surveys.

The China Space Station Telescope (CSST) is a forthcoming Stage IV galaxy survey. It will simultaneously undertake the photometric redshift (photo-z) and slitless spectroscopic redshift (spec-z) surveys mainly for weak lensing and galaxy clustering studies. The two surveys cover the same sky area and overlap on the redshift range. Due to the sparse number density of the spec-z sample at $z>1$, it limits the constraints on the scale of Baryon Acoustic Oscillations (BAO). By cross-correlating the spec-z sample with the high-density photo-z sample, we can effectively enhance the constraints on the angular diameter distances from BAO. We estimate a greater than 35 per cent improvement utilising the Fisher matrix formalism. Such improvement is robust against different systematic effects including the systematic noise and the redshift success rate of the spec-z survey, as well as the photo-z error. Our study can be a reference for future BAO analysis on real CSST data. The methodology can be applied to other surveys with spec-z and photo-z data in the same survey volume.

Asymmetric features in exoplanet transit light curves are often interpreted as a gravity darkening effect especially if there is spectroscopic evidence of a spin-orbit misalignment. Since other processes can also lead to light curve asymmetries this may lead to inaccurate gravity darkening parameters. Here we investigate the case of non-radial pulsations as possible sources of asymmetry and likely source of misinterpreted parameters through simulations. We obtained a series of simulated transit light curves of a hypothetical exoplanet-star system: a host star with no gravity darkening exhibiting small amplitude pulsations, and a typical hot Jupiter in a circular, edge-on orbit. A number of scenarios of pulsations of various amplitudes were considered, and a proper account of the obscuring effect of transits on all the surface intensity components was made. The magnitude of amplitude and phase modulations of nonradial pulsations during transits was also also investigated. We then fitted both a non-gravity-darkened, and a gravity-darkened, free spin-orbit axis model on the data. The Akaike and Bayesian Information Criteria were used for an objective selection of the most plausible model. We then explored the dependence of the parameter deviations on the pulsation properties, in order to identify configurations that can lead to falsely misaligned solutions. Low-amplitude pulsations in general do not affect the determination of the system parameters beyond their noise nature. However, frequencies close to multiples of the orbital frequency are found to cause distortions leading to solutions with a side tilted stellar rotational axis, they are therefore preferable to clean beforehand for the sake of a correct analysis. Additionally, for cases with higher-amplitude pulsations, it is recommended to pre-process and clean the pulsations before analysis.

In 2022 we carried out an upgrade of the speckle polarimeter (SPP) -- the facility instrument of the 2.5-m telescope of the Caucasian Observatory of the SAI MSU. During the overhaul, CMOS Hamamatsu ORCA-Quest qCMOS C15550-20UP was installed as the main detector, some drawback of the previous version of the instrument were eliminated. In this paper, we present a description of the instrument, as well as study some features of the CMOS detector and ways to take them into account in speckle interferometric processing. Quantitative comparison of CMOS and EMCCD in the context of speckle interferometry is performed using numerical simulation of the detection process. Speckle interferometric observations of 25 young variable stars are given as an example of astronomical result. It was found that BM And is a binary system with a separation of 273 mas. The variability of the system is dominated by the brightness variations of the main component. A binary system was also found in NSV 16694 (TYC 120-876-1). The separation of this system is 202 mas.

It is a common belief that for magnetic fields typical for accreting neutron stars in High-Mass X-ray Binaries vacuum polarization only affects the propagation of polarized emission in the neutron star magnetosphere. We show that vacuum resonances can significantly alter the emission from the poles of accreting neutron stars. The effect is similar to vacuum polarization in the atmospheres of isolated neutron stars and can result in suppression of the continuum and the cyclotron lines. It is enhanced by magnetic Comptonization in the hot plasma and proximity to the electron cyclotron resonance. We present several models to illustrate the vacuum polarization effect for various optically thick media and discuss how the choice of polarization modes affects the properties of the emergent radiation by simulating polarized energy- and angle-dependent radiative transfer. Polarization effects, including vacuum polarization, crucially alter the emission properties. Together with strongly angle- and energy- dependent magnetic Comptonization, they result in a complex spectral shape, which can be described by dips and humps on top of a power-law-like continuum with high-energy cutoff. These effects provide a possible explanation for the common necessity of additional broad Gaussian components and two-component Comptonization models that are used to describe spectra of accreting X-ray pulsars. We also demonstrate the character of depolarization introduced by the radiation field's propagation inside the inhomogeneous emission region.

In the interstellar medium (ISM), the complex organic molecules that contain the thiol group ($-$SH) play an important role in the polymerization of amino acids. We look for SH-bearing molecules in the chemically rich solar-type protostar IRAS 16293-2422. After the extensive spectral analysis using the local thermodynamic equilibrium (LTE) model, we have detected the rotational emission lines of trans-isomer monothioformic acid (t-HC(O)SH) towards the IRAS 16293 B using the Atacama Large Millimeter/Submillimeter Array (ALMA). We did not observe any evidence of cis-isomer monothioformic acid (c-HC(O)SH) towards the IRAS 16293 B. The column density of t-HC(O)SH towards the IRAS 16293 B was (1.02$\pm$0.6)$\times$10$^{15}$ cm$^{-2}$ with an excitation temperature of 125$\pm$15 K. The fractional abundance of t-HC(O)SH with respect to H$_{2}$ towards the IRAS 16293 B is 8.50$\times$10$^{-11}$. The column density ratio of t-HC(O)SH/CH$_{3}$SH towards the IRAS 16293 B is 0.185. We compare our estimated abundance of t-HC(O)SH towards the IRAS 16293 B with the abundance of t-HC(O)SH towards the galactic center quiescent cloud G+0.693-0.027 and hot molecular core G31.41+0.31. After the comparison, we found that the abundance of t-HC(O)SH towards the IRAS 16293 B is several times of magnitude lower than G+0.693-0.027 and G31.41+0.31. We also discuss the possible formation mechanism of t-HC(O)SH in the ISM.

Cool stars, such as M giants, can only be analysed in the near-infrared (NIR) regime due to the ubiquitous TiO features in optical spectra of stars with Teff < 4000 K. In dust obscured regions, like the inner bulge and Galactic Center, the intrinsically bright M giants observed in the NIR is an optimal option to determine their stellar abundances. Due to uncertainties in photometric methods, a method to determine the stellar parameters for M giants from the NIR spectra themselves is needed. We have carried out new observations of 44 M giant stars (also in APOGEE DR17) with IGRINS (R=45,000) mounted on the Gemini South telescope. We also obtained HK band IGRINS spectra of six nearby well-studied M giants from the IGRINS spectral library. Using this sample, we have developed a method to determine the stellar parameters for M giants from the NIR spectra by spectral synthesis using SME. The method is validated using the six nearby well-studied M-giants. We demonstrate the accuracy and precision by determining stellar parameters and $\alpha$-element trends versus metallicity for solar neighbourhood M giants. The effective temperatures that we derive (tested for 3400$\lesssim$ Teff $\lesssim$4000\,K) agree excellently with the six nearby M giants which indicates that the accuracy is indeed high. For the 43 solar neighborhood M giants, our Teff, logg, [Fe/H], $\xi_\mathrm{micro}$, [C/Fe], [N/Fe], and [O/Fe] are in unison with APOGEE with mean differences and scatter (our method - APOGEE) of -67$\pm$33 K, -0.31$\pm$0.15 dex, 0.02$\pm$0.05 dex, 0.22$\pm$0.13 km/s, -0.05$\pm$0.06 dex, 0.06$\pm$0.06 dex, and 0.02$\pm$0.09 dex, respectively. The $\alpha$-element trends versus metallicity for Mg, Si, Ca and Ti are consistent with both APOGEE DR17 trends for the same stars as well as with the GILD optical trends. We also find clear enhancement in abundances for thick disc stars.

We use the SIMBA galaxy formation simulation suite to explore anisotropies in the properties of circumgalactic gas that result from accretion and feedback processes. We particularly focus on the impact of bipolar active galactic nuclei (AGN) jet feedback as implemented in SIMBA, which quenches galaxies and has a dramatic effect on large-scale gas properties. We show that jet feedback at low redshifts is most common in the stellar mass range $(1-5)\times 10^{10}M_\odot$, so we focus on galaxies with active jets in this mass range. In comparison to runs without jet feedback, jets cause lower densities and higher temperatures along the galaxy minor axis (SIMBA jet direction) at radii >=$0.5r_{200c}-4r_{200c}$ and beyond. This effect is less apparent at higher or lower stellar masses, and is strongest within green valley galaxies. The metallicity also shows strong anisotropy out to large scales, driven by star formation feedback. We find substantially stronger anisotropy at <=$0.5r_{200c}$, but this also exists in runs with no explicit feedback, suggesting that it is due to anisotropic accretion. Finally, we explore anisotropy in the bulk radial motion of the gas, finding that both star formation and AGN wind feedback contribute to pushing the gas outwards along the minor axis at <=1 Mpc, but AGN jet feedback further causes bulk outflow along the minor axis out to several Mpc, which drives quenching via gas starvation. These results provide observational signatures for the operation of AGN feedback in galaxy quenching.

In this work we revisit power law, $\frac{1}{M^2}R^\beta$, inflation to find the deviations from $R^2$ inflation allowed by current CMB and LSS observations. We compute the power spectra for scalar and tensor perturbations numerically and perform MCMC analysis to put constraints on parameters $M$ and $\beta$ from Planck-2018, BICEP3 and other LSS observations. We consider general reheating scenario and also vary the number of e-foldings during inflation, $N_{pivot}$, along with the other parameters. We find $\beta = 1.966^{+0.035}_{-0.042}$, $M= \left(3.31^{+5}_{-2}\right)\times 10^{-5}$ and $N_{pivot} = 41^{+10}_{-10}$ with $95\%\, C.\, L.$. This indicates that the current observations allow deviation from Starobinsky inflation. The scalar spectral index, $n_s$, and tensor-to-scalar ratio, $r$, derived from these parameters, are consistent with the Planck and BICEP3 observations.

FU Orionis (FUor) and EX Lupi (EXor) type objects are two groups of peculiar and rare pre-main sequence low-mass stars that are undergoing powerful accretion outbursts during their early stellar evolution. Water masers are widespread in star forming regions and are powerful probes of mass accretion and ejection, but little is known about the prevalence of them toward FUors/EXors. We perform the first systematic search for the 22.2 GHz water maser line in FUors/EXors to determine its overall incidence to perform follow-up high angular resolution observations. We used the Effelsberg 100-m radio telescope to observe the 22.2 GHz H2O maser toward a sample of 51 objects. We detect 5 water masers; 3 are associated with eruptive stars, resulting in a 6% detection rate for eruptive sources. These detections include one EXor, V512 Per (also known as SVS 13 or SVS 13A), and two FUors, Z CMa and HH 354 IRS. This is the first reported detection of water maser emission towards HH 354 IRS. We detect water maser emission in our pointing towards the FUor binary RNO 1B/1C, which most likely originates from the nearby deeply embedded source IRAS 00338+6312 (~4'', from RNO 1B/1C). Emission was also detected from H$_2$O(B) (also known as SVS 13C), a Class 0 source ~30'', from the EXor V512 Per. The peak flux density of H$_2$O(B) in our observations, 498.7 Jy, is the highest observed to date. In addition to the two non-eruptive Class 0 sources (IRAS 00338+6312 and H$_2$O(B) /SVS 13C), we detect maser emission towards one Class 0/I (HH 354 IRS) and two Class I (V512 Per and Z CMa) eruptive stars. We demonstrate the presence of 22.2 GHz water maser emission in FUor/EXor systems, opening the way to radio interferometric observations to study these eruptive stars on small scales. Comparing our data with historical observations suggest that multiple water maser flares have occurred in both V512 Per and H$_2$O(B).

Observations across the electromagnetic spectrum of radiative processes involving interstellar dust -- emission, extinction, and scattering -- are used to constrain the parameters of dust models and more directly to aid in foreground removal of dust for extragalactic and cosmology observations. The more complementary observations, the better. Here, we quantify the relationship between scattered light and thermal emission from dust in a diffuse (cirrus) intermediate latitude cloud, Spider, using data from the Dragonfly Telephoto Array and the Herschel Space Observatory. A challenge for optical observations of faint cirrus is accurate removal of a contaminating spatially varying sky background. We present a technique to analyse two images of the same cirrus field concurrently, correlating pixel values to capture the relationship and simultaneously fitting the sky background as a complex non-correlating additive component. For the Spider, we measure a $g-r$ color of 0.644$\pm 0.024$ and a visible wavelength to 250 $\mu$m intensity ratio of $10^{-3} \times (0.855 \pm0.025)$ and $10^{-3} \times (1.55\pm0.08)$ for $g$ and $r$-band respectively. We show how to use any dust model that matches the thermal dust emission to predict an upper limit to the amount of scattered light. The actual brightness of the cirrus will be fainter than this limit because of anisotropic scattering by the dust combined with anisotropy of the incident interstellar radiation field (ISRF). Using models of dust and the ISRF in the literature we illustrate that the predicted brightness is indeed lower, though not as faint as the observations indicate.

Current measurements of planet population as a function of stellar mass show three seemingly contradictory signatures: close-in super-Earths are more prevalent around M dwarfs than FGK dwarfs; inner super-Earths are correlated with outer giants; and outer giants are less common around M dwarfs than FGK dwarfs. Here, we build a simple framework that combines the theory of pebble accretion with the measurements of dust masses in protoplanetary disks to reconcile all three observations. First, we show that cooler stars are more efficient at converting pebbles into planetary cores at short orbital periods. Second, when disks are massive enough to nucleate a heavy core at 5 AU, more than enough dust can drift in to assemble inner planets, establishing the correlation between inner planets and outer giants. Finally, while stars of varying masses are similarly capable of converting pebbles into cores at long orbital periods, hotter stars are much more likely to harbor more massive dust disks so that the giant planet occurrence rate rises around hotter stars. Our results are valid over a wide range of parameter space for a disk accretion rate that follows $\dot{M}_\star \sim 10^{-8}\,M_\odot\,{\rm yr}^{-1}(M_\star/M_\odot)^2$. We predict a decline in mini-Neptune population (but not necessarily terrestrial planets) around stars lighter than $\sim 0.3-0.5 \, M_\odot$. Cold giants ($\gtrsim$5 AU), if they exist, should remain correlated with inner planets even around lower mass stars.

In exoplanetary systems, the interaction between the central star and the planet can trigger Auroral Radio Emission (ARE), due to the Electron Cyclotron Maser mechanism. The high brightness temperature of this emission makes it visible at large distances, opening new opportunities to study exoplanets and to search for favourable conditions for the development of extra-terrestrial life, as magnetic fields act as a shield that protects life against external particles and influences the evolution of the planetary atmospheres. In the last few years, we started an observational campaign to observe a sample of nearby M-type stars known to host exoplanets with the aim to detect ARE. We observed YZ Ceti with the upgraded Giant Metrewave Radio Telescope (uGMRT) in band 4 (550-900 MHz) nine times over a period of five months. We detected radio emission four times, two of which with high degree of circular polarization. With statistical considerations we exclude the possibility of flares due to stellar magnetic activity. Instead, when folding the detections to the orbital phase of the closest planet YZ Cet b, they are at positions where we would expect ARE due to star-planet interaction (SPI) in sub-Alfvenic regime. With a degree of confidence higher than 4.37 sigma, YZ Cet is the first extrasolar systems with confirmed SPI at radio wavelengths. Modelling the ARE, we estimate a magnetic field for the star of about 2.4 kG and we find that the planet must have a magnetosphere. The lower limit for the polar magnetic field of the planet is 0.4 G.

The spin-orbit angle, or obliquity, is a powerful observational marker that allows us to access the dynamical history of exoplanetary systems. Here, we have examined the distribution of spin-orbit angles for close-in exoplanets and put it in a statistical context of tidal interactions between planets and their stars. We confirm the observed trends between the obliquity and physical quantities directly connected to tides, namely the stellar effective temperature, the planet-to-star mass ratio, and the scaled orbital distance. We further devised a tidal efficiency factor combining critical parameters that control the strength of tidal effects and used it to corroborate the strong link between the spin-orbit angle distribution and tidal interactions. In particular, we developed a readily usable formula to estimate the probability that a system is misaligned, which will prove useful in global population studies. By building a robust statistical framework, we reconstructed the distribution of the three-dimensional spin-orbit angles, allowing for a sample of nearly 200 true obliquities to be analyzed for the first time. This realistic distribution maintains the sky-projected trends, and additionally hints toward a striking pileup of truly aligned systems. The comparison between the full population and a pristine subsample unaffected by tidal interactions suggests that perpendicular architectures are resilient toward tidal realignment, providing evidence that orbital misalignments are sculpted by disruptive dynamical processes that preferentially lead to polar orbits. On the other hand, star-planet interactions seem to efficiently realign or quench the formation of any tilted configuration other than for polar orbits, and in particular for antialigned orbits.

The BlueWalker 3 (BW3) satellite was folded into a compact object when launched on 2022 September 11. The spacecraft's apparent visual magnitude initially ranged from about 4 to 8. Observations on November 11 revealed that the brightness increased by 4 magnitudes which indicated that the spacecraft had deployed into a large flat-panel shape. The satellite then faded by several magnitudes in December before returning to its full luminosity; this was followed by additional faint periods in 2023 February and March. We discuss the probable cause of the dimming phenomena and identify a geometrical circumstance where the satellite is abnormally bright. The luminosity of BW3 can be represented with a brightness model which is based on the satellite shape and orientation as well as a reflection function having Lambertian and pseudo-specular components. Apparent magnitudes are most frequently between 2.0 and 3.0. When BW3 is near zenith the magnitude is about 1.4.

Gamma-ray bursts (GRBs) 201015A and 201216C are valuable cases with detection of very high energy (VHE) gamma-ray afterglows. By analysing their prompt emission data, we find that GRB 201216C is an extremely energetic long GRB with a hard gamma-ray spectrum, while GRB 201015A is a relative sub-energetic, soft spectrum GRB. Attributing their radio-optical-X-ray afterglows to the synchrotron radiation of the relativistic electrons accelerated in their jets, we fit their afterglow lightcurves with the standard external shock model and infer their VHE afterglows from the synchrotron self-Compton scattering process of the electrons. It is found that the jet of GRB 201015A is mid-relativistic ($\Gamma_0=44$) surrounded by a very dense medium ($n=1202$ cm$^{-3}$) and the jet of GRB 201216C is ultra-relativistic ($\Gamma_0=331$) surrounded by a moderate dense medium ($n=5$ cm$^{-3}$). The inferred peak luminosity of the VHE gamma-ray afterglows of GRB 201216C is approximately $10^{-9}$ erg cm$^{-2}$ s$^{-1}$ at $57-600$ seconds after the GRB trigger, making it can be detectable with the MAGIC telescopes at a high confidence level, even the GRB is at a redshift of 1.1. Comparing their intrinsic VHE gamma-ray lightcurves and spectral energy distributions with GRBs~180720B, 190114C, and 190829A, we show that their intrinsic peak luminosity of VHE gamma-ray afterglows at $10^{4}$ seconds post the GRB trigger is variable from $10^{45}$ to $5\times 10^{48}$ erg s$^{-1}$, and their kinetic energy, initial Lorentz factor, and medium density are diverse among bursts.

I discuss Local Group galaxies from the perspective of external galaxies that define benchmark scaling relations. Making use of this information leads to a model for the Milky Way that includes bumps and wiggles due to spiral arms. This model reconciles the terminal velocities observed in the interstellar medium with the rotation curve derived from stars, correctly predicts the gradual decline of the outer rotation curve ($dV/dR = -1.7\;\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}^{-1}$), and extrapolates well out to 50 kpc. Rotationally supported Local Group galaxies are in excellent agreement with the baryonic Tully-Fisher relation. Pressure supported dwarfs that are the most likely to be in dynamical equilibrium also align with this relation. Local Group galaxies thus appear to be normal members of the low redshift galaxy population. There is, however, a serious tension between the dynamical masses of the Milky Way and M31 ($M_{200} \approx 1.4$ and $1.6 \times 10^{12}\;\mathrm{M}_{\odot}$, respectively) and those expected from the stellar mass-halo mass relation of abundance matching ($M_{200} \approx 3$ and $20 \times 10^{12}\;\mathrm{M}_{\odot}$, respectively).

The abundances of CHNOS are crucial for the composition of planets. At the onset of planet formation, large amounts of these elements are stored in ices on dust grains in planet-forming disks. The evolution of this ice is affected by dynamical transport, collisional processes, and the formation and sublimation of ice. We aim to constrain the disk regions where these processes are fully coupled, and develop a flexible modelling approach that is able to predict the effects of these processes acting simultaneously on the CHNOS budgets of the dust in these regions. We compared timescales associated with these disk processes to constrain the disk regions where this approach is necessary, and developed the SHAMPOO code, which tracks the CHNOS abundances in the ice mantle of a single monomer dust particle, embedded in a larger aggregate and undergoing these processes simultaneously. The adsorption and photodesorption of monomer ices depend on the depth of the monomer in the aggregate. We investigated the effect of fragmentation velocity and aggregate filling factor on the amount of ice on monomers residing at r = 10 AU. The locations where disk processes are fully coupled depend on both grain size and ice species. Monomers embedded in aggregates with fragmentation velocities of 1 m/s are able to undergo adsorption and photodesorption more often compared to a fragmentation velocity of 5 m/s or 10 m/s. Aggregates with a filling factor of $10^{-3}$ are able to accumulate ice 22 times faster on average than aggregates with a filling factor of 1. As different grain sizes are coupled through collisions and the grain ice consists of multiple ice species, it is difficult to isolate the locations where disk processes are fully coupled, necessitating the development of the SHAMPOO code. The processing of ice may not be spatially limited to dust aggregate surfaces for either fragile or porous aggregates.

Planets orbiting M-dwarf stars are prime targets in the search for rocky exoplanet atmospheres. The small size of M dwarfs renders their planets exceptional targets for transmission spectroscopy, facilitating atmospheric characterization. However, it remains unknown whether their host stars' highly variable extreme-UV radiation environments allow atmospheres to persist. With JWST, we have begun to determine whether or not the most favorable rocky worlds orbiting M dwarfs have detectable atmospheres. Here, we present a 2.8-5.2 micron JWST NIRSpec/G395H transmission spectrum of the warm (700 K, 40.3x Earth's insolation) super-Earth GJ 486b (1.3 R$_{\oplus}$ and 3.0 M$_{\oplus}$). The measured spectrum from our two transits of GJ 486b deviates from a flat line at 2.2 - 3.3 $\sigma$, based on three independent reductions. Through a combination of forward and retrieval models, we determine that GJ 486b either has a water-rich atmosphere (with the most stringent constraint on the retrieved water abundance of H2O > 10% to 2$\sigma$) or the transmission spectrum is contaminated by water present in cool unocculted starspots. We also find that the measured stellar spectrum is best fit by a stellar model with cool starspots and hot faculae. While both retrieval scenarios provide equal quality fits ($\chi^2_\nu$ = 1.0) to our NIRSpec/G395H observations, shorter wavelength observations can break this degeneracy and reveal if GJ 486b sustains a water-rich atmosphere.

We study the extent to which Milne-Eddington inversions are able to retrieve and characterize the magnetic landscape of the solar poles from observations by the spectropolarimeter onboard Hinode. In particular, we evaluate whether a variable magnetic filling factor is an adequate modeling technique for retrieving the intrinsic magnetic properties from every pixel in the polar field of view. We first generate synthetic spectra emerging from a numerical simulation of a "plage" region at an inclined line of sight of 65$^{\circ}$, and degrade the data to emulate real observations. Then, we invert the synthetic spectra with two Milne-Eddington inversion codes that feature different treatments of the magnetic filling factor, and relate the retrieved magnetic quantities back to their original values in the simulation cube. We find that while the apparent retrieved magnetic properties map well the spatially-degraded simulation, the intrinsic magnetic quantities bear little relation to the magnetic field at the native resolution of the simulation. We discuss the systematic biases caused by line-of-sight foreshortening, spatial degradation, photon noise and modeling assumptions embedded in the inversion algorithm.

We demonstrate a prototype kinetic inductance detector (KID) readout system that uses less than 10 mW per pixel. The CCAT-prime RFSoC based readout is capable of reading four independent detector networks of up to 1000 KIDs each. The power dissipation was measured to be less than 40 W while running multi-tone combs on all four channels simultaneously. The system was also used for the first time to perform sweeps and resonator identification on a prototype 280 GHz array.

High-resolution observations of several debris disks reveal structures such as gaps and spirals, suggestive of gravitational perturbations induced by underlying planets. Most existing studies of planet--debris disk interactions ignore the gravity of the disk, treating it as a reservoir of massless planetesimals. In this paper, we continue our investigation into the long-term interaction between a single eccentric planet and an external, massive debris disk. Building upon our previous work, here we consider not only the axisymmetric component of the disk's gravitational potential, but also the non-axisymmetric torque that the disk exerts on the planet (ignoring for now only the non-axisymmetric component of the disk \textit{self}-gravity). To this goal, we develop and test a semi-analytic `$N$-ring' framework that is based on a generalized (softened) version of the classical Laplace--Lagrange secular theory. Using this tool, we demonstrate that even when the disk is less massive than the planet, not only can a secular resonance be established within the disk that leads to the formation of a wide non-axisymmetric gap (akin to those observed in HD 107146, HD 92945, and HD 206893), but that the very same resonance also damps the planetary eccentricity via a process known as resonant friction. We also develop analytic understanding of these findings, finding good quantitative agreement with the outcomes of the $N$-ring calculations. Our results may be used to infer both the dynamical masses of gapped debris disks and the dynamical history of the planets interior to them, as we exemplify for HD 206893.

Emergent bulk properties of matter governed by the strong nuclear force give rise to physical phenomena across vastly different scales, ranging from the shape of atomic nuclei to the masses and radii of neutron stars. They can be accessed on Earth by measuring the spatial extent of the outer skin made of neutrons that characterises the surface of heavy nuclei. The isotope $^{208}$Pb, owing to its simple structure and neutron excess, has been in this context the target of many dedicated efforts. Here, we determine the neutron skin from measurements of particle distributions and their collective flow in $^{208}$Pb+$^{208}$Pb collisions at ultrarelativistic energy performed at the Large Hadron Collider, which are sensitive to the overall size of the colliding $^{208}$Pb ions. By means of state-of-the-art global analysis tools within the hydrodynamic model of heavy-ion collisions, we infer a neutron skin $\Delta r_{np}=0.217\pm0.058$ fm, consistent with nuclear theory predictions, and competitive in accuracy with a recent determination from parity-violating asymmetries in polarised electron scattering. We establish thus a new experimental method to systematically measure neutron distributions in the ground state of atomic nuclei.

We argue that demanding a consistent cosmological history, including the absence of domain walls and strongly interacting relics at the Peccei-Quinn scale, singles out two concrete realizations of hadronic QCD axions as viable dark matter models. These realizations generally feature flavor-violating axion couplings to Standard Model quarks that are unsuppressed at low energies. As a consequence, experiments looking for flavor-violating hadronic processes involving the axion can be sensitive probes of QCD axion dark matter models. In particular, we show that the NA62 and KOTO experiments could detect the $K\rightarrow\pi + a$ decay for axions consistent with the observed dark matter abundance via the post-inflationary misalignment mechanism.

We investigate the reach of future gravitational wave (GW) detectors in probing inflaton couplings with visible sector particles that can either be bosonic or fermionic in nature. Assuming reheating takes place through perturbative quantum production from vacuum in presence of classical inflaton background field, we find that the spectral energy density of the primordial GW generated during inflation becomes sensitive to inflaton-matter coupling. We conclude, obeying bounds from Big Bang Nucleosysthesis and Cosmic Microwave Background, that, e.g., inflaton-scalar couplings of the order of $\sim\mathcal{O}(10^{-20})$ GeV fall within the sensitivity range of several proposed GW detector facilities. However, this prediction is sensitive to the size of the inflationary scale, nature of the inflaton-matter interaction and shape of the potential during reheating. Having found the time-dependent effective inflaton decay width, we also discuss its implications for dark matter (DM) production from the thermal plasma via UV freeze-in during reheating. It is shown, that one can reproduce the observed DM abundance for its mass up to several PeVs, depending on the dimension of the operator connecting DM with the thermal bath and the associated scale of the UV physics. Thus we promote primordial GW to observables sensitive to feebly coupled inflaton, which is very challenging if not impossible to test in conventional particle physics laboratories or astrophysical measurements.

For a cosmological first-order phase transition in the early Universe, the associated stochastic gravitational wave background is usually dominated by sound waves from plasma fluid motions, which have been analytically modeled as a random superposition of freely propagating sound shells but with the force by the scalar field that produces the self-similar profile removed. In this Letter, we propose a new analytic sound shell model by focusing on the forced propagating contribution from the initial collision stage of sound shells when their self-similar profiles are still maintained by the moving bubble walls. We reproduce the causal $k^3$-scaling in the infrared consistent with numerical simulations, and also recover the broad dome in the power spectrum first observed in numerical simulations. The total sound waves should contain both contributions from forced collisions and free propagation of sound shells at early and late stages of the phase transition, respectively.

In the present paper, we perform a sub-Planckian quantum mode analysis of linear cosmological perturbation in the inflaton field over a classical quasi de-Siter metric background by dynamical horizon exit (DHE) method. In this way, we probe the inflationary regime of a quintessential $\alpha$-attractor model by analysing the COBE/Planck normalized power spectra, spectral indices, tensor to scalar ratio, number of e-folds, running of the spectral index and inflationary Hubble parameter in $k$-space. We compare our results with ordinary $\alpha$-attractor $E$ and $T$ models and with that of Planck-2018 results. Our estimated values of $n_s$ and $r$ lie within $68\%$ CL with respect to Planck data for $k=0.001 - 0.009$ Mpc$^{-1}$ for all values of $\alpha$. The $\alpha$ values, obtained in our calculations satisfy various post inflationary constraints regarding preheating and reheating, reported in current literature. We observe that quintessence sets an upper bound of $\alpha=4.3$ and thereby restricts the model from becoming of the power law type, making it more efficacious than ordinary $\alpha$-attractors in explaining both inflation and dark energy. A striking observation in our analyses is that, unlike in our previous study, we find a continuous values of $\alpha$ within $\frac{1}{10}\leq \alpha\leq 4.3$ for the specified $k$ range. At the end, we have shown that the model parameters constrained in this work give a very small vacuum density $\sim 10^{-117}-10^{-115} M_P^4$ which is an essential criterion for current and future dark energy observations of the universe.

Lithium-drifted silicon [Si(Li)] has been used for decades as an ionizing radiation detector in nuclear, particle, and astrophysical experiments, though such detectors have frequently been limited to small sizes (few cm$^2$) and cryogenic operating temperatures. The 10-cm-diameter Si(Li) detectors developed for the General Antiparticle Spectrometer (GAPS) balloon-borne dark matter experiment are novel particularly for their requirements of low cost, large sensitive area (~10 m$^2$ for the full 1440-detector array), high temperatures (near -40$\,^\circ$C), and energy resolution below 4 keV FWHM for 20--100-keV x-rays. Previous works have discussed the manufacturing, passivation, and small-scale testing of prototype GAPS Si(Li) detectors. Here we show for the first time the results from detailed characterization of over 1100 flight detectors, illustrating the consistent intrinsic low-noise performance of a large sample of GAPS detectors. This work demonstrates the feasibility of large-area and low-cost Si(Li) detector arrays for next-generation astrophysics and nuclear physics applications.

Accurate nuclear reaction rates for 26P(p,{\gamma})27S are pivotal for a comprehensive understanding of rp-process nucleosynthesis path in the region of proton-rich sulfur and phosphorus isotopes. However, large uncertainties still exist in the current rate of 26P(p,{\gamma})27S because of the lack of the nuclear mass and the energy level structure information of 27S. We reevaluate this reaction rate using the experimentally constrained 27S mass, together with the shell-model predicted level structure. It is found that the 26P(p,{\gamma})27S reaction rate is dominated by a direct-capture (DC) reaction mechanism despite the presence of three resonances at E = 1.104, 1.597, 1.777 MeV above the proton threshold in 27S. The new rate is overall smaller than the other previous rates from Hauser-Feshbach statistical model by at least one order of magnitude in the temperature range of X-ray burst interest. In addition, we consistently update the photodisintegration rate using the new 27S mass. The influence of new rates of forward and reverse reaction in the abundances of isotopes produced in rp-process is explored by post-processing nucleosynthesis calculations. The final abundance ratio of 27S/26P obtained using the new rates is only 10% of that from the old rate. The abundance flow calculations show the reaction path 26P(p,{\gamma})27S(\b{eta}+,{\nu})27P is not as important as thought previously for producing 27P. The adoption of the new reaction rates for 26P(p,{\gamma})27S only reduces the final production of aluminum by 7.1%, and has no discernible impact on the yield of other elements.

We implement the Bayesian inference to retrieve energy spectra of all neutrinos from a galactic core-collapse supernova (CCSN). To achieve high statistics and full sensitivity to all flavours of neutrinos, we adopt a combination of several reaction channels from different large-scale neutrino observatories, namely inverse beta decay on proton and elastic scattering on electron from Hyper-Kamiokande (Hyper-K), charged current absorption on Argon from Deep Underground Neutrino Experiment (DUNE) and coherent elastic scattering on Lead from RES-NOVA. Assuming no neutrino oscillation or specific oscillation models, we obtain mock data for each channel through Poisson processes with the predictions, for a typical source distance of 10 kpc in our Galaxy, and then evaluate the probability distributions for all spectral parameters of theoretical neutrino spectrum model with Bayes' theorem. Although the results for either the electron-neutrinos or electron-antineutrinos reserve relatively large uncertainties (according to the neutrino mass hierarchy), a precision of a few percent (i.e., $\pm 1 \% \sim \pm 4 \%$ at a credible interval of $2 \sigma$) is achieved for primary spectral parameters (e.g., mean energy and total emitted energy) of other neutrino species. Moreover, the correlation coefficients between different parameters are computed as well and interesting patterns are found. Especially, the mixing-induced correlations are sensitive to the neutrino mass hierarchy, which potentially makes it a brand new probe to determine the neutrino mass hierarchy in the detection of galactic supernova neutrinos. Finally, we discuss the origin of such correlation patterns and perspectives for further improvement on our results.

Models of induced-gravity inflation are formulated within Supergravity employing as inflaton the Higgs field which leads to a spontaneous breaking of a U(1)_{B-L} symmetry at Mgut=2x10^16 GeV. We use a renormalizable superpotential, fixed by a U(1) R symmetry, and logarithmic or semi-logarithmic Kahler potentials with integer prefactors which exhibit a quadratic non-minimal coupling to gravity. We find inflationary solutions of Starobinsky type in accordance with the observations. The inflaton mass is predicted to be of the order of 10^13 GeV. The model can be nicely linked to MSSM offering an explanation of the magnitude of the mu parameter consistently with phenomenological data. Also it allows for baryogenesis via non-thermal leptogenesis, provided that the gravitino is heavier than about 10 TeV.

In big bang nucleosynthesis (BBN), the deuterium-tritium (DT) fusion reaction, D(T,n)$\alpha$, enhanced by the 3/2$^+$ resonance, is responsible for 99% of primordial $^4$He. This has been known for decades and has been well documented in the scientific literature. However, following the tradition adopted by authors of learned articles, it was stated in a matter-of-fact manner and not emphasized; for most people, it has remained unknown. This helium became a source for the subsequent creation of $\geq$25\% of the carbon and other heavier elements and, thus, a substantial fraction of our human bodies. (To be more precise than $\geq$25\% will require future simulation studies on stellar nucleosynthesis.) Also, without this resonance, controlled fusion energy would be beyond reach. For example, for inertial confinement fusion (ICF), laser energy delivery for the National Ignition Facility (NIF) would have to be approximately 70 times larger for ignition. Because the resonance enhances the DT fusion cross section a hundredfold, we propose that the 3/2$^+$ $^5$He excited state be referred to as the "Bretscher state" in honor of the Manhattan Project scientist who discovered it, in analogy with the well-known 7.6 MeV "Hoyle state" in $^{12}$C that allows for the resonant 3$\alpha$ formation.

Hawking's quasi-local energy definition quantifies the energy enclosed by a spacelike 2-sphere in terms of the amount of lightbending on the sphere caused by the energy distribution inside the sphere. This paper establishes for the first time a direct connection between the formal mathematical definition of a quasi-local energy and observations, in the context of cosmological perturbation theory. This is achieved by studying the Hawking Energy of spherical sections of the past lightcone of a cosmic observer in a perturbed Friedmann-Lema\^{i}tre spacetime. We express the Hawking Energy in terms of gauge-invariant perturbation variables and comment on the cosmic observables needed to in principle measure it. We then calculate its angular power spectrum and interpret its contributions.

Gravitational waves can generate electromagnetic effects inside a strong electric or magnetic field within the Standard Model and general relativity. Here we propose using a quarterly split cavity and LC-resonance circuit to detect a high-frequency gravitational wave from 0.1 MHz to GHz. We perform a full 3D simulation of the cavity's signal for sensitivity estimate. Our sensitivity depends on the coherence time scale of the high-frequency gravitational wave sources and the volume size of the split cavity. We discuss the resonant measurement schemes for narrow-band gravitational wave sources and also a non-resonance scheme for broadband signals. For a meter-sized split cavity under a 14 Tesla magnetic field, the LC resonance enhanced sensitivity to the gravitational wave strain is expected to reach $h\sim 10^{-20}$ around $10$ MHz.