Papers for Wednesday, Oct 05 2022

The HectoMAP survey provides a complete, mass-limited sample of 30,231 quiescent galaxies with $i-$band Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) imaging that spans the redshift range $0.2 <z < 0.6$. We combine half-light radii based on HSC SSP imaging with redshifts and D$_n4000$ to explore the size - mass relation, $R_{e} = A \times M_{*}^{\alpha}$, and its evolution for the entire HectoMAP quiescent population and for two subsets of the data. Newcomers with $1.5 < \mathrm{D}_n4000 < 1.6$ at each redshift show a steeper increase in $A$ as the universe ages than the population that descends from galaxies that are already quiescent at the survey limit, $z \sim 0.6$ (the resident population). In broad agreement with previous studies, evolution in the size - mass relation both for the entire HectoMAP sample and for the resident population (but not for the newcomers alone) is consistent with minor merger driven growth. For the resident population, the evolution in the size - mass relation is independent of the population age at $z \sim 0.6$. The contrast between the sample of newcomers and the resident population provides insight into the role of commonly termed "progenitor bias" on the evolution of the size - mass relation.

We aim to directly trace the vertical location of the emitting surface of multiple molecular tracers in protoplanetary disks. Our sample of disks includes Elias 2-27, WaOph 6 and the sources targeted by the MAPS ALMA Large Program. The set of molecules studied include CO isotopologues in various transitions, HCN, CN, H2CO, HCO+, C2H and c-C3H2. The vertical emitting region is determined directly from the channel maps, implementing accurate masking of the channel emission to recover the vertical location of the emission surface even at large radial distances from the star and for low-SNR lines. The vertical location of the emitting layer is obtained for 4-10 lines in each disk. IM Lup, HD163296 and MWC 480 12CO and 13CO show vertical modulations, which are coincident with dust gaps and kinematical perturbations. We also present estimates of the gas pressure scale height in the disks from the MAPS sample. Compared to physical-chemical models we find good agreement with the vertical location of CO isotopologues. In HD 163296 CN and HCN trace a similar intermediate layer, for the other disks, the UV flux tracers and the vertical profiles of HCN and C2H are lower than predicted in theoretical models. HCN and H2CO show a highly structured vertical profile, possibly indicative of different formation pathways. It is possible to trace the vertical locations of multiple molecular species that trace a wide variety of physical and chemical disk properties. The distribution of CO isotopologues are found at a wide range of vertical heights $z/r = $ 0.5-0.05. Other molecular lines are mostly found at $z/r \leq $0.15. The vertical layering of molecules is in agreement with theory in some systems, but not in all, therefore dedicated chemical-physical models are needed to further study and understand the emission surfaces.

The cuspy central density profiles of cold dark matter (CDM) haloes make them highly resilient to disruption by tides. Self-interactions between dark matter particles, or the cycling of baryons during galaxy formation, may result in the formation of a constant density core which would make haloes more susceptible to tidal disruption. We use N-body simulations to study the evolution of NFW-like "cored" subhaloes in the tidal field of a massive host, and identify the criteria and timescales for full tidal disruption. Applied to the Milky Way (MW), our results imply that the survival of MW satellites places interesting constraints on core formation. Indeed, we find that no subhaloes with cores larger than 1 per cent of their initial NFW scale radius can survive for a Hubble time on orbits with pericentres <10 kpc. A satellite like Tucana 3, with pericentre ~3.5 kpc, must have a core size smaller than ~2 pc to survive just three orbital periods on its current orbit. The core sizes expected in self-interacting dark matter (SIDM) models with a velocity-independent cross section of 1 cm^2/g seem incompatible with ultra-faint satellites with small pericentric radii, such as Tuc 3, Seg 1, Seg 2, Wil 1, as these would fully disrupt in less than 10 Gyr after infall. These results suggest that many satellites have vanishingly small core sizes, consistent with CDM cusps. The discovery of further Milky Way satellites on orbits with small pericentric radii would strengthen these conclusions and allow for stricter upper limits on the core sizes implied by the survival of Milky Way satellites.

Using 3D radiation-hydrodynamical simulations, we study the effects of ionising radiation on the formation of second-generation (SG) stars in Globular Clusters (GCs) with multiple stellar populations. In particular, we focus on massive ($10^7 \mathrm{M}_{\odot}$) and young (40-Myr old) GCs. We consider stellar winds from asymptotic giant branch (AGB) stars, ram pressure, gas accretion onto the cluster, and photoionisation feedback of binary stars. We find that the stellar luminosity is strong enough to warm and ionise the intracluster medium, but it does not lead to a significant gas expulsion. The cluster can thus retain the ejecta of AGB stars and the accreted pristine gas. In addition, efficient cooling occurs in the central region of the cluster within $50 \mathrm{Myr}$ from the formation of first generation stars, leading to the formation of SG stars. Our results indicate that the inclusion of photoionisation does not suppress SG formation, but rather delays it by about $\sim10 \mathrm{Myr}$. The time delay depends on the density of the pristine gas, so that a denser medium exhibits a shorter delay in star formation. Moreover, photoionisation leads to a modest decrease in the total SG mass, compared to a model without it.

Gyrochronology can yield useful ages for field main-sequence stars, a regime where other techniques are problematic. Typically, gyrochronology relations are calibrated using young ($\lesssim$ 2 Gyr) clusters, but the constraints at older ages are scarce, making them potentially inaccurate and imprecise. In order to test the performance of current gyrochronology relations, we construct samples of stellar pairs with coeval components, for a range of ages and with available rotation periods. These include both randomly paired stars in clusters, and wide binaries in the Kepler field. We design a parameter, $\Delta P_{rot, gyro}$, that quantifies the level of agreement between the components of coeval stellar pairs and the gyrochronology calibrations under scrutiny. Our results show that wide binaries and cluster members are in better concordance with gyrochronology than samples of randomly paired field stars (which we use as a control group), confirming that the relations do have predicting power in identifying coeval stars. However, the agreement with the examined relations decreases for older stars, suggesting a degradation of the current gyrochronology with age, in agreement with recent literature results. We argue for the need of novel gyrochronology constraints at older ages that may allow revised calibrations. Ultimately, using coeval stars to test gyrochronology poses the advantage of circumventing the need for age determinations while simultaneously exploiting larger samples at older ages. In addition, taking gyrochronology at face value, we note that our results provide new empirical evidence that the components of field wide binaries are indeed coeval.

We present the first comprehensive release of photometric redshifts (photo-z's) from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) team. We use statistics based upon the Quantile-Quantile (Q--Q) plot to identify biases and signatures of underestimated or overestimated errors in photo-z probability density functions (PDFs) produced by six groups in the collaboration; correcting for these effects makes the resulting PDFs better match the statistical definition of a PDF. After correcting each group's PDF, we explore three methods of combining the different groups' PDFs for a given object into a consensus curve. Two of these methods are based on identifying the minimum f-divergence curve, i.e., the PDF that is closest in aggregate to the other PDFs in a set (analogous to the median of an array of numbers). We demonstrate that these techniques yield improved results using sets of spectroscopic redshifts independent of those used to optimize PDF modifications. The best photo-z PDFs and point estimates are achieved with the minimum f-divergence using the best 4 PDFs for each object (mFDa4) and the Hierarchical Bayesian (HB4) methods, respectively. The HB4 photo-z point estimates produced $\sigma_{\rm NMAD} = 0.0227/0.0189$ and $|\Delta z/(1+z)| > 0.15$ outlier fraction = 0.067/0.019 for spectroscopic and 3D-HST redshifts, respectively. Finally, we describe the structure and provide guidance for the use of the CANDELS photo-z catalogs, which are available at https://archive.stsci.edu/hlsp/candels.

We present the analysis of existing optical photometry and new optical spectroscopy of the candidate cataclysmic variable star OGLE-BLG504.12.201843. As was shown previously, this object has an orbital period of 0.523419 days and exhibits year-long outbursts with a mean period of 973 days. Using digitized photographic archives, we show that the earliest recorded outburst occurred in 1910. We propose that this object is a U Gem-type dwarf nova with extreme properties. The orbital variability of the system in outburst shows clear signs of an accretion disc, from which the outburst likely originates. During quiescence, the object slowly brightens by up to $0.75$ mag in the $I$ band over 600 days before the outburst and exhibits small flares with amplitude $\lesssim 0.2$ mag in the $I$ band. We interpret the gradual brightening as an increase in the luminosity and temperature of the accretion disc, which is theoretically predicted but only rarely seen in DNe. The origin of small flares remains unexplained. The spectra shows Balmer absorption lines both in quiescence and outburst, which can be associated with a bright secondary star or a cold accretion disc. During outbursts, emission lines with FWHM of about 450 km s$^{-1}$ appear, but they lack typical double-peaked profiles. We suggest that either these lines originate in the disc winds or the orbital inclination is low, the latter being consistent with constrains obtained from the orbital variability of the system. Due to its extreme properties and peculiarities, OGLE-BLG504.12.201843 is an excellent object for further follow-up studies.

We perform 3D hydrodynamics simulations of disc-disc stellar flybys with on-the-fly Monte Carlo radiative transfer. We show that pre-existing circumstellar discs around both stars result in fast rising ($\sim$yrs) outbursts lasting 2-5 times longer than for a star-disc flyby. The perturber always goes into outburst ($\dot{M}>10^{-5}~{\rm M_{\odot}~ yr^{-1}}$). Whereas we find that the primary goes into a decades long outburst only when the flyby is retrograde to the circumprimary disc rotation. High accretion rates during the outburst are triggered by angular momentum cancellation in misaligned material generated by the encounter. A large fraction of accreted material is alien.

The nature of the newly discovered fast blue optical transients (FBOTs) is still puzzling astronomers. In this paper we carry out a comprehensive analysis of the molecular gas, ionized gas and stellar populations in the environment of the nearby FBOT AT2018cow based on ALMA, VLT/MUSE and HST/WFC3 observations. A prominent molecular concentration of 6 ($\pm$ 1) $\times$ 10$^6$ $M_\odot$ is found in the vicinity of AT2018cow, which has given rise to two active star-forming complexes with ages of 4 $\pm$ 1 Myr and $\lesssim$2.5 Myr, respectively. Each star-forming complex has a stellar mass of 3 $\times$ 10$^5$ $M_\odot$ and has photoionized a giant H II region with H$\alpha$ luminosity even comparable to that of the 30 Dor mini-starburst region. AT2018cow is spatially coincident with one of the star-forming complexes; however, it is most likely to reside in its foreground since it has a much smaller extinction than the complex. Its progenitor could have been formed at an earlier epoch in this area; if it were from a major star-forming event, the non-detection of the associated stellar population constrains the progenitor's age to be $\gtrsim$10 Myr and initial mass to be $\lesssim$ 20 $M_\odot$. We further find the late-time brightness of AT2018cow is unlikely to be a stellar object. Its brightness has slightly declined from 2 yr to 4 yr after explosion and is most likely to originate from AT2018cow itself due to some powering mechanism still working at such late times.

Misalignments between the rotation axis of stars and gas are an indication of external processes shaping galaxies throughout their evolution. Using observations of 3068 galaxies from the SAMI Galaxy Survey, we compute global kinematic position angles for 1445 objects with reliable kinematics and identify 169 (12%) galaxies which show stellar-gas misalignments. Kinematically decoupled features are more prevalent in early-type/passive galaxies compared to late-type/star-forming systems. Star formation is the main source of gas ionisation in only 22% of misaligned galaxies; 17% are Seyfert objects, while 61% show Low-Ionisation Nuclear Emission-line Region features. We identify the most probable physical cause of the kinematic decoupling and find that, while accretion-driven cases are dominant, for up to 8% of our sample, the misalignment may be tracing outflowing gas. When considering only misalignments driven by accretion, the acquired gas is feeding active star formation in only $\sim$1/4 of cases. As a population, misaligned galaxies have higher S\'ersic indices and lower stellar spin & specific star formation rates than appropriately matched samples of aligned systems. These results suggest that both morphology and star formation/gas content are significantly correlated with the prevalence and timescales of misalignments. Specifically, torques on misaligned gas discs are smaller for more centrally concentrated galaxies, while the newly accreted gas feels lower viscous drag forces in more gas-poor objects. Marginal evidence of star formation not being correlated with misalignment likelihood for late-type galaxies suggests that such morphologies in the nearby Universe might be the result of preferentially aligned accretion at higher redshifts.

The near-ultraviolet (NUV) spectral region is a useful diagnostic for stellar flare physics and assessing the energy environment of young exoplanets, especially as relates to prebiotic chemistry. We conducted a pilot NUV spectroscopic flare survey of the young M dwarf AU Mic with the Neil Gehrels Swift Observatory's UltraViolet and Optical Telescope. We detected four flares and three other epochs of significantly elevated count rates during the 9.6 hours of total exposure time, consistent with a NUV flare rate of $\sim$0.5 hour$^{-1}$. The largest flare we observed released a minimum energy of 6$\times$10$^{33}$ erg between 1730-5000 \r{A}. All flares had durations longer than the $\sim$14-17 minute duration of each Swift visit, making measuring total flare energy and duration infeasible.

We present an updated model for the extragalactic background light (EBL) from stars and dust, over wavelengths approximately 0.1 to 1000 $\mu$m. This model uses accurate theoretical stellar spectra, and tracks the evolution of star formation, stellar mass density, metallicity, and interstellar dust extinction and emission in the universe with redshift. Dust emission components are treated self-consistently, with stellar light absorbed by dust reradiated in the infrared as three blackbody components. We fit our model, with free parameters associated with star formation rate and dust extinction and emission, to a wide variety of data: luminosity density, stellar mass density, and dust extinction data from galaxy surveys; and $\gamma$-ray absorption optical depth data from $\gamma$-ray telescopes. Our results strongly constraint the star formation rate density and dust photon escape fraction of the universe out to redshift $z=10$, about 90% of the history of the universe. We find our model result is, in some cases, below lower limits on the $z=0$ EBL intensity, and below some low-$z$ $\gamma$-ray absorption measurements.

This work investigates the evolution of the distribution of charged particles (cosmic rays) due to the mechanism of stochastic turbulent acceleration (STA) in presence of small-scale turbulence with a mean magnetic field. STA is usually modelled as a biased random walk process in the momentum space of the non-thermal particles. This results in an advection-diffusion type transport equation for the non-thermal particle distribution function. Under quasilinear approximation, and by assuming a turbulent spectrum with single scale injection at sub-gyroscale, we find that the Fokker-Planck diffusion coefficients $D_{\gamma\gamma}$ and $D_{\mu\mu}$ scale with the Lorentz factor $\gamma$ as: $D_{\gamma\gamma}\propto\gamma^{-2/3}$ and $D_{\mu\mu}\propto\gamma^{-8/3}$. Furthermore, with the calculated transport coefficients, we numerically solve the advection-diffusion type transport equation for the non-thermal particles. We demonstrate the interplay of various mircophysical processes such as STA, synchrotron loss and particle escape on the particle distribution by systematically varying the parameters of the problem.

Relativistic plasma flows from the jets of black hole binary systems consist the environment of multiple particle production and radiation emission including neutrinos and gamma-rays. We implement a hadronic model based on $p-p$ interactions with the purpose of predicting the produced secondary particle distributions inside the jet. Our ultimate goal is the neutrino and gamma-ray intensities calculation while taking into account the most important gamma-ray absorption processes in order to present more realistic results.

From 2018 to 2020, the Scientific Balloon Roadmap Program Analysis Group (Balloon Roadmap PAG) served as an community-based, interdisciplinary forum for soliciting and coordinating community analysis and input in support of the NASA Scientific Balloon Program. The Balloon Roadmap PAG was tasked with articulating and prioritizing the key science drivers and needed capabilities of the Balloon Program for the next decade. Additionally, the Balloon Roadmap PAG was asked to evaluate the potential for achieving science goals and maturing technologies of the Science Mission Directorate, evaluate the Balloon Program goals towards community outreach, and asses commercial balloon launch opportunities. The culmination of this work has been a written report submitted to the NASA Astrophysics Division Director.

Halo abundance and structures are critical to understand the small scale challenges for $\Lambda$CDM cosmology. We present a unified theory and analytical models for both halo mass function and halo density profile based on a random walk of halos and dark matter particles. The position dependent waiting time in random walk leads to a stretched Gaussian for mass function and density with a power-law on small scale and an exponential decay on large scale. Both Press-Schechter mass function and Einasto density profile can be easily recovered. This new perspective provides a simple theory for universal halo mass function and density profile.

We investigate general relativistic magnetohydrodynamic simulations (GRMHD) to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterising the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnetic field structure in the emitting region. Moreover, the dependence of this twisty pattern on spin can be attributed to changes in the magnetic field geometry that occur due to the frame dragging. By studying an analytic ring model, we find that the roles of Doppler boosting and lensing are subdominant. Faraday rotation may cause a systematic shift in the linear polarization pattern, but we find that its impact is subdominant for models with strong magnetic fields and modest ion-to-electron temperature ratios. Models with weaker magnetic fields are much more strongly affected by Faraday rotation and have more complicated emission geometries than can be captured by a ring model. However, these models are currently disfavoured by the recent EHT observations of M87*. Our results suggest that linear polarization maps can provide a probe of the underlying magnetic field structure around a black hole, which may then be usable to indirectly infer black hole spins. The generality of these results should be tested with alternative codes, initial conditions, and plasma physics prescriptions.

The observation of small bodies in the Space Environment is an ongoing important task in astronomy. While nowadays new objects are mostly detected in larger sky surveys, several follow-up observations are usually needed for each object to improve the accuracy of orbit determination. In particular objects orbiting close to Earth, so called Near-Earth Objects are of special concern as a small but not negligible fraction of them can have a non-zero impact probability with Earth. Telescopes are often hosted by amateur observatories. With upcoming new NEO search campaigns by very wide field of view telescopes, like the Vera C. Rubin Observatory, NASA's NEO surveyor space mission and ESA's Flyeye telescopes, the number of NEO discoveries will increase dramatically. This will require an increasing number of useful telescopes for follow-up observations at different geographical locations. While well-equipped amateur astronomers often host instruments which might be capable of creating useful measurements, both observation planning and scheduling, and also analysis are still a major challenge for many observers. In this work we present a fully robotic planning, scheduling and observation pipeline that extends the widely used open-source cross-platform software KStars/Ekos for INDI devices. The method consists of algorithms which automatically select NEO candidates with priority according to ESA's NEOCC. It then analyses detectable objects (based on limiting magnitudes, geographical position, and time) with preliminary ephemeris from the Minor Planet Center. Optimal observing slots during the night are calculated and scheduled. Immediately before the measurement the accurate position of the minor body is recalculated and finally the images are taken. Besides the detailed description of all components, we will show a complete robotic hard- and software solution based on our methods.

We have collected a new deep {\it Chandra X-ray Observatory} ({\it CXO}) exposure of PSR B2224+65 and the `Guitar Nebula', mapping the complex X-ray structure. This is accompanied by a new {\it HST} H$\alpha$ image of the head of the Guitar. Comparing the {\it HST} and {\it CXO} structures in 4 epochs over 25 years, we constrain the evolution of the TeV particles that light up the filament. Cross-field diffusion appears to be enhanced, likely by the injected particles, behind the filament's sharp leading edge, explaining the filament width and its evolving surface brightness profile.

The purpose of this white paper is to put together a coherent vision for the role of helioseismic monitoring of magnetic activity in the Sun's far hemisphere that will contribute to improving space weather forecasting as well as fundamental research in the coming decade. Our goal fits into the broader context of helioseismology in solar research for any number of endeavors when helioseismic monitors may be the sole synoptic view of the Sun's far hemisphere. It is intended to foster a growing understanding of solar activity, as realistically monitored in both hemispheres, and its relationship to all known aspects of the near-Earth and terrestrial environment. Some of the questions and goals that can be fruitfully pursued through seismic monitoring of farside solar activity in the coming decade include: What is the relationship between helioseismic signatures and their associated magnetic configurations, and how is this relationship connected to the solar EUV irradiance over the period of a solar rotation?; How can helioseismic monitoring contribute to data-driven global magnetic-field models for precise space weather forecasting?; What can helioseismic monitors tell us about prospects of a flare, CME or high-speed stream that impacts the terrestrial environment over the period of a solar rotation?; How does the inclusion of farside information contribute to forecasts of interplanetary space weather and the environments to be encountered by human crews in interplanetary space? Thus, it is crucial for the development of farside monitoring of the Sun be continued into the next decade either through ground-based or space-borne observations.

Cosmic dust plays a dominant role in the universe, especially in the formation of stars and planetary systems. Furthermore, the surface of cosmic dust grains is the bench-work where molecular hydrogen and simple organic compounds are formed. We manipulate individual dust particles in water solution by contactless and non-invasive techniques such as standard and Raman tweezers, to characterize their response to mechanical effects of light (optical forces and torques) and to determine their mineral compositions. Moreover, we show accurate optical force calculations in the T-matrix formalism highlighting the key role of composition and complex morphology in optical trapping of cosmic dust particles.This opens perspectives for future applications of optical tweezers in curation facilities for sample return missions or in extraterrestrial environments.

With respect to the $\rm H\beta$ full width at half-maximum ($\rm FWHM_{H\beta}$), the broad $\rm H\beta$ line dispersion ($\sigma_{\rm H\beta}$) was preferred as a velocity tracer to calculate the single-epoch supermassive black hole mass ($M_{\rm BH}$) suggested by \cite{Yu2020b}. For a compiled sample of 311 broad-line active galactic nuclei (AGN) with measured hard X-ray photon index ($z<0.7$), $\sigma_{\rm H\beta}$ and the optical Fe II relative strength ($R_{\rm Fe}$) are measured from their optical spectra, which are used to calculate $\sigma_{\rm H\beta}$-based virial $M_{\rm BH}$ and dimensionless accretion rate ($\dot{\mathscr{M}}$). With respect to $\rm FWHM_{\rm H\beta}$, it is found that the mean value of $\sigma_{\rm H\beta}$-based $M_{\rm BH}$ is on average larger by 0.26 dex, and the mean value of $\sigma_{\rm H\beta}$-based $\dot{\mathscr{M}}$ is on average smaller by 0.51 dex. It is found that there exists a non-linear relationship between the Eddington ratio ($L_{\rm Bol}/L_{\rm Edd}$) and $\dot{\mathscr{M}}$, i.e., $L_{\rm Bol}/L_{\rm Edd} \propto \dot{\mathscr{M}}^{0.56\pm 0.01}$. This non-linear relationship comes from the accretion efficiency $\eta$, which is smaller for AGN with higher $\dot{\mathscr{M}}$. We find a strong bivariate correlation of the fraction of energy released in the corona $F_{\rm X}$ with $\dot{\mathscr{M}}$ and \mbh, $F_{\rm X} \propto \dot{\mathscr{M}}^{-0.57\pm 0.05} M_{\rm BH}^{-0.54\pm 0.06}$. The flat slope of $-0.57\pm 0.05$ favours the shear stress tensor of the accretion disk being proportional to the geometric mean of gas pressure and total pressure. We find a strong bivariate relation of $\Gamma$ with $\dot{\mathscr{M}}$ and $F_{\rm X}$, $\Gamma \propto \dot{\mathscr{M}}^{-0.21\pm 0.02}F_{\rm X}^{0.02\pm 0.04}$. The hard X-ray spectrum becomes softer with increasing of $F_{\rm X}$, although the scatter is large.

The shape of emission lines in the optical spectra of star-forming galaxies reveals the kinematics of the diffuse gaseous component. We analyse the shape of prominent emission lines in a sample of ~53,000 star-forming galaxies from the Sloan Digital Sky Survey, focusing on departures from gaussianity. Departures from a single gaussian profile allow us to probe the motion of gas and to assess the role of outflows. The sample is divided into groups according to their stellar velocity dispersion and star formation rate. The spectra within each group are stacked to improve the signal-to-noise ratio of the emission lines, to remove individual signatures, and to enhance the effect of star formation rate on the shapes of the emission lines. The moments of the emission lines, including kurtosis and skewness, are determined. We find that most of the emission lines in strong star-forming systems unequivocally feature negative kurtosis. This signature is present in H$\beta$, H$\alpha$, [N II] and [S II] in massive galaxies with high star formation rates. We attribute it as evidence of radial outflows of ionised gas driven by the star formation of the galaxies. Also, most of the emission lines in low-mass systems with high star formation rates feature negative skewness, and we interpret it as evidence of dust obscuration in the galactic disk. These signatures are however absent in the [O III] line, which is believed to trace a different gas component. The observed trend is significantly stronger in face-on galaxies, indicating that star formation drives the outflows along the galactic rotation axis, presumably the path of least resistance. The data suggest that outflows driven by star formation exert accumulated impacts on the interstellar medium, and the outflow signature is more evident in older galaxies as they have experienced a longer total duration of star formation.

We report the Fermi Large Area Telescope (Fermi-LAT) detection of the $\gamma$-ray emission toward the massive star forming region of Carina Nebula Complex (CNC). Using the latest source catalog and diffuse background models, we found that the GeV $\gamma$-ray emission in this region can be resolved into three different components. The GeV $\gamma$-ray emission from the central point source is considered to originate from the Eta Carina ($\eta$ Car). We further found the diffuse GeV $\gamma$-ray emission around the CNC which can be modelled by two Gaussian disks with radii of 0.4\deg (region A) and 0.75\deg (region B), respectively. The GeV $\gamma$-ray emission from both the regions A and B have good spatial consistency with the derived molecular gas in projection on the sky. The GeV $\gamma$-ray emission of region A reveals a characteristic spectral shape of the pion-decay process, which indicates that the $\gamma$-rays are produced by the interactions of hadronic cosmic rays with ambient gas. The $\gamma$-rays spectrum of region B has a hard photon index of 2.12 $\pm$ 0.02, which is similar to other young massive star clusters. We argue that the diffuse GeV $\gamma$-ray emission in region A and region B likely originate from the interaction of accelerated protons in clusters with the ambient gas.

Optical SETI (Search for Extraterrestrial Intelligence) instruments that can explore the very fast time domain, especially with large sky coverage, offer an opportunity for new discoveries that can complement multimessenger and time domain astrophysics. The Panoramic SETI experiment (PANOSETI) aims to observe optical transients with nanosecond to second duration over a wide field-of-view ($\thicksim$2,500 sq.deg.) by using two assemblies of tens of telescopes to reject spurious signals by coincidence detection. Three PANOSETI telescopes, connected to a White Rabbit timing network used to synchronize clocks at the nanosecond level, have been deployed at Lick Observatory on two sites separated by a distance of 677 meters to distinguish nearby light sources (such as Cherenkov light from particle showers in the Earth's atmosphere) from astrophysical sources at large distances. In parallel to this deployment, we present results obtained during four nights of simultaneous observations with the four 12-meter VERITAS gamma-ray telescopes and two PANOSETI telescopes at the Fred Lawrence Whipple Observatory. We report PANOSETI's first detection of astrophysical gamma rays, comprising three events with energies in the range between $\thicksim$15 TeV and $\thicksim$50 TeV. These were emitted by the Crab Nebula, and identified as gamma rays using joint VERITAS observations.

Titan has a diverse range of materials in its atmosphere and on its surface: the simple organics that reside in various phases (gas, liquid, ice) and the solid complex refractory organics that form Titan's haze layers. These materials all actively participate in various physical processes on Titan, and many material properties are found to be important in shaping these processes. Future in-situ exploration on Titan would likely encounter a range of materials, and a comprehensive database to archive the material properties of all possible material candidates will be needed. Here we summarize several important material properties of the organic liquids, ices, and the refractory hazes on Titan that are available in the literature and/or that we have computed. These properties include thermodynamic properties (phase change points, sublimation and vaporization saturation vapor pressure, and latent heat), physical property (density), and surface properties (liquid surface tensions and solid surface energies). We have developed a new database to provide a repository for these data and make them available to the science community. These data can be used as inputs for various theoretical models to interpret current and future remote sensing and in-situ atmospheric and surface measurements on Titan. The material properties of the simple organics may also be applicable for giant planets and icy bodies in the outer solar system, interstellar medium, and protoplanetary disks.

Although a large number of Galactic open clusters (OCs) have been identified, the internal kinematic properties (e.g., rotation) of almost all the known OCs are still far from clear. With the high-precision astrometric data of Gaia EDR3, we have developed a methodology to unveil the rotational properties of the Praesepe cluster. Statistics of the three-dimensional residual motions of the member stars reveal the presence of Praesepe's rotation and determine its spatial rotation axis. The mean rotation velocity of the Praesepe cluster within its tidal radius is estimated to be 0.2 $\pm$ 0.05 km s$^{-1}$, and the corresponding rotation axis is tilted in relation to the Galactic plane with an angle of 41 $\pm$ 12 degree. We also analysed the rms rotational velocity of the member stars around the rotation axis, and found that the rotation of the member stars within the tidal radius of Praesepe probably follows the Newton's classical theorems.

Wide-field time domain facilities detect transient events in large numbers through difference imaging. For example, Zwicky Transient Facility produces alerts for hundreds of thousands of transient events per night, a rate set to be dwarfed by the upcoming Vera Rubin Observatory. The automation provided by Machine Learning (ML) is, therefore, necessary to classify these events and select the most interesting sources for follow-up observations. Cataclysmic Variables (CVs) are a transient class that are numerous, bright, and nearby, providing excellent laboratories for the study of accretion and binary evolution. Here we focus on our use of ML to identify CVs from photometric data of transient sources published by the Gaia Science Alerts program (GSA) - a large, easily accessible resource, not fully explored with ML. The use of light curve feature extraction techniques and source metadata from the Gaia survey resulted in a Random Forest model capable of distinguishing CVs from supernovae, Active Galactic Nuclei, and Young Stellar Objects with a 92\% precision score and an 85\% hit rate. Of 13,280 sources within GSA without an assigned transient classification our model predicts the CV class for $\sim$2800. Spectroscopic observations are underway to classify a statistically significant sample of these targets to validate the performance of the model. This work puts us on a path towards the classification of rare CV subtypes from future wide-field surveys such as the Legacy Survey of Space and Time.

We present a dynamical analysis of the fragmented active asteroid 331P/Gibbs. Using archival images taken by the Hubble Space Telescope from 2015 to 2018, we measured the astrometry of the primary and the three brightest (presumably the largest) components. Conventional orbit determination revealed a high-degree of orbital similarity between the components. We then applied a fragmentation model to fit the astrometry, obtaining key parameters including the fragmentation epochs and separation velocities. Our best-fit models show that Fragment B separated from the primary body at a speed of $\sim$1 cm s$^{-1}$ between 2011 April and May, whereas two plausible scenarios were identified for Fragments A and C. The former split either from the primary or from Fragment B, in 2011 mid-June at a speed of $\sim$8 cm s$^{-1}$, and the latter split from Fragment B either in late 2011 or between late 2013 and early 2014, at a speed of $\sim$0.7-0.8 cm s$^{-1}$. The results are consistent with rotational disruption as the mechanism causing the cascading fragmentation of the asteroid, as suggested by the rapid rotation of the primary. The fragments constitute the youngest known asteroid cluster, providing us with a great opportunity to study asteroid fragmentation and formation of asteroid clusters.

The Atacama Large Millimeter/submillimeter Array with the planned electronic upgrades will deliver an unprecedented amount of deep and high resolution observations. Wider fields of view are possible with the consequential cost of image reconstruction. Alternatives to commonly used applications in image processing have to be sought and tested. Advanced image reconstruction methods are critical to meet the data requirements needed for operational purposes. Astrostatistics and astroinformatics techniques are employed. Evidence is given that these interdisciplinary fields of study applied to synthesis imaging meet the Big Data challenges and have the potentials to enable new scientific discoveries in radio astronomy and astrophysics.

The probability density function of accretion disc luminosity fluctuations at high observed energies (i.e., energies larger than the peak temperature scale of the disc) is derived, under the assumption that the temperature fluctuations are log-normally distributed. Thin disc theory is used throughout. While log-normal temperature fluctuations would imply that the disc's bolometric luminosity is also log-normal, the observed Wien-like luminosity behaves very differently. For example, in contrast to a log-normal distribution, the standard deviation of the derived distribution is not linearly proportional to its mean. This means that these systems do not follow a linear rms-flux relationship. Instead they exhibit very high intrinsic variance, and undergo what amounts to a phase transition, in which the mode of the distribution (in the statistical sense) ceases to exist, even for physically reasonable values of the underlying temperature variance. The moments of this distribution are derived using asymptotic expansion techniques. A result that is important for interpreting observations is that the theory predicts that the fractional variability of these disc systems should increase as the observed frequency is increased. The derived distribution will be of practical utility in quantitatively understanding the variability of disc systems observed at energies above their peak temperature scale, including X-ray observations of tidal disruption events.

We probe how common extremely rapid rotation is among massive stars in the early universe by measuring the OBe star fraction in nearby metal-poor dwarf galaxies. We apply a new method that uses broad-band photometry to measure the galaxy-wide OBe star fractions in the Magellanic Clouds and three more distant, more metal-poor dwarf galaxies. We find OBe star fractions of ~20% in the Large Magellanic Cloud (0.5 Z_Solar), and ~30% in the Small Magellanic Cloud (0.2 Z_Solar) as well as in the so-far unexplored metallicity range from 0.1 Z_solar to 0.2 Z_solar occupied by the other three dwarf galaxies. Our results imply that extremely rapid rotation is common among massive stars in metal-poor environments such as the early universe.

POLAR-2 is a space-borne polarimeter, built to investigate the polarization of Gamma-Ray Bursts and help elucidate their mechanisms. The instrument is targeted for launch in 2024 or 2025 aboard the China Space Station and is being developed by a collaboration between institutes from Switzerland, Germany, Poland and China. The instrument will orbit at altitudes between 340km and 450km with an inclination of 42$^{\circ}$ and will be subjected to background radiation from cosmic rays and solar events. It is therefore pertinent to better understand the performance of sensitive devices under space-like conditions. In this paper we focus on the radiation damage of the silicon photomultiplier arrays S13361-6075NE-04 and S14161-6050HS-04 from Hamamatsu. The S13361 are irradiated with 58MeV protons at several doses up to 4.96Gy, whereas the newer series S14161 are irradiated at doses of 0.254Gy and 2.31Gy. Their respective performance degradation due to radiation damage are discussed. The equivalent exposure time in space for silicon photomultipliers inside POLAR-2 with a dose of 4.96Gy is 62.9 years (or 1.78 years when disregarding the shielding from the instrument). Primary characteristics of the I-V curves are an increase in the dark current and dark counts, mostly through cross-talk events. Annealing processes at $25^{\circ}C$ were observed but not studied in further detail. Biasing channels while being irradiated have not resulted in any significant impact. Activation analyses showed a dominant contribution of $\beta^{+}$ particles around 511keV. These resulted primarily from copper and carbon, mostly with decay times shorter than the orbital period.

Atmospheres play a crucial role in planetary habitability. Around M dwarfs and young Sun-like stars, planets receiving the same insolation as the present-day Earth are exposed to intense stellar X-rays and extreme-ultraviolet (XUV) radiation. This study explores the fundamental question of whether the atmosphere of present-day Earth could survive in such harsh XUV environments. Previous theoretical studies suggest that stellar XUV irradiation is sufficiently intense to remove such atmospheres completely on short timescales. In this study, we develop a new upper-atmospheric model and re-examine the thermal and hydrodynamic responses of the thermospheric structure of an Earth-like N2-O2 atmosphere, on an Earth-mass planet, to an increase in the XUV irradiation. Our model includes the effects of radiative cooling via electronic transitions of atoms and ions, known as atomic line cooling, in addition to the processes accounted for by previous models. We demonstrate that atomic line cooling dominates over the hydrodynamic effect at XUV irradiation levels greater than several times the present level of the Earth. Consequentially, the atmosphere's structure is kept almost hydrostatic, and its escape remains sluggish even at XUV irradiation levels up to a thousand times that of the Earth at present. Our estimates for the Jeans escape rates of N2-O2 atmospheres suggest that these 1 bar atmospheres survive in early active phases of Sun-like stars. Even around active late M dwarfs, N2-O2 atmospheres could escape significant thermal loss on timescales of gigayears. These results give new insights into the habitability of terrestrial exoplanets and the Earth's climate history.

The post-common envelope binary V471 Tauri has been an object of interest for decades. V471 Tau shows various phenomena due to its evolutionary state and unique properties, e.g. its magnetic accretion and eclipse timing variation (ETV). Previous authors explained the ETVs by different, sometimes contradictory theories. In this paper, we present and analyse the variability of the eclipse timing of this star. We observed V471 Tauri over the last ten years and covered the second cycle of its period variation. Based on our analysis of the presented data, we assess the possible existence of a brown dwarf in this system and derive its orbital parameters. We compare the results of our dynamical modelling to the solution predicted by Applegate-mechanism theories, which have been developed in recent studies. We found that the observed ETV cannot be explained only by the presence of additional components to the binary.

Understanding the radiative-dynamical coupling between upper photosphere and deeper atmosphere is a key in understanding the abnormal large radii of hot Jupiters. One needs very long integration times of 3D general circulation models (GCMs) with self consistent radiative transfer to achieve a better understanding of the feedback process between dynamics and radiation. We here present the longest 3D non-gray GCM study (86000 d) of an ultra hot Jupiter (WASP-76 b) published to this date that reached a final converged state. Furthermore, we present a method that can be used to accelerate the path towards temperature convergence in the deep atmospheric layers. We find that the final converged temperature profile is cold in the deep atmospheric layers, lacking any sign of vertical transport of potential temperature by large scale atmospheric motions. We thus conclude that the coupling between radiation and dynamics alone is not sufficient to explain the abnormal large radii of inflated hot gas giants.

Observations and simulations of barred spiral galaxies have shown that, in general, the spiral arms rotate at a different pattern speed to that of the bar. The main conclusion from the bibliography is that the bar rotates faster than the spiral arms with a double or even a triple value of angular velocity. The theory that prevails in explaining the formation of the spiral arms in the case of a barred spiral galaxy with two pattern speeds is the manifold theory, where the orbits that support the spiral density wave are chaotic, and are related to the manifolds emanating from the Lagrangian points L_1 and L_2 at the end of the bar. In the present study, we consider an alternative scenario in the case where the bar rotates fast enough in comparison with the spiral arms and the bar potential can be considered as a perturbation of the spiral potential. In this case, the stable elliptical orbits that support the spiral density wave (in the case of grand design galaxies) are transformed into quasiperiodic orbits (or 2D tori) with a certain thickness. The superposition of these perturbed preccesing ellipses for all the energy levels of the Hamiltonian creates a slightly perturbed symmetrical spiral density wave.

Forming massive stars launch outflows of magnetic origin, which in fact serve as a marker for finding sites of massive star formation. However, both the theoretical and observational study of the mechanisms that intervene in the formation and propagation of such outflows has been possible only until recent years. With this work, we aim to study in detail the mechanisms that drive highly collimated outflows from early stages of the formation of a massive star, and how those processes are impacted by the properties of the natal environment of the forming massive star. We perform a series of 31 simulations with the aim of building a unified theoretical picture of these mechanisms, and see how the impact of different environments alter their morphology and momentum output. The magnetohydrodynamical simulations consider also Ohmic dissipation as a nonideal effect, self-gravity, and diffusive radiation transport for thermal absorption and emission by the dust and gas. We start from a collapsing cloud core that is threaded by an initially-uniform magnetic field and which is slowly rotating. We utilize a two-dimensional axisymmetric grid in spherical coordinates. In the simulations, we can clearly distinguish a fast, magneto-centrifugally launched and collimated jet (of speeds > 100 km/s), from a wider magnetic tower flow driven by magnetic pressure which broadens in time. We analyze in detail the acceleration of the flow, and its re-collimation by magnetic forces happening at distances of several hundreds of astronomical units. We quantify the impact of magnetic braking in the outflows, which narrows the outflow cavity for the late evolution of the system. We observe the presence of the same jet-driving mechanisms for a wide range of assumptions on the natal environment of the massive protostar, but with changes to their morphology and mechanical feedback into larger scales over time.

The latest Gaia Data Release 3 provides an opportunity to expand the census of Galactic open clusters harboring classical Cepheid variables, thereby bolstering the cosmic distance scale. A comprehensive analysis yielded a total of 50 classical Cepheids associated with 45 open clusters, of which 39 open cluster-classical Cepheid pairs are considered probable, with the remaining 11 pairs considered improbable but worth following up. Two previously identified clusters by us possibly host classical Cepheids (OC-0125/V1788 Cyg and OC-0675/OGLE-BLG-CEP-114). In addition, we identify 38 new open cluster candidates within the Galactic disk.

We present the first survey of quiet Sun features observed in hard X-rays (HXRs), using the the Nuclear Spectroscopic Telescope ARray (NuSTAR), a HXR focusing optics telescope. The recent solar minimum combined with NuSTAR's high sensitivity has presented a unique opportunity to perform the first HXR imaging spectroscopy on a range of features in the quiet Sun. By studying the HXR emission of these features we can detect or constrain the presence of high temperature (>5 MK) or non-thermal sources, to help understand how they relate to larger more energetic solar phenomena, and determine their contribution to heating the solar atmosphere. We report on several features observed in the 28 September 2018 NuSTAR full-disk quiet Sun mosaics, the first of the NuSTAR quiet Sun observing campaigns, which mostly include steady features of X-ray bright points and an emerging flux region which later evolved into an active region, as well as a short-lived jet. We find that the features' HXR spectra are well fitted with isothermal models with temperatures ranging between 2.0-3.2 MK. Combining the NuSTAR data with softer X-ray emission from Hinode/XRT and EUV from SDO/AIA we recover the differential emission measures, confirming little significant emission above 4 MK. The NuSTAR HXR spectra allow us to constrain the possible non-thermal emission that would still be consistent with a null HXR detection. We found that for only one of the features (the jet) was there a potential non-thermal upper limit capable of powering the heating observed. However, even here the non-thermal electron distribution had to be very steep (effectively mono-energetic) with a low energy cut-off between 3-4 keV. The higher temperature or non-thermal sources in the typical quiet Sun features found in this September 2018 data are therefore found to be very weak, if present at all.

Nuclear star clusters (NSCs) are massive star clusters found in all types of galaxies from dwarfs to massive galaxies. Recent studies show that while low-mass NSCs in dwarf galaxies ($M_\text{gal} < 10^{9} M_\odot$) form predominantly out of the merger of globular clusters (GCs), high-mass NSCs in massive galaxies have assembled most of their mass through central enriched star formation. So far, these results of a transition in the dominant NSC formation channel have been based on studies of early-type galaxies and massive late-type galaxies. Here, we present the first spectroscopic analysis of a sample of nine nucleated late-type dwarf galaxies with the aim of identifying the dominant NSC formation pathway. We use integral-field spectroscopy data obtained with the Multi Unit Spectroscopic Explorer (MUSE) instrument to analyse the ages, metallicities, star formation histories, and star formation rates of the NSCs and their surroundings. Our sample includes galaxies with stellar masses $M_\text{gal} = 10^7 - 10^9 M_\odot$ and NSC masses $M_\text{NSC} = 6 \times 10^4 - 6 \times 10^{6} M_\odot$. Although all NSC spectra show emission lines, this emission is not always connected to star formation within the NSC, but rather to other regions along the line-of-sight. The NSC star formation histories reveal that metal-poor and old populations dominate the stellar populations in five NSCs, possibly stemming from the inspiral of GCs. The NSCs of the most massive galaxies in our sample show significant contributions from young and enriched populations that indicate additional mass growth through central star formation. Our results support previous findings of a transition in the dominant NSC formation channel with galaxy mass, showing that the NSCs in low-mass galaxies predominantly grow through the inspiral of GCs, while central star formation can contribute to NSC growth in more massive galaxies.

It is well known that the solar gravitational field can be considered as a telescope with a prime focus at locations beyond 550 au. In this work we present a new derivation of the wave-optical properties of the system, by adapting the arrival-time formalism from gravitational lensing. At the diffraction limit the angular resolution is similar to that of a notional telescope with the diameter of the Sun, and the maximum light amplification is $8{\pi}4GM /(c^2{\lambda})$, enough to detect a 1 W laser on Proxima Centauri b pointed in the general direction of the Sun. Extended sources, however, would be blurred by the wings of the point spread function into the geometrical-optics regime of gravitational lensing. Broad-band sources would have to further contend with the solar corona. Imaging an exoplanet surface as advocated in the literature, without attempting to reach the diffraction limit, appears achievable. For diffraction-limited imaging (sub-km scales from 100 pc) nearby neutron stars appear to be most plausible targets.

We present observations of ammonia emission lines toward the interstellar filament WB~673 hosting the dense clumps WB~673, WB~668, S233-IR and G173.57+2.43. LTE analysis of the lines allows us to estimate gas kinetic temperature ($\lesssim$ 30~K in all the clumps), number density ($7-17\times10^3$~cm$^{-3}$), and ammonia column density ($\approx 1-1.5\times 10^{15}$~cm$^{-2}$) in the dense clumps. We find signatures of collapse in WB 673 and presence of compact spatially unresolved dense clumps in S233-IR. We reconstruct 1D density and temperature distributions in the clumps and estimate their ages using astrochemical modelling. Considering CO, CS, NH$_3$ and N$_2$H$^+$ molecules (plus HCN and HNC for WB~673), we find a chemical age of $t_{\rm chem}=1-3\times 10^5$~yrs providing the best agreement between the simulated and observed column densities in all the clumps. Therefore, we consider $t_{\rm chem}$ as the chemical age of the entire filament. A long preceding low-density stage of gas accumulation in the astrochemical model would break the agreement between the simulated and observed column densities. We suggest that rapid star formation over a $\sim 10^5$~yrs timescale take place in the filament.

We simulate the evolution of relativistic electrons injected into the medium of a small galaxy cluster by a central radio galaxy, studying how the initial jet power affects the dispersal and the emission properties of radio plasma. By coupling passive tracer particles to adaptive-mesh cosmological MHD simulations, we study how cosmic-ray electrons are dispersed as a function of the input jet power. We also investigate how the latter affects the thermal and non-thermal properties of the intracluster medium, with differences discernible up to $\sim$ Gyr after the start of the jet. We evolved the energy spectra of cosmic-ray electrons, subject to energy losses that are dominated by synchrotron and inverse Compton emission as well as energy gains via re-acceleration by shock waves and turbulence. We find that in the absence of major mergers the amount of re-acceleration experienced by cosmic-ray electrons is not enough to produce long-lived detectable radio emissions. However, for all simulations the role of re-acceleration processes is crucial to maintain a significant and volume-filling reservoir of fossil electrons ($\gamma \sim 10^3$) for several Gyrs after the first injection by jets. This is important to possibly explain recent discoveries of cluster-wide emission and other radio phenomena in galaxy clusters.

We present a spectroscopic survey of 248 white dwarf candidates within 40 pc of the Sun; of these 244 are in the southern hemisphere. Observations were performed mostly with the Very Large Telescope (X-Shooter) and Southern Astrophysical Research Telescope. Almost all candidates were selected from $\textit{Gaia}$ Data Release 3 (DR3). We find a total of 246 confirmed white dwarfs, 209 of which had no previously published spectra, and two main-sequence star contaminants. Of these, 100 white dwarfs display hydrogen Balmer lines, 69 have featureless spectra, and two show only neutral helium lines. Additionally, 14 white dwarfs display traces of carbon, while 37 have traces of other elements that are heavier than helium. We observe 36 magnetic white dwarfs through the detection of Zeeman splitting of their hydrogen Balmer or metal spectral lines. High spectroscopic completeness (> 97 per cent) has now been reached, such that we have 1058 confirmed $\textit{Gaia}$ DR3 white dwarfs out of 1083 candidates within 40 pc of the Sun at all declinations.

Context. The magnetic field in the solar atmosphere continually reconnects and accelerates charged particles to high energies. Simulations of the atmosphere in three dimensions that include the effects of accelerated particles can aid our understanding of the interplay between energetic particle beams and the environment where they emerge and propagate. We presented the first attempt at such a simulation in a previous paper, emphasising the physical model of particle beams. However, the numerical implementation of this model is not straightforward due to the diverse conditions in the atmosphere and the way we must distribute computation between multiple CPU cores. Aims. Here, we describe and verify our numerical implementation of energy transport by electron beams in a 3D magnetohydrodynamics code parallelised by domain decomposition. Methods. We trace beam trajectories using a Runge-Kutta scheme with adaptive step length control and integrate deposited beam energy along the trajectories with a hybrid analytical and numerical approach. To parallelise this, we coordinate beam transport across subdomains owned by separate processes using a buffering system designed to optimise data flow. Results. Using an ad hoc magnetic field with analytical field lines as a test scenario, we show that our parallel implementation of adaptive tracing efficiently follows a challenging trajectory with high precision. By timing executions of electron beam transport with different numbers of processes, we found that the processes communicate with minimal overhead but that the parallel scalability is still sublinear due to workload imbalance caused by the uneven spatial distribution of beams.

Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers (M_son < 0.01 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant-Friedrichs-Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume Seven-League Hydro (SLH) code, and it makes use of a low-Mach version of the five-wave Harten-Lax-van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit-explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at M_son ~ 0.001 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.

The existence of high-energy astrophysical neutrinos has been unambiguously demonstrated, but their sources remain elusive. IceCube reported an association of a 290-TeV neutrino with a gamma-ray flare of TXS 0506+056, an active galactic nucleus with a compact radio jet pointing to us. Later, radio blazars were shown to be associated with IceCube neutrino events with high statistical significance. These associations remained unconfirmed with the data of independent experiments. Here we report on the detection of a rare neutrino event with the estimated energy of 224 +- 75 TeV from the direction of TXS 0506+056 by the new Baikal-GVD neutrino telescope in April 2021 followed by a radio flare observed by RATAN-600. This event is the highest-energy cascade detected so far by Baikal-GVD from a direction below horizon. The result supports previous suggestions that radio blazars in general, and TXS 0506+056 in particular, are the sources of high-energy neutrinos, and opens up the cascade channel for the neutrino astronomy.

Optical astronomical images are strongly affected by the point spread function (PSF) of the optical system and the atmosphere (seeing) which blurs the observed image. The amount of blurring depends both on the observed band, and more crucially, on the atmospheric conditions during observation. A typical astronomical image will therefore have a unique PSF that is non-circular and different in different bands. Observations of known stars give us a determination of this PSF. Therefore, any serious candidate for production analysis of astronomical images must take the known PSF into account during the image analysis. So far the majority of applications of neural networks (NN) to astronomical image analysis have ignored this problem by assuming a fixed PSF in training and validation. We present a neural network based deconvolution algorithm based on Deep Wiener Deconvolution Network (DWDN) that takes the PSF shape into account when performing deconvolution as an example of one possible approach to enabling neural network to use the PSF information. We study the performance of several versions of this algorithm under realistic observational conditions in terms of recovery of most relevant astronomical quantities such as colors, ellipticities and orientations. We also investigate the performance of custom loss functions and find that they cause modest improvements in the recovery of astronomical quantities.

Nitrous oxide (N2O) -- a product of microbial nitrogen metabolism -- is a compelling exoplanet biosignature gas with distinctive spectral features in the near- and mid-infrared, and only minor abiotic sources on Earth. Previous investigations of N2O as a biosignature have examined scenarios using Earthlike N2O mixing ratios or surface fluxes, or those inferred from Earth's geologic record. However, biological fluxes of N2O could be substantially higher, due to a lack of metal catalysts or if the last step of the denitrification metabolism that yields N2 from N2O had never evolved. Here, we use a global biogeochemical model coupled with photochemical and spectral models to systematically quantify the limits of plausible N2O abundances and spectral detectability for Earth analogs orbiting main-sequence (FGKM) stars. We examine N2O buildup over a range of oxygen conditions (1%-100% present atmospheric level) and N2O fluxes (0.01-100 teramole per year; Tmol = 10^12 mole) that are compatible with Earth's history. We find that N2O fluxes of 10 [100] Tmol yr$^{-1}$ would lead to maximum N2O abundances of ~5 [50] ppm for Earth-Sun analogs, 90 [1600] ppm for Earths around late K dwarfs, and 30 [300] ppm for an Earthlike TRAPPIST-1e. We simulate emission and transmission spectra for intermediate and maximum N2O concentrations that are relevant to current and future space-based telescopes. We calculate the detectability of N2O spectral features for high-flux scenarios for TRAPPIST-1e with JWST. We review potential false positives, including chemodenitrification and abiotic production via stellar activity, and identify key spectral and contextual discriminants to confirm or refute the biogenicity of the observed N2O.

Scandium is a key element of the Am star phenomenon since its surface under-abundance is one of the criteria that characterise such stars. Thanks to the availability of a sufficiently complete set of theoretical atomic data for this element, reliable radiative accelerations for Sc can now be computed, which allows its behaviour under the action of atomic diffusion to be modelled. We explore the required conditions, in terms of mixing processes or mass loss, for our models to reproduce the observed surface abundances of Sc in Am stars. The models are computed with the Toulouse-Geneva evolution code, which uses the parametric single-valued parameter method for the calculation of radiative accelerations. Fingering mixing is included, using a prescription that comes from 3D hydrodynamical simulations. Other parameter-dependent turbulent mixing processes are also considered. A global mass loss is also implemented. When no mass loss is considered, the observed abundances of Sc are rather in favour of the models whose superficial layers are fully mixed down to the iron accumulation zone, although other mixing prescriptions are also able to reproduce the observations for the most massive model presented here ($2.0 M_\odot$). The models including mass loss with rates in the range of $[10^{-13};10^{-14}] M_\odot$/yr are compatible with some of the observations, while other observations suggest that the mass-loss rate could be lower. The constraints brought by the modelling of Sc are consistent with those derived using other chemical elements.

Context: The high quality of the Gaia mission data is allowing to study the internal kinematics of the Large Magellanic Cloud (LMC) in unprecedented detail, providing insights on the non-axisymmetric structure of its disc. Aims: To define and validate an improved selection strategy to distinguish the LMC stars from the Milky Way foreground. To check the possible biases that assumed parameters or sample contamination from the Milky Way can introduce in the analysis of the internal kinematics of the LMC using Gaia data. Methods: Our selection is based on a supervised Neural Network classifier using as much as of the Gaia DR3 data as possible. We select three samples of candidate LMC stars with different degrees of completeness and purity; we validate them using different test samples and we compare them with the Gaia Collaboration paper sample. We analyse the resulting velocity profiles and maps, and we check how these results change when using also the line-of-sight velocities, available for a subset of stars. Results: The contamination in the samples from Milky Way stars affects basically the results for the outskirts of the LMC, and the absence of line-of-sight velocities does not bias the results for the kinematics in the inner disc. For the first time, we perform a kinematic analysis of the LMC using samples with the full three dimensional velocity information from Gaia DR3. Conclusions: The dynamics in the inner disc is mainly bar dominated; the kinematics on the spiral arm over-density seem to be dominated by an inward motion and a rotation faster than that of the disc in the piece of the arm attached to the bar; contamination of MW stars seem to dominate the outer parts of the disc and mainly affects old evolutionary phases; uncertainties in the assumed disc morphological parameters and line-of-sight velocity of the LMC can in some cases have significant effects. [ABRIDGED]

Supernova (SN) 2021ocs was discovered in the galaxy NGC 7828 within the interacting system Arp 144, and subsequently classified as a normal type-Ic SN around peak brightness. VLT/FORS2 observations in the nebular phase at 150 d reveal that the spectrum is dominated by oxygen and magnesium emission lines of different transitions and ionization states: O I, [O I], [O II], [O III], Mg I, and Mg II. Such a spectrum has no counterpart in the literature, though it bears a few features similar to those of some interacting type Ibn and Icn SNe. Additionally, SN 2021ocs showed a blue color, $(g-r) \lesssim -0.5$ mag, after the peak, atypical for a type-Ic SN. Together with the nebular spectrum, this suggests that SN 2021ocs underwent late-time interaction with an H/He-poor circumstellar medium (CSM), resulting from the pre-SN progenitor mass loss during its final $\sim$1000 days. The strong O and Mg lines and the absence of strong C and He lines suggest that the progenitor star's O-Mg layer is exposed, which places SN 2021ocs as the most extreme case of massive progenitor star's envelope stripping in interacting SNe, followed by type-Icn (stripped C-O layer) and Ibn (stripped He-rich layer) SNe. This is the first time such case is reported in the literature. SN 2021ocs emphasizes the importance of late-time spectroscopy of even seemingly normal SNe, which reveals the inner ejecta and progenitor star's mass loss history.

JWST was designed to peer into the distant universe and study galaxies nearer the beginning of time than previously. Here we report the discovery of 12 galaxy candidates observed 300-600 Myr after the Big Bang with photometric redshifts between z ~ 8.5-13 measured using JWST NIRCam imaging of the galaxy cluster WHL0137 observed in 8 filters spanning 0.8-5.0 $\mu$m, plus 9 HST filters spanning 0.4-1.7 $\mu$m. Three of these candidates are gravitationally lensed by the foreground galaxy cluster and have magnifications of $\mu \sim 3 - 8$. The remaining nine candidates are located in a second JWST NIRCam module, centered ~29' from the cluster center, with expected magnifications of $\mu$ <~ 1.1. Our sample of high-redshift candidates have observed F200W AB magnitudes between 25.9 and 28.1 mag and intrinsic F200W AB magnitudes between 26.4 and 29.7 mag ($M_{UV}$ = -22.5 to -17). We find the stellar masses of these galaxies are in the range $\log M_{*}/M_{\odot}$ = 8 - 9, and down to 7.5 for the lensed galaxies. All are young with mass-weighted ages < 100 Myr, low dust content $A_V$ < 0.15 mag, and high specific star formation rates sSFR ~10-50 Gyr$^{-1}$ for most. One z ~ 9 candidate is consistent with an age < 5 Myr and a sSFR ~250 Gyr$^{-1}$, as inferred from a strong F444W excess, implying [OIII]+H-beta rest-frame equivalent width ~2000 Angstrom, although an older and redder z~ 10 object is also allowed. Another z~9 candidate ID9356 is lensed into an arc 2.6" long by the effects of strong gravitational lensing ($\mu$~8), and has at least two bright knots of unevenly distributed star formation. This arc is the most spatially-resolved galaxy at z~9 known to date, revealing structures ~30 pc across. Follow-up spectroscopy of WHL0137 with JWST/NIRSpec is planned for later this year, which will validate some of these candidates and study their physical properties in more detail.

While previous studies have shown a strong preference for a future mid-infrared nulling interferometer space mission to detect habitable zone planets around M-dwarfs, we here focus on a more conservative approach towards the concept of habitability and present yield estimates for two stellar samples consisting of nearby (d<20 pc) Sun-like stars (4800-6300 K) and nearby FGK-type stars (3940-7220 K) accessible to such a mission. Our yield estimates are based on recently derived rocky planet occurrence rates from the Kepler mission and our LIFE exoplanet observation simulation tool LIFEsim, which includes all major astrophysical noise sources, but no instrumental noise sources, yet. Depending on a pessimistic or optimistic extrapolation of the Kepler results, we find that during a 2.5-year search phase LIFE could detect between ~10-16 (average) or ~5-34 (including 1-$\sigma$ uncertainties) rocky planets (0.5-1.5 R${}_\rm{Earth}$) within the optimistic HZ of Sun-like stars and between ~4-6 (average) or ~1-13 (including 1-$\sigma$ uncertainties) exo-Earth candidates (EECs) assuming four collector spacecraft equipped with 2 m mirrors and a conservative instrument throughput of 5%. With D=3.5 m or 1 m mirrors, the yield $Y$ changes strongly following approximately $Y \propto D^{3/2}$. With the larger sample of FGK-type stars, the yield increases to ~16-22 (average) rocky planets within the optimistic HZ and ~5-8 (average) EECs. Furthermore, we find that, besides the mirror diameter, the yield depends strongly on the total throughput, but only weakly on the exozodiacal dust level and the accessible wavelength range of the mission. When focusing entirely on Sun-like stars, larger mirrors (~3 m with 5% total throughput) or a better total throughput (~20% with 2 m mirrors) are required to detect a statistically relevant sample of ~30 rocky planets within the optimistic HZ.

The majority of brown dwarfs show some level of photometric or spectro-photometric variability in different wavelength ranges. This variability allow us to trace the 3D atmospheric structures of variable brown dwarfs and directly-imaged exoplanets with radiative-transfer models and mapping codes. Nevertheless, to date, we do not have an informed method to pre-select the brown dwarfs that might show a higher variability amplitude for a thorough variability study. In this work, we designed and tested near-infrared spectral indices to pre-select the most likely variable mid- and late-T dwarfs, which overlap in effective temperatures with directly-imaged exoplanets. We used time-resolved near-infrared Hubble Space Telescope Wide Field Camera 3 spectra of a T6.5 dwarf, 2MASS J22282889--431026, to design our novel spectral indices. We tested these spectral indices on 26 T5.5--T7.5 near-infrared SpeX/IRTF spectra, and we provided eight new mid- and late-T variable candidates. We estimated the variability fraction of our sample in $38^{+4}_{-30}$%, which agrees with the variability fractions provided by Metchev, et. al, (2015) for mid-to late-T dwarfs. In addition, two of the three previously known variables in our SpeX spectra sample are flagged as variable candidates by our indices. Similarly, all seven known non-variable in our sample are flagged as non-variable objects by our indices. These results suggest that our spectral indices might be used to find variable mid- and late-T brown dwarf variables. These indices may be crucial in the future to select cool directly-imaged exoplanets for variability studies.

We expand our analysis of Newtonian Fractional-Dimension Gravity (NFDG), an extension of the classical laws of Newtonian gravity to lower dimensional spaces, including those with fractional (i.e., non-integer) dimension. We apply our model to four rotationally supported galaxies (NGC 5033, NGC 6674, NGC 5055, NGC 1090), in addition to other three galaxies (NGC 7814, NGC 6503, NGC 3741) which were analyzed in previous studies. NFDG is able to fit the rotation curves of all these galaxies without any dark matter component. We also investigate the possible violation of the strong equivalence principle, in relation to the External Field Effect (EFE), i.e., the dependence of the internal motion of a self-gravitating system under freefall on an external gravitational field. This effect is not present in Newtonian or Einstein gravity, but is predicted by some alternative theories of gravity. On the contrary, we show that NFDG does not imply the EFE, at least for values of the fractional dimension in the range $1 \leq D \leq 3$. Using improved NFDG numerical computations, we analyze the rotation curves of the aforementioned galaxies and obtain perfect fits to the experimental data, by using a fractional-dimension function $D\left (R\right )$ which characterizes each individual galaxy. In the galactic sample studied here with NFDG methods, we do not detect any significant differences between galaxies that are supposed to show/not show the EFE according to other alternative theories. A larger sample of galaxies will be needed to fully determine the absence of any external field effect in NFDG.

Laboratory experiments have confirmed that the radiolytic decay rate of astrochemical ice analogues is dependent upon the solid phase of the target ice, with some crystalline molecular ices being more radio-resistant than their amorphous counterparts. The degree of radio-resistance exhibited by crystalline ice phases is dependent upon the nature, strength, and extent of the intermolecular interactions that characterise their solid structure. For example, it has been shown that crystalline CH3OH decays at a significantly slower rate when irradiated by 2 keV electrons at 20 K than does the amorphous phase due to the stabilising effect imparted by the presence of an extensive array of strong hydrogen bonds. These results have important consequences for the astrochemistry of interstellar ices and outer Solar System bodies, as they imply that the chemical products arising from the irradiation of amorphous ices (which may include prebiotic molecules relevant to biology) should be more abundant than those arising from similar irradiations of crystalline phases. In this present study, we have extended our work on this subject by performing comparative energetic electron irradiations of the amorphous and crystalline phases of the sulphur-bearing molecules H2S and SO2 at 20 K. We have found evidence for phase-dependent chemistry in both these species, with the radiation-induced exponential decay of amorphous H2S being more rapid than that of the crystalline phase, similar to the effect that has been previously observed for CH3OH. For SO2, two fluence regimes are apparent: a low-fluence regime in which the crystalline ice exhibits a rapid exponential decay while the amorphous ice possibly resists decay, and a high-fluence regime in which both phases undergo slow exponential-like decays.

We present a fit to observational data in an asymmetric self-interacting dark matter model using our recently calculated cross sections that incorporate both $t$-channel and $u$-channel exchanges in the scattering of identical particles. We find good fits to the data ranging from dwarf galaxies to galaxy clusters, and equivalent relative velocities from $\sim 20$ km/sec to $\gtrsim 10^3$ km/s. We compare our results with previous fits that used only $t$-channel exchange contributions to the scattering.

We establish a generic, fully-relativistic formalism to study gravitational-wave emission by extreme-mass-ratio systems in spherically-symmetric, non-vacuum black-hole spacetimes. The potential applications to astrophysical setups range from black holes accreting baryonic matter to those within axionic clouds and dark matter environments, allowing to assess the impact of the galactic potential, of accretion, gravitational drag and halo feedback on the generation and propagation of gravitational-waves. We apply our methods to a black hole within a halo of matter. We find fluid modes imparted to the gravitational-wave signal (a clear evidence of the black hole fundamental mode instability) and the tantalizing possibility to infer galactic properties from gravitational-wave measurements by sensitive, low-frequency detectors.

We obtain the analytic solution of the Friedmann equation for fully realistic cosmologies including radiation, non-relativistic matter, a cosmological constant $\lambda$ and arbitrary spatial curvature $\kappa$. The general solution for the scale factor $a(\tau)$, with $\tau$ the conformal time, is an elliptic function, meromorphic and doubly periodic in the complex $\tau$-plane, with one period along the real $\tau$-axis, and the other along the imaginary $\tau$-axis. The periodicity in imaginary time allows us to compute the thermodynamic temperature and entropy of such spacetimes, just as Gibbons and Hawking did for black holes and the de Sitter universe. The gravitational entropy favors universes like our own which are spatially flat, homogeneous, and isotropic, with a small positive cosmological constant.

We revisit the original proposal of cosmological relaxation of the electroweak scale by Graham, Kaplan and Rajendran in which the Higgs mass is scanned during inflation by an axion field, the relaxion. We investigate the regime where the relaxion is subject to large fluctuations during inflation. The stochastic dynamics of the relaxion is described by means of the Fokker-Planck formalism. We derive a new stopping condition for the relaxion taking into account transitions between the neighboring local minima of its potential. Relaxion fluctuations have important consequences even in the "classical-beats-quantum" regime. We determine that for a large Hubble parameter during inflation, the random walk prevents the relaxion from getting trapped at the first minimum. The relaxion stops much further away, where the potential is less shallow. Interestingly, this essentially jeopardises the "runaway relaxion" threat from finite-density effects, restoring most of the relaxion parameter space. We also explore the "quantum-beats-classical" regime, opening large new regions of parameter space. We investigate the consequences for both the QCD and the non-QCD relaxion. The misalignment of the relaxion due to fluctuations around its local minimum opens new phenomenological opportunities.

Axion clumps in an external magnetic field can emit electromagnetic radiation which causes them to decay. In the presence of a plasma, such radiation can become resonant if the clump frequency matches the plasma frequency. Typically, the decay or back-reaction of the clump is ignored in the literature when analyzing such radiation. In this paper we present a self consistent, semi-analytic approach which captures axion back-reaction using energy conservation. We find that inclusion of back-reaction changes the clump frequency over time enabling clumps with a range of different initial frequencies to become resonant at some point in their time evolution.

While cosmological and astrophysical probes suggest that dark matter would make up for 85% of the total matter content of the Universe, the determination of its nature remains one of the greatest challenges of fundamental physics. Assuming the $\Lambda$CDM cosmological model, Weakly Interacting Massive Particles would annihilate into Standard Model particles, yielding $\gamma$-rays, which could be detected by ground-based telescopes. Dwarf spheroidal galaxies represent promising targets for such indirect searches as they are assumed to be highly dark matter dominated with the absence of astrophysical sources nearby. Previous studies have led to upper limits on the annihilation cross-section assuming single exclusive annihilation channels. In this work, we consider a more realistic situation and take into account the complete annihilation pattern within a given particle physics model. This allows us to study the impact on the derived upper limits on the dark matter annihilation cross-section from a full annihilation pattern compared to the case of a single annihilation channel. We use mockdata for the Cherenkov Telescope Array simulating the observations of the promising dwarf spheroidal galaxy Sculptor. We show the impact of considering the full annihilation pattern within a simple framework where the Standard Model of particle physics is extended by a singlet scalar. Such a model shows new features in the shape of the predicted upper limit which reaches a value of $\langle \sigma v \rangle = 3.8\times10^{-24}~\rm{cm}^{-3}\rm{s}^{-1}$ for a dark matter mass of 1 TeV at 95% confidence level. We suggest to consider the complete particle physics information in order to derive more realistic limits.

Spectral shaping is critical to many fields of science. In astronomy for example, the detection of exoplanets via the Doppler effect hinges on the ability to calibrate a high resolution spectrograph. Laser frequency combs can be used for this, but the wildly varying intensity across the spectrum can make it impossible to optimally utilize the entire comb, leading to a reduced overall precision of calibration. To circumvent this, astronomical applications of laser frequency combs rely on a bulk optic setup which can flatten the output spectrum before sending it to the spectrograph. Such flatteners require complex and expensive optical elements like spatial light modulators and have non-negligible bench top footprints. Here we present an alternative in the form of an all-photonic spectral shaper that can be used to flatten the spectrum of a laser frequency comb. The device consists of a circuit etched into a silicon nitride wafer that supports an arrayed-waveguide grating to disperse the light over hundreds of nanometers in wavelength, followed by Mach-Zehnder interferometers to control the amplitude of each channel, thermo-optic phase modulators to phase the channels and a second arrayed-waveguide grating to recombine the spectrum. The demonstrator device operates from 1400 to 1800 nm (covering the astronomical H band), with twenty 20 nm wide channels. The device allows for nearly 40 dBs of dynamic modulation of the spectrum via the Mach-Zehnders , which is greater than that offered by most spatial light modulators. With a superluminescent diode, we reduced the static spectral variation to ~3 dB, limited by the properties of the components used in the circuit and on a laser frequency comb we managed to reduce the modulation to 5 dBs, sufficient for astronomical applications.

We examined soft X-ray emission by the solar wind charge-exchange process around the Earth's magnetosphere using a global magnetohydrodynamic simulation model. The dayside magnetopause reconnection heats and accelerates the plasma whereby the X-ray emission becomes as bright as $\sim 6 \times 10^{-6} {\rm\ eV}\ {\rm cm}^{-3}\ {\rm s}^{-1}$ under the southward interplanetary magnetic field conditions. In particular, under low plasma-$\beta$ solar wind conditions, we found that the X-ray intensity reflects the bulk motion of outflows from the reconnection region. We propose that this particular solar wind condition would allow visualization of the mesoscale magnetopause reconnection site, as observed in the solar corona.

It is possible that the strongest interactions between dark matter and the Standard Model occur via the neutrino sector. Unlike gamma rays and charged particles, neutrinos provide a unique avenue to probe for astrophysical sources of dark matter, since they arrive unimpeded and undeflected from their sources. Previously, we reported on annihilations of dark matter to neutrinos; here, we review constraints on the decay of dark matter into neutrinos over a range of dark matter masses from MeV to ZeV, compiling previously reported limits, exploring new electroweak corrections and computing constraints where none have been computed before. We examine the expected contributions to the neutrino flux at current and upcoming neutrino experiments as well as photons from electroweak emission expected at gamma-ray telescopes, leading to constraints on the dark matter decay lifetime, which ranges from $\tau \sim 1.2\times10^{21}$ s at 10 MeV to $1.5\times10^{29}$s at 1 PeV.

Since the first in-situ measurements of the altitude profile of upper atmospheric density and composition were carried out by the Viking lander missions in 1976, similar data are continuously gathered by MAVEN and MOM spacecraft orbiting Mars since their launch in September 2014 with mass spectrometers and other related payloads. Using near-simultaneous observations by the two orbiters, it is seen that both data sets indicate significant day-to-day variations of Argon density profiles in the thermosphere-exosphere, 150-300 km region, during the period 1-15, June 2018, when the solar EUV radiation did not show any appreciable change but the solar wind energetic particle fluxes did so. Extending this study to include the other parent atmospheric constituents carbon dioxide, helium, nitrogen and their photochemical products atomic oxygen, and carbon monoxide during the same period it is found that the density profiles of carbon dioxide and atomic oxygen also show similar variations with carbon dioxide densities showing an increasing trend similar to Argon, but a reversal of this trend for atomic oxygen densities. Using insitu and near simultaneous measurements of solar EUV fluxes and the solar wind plasma velocities and densities near MAVEN periapsis it is noted that, unlike the solar EUV radiation, solar wind parameters showed a decrease by a factor of 2-3. Hence, it is inferred that the energetic and penetrating solar wind charged particle impact-driven dissociation, ionisation and ion-chemical processes could decrease the carbon dioxide densities leading to an increase in atomic oxygen densities. This result is also discussed from the considerations of the proton gyro radius effect, pickup ions, sputtering, energetic neutral atoms driven ionisation and ion losses. Further data and modelling efforts would be necessary to confirm this finding.

In this article we discuss some consequences of the well-known proposition of Fritz Zwicky [1], published in the nineteen thirties, that Dark Matter `mimics' the inertia-gravitational behaviour of usual matter. In particular, we consider some special dynamical regions such as those of the Ring Systems of the gaseous giants at the edge of the Planetary System. This article is a continuation of an earlier paper [2], where it was shown that gravitationally interacting particles may remain near the Lagrange Points L4 and L5 for many thousands of years. This provides enough time for the Dark Matter, if present there, to interact with the usual matter. We discuss also a number of questions related to places which might be considered singular in the mathematical sense.

The closed-form solution of the 1.5 post-Newtonian (PN) accurate binary black hole (BBH) Hamiltonian system has proven to be difficult to obtain for a long time since its introduction in 1966. Closed-form solutions of the PN BBH systems with arbitrary parameters (masses, spins, eccentricity) are required for modeling the gravitational waves (GWs) emitted by them. Accurate models of GWs are crucial for their detection by LIGO/Virgo and LISA. Only recently, two solution methods for solving the BBH dynamics were proposed in arXiv:1908.02927 (without using action-angle variables), and arXiv:2012.06586, arXiv:2110.15351 (action-angle based). This paper combines the ideas laid out in the above articles, fills the missing gaps and provides the two solutions which are fully 1.5PN accurate. We also present a public Mathematica package BBHpnToolkit which implements these two solutions and compares them with a fully numerical treatment. The level of agreement between these solutions provides a numerical verification for all the five actions constructed in arXiv:2012.06586, and arXiv:2110.15351. This paper hence serves as a stepping stone for pushing the action-angle-based solution to 2PN order via canonical perturbation theory.

We present a convolutional neural network, designed in the auto-encoder configuration that can detect and denoise astrophysical gravitational waves from merging black hole binaries, orders of magnitude faster than the conventional matched-filtering based detection that is currently employed at advanced LIGO (aLIGO). The Neural-Net architecture is such that it learns from the sparse representation of data in the time-frequency domain and constructs a non-linear mapping function that maps this representation into two separate masks for signal and noise, facilitating the separation of the two, from raw data. This approach is the first of its kind to apply machine learning based gravitational wave detection/denoising in the 2D representation of gravitational wave data. We applied our formalism to the first gravitational wave event detected, GW150914, successfully recovering the signal at all three phases of coalescence at both detectors. This method is further tested on the gravitational wave data from the second observing run ($O2$) of aLIGO, reproducing all binary black hole mergers detected in $O2$ at both the aLIGO detectors. The Neural-Net seems to have uncovered a pattern of 'ringing' after the ringdown phase of the coalescence, which is not a feature that is present in the conventional binary merger templates. This method can also interpolate and extrapolate between modeled templates and explore gravitational waves that are unmodeled and hence not present in the template bank of signals used in the matched-filtering detection pipelines. Faster and efficient detection schemes, such as this method, will be instrumental as ground based detectors reach their design sensitivity, likely to result in several hundreds of potential detections in a few months of observing runs.

Black holes below Chandrasekhar mass limit (1.4 $M_{\odot}$) can not be produced via any standard stellar evolution. Recently, gravitational wave experiments have also discovered unusually low mass black holes whose origin is yet to be known. We propose a simple yet novel formation mechanism of such low mass black holes. Non-annihilating particle dark matter, owing to their interaction with stellar nuclei, can gradually accumulate inside compact stars, and eventually swallows them to low mass black holes, ordinarily impermissible by the Chandrasekhar limit. We point out several avenues to test this proposal, concentrating on the cosmic evolution of the binary merger rates.

The fundamental nature of Dark Matter is a central theme of the Snowmass 2021 process, extending across all frontiers. In the last decade, advances in detector technology, analysis techniques and theoretical modeling have enabled a new generation of experiments and searches while broadening the types of candidates we can pursue. Over the next decade, there is great potential for discoveries that would transform our understanding of dark matter. In the following, we outline a road map for discovery developed in collaboration among the frontiers. A strong portfolio of experiments that delves deep, searches wide, and harnesses the complementarity between techniques is key to tackling this complicated problem, requiring expertise, results, and planning from all Frontiers of the Snowmass 2021 process.