Papers for Friday, May 12 2023

We present a spatially-resolved study of the ionised gas in the central 2 kpc of the Seyfert 2 galaxy NGC 2110 and investigate the role of its moderate luminosity radio jet (kinetic radio power of $P_\mathrm{jet} = 2.3 \times 10^{43}\mathrm{erg\ s^{-1}}$). We use new optical integral-field observations taken with the MEGARA spectrograph at GTC. We fit the emission lines with a maximum of two Gaussian components, except at the AGN position where we used three. Aided by existing stellar kinematics, we use the observed velocity and velocity dispersion of the emission lines to classify the different kinematic components. The disc component is characterised by lines with $\sigma \sim 60-200\ \mathrm{km\ s^{-1}}$. The outflow component has typical values of $\sigma \sim 700\ \mathrm{km\ s^{-1}}$ and is confined to the central 400 pc, which is coincident with linear part of the radio jet detected in NGC 2110. At the AGN position, the [O III]$\lambda$5007 line shows high velocity components reaching at least $1000\ \mathrm{km\ s^{-1}}$. This and the high velocity dispersions indicate the presence of outflowing gas outside the galaxy plane. Spatially-resolved diagnostic diagrams reveal mostly LI(N)ER-like excitation in the outflow and some regions in the disc, which could be due to the presence of shocks. However, there is also Seyfert-like excitation beyond the bending of the radio jet, probably tracing the edge of the ionisation cone that intercepts with the disc of the galaxy. NGC 2110 follows well the observational trends between the outflow properties and the jet radio power found for a few nearby Seyfert galaxies. All these pieces of information suggest that part of observed ionised outflow in NGC 2110 might be driven by the radio jet. However, the radio jet was bent at radial distances of 200 pc (in projection) from the AGN, and beyond there, most of the gas in the galaxy disc is rotating.

We present Keck/DEIMOS spectroscopy of the first complete sample of ultra-diffuse galaxies (UDGs) in the Virgo cluster. We select all UDGs in Virgo that contain at least 10 globular cluster (GC) candidates and are more than $2.5\sigma$ outliers in scaling relations of size, surface brightness, and luminosity (a total of 10 UDGs). We use the radial velocity of their GC satellites to measure the velocity dispersion of each UDG. We find a mixed bag of galaxies: from one UDG that shows no signs of dark matter, to UDGs that follow the luminosity-dispersion relation of early-type galaxies, to the most extreme examples of heavily dark matter dominated galaxies that break well-known scaling relations such as the luminosity-dispersion or the U-shaped total mass-to-light ratio relations. This is indicative of a number of mechanisms at play forming these peculiar galaxies. Some of them may be the most extended version of dwarf galaxies, while others are so extreme that they seem to populate dark matter halos consistent with that of the Milky-Way or even larger. Even though Milky-Way stars and other GC interlopers contaminating our sample of GCs cannot be fully ruled-out, our assessment of this potential problem and simulations indicate that the probability is low and, if present, unlikely to be enough to explain the extreme dispersions measured. Further confirmation from stellar kinematics studies in these UDGs would be desirable. The lack of such extreme objects in any of the state-of-the-art simulations, opens an exciting avenue of new physics shaping these galaxies.

One of the intriguing puzzles concerning Swift J1644+57, the first jetted tidal disruption event (TDE) discovered, is the constant increase in its jet energy, as implied by radio observations. During the first two hundred days the jet energy has increased by an order of magnitude. We suggest that the jet was viewed slightly off-axis. In this case, the apparent energy increase arises due to the slowing down of the jet and the corresponding broadening of its beaming cone. Using equipartition analysis, we infer an increasing jet energy as a larger region of the jet is observed. A simple off-axis model accounts nicely for the multi-wavelength radio observations, resolving this long-standing puzzle. The model allows us to self-consistently evolve the synchrotron signature from an off-axis jet as a function of time. It also allows us to estimate, for the first time, the beaming angle of the jet, $\theta_0 \approx 21^{\circ}$. This implies that the prompt phase of Swift J1644+57 involved super Eddington jet luminosity. We also present a closure relation between the spectral and temporal flux for off-axis jets, which can be used to test whether a given radio transient is off-axis or not.

It is often stated that the observation of high-energy neutrinos from an astrophysical source would constitute a smoking gun for the acceleration of hadronic cosmic rays. Here, we point out that there exists a purely leptonic mechanism to produce TeV-scale neutrinos in astrophysical environments. In particular, very high energy synchrotron photons can scatter with X-rays, exceeding the threshold for muon-antimuon pair production. When these muons decay, they produce neutrinos without any cosmic-ray protons or nuclei being involved. In order for this mechanism to be efficient, the source in question must feature both kG-scale magnetic fields and a high density of keV-scale photons. As an example, we consider the active galaxy NGC 1068, which IceCube has recently detected as a source of TeV-scale neutrinos. We find that the neutrino emission observed from this source could potentially be generated through muon pair production for reasonable choices of physical parameters.

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver.

Stars with masses above 1.6 solar masses generally possess convective cores and radiative envelopes, which allows the propagation of outward-travelling internal gravity waves. We have studied the generation and propagation of IGWs in such stars using two-dimensional, fully nonlinear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics. Compared to previous similar works, this study utilises radius-dependent thermal diffusivity profiles for 5 different stellar masses at the middle of main sequence: 3 - 13 solar masses. From the simulations, we find that the surface perturbations are larger for higher masses, but no noticeable trends are observed for the frequency slopes with different stellar masses. The slopes are also similar to the results from previous works. We compare our simulation results with stellar photometric data from a recent survey and find that for frequency intervals above 8 microHz, there is a good agreement between the temperature frequency slopes from the simulations and the surface brightness variations of these observed stars, indicating that the brightness variations are caused by core-generated IGWs.

Galactic bars, made up of elongated and aligned stellar orbits, can lose angular momentum via resonant torques with dark matter particles in the halo and slow down. Here we show that if a stellar bar is decelerated to zero rotation speed, it can flip the sign of its angular momentum and reverse rotation direction. We demonstrate this in a collisionless N-body simulation of a galaxy in a live counter-rotating halo. Reversal begins at small radii and propagates outward. The flip generates a kinematically-decoupled core both in the visible galaxy and in the dark matter halo, and counter-rotation generates a large-scale warp of the outer disk with respect to the bar.

Intracluster light is thought to originate from stars that were ripped away from their parent galaxies by gravitational tides and galaxy interactions during the build up of the cluster. The stars from such interactions will accumulate over time, so semi-analytic models suggest that the abundance of intracluster stars is negligible in young proto-clusters at z~2 and grows to around a quarter of the stellar mass in the oldest, most mature clusters. In contrast to these theoretical expectations, we report on the detection of intracluster light within two proto-clusters at z=2 using deep HST images. We use the colour of the intracluster light to estimate its mass-to-light ratio in annuli around the brightest cluster galaxies (BCG), up to a radius of 100 kpc. We find that $54\pm5$\% and $71\pm3$\% of the stellar mass in these regions is located more than 10 kpc away from the BCGs in the two proto-clusters. This low concentration is similar to BCGs in lower redshift clusters, and distinct from other massive proto-cluster galaxies. This suggests that intracluster stars are already present within the core 100 kpc of proto-clusters. We compare these observations to the Hydrangea hydrodynamical galaxy cluster simulations and find that intracluster stars are predicted to be a generic feature of group-sized halos at z=2. These intracluster stars will gradually move further away from the BCG as the proto-cluster assembles into a cluster.

We analyze the dust properties of 57 870 $\mu$m selected dusty star-forming galaxies in the GOODS-S using new deep ALMA 1.2 mm, 2 mm, and 3 mm continuum imaging together with other far-infrared through millimeter data. We fit the spectral energy distributions (SEDs) with optically thin modified blackbodies to constrain the emissivity indices and effective dust temperatures, finding a median emissivity index of $\beta = 1.78^{+0.43}_{-0.25}$ and a median temperature of $T_d = 33.6^{+12.1}_{-5.4}$ K. We observe a negative correlation between $\beta$ and $T_d$. By testing several SED models, we determine that the derived emissivity indices can be influenced by opacity assumptions. Our temperature measurements are consistent with no evolution in dust temperature with redshift.

Nowadays, efforts are being devoted to the study of alternative cosmological scenarios, in which, modifications of General Relativity (GR) theory have been proposed to explain the late cosmic acceleration, without assuming the existence of the dark energy (DE) component. We investigate the $R^2$-corrected Appleby-Battye model, or $R^2$-AB model, which consists of an $f(R)$ model with only one extra free parameter $b$, besides the cosmological parameters of the flat-$\Lambda$CDM model: $H_0$ and $\Omega_{m,0}$. Regarding this model, it was already shown that a positive value for $b$ is required for the model to be consistent with Solar System tests, moreover, the condition for the existence of a de~Sitter state requires $b \ge 1.6$. To impose observational constraints on the $R^2$-AB model, we consider in our analyses two data sets: cosmic chronometer $H(z)$ data for the background level, and $[f\sigma_8](z)$ data, for the perturbative level. The first one provides $b = 1.6^{+3.1}_{-0.0}$ and the cosmological parameters $\{H_0 ,\Omega_{m,0}\}$ in agreement to Planck values, while the second one, indicates $b = 1.76^{+2.91}_{-0.15}$ and the parameters $\{\Omega_{m,0},\sigma_{8,0} \}$ also in agreement to Planck values; in the last case the data was marginalized over the parameter $H_0$. Additionally, we perform illustrative analyses that compare this $f(R)$ model with the flat-$\Lambda$CDM model, considering several values of the parameter $b$, for diverse cosmological functions like the Hubble function $H(z)$, the equation of state $w_{eff}(z)$, the parametrized growth rate of cosmic structures $[f \sigma_8](z)$, and $\sigma_8(z)$. The overall conclusion is that the $R^2$-AB model is a promising $f(R)$ model that deserves to continue being tested with diverse cosmological data.

LOFAR (LOw Frequency ARray) has previously detected bursts from the periodically active, repeating fast radio burst (FRB) source FRB 20180916B down to unprecedentedly low radio frequencies of 110 MHz. Here we present 11 new bursts in 223 more hours of continued monitoring of FRB 20180916B in the 110-188 MHz LOFAR band. We place new constraints on the LOFAR activity width $w = 4.3^{+0.7}_{-0.2}$ day, and phase centre $\phi_{\mathrm{c}}^{\mathrm{LOFAR}} = 0.67^{+0.03}_{-0.02}$ in its 16.33-day activity cycle, strengthening the evidence for its frequency-dependent activity cycle. Propagation effects like Faraday rotation and scattering are especially pronounced at low frequencies and constrain properties of FRB 20180916B's local environment. We track variations in scattering and time-frequency drift rates, and find no evidence for trends in time or activity phase. Faraday rotation measure (RM) variations seen between June 2021 and August 2022 show a fractional change $>$50% with hints of flattening of the gradient of the previously reported secular trend seen at 600 MHz. The frequency-dependent window of activity at LOFAR appears stable despite the significant changes in RM, leading us to deduce that these two effects have different causes. Depolarization of and within individual bursts towards lower radio frequencies is quantified using LOFAR's large fractional bandwidth, with some bursts showing no detectable polarization. However, the degree of depolarization seems uncorrelated to the scattering timescales, allowing us to evaluate different depolarization models. We discuss these results in the context of models that invoke rotation, precession, or binary orbital motion to explain the periodic activity of FRB 20180916B.

Abundances of s-capture process elements in stars with exquisite asteroseismic, spectroscopic, and astrometric constraints offer a novel opportunity to study stellar evolution, nucleosynthesis, and Galactic chemical evolution. We aim to investigate one of the least studied s-process elements in the literature, Ce, using stars with asteroseismic constraints from the Kepler, K2 and TESS missions. We combine the global asteroseismic parameters derived from precise light curves obtained by the Kepler, K2 and TESS missions with chemical abundances from the APOGEE DR17 survey and astrometric data from the Gaia mission. Finally, we compute stellar ages using the code PARAM. We investigate the different trends of [Ce/Fe] as a function of [Fe/H], [alpha/Fe] and age considering the dependence on the radial position, specially in the case of K2 targets which cover a large Galactocentric range. We, finally, explore the [Ce/alpha] ratios as a function of age in different Galactocentric intervals. The studied trends display a strong dependence of the Ce abundances on [Fe/H] and star formation history. Indeed, the [Ce/Fe] ratio shows a non-monotonic dependence on [Fe/H] with a peak around -0.2 dex. Moreover, younger stars have higher [Ce/Fe] and [Ce/alpha] ratios than older stars, confirming the latest contribution of low- and intermediate-mass asymptotic giant branch stars to the Galactic chemical enrichment. In addition, the trends of [Ce/Fe] and [Ce/alpha] with age become steeper moving towards the outer regions of the Galactic disc, demonstrating a more intense star formation in the inner regions than in the outer regions. Ce is thus a potentially interesting element to help constraining stellar yields and the inside-out formation of the Milky Way disc. However, the large scatter in all the relations studied here, suggests that spectroscopic uncertainties for this element are still too large.

We analyze the physical properties of 8 X-ray selected active galactic nuclei (AGN) and one candidate protoquasar system (ADF22A1) in the $z = 3.09$ SSA22 protocluster by fitting their X-ray-to-IR spectral energy distributions (SEDs) using our SED fitting code, Lightning. We recover star formation histories (SFH) for 7 of these systems which are well-fit by composite stellar population plus AGN models. We find indications that 4/9 of the SSA22 AGN systems we study have host galaxies below the main sequence, with $\rm SFR/SFR_{MS} \leq -0.4$. The remaining SSA22 systems, including ADF22A1, are consistent with obscured supermassive black hole (SMBH) growth in star forming galaxies. We estimate the SMBH accretion rates and masses, and compare the properties and SFH of the 9 protocluster AGN systems with X-ray detected AGN candidates in the Chandra Deep Fields (CDF), finding that the distributions of SMBH growth rates, star formation rates, SMBH masses, and stellar masses for the protocluster AGN are consistent with field AGN. We constrain the ratio between the sample-averaged SSA22 SMBH mass and CDF SMBH mass to $<1.41$. While the AGN are located near the density peaks of the protocluster, we find no statistically significant trends between the AGN or host galaxy properties and their location in the protocluster. We interpret the similarity of the protocluster and field AGN populations together with existing results as suggesting that the protocluster and field AGN co-evolve with their hosts in the same ways, while AGN-triggering events are more likely in the protocluster.

Wide-field and deep DECam multi-band photometry, combined with HST data for the core of the Galactic globular cluster NGC 2808, allowed us to study the distribution of various stellar sub-populations and stars in different evolutionary phases out to the cluster tidal radius. We used the C_ugi = (u-g)-(g-i) index to identify three chemically distinct sub-populations along the red giant branch and compared their spatial distributions. The most light-element enriched sub-population (P3) is more centrally concentrated; however, it shows a more extended distribution in the external regions of the cluster compared to the primordial (P1) and intermediate (P2) composition populations. Furthermore, the P3 sub-population centroid is off-center relative to those of the P1 and P2 groups. We also analyzed the spatial distribution of horizontal branch stars and found that the relative fraction of red horizontal branch stars increases for radial distances larger than ~ 1.5' while that of the blue and hotter stars decreases. These new observations, combined with literature spectroscopic measurements, suggest that the red horizontal branch stars are the progeny of all the stellar sub-populations in NGC 2808, i.e. primordial and light-element enhanced, while the blue stars are possibly the result of a combination of the "hot-flasher" and the "helium-enhanced" scenarios. A similar distribution of different red giant branch sub-populations and horizontal branch stars was also found for the most massive Galactic globular cluster, omega Cen, based on combined DECam and HST data, which suggests the two may share a similar origin.

Primordial magnetic fields (PMFs) offer a compelling explanation for the origin of observed magnetic fields, especially on extragalactic scales. Such PMFs give rise to excess of power in small scale matter perturbations that could strongly influence structure formation. We study the impact of the magnetically enhanced matter power spectrum on the signal that will be observed by line-intensity mapping (LIM) surveys targeting carbon monoxide (CO) emission from star-forming galaxies at high redshifts. Specifically, the voxel intensity distribution of intensity maps provides access to small-scale information, which makes it highly sensitive to signatures of PMFs on matter overdensities. We present forecasts for future LIM CO surveys, finding that they can constrain PMF amplitudes as small as $\sigma_{B,0}\sim0.04-1\,{\rm nG}$, depending on the magnetic spectral index and the targeted redshifts.

Due to the non-linear ionizing and heating processes, the 21-cm signals from epoch of reionization (EoR) are expected to have strong non-Gaussian fluctuations. In this paper, we use the semi-numerical simulations to study the non-Gaussian statistics i.e. skew spectrum and smoothed skewness of the 21-cm signals from EoR. We find the 21-cm skew spectrum and smoothed skewness have similar evolution features with the 21-cm bispectrum. All of them are sensitive to the EoR models, while not too much to the cosmic volume applied. With the SKA1-low telescope as reference, we find both the skew spectrum and smoothed skewness have much higher S/N ratios than the 21-cm bispectrum.

Sunspot Cycle 25 is now over 3 years past the cycle minimum of December 2019. At this point in the cycle, curve-fitting to the activity becomes reliable and now consistently indicates a maximum sunspot number of 135 +/- 10 - slightly larger than Cycle 24's maximum of 116.4, but well below the average of 179. A geomagnetic precursor, the minimum in the aa-index, and the Sun's magnetic precursors, the Sun's polar field strength and its axial dipole moment at the time of minimum, are often used to predict the amplitude of the cycle at (or before) the onset of the cycle. We examine Cycle 25 predictions produced by these precursors. The geomagnetic precursor indicated a Cycle 25 slightly stronger that Cycle 24, with a maximum of 132 +/- 8. The Sun's magnetic precursors indicated that Cycle 25 would be more similar to Cycle 24, with a maximum sunspot number of 120 +/- 10 or 114 +/- 15. Combining the curve-fitting results with the precursor predictions, we conclude that Cycle 25 will have a maximum smoothed sunspot number of 134 +/- 8 with maximum occurring late in the fall of 2024. Models for predicting the Sun's magnetic field ahead of minimum, were generally successful at predicting the polar precursors years in advance. The fact that Sun's magnetic precursors at cycle minimum were successfully predicted years before minimum and that the precursors are consistent with the size of Cycle 25 suggests that we can now reliably predict the solar cycle.

The variability of submillimeter emission provides a useful tool to probe the accretion physics in low-luminosity active galactic nuclei. We accumulate four years of observations using Submillimeter Array for Centaurus A, NGC 4374, NGC 4278, and NGC 5077 and one year of observations for NGC 4552 and NGC 4579. All sources are variable. We measure the characteristic timescale at which the variability is saturated by modeling these sources' light curve as a damped random walk. We detect a timescale for all the sources except NGC 4552. The detected timescales are comparable to the orbital timescale at the event horizon scale for most sources. Combined with previous studies, we show a correlation between the the timescale and the black hole mass over three orders of magnitude. This discovery suggests the sub-mm emission is optically thin with the emission originating from the event horizon. The mass scaling relationship further suggests that a group of radio sources with a broadband spectrum that peaks at submillimeter wavelengths have similar inner accretion physics. Sources that follow this relationship may be good targets for high-resolution imaging with the Event Horizon Telescope.

The formation path to ultra-compact X-ray binaries (UCXBs) with black hole (BH) accretors is still unclear. In the classical formation scenario, it is difficult to eject the massive envelope of the progenitor star of the BH via common envelope process. Given that some neutron stars (NSs) in binary systems evidently have birth masses close to $\sim 2.0\;M_\odot$, we explore here the possibility that BH-UCXBs may form via accretion-induced collapse (AIC) of accreting NSs, assuming that these previously evolved in LMXBs to masses all the way up to the maximum limit of a NS. We demonstrate this formation path by modelling a few cases of NS-UCXBs with initial NS masses close to the maximum mass of a NS that evolve into BH-UCXBs after the NS accretes material from its He~WD companion. We follow the evolution of the post-AIC BH-UCXB and, based on simple arguments, we anticipate that there is about one BH-UCXB with an AIC origin and a He~WD donor within the current sample of known UCXBs and that 2--5 such BH-UCXBs may be detected in gravitational waves by LISA. In addition, we find that the X-ray luminosity of NS-UCXBs near their orbital period minimum exceeds $\sim 10^{39}\;{\rm erg\;s^{-1}}$ and thus such systems may appear as ultraluminous X-ray sources.

Low mass stars ($\leq 1$ M$_{\odot}$) are some of the best candidates for hosting planets with detectable life because of these stars' long lifetimes and relative planet to star mass and radius ratios. An important aspect of these stars to consider is the amount of ultraviolet (UV) and X-ray radiation incident on planets in the habitable zones due to the ability of UV and X-ray radiation to alter the chemistry and evolution of planetary atmospheres. In this work, we build on the results of the HAZMAT I (Shkolnik & Barman 2014) and HAZMAT III (Schneider & Shkolnik 2018) M star studies to determine the intrinsic UV and X-ray flux evolution with age for M stars using Gaia parallactic distances. We then compare these results to the intrinsic fluxes of K stars adapted from HAZMAT V (Richey-Yowell et al. 2019). We find that although the intrinsic M star UV flux is 10 to 100 times lower than that of K stars, the UV fluxes in their respective habitable zone are similar. However, the habitable zone X-ray flux evolutions are slightly more distinguishable with a factor of 3 -- 15 times larger X-ray flux for late-M stars than for K stars. These results suggest that there may not be a K dwarf advantage compared to M stars in the UV, but one may still exist in the X-ray.

We inspect the microlensing data of the KMTNet survey collected during the 2018--2020 seasons in order to find lensing events produced by binaries with brown-dwarf companions. In order to pick out binary-lens events with candidate BD lens companions, we conduct systematic analyses of all anomalous lensing events observed during the seasons. By applying the selection criterion with mass ratio between the lens components of $0.03\lesssim q\lesssim 0.1$, we identify four binary-lens events with candidate BD companions, including KMT-2018-BLG-0321, KMT-2018-BLG-0885, KMT-2019-BLG-0297, and KMT-2019-BLG-0335. For the individual events, we present the interpretations of the lens systems and measure the observables that can constrain the physical lens parameters. The masses of the lens companions estimated from the Bayesian analyses based on the measured observables indicate that the probabilities for the lens companions to be in the brown-dwarf mass regime are high: 59\%, 68\%, 66\%, and 66\% for the four events respectively.

We present new results from INTEGRAL and Swift observations of the hitherto poorly studied and unidentified X-ray source XTE J1906+090. A bright hard X-ray outburst (luminosity of $\sim$10$^{36}$ erg s$^{-1}$ above 20 keV) has been discovered with INTEGRAL observations in 2010, this being the fourth outburst ever detected from the source. Such events are sporadic, the source duty cycle is in the range (0.8--1.6)% as inferred from extensive INTEGRAL and Swift monitoring in a similar hard X-ray band. Using five archival unpublished Swift/XRT observations, we found that XTE J1906+090 has been consistently detected at a persistent low X-ray luminosity value of $\sim$10$^{34}$ erg s$^{-1}$, with limited variability (a factor as high as 4). Based on our findings, we propose that XTE J1906+090 belongs to the small and rare group of persistent low luminosity Be X-ray Binaries.

Recently, the Fermi-LAT gamma-ray satellite has detected six Giant Molecular Clouds (GMCs) located in the Gould Belt and the Aquila Rift regions. In half of these objects (Taurus, Orion A, Orion B), the observed gamma-ray spectrum can be explained using the Galactic diffused Cosmic Ray (CR) interactions with the gas environments. In the remaining three GMCs (Rho Oph, Aquila Rift, Cepheus), the origin of the gamma-ray spectrum is still not well established. We use the GEometry ANd Tracking (GEANT4) simulation framework in order to simulate gamma-ray emission due to CR/GMC interaction in these three objects, taking into account the gas density distribution inside the GMCs. We find that propagation of diffused Galactic CRs inside these GMCs can explain the Fermi-LAT detected gamma-ray spectra. Further, our estimated TeV-PeV fluxes are consistent with the HAWC upper limits, available for the Aquila Rift GMC. As last step, we compute the total neutrino flux estimated for these GMCs and compare it with the IceCube detection sensitivity.

We present a design study for a new production technology for ultra-high vacuum pipes. The pipes are produced in a fully automatised process in sections of hundreds of meters directly in the later location of usage. We estimate the effort for such a production and show that it might be substantially lower than the effort for an off-site production of transportable sections.

We analyse VLBI and optical images of AGNs and their host galaxies and look for statistical correlations between the shape and orientation of the galaxy and the direction of the jet. We utilise the Astrogeo catalogue, which has over 9000 VLBI sources, many of those with a clear core-jet like structure that allows for the jet position angle to be reliably determined. We then use the VLBI source positions to search for optical counterparts within various optical surveys. In order to parameterise the orientation and shape of the host galaxy, we fitted a Gaussian elliptical model to the optical image, taking the PSF into account. We check our own shape parameters from this fit against the ones provided by the optical surveys. As of yet, no clear correlation between the galaxy morphology and the jet direction is seen.

Tensions between the diffuse gamma-ray sky observed by the Fermi Large Area Telescope (LAT) and the diffuse high-energy neutrino sky detected by the IceCube South Pole Neutrino Observatory question our knowledge about high-energy neutrino sources in the gamma-ray regime. While blazars are among the most energetic persistent particle accelerators in the Universe, studies suggest that they could account for up to for 10-30% of the neutrino flux measured by IceCube. Our recent results highlighted that the associated IceCube neutrinos arrived in a local gamma-ray minimum (dip) of three strong neutrino point-source candidates. We increase the sample of neutrino-source candidates to study their gamma-ray light curves. We generate the one-year Fermi-LAT light curve for 8 neutrino source candidate blazars (RBS 0958, GB6 J1040+0617, PKS 1313-333, TXS 0506+056, PKS 1454-354, NVSS J042025-374443, PKS 0426-380 and PKS 1502+106), centered on the detection time of the associated IceCube neutrinos. We apply the Bayesian block algorithm on the light curves to characterize their variability. Our results indicate that GB6 J1040+0617 was in the phase of high gamma-ray activity, while none of the other 7 neutrino source candidates were statistically bright during the detection of the corresponding neutrinos and that indeed even most of the times neutrinos arrived in a faint gamma-ray phase of the light curves. This suggests that the 8 source-candidate blazars (associated with 7 neutrino events) in our reduced sample are either not the sources of the corresponding IceCube neutrinos, or that an in-source effect (e.g. suppression of gamma rays due to high gamma-gamma opacity) complicates the multimessenger scenario of neutrino emission for these blazars.

In N-body simulations, nearly radial mergers can form shell-like overdensities in the sky position and phase space ($r-v_r$) due to the combination of dynamical friction and tidal stripping. The merger event of Gaia-Sausage-Enceladus has provided a unique opportunity to study the shells in the phase space. To search for them, we integrate the orbits of 5949 GSE-related halo K giants from the LAMOST survey and record their positions at all time intervals in $r-v_r$ diagram. After the subtraction of a smoothed background, we find six significant and complete thin chevron-like overdensities. The apocenters $r_\mathrm{apo}$ of stars in the six chevrons are around 6.75, 12.75, 18.75, 25.25, 27.25, and 30.25 kpc. These chevrons reveal the multiple pile-ups of GSE stars at different apocenters. The application of a different Milky Way mass $M_\mathrm{vir}$ will change the opening angles of these chevrons, while leave their apocenters almost unchanged. By comparing with a recent study of the phase space overdensities of local halo stars from Gaia RVS survey, our results are more inclined to a medium $M_\mathrm{vir}$ of $10^{12}\,M_\odot$. The application of a non-axisymmetric Galactic potential with a steadily rotating bar has a blurring effect on the appearance of these chevron-like overdensities, especially for the chevrons with $r_\mathrm{apo} > 20$ kpc.

Anomalous Microwave Emission (AME) is an important emission component between 10 and 60 GHz that is not yet fully understood. It seems to be ubiquituous in our Galaxy and is observed at a broad range of angular scales. Here we use the new QUIJOTE-MFI wide survey data at 11, 13, 17 and 19 GHz to constrain the AME in the Galactic plane ($|b|<10^\circ$) on degree scales. We built the spectral energy distribution between 0.408 and 3000 GHz for each of the 5309 0.9$^\circ$, pixels in the Galactic plane, and fitted a parametric model by considering five emission components: synchrotron, free-free, AME, thermal dust and CMB anisotropies. We show that not including QUIJOTE-MFI data points leads to the underestimation (up to 50 %) of the AME signal in favour of free-free emission. The parameters describing these components are then intercompared, looking for relations that help to understand AME physical processes. We find median values for the AME width, $W_{\rm AME}$, and for its peak frequency, $\nu_{\rm AME}$, respectively of $0.560^{+0.059}_{-0.050}$ and $20.7^{+2.0}_{-1.9}$ GHz, slightly in tension with current theoretical models. We find spatial variations throughout the Galactic plane for $\nu_{\rm AME}$, but only with reduced statistical significance. We report correlations of AME parameters with certain ISM properties, such as that between the AME emissivity (which shows variations with the Galactic longitude) and the interstellar radiation field, and that between the AME peak frequency and dust temperature. Finally, we discuss the implications of our results on the possible molecules responsible for AME.

We present a novel test of the cosmological principle: the idea that, on sufficiently large scales, the universe should appear homogeneous and isotropic to observers comoving with the Hubble flow. This is a fundamental assumption in modern cosmology, underpinning the use of the Friedmann-Lema\^itre-Robertson-Walker metric as part of the concordance $\Lambda$CDM paradigm. However, the observed dipole imprinted on the Cosmic Microwave Background (CMB) is interpreted as our departure from the Hubble flow, and such a proper motion will induce a directionally-dependent time dilation over the sky. We illustrate the feasibility of detection of this 'time dilation dipole' and sketch the practical steps involved in its extraction from a catalogue of sources with intrinsic time-scales. In essence, whilst the scale of this dilation is small, being of order of 0.1%, it will in principle be detectable in large scale surveys of variable cosmological sources, such as quasars and supernovae. The degree of alignment of the time dilation dipole with the kinematic dipole derived from the CMB will provide a new assessment of the cosmological principle, and address the tension in dipole measures from other observations.

Neptune's tropospheric winds are among the most intense in the Solar System, but the dynamical mechanisms that produce them remain uncertain. Measuring wind speeds at different pressure levels may help understand the atmospheric dynamics of the planet. The goal of this work is to directly measure winds in Neptune's stratosphere with ALMA Doppler spectroscopy. We derived the Doppler lineshift maps of Neptune at the CO(3-2) and HCN(4-3) lines at 345.8 GHz ($\lambda$~0.87 mm) and 354.5 GHz (0.85 mm), respectively. For that, we used spectra obtained with ALMA in 2016 and recorded with a spatial resolution of ~0.37" on Neptune's 2.24" disk. After subtracting the planet solid rotation, we inferred the contribution of zonal winds to the measured Doppler lineshifts at the CO and HCN lines. We developed an MCMC-based retrieval methodology to constrain the latitudinal distribution of wind speeds. We find that CO(3-2) and HCN(4-3) lines probe the stratosphere of Neptune at pressures of $2^{+12}_{-1.8}$ mbar and $0.4^{+0.5}_{-0.3}$ mbar, respectively. The zonal winds at these altitudes are less intense than the tropospheric winds based on cloud tracking from Voyager observations. We find equatorial retrograde (westward) winds of $-180^{+70}_{-60}$ m/s for CO, and $-190^{+90}_{-70}$ m/s for HCN. Wind intensity decreases towards mid-latitudes, and wind speeds at 40$^\circ$S are $-90^{+50}_{-60}$ m/s for CO, and $-40^{+60}_{-80}$ m/s for HCN. Wind speeds become 0 m/s at about 50$^\circ$S, and we find that the circulation reverses to a prograde jet southwards of 60$^\circ$S. Overall, our direct stratospheric wind measurements match previous estimates from stellar occultation profiles and expectations based on thermal wind equilibrium. These are the first direct Doppler wind measurements performed on the Icy Giants, opening a new method to study and monitor their stratospheric dynamics.

Accurate quantification of the signal-to-noise ratio (SNR) of a given observational phenomenon is central to associated calculations of sensitivity, yield, completeness and occurrence rate. Within the field of exoplanets, the SNR of a transit has been widely assumed to be the formula that one would obtain by assuming a boxcar light curve, yielding an SNR of the form $(\delta/\sigma_0) \sqrt{D}$. In this work, a general framework is outlined for calculating the SNR of any analytic function and it is applied to the specific case of a trapezoidal transit as a demonstration. By refining the approximation from boxcar to trapezoid, an improved SNR equation is obtained that takes the form $(\delta/\sigma_0) \sqrt{(T_{14}+2T_{23})/3}$. A solution is also derived for the case of a trapezoid convolved with a top-hat, corresponding to observations with finite integration time, where it is proved that SNR is a monotonically decreasing function of integration time. As a rule of thumb, integration times exceeding $T_{14}/3$ lead to a 10% loss in SNR. This work establishes that the boxcar transit is approximate and it is argued that efforts to calculate accurate completeness maps or occurrence rate statistics should either use the refined expression, or even better numerically solve for the SNR of a more physically complete transit model.

Based on the DECaLS shear catalog, we study the scaling relations between halo mass($M_{\rm h}$) and various proxies for SDSS central galaxies, including stellar mass($M_*$), stellar velocity dispersion($\sigma_*$), abundance matching halo mass($M_{\rm AM}$) and satellite velocity dispersion($\sigma_{\rm s}$), and their dependencies on galaxy and group properties. In general, they are all good proxies of $M_{\rm h}$, and their scaling relations are consistent with previous studies. We find that the $M_{\rm h}$-$M_*$ and $M_{\rm h}$-$\sigma_*$ relations depend strongly on group richness($N_{\rm sat}$), while the $M_{\rm h}$-$M_{\rm AM}$ and $M_{\rm h}$-$\sigma_{\rm s}$ relations are independent of it. Moreover, the dependence on star formation rate(SFR) is rather weak in the $M_{\rm h}$-$\sigma_*$ and $M_{\rm h}$-$\sigma_{\rm s}$ relations, but very prominent in the other two. $\sigma_{\rm s}$ is thus the best proxy among them, and its scaling relation is in good agreement with hydro-dynamical simulations. However, estimating $\sigma_{\rm s}$ accurately for individual groups/clusters is challenging because of interlopers and the requirement for sufficient satellites. We construct new proxies by combining $M_*$, $\sigma_*$, and $M_{\rm AM}$, and find the proxy with 30\% contribution from $M_{\rm AM}$ and 70\% from $\sigma_*$ can minimize the dependence on $N_{\rm sat}$ and SFR. We obtain the $M_{\rm h}$-supermassive black hole(SMBH) mass relation via the SMBH scaling relation and find indications for rapid and linear growth phases for SMBH. We also find that correlations among $M_{\rm h}$, $M_*$ and $\sigma_*$ change with $M_*$, indicating that different processes drive the growth of galaxies and SMBH at different stages.

To decipher complex patterns of gravity-mode period spacings observed for intermediate-mass main-sequence stars is an important step toward the better understanding of the structure and dynamics in the deep radiative region of the stars. In this study, we apply JWKB approximation to derive a semi-analytical expression of the g-mode period spacing pattern, for which the gradient in the Brunt-V\"ais\"al\"a frequency is taken into account. The formulation includes a term $P^{-1} B_{\star}$, where $P$ and $B_{\star}$ represent the g-mode period and degree of the structural variation, the latter of which especially is related to the steepness of the gradient of the Brunt-V\"ais\"al\"a frequency. Tests with 1-dimensional stellar models show that the semi-analytical expression derived in this study is useful for inferring the degree of the structural variation $B_{\star}$ with accuracy of $\sim 10\,\%$ in the case of relatively massive intermediate-mass models with the mass $M$ larger than $3 \,M_{\odot}$. The newly formulated expression will possibly allow us to put further constraints on, e.g., mixing processes inside intermediate-mass main-sequence g-mode pulsators such as $\beta$ Cep, SPB, and $\gamma$ Dor stars that have been principal targets in asteroseismology.

Massive close binary stars with extremely small separations have been observed, and they are possible progenitors of gravitational-wave sources. The evolution of massive binaries in the protostellar accretion stage is key to understanding their formation process. We, therefore, investigate how close the protostars, consisting of a high-density core and a vast low-density envelope, can approach each other but not coalesce. To investigate the coalescence conditions, we conduct smoothed particle hydrodynamics simulations following the evolution of equal-mass binaries with different initial separations. Since Population (Pop) I and III protostars have similar interior structures, we adopt a specific Pop~III model with the mass and radius of $7.75\;M_{\odot}$ and $61.1\;R_{\odot}$ obtained by the stellar evolution calculations. Our results show that the binary separation decreases due to the transport of the orbital angular momentum to spin angular momentum. If the initial separation is less than about 80 per~cent of the sum of the protostellar radius, the binary coalesces in a time shorter than the tidal lock timescale. The mass loss up to the merging is $\lesssim 3$ per~cent. After coalescence, the star rotates rapidly, and its interior structure is independent of the initial separation. We conclude that there must be some orbital shrinking mechanism after the protostars contract to enter the zero-age main-sequence stage.

We aim to characterise the small-scale magnetic fields for a sample of 16 Sun-like stars and investigate the capabilities of the newly upgraded near-infrared (NIR) instrument CRIRES$^+$ at the VLT in the context of small-scale magnetic field studies. Our targets also had their magnetic fields studied in the optical, which allows us to compare magnetic field properties at different spatial scales on the stellar surface and to contrast small-scale magnetic field measurements at different wavelengths. We analyse the Zeeman broadening signature for six magnetically sensitive and insensitive \ion{Fe}{I} lines in the H-band to measure small-scale magnetic fields on the stellar surface. We use polarised radiative transfer modelling and NLTE departure coefficients in combination with MCMC to determine magnetic field characteristics together with non-magnetic stellar parameters. We use two different approaches to describe small-scale magnetic fields. The first is a two-component model with a single magnetic region and a free magnetic field strength. The second model contains multiple magnetic components with fixed magnetic field strengths. We find average magnetic field strengths ranging from $\sim 0.4$ kG down to $<0.1$ kG. The results align closely with other results from high resolution NIR spectrographs such as SPIRou. We find that the small-scale fields correlate with the large-scale fields and that the small-scale fields are at least 10 times stronger than the large-scale fields inferred with Zeeman Doppler imaging. The two- and multi-component models produce systematically different results as the strong fields from the multi-component model increase the obtained mean magnetic field strength. When comparing our results with the optical measurements of small-scale fields we find a systematic offset of 2--3 times stronger fields in the optical.

Solar coronal mass ejections (CMEs) can catch up and interact with preceding CMEs and solar wind structures to undergo rotation and deflection during their propagation. We aim to show how interactions undergone by a CME in the corona and heliosphere can play a significant role in altering its geoeffectiveness predicted at the time of its eruption. We consider a case study of two successive CMEs launched from the active region NOAA 12158 in early September 2014. The second CME was predicted to be extensively geoeffective based on the remote-sensing observations of the source region. However, in situ measurements at 1~au recorded only a short-lasting weak negative Bz component followed by a prolonged positive Bz component. The EUropean Heliosphere FORecasting Information Asset (EUHFORIA) is used to perform a self-consistent 3D MHD simulation of the two CMEs in the heliosphere. The initial conditions of the CMEs are determined by combining observational insights near the Sun, fine-tuned to match the in situ observations near 1~au, and additional numerical experiments of each individual CME. By introducing CME1 before CME2 in the EUHFORIA simulation, we modelled the negative Bz component in the sheath region ahead of CME2 whose formation can be attributed to the interaction between CME1 and CME2. To reproduce the positive Bz component in the magnetic ejecta of CME2, we had to initialise CME2 with an orientation determined at 0.1~au and consistent with the orientation interpreted at 1~au, instead of the orientation observed during its eruption. EUHFORIA simulations suggest the possibility of a significant rotation of CME2 in the low corona in order to explain the in situ observations at 1~au. Coherent magnetic field rotations, potentially geoeffective, can be formed in the sheath region as a result of CME-CME interactions in the heliosphere even if the individual CMEs are not geoeffective.

Microstreams are fluctuations in the solar wind speed and density associated with polarity-reversing folds in the magnetic field (also denoted switchbacks). Despite their long heritage, the origin of these microstreams/switchbacks remains poorly understood. For the first time, we investigated periodicities in microstreams during Parker Solar Probe (PSP) Encounter 10 to understand their origin. Our analysis was focused on the inbound corotation interval on 2021 November 19-21, while the spacecraft dove toward a small area within a coronal hole (CH). Solar Dynamics Observatory remote-sensing observations provide rich context for understanding the PSP in-situ data. Extreme ultraviolet images from the Atmospheric Imaging Assembly reveal numerous recurrent jets occurring within the region that was magnetically connected to PSP during intervals that contained microstreams. The periods derived from the fluctuating radial velocities in the microstreams (approximately 3, 5, 10, and 20 minutes) are consistent with the periods measured in the emission intensity of the jetlets at the base of the CH plumes, as well as in larger coronal jets and in the plume fine structures. Helioseismic and Magnetic Imager magnetograms reveal the presence of myriad embedded bipoles, which are known sources of reconnection-driven jets on all scales. Simultaneous enhancements in the PSP proton flux and ionic ($^3$He, $^4$He, Fe, O) composition during the microstreams further support the connection with jetlets and jets. In keeping with prior observational and numerical studies of impulsive coronal activity, we conclude that quasiperiodic jets generated by interchange/breakout reconnection at CH bright points and plume bases are the most likely sources of the microstreams/switchbacks observed in the solar wind.

Atmospheric mass loss plays a major role in the evolution of exoplanets. This process is driven by the stellar high-energy irradiation, especially in the first hundreds of millions of years after dissipation of the proto-planetary disk. A major source of uncertainty in modeling atmospheric photo-evaporation and photo-chemistry is due to the lack of direct measurements of the stellar flux at EUV wavelengths. Several empirical relationships have been proposed in the past to link EUV fluxes to emission levels in X-rays, but stellar samples employed for this aim are heterogeneous, and available scaling laws provide significantly different predictions, especially for very active stars. We present new UV and X-ray observations of V1298 Tau with HST/COS and XMM-Newton, aimed to determine more accurately the XUV emission of this solar-mass pre-Main Sequence star, which hosts four exoplanets. Spectroscopic data were employed to derive the plasma emission measure distribution vs.\ temperature, from the chromosphere to the corona, and the possible variability of this irradiation on short and year-long time scales, due to magnetic activity. As a side result, we have also measured the chemical abundances of several elements in the outer atmosphere of V1298 Tau. We employ our results as a new benchmark point for the calibration of the X-ray to EUV scaling laws, and hence to predict the time evolution of the irradiation in the EUV band, and its effect on the evaporation of exo-atmospheres.

The diffuse gamma-ray emission between 10 and 1000 TeV from the Galactic plane was recently measured precisely by the Large High Altitude Air Shower Observatory (LHAASO), which is very useful in constraining the propagation and interaction of cosmic rays in the Milky Way. On the other hand, new measurements of CR spectra reach a very high precision up to 100 TeV energies, revealing multiple spectral structures of various species. In this work, we confront the model prediction of the diffuse gamma-ray emission, based on up-to-date measurements of the local cosmic ray spectra and simplified propagation setup, with the measurements of diffuse gamma-rays. To better constrain the low-energy part of the model, we analyze Fermi-LAT data to extract the diffuse emission between 1 and 500 GeV from the same sky regions of LHAASO. Compared with the prediction, we find that clear excesses between several GeV and ~60 TeV of the diffuse emission exist. Possible reasons to explain the excesses may include unresolved sources or more complicated propagation models. We illustrate that an exponential-cutoff-power-law component with an index of -2.40 and cutoff energy of ~30 TeV is able to account for such excesses.

We study the use of U-Nets in reconstructing the linear dark matter density field and its consequences for constraining cosmological parameters, in particular primordial non-Gaussianity. Our network is able to reconstruct the initial conditions of redshift $z=0$ density fields from N-body simulations with $90\%$ accuracy out to $k \leq 0.4$ h/Mpc, competitive with state-of-the-art reconstruction algorithms at a fraction of the computational cost. We study the information content of the reconstructed $z=0$ density field with a Fisher analysis using the QUIJOTE simulation suite, including non-Gaussian initial conditions. Combining the pre- and post-reconstructed power spectrum and bispectrum data up to $k_{\rm max} = 0.52$ h/Mpc, we find significant improvements on all parameters. Most notably, we find a factor $3.65$ (local), $3.54$ (equilateral) and $2.90$ (orthogonal) improvement on the marginalized errors of $f_{\rm NL}$ as compared to only using the pre-reconstructed data. We show that these improvements can be attributed to a combination of reduced data covariance and parameter degeneracy. The results constitute an important step towards more optimal inference of primordial non-Gaussianity from non-linear scales.

The Macquart relation describes the correlation between the dispersion measure (DM) of fast radio bursts (FRBs) and the redshift $z$ of their host galaxies. The scatter of the Macquart relation is sensitive to the distribution of baryons in the intergalactic medium (IGM) including those ejected from galactic halos through feedback processes. The width of the distribution in DMs from the cosmic web (${\rm DM}_{\rm cosmic}$) is parameterized by a fluctuation parameter $F$, which is related to the cosmic DM variance by $\sigma_{\rm DM}= F z^{-0.5}$. In this work, we present a new measurement of $F$ using 78 FRBs of which 21 have been localized to host galaxies. Our analysis simultaneously fits for the Hubble constant $H_0$ and the DM distribution due to the FRB host galaxy. We find that the fluctuation parameter is degenerate with these parameters, most notably $H_0$, and use a uniform prior on $H_0$ to measure $\log_{10} F > -0.89$ at the $3\sigma$ confidence interval and a new constraint on the Hubble constant $H_0 = 85.3_{-8.1}^{+9.4} \, {\rm km \, s^{-1} \, Mpc^{-1}}$. Using a synthetic sample of 100 localized FRBs, the constraint on the fluctuation parameter is improved by a factor of $\sim 2$. Comparing our $F$ measurement to simulated predictions from cosmological simulation (IllustrisTNG), we find agreement between $0.4 < z < 2$. However, at $z < 0.4$, the simulations underpredict $F$ which we attribute to the rapidly changing extragalactic DM excess distribution at low redshift.

We propose that the dark matter of our universe could be sterile neutrinos which reside within the twin sector of a mirror twin Higgs model. In our scenario, these particles are produced through a version of the Dodelson-Widrow mechanism that takes place entirely within the twin sector, yielding a dark matter candidate that is consistent with X-ray and gamma-ray line constraints. Furthermore, this scenario can naturally avoid the cosmological problems that are typically encountered in mirror twin Higgs models. In particular, if the sterile neutrinos in the Standard Model sector decay out of equilibrium, they can heat the Standard Model bath and reduce the contributions of the twin particles to $N_\mathrm{eff}$. Such decays also reduce the effective temperature of the dark matter, thereby relaxing constraints from large-scale structure. The sterile neutrinos included in this model are compatible with the seesaw mechanism for generating Standard Model neutrino masses.

We present a new divergence-free and well-balanced hybrid FV/FE scheme for the incompressible viscous and resistive MHD equations on unstructured mixed-element meshes in 2 and 3 space dimensions. The equations are split into subsystems. The pressure is defined on the vertices of the primary mesh, while the velocity field and the normal components of the magnetic field are defined on an edge-based/face-based dual mesh in two and three space dimensions, respectively. This allows to account for the divergence-free conditions of the velocity field and of the magnetic field in a rather natural manner. The non-linear convective and the viscous terms are solved at the aid of an explicit FV scheme, while the magnetic field is evolved in a divergence-free manner via an explicit FV method based on a discrete form of the Stokes law in the edges/faces of each primary element. To achieve higher order of accuracy, a pw-linear polynomial is reconstructed for the magnetic field, which is guaranteed to be divergence-free via a constrained L2 projection. The pressure subsystem is solved implicitly at the aid of a classical continuous FE method in the vertices of the primary mesh. In order to maintain non-trivial stationary equilibrium solutions of the governing PDE system exactly, which are assumed to be known a priori, each step of the new algorithm takes the known equilibrium solution explicitly into account so that the method becomes exactly well-balanced. This paper includes a very thorough study of the lid-driven MHD cavity problem in the presence of different magnetic fields. We finally present long-time simulations of Soloviev equilibrium solutions in several simplified 3D tokamak configurations even on very coarse unstructured meshes that, in general, do not need to be aligned with the magnetic field lines.

We develop a method to perform an untargeted Bayesian search for anisotropic gravitational-wave backgrounds that can efficiently and accurately reconstruct the background intensity map. Our method employs an analytic marginalization of the posterior of the spherical-harmonic components of the intensity map, without assuming the background possesses any specific angular structure. The key idea is to realize that the likelihood function is a multivariable Gaussian of the spherical-harmonic components of the energy spectrum of the gravitational-wave background. If a uniform and wide prior of these spherical-harmonic components is prescribed, the marginalized posterior and the Bayes factor can be well approximated by a high-dimensional Gaussian integral. The analytical marginalization allows us to regard the spherical-harmonic components of the intensity map of the background as free parameters, and to construct their individual marginalized posterior distribution in a reasonable time, even though many spherical-harmonic components are required. The marginalized posteriors can, in turn, be used to accurately construct the intensity map of the background. By applying our method to mock data, we show that we can recover precisely the angular structures of various simulated anisotropic backgrounds, without assuming prior knowledge of the relation between the spherical-harmonic components predicted by a given model. Our method allows us to bypass the time-consuming numerical sampling of a high-dimensional posterior, leading to a more model-independent and untargeted Bayesian measurement of the angular structures of the gravitational-wave background.

Putting constraints on a possible Lorentz Invariance Violation (LIV) from astrophysical sources such as gamma-ray bursts (GRBs) is an essential tool for finding evidences of new theories of quantum gravity (QG) that predict energy-dependent speed of light. Such a search has its own difficulties, so usually, the effect of the cosmological model is understudied and the default model is a fixed-parameters $\Lambda$CDM. In this work, we use different astrophysical datasets to study the effect of a number of dark energy models on the LIV constrains. To this end, we combine two public time-delay GRB datasets with the supernovae Pantheon dataset, a number of angular baryonic acoustic oscillations (BAO), the cosmic microwave background (CMB) distance prior and a GRB or quasars dataset. We find for $\alpha$ the expected average value of $\sim 4 \times 10^{-4}$, corresponding to $E_{QG}\ge 10^{17}$ GeV for both time-delay (TD) datasets, with the second one being more sensitive to the cosmological model. We find that the cosmology amounts to at least 20\% deviation in our constraints on the energy. Also interestingly, adding the TD points makes the DE models less-preferable statistically and shifts the value of the parameter $c/(H_0 r_d)$ down, making it smaller than the expected value. We see that possible LIV measurements depend critically on the transparency of the assumptions behind the published data with respect to cosmology and that taking it into account may be important contribution in the case of possible detection.

Neutron stars are one of the most mysterious wonders in the Universe. Their extreme densities hint at new and exotic physics at work within. Gravitational waves could be the key to unlocking their secrets. In particular, a first detection of gravitational waves from rapidly-spinning, deformed neutron stars could yield new insights into the physics of matter at extreme densities and under strong gravity. Once a first detection is made, a critical challenge will be to robustly extract physically interesting information from the detected signals. In this essay, we describe initial research towards answering this challenge, and thereby unleashing the full power of gravitational waves as an engine for the discovery of new physics.

Onboard electrostatic suspension inertial sensors are important applications for gravity satellites and space gravitational wave detection missions, and it is important to suppress noise in the measurement signal. Due to the complex coupling between the working space environment and the satellite platform, the process of noise generation is extremely complex, and traditional noise modeling and subtraction methods have certain limitations. With the development of deep learning, applying it to high-precision inertial sensors to improve the signal-to-noise ratio is a practically meaningful task. Since there is a single noise sample and unknown true value in the measured data in orbit, odd-even sub-samplers and periodic sub-samplers are designed to process general signals and periodic signals, and adds reconstruction layers consisting of fully connected layers to the model. Experimental analysis and comparison are conducted based on simulation data, GRACE-FO acceleration data and Taiji-1 acceleration data. The results show that the deep learning method is superior to traditional data smoothing processing sol

The Electron, Proton and Alpha Monitor, EPAM, located at the L1 Position approximately 1-million miles from the earth in the direction of the sun, was designed to detect fluctuations in solar output through counting the numbers of various particles hitting the detector. The EPAM detector is part of an early warning system that can alert the earth to coronal mass ejection events that can damage our electronic grids and satellite equipment. EPAM gives a real-time estimate of changes in the local solar magnetic field directed towards the earth, recorded in the fluctuations of solar particles being ejected. This paper presents an analysis of fluctuations in data taken by the Geological Survey of Israel, GSI, compared to the changes in detected numbers of protons as seen by EPAM. Surprisingly, the GSI and EPAM detectors show an unexpected correlation between the variation in count rate detected by the GSI detectors and an increased numbers of protons seen at EPAM; well above statistical significance of 5-sigma, indicating a non-random connection between the data sets. The statistically significant overlap between data taken by these two detectors, subject to very different conditions, may hint at a Primakoff mechanism whereby exotic particles, e.g. galactic Dark Matter, couple through magnetic fields to both photons and even nuclei. This work builds on an earlier paper on the observations of Radon decay and their implications for particle physics.