Papers for Friday, Oct 01 2021

The recent discovery of the binary black hole (BBH) merger event GW190521, between two black holes (BHs) of $\sim100M_\odot$, and as well as other massive BBH merger events involving BHs within the pair-instability supernova (PSN) mass gap have sparked widespread debate on the origin of such extreme gravitational-wave (GW) events. In this study, I investigate whether dynamical interactions in young massive clusters (YMCs) serves as a viable scenario for assembling PSN-gap BBH mergers. To that end, I explore a grid of 24 new evolutionary models of a representative YMC of initial mass $M_{\rm cl}=7.5\times10^4M_\odot$ ($N\approx1.28\times10^5$) and size $r_h=2$ pc, with all BH progenitor stars being initially in primordial binaries. All cluster models are evolved with the direct, relativistic N-body code NBODY7 incorporating up to date remnant formation, BH natal spin, and general-relativistic (GR) merger recoil schemes. The BBH mergers from these model cluster computations agree well with the GWTC-2 and GWTC-2.1 events' masses and effective spin parameters. In particular, GW190521-like, i.e., $\sim200M_\odot$, low aligned spin events are produced via dynamical merger among BHs derived from star-star merger products. GW190403-like, i.e., PSN-gap, highly asymmetric and aligned events result from mergers involving BHs that are spun up via matter accretion or binary interaction. The resulting differential merger rate density within the PSN gap well accommodates that from GWTC-2. Particularly, the models well reproduce the LVK-estimated merger rate density of GW190521-like events. This study demonstrates that, subject to model uncertainties, the tandem of massive binary evolution and dynamical interactions in low metallicity YMCs in the Universe can plausibly produce GR mergers involving PSN-gap BHs and in rates consistent with that from to-date GW observations. [Abridged]

We present an analysis of 991 heartbeat stars (HBSs) from the OGLE Collection of Variable Stars (OCVS). The sample consists of 512 objects located toward the Galactic bulge (GB), 439 in the Large Magellanic Cloud (LMC) and 40 in the Small Magellanic Cloud (SMC). We model the $I$-band OGLE light curves using an analytical model of flux variations, reflecting tidal deformations between stars. We present distributions of the model parameters that include the eccentricity, orbital inclination, and argument of the periastron, but also the period-amplitude diagrams. On the Hertzsprung-Russell (HR) diagram, our HBS sample forms two separate groups of different evolutionary status. The first group of about 90 systems, with short orbital periods ($P\lesssim50$~days), consists of an early-type primary star lying on (or close to) the main sequence (MS). The second group of about 900 systems, with long orbital periods ($P\gtrsim100$~days), contains a red giant (RG). The position of RG HBSs on the period-luminosity diagram strongly indicates their binary nature. They appear to be a natural extension of confirmed binary systems that include the OGLE ellipsoidal and Long Secondary Period (LSP) variables. We also present a time-series analysis leading to detection of tidally-excited oscillations (TEOs). We identify such pulsations in about 5\% of stars in the sample with a total number of 78 different modes. This first relatively large homogeneous sample of TEOs allowed us to construct a diagram revealing the correlation between the TEO's orbital harmonic number and the eccentricity of the host binary system.

Time-delay strong lensing (TDSL) is a powerful probe of the current expansion rate of the Universe. However, in light of the discrepancies between early and late-time cosmological studies, efforts revolve around the characterisation of systematic uncertainties in the methods. Here, we focus on the mass-sheet degeneracy (MSD), which is considered a significant source of systematics in TDSL, and aim to assess the constraining power provided by IFU stellar kinematics. We approximate the MSD with a cored, two-parameter extension to the lensing mass profiles (with core radius $r_{\rm c}$ and mass-sheet parameter $\lambda_{\rm int}$). In addition, we utilise mock IFU stellar kinematics of time-delay strong lenses, given the prospects of obtaining such data with JWST. We construct joint strong lensing and stellar dynamical models, where the time delays, mock imaging and IFU observations are used to constrain the mass profile of lens galaxies, and yield joint constraints on the time-delay distance ($D_{\Delta t}$) and angular diameter distance ($D_{\rm d}$) to the lens. We find that mock JWST-like stellar kinematics constrain the internal mass sheet and limit its contribution to the uncertainties of $D_{\Delta t}$ and $D_{\rm d}$, each at the < 4% level, without assumptions on the background cosmological model. These distance constraints would translate to a < 4% precision measurement on $H_{\rm 0}$ in flat $\Lambda CDM$ for a single lens. Our study shows that IFU stellar kinematics of time-delay strong lenses will be key in lifting the MSD on a per lens basis, assuming reasonable core sizes. However, even in the limit of infinite $r_{\rm c}$, where $D_{\Delta t}$ is degenerate with $\lambda_{\rm int}$, stellar kinematics of the deflector, time delays and imaging data will provide powerful constraints on $D_{\rm d}$, which becomes the dominant source of information in the cosmological inference.

We present a collection of 991 heartbeat star (HBS) candidates found in the Optical Gravitational Lensing Experiment (OGLE) project data archive. We discuss the selection process of the HBS candidates and the structure of the catalog itself. It consists of 512 stars located toward the Galactic bulge (GB), 439 stars located in the Large Magellanic Cloud (LMC), and 40 in the Small Magellanic Cloud (SMC). The collection contains two large groups of HBSs with different physical properties. The main distinction between the two groups is the evolutionary status of the primary star. The first group of about 100 systems contains a hot main-sequence (MS) or a Hertzsprung-gap primary star, while the second group of about 900 systems includes a red giant (RG). For each star, we provide two-decade-long time-series photometry, in the Cousins $I$- and Johnson $V$-band filters, obtained by the OGLE project. We also present basic observational information as well as orbital parameters derived from the light curve modeling.

The majority of existing results for the kilonova (or macronova) emission from material ejected during a neutron-star (NS) merger is based on (quasi-)one-zone models or manually constructed toy-model ejecta configurations. In this study we present a kilonova analysis of the material ejected during the first ~10ms of a NS merger, called dynamical ejecta, using directly the outflow trajectories from general relativistic smoothed-particle hydrodynamics simulations including a sophisticated neutrino treatment and the corresponding nucleosynthesis results, which have been presented in Part I of this study. We employ a multi-dimensional two-moment radiation transport scheme with approximate M1 closure to evolve the photon field and use a heuristic prescription for the opacities found by calibration with atomic-physics based reference results. We find that the photosphere is generically ellipsoidal but augmented with small-scale structure and produces emission that is about 1.5-3 times stronger towards the pole than the equator. The kilonova typically peaks after 0.7-1.5days in the near-infrared frequency regime with luminosities between 3-7x10^40erg/s and at photospheric temperatures of 2.2-2.8x10^3K. A softer equation of state or higher binary-mass asymmetry leads to a longer and brighter signal. Significant variations of the light curve are also obtained for models with artificially modified electron fractions, emphasizing the importance of a reliable neutrino-transport modeling. None of the models investigated here, which only consider dynamical ejecta, produces a transient as bright as AT2017gfo. The near-infrared peak of our models is incompatible with the early blue component of AT2017gfo.

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, $ l $, of the collapsing material and the magnetic field, $ B $, on the BH horizon: (i) Subrelativistic outflow (low $ l $ and high $ B $), (ii) Stationary accretion shock instability (SASI; high $ l $ and low $ B $), (iii) Relativistic jets (high $ l $ and high $ B $). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices ($ \gtrsim 2 $), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the lightcurve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova+GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disk-jet flows, even in the absence of any tilt in the progenitor.

We introduce the public version of the BAyesian STellar Algorithm (BASTA), an open-source code written in {\tt Python} to determine stellar properties based on a set of astrophysical observables. BASTA has been specifically designed to robustly combine large datasets that include asteroseismology, spectroscopy, photometry, and astrometry. We describe the large number of asteroseismic observations that can be fit by the code and how these can be combined with atmospheric properties (as well as parallaxes and apparent magnitudes), making it the most complete analysis pipeline available for oscillating main-sequence, subgiant, and red giant stars. BASTA relies on a set of pre-built stellar isochrones or a custom-designed library of stellar tracks which can be further refined using our interpolation method (both along and across stellar tracks/isochrones). We perform recovery tests with simulated data that reveal levels of accuracy at the few percent level for radii, masses, and ages when individual oscillation frequencies are considered, and show that asteroseismic ages with statistical uncertainties below 10% are within reach if our stellar models are reliable representations of stars. BASTA is extensively documented and includes a suite of examples to support easy adoption and further development by new users.

Filaments and clusters of the cosmic web have an impact on the properties of galaxies, switching off their star-formation, contributing to the build-up of their stellar mass, and influencing the acquisition of their angular momentum. In this work we make use of the IllustrisTNG simulation, coupled with the DisPerSE cosmic web extraction algorithm, to test which is the galaxy property most affected by the cosmic web and, conversely, to assess the differential impact of the various cosmic web features on a given galaxy property. Our aim is to use this information to better understand galaxy evolution and to identify on which galaxy property future efforts should focus to detect the cosmic web from the galaxy distribution. We provide a comprehensive analysis of the relation between galaxy properties and cosmic web features. We also perform extensive tests in which we try to disentangle the effect of local overdensities of galaxies on their properties from the effect of the large scale structure environment. Our results show that star-formation is the quantity that shows the strongest variation with the distances from the cosmic web features, but it is also the one that shows the strongest relation to the local environment of galaxies. On the other hand, the direction of the angular momentum of galaxies is the property that shows the weakest trends with distance from cosmic web features, while also being more independent from the local environment of galaxies. We conclude that the direction of the angular momentum of galaxies and its use to improve our detection of the cosmic web features could be the focus of futures studies benefitting from larger statistical samples.

We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of HI shells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows near $z \sim 2$. We find that high-redshift galaxies may 'capture' 20-50% of their feedback momentum in the dense ISM (with the remainder escaping into the nearby CGM), whereas local galaxies may contain $\lesssim$10% of their feedback momentum from the central starburst. Using azimuthally averaged galaxy properties, we predict that most superbubbles stall and fragment within the ISM, and that this occurs at, or near, the gas scale height. We find good agreement when comparing our predicted bubble radii and velocities with observed HI bubbles/holes and find that most will fragment very near the gas scale height or are able to break-out and drive local galactic fountains. Additionally, we demonstrate that models with constant star cluster formation efficiency per Toomre mass are inconsistent with the occurrence of outflows from high-$z$ starbursts and local circumnuclear regions.

In core-collapse supernovae or compact binary merger remnants, neutrino-neutrino refraction can spawn fast pair conversion of the type $\nu_e \bar\nu_e \leftrightarrow \nu_x \bar\nu_x$ (with $x=\mu, \tau$), governed by the angle-dependent density matrices of flavor lepton number. All angle modes evolve coherently for a homogeneous and axially symmetric two-flavor system, and the nonlinear equations of motion are equivalent to those of a gyroscopic pendulum. Our main innovation is to identify the real part of the eigenfrequency resulting from the linear normal-mode analysis of the neutrino system as the "pendulum spin." This is the elusive characteristic of the lepton-number angle distribution that determines the depth of conversion, without solving the nonlinear evolution equations.

We investigate the relationship between dust attenuation and stellar mass ($M_*$) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey at $z\sim2.3$ and the Sloan Digital Sky Survey (SDSS) at $z\sim0$, augmenting the latter optical dataset with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-field Infrared Survey Explorer (WISE) photometry from the GALEX-SDSS-WISE Catalog. We quantify dust attenuation using both spectroscopic measurements of H$\alpha$ and H$\beta$ emission lines, and photometric measurements of the rest-UV stellar continuum. The H$\alpha$/H$\beta$ ratio is used to determine the magnitude of extinction at the wavelength of H$\alpha$, $A_{{\rm H}\alpha}$. Rest-UV colors and spectral-energy-distribution fitting are used to estimate $A_{1600}$, the magnitude of extinction at a rest wavelength of 1600\AA. As in previous work, we find a lack of significant evolution in the relation between dust attenuation and $M_*$ over the redshift range $z\sim0$ to $z\sim2.3$. Folding in the latest estimates of the evolution of $M_{{\rm dust}}$, $({M_{{\rm dust}}}/{M_{{\rm gas}}})$, and gas surface density at fixed $M_*$, we find that the expected $M_{{\rm dust}}$ and dust mass surface density are both significantly higher at $z\sim2.3$ than at $z\sim0$. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancy in attenuation vs. $M_*$, it is essential to determine the relationship between metallicity and $({M_{{\rm dust}}}/{M_{{\rm gas}}})$, the dust mass absorption coefficient, and dust geometry, and the evolution of these relations and quantities from $z\sim0$ to $z\sim2.3$.

Broadband spectral studies of Type-I X-ray bursts can put strong constraints on the physics of burst spectra as well as their interaction with the environment. We present the results obtained from the broadband time-resolved spectroscopy of 15 thermonuclear bursts detected simultaneously from the neutron star atoll source 4U 1636-536 using LAXPC and SXT onboard AstroSat. During the observations with AstroSat, the Low mass X-ray binary (LMXB) 4U 1636-536 is observed to show a modest spectral evolution within the island state. The broadband burst spectra are observed to show an excess in addition to the thermal emission from the neutron star surface near the peak of the bursts. We investigate the interpretation of the excess observed near the peak of the burst as re-emission/reprocessing of the photons by the accretion disk/corona or scattering of the photons in the neutron star atmosphere or the enhanced persistent emission due to Poynting-Robertson drag. This is the first reported broadband simultaneous study of Type-I bursts using LAXPC and SXT onboard AstroSat. This kind of study may provide a better understanding of the burst-accretion interaction and how the bursts influence the overall accretion process contributed by the accretion disk as well as the corona.

Mass measurements and absorption line studies indicate that the stellar initial mass function (IMF) is bottom-heavy in the central regions of many early-type galaxies, with an excess of low mass stars compared to the IMF of the Milky Way. Here we test this hypothesis using a method that is independent of previous techniques. Low mass stars have strong chromospheric activity characterized by non-thermal emission at short wavelengths. Approximately half of the UV flux of M dwarfs is contained in the $\lambda{}1215.7$ Ly$\alpha$ line, and we show that the total Ly$\alpha$ emission of an early-type galaxy is a sensitive probe of the IMF with a factor of $\sim 2$ flux variation in response to plausible variations in the number of low mass stars. We use the Cosmic Origins Spectrograph on the Hubble Space Telescope to measure the Ly$\alpha$ line in the centers of the massive early-type galaxies NGC1407 and NGC2695. We detect Ly$\alpha$ emission in both galaxies and demonstrate that it originates in stars. We find that the Ly$\alpha$ to i-band flux ratio is a factor of 2.0$\pm$0.4 higher in NGC1407 than in NGC2695, in agreement with the difference in their IMFs as previously determined from gravity-sensitive optical absorption lines. Although a larger sample of galaxies is required for definitive answers, these initial results support the hypothesis that the IMF is not universal but varies with environment.

We analyse VLT/MUSE observations of PGC 046832, the brightest cluster galaxy of Abell 3556. The velocity structure of this galaxy is startling and shows two reversals in sign along the minor axis, and one along the major axis. We use triaxial Schwarzschild models to infer the intrinsic shape, central black hole mass and orbit distribution of this galaxy. The shape determination suggests that the galaxy is highly triaxial in the centre (almost prolate) but has a low triaxiality (almost oblate) in the outer parts. The orbit distribution of the best-fit Schwarzschild model shows that the kinematic reversal along the projected minor axis is driven by a slight asymmetry in the distribution and amount of long axis tubes in the inner parts. The kinematic reversals along the projected major axis are driven by a high fraction of counter-rotating orbits at intermediate radii in the galaxy. Using chemical tagging of orbits in the Schwarzschild model, we do not find evidence for any association of these orbits with specific stellar population parameters. Although the inner part of the galaxy almost certainty formed through one or more dry mergers producing the prolate shape, the outer parts are consistent with both accretion and in situ formation. While axisymmetric models suggests the presence of a supermassive black hole with mass $\sim 6 \times 10^9$M$_{\odot}$ and $\sim 10^{10}$M$_{\odot}$ (with Schwarzschild and Jeans modelling, resp), triaxial Schwarzschild models provide only an upper limit of $\sim 2 \times 10^9$M$_{\odot}$.

We present optical and near-infrared photometric and spectroscopic observations of the fast-declining Type Ia Supernova (SN) 2015bo. SN 2015bo is under-luminous (M$_B$ = -17.50 $\pm$ 0.13 mag) and has a fast-evolving light curve ($\Delta \mathrm{m}_{15}\mathrm{(B)}$ = 1.91 $\pm$ 0.01 mag and $s_{BV}$ = 0.48 $\pm$ 0.01). It has a unique morphology in the $V-r$ color curve, where it is bluer than all other SNe in the comparison sample. A $^{56}$Ni mass of 0.17 $\pm$ 0.03 $M_{\odot}$ was derived from the peak bolometric luminosity, which is consistent with its location on the luminosity-width relation. Spectroscopically, SN 2015bo is a Cool SN in the Branch classification scheme. The velocity evolution measured from spectral features is consistent with 1991bg-like SNe. SN 2015bo has a SN twin (similar spectra) and sibling (same host galaxy), SN 1997cn. Distance modulii of $\mu$ = 34.36 $\pm$ 0.01 (stat) $\pm$0.13 (sys) mag and $\mu$ = 34.37 $\pm$ 0.04 (stat) $\pm$ 0.12 (sys) mag were derived for SN 2015bo and SN 1997cn, respectively. These distances are consistent at the 0.06-$\sigma$ level with each other, and are also consistent with distances derived using surface-brightness fluctuations and redshift-corrected cosmology. This suggests that fast-declining SNe could be accurate distance indicators which should not be excluded from future cosmological analyses.

Observations of massive galaxies at low redshift have revealed approximately linear scaling relations between the mass of a supermassive black hole (SMBH) and properties of its host galaxy. How these scaling relations evolve with redshift and whether they extend to lower-mass galaxies however remain open questions. Recent galaxy formation simulations predict a delayed, or "two-phase", growth of SMBHs: slow, highly intermittent BH growth due to repeated gas ejection by stellar feedback in low-mass galaxies, followed by more sustained gas accretion that eventually brings BHs onto the local scaling relations. The predicted two-phase growth implies a steep increase, or "kink", in BH-galaxy scaling relations at a stellar mass $M_{*}\sim5\times10^{10} M_{\odot}$. We develop a parametric, semi-analytic model to compare different SMBH growth models against observations of the quasar luminosity function (QLF) at $z\sim0.5-4$. We compare models in which the relation between SMBH mass and galaxy mass is purely linear versus two-phase models. The models are anchored to the observed galaxy stellar mass function, and the BH mass functions at different redshifts are consistently connected by the accretion rates contributing to the QLF. The best fits suggest that two-phase evolution is significantly preferred by the QLF data over a purely linear scaling relation. Moreover, when the model parameters are left free, the two-phase model fits imply a transition mass consistent with that predicted by simulations. Our analysis motivates further observational tests, including measurements of BH masses and AGN activity at the low-mass end, which could more directly test two-phase SMBH growth.

We present an analytic model for clustered supernovae (SNe) feedback in galaxy disks, incorporating the dynamical evolution of superbubbles formed from spatially overlapping SNe remnants. We propose two realistic outcomes for the evolution of superbubbles in galactic disks: (1) the expansion velocity of the shock front falls below the turbulent velocity dispersion of the ISM in the galaxy disk, whereupon the superbubble stalls and fragments, depositing its momentum entirely within the galaxy disk, or (2) the superbubble grows in size to reach the gas scale height, breaking out of the galaxy disk and driving galactic outflows/fountains. In either case, we find that superbubble breakup/breakout almost always occurs before the last Type-II SN ($\lesssim$40 Myr) in the recently formed star cluster, assuming a standard high-end IMF slope, and scalings between stellar lifetimes and masses. The threshold between these two cases implies a break in the effective strength of feedback in driving turbulence within galaxies, and a resulting change in the scalings of, e.g., star formation rates with gas surface density (the Kennicutt-Schmidt relation) and the star formation efficiency in galaxy disks.

Most major scientific results produced by ground-based gamma-ray telescopes in the last 30 years have been obtained by expert members of the collaborations operating these instruments. This is due to the proprietary data and software policies adopted by these collaborations. However, the advent of the next generation of telescopes and their operation as observatories open to the astronomical community, along with a generally increasing demand for open science, confront gamma-ray astronomers with the challenge of sharing their data and analysis tools. As a consequence, in the last few years, the development of open-source science tools has progressed in parallel with the endeavour to define a standardised data format for astronomical gamma-ray data. The latter constitutes the main topic of this review. Common data specifications provide equally important benefits to the current and future generation of gamma-ray instruments: they allow the data from different instruments, including legacy data from decommissioned telescopes, to be easily combined and analysed within the same software framework. In addition, standardised data accessible to the public, and analysable with open-source software, grant fully-reproducible results. In this article we provide an overview of the evolution of the data format for gamma-ray astronomical data, focusing on its progression from private and diverse specifications to prototypical open and standardised ones. The latter have already been successfully employed in a number of publications paving the way to the analysis of data from the next generation of gamma-ray instruments, and to an open and reproducible way of conducting gamma-ray astronomy.

This work communicates the discovery of a binary open cluster within the Galaxy. NGC 1605 presents an unusual morphology with a sparse stellar distribution and a double core in close angular proximity. The 2MASS and Gaia-EDR3 field-star decontaminated colour-magnitude diagrams (CMDs) show two distinct stellar populations located at the same heliocentric distance of $\sim2.6$ kpc suggesting that there are two clusters in the region, NGC 1605a and NGC 1605b, with ages of $2$ Gyr and $600$ Myr, respectively. Both Gaia parallax and PM distributions are compact and very similar indicating that they are open clusters (OCs) and share the same kinematics. The large age difference, 1.4 Gyr, points to a formation by tidal capture during a close encounter and the close spatial proximity and similar kinematics suggest an ongoing merger event. There are some prominent tidal debris that appear to trace the clusters' orbits during the close encounter and, unexpectedly, some of them appear to be bound structures, which may suggest that additionaly to the evaporation the merging clusters are being broken apart into smaller structures by the combination of Galactic disk, Perseus arm, and mutual tidal interactions. In this sense, the newly found binary cluster may be a key object on the observational validation of theoretical studies on binary cluster pairs formation by tidal capture as well as in the formation of massive clusters by merging, and tidal disruption of stellar systems.

Most Ultraluminous X-ray sources (ULXs) are thought to be powered by super-Eddington accretion onto stellar-mass compact objects. Accretors in this extreme regime are naturally expected to ionise copious amounts of plasma in their vicinity and launch powerful radiation-driven outflows from their discs. High spectral resolution X-ray observations (with RGS gratings onboard XMM-Newton) of a few ULXs with the best datasets indeed found complex line spectra and confirmed such extreme (0.1-0.3c) winds. However, a search for plasma signatures in a large ULX sample with a rigorous technique has never been performed, thereby preventing us from understanding their statistical properties such as the rate of occurrence, to constrain the outflow geometry and its duty cycle. We developed a fast method for automated line detection in X-ray spectra and applied it to the full RGS ULX archive, rigorously quantifying the statistical significance of any candidate lines. Collecting the 135 most significant features detected in 89 observations of 19 objects, we created the first catalogue of spectral lines detected in soft X-ray ULX spectra. We found that the detected emission lines are concentrated around known rest-frame elemental transitions and thus originate from low-velocity material. The absorption lines instead avoid these transitions, suggesting they were imprinted by blueshifted outflows. Such winds therefore appear common among the ULX population. Additionally, we found that spectrally hard ULXs show fewer line detections than soft ULXs, indicating some difference in their accretion geometry and orientation, possibly causing over-ionisation of plasma by the harder spectral energy distributions of harder ULXs.

We propose that the origin of faint, broad emission-line wings in the Green Pea (GP) analog Mrk 71 is a clumpy, LyC and/or Ly$\alpha$-driven superwind. Our spatially-resolved analysis of Gemini-N/GMOS-IFU observations shows that these line wings with terminal velocity $>3000~\rm{km~s^{-1}}$ originate from the super star cluster (SSC) Knot A, and propagate to large radii. The object's observed ionization parameter and stellar surface density are close to their theoretical maxima, and radiation pressure dominates over gas pressure. Together with a lack of evidence for supernova feedback, these imply a radiation-dominated environment. We demonstrate that a clumpy, radiation-driven superwind from Knot A is a viable model for generating the extreme velocities, and in particular, that Lyman continuum and/or Ly$\alpha$ opacity must be responsible. We find that the Mrk 71 broad wings are best fitted with power laws, as are those of a representative extreme GP and a luminous blue variable star, albeit with different slopes. This suggests that they may share a common wind-acceleration mechanism. We propose that high-velocity, power-law wings may be a distinctive signature of radiation feedback, and of radiatively-driven winds, in particular.

Type-Ia supernovae (SN Ia) are powerful stellar explosions that provide important distance indicators in cosmology. Recently, we proposed a new SN Ia mechanism that involves a nuclear fission chain-reaction in an isolated white dwarf [PRL 126, 1311010]. Here we perform novel reaction network simulations of the actinide-rich first solids in a cooling white dwarf. The network includes neutron-capture and fission reactions on a range of U and Th isotopes with various possible values for U-235 enrichment. We find, for modest U-235 enrichments, neutron-capture on U-238 and Th-232 can breed additional fissile nuclei so that a significant fraction of all U and Th nuclei may fission during the chain-reaction. The resulting large energy release could ignite thermonuclear carbon burning and possibly trigger a SN Ia.

We present a simple, unified model that can explain two of the brightest, large-scale, diffuse, polarized radio features in the sky, the North Polar Spur (NPS) and the Fan Region, along with several other prominent loops. We suggest that they are long, magnetized, and parallel filamentary structures that surround the Local arm and/or Local Bubble, in which the Sun is embedded. We show this model is consistent with the large number of observational studies on these regions, and is able to resolve an apparent contradiction in the literature that suggests the high latitude portion of the NPS is nearby, while lower latitude portions are more distant. Understanding the contributions of this local emission is critical to developing a complete model of the Galactic magnetic field. These very nearby structures also provide context to help understand similar non-thermal, filamentary structures that are increasingly being observed with modern radio telescopes.

We analyze the distribution of rest-frame U-V and V-J colors for star-forming galaxies at 0.5 < z < 2.5. Using stellar population synthesis, stochastic star formation histories, and a simple prescription for the dust attenuation that accounts for the shape and inclination of galaxies, we construct a model for the distribution of galaxy colors. With only two free parameters, this model is able to reproduce the observed galaxy colors as a function of redshift and stellar mass remarkably well. Our analysis suggests that the wide range of dust attenuation values measured for star-forming galaxies at a given redshift and stellar mass is almost entirely due to the effect of inclination; if all galaxies were observed edge-on, they would show very similar dust attenuation. This result has important implications for the interpretation of dust attenuation measurements, the treatment of UV and IR luminosity, and the comparison between numerical simulations and observations.

Small-scale brightenings in solar atmospheric observations are a manifestation of heating and/or energy transport events. We present statistical characteristics of brightenings from a new detection method applied to 1330, 1400, and 2796 \AA\ IRIS slitjaw image time series. 2377 events are recorded which coexist in all three channels, giving high confidence that they are real. $\approx$1800 of these are spatially coherent, equating to event densities of $\sim9.7\times10^{-5}$arcsec$^{-2}$s$^{-1}$ within a $90\arcsec\times100\arcsec$ FOV over 34.5 minutes. Power Law indices estimates are determined for total brightness ($2.78<\alpha<3.71$), maximum brightness ($3.84<\alpha<4.70$), and average area ($4.31<\alpha<5.70$) distributions. Duration and speed distributions do not obey a power law. A correlation is found between the events' spatial fragmentation, area, and duration, and a weak relationship with total brightness, showing that larger/longer-lasting events are more likely to fragment during their lifetime. Speed distributions show that all events are in motion, with an average speed of $\sim7$\kms. The events' spatial trajectories suggest that cooler 2796 \AA\ events tend to appear slightly later, and occupy a different position/trajectory to the hotter channel results. This suggests that either many of these are impulsive events caused by reconnection, with subsequent rapid cooling, or that the triggering event occurs near the TR, with a subsequent propagating disturbance to cooler atmospheric layers. The spatial distribution of events is not uniform, with broad regions devoid of events. A comparison of spatial distribution with properties of other atmospheric layers shows a tentative connection between high magnetic field strength, the corona's multithermality, and high IRIS brightening activity.

GW200115 was the second merger of a black hole and a neutron star confidently detected through gravitational waves. Inference on the signal allows for a large black hole spin misaligned with the orbital angular momentum, but shows little support for aligned spin values. We show that this is a natural consequence of measuring the parameters of a black hole -- neutron star binary with non-spinning components while assuming the priors used in the LIGO-Virgo-KAGRA analysis. We suggest that, a priori, a non-spinning binary is more consistent with current astrophysical understanding.

The Space-Weather ANalytics for Solar Flares (SWAN-SF) is a multivariate time series benchmark dataset recently created to serve the heliophysics community as a testbed for solar flare forecasting models. SWAN-SF contains 54 unique features, with 24 quantitative features computed from the photospheric magnetic field maps of active regions, describing their precedent flare activity. In this study, for the first time, we systematically attacked the problem of quantifying the relevance of these features to the ambitious task of flare forecasting. We implemented an end-to-end pipeline for preprocessing, feature selection, and evaluation phases. We incorporated 24 Feature Subset Selection (FSS) algorithms, including multivariate and univariate, supervised and unsupervised, wrappers and filters. We methodologically compared the results of different FSS algorithms, both on the multivariate time series and vectorized formats, and tested their correlation and reliability, to the extent possible, by using the selected features for flare forecasting on unseen data, in univariate and multivariate fashions. We concluded our investigation with a report of the best FSS methods in terms of their top-k features, and the analysis of the findings. We wish the reproducibility of our study and the availability of the data allow the future attempts be comparable with our findings and themselves.

Measurements of cosmic flows enable us to test whether cosmological models can accurately describe the evolution of the density field in the nearby Universe. In this paper, we measure the low-order kinematic moments of the cosmic flow field, namely bulk flow and shear moments, using the Cosmicflows-4 Tully-Fisher catalogue (CF4TF). To make accurate cosmological inferences with the CF4TF sample, it is important to make realistic mock catalogues. We present the mock sampling algorithm of CF4TF. These mock can accurately realize the survey geometry and luminosity selection function, enabling researchers to explore how these systematics affect the measurements. These mocks can also be further used to estimate the covariance matrix and errors of power spectrum and two-point correlation function in future work. In this paper, we use the mocks to test the cosmic flow estimator and find that the measurements are unbiased. The measured bulk flow in the local Universe is 376 $\pm$ 23 (error) $\pm$ 183 (cosmic variance) km s$^{-1}$ at depth $d_{\text{MLE}}=35$ Mpc $h^{-1}$, to the Galactic direction of $(l,b)=(298\pm 3^{\circ}, -6\pm 3^{\circ})$. Both the measured bulk and shear moments are consistent with the concordance $\Lambda$ Cold Dark Matter cosmological model predictions.

Helium star - carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code $\texttt{MESA}$, we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the time of explosion. We convert these model properties to NUV/optical magnitudes assuming a blackbody spectrum and support this approach using a matched stellar atmosphere model. These models will be valuable to compare with pre-explosion imaging for future supernovae, though we emphasize the observational difficulty of detecting extremely blue companions. The pre-explosion source detected in association with SN 2012Z has been interpreted as a helium star binary containing an initially ultra-massive WD in a multi-day orbit. However, extending our binary models to initial CO WD masses of up to $1.2\,M_{\odot}$, we find that these systems undergo off-center carbon ignitions and thus are not expected to produce thermonuclear supernovae. This tension suggests that, if SN 2012Z is associated with a helium star - WD binary, then the pre-explosion optical light from the system must be significantly modified by the binary environment and/or the WD does not have a carbon-rich interior composition.

We report results from new and archival observations of the newly discovered active asteroid (248370) 2005 QN_137, which has been determined to be a likely main-belt comet based on a subsequent discovery that it is recurrently active near perihelion. From archival data analysis, we estimate g'-, r'-, i'-, and z'-band absolute magnitudes for the nucleus of H_g=16.62+/-0.13, H_r=16.12+/-0.10, H_i=16.05+/-0.11, and H_z=15.93+/-0.08, corresponding to nucleus colors of g'-r'=0.50+/-0.16, r'-i'=0.07+/-0.15, and i'-z'=0.12+/-0.14, an equivalent V-band absolute magnitude of H_V=16.32+/-0.08, and a nucleus radius of r_n=1.6+/-0.2 km (using a V-band albedo of p_V=0.054+/-0.012). Meanwhile, we find mean near-nucleus coma colors when 248370 was active of g'-r'=0.47+/-0.03, r'-i'=0.10+/-0.04, and i'-z'=0.05+/-0.05, and similar mean dust tail colors, suggesting that no significant gas coma is present. We find approximate ratios between the scattering cross-sections of near-nucleus dust (within 5000 km of the nucleus) and the nucleus of A_d/A_n=0.7+/-0.3 on 2016 July 22, and 1.8<A_d/A_n<2.9 in 2021 July and August. During the 2021 observation period, the coma declined in intrinsic brightness by ~0.35 mag (or ~25%) in 37 days, while the surface brightness of the dust tail remained effectively constant over the same period. Constraints derived from the sunward extent of the coma suggest that terminal velocities of ejected dust grains are extremely slow (~1 m/s for 1 micron particles), indicating that the observed dust emission may have been aided by rapid rotation of the nucleus lowering the effective escape velocity.

It is noted that effective very recently, the orbital motion of Encke's comet has been found to be affected by a very slight nongravitational deceleration. Soon after J. F. Encke established in the early 19th century that the comet was returning to perihelion every 3.3 years, he also discovered that the object was notorious for returning to perihelion a little earlier than predicted by the Newtonian theory. The acceleration persisted over a period of two centuries, but its rate was gradually decreasing, Generations of cometary astronomers were curious to know whether or not the comet would eventually move in purely gravitational orbit. A model based on the assumption of a precession of the comet's nucleus, which predicted that the acceleration would change into a deceleration, was not published until 1979. This transition has now been documented by two independent, highly-accurate orbit determinations. The era of the comet's persevering acceleration is finally over.

We study the formation of primordial black hole (PBH) dark matter and the generation of scalar induced secondary gravitational waves (SIGWs) in a non-supersymmetric model of hybrid inflation with chaotic (polynomial-like) potential, including one-loop radiative corrections. A radiatively corrected version of these models is entirely consistent with Planck's data. By adding non-canonical kinetic energy term in the lagrangian, the inflaton experiences a period of ultra-slow-roll, and the amplitude of primordial power spectrum is enhanced to $O(10^{-2})$. The enhanced power spectra of primordial curvature perturbations can have both sharp and broad peaks. A wide mass range of PBH is realized in our models, and the frequencies of scalar induced gravitational waves are ranged from nHz to Hz. We present several benchmark points where the PBH mass generated during inflation is around $(1 - 100) \, M_{\odot}$, $(10^{-9} - 10^{-7}) \, M_{\odot}$ and $(10^{-16} - 10^{-11}) \, M_{\odot}$. The PBHs can make up most of the dark matter with masses around $(10^{-16} - 10^{-11}) \, M_{\odot}$ and $(1 - 100) \, M_{\odot}$, and their associated SIGWs can be probed by the upcoming space-based gravitational wave (GW) observatories. The formation of $~10 \, M_{\odot}$ PBH with wide peaks of SIGWs can be used to interpret the stochastic GW signal in the nHz band, detected by the North American Nano hertz Observatory (NANO-Grav) for Gravitational Waves. These signals may also be tested by future interferometer-type GW observations of EPTA, SKA, LISA, TaiJi, TianQin and Einstein Telescope (ET).

We study the stellar binary black holes (BBHs) inspiralling/merging in galactic nuclei based on our numerical method GNC. We find that $3-40\%$ of all new born BBHs will finally merge due to various dynamical effects. In a five year's mission, up to $10^4$, $10^5$, $\sim100$ of BBHs inspiralling/merging in galactic nuclei can be detected with SNR$>8$ in aLIGO, Einstein/DECIGO, TianQin/LISA/TaiJi, respectively. About tens are detectable in both LISA/TaiJi/TianQin and aLIGO. These BBHs have two unique characteristics: (1) Significant eccentricities. $1-3\%$, $2-7\%$, or $30-90\%$ of them is with $e_i>0.1$ when they enter into aLIGO, Einstein, or space observatories, respectively. Such high eccentricities provide a possible explanation for that of GW 190521. Most highly-eccentric BBHs are not detectable in LISA/Tianqin/TaiJi before entering into aLIGO/Einstein as their strain become significant only at $f_{\rm GW}\gtrsim0.1$ Hz. DECIGO become an ideal observatory to detect those events as it can fully cover the rising phase. (2) Up to $2\%$ of BBHs can inspiral/merge at distances $\lesssim10^3 r_{\rm SW}$ from the massive black hole (MBH), with significant accelerations, such that the Doppler phase drift of $\sim10-10^5$ of them can be detectable with SNR$>8$ in space observatories. The energy density of the gravitational wave backgrounds (GWB) contributed by these BBHs deviate from the powerlaw slope of $2/3$ at $f_{\rm GW}\lesssim 1$mHz. The high eccentricity, significant accelerations and different profile of GWB of these sources make them distinguishable, thus interesting for future GW detections and tests of relativities.

We provide a model-independent reconstruction of dark energy from $z=0$ to $ \gtrsim 10^5$. We parameterise the model by a perfect fluid with a series of physically well-motivated bins in energy-density, such that the equation of state is always $-1 \le w \le 1$. Our method is capable of describing a range of theoretical models with smooth modifications to the expansion history. Combining the latest CMB, BAO, SN and local $H_0$ measurements, we obtain a large improvement of $\Delta \chi^2=41.3$ over LCDM, at the expense of 33 additional parameters in the fit, with dark energy contributing significantly between $z \sim 10^4 - 10^5$, and intriguingly with a sound speed $c_s^2 \sim 1/3$. A significant part of the $\Delta \chi^2$ improvement comes from \Planck\ + Atacama Cosmology Telescope (\textsc{Act}) data, alleviating tension between them within LCDM. We apply a correlation prior to penalise models with unnecessary degrees of freedom and find no preference for deviations from LCDM at late-times, but moderate Bayesian evidence of an early dark energy (EDE) component. Although the model has a large amount of freedom, it is unable to reduce $S_8 \equiv \sigma_8 (\Omega_\mathrm{c} / 0.3)^{0.5}$ below that of \lcdm, to bring about full concordance with large-scale structure data.

We present the first sub-mJy ($\approx0.7$ mJy beam$^{-1}$) survey to be completed below 100 MHz, which is over an order of magnitude deeper than previously achieved for widefield imaging of any field at these low frequencies. The high resolution ($15 \times 15$ arcsec) image of the Bo\"otes field at 34-75 MHz is made from 56 hours of observation with the LOw Frequency ARray (LOFAR) Low Band Antenna (LBA) system. The observations and data reduction, including direction-dependent calibration, are described here. We present a radio source catalogue containing 1,948 sources detected over an area of $23.6$ deg$^2$, with a peak flux density threshold of $5\sigma$. Using existing datasets, we characterise the astrometric and flux density uncertainties, finding a positional uncertainty of $\sim1.2$ arcsec and a flux density scale uncertainty of about 5 per cent. Using the available deep 144-MHz data, we identified 144-MHz counterparts to all the 54-MHz sources, and produced a matched catalogue within the deep optical coverage area containing 829 sources. We calculate the Euclidean-normalised differential source counts and investigate the low-frequency radio source spectral indices between 54 and 144 MHz, both of which show a general flattening in the radio spectral indices for lower flux density sources, from $\sim-0.75$ at 144-MHz flux densities between 100-1000 mJy to $\sim-0.5$ at 144-MHz flux densities between 5-10 mJy, due to a growing population of star forming galaxies and compact core-dominated AGN.

We examine the possibility for the existence of gravitomagnetic monopole ($n_*$) in M87* by using the results obtained from its first Event Horizon Telescope image. By numerically deducing the shadow sizes in Kerr-Taub-NUT (KTN) spacetime, we show that the shadow size increases with increasing $|n_*|$ for a fixed Kerr parameter $|a_*|$ in case of the KTN black hole, whereas for a KTN naked singularity it increases with increasing $n_*$ for a fixed $a_* > 0$ if $n_* > -\cot 17^{\circ}$ . In general, the asymmetry of shadow shape increases if the central dark object in M87 is a KTN/Kerr naked singularity instead of a KTN/Kerr black hole. We find that a non-zero gravitomagnetic monopole is still compatible with the current EHT observations, in which case the upper limit of $n_*$ cannot be greater than $1.1$, i.e., $n_* \lesssim 1.1$ for the prograde rotation ($a_* > 0$), and the lower limit of $n_*$ cannot be less than $-1.1$, i.e., $ n_* \gtrsim -1.1$ for the retrograde rotation ($a_* < 0$). Moreover, if the circularity of the shadow can be measured on a precision of $\lesssim 1\%$, the Kerr and KTN naked singularities can be falsified for M87*.

The Fourier transformation is an effective and efficient operation of Gaussianization at the one-point level. Using a set of N-body simulation data, we verified that the one-point distribution functions of the dark matter momentum divergence and density fields closely follow complex Gaussian distributions. Statistical theories about the one-point distribution function of the quotient of two complex Gaussian random variables are introduced and applied to model one-point statistics about the growth of individual Fourier mode of the dark matter density field, namely the mode-dependent growth function and the mode growth rate, which can be obtained by the ratio of two Fourier transformed cosmic fields. Our simulation results proved that the models based on the Gaussian approximation are impressively accurate, and our analysis revealed many interesting aspects about the growth of dark matter's density fluctuation in Fourier space.

The evaporations of Primordial Black Holes (PBH) (via Hawking radiation) can produce electrons/positrons ($e^-/e^+$) in the Galactic Centre (GC) region which under the influence of the magnetic field of Centre region can emit synchrotron radiation. These $e^-/e^+$ can also induce Inverse Compton radiation due to the scattering with ambient photons. In this work three different PBH mass distributions namely, monochromatic, power law and lognormal distributions are considered to calculate such radiation fluxes. On the other hand, annihilation or decay of dark matter in the Galactic Centre region can also yield $e^-/e^+$ as the end product which again may emit synchrotron radiation in the Galactic magnetic field and also induce Inverse Compton scattering. In this work a comparative study is made for these radiation fluxes from both PBH evaporations and from dark matter origins and their detectabilities are addressed in various ongoing and other telescopes as well as in upcoming telescopes such as SKA. The variations of these radiation fluxes with the distance from the Galactic Centre are also computed and it is found that such variations could be a useful probe to determine the mass of PBH or the mass of dark matter.

In the field of multi-messenger astronomy, Bayesian inference is commonly adopted to compare the compatibility of models given the observed data. However, to describe a physical system like neutron star mergers and their associated gamma-ray burst (GRB) events, usually more than ten physical parameters are incorporated in the model. With such a complex model, likelihood evaluation for each Monte Carlo sampling point becomes a massive task and requires a significant amount of computational power. In this work, we perform quick parameter estimation on simulated GRB X-ray light curves using an interpolated physical GRB model. This is achieved by generating a grid of GRB afterglow light curves across the parameter space and replacing the likelihood with a simple interpolation function in the high-dimensional grid that stores all light curves. This framework, compared to the original method, leads to a $\sim$90$\times$ speedup per likelihood estimation. It will allow us to explore different jet models and enable fast model comparison in the future.

Context. The magnetic Kelvin-Helmholtz instability (KHI) has been proposed as a means of generating magnetohy- drodynamic turbulence and encouraging wave energy dissipation in the solar corona, particularly within transversely oscillating loops. Aims. Our goal is to determine whether the KHI encourages magnetic reconnection in oscillating flux tubes in the solar corona. This will establish whether the instability enhances the dissipation rate of energy stored in the magnetic field. Methods. We conducted a series of three-dimensional magnetohydrodynamic simulations of the KHI excited by an oscillating velocity shear. We investigated the effects of numerical resolution, field line length, and background currents on the growth rate of the KHI and on the subsequent rate of magnetic reconnection. Results. The KHI is able to trigger magnetic reconnection in all cases, with the highest rates occurring during the initial growth phase. Reconnection is found to occur preferentially along the boundaries of Kelvin-Helmholtz vortices, where the shear in the velocity and magnetic fields is greatest. The estimated rate of reconnection is found to be lowest in simulations where the KHI growth rate is reduced. For example, this is the case for shorter field lines or due to shear in the background field. Conclusions. In non-ideal regimes, the onset of the instability causes the local reconnection of magnetic field lines and enhances the rate of coronal wave heating. However, we found that if the equilibrium magnetic field is sheared across the Kelvin-Helmholtz mixing layer, the instability does not significantly enhance the rate of reconnection of the background field, despite the free energy associated with the non-potential field.

Intensity bursts in ultraviolet (UV) to X-ray wavelengths, and plasma jets are typical signatures of magnetic reconnection and the associated impulsive heating of the solar atmospheric plasma. To gain new insights into the process, high-cadence observations are required to capture the rapid response of plasma to magnetic reconnection as well as the highly dynamic evolution of jets. Here we report 2\,s cadence extreme-UV observations recorded by the 174\,\AA\ High Resolution Imager of the Extreme Ultraviolet Imager onboard the Solar Orbiter mission. These observations, covering a quiet-Sun coronal region, reveal the onset signatures of magnetic reconnection as localized heating events. These localized sources then exhibit repeated plasma eruptions or jet activity. Our observations show that this spatial morphological change from localized sources to jet activity could occur rapidly on timescales of about 20\,s. The jets themselves are intermittent and are produced from the source region on timescales of about 20\,s. In the initial phases of these events, plasma jets are observed to exhibit speeds, as inferred from propagating intensity disturbances, in the range of 100\,km\,s$^{-1}$ to 150\,km\,s$^{-1}$. These jets then propagate to lengths of about 5\,Mm. We discuss examples of bidirectional and unidirectional jet activity observed to be initiated from the initially localized bursts in the corona. The transient nature of coronal bursts and plasma jets/flows and their dynamics could provide a benchmark for magnetic reconnection models of coronal bursts and jets.

Strange quark matter (SQM) may be the true ground state of matter. According to this SQM hypothesis, the observed neutron stars actually should all be strange quark stars. But distinguishing between neutron stars and strange quark stars by means of obser- vations is extremely difficult. It is interesting to note that under the SQM hypothesis, less massive objects such as strange quark planets and strange dwarfs can also stably exist. The extremely high density and small radius of strange quark planets give us some new perspectives to identify SQM objects and to test the SQM hypothesis. First, the tidal disruption radius of strange quark planets is much smaller than normal planets, so, very close-in exoplanets can be safely identified as candidates of SQM objects. Sec- ond, gravitational waves (GW) from mergers of strange quark star-strange quark planet systems are strong enough to be detected by ground-based GW detectors. As a result, GW observation will be a powerful tool to probe SQM stars. At the same time, the tidal deformability of SQM planets can be measured to further strengthen the result.

Filaments play an important role in star formation, but the formation process of filaments themselves is still unclear. The high-mass star forming clump G286.21+0.17 (G286 for short) that contains an "L" type filament was thought to undergo global collapse. Our high resolution ALMA band 3 observations resolve the gas kinematics of G286 and reveal two sub-clumps with very different velocities inside it. We find that the "blue profile" (an indicator of gas infall) of HCO+ lines in single dish observations of G286 is actually caused by gas emission from the two sub-clumps rather than gas infall. We advise great caution in interpreting gas kinematics (e.g., infall) from line profiles toward distant massive clumps in single dish observations. Energetic outflows are identified in G286 but the outflows are not strong enough to drive expansion of the two sub-clumps. The two parts of the "L" type filament ("NW-SE" and "NE-SW" filaments) show prominent velocity gradients perpendicular to their major axes, indicating that they are likely formed due to large-scale compression flows. We argue that the large-scale compression flows could be induced by the expansion of nearby giant HII regions. The "NW-SE" and "NE-SW" filaments seem to be in collision, and a large amount of gas has been accumulated in the junction region where the most massive core G286c1 forms.

The existence of cosmic accelerators able to emit charged particles up to EeV energies has been confirmed by the observations made in the last years by experiments such as Auger and Telescope Array. The interaction of such energetic cosmic-rays with gas or low energy photons, surrounding the astrophysical sources or present in the intergalactic medium, guarantee an ultra-high-energy neutrino related emission. When these energetic neutrinos interact in a medium produce a thermo-acoustic process where the energy of generated particle cascades can be conveyed in a pressure pulse propagating into the same medium. The kilometric attenuation length as well as the well-defined shape of the expected pulse suggest a large-area-undersea-array of acoustic sensors as an ideal observatory. For this scope, we propose the recycle of ENI offshore (oil rigs) powered platforms in the Adriatic sea as the main infrastructure to build an acoustic submarine array of dedicated hydrophones. In this work we describe the advantages of this detector concept using a new wave-tracing technique as well as the scientific goals linked to the challenging purpose of observing for the first time ultra-high-energy cosmic neutrinos. This observatory will be complementary to the dedicated radio array detectors with the advantages of avoiding any possible thermo-acoustic noise from the atmospheric muons.

In this contribution we present the first numerical simulations of a relativistic outflow propagating through the inner hundreds of parsecs of its host galaxy, including atomic and ionised hydrogen, and the cooling effects of ionisation. Our results are preliminary, but we observe efficient shock ionization of atomic hydrogen in interstellar clouds. The mean density of the interstellar medium in these initial simulations is lower than that expected in typical galaxies, which makes cooling times longer and thus no recombination is observed inside the shocked region. The velocities achieved by the shocked gas in the simulations are in agreement with observational results, although with a wide spectrum of values.

We present $\texttt{matryoshka}$, a suite of neural network based emulators and accompanying Python package that have been developed with the goal of producing fast and accurate predictions of the nonlinear galaxy power spectrum. The suite of emulators consists of four linear component emulators, from which fast linear predictions of the power spectrum can be made, allowing all nonlinearities to be included in predictions from a nonlinear boost component emulator. The linear component emulators includes an emulator for the matter transfer function that produces predictions in $\sim 0.0004 \ \mathrm{s}$, with an error of $<0.08\%$ (at $1\sigma$ level) on scales $10^{-4} \ h \ \mathrm{Mpc}^{-1}<k<10^1 \ h \ \mathrm{Mpc}^{-1}$. In this paper we demonstrate $\texttt{matryoshka}$ by training the nonlinear boost component emulator with analytic training data calculated with HALOFIT, that has been designed to replicate training data that would be generated using numerical simulations. Combining all the component emulator predictions we achieve an accuracy of $< 0.75\%$ (at $1\sigma$ level) when predicting the real space nonlinear galaxy power spectrum on scales $0.0025 \ h \ \mathrm{Mpc}^{-1}<k<1 \ h \ \mathrm{Mpc}^{-1}$. We use $\texttt{matryoshka}$ to investigate the impact of the analysis setup on cosmological constraints by conducting several full shape analyses of the real space galaxy power spectrum. Specifically we investigate the impact of the minimum scale (or $k_\mathrm{max}$), finding an improvement of $\sim 1.8\times$ in the constraint on $\sigma_8$ by pushing $k_\mathrm{max}$ from $k_\mathrm{max}=0.25 \ h \ \mathrm{Mpc}^{-1}$ to $k_\mathrm{max}=0.85 \ h \ \mathrm{Mpc}^{-1}$, highlighting the potential gains when using clustering emulators such as $\texttt{matryoshka}$ in cosmological analyses. $\texttt{matryoshka}$ is publicly available at https://github.com/JDonaldM/Matryoshka.

We present wideband polarimetric observations, obtained with the Karl G. Jansky Very Large Array (VLA), of the merging galaxy cluster MACS J0717.5+3745, which hosts one of the most complex known radio relic and halo systems. We use both Rotation Measure Synthesis and QU-fitting, and find a reasonable agreement of the results obtained with these methods, in particular, when the Faraday distribution is simple and the depolarization is mild. The relic is highly polarized over its entire length reaching a fractional polarization ${>}30\%$ in some regions. We also observe a strong wavelength-dependent depolarization for some regions of the relic. The northern part of the relic shows a complex Faraday distribution suggesting that this region is located in or behind the intracluster medium (ICM). Conversely, the southern part of the relic shows a Rotation Measure very close to the Galactic foreground, with a rather low Faraday dispersion, indicating very little magnetoionic material intervening the line-of-sight. From spatially resolved polarization analysis, we find that the scatter of Faraday depths correlates with the depolarization, indicating that the tangled magnetic field in the ICM causes the depolarization. At the position of a well known narrow-angle-tailed galaxy (NAT), we find evidence of two components clearly separated in Faraday space. The high Faraday dispersion component seems to be associated with the NAT, suggesting the NAT is embedded in the ICM while the southern part of the relic lies in front of it. The magnetic field orientation follows the relic structure indicating a well-ordered magnetic field. We also detect polarized emission in the halo region; however the absence of significant Faraday rotation and a low value of Faraday dispersion suggests the polarized emission, previously considered as the part of the halo, has a shock(s) origin.

The Cygnus Loop (G74.0-8.5) is a very well-known nearby supernova remnant (SNR) in our Galaxy. Thanks to its large size, brightness, and angular offset from the Galactic plane, it has been studied in detail from radio to $\gamma$-ray emission. The $\gamma$ -rays probe the populations of energetic particles and their acceleration mechanisms at low shock speeds. We present an analysis of the $\gamma$-ray emission detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope over 11 years in the region of the Cygnus Loop. We performed detailed morphological and spectral studies of the $\gamma$-ray emission toward the remnant from 100 MeV to 100 GeV and compared it with X-ray, UV, optical, and radio images. The higher statistics with respect to the previous studies enabled us to decompose the emission from the remnant into two morphological components to model its nonthermal multiwavelength emission. The extended $\gamma$-ray emission is well correlated with the thermal X-ray and UV emission of the SNR. Our morphological analysis reveals that a model considering two contributions from the X-ray and the UV emission regions is the best description of the $\gamma$-ray data. Both components show a curved spectrum, but the X-ray component is softer and more curved than the UV component, suggesting a different physical origin. The multiwavelength modeling of emission toward the SNR suggests that the nonthermal radio and $\gamma$-ray emission associated with the UV component is mostly due to the reacceleration of preexisting cosmic rays by radiative shocks in the adjacent clouds, while the nonthermal emission associated with the X-ray component arises from freshly accelerated cosmic rays.

The 182Hf-182W chronology of iron meteorites provides crucial information on the timescales of accretion and differentiation of some of the oldest planetesimals of the Solar System. Determining accurate Hf-W model ages of iron meteorites requires correction for cosmic ray expo-sure (CRE) induced modifications of W isotope compositions, which can be achieved using in-situ neutron dosimeters such as Pt isotopes. Until now it has been assumed that all Pt isotope variations in meteorites reflect CRE, but here we show that some ungrouped iron meteorites display small nucleosynthetic Pt isotope anomalies. These provide the most appropriate starting composition for the correction of CRE-induced W isotope variations in iron meteorites from all major chemical groups, which leads to a ~1 Ma upward revision of previously reported Hf-W model ages. The revised ages indicate that core formation in non-carbonaceous (NC) iron meteorite parent bodies occurred at ~1-2 Ma after CAI formation, whereas most carbonaceous (CC) iron meteorite parent bodies underwent core formation ~2 Ma later. We show that the younger CC cores have lower Fe/Ni ratios than the earlier-formed NC cores, indicating that core formation under more oxidizing conditions occurred over a more protracted timescale. Thermal modeling of planetesimals heated by 26Al-decay reveals that this protracted core formation timescale is consistent with a higher fraction of water ice in CC compared to NC planetesimals, implying that in spite of distinct core formation timescales, NC and CC iron meteorite parent bodies accreted about contemporaneously within ~1 Ma after CAI formation, but at different radial locations in the disk.

In the solar corona and solar wind, electron heat conduction is an important process that transports energy over large distances and helps determine the spatial variation of temperature. High-density regions undergoing rapid particle-particle collisions exhibit a heat flux described well by classical Spitzer-Harm theory. However, much of the heliosphere is closer to a more collisionless state, and there is no standard description of heat conduction for fluid-based (e.g., magnetohydrodynamic) models that applies generally. Some proposed models rely on electron velocity distributions that exhibit negative values of the phase-space density. In this paper, we explore how positive-definite velocity distributions can be used in fluid-based conservation equations for the electron heat flux along magnetic-field lines in the corona and solar wind. We study both analytic forms of skewed distributions (e.g., skew-normal distributions, two-sided bi-Maxwellians, and constant-collision-time electrostatic solutions) and empirical fits to measurements of core, halo, and strahl electrons in interplanetary space. We also present example solutions to a generalized conservation equation for the heat flux in the solar wind, with some limiting cases found to resemble known free-streaming approximations. The resulting values of the electron heat flux vary as a function of radial distance and Knudsen number in ways that resemble observed data. We note that this model does not include the effects of kinetic instabilities (which may impose saturation limits when active), so for now its regime of applicability is limited to collisionless heat-flux evolution away from the known instability boundaries in parameter space.

We analyze the structure and evolution of ribbons from the M7.3 SOL2014-04-18T13 flare using ultraviolet (UV) images from IRIS and SDO/AIA, magnetic data from SDO/HMI, hard X-ray (HXR) images from RHESSI, and light curves from Fermi/GBM, in order to infer properties of coronal magnetic reconnection. As the event progresses, two flare ribbons spread away from the magnetic polarity inversion line. The width of the newly brightened front along the extension of the ribbon is highly intermittent in both space and time, presumably reflecting non-uniformities in the structure and/or dynamics of the flare current sheet. Furthermore, the ribbon width grows most rapidly in regions exhibiting concentrated non-thermal HXR emission, with sharp increases slightly preceding the HXR bursts. The light curve of the ultraviolet emission matches the HXR light curve at photon energies above 25 keV. In other regions the ribbon-width evolution and light curves do not temporally correlate with the HXR emission. This indicates that the production of non-thermal electrons is highly non-uniform within the flare current sheet. Our results suggest a strong connection between the production of non-thermal electrons and the locally enhanced perpendicular extent of flare ribbon fronts, which in turn reflects inhomogeneous structure and/or reconnection dynamics of the current sheet. Despite this variability, the ribbon fronts remain nearly continuous, quasi-one-dimensional features. Thus, although the reconnecting coronal current sheets are highly structured, they remain quasi-two-dimensional and the magnetic energy release occurs systematically, rather than stochastically, through the volume of reconnecting magnetic flux.

We fabricated a 302 mm diameter low-pass filter made of alumina that has an anti-reflection coating (ARC) made with laser-ablated sub-wavelength structures (SWS). The filter has been integrated into and is operating with the MUSTANG2 instrument, which is coupled to the Green Bank Telescope. The average transmittance of the filter in the MUSTANG2 operating band between 75 and 105 GHz is 98%. Reflective loss due to the ARC is 1%. The difference in transmission between the s- and p-polarization states is less than 1%. To within 1% accuracy we observe no variance in these results when transmission is measured in six independent filter spatial locations. The alumina filter has replaced a prior MUSTANG2 Teflon filter. Data taken with the filter heat sunk to its nominal 40 K stage show performance consistent with expectations: a reduction of about 50% in filters-induced optical power load on the 300 mK stage, and in in-band optical loading on the detectors. This is the first report of an alumina filter with SWS ARC deployed with an operating instrument, and the first demonstration of a large area fabrication of SWS with laser ablation.

A number of popular extensions of the Standard Model of particle physics predict the existence of doubly charged scalar particles $X^{\pm\pm}$. Such particles may be long-lived or even stable. If exist, $X^{--}$ could form atomic bound states with light nuclei and catalyze their fusion by essentially eliminating the Coulomb barrier between them. Such an $X$-catalyzed fusion ($X$CF) process does not require high temperatures or pressure and may have important applications for energy production. A similar process of muon-catalyzed fusion ($\mu$CF) has been shown not to be a viable source of energy because of the sticking of negative muons to helium nuclei produced in the fusion of hydrogen isotopes, which stops the catalytic process. We analyze $X$CF in deuterium environments and show that the $X$-particles can only stick to $^6$Li nuclei, which are produced in the third-stage reactions downstream the catalytic cycle. The corresponding sticking probability is very low, and, before getting bound to $^6$Li, each $X$-particle can catalyze $\sim 3.5\cdot 10^{9}$ fusion cycles, producing $\sim 7\cdot 10^{4}$ TeV of energy. We also discuss the ways of reactivating the $X$-particles from the Coulomb-bound (${\rm ^6Li}X$) states, which would allow re-using them in $X$CF reactions.

If the QCD axion solves the strong CP problem then light axion-like-particles (ALPs) are expected to be ubiquitous in string theory - the string axiverse. Such ALPs can be the QCD axion and constitute dark matter (DM) or radiation, quintessence, and lead to new forces. String ALPs are also expected to give rise to a multiplicity of cosmologically important global axion strings. We study the properties of these axiverse cosmic strings including the vital effects of moduli stabilization, and find that the string cores provide `portals' to different decompactifications - to be precise, the cores explore the large K\"ahler or complex structure boundary of moduli space. As usual for global strings the tension $T_1\sim \Lambda^2 \log(L\Lambda)$ with inter-string separation, $L$, while $\Lambda$ can be small $\ll M_{\rm pl}$. At long distances from the string there are potential new signatures involving variations in Standard Model (SM) parameters (Yukawa couplings, gauge couplings, masses) and equivalence principle violations.

Originally introduced in connection with general relativistic Coriolis forces, the term $\textit{frame-dragging}$ is associated today with a plethora of effects related to the off-diagonal element of the metric tensor. It is also frequently the subject of misconceptions leading to incorrect predictions, even of nonexistent effects. We show that there are three different levels of frame-dragging corresponding to three distinct gravitomagnetic objects: gravitomagnetic potential 1-form, field, and tidal tensor, whose effects are independent, and sometimes opposing. It is seen that, from the two analogies commonly employed, the analogy with magnetism holds strong where it applies, whereas the fluid-dragging analogy (albeit of some use, qualitatively, in the first level) is in general misleading. Common misconceptions (such as viscous-type "body-dragging") are debunked. Applications considered include rotating cylinders (Lewis-Weyl metrics), Kerr, Kerr-Newman and Kerr-dS spacetimes, black holes surrounded by disks/rings, and binary systems.

Recently the Baksan Experiment on Sterile Transitions (BEST) has presented results confirming the gallium anomaly -- lacks of electron neutrinos $\nu_e$ at calibrations of SAGE and GALLEX -- at the statistical significance exceeding 5$\,\sigma$. This result is consistent with explanation of the gallium anomaly as electron neutrino oscillations into sterile neutrino, $\nu_s$. Within this explanation the BEST experiment itself provides the strongest evidence for the sterile neutrino among all the previous anomalous results in the neutrino sector. We combine the results of gallium experiments with searches for sterile neutrinos in reactor antineutrino experiments (assuming CPT-conservation in the $3+1$ neutrino sector). While the 'gallium' best fit point in the model parameter space (sterile neutrino mass squared $m_{\nu_s}^2\approx 1.25\,$eV$^2$, sterile-electron neutrino mixing $\sin^22\theta\approx 0.34$) is excluded by these searches, a part of the BEST-favored 2$\,\sigma$ region with $m^2_{\nu_s}>5\,$eV$^2$ is consistent with all of them. Remarkably, the regions advertised by anomalous results of the NEUTRINO-4 experiment overlap with those of the BEST experiment: the best fit point of the joint analysis is $\sin^22\theta\approx 0.38$, $m_{\nu_s}^2\approx7.3\,$eV$^2$, the favored region will be explored by the KATRIN experiment. The sterile neutrino explanation of the BEST results would suggest not only the extension of the Standard Model of particle physics, but also either serious modifications of the Standard Cosmological Model and Solar Model, or a specific modification of the sterile sector needed to suppress the sterile neutrino production in the early Universe and in the Sun.

We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (tensor transverse, TT) correlations. Specifically, we find ST correlations with a signal-to-noise ratio of 2.8 that are preferred over TT correlations (Hellings and Downs correlations) with Bayesian odds of about 20:1. However, the significance of ST correlations is reduced dramatically when we include modeling of the Solar System ephemeris systematics and/or remove pulsar J0030$+$0451 entirely from consideration. Even taking the nominal signal-to-noise ratios at face value, analyses of simulated data sets show that such values are not extremely unlikely to be observed in cases where only the usual TT modes are present in the GWB. In the absence of a detection of any polarization mode of gravity, we place upper limits on their amplitudes for a spectral index of $\gamma = 5$ and a reference frequency of $f_\text{yr} = 1 \text{yr}^{-1}$. Among the upper limits for eight general families of metric theories of gravity, we find the values of $A^{95\%}_{TT} = (9.7 \pm 0.4)\times 10^{-16}$ and $A^{95\%}_{ST} = (1.4 \pm 0.03)\times 10^{-15}$ for the family of metric spacetime theories that contain both TT and ST modes.

First-order phase transitions exist in many models beyond the Standard Model and can generate detectable stochastic gravitational waves for a strong one. Using the cosmological observables in big bang nucleosynthesis and cosmic microwave background, we derive constraints on the phase transition temperature and strength parameter in a model-independent way. For a strong phase transition, we find that the phase transition temperature should be above around 2 MeV for both reheating photon and neutrino cases. For a weak one with the temperature below 1 MeV, the phase transition strength parameter is constrained to be smaller than around 0.1. Implications for using a first-order phase transition to explain the NANOGrav observed gravitational wave signal are also discussed.

In this work, we discuss the deconfinement phase transition to quark matter in hot/dense matter. We {examine} the effect that different charge fractions, isospin fractions, net strangeness, and chemical equilibrium with respect to leptons have on the position of the coexistence line between different phases. In particular, we investigate how different sets of conditions that describe matter in neutron stars and their mergers, or matter created in heavy-ion collisions affect the position of the critical end point, namely where the first-order phase transition becomes a crossover. We also present an introduction to the topic of critical points, including a review of recent {advances} concerning QCD critical points.

In recent years, stable and unstable manifolds of invariant objects (such as libration points and periodic orbits) have been increasingly recognized as an efficient tool for designing transfer trajectories in space missions. However, most methods currently used in mission design rely on using eigenvectors of the linearized dynamics as local approximations of the manifolds. Since such approximations are not accurate except very close to the base invariant object, this requires large amounts of numerical integration to globalize the manifolds and locate intersections. In this paper, we study hyperbolic resonant periodic orbits in the planar circular restricted 3-body problem, and transfer trajectories between them, by: 1) determining where to search for resonant periodic orbits; 2) developing and implementing a parameterization method for accurate computation of their invariant manifolds as Taylor series; and 3) developing a procedure to compute intersections of the computed stable and unstable manifolds. We develop and implement algorithms that accomplish these three goals, and demonstrate their application to the problem of transferring between resonances in the Jupiter-Europa system.

When the planar circular restricted 3-body problem is periodically perturbed, most unstable periodic orbits become invariant tori. However, 2D Poincar\'e sections no longer work to find their manifolds' intersections; new methods are needed. In this study, we first review a method of restricting the intersection search to only certain manifold subsets. We then implement this search using Julia and OpenCL, representing the manifolds as triangular meshes and gaining a 30x speedup using GPUs. We finally show how to use manifold parameterizations to refine the approximate connections found in the mesh search. We demonstrate the tools on the planar elliptic RTBP.

Many unstable periodic orbits of the planar circular restricted 3-body problem (PCRTBP) persist as invariant tori when a periodic forcing is added to the equations of motion. In this study, we compute tori corresponding to exterior Jupiter-Europa and interior Jupiter-Ganymede PCRTBP resonant periodic orbits in a concentric circular restricted 4-body problem (CCR4BP). Motivated by the 2:1 Laplace resonance between Europa and Ganymede's orbits, we then attempt the continuation of a Jupiter-Europa 3:4 resonant orbit from the CCR4BP into the Jupiter-Ganymede PCRTBP. We strongly believe that the resulting dynamical object is a KAM torus lying near but not on the 3:2 Jupiter-Ganymede resonance.