Papers for Monday, Dec 05 2022

The key to the exploration of the Ice Giant planets is avoiding cutting edge technology. Complexity produces delay and financial roadblocks. Simple robot scouts can be launched in time to utilize gravity assists from Jupiter in the early 2030s. Demands on NASA's budget from large missions, such as Mars sample return, will not allow Flagship missions to Uranus and Neptune in the near term. The science goals of Ice Giant exploration can be accomplished by a series of fast, simple, affordable (FSA) craft. Separate lines of cost-capped Orbiters and Probes would be launched at a cadence dictated by trajectories and funding. Contractors would be selected using competitive Announcements of Opportunity (AO). The march of progress in spacecraft technology offers hope and a path forward. The key is to start small and keep it affordable.

The study of flaring astrophysical events in the multi-messenger approach requires instantaneous follow-up observations to better understand the nature of these events through complementary observational data. We present Astro-COLIBRI as a meta-platform for the patchwork of different specific tools in the real-time multi-messenger ecosystem. The Astro-COLIBRI platform bundles and evaluates alerts about transients from various channels. It further automates the coordination of follow-up observations by providing and linking detailed information through its comprehensible graphical user interface. We present the functionalities using documented examples of Astro-COLIBRI usage through the community since its public release in August 2021. We highlight the use cases of Astro-COLIBRI for planning follow-up observations by professional and amateur astronomers, as well as checking predictions from theoretical models.

The majority of the Milky Way's stellar halo consists of debris from our Galaxy's last major merger, the Gaia-Sausage-Enceladus (GSE). In the past few years, stars from GSE have been kinematically and chemically studied in the inner $30$ kpc of our Galaxy. However, simulations predict that accreted debris could lie at greater distances, forming substructures in the outer halo. Here we derive metallicities and distances using Gaia DR3 XP spectra for an all-sky sample of luminous red giant stars, and map the outer halo with kinematics and metallicities out to $100$ kpc. We obtain follow-up spectra of stars in two strong overdensities - including the previously identified Outer Virgo Overdensity - and find them to be relatively metal-rich and on predominantly retrograde orbits, matching predictions from simulations of the GSE merger. We argue that these are apocentric shells of GSE debris, forming $60-90$ kpc counterparts to the $15-20$ kpc shells that are known to dominate the inner stellar halo. Extending our search across the sky with literature radial velocities, we find evidence for a coherent stream of retrograde stars encircling the Milky Way from $50-100$ kpc, in the same plane as the Sagittarius stream but moving in the opposite direction. These are the first discoveries of distant and structured imprints from the GSE merger, cementing the picture of an inclined and retrograde collision that built up our Galaxy's stellar halo.

Binary systems are ubiquitous and their formation requires two-body interaction and dissipation. In gaseous media, interactions between two initially unbound objects could result in gas-assisted binary formation, induced by a loss of kinetic energy to the ambient gas medium. Here we use analytic tools to derive the criteria for gas-assisted binary capture through gas dynamical friction dissipation. We validate them with few-body simulations and explore this process in different gas-rich environments, including gas-embedded stars in star forming regions (SFR), gas-enriched globular clusters, AGN disks and gaseous protoplanetary disks. We find that gas-assisted binary capture could be highly efficient in SFRs, potentially providing the main channel for the formation of stellar binaries. It could also operate under certain conditions in gas-enriched globular clusters. AGN disks could also provide a fertile ground for gas-assisted binary capture and in particular the formation of black-hole and other compact object binaries and the production of gravitational-wave (GW) and other high-energy transients, under the assumption of very thin AGN disks. Large scale gaseous disks might be too thick to enable gas-assisted binary capture and previous estimates of the production of GW-sources could be significantly overestimated, and in any case, sensitive to the specific condition and structure of the disks. In protoplanetary disks, however, while gas-assisted binary capture can produce binary KBOs, dynamical friction by small planetsimals is likely to be more efficient. Overall, we show that gas-assisted binary formation is robust and can contribute significantly to the binary formation rate in many environments. In fact, the gas-assisted binary capture rates are sufficiently high such that they will lead to multicaptures, and the formation of higher multiplicity systems.

The ultra-faint dwarf galaxy Reticulum II (Ret II) exhibits a unique chemical evolution history, with 72 +10/-12% of its stars strongly enhanced in r-process elements. We present deep Hubble Space Telescope photometry of Ret II and analyze its star formation history. As in other ultra-faint dwarfs, the color-magnitude diagram is best fit by a model consisting of two bursts of star formation. If we assume that the bursts were instantaneous, then the older burst occurred around the epoch of reionization and formed ~80% of the stars in the galaxy, while the remainder of the stars formed ~3 Gyr later. When the bursts are allowed to have nonzero durations we obtain slightly better fits. The best-fitting model in this case consists of two bursts beginning before reionization, with approximately half the stars formed in a short (100 Myr) burst and the other half in a more extended period lasting 2.6 Gyr. Considering the full set of viable star formation history models, we find that 28% of the stars formed within 500 +/- 200 Myr of the onset of star formation. The combination of the star formation history and the prevalence of r-process-enhanced stars demonstrates that the r-process elements in Ret II must have been synthesized early in its initial star-forming phase. We therefore constrain the delay time between the formation of the first stars in Ret II and the r-process nucleosynthesis to be less than 500 Myr. This measurement rules out an r-process source with a delay time of several Gyr or more such as GW170817.

The first JWST observations of nearby galaxies have unveiled a rich population of bubbles that trace the stellar feedback mechanisms responsible for their creation. Studying these bubbles therefore allows us to chart the interaction between stellar feedback and the interstellar medium, and the larger galactic flows needed to regulate star formation processes globally. We present the first catalog of bubbles in NGC628, visually identified using MIRI F770W PHANGS-JWST observations, and use them to statistically evaluate bubble characteristics. We classify 1694 structures as bubbles with radii between 6-552 pc. Of these, 31% contain at least one smaller bubble at their edge, indicating that previous generations of star formation have a local impact on where new stars form. On large scales, most bubbles lie near a spiral arm, and their radii increase downstream compared to upstream. Furthermore, bubbles are elongated in a similar direction to the spiral arm ridge-line. These azimuthal trends demonstrate that star formation is intimately connected to the spiral arm passage. Finally, the bubble size distribution follows a power-law of index $p=-2.2\pm0.1$, which is slightly shallower than the theoretical value by 1-3.5$\sigma$ that did not include bubble mergers. The fraction of bubbles identified within the shells of larger bubbles suggests that bubble merging is a common process. Our analysis therefore allows us to quantify the number of star-forming regions that are influenced by an earlier generation, and the role feedback processes have in setting the global star formation rate. With the full PHANGS-JWST sample, we can do this for more galaxies.

We present a high-resolution view of bubbles within The Phantom Galaxy (NGC 628); a nearby (~10Mpc), star-forming (~2Msun/yr), face-on (i~9deg) grand-design spiral galaxy. With new data obtained as part of the PHANGS-JWST treasury program, we perform a detailed case-study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (The Phantom Void; over 1kpc in diameter), and the other being a smaller region that may be the precursor to such a large bubble (The Precursor Phantom Void). When comparing to matched resolution Halpha observations from the Hubble Space Telescope (HST), we see that the ionized gas is brightest in the shells of both bubbles, and is coincident with the youngest (~1Myr) and most massive (~100,000Msun) stellar associations. We also find an older generation (~20Myr) of stellar associations is present within the bubble of The Phantom Void. From our kinematic analysis of the HI, H2 (CO) and HII gas across The Phantom Void, we infer a high expansion speed of around 15 to 50km/s. The large size and high expansion speed of The Phantom Void suggest that the driving mechanism is sustained stellar feedback due to multiple mechanisms, where early feedback first cleared a bubble (as we observe now in The Precursor Phantom Void), and since then SNe have been exploding within the cavity, and have accelerated the shell. Finally, comparison to simulations shows a striking resemblance to our JWST observations, and suggests that such large-scale stellar feedback-driven bubbles should be common within other galaxies.

New HCN, HNC, and HCO+ measurements of 46 normal galaxies in transitions up to J=4-3 are included in a multitransition database covering HCN and HCO+ (130 galaxies) and HNC (94 galaxies). The near-linear luminosity relations are dominated by distance effects and do not reflect galaxy physical properties. Individual galaxies show significant dispersion in both their luminosity and line ratio. Line ratios do not correlate well with either luminosities or other line ratios. Only the normalized J-transition ladders of HCN and HCO+ and the J=1-0 12CO/13CO isotopologue ratio are positively correlated with CO and far infrared (FIR) luminosity. The HCN and HCO+ molecules have very similar intensities and trace the same gas. In galaxies dominated by an active nucleus, HCO+ intensities are depressed relative to HCN intensities. Only a small fraction of the CO emission is directly associated with gas emitting in HCN and HCO, yet a significant fraction of even that gas appears to be translucent molecular gas. In the observed galaxy centers, the HCN/CO line intensity ratio is not a proxy for the dense gas fraction. Likewise, the FIR/HCN and FIR/CO ratios are not proxies for the star formation efficiency. The observed molecular line emission is fully consistent with UV-photon heating boosted by significant mechanical heating. The molecular gas sampled by HCN and HCO+ has low kinetic temperatures T(kin)=10-50 K, low densities n(H)=10^4-10^5 cm^(-3), and low optical depths in the ground-state lines. Most of the gas sampled by CO has densities lower by one to two orders of magnitude. For a mechanical heating fraction of 0.5, a modest energy input of only G=300 Go is required. A proper understanding of star formation requires a more appropriate determination of the gas mass than provided by the intensities of individual HCN or CO transitions.

We present a variety of analyses of the turbulent dynamics of the boundary of a photo-dissociation region (PDR) in the Carina Nebula using high resolution ALMA observations. Using Principal Component Analysis we suggest that the turbulence in this molecular cloud is driven at large scales. Analysis of the centroid velocity structure functions indicate that the turbulence is dominated by shocks rather than local (in k-space) transport of energy. We further find that length-scales in the range 0.02 - 0.03 pc are important in the dynamics of this cloud and this finding is supported by analysis of the dominant emission structure length-scale. These length-scales are well resolved by the observational data and we conclude that the apparent importance of this range of scales is physical in origin. Given that it is also well within the range strongly influenced by ambipolar diffusion, we conclude that it is not primarily a product of turbulence alone, but is more likely to be a result of the interplay between gravity and turbulence. Finally, through comparison of these results with previous observations of H2 emission from the Western Wall we demonstrate that observations of a PDR can be used to probe the internal structure of the undisturbed portion of a molecular cloud.

We have performed detailed high-resolution spectroscopic analysis on seven metal-poor stars (BD+75 348, BD+09 3019, HD238020, HE0319-0215, HE0507-1653, HE0930-0018, HE1023-1504) and derived their atmospheric parameters T$_{eff}$, log$g$, [Fe/H], and mictroturbulent velocity ($\xi$). The metallicity range is found to be -2.57$<$[Fe/H]$<$-0.42. The elemental abundances of 17 light elements and 12 heavy elements are estimated. We have classified BD+75 348 and BD+09 3019 as strong Ba stars, HD238020 as a mild Ba star, and the remaining four objects as CEMP-s stars. We have estimated the masses of the stars from Hertzsprung-Russel (HR) diagram, and, compiling the data of 205 Ba stars from literature, estimated the mass distribution of Ba stars. We have also estimated the initial masses of the companion AGBs of the program stars as well as the masses of the companion AGBs of 159 Ba and 36 CEMP-s stars from literature, with the help of a parametric-model-based analysis using FRUITY models. While the primary mass distribution of mild Ba stars peaks at 3.7 M$_{\odot}$, for the strong Ba stars the peak appears at 2.5 M$_{\odot}$. We, therefore, propose that the initial masses of the progenitor AGBs dominantly control the formation of mild and strong Ba stars. However, a clear overlap, in the range 1.3-4.0 M$_{\odot}$, noticed between the progenitor masses of both the subclasses of Ba stars, may indicate that other factors, such as metallicities and the orbital periods, may also have significant contributions. The progenitor AGBs' mass distribution of CEMP-s stars is found to peak at 2.03 M$_{\odot}$.

The fireball of the first interstellar meteor, CNEOS 2014-01-08 (IM1) (Siraj & Loeb 2019), was detected off the northern coast of Papua New Guinea. A recently announced ocean expedition will retrieve any extant fragments by towing a magnetic sled across a 10 km x 10 km area of ocean floor approximately 300 km north of Manus Island (Siraj, Loeb, & Gallaudet 2022). We formulate a model that includes both the probabilistic mass distribution of meteor fragments immediately after the fragmentation event, the ablation of the fragments, and the geographic distribution of post-ablation fragments along the ground track trajectory of the bulk fragment cloud. We apply this model to IM1 to provide a heuristic estimate of the impactor's post-ablation fragment mass distribution, constructed through a Monte Carlo simulation. We find between ~14% and ~36% of IM1 fragments are expected to survive ablation with a mass $\geq$ .001 g, and also provide an estimation for the geographic distribution of post-ablation fragments.

Observational studies have widely demonstrated that galaxy physical properties are strongly affected by the surrounding environment. On one side, gas inflows provide galaxies with new fuel for star formation. On the other side, the high temperatures and densities of the medium are expected to induce quenching in the star formation. Observations of large structures, in particular filaments at the cluster outskirts (r>2r$_{200}$), are currently limited to the low redshift Universe. We present a multi-band dataset for the cluster MACS J0416.1-2403 (z=0.397), observed in the context of the Galaxy Assembly as a function of Mass and Environment with VST (VST-GAME) survey. The project aims at gathering deep ($r$<24.4) and wide (20x20Mpc$^2$) observations at optical wavelengths for six massive galaxy clusters at 0.2<z<0.6, complemented with near infrared data. This work describes the photometric analysis of the cluster, defines a density field and studies galaxy properties in the cluster outskirts. We extract sources paying particular attention to recover the faintest ones. We combine all the extractions in a multi-band catalog and compute photometric redshifts. We then define cluster memberships up to 5r$_{200}$ from the cluster core and measure the density field, comparing galaxy properties in different environments. We found that the $g-r$ colors show bimodal behaviours in all the environments, but the peak of the distribution of red galaxies shifts toward redder colors with increasing density and the fraction of galaxies in the blue cloud increases with decreasing density. We also found 3 overdense regions in the cluster outskirts at r$\sim$5r$_{200}$. The color of galaxies suggests the presence of evolved galaxy populations, an insight for pre-processing phenomena over these substructures. We release the multi-band catalog, down to the completeness limit $r$<24.4 mag.

We model substructure in the protoplanetary disks around DoAr 44 and HD 163296 in order to better understand the conditions under which planets may form. We match archival millimeter-wavelength thermal emission against models of the disks' structure that are in radiation balance with the starlight heating and in vertical hydrostatic equilibrium, and then compare to archival polarized scattered near-infrared images of the disks. The millimeter emission arises in the interior, while the scattered near-infrared radiation probes the disks' outer layers. Our best model of the HD 163296 disk has dust masses $81\pm 13$ $M_\oplus$ in the inner ring at 68 au and $82^{+26}_{-16}$ $M_\oplus$ in the outer ring at 102 au, both falling within the range of estimates from previous studies. Our DoAr 44 model has total dust mass $84^{+7.0}_{-3.5}$ $M_\oplus$. Unlike HD 163296, DoAr 44 as of yet has no detected planets. If the central cavity in the DoAr 44 disk is caused by a planet, the planet's mass must be at least 0.5 $M_J$ and is unlikely to be greater than 1.6 $M_J$. We demonstrate that the DoAr 44 disk's structure with a bright ring offset within a fainter skirt can be formed by dust particles drifting through a plausible distribution of gas.

At the 450 yr anniversary of its observation, the supernova named after Tycho Brahe, SN 1572, can be explained in the terms used nowadays to characterize Type supernovae (SNe Ia). By assembling the records of the observations made in 1572--74 and evaluating their uncertainties, it is possible to recover the light curve and the color evolution of this supernova. It is found that, within the SNe Ia family, the event should have been a SN Ia with a normal rate of decline. Concerning the color evolution of SNe Ia, the most recently recovered records reaffirm previous findings of its being a normal SN Ia. Since 2004, an exploration for the progenitor of this SN Ia has tentatively suggested that the binary path to the explosion was that of a single degenerate. However, this is not conclusive at this point. The abundance studies from X--ray spectroscopy of the whole remnant point to a nuclear burning of the kind of a delayed detonation explosion of a Chandrasekhar--mass white dwarf. This was the first supernova studied by astronomers and it is still subject of very active research.

Planets and moons reorient in space due to mass redistribution associated with various types of internal and external processes. While the equilibrium orientation of a tidally locked body is well understood, much less explored are the dynamics of the reorientation process (or true polar wander, TPW, used here for the motion of either the rotation or the tidal pole). This is despite their importance for predicting the patterns of TPW-induced surface fractures, and for assessing whether enough time has passed for the equilibrium orientation to be reached. The only existing, and relatively complex numerical method for an accurate evaluation of the reorientation dynamics of a tidally locked body was described in a series of papers by Hu et al. (2017a,b, 2019). Here we demonstrate that an identical solution can be obtained with a simple approach, denoted as ow||mMIA, because, contrary to previous claims, during TPW the tidal and the rotation axes closely follow respectively the minor and the major axes of the total, time-evolving inertia tensor. Motivated by the presumed reorientation of Pluto, the use of the ow||mMIA method is illustrated on several test examples. In particular, we analyze whether reorientation paths are curved or straight when the load sign and the mass of the host body are varied. When tidal forcing is relatively small, the paths of negative anomalies (e.g.~basins) towards the rotation pole are highly curved, while positive loads may reach the sub- or anti-host point in a straightforward manner. Our results suggest that the Sputnik Planitia basin cannot be a negative anomaly at present day, and that the remnant figure of Pluto must have formed prior to the reorientation. Finally, the presented method is complemented with an energy balance that can be used to test the numerical solution.

Gamma rays are produced by cosmic ray (CR) protons interacting with the particles at solar photosphere and by cosmic ray electrons and positrons (CRes) via inverse Compton scattering of solar photons. The former come from the solar disk while the latter extend beyond the disk. Evaluation of these emissions requires the flux and spectrum of CRs in the vicinity of the Sun, while most observations provide flux and spectra near the Earth, at around 1 AU from the Sun. Past estimates of the quiet Sun gamma-ray emission use phenomenological modulation procedures to estimate spectra near the Sun (see review by Orlando and Strong 2021 and references therein). We show that CRe transport in the inner heliosphere requires a kinetic approach and use a novel approximation to determine the variation of CRe flux and spectrum from 1 AU to the Sun including effects of (1) the structure of large scale magnetic field, (2) small scale turbulence in the solar wind from several in situ measurements, in particular, those by Parker Solar Probe that extend this information to 0.1 AU, and (3) most importantly, energy losses due to synchrotron and inverse Compton processes. We present results on the flux and spectrum variation of CRes from 1 AU to the Sun for several transport models. In forthcoming papers we will use these results for a more accurate estimate of quiet Sun inverse Compton gamma-ray spectra, and, for the first time, the spectrum of extreme ultraviolet to hard X-ray photons produced by synchrotron emission. These can be compared with the quiet Sun gamma-ray observation by Fermi (see, e.g.~Fermi-LAT Collaboration, 2011) and X-ray upper limits set by RHESSI (Hannah et al., 2010).

The butterfly diagram of the solar cycle is the equatorward migration of the emergence latitudes of sunspots as the solar cycle evolves. Revealing the mechanism for the butterfly diagram is essential for understanding the solar and stellar dynamo. The equatorward meridional flow at the base of the convection zone (CZ) was believed to be responsible for the butterfly diagram. However, helioseismological studies indicate controversial forms of the flow, and even present poleward flow at the base of the CZ, which poses a big challenge to the widely accepted mechanism. This motivates us to propose a new mechanism in this study. Using a data-driven Babcock-Leighton-type dynamo model, we carry out numerical simulations to explore how the latitude-dependent radial flux transport affects the latitudinal migration of the toroidal field, under different meridional flow profiles. The results indicate that when the radial transport of the poloidal field at higher latitudes is sufficiently faster, the toroidal fields of a new cycle at higher latitudes are generated earlier than that at lower latitudes, and vice versa. Thus, the butterfly diagram is suggested to correspond to the time- and latitude-dependent regeneration of the toroidal field due to the latitude-dependent radial transport of the poloidal flux.

Stellar streams are sensitive probes of the Galactic potential. The likelihood of a stream model given stream data is often assessed using simulations. However, comparing to simulations is challenging when even the stream paths can be hard to quantify. Here we present a novel application of Self-Organizing Maps and first-order Kalman Filters to reconstruct a stream's path, propagating measurement errors and data sparsity into the stream path uncertainty. The technique is Galactic-model independent, non-parametric, and works on phase-wrapped streams. With this technique, we can uniformly analyze and compare data with simulations, enabling both comparison of simulation techniques and ensemble analysis with stream tracks of many stellar streams. Our method is implemented in the public Python package TrackStream, available at https://github.com/nstarman/trackstream.

The bright end of the rest-frame UV luminosity function (UVLF) of high-redshift galaxies is modified by gravitational lensing magnification bias. Motivated by recent discoveries of very high-z galaxies with JWST, we study the dependence of magnification bias on the finite size of sources at $6<z<14$. We calculate the magnification probability distributions and use these to calculate the magnification bias assuming a rest-frame Schechter UVLF for galaxies at redshift $6<z<14$. We find that the finite size of bright high-redshift galaxies together with lens ellipticity significantly suppresses magnification bias, producing an observed bright end which declines more sharply than the power-law resulting from assumption of point sources. By assuming a luminosity-size relation for the source population and comparing with the observed $z=6$ galaxy luminosity function from Harikane+(2022), we show that the UVLF can be used to set mild constraints on the galaxies intrinsic size, favoring smaller galaxies compared to the fiducial luminosity-size relation. In the future, wide surveys using Euclid and Roman Space Telescope will place stronger constraints. We also tabulate the maximum magnification possible as a function of source size and lens ellipticity.

The identification and localization of Fast Radio Bursts to their host galaxies has revealed important details about the progenitors of these mysterious, millisecond-long bursts of coherent radio emission. In this work we study the most probable host galaxy of the apparently non-repeating CHIME/FRB event FRB 20190425A -- a particularly high luminosity, low dispersion measure event that was demonstrated by Moroianu et al. 2022 to be temporally and spatially coincident with the LIGO-Virgo-KAGRA binary neutron star merger GW190425, suggesting an astrophysical association (p-value 0.0052). In this paper we remain agnostic to this result, and we confirm UGC10667 as the most probable host galaxy of FRB 20190425A, demonstrating that the host galaxies of low dispersion measure, one-off CHIME FRBs can be plausibly identified. We then perform multi-wavelength observations to characterize the galaxy and search for any afterglow emission associated with the FRB and its putative GW counterpart. We find no radio or optical transient emission in our observations $2.5\,\mathrm{yr}$ post-burst. UGC10667 is a spiral galaxy at $z\sim0.03$, dominated by an old stellar population. We find no evidence of a large population of young stars, with nebular emission dominated by star formation at a rate of $1-2\,\mathrm{M_\odot\,yr^{-1}}$. While we cannot rule out a young magnetar as the origin of FRB 20190425A, our observations are consistent with an origin in a long delay-time neutron star binary merger as posited by Moroianu et al. 2022.

Flare occurrence on the Sun is highly variable, exhibiting both short term variation due to the emergence and evolution of active regions, and long-term variation from the solar cycle. On solar-like stars, much larger stellar flares (superflares) have been observed, and it is of interest to determine whether observed rates of superflare occurrence exhibit similar variability to solar flares. We analyse 274 G-type stars using data from the Transiting Exoplanet Survey Satellite (TESS) and identify seven stars which exhibit statistically significant changes in the rate of superflare occurrence by fitting a piecewise constant-rate model with the Bayesian Blocks algorithm (Scargle et al 2012; arXiv:1207.5578). We investigate the properties of these stars and their flaring rates, and discuss the possible reasons for the low number of stars with detectable rate variation.

Recent observations demonstrate that the symbiotic X-ray binary (SyXB) IGR J17329-2731 contains a highly magnetized neutron star (NS) which accretes matter through the wind from its giant star companion, and suggest that 4U 1700+24 may also have a highly magnetized NS. Accretion-induced collapse (AIC) from oxygen-neon-magnesium white dwarf (ONeMg WD) + red giant (RG) star binaries is one promising channel to form these SyXBs, while other long standing formation channels have difficulties to produce these SyXBs. By considering non-magnetic and magnetic ONeMg WDs, I investigate the evolution of ONeMg WD + RG binaries with the MESA stellar evolution code for producing SyXBs with non-magnetic or magnetized NSs. In the pre-AIC evolution with magnetic confinement, the mass accumulation efficiency of the accreting WD is increased at low mass transfer rate compared to the non-magnetic case. The newborn NSs formed via AIC of highly magnetized WDs could inherit the large magnetic field through conservation of magnetic flux, and the systems could have a long age compatible with that of the red giant companions. These young and highly magnetized NSs could accrete matters from the stellar wind of the giant companions to that shine as those observed SyXBs, and could preserve their high magnetic field during this time. The MESA calculation results show that the initial parameter (initial RG mass and orbital period) space for the AIC with magnetic confinement to form SyXBs with highly magnetized NSs shifts to be lower and narrower compared to that of the no magnetic confinement case.

Future launches are projected to significantly increase both the number of active satellites and aggregate collision risk in Low Earth Orbit (LEO). In this paper, a dynamical systems theory approach is used to analyze the effect of launch rate distribution on the stability of the LEO environment. A multi-shell, three-species source-sink model of the LEO environment, referred to as MOCAT-3 for MIT Orbital Capacity Assessment Tool 3 Species, is used to study the evolution of the species populations. The three species included in the model are active satellites, derelict satellites, and debris. The model's coefficients represent atmospheric drag, collision rate, mean satellite lifetime, post-mission disposal probability, and active debris removal rate. Solutions of the system of differential equations are computed, and an analysis of the stability of the equilibrium points is conducted for numerous launch rate distributions. The stability of the equilibrium points is used to test the sensitivity of the environment to run-away debris growth, known as Kessler syndrome, that occurs at the instability threshold. An analysis of the environment's response to perturbations in launch rate and debris population is conducted. The maximum perturbation in the debris population from the equilibrium state, for which the system remains in a stable configuration, is calculated. Plots of the phase space about the equilibrium points are generated. The results will help to better understand the orbital capacity of LEO and the stability of the space environment, as well as provide improved guidelines on future launch plans to avoid detrimental congestion of LEO.

We report the results of broadband timing and spectral analysis of data from an AstroSat observation of the High Mass X-ray binary LMC X$-$4. The Large Area X-ray Proportional Counter (LAXPC) and Soft X-ray Telescope (SXT) instruments on-board the AstroSat observed the source in August 2016. A complete X-ray eclipse was detected with the LAXPC. The 3$-$40 keV power density spectrum showed the presence of coherent pulsations along with a $\sim 26$ mHz quasi-periodic oscillation feature. The spectral properties of LMC X$-$4 were derived from a joint analysis of the SXT and LAXPC spectral data. The 0.5$-$25 keV persistent spectrum comprised of an absorbed high energy cutoff power law with photon index of $\Gamma \sim$ 0.8 and cutoff at $\sim$16 keV, a soft thermal component with kT$_{BB} \sim$ 0.14 keV and Gaussian components corresponding to Fe K$_\alpha$, Ne \textsc{ix} and Ne \textsc{x} emission lines. Assuming a source distance of 50 kpc, we determined 0.5--25 keV luminosity to be $\sim 2 \times 10^{38}$ erg s$^{-1}$.

Starting from the 1970s, some relations connecting dark matter and baryons were discovered, such as the Tully-Fisher relation. However, many of the relations found in galaxies are quite different from that found in galaxy clusters. Here, we report two new mysterious universal relations connecting dark matter and baryons in both galaxies and galaxy clusters. The first relation indicates that the total dynamical mass of a galaxy or a galaxy cluster $M_{500}$ has a power-law relation with its total baryonic mass $M_b$ within the `virial region': $M_{500} \propto M_b^a$, with $a \approx 3/4$. The second relation indicates that the enclosed dynamical mass $M_d$ is almost directly proportional to the baryonic mass for galaxies and galaxy clusters within the central baryonic region: $M_d \propto M_b$. The close relations between dark matter and baryons in both galaxies and galaxy clusters suggest that some unknown interaction or interplay except gravity might exist between dark matter and baryons.

Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially sub-$\mu m$ dust particles in self-gravitating young protoplanetary disks with various strengths of turbulent viscosity. We aim to understand the physical conditions that determine the formation and spatial distribution of pebbles when both disk self-gravity and turbulent viscosity can be concurrently at work. Methods. We perform the thin-disk hydrodynamics simulations of self-gravitating protoplanetary disks over an initial time period of 0.5 Myr using the FEOSAD code. Turbulent viscosity is parameterized in terms of the spatially and temporally constant $\alpha$-parameter, while the effects of gravitational instability on dust growth is accounted for by calculating the effective parameter $\alpha_{\rm GI}$. We consider the evolution of dust component including momentum exchange with gas, dust self-gravity, and also a simplified model of dust growth. Results. We find that the level of turbulent viscosity strongly affects the spatial distribution and total mass of pebbles in the disk. The $\alpha=10^{-2}$ model is viscosity-dominated, pebbles are completely absent, and dust-to-gas mass ratio deviates from the reference 1:100 value no more than by 30\% throughout the disk extent. On the contrary, the $\alpha=10^{-3}$ model and, especially, the $\alpha=10^{-4}$ model are dominated by gravitational instability. The effective parameter $\alpha+\alpha_{\rm GI}$ is now a strongly varying function of radial distance. As a consequence, a bottle neck effect develops in the innermost disk regions, which makes gas and dust accumulate in a ring-like structure. Abridged.

We investigate the propagation of MHD fast waves into a cylindrical coronal loop through an inhomogeneous stationary flow region. The background flow is assumed to have a small, spatially periodic structure in addition to a constant speed. We focus on the absorption of the wave energy in Alfv\'{e}n resonance, comparing with the constant flow case. A new flow (absorption) regime is induced by the periodic flow structure which enhances the absorption for the antiparallel flow and inverse absorption (overreflection) for the parallel flow with respect to the axial wave vector, depending on the transitional layer and flow profiles. A giant overreflection and anomalous absorption behavior arise for some flow configurations. In the other flow regimes, its effect on the absorption is shown to be weak.

Isotropic diffuse gamma-ray flux in the PeV energy band is an important tool for multimessenger tests of models of the origin of high-energy astrophysical neutrinos and for new-physics searches. So far, this flux has not yet been observed. Carpet-2 is an air-shower experiment capable of detecting astrophysical gamma rays with energies above 0.1 PeV. Here we report the upper limits on the isotropic gamma-ray flux from Carpet-2 data obtained in 1999-2011 and 2018-2022. These results, obtained with the new statistical method based on the shape of the muon-number distribution, summarize Carpet-2 observations as the upgraded installation, Carpet-3, starts its operation.

On August 17, 2017, Advanced LIGO and Virgo observed GW170817, the first gravitational-wave (GW) signal from a binary neutron star merger. It was followed by a short-duration gamma-ray burst, GRB 170817A, and by a non-thermal afterglow emission. In this work, a combined simultaneous fit of the electromagnetic (EM, specifically, afterglow) and GW domains is implemented, both using the posterior distribution of a GW standalone analysis as prior distribution to separately process the EM data, and fitting the EM and GW domains simultaneously. These approaches coincide mathematically, as long as the actual posterior of the GW analysis, and not an approximation, is used as prior for the EM analysis. We treat the viewing angle, $\theta_v$, as shared parameter across the two domains. In the afterglow modelling this parameter and the jet core angle, $\theta_c$, are correlated, leading to high uncertainties on their values. The joint EM+GW analysis relaxes this degeneracy, reducing the uncertainty compared to an EM-only fit. We also apply our methodology to a hypothetical GW170817-like event occurring in the next GW observing run at $\sim$140 Mpc, so that the afterglow flux is one order of magnitude fainter, leaving only the light curve peak visible. The EM-only fit cannot constrain $\theta_v$ nor $\theta_c$, while folding the GW data into the analysis leads to tighter constraints on $\theta_v$. Although the geometry of the jet remains unconstrained because of the lack of data in the afterglow rising phase, in all fits the jet total width results to be smaller than $\theta_v$ (off-axis case).

Star formation can produce bubbles and outflows, as a result of stellar feedback. Outflows and bubbles inject momentum and energy into the surrounding interstellar medium, and so are related to the overall energy balance of the molecular cloud. Molecular bubbles can be resolved by higher-resolution radio telescopes to quantify the effect of star formation on molecular clouds. We report here the identification of a new molecular bubble with an outflow, and an Herbig Haro object, HH319, located at the bubble center. Multi-wavelength data have been utilized to study its spatial structure, energy injection, and dynamical timescale. This bubble has a kinetic energy of $\rm 5.8 \times 10^{43}$ erg within the smallest radius of a bubble in Taurus, 0.077 pc. The bubble formed $\sim$70,000 years ago. According to the proper motion velocities of protostars from $Gaia$ EDR3, the T Tauri binary stars (FY Tau and FZ Tau) at the southwest edge of the bubble may have produced the outflow-bubble structure. This is an unusual new structure found in low- and intermediate-mass star formation regions. Only a bubble in Orion A, driven by V380 Ori, has a similar structure. The bubble-outflow structure provides additional observational evidence for the theory of stellar wind from T Tauri stars. It enhances our understanding of how stellar feedback acts on molecular clouds.

Relativistic reflection observed in the hard states of accreting black holes often shows a weak amplitude relative to the main Comptonization component, which may result from either a disc truncation or a non-isotropy of the X-ray source, e.g. due to a motion away from the reflector. We investigate here the latter case, assuming that the X-ray source is located on the symmetry axis of the Kerr black hole. We discuss effects relevant to a proper computation of the reflected radiation and we implement them in the model for data analysis, reflkerrV. We apply it to the simultaneous Suzaku and NuSTAR observation of Cyg X-1 in the hard state and we find a good fit for an untruncated disc irradiated by the source moving away from it at 0.36c. However, we find a slightly better solution in a geometry closely approximating the truncated disc irradiated by an inner hot flow. In this solution we either still need a subrelativistic outflow or the source opposite to the observer must contribute to the directly observed radiation. We also discuss differences between the implementation of the outflow effect in reflkerrV and in relxilllpCp.

Physics-based solar cycle predictions provide an effective way to verify our understanding of the solar cycle. Before the start of cycle 25, several physics-based solar cycle predictions were developed. These predictions use flux transport dynamo (FTD) models, surface flux transport (SFT) models, or a combination of the two kinds of models. The common physics behind these predictions is that the surface poloidal fields around cycle minimum dominate the subsequent cycle strength. In the review, we first give short introductions to SFT and FTD models. Then we compare 7 physics-based prediction models from 4 aspects, which are what the predictor is, how to get the predictor, how to use the predictor, and what to predict. Finally, we demonstrate the large effect of assimilated magnetograms on predictions by two SFT numerical tests. We suggest that uncertainties in both initial magnetograms and sunspot emergence should be included in such physics-based predictions because of their large effects on the results. In addition, in the review we put emphasis on what we can learn from different predictions, rather than an assessment of prediction results.

We reconsider the gravitational wave spectrum induced by scalar perturbations in spatially flat Friedmann-Lema\^itre-Robertson-Walker spacetimes, focusing on the matter- and $\Lambda$-dominated epochs. During matter domination, sub-horizon modes are not free and a commonly applied approximation for the derivative of the tensor perturbation is flawed. We show analytically that this leads to a significant overestimation of the energy density spectrum. In addition, we demonstrate that gauge-dependent non-oscillating tensor perturbations appear in the presence of a cosmological constant. Complementing the analytical calculations, we compute the according present-day spectrum numerically for a Planck-like $\Lambda$CDM cosmology, finding that non-oscillating growing modes appear during the transition between matter and $\Lambda$ domination in conformal Newtonian gauge.

Context: Outflows from low-mass star-forming galaxies are a fundamental ingredient for models of galaxy evolution and cosmology. Aims: The onset of kpc-scale ionised filaments in the halo of the metal-poor compact dwarf SBS 0335-052E was previously not linked to an outflow. We here we investigate whether these filaments provide evidence for an outflow. Methods: We obtained new VLT/MUSE WFM and deep NRAO/VLA B-configuration 21cm data of the galaxy. The MUSE data provide morphology, kinematics, and emission line ratios H$\beta$/H$\alpha$ and [\ion{O}{iii}]$\lambda5007$/H$\alpha$ of the low surface-brightness filaments, while the VLA data deliver morphology and kinematics of the neutral gas in and around the system. Both datasets are used in concert for comparisons between the ionised and the neutral phase. Results: We report the prolongation of a lacy filamentary ionised structure up to a projected distance of 16 kpc at $\mathrm{SB}_\mathrm{H\alpha} = 1.5\times10^{-18}$erg s$^{-1}$ cm$^{-2}$arcsec$^{-2}$. The filaments exhibit unusual low H$\alpha$/H$\beta \approx 2.4$ and low [\ion{O}{iii}]/H$\alpha \sim 0.4 - 0.6$ typical of diffuse ionised gas. They are spectrally narrow ($\sim 20$ km s$^{-1}$) and exhibit no velocity sub-structure. The filaments extend outwards of the elongated \ion{H}{I} halo. On small scales the $N_\mathrm{HI}$ peak is offset from the main star-forming sites. Morphology and kinematics of \ion{H}{I} and \ion{H}{II} reveal how star-formation driven feedback interacts differently with the ionised and the neutral phase. Conclusions: We reason that the filaments are a large scale manifestation of star-formation driven feedback, namely limb-brightened edges of a giant outflow cone that protrudes through the halo of this gas-rich system. A simple toy model of such a conical-structure is found to be commensurable with the observations.

Intensity mapping with the redshifted 21-cm line is an emerging tool in cosmology. Drift scan observations, where the antennas are fixed to the ground and the telescope's pointing center (PC) changes continuously on the sky due to earth's rotation, provide broad sky coverage and sustained instrumental stability needed for 21-cm intensity mapping. Here we present the Tracking Tapered Grided Estimator (TTGE) to quantify the power spectrum of the sky signal estimated directly from the visibilities measured in drift scan radio interferometric observations. The TTGE uses the data from the different PC to estimate the power spectrum of the signal from a small angular region located around a fixed tracking center (TC). The size of this angular region is decided by a suitably chosen tapering window function which serves to reduce the foreground contamination from bright sources located at large angles from the TC. It is possible to cover the angular footprint of the drift scan observations using multiple TC, and combine the estimated power spectra to increase the signal to noise ratio. Here we have validated the TTGE using simulations of $154 \, {\rm MHz}$ MWA drift scan observations. We show that the TTGE can recover the input model angular power spectrum $C_{\ell}$ within $20 \%$ accuracy over the $\ell$ range $40 < \ell < 700$.

Mrk 477 is the closest type-2 quasar (QSO2), at a distance of 163 Mpc. This makes it an ideal laboratory for studying the interplay between nuclear activity and star formation with a great level of detail and signal-to-noise. In this Letter we present new mid-infrared (mid-IR) imaging and spectroscopic data with an angular resolution of 0.4 arcsec (~300 pc) obtained with the Gran Telescopio Canarias (GTC) instrument CanariCam. The N-band (8-13 micron) spectrum of the central ~400 pc of the galaxy reveals [S IV]10.51 micron emission, but no 8.6 or 11.3 micron PAH features, which are commonly used as tracers of recent star formation. This is in stark contrast with the presence of a nuclear starburst with age of about 6 Myr and ~300 pc in size, as constrained from ultraviolet HST observations. Considering this, we argue that even the more resilient, neutral molecules that mainly produce the 11.3 micron PAH band are most likely being destroyed in the vicinity of the active nucleus despite the relatively large X-ray column density, of log N_H=23.5 cm^-2, and modest X-ray luminosity, of 1.5x10^43 erg/s. This highlights the importance of being cautious when using PAH features as star formation tracers in the central region of galaxies to evaluate the impact of active galactic nuclei (AGN) feedback.

We investigate the ionization of the diffuse interstellar medium by cosmic rays by modeling their propagation along the wandering magnetic fields using a Monte Carlo method. We study how low-energy cosmic rays propagate in turbulent, translucent molecular clouds, and how they regulate the ionization and both lose and gain energy from the medium. As a test case, we use high spatial resolution (0.03 pc) CO maps of a well-studied high latitude translucent cloud, MBM 3, to model turbulence. The propagation problem is solved with a modified Monte Carlo procedure that includes trapping, energization, and ionization losses. In a homogeneous medium, trapping and re-energization do not produce a significant effect. In a nonuniform medium, particles can be trapped for a long time inside the cloud. This modifies the cosmic ray distribution due to stochastic acceleration at the highest energies (about 100 MeV). At lower energies, the re-energization is too weak to produce an appreciable effect. The change in the energy distribution does not significantly affect the ionization losses, so ionization changes are due to trapping effects.

Observations of protoplanetary disks provide information on planet formation and the reasons for the diversity of planetary systems. The key to understanding planet formation is the study of dust evolution from small grains to pebbles. Smaller grains are well-coupled to the gas dynamics, and their distribution is significantly extended above the disk midplane. Larger grains settle much faster and are efficiently formed only in the midplane. By combining near-infrared polarized light and millimeter observations, it is possible to constrain the spatial distribution of both the small and large grains. We aim to construct detailed models of the size distribution and vertical/radial structure of the dust particles in protoplanetary disks based on observational data. In particular, we are interested in recovering the dust distribution in the IM Lup protoplanetary disk. We create a physical model for the dust distribution of protoplanetary disks and simulate the radiative transfer of the millimeter continuum and the near-infrared polarized radiation. Using a Markov chain Monte Carlo method, we compare the derived images to the observations available for the IM Lup disk to constrain the best physical model for IM Lup and to recover the vertical grain size distribution. The millimeter and near-infrared emission tightly constrain the dust mass and grain size distribution of our model. We find size segregation in the dust distribution, with millimeter-sized grains in the disk midplane. These grains are efficiently formed in the disk, possibly by sedimentation-driven coagulation, in accord with the short settling timescales predicted by our model. This also suggests a high dust-to-gas ratio at smaller radii in the midplane, possibly triggering streaming instabilities and planetesimal formation in the inner disk.

The code FADO is the first publicly available population spectral synthesis tool that treats the contribution from ionised gas to the observed emission self-consistently. We study the impact of the nebular contribution on the determination of the star formation rate (SFR), stellar mass, and consequent effect on the star-forming main sequence (SFMS) at low redshift. We applied FADO to the spectral database of the SDSS to derive the physical properties of galaxies. As a comparison, we used the data in the MPA-JHU catalogue, which contains the properties of SDSS galaxies derived without the nebular contribution. We selected a sample of SF galaxies with H$\alpha$ and H$\beta$ flux measurements, and we corrected the fluxes for the nebular extinction through the Balmer decrement. We then calculated the H$\alpha$ luminosity to estimate the SFR. Then, by combining the stellar mass and SFR estimates from FADO and MPA-JHU, the SFMS was obtained. The H$\alpha$ flux estimates are similar between FADO and MPA-JHU. Because the H$\alpha$ flux was used as tracer of the SFR, FADO and MPA-JHU agree in their SFR. The stellar mass estimates are slightly higher for FADO than for MPA-JHU on average. However, considering the uncertainties, the differences are negligible. With similar SFR and stellar mass estimates, the derived SFMS is also similar between FADO and MPA-JHU. Our results show that for SDSS normal SF galaxies, the additional modelling of the nebular contribution does not affect the retrieved fluxes and consequentially also does not influence SFR estimators based on the extinction-corrected H$\alpha$ luminosity. For the stellar masses, the results point to the same conclusion. These results are a consequence of the fact that the vast majority of normal SF galaxies in the SDSS have a low nebular contribution.

Recent news reports claim that China is overtaking the United States and all other countries in scientific productivity and scientific impact. A straightforward analysis of high-impact papers in astronomy reveals that this is not true in our field. In fact, the United States continues to host, by a large margin, the authors that lead high-impact papers. Moreover, this analysis shows that 90% of all high-impact papers in astronomy are led by authors based in North America and Europe. That is, only about 10% of countries in the world host astronomers that publish "astronomy's greatest hits".

The prevalent view that the radio-loud gamma-ray pulsars have gamma-ray luminosities that exceed their radio luminosities by several orders of magnitude is based on the assumption that the decay with distance of their gamma-ray fluxes obeys the inverse-square law as does that of their radio fluxes. The results presented here, of testing the hypothesis of independence of luminosities and distances of gamma-ray pulsars by means of the Efron-Petrosian statistic, do not uphold this assumption however: they imply that the observational data in the Fermi-LAT 12-Year Catalog are consistent with the dependence $S\propto D^{-3/2}$ of the flux densities $S$ of the gamma-ray pulsars on their distances $D$ at substantially higher levels of significance than they are with the dependence $S\propto D^{-2}$. These results, which were theoretically predicted in Ardavan (2021, MNRAS, 507, 4530), are not incompatible with the requirements of the conservation of energy because the radiation process by which the superluminally moving current sheet in the magnetosphere of a neutron star has been shown to generate the slowly decaying gamma-ray pulses is intrinsically transient: the difference in the fluxes of power across any two spheres centred on the star is balanced by the change with time of the energy contained inside the shell bounded by those spheres. Once the over-estimation of their values is rectified, the luminosities of gamma-ray pulsars turn out to have the same range of values as do the luminosities of radio pulsars. This conclusion agrees with that reached earlier on the basis of the smaller data set in the Second Fermi-LAT Catalog of Gamma-ray Pulsars.

Many aspects of the three-dimensional (3-D) structure and evolution of interplanetary coronal mass ejections (ICMEs) remain unexplained. Here, we investigate two main topics: (1) the coherence scale of magnetic fields inside ICMEs, and (2) the dynamic nature of ICME magnetic complexity. We simulate ICMEs interacting with different solar winds using the linear force-free spheromak model incorporated into the EUHFORIA model. We place a swarm of ~20000 spacecraft in the 3-D simulation domain and characterize ICME magnetic complexity and coherence at each spacecraft based on simulated time series. Our simulations suggest that ICMEs retain a lower complexity and higher coherence along their magnetic axis, but that a characterization of their global complexity requires crossings along both the axial and perpendicular directions. For an ICME of initial half angular width of $45^\circ$ that does not interact with other large-scale solar wind structures, global complexity can be characterized by as little as 7-12 spacecraft separated by $25^\circ$, but the minimum number of spacecraft rises to 50-65 (separated by $10^\circ$) if interactions occur. Without interactions, ICME coherence extends for $45^\circ$, $20^\circ$-$30^\circ$, $15^\circ$-$30^\circ$, and $0^\circ$-$10^\circ$ for $B$, $B_\phi$, $B_\theta$, and $B_r$, respectively. Coherence is also lower in the ICME west flank compared to the east flank due to Parker spiral effects. Moreover, coherence is reduced by a factor of 3-6 by interactions with solar wind structures. Our findings help constrain some of the critical scales that control the evolution of ICMEs and aid in the planning of future dedicated multi-spacecraft missions.

Pseudo two-colour diagrams or Chromosome maps (ChM) indicate that NGC 2808 host five different stellar populations. The existing ChMs have been derived by the Hubble Space Telescope photometry, and comprise of stars in a small field of view around the cluster centre. To overcome these limitations, we built a ChM with U,B,I photometry from ground-based facilities that disentangle the multiple stellar populations of NGC 2808 over a wider field of view. We used spectra collected by GIRAFFE@VLT in a sample of 70 red giant branch (RGB) and seven asymptotic giant branch (AGB) stars to infer the abundances of C, N, O, Al, Fe, and Ni, which combined with literature data for other elements (Li, Na, Mg, Si, Ca, Sc, Ti, Cr and Mn), and together with both the classical and the new ground-based ChMs, provide the most complete chemical characterisation of the stellar populations in NGC 2808 available to date. As typical of the multiple population phenomenon in globular clusters, the light elements vary from one stellar population to another; whereas the iron peak elements show negligible variation between the different populations (at a level of $\lesssim0.10$~dex). Our AGB stars are also characterised by the chemical variations associated with the presence of multiple populations, confirming that this phase of stellar evolution is affected by the phenomenon as well. Intriguingly, we detected one extreme O-poor AGB star (consistent with a high He abundance), challenging stellar evolution models which suggest that highly He-enriched stars should avoid the AGB phase and evolve as AGB-manqu\'e star.

It has been pointed out in arXiv:2211.17057 that our recent (published) paper might be revised, due to an incorrect evaluation of the diffusion coefficients, $D(E)$, employed in the calculations. Unfortunately, there is no indication as to where the incorrectness might be. Here we offer the opportunity to be more specific, by providing the community with the whole description of the equations involved in the calculation of $D(E)$, which is missing in the {\tt arXiv} note. In this context, we mention that no \textit{ad hoc} parameterisation has been used in our paper. Furthermore, any assumption on the injection mechanisms is explicitly described in the paper as an input factor and is obviously part of the modelisation procedure, hence the final outcome is subject to it. Finally, we discuss in a more broad context what, in this calculation of the diffusion coefficient, we believe is the key message.

The purpose of this work is to extend a sample of accurately modeled, benchmark-grade eclipsing binaries with accurately determined masses and radii. We select four "well-behaved" Kepler binaries, KIC2306740, KIC4076952, KIC5193386 and KIC5288543, each with at least 8 double-lined spectra from the apogee instrument that is part of the Sloan Digital Sky Surveys III and IV, and from the Hobby-Eberly High Resolution Spectrograph. We obtain masses and radii with uncertainties of 2.5% or less for all four systems. Three of these systems have orbital periods longer than 9 days, and thus populate an under-sampled region of the parameter space for extremely well-characterized detached eclipsing binaries. We compare the derived masses and radii against MESA MIST isochrones to determine the ages of the systems. All systems were found to be coeval, showing that the results are consistent across MESA MIST and PHOEBE.

The quiet Sun, i.e. in its non-flaring state or non-flaring regions, emits thermal radiation from radio to ultraviolet. The quiet Sun produces also non-thermal radiation observed in gamma rays due to interactions of Galactic Cosmic Rays (GCR) with the solar gas and photons. We report on a new component: the synchrotron emission by GCR electrons in the solar magnetic field. To the best of our knowledge this is the first time this emission has been theoretically claimed and modeled. We find that the measured GCR electrons with energies from tens of GeV to a few TeV produce synchrotron emission in X-rays, which is a few orders of magnitude lower than current upper limits of the quiet Sun set by RHESSI and FOXSI. For a radially decreasing solar magnetic field we find the expected synchrotron intensity to be almost constant in the solar disk, to peak in the close proximity of the Sun, and to quickly drop away from the Sun. We also estimate the synchrotron emission from radio to gamma rays and we compare it with current observations, especially with LOFAR. While it is negligible from radio to UV compared to the solar thermal radiation, this emission can potentially be observed at high energies with NuSTAR and more promising future FOXSI observations. This could potentially allow for constraining CR densities and magnetic-field intensities at the Sun. This study provides a more complete description and a possible new way for understanding the quite Sun and its environment.

A variety of scenarios for early-universe cosmology give rise to a population of primordial black holes (PBHs) with a broad spectrum of masses. In this paper we demonstrate that the evaporation of PBHs in such scenarios has the potential to place the universe into an extended period of "stasis" during which the abundances of matter and radiation remain absolutely constant despite cosmological expansion. This surprising phenomenon can give rise to new possibilities for early-universe dynamics and lead to distinctive signatures of the evaporation of such PBHs.

We present laboratory characterization of kilo-pixel, filled backshort-under-grid (BUG) transition-edge sensor (TES) arrays developed for the Primordial Inflation Polarization ExploreR (PIPER) balloon-borne instrument. PIPER is designed to map the polarization of the CMB on the largest angular scales and characterize dust foregrounds by observing a large fraction of the sky in four frequency bands in the range 200 to 600 GHz. The BUG TES arrays are read out by planar SQUID-based time division multiplexer chips (2dMUX) of matching form factor and hybridized directly with the detector arrays through indium bump bonding. Here, we discuss the performance of the 2dMUX and present measurements of the TES transition temperature, thermal conductance, saturation power, and preliminary noise performance. The detectors achieve saturation power below 1 pW and phonon noise equivalent power (NEP) on the order of a few aW/rtHz. Detector performance is further verified through pre-flight tests in the integrated PIPER receiver, performed in an environment simulating balloon float conditions.

We propose a solution to the problem of accretion disc sizes in active galactic nuclei being larger when measured by reverberation mapping than predicted by theory. Considering blackbody reprocessing on a disc with thickness profile $H(r)$, our solution invokes a steep rim or rippled structures irradiated by the central lamp-post. We model the continuum lags and the faint and bright disc spectral energy distribution (SED) in the best-studied case NGC 5548 (black hole mass $M = 7\times10^{7} M_\odot$, disc inclination $i=45^\circ$). With the lamp-post off, the observed disc SED requires a low accretion rate ($\dot{M} \sim 0.0014 M_\odot$/yr) and high prograde black hole spin ($a \sim 0.93$). Reprocessing on the thin disc gives time lags increasing with wavelength but 3 times smaller than observed. Introducing a steep $H(r)$ rim, or multiple crests, near $r = 5$ light days, reprocessing on their steep centre-facing slopes increases temperatures from $\sim1500$ K to $\sim6000$ K and this increases optical lags to match the lag data. Most of the disc surface maintains the cooler $T\propto r^{-3/4}$ temperature profile that matches the SED. The bright lamp-post may be powered by magnetic links tapping the black hole spin. The steep rim occurs near the sublimation radius for dust in the disc, as in the "failed disc wind model" for broad-line clouds. Lens-Thirring torques aligning the disc and black hole spin may also raise a warp and associated waves. In both scenarios, the small density scale height implied by the inferred value of $H(r)$ suggests possible marginal gravitational instability in the disc.

We present the final results of an imaging and spectroscopic search for stars in the Large Magellanic Cloud (LMC) with C II 7231, 7236 emission lines. The goal is to identify and study [WC11] stars, the coolest of the low-mass Wolf-Rayet sequence, and a subset of central stars of planetary nebulae where the C II lines are known to be especially prominent. A recent serendipitous discovery of an LMC [WC11] raised the possibility that these objects, although difficult to identify, might in fact be more common than previously believed. Several new members of this rare class have been found in this survey. It now seems clear, however, that a significant number of these stars are not hiding amongst the general [WC] population. We point out that the C II doublet intensity ratio observed in our spectra proves to neatly divide the objects into two distinct groups, with the C II emission likely originating from either the stellar wind or a surrounding nebula. The physics of the C II emission mechanism correctly explains this bifurcation. Spectral subtypes are suggested for most of the objects. The numerous spectroscopic clues now available for these objects should facilitate future detailed modeling.

We investigate the energy pathways between the velocity and the magnetic fields in a rotating plane layer dynamo driven by Rayleigh-B\'enard convection using direct numerical simulations. The kinetic and magnetic energies are divided into mean and turbulent components to study the production, transport, and dissipation associated with large and small-scale dynamos. This energy balance-based characterization reveals distinct mechanisms for large- and small-scale magnetic field generation in dynamos, depending on the nature of the velocity field and the conditions imposed at the boundaries.

In the past few decades, various versions of modified gravity theories were proposed to mimic the effect of dark matter. Compared with the conventional Newtonian or relativistic dynamics, these theories contain some extra apparent force terms in the dynamical equations to replace the role of dark matter. Generally speaking, the extra apparent force terms usually scale with radius so that the effect would be significant only on large scale to explain the missing mass in galaxies or galaxy clusters. Nevertheless, the apparent effect may still be observable in small structures like the solar system. In this article, we derive analytic general formulae to represent the contribution of the precession angle of the planets in the solar system due to the general modified gravity theories, in which the extra apparent force terms can be written in a power law of radius $r$ or an exponential function in $r$. We have tested three popular modified gravity theories, the Modified Newtonian Dynamics (MOND), the Emergent Gravity (EG), and the Modified Gravity (MOG). In particular, based on the solar system data, we have constrained the parameters involved for two popular general interpolating functions used in MOND. Our results can be generally applied to both of the modified inertia and modified gravity versions of MOND.

High resolution imaging is achieved using increasingly larger apertures and successively shorter wavelengths. Optical aperture synthesis is an important high-resolution imaging technology used in astronomy. Conventional long baseline amplitude interferometry is susceptible to uncontrollable phase fluctuations, and the technical difficulty increases rapidly as the wavelength decreases. The intensity interferometry inspired by HBT experiment is essentially insensitive to phase fluctuations, but suffers from a narrow spectral bandwidth which results in a lack of detection sensitivity. In this study, we propose optical synthetic aperture imaging based on spatial intensity interferometry. This not only realizes diffraction-limited optical aperture synthesis in a single shot, but also enables imaging with a wide spectral bandwidth. And this method is insensitive to the optical path difference between the sub-apertures. Simulations and experiments present optical aperture synthesis diffraction-limited imaging through spatial intensity interferometry in a 100 $nm$ spectral width of visible light, whose maximum optical path difference between the sub-apertures reach $69.36\lambda$. This technique is expected to provide a solution for optical aperture synthesis over kilometer-long baselines at optical wavelengths.

The Hot Big Bang model predicts the existence of a \emph{cosmic neutrino background}. The number of particles and anti-particles in this primordial bath of neutrinos can be different -- a memory of processes that took place at very early epochs. If neutrinos were massless, this asymmetry would not change once neutrinos froze out. However, in the case of massive particles, the asymmetry is not protected by conservation laws and can get erased via helicity-flipping scatterings off matter inhomogeneities. We evaluate this helicity-flipping rate and demonstrate that if relic lepton asymmetry ever existed, it would remain largely intact in the Earth's neighborhood for realistic values of neutrino masses.

Lattice Field Theory can be used to study finite temperature first-order phase transitions in new, strongly-coupled gauge theories of phenomenological interest. Metastable dynamics arising in proximity of the phase transition can lead to large, uncontrolled numerical errors when analysed with standard methods. In this contribution, we discuss a prototype lattice calculation in which the first-order deconfinement transition in the strong Yang-Mills sector of the standard model is analysed using a novel lattice method, the logarithmic linear relaxation algorithm. This method provides a determination of the density of states of the system with exponential error suppression. Thermodynamic observables can be reconstructed with a controlled error, providing a promising direction for accurate numerical predictions.

We reexamine the regularization of the two-point function of a scalar field in a Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime. Adiabatic regularization provides a set of subtraction terms in momentum space that successfully remove its ultraviolet divergences at coincident points, but can significantly distort the power spectrum at infrared scales, especially for light fields. In this work we propose, by using the intrinsic ambiguities of the renormalization program, a new set of subtraction terms that minimize the distortions for scales $k \lesssim M$, with $M$ an arbitrary mass scale. Our method is consistent with local covariance and equivalent to general regularization methods in curved spacetime. We apply our results to the regularization of the power spectrum in de Sitter space: while the adiabatic scheme yields exactly $\Delta_{\phi}^{\rm (reg)} = 0$ for a massless field, our proposed prescription recovers the standard scale-invariant result $\Delta_{\phi}^{\rm (reg)} \simeq H^2 /(4\pi^2)$ at super-horizon scales.

We review the physics of hyperons and $\Delta$-resonances in the dense matter in compact stars. The covariant density functional approach to the equation of state and composition of dense nuclear matter in the mean-field Hartree and Hartree-Fock approximation is presented, with regimes covering cold $\beta$-equilibrated matter, hot and dense matter with and without neutrinos relevant for the description of supernovas and binary neutron star mergers, as well as dilute expanding nuclear matter in collision experiments. We discuss the static properties of compact stars with hyperons $\Delta$-resonances in light of constraints placed in recent years by the multimessenger astrophysics of compact stars on the compact stars' masses, radii, and tidal deformabilities. The effects of kaon condensation and strong magnetic fields on the composition of hypernuclear stars are also discussed. The properties of rapidly rotating compact hypernuclear stars are discussed and confronted with the observations of 2.5-2.8 solar mass compact objects found in gravitational wave events. We further discuss the cooling of hypernuclear stars, neutrino emission mechanisms hyperonic pairing, and the mass hierarchy in the cooling curves that arise due to the onset of hyperons. The effects of hyperons and $\Delta$-resonances on the equation of state of hot nuclear matter in the dense regime, relevant for the transient astrophysical event and in the dilute regime relevant to the collider physics is discussed. The review closes with a discussion of universal relations among the integral parameters of hot and cold hypernuclear stars and their implications for the analysis of binary neutron star merger events.

Non-singular bouncing cosmologies are well--motivated models for the early universe. Recent observational data are consistent with positive spatial curvature and allow for a natural collapsing and bouncing phase in the very early universe. Additionally, bouncing cosmologies have the potential to rectify conceptual shortcomings identified in the theory of inflation, such as the singularity problem. In this paper we present a classical bouncing model in the context of modified gravity, including an $R^2$-term in the action. We show that after the bounce, the universe enters naturally a period of inflation, driven by the $R^2$--term. We analyse the stability of the model and find that the scalaron assists the stability of the model.

The possible existence of hybrid stars is studied using several multi-quark interaction channels. The hadronic phase consists of an equation of state (EoS) with presently accepted nuclear matter properties and the quark model is constrained by the vacuum properties of several light mesons. The dependence of several NS properties on the different quark interactions is analyzed. We show that the present constraints from neutron star observations allow for the existence of hybrid stars with large strangeness content and large quark cores.

In this paper, we evaluate the energy loss rate of supernovae induced by the axion emission process $\pi^- + p \to n + a$ with the $\Delta(1232)$ resonance in the heavy baryon chiral perturbation theory for the first time. Given the axion-nucleon-$\Delta$ interactions, we include the previously ignored $\Delta$-mediated graphs to $\pi^- + p \to n + a$ process. In particular, the $\Delta^0$-mediated diagram can give a resonance contribution to the supernova axion emission rate when the center-of-mass energy of the pion and proton approaches the $\Delta(1232)$ mass. With these new contributions, we find that for the typical supernova temperatures, compared to the earlier work, the supernova axion emission rate can be enhanced by a factor of $\sim$ 4 in the KSVZ model and up to a factor of $\sim$5 in the DFSZ model with small $\tan\beta$ values. Remarkably, we notice that the $\Delta(1232)$ resonance gives a destructive contribution to the supernova axion emissivity at high supernova temperatures, which is a nontrivial result in this study.

The propagation of perturbations is studied with generalized holonomy corrections in a fully consistent way, ensuring that the deformed algebra of constraints remain closed. The primordial cosmological power spectra are calculated. It is shown that, although the detailed form of the correction does unavoidably impact the observables, the main known results of loop quantum cosmology are robust in this respect.

The lunar poles are unique environments of both great scientific and, increasingly, commercial interest. Consequently, a tension exists between the twin objectives of (a) Exploring the lunar poles for both scientific and commercial purposes and ultimately supporting a lunar economy; and (b) Minimising the environmental impacts on the lunar polar regions so as to preserve them for future scientific investigations. We suggest that the best compromise between these equally valuable objectives would be to restrict scientific and commercial activities to the lunar South Pole, while placing a moratorium on activities at the North Pole until the full consequences of human activities at the South Pole are fully understood and mitigation protocols established. Depending on the pace at which lunar exploration proceeds, such a moratorium might last for several decades in order to properly assess the effects of exploration and commercial activities in regions surrounding the South Pole. A longer term possibility might be to consider designating the lunar North Polar region as a (possibly temporary) Planetary Park. Similar protected status might also be desirable for other unique lunar environments, and, by extension, other scientifically important localities elsewhere in the Solar System.