Daily Astro papers on the arXiv in all categories.
The dynamical state of the globular clusters Rup 106 and IC 4499 by Giacomo Beccari et al. on Wednesday 30 November The dynamical evolution of globular clusters is theoretically described by a series of well known events typical of N-body systems. Still, the identification of observational signatures able to empirically describe the stage of dynamical evolution of a stellar system of the density typical of a globular cluster, represents a challenge. In this paper we study the dynamical age of the globular clusters Rup 106 and IC 4499. To this aim, we study the radial distribution of the Blue Straggler Stars via the A+ parameter and of the slope of the Main Sequence Mass Function. Both tracers show that Rup 106 and IC 4499 are dynamically young clusters where dynamical friction has just started to segregate massive stars towards the clusters' centre. In fact, we observe that the Blue Straggler stars are more centrally concentrated in both clusters than the reference population. On the same line, we find that in both cases the slope of the mass function significantly decreases as a function of the cluster-centric distances. This result provides additional support for the use of the the radial distribution of the blue stragglers as a powerful observationally convenient indicator of the cluster dynamical age. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.12884v2
Numerical Modeling of Energetic Electron Acceleration, Transport, and Emission in Solar Flares: Connecting Loop-top and Footpoint Hard X-Ray Sources by Xiangliang Kong et al. on Wednesday 30 November The acceleration and transport of energetic electrons during solar flares is one of the outstanding topics in solar physics. Recent X-ray and radio imaging and spectroscopy observations have provided diagnostics of the distribution of nonthermal electrons and suggested that, in certain flare events, electrons are primarily accelerated in the loop-top and likely experience trapping and/or scattering effects. By combining the focused particle transport equation with magnetohydrodynamic (MHD) simulations of solar flares, we present a macroscopic particle model that naturally incorporates electron acceleration and transport. Our simulation results indicate that the physical processes such as turbulent pitch-angle scattering can have important impacts on both electron acceleration in the loop-top and transport in the flare loop, and their influences are highly energy dependent. A spatial-dependent turbulent scattering with enhancement in the loop-top can enable both efficient electron acceleration to high energies and transport of abundant electrons to the footpoints. We further generate spatially resolved synthetic hard X-ray (HXR) emission images and spectra, revealing both the loop-top and footpoint HXR sources. Similar to the observations, we show that the footpoint HXR sources are brighter and harder than the loop-top HXR source. We suggest that the macroscopic particle model provides new insights into understanding the connection between the observed loop-top and footpoint nonthermal emission sources by combining the particle model with dynamically evolving MHD simulations of solar flares. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15333v1
Electromagnetic draping of merging neutron stars by Maxim Lyutikov. on Wednesday 30 November We first derive a set of equations describing general stationary configurations of relativistic force-free plasma, without assuming any geometric symmetries. We then demonstrate that electromagnetic interaction of merging neutron stars is necessarily dissipative due to the effect of electromagnetic draping - creation of dissipative regions near the star (in the single-magnetized case) or at the magnetospheric boundary (in the double-magnetized case). Our results indicate that even in the single magnetized case we expect that relativistic jets (or ``tongues'') are produced, with correspondingly beamed emission pattern. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.14433v2
Stellar space weather effects on potentially habitable planets by A. A. Vidotto. on Wednesday 30 November Stellar activity can reveal itself in the form of radiation (eg, enhanced X-ray coronal emission, flares) and particles (eg, winds, coronal mass ejections). Together, these phenomena shape the space weather around (exo)planets. As stars evolve, so do their different forms of activity -- in general, younger solar-like stars have stronger winds, enhanced flare occurrence and likely more frequent coronal mass ejections. Altogether, these effects can create harsher particle and radiation environments for habitable-zone planets, in comparison to Earth, in particular at young ages. In this article, I will review some effects of these harsher environments on potentially habitable exoplanets. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15396v1
The AstroSat UV Deep Field North: the far and near ultraviolet photometric catalog by Chayan Mondal et al. on Wednesday 30 November We present deep UV imaging observations of the Great Observatories Origins Survey Northern (GOODS-N) field with AstroSat/UVIT (AstroSat UV Deep Field north - AUDFn), using one far-UV (F154W, 34.0 kilosec) and two near-UV filters (N242W, 19.2 kilosec; N245M, 15.5 kilosec). The nature of the UV sky background was explored across the UVIT field and a global mean and rms was estimated for each filter. We reach 3$sigma$ detection limits of $m_{rm AB}$ $sim$ 27.35 mag, 27.28 mag and 27.02 mag for a point source in the F154W, N242W and N245M bands respectively. The 50% completeness limits of the FUV and NUV images are $m_{rm AB}=$ 26.40 mag and 27.05 mag respectively. We constructed PSFs for each band and estimated their FWHM, which were found to be almost the same: 1.18" in F154W, 1.11" in N242W, and 1.24" in N245M. We used SExtractor to separately identify sources in the FUV and NUV filters and produce the UV source catalog of the entire AUDFn field. The source count slope estimated in FUV and NUV is 0.57 dex mag$^{-1}$ (between 19 - 25 mag) and 0.44 dex mag$^{-1}$ (between 18 - 25 mag), respectively. The catalog contains 6839 and 16171 sources (brighter than the 50% completeness limit) in the FUV and NUV, respectively. Our FUV and NUV flux measurements of the identified sources complement existing multi-band data in the GOODS-N field, and enable us to probe rest-frame FUV properties of galaxies at redshift $z < 1$ and search for candidate Lyman continuum leakers at redshift $z > 0.97$. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16923v1
Status and prospects of the CORSIKA 8 air shower simulation framework by Alexander Sandrock. on Wednesday 30 November The Fortran-versions of the CORSIKA air shower simulation code have been at the core of simulations for many astroparticle physics experiments for the last 30 years. Having grown over decades into an ever more complex software, maintainability of CORSIKA has become increasingly difficult, though its performance is still excellent. In 2018, therefore a complete rewrite of CORSIKA has begun in modern modular C++. Today, CORSIKA 8 has reached important milestones with a full-fledged implementation of both the hadronic and electromagnetic cascades, the ability to simulate radio and Cherenkov-light emission from air showers and an unprecedented flexibility to configure simulation media and their geometries. This presentation will discuss the current status of CORSIKA 8, highlight the new possibilities already available, and future prospects of this new air shower simulation framework. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16057v1
RSD measurements from BOSS galaxy power spectrum using the halo perturbation theory model by Byeonghee Yu et al. on Wednesday 30 November We present growth of structure constraints from the cosmological analysis of the power spectrum multipoles of SDSS-III BOSS DR12 galaxies. We use the galaxy power spectrum model of Hand et al. (2017), which decomposes the galaxies into halo mass bins, each of which is modeled separately using the relations between halo biases and halo mass. The model combines Eulerian perturbation theory and halo model calibrated on $N$-body simulations to model the halo clustering. In this work, we also generate the covariance matrix by combining the analytic disconnected part with the empirical connected part: we smooth the connected component by selecting a few principal components and show that it achieves good agreement with the mock covariance. We find tight constraints on $fsigma_8$: $fsigma_8(z_{mathrm{eff}}=0.38)=0.489 pm 0.036$ and $fsigma_8(z_{mathrm{eff}}=0.61)=0.455 pm 0.026$ at $k_{mathrm{max}} = 0.2 h$Mpc$^{-1}$, in good agreement with Planck amplitude. This corresponds to $S_8 = 0.821 pm 0.037$ or an overall amplitude error of 4%, within 0.3 sigma of Planck's $S_8 = 0.832 pm 0.013$. We discuss the sensitivity of cosmological parameter estimation to the choice of scale cuts, covariance matrix, and the inclusion of hexadecapole $P_4(k)$. We show that with $k_{mathrm{max}} = 0.4 h$Mpc$^{-1}$ the constraints improve considerably to an overall 2.7% amplitude error (with $S_8 = 0.786 pm 0.021$), but there is some evidence of model misspecification on MultiDark-PATCHY mocks. Choosing $k_{mathrm{max}}$ consistently and reliably remains the main challenge of RSD analysis methods. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16794v1
Star & planet formation: Upcoming opportunities in the space-based infrared by Floris van der Tak. on Wednesday 30 November While ALMA and JWST are revolutionizing our view of star and planet formation with their unprecedented sensitivity and resolution at submillimeter and near-IR wavelengths, many outstanding questions can only be answered with observations in the thermal (mid- and far-) infrared domain. Many of these questions require space-based observations, to achieve the necessary sensitivity and/or wavelength coverage. In particular, how do interstellar clouds develop filamentary structures and dense cores? What are the masses and luminosities of objects at the earliest stages of star formation? What are the gas masses of planet-forming disks, and how do these disks disperse during planet formation? How is refractory and volatile material distributed within the disks, and how does this evolve with time? This article reviews how upcoming and planned balloon-borne and space-based telescopes for the mid- and far-infrared will address these questions, and outlines which further missions will be needed beyond 2030, when the ELTs will be in full operation. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15149v2
How has the solar wind evolved to become what it is today? by A. A. Vidotto. on Wednesday 30 November In this contribution, I briefly review the long-term evolution of the solar wind (its mass-loss rate), including the evolution of observed properties that are intimately linked to the solar wind (rotation, magnetism and activity). I also briefly discuss implications of the evolution of the solar wind on the evolving Earth. I argue that studying exoplanetary systems could open up new avenues for progress to be made in our understanding of the evolution of the solar wind. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15400v1
SILCC VII -- Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities by Tim-Eric Rathjen et al. on Wednesday 30 November We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities $Sigma_mathrm{gas} =$ 10, 30, 50, and 100 $mathrm{M_odot,pc^{-2}}$. The SILCC project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionising UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between $Sigma_mathrm{gas}$ and the star formation rate surface density $Sigma_mathrm{SFR}$. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher $Sigma_mathrm{gas}$ and the mass loading factors of the CR-supported outflows are of order unity independent of $Sigma_mathrm{SFR}$. Similar to observations, vertical velocity dispersions of the warm ionised medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $sigma_mathrm{z} propto Sigma_mathrm{SFR}^a$, with $a sim 0.20$. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor $sim 2.2$ higher than that of the CNM, in agreement with local observations. For $Sigma_mathrm{SFR} gtrsim 1.5 times 10^{-2},mathrm{M}_odot,mathrm{yr}^{-1},mathrm{kpc}^{-2}$ the WIM motions become supersonic. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15419v1
Simultaneous gas accretion onto a pair of giant planets: Impact on their final mass and on the protoplanetary disk structure by Camille Bergez-Casalou et al. on Wednesday 30 November Several planetary systems are known to host multiple giant planets. However, when two giant planets are accreting from the same disk, it is unclear what effect the presence of the second planet has on the gas accretion process of both planets. In this paper we perform long-term 2D isothermal hydrodynamical simulations (over more than 0.5 Myrs) with the FARGO-2D1D code, considering two non-migrating planets accreting from the same gaseous disk. We find that the evolution of the planets' mass ratio depends on gap formation. However, in all cases, when the planets start accreting at the same time, they end up with very similar masses (0.9 $
ISPY-NACO Imaging Survey for Planets around Young stars The demographics of forming planets embedded in protoplanetary disks by Gabriele Cugno et al. on Wednesday 30 November We present the statistical analysis of a subsample of 45 young stars surrounded by protoplanetary disks (PPDs). This is the largest imaging survey uniquely focused on PPDs to date. Our goal is to search for young forming companions embedded in the disk material and to constrain their occurrence rate in relation to the formation mechanism. We used principal component analysis based point spread function subtraction techniques to reveal young companions forming in the disks. We calculated detection limits for our datasets and adopted a black-body model to derive temperature upper limits of potential forming planets. We then used Monte Carlo simulations to constrain the population of forming gas giant companions and compare our results to different types of formation scenarios. Our data revealed a new binary system (HD38120) and a recently identified triple system with a brown dwarf companion orbiting a binary system (HD101412), in addition to 12 known companions. Furthermore, we detected signals from 17 disks, two of which (HD72106 and TCrA) were imaged for the first time. We reached median detection limits of L =15.4 mag at 2.0 arcsec, which were used to investigate the temperature of potentially embedded forming companions. We can constrain the occurrence of forming planets with semi-major axis a in [20 - 500] au and Teff in [600 - 3000] K, in line with the statistical results obtained for more evolved systems from other direct imaging surveys. The NaCo-ISPY data confirm that massive bright planets accreting at high rates are rare. More powerful instruments with better sensitivity in the near- to mid-infrared are likely required to unveil the wealth of forming planets sculpting the observed disk substructures. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15434v1
Observations of planetary winds and outflows by Leonardo A. Dos Santos. on Wednesday 30 November We have recently hit the milestone of 5,000 exoplanets discovered. In stark contrast with the Solar System, most of the exoplanets we know to date orbit extremely close to their host stars, causing them to lose copious amounts of gas through atmospheric escape at some stage in their lives. In some planets, this process can be so dramatic that they shrink in timescales of a few million to billions of years, imprinting features in the demographics of transiting exoplanets. Depending on the transit geometry, ionizing conditions, and atmospheric properties, a planetary outflow can be observed using transmission spectroscopy in the ultraviolet, optical or near-infrared. In this review, we will discuss the main techniques to observe evaporating exoplanets and their results. To date, we have evidence that at least 28 exoplanets are currently losing their atmospheres, and the literature has reported at least 42 non-detections. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16243v1
A Bayesian approach to RFI mitigation by S. A. K. Leeney et al. on Wednesday 30 November Radio Frequency Interference (RFI) is an endemic problem in radio astronomy. Information in contaminated frequency channels is lost and it can lead to significant systematic error if not properly modelled. In this paper we propose a RFI mitigation methodology that takes a Bayesian approach, where contaminated data is both flagged and managed as part of a single step fitting process. To the authors knowledge, our approach is first-of-its-kind. We provide a proof-of-concept of our methods which are also tested and shown to be effective on both a simple toy model and when incorporated into the Bayesian data analysis pipeline for REACH (a modern cosmological radio experiment). arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15448v1
EIGER III JWST NIRCam observations of the ultra-luminous high-redshift quasar J0100+2802 by Anna-Christina Eilers et al. on Wednesday 30 November We present the first rest-frame optical spectrum of a high-redshift quasar observed with JWST/NIRCam in Wide Field Slitless (WFSS) mode. The observed quasar, J0100+2802, is the most luminous quasar known at $z>6$. We measure the mass of the central supermassive black hole (SMBH) by means of the rest-frame optical H$beta$ emission line, and find consistent mass measurements of the quasar's SMBH of $M_bulletapprox10^{10},M_odot$ when compared to the estimates based on the properties of rest-frame UV emission lines CIV and MgII, which are accessible from ground-based observatories. To this end, we also present a newly reduced rest-frame UV spectrum of the quasar observed with X-Shooter/VLT and FIRE/Magellan for a total of 16.8 hours. We readdress the question whether this ultra-luminous quasar could be effected by strong gravitational lensing making use of the diffraction limited NIRCam images in three different wide band filters (F115W, F200W, F356W), which improves the achieved spatial resolution compared to previous images taken with the Hubble Space Telescope by a factor of two. We do not find any evidence for a foreground deflecting galaxy, nor for multiple images of the quasar, and determine the probability for magnification due to strong gravitational lensing with image separations below the diffraction limit of $Deltathetalesssim 0.05''$ to be $lesssim 2.2times 10^{-3}$. Our observations therefore confirm that this quasar hosts a ten billion solar mass black hole less than $1$ Gyr after the Big Bang, which is challenging to explain with current black hole formation models. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16261v1
Assessing the physical reality of Milky Way open cluster candidates by A. E. Piatti et al. on Wednesday 30 November We report results on the analysis of eleven new Milky Way open cluster candidates, recently discovered from the detection of stellar overdensities in the Vector Point diagram, by employing extreme deconvolution Gaussian mixture models. We treated these objects as real open clusters and derived their fundamental properties with their associated intrinsic dispersions by exploring the parameter space through the minimization of likelihood functions on generated synthetic colour-magnitude diagrams (CMDs). The intrinsic dispersions of the resulting ages turned out to be much larger than those usually obtained for open clusters. Indeed, they resemble those of ages and metallicities of composite star field populations. We also traced their stellar number density profiles and mass functions, derived their total masses, Jacobi and tidal radii, which helped us as criteria while assessing their physical nature as real open clusters. Because the eleven candidates show a clear gathering of stars in the proper motion plane and some hint for similar distances, we concluded that they are possibly sparse groups of stars. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15483v1
The Complex Structure of the Bulge of M31 by Denis Leahy et al. on Wednesday 30 November The bulge of M31 is of interest in the context of the nature of galactic bulges and how their structure relates to bulge formation mechanisms and their subsequent evolution. With the UVIT instrument on AstroSat, we have observed the bulge of M31 in five far ultraviolet (FUV) and near ultraviolet (NUV) filters at 1" spatial resolution. Models for the luminosity distribution of the bulge are constructed using the UVIT data and the galaxy image-fitting algorithm GALFIT. We fit the bulge without nuclear region with a Sersic function the five images and find Sersic indices ($simeq2.1$ to 2.5) similar to previous studies but smaller $R_e$ values ($simeq0.5$ to 0.6 kpc). When the images include the nuclear region, a multicomponent model is used to find the best-fit. We use an 8-component fit for the FUV 148nm image, which has the highest sensitivity. The other images (169 to 279 nm) are fit with 4-component models. The dust lanes in the bulge region are recovered in the residual image after subtraction of the bright bulge light using the multicomponent model. The dust lanes show that M31's nuclear spiral is visible in absorption at NUV and FUV wavelengths. The bulge images show boxy contours in all five UVIT wavebands, which is confirmed by fitting using GALFIT. The Sersic indices of $sim$2.1-2.5 are intermediate between expected values for a classical bulge and for a pseudobulge. The boxiness of the bulge provides further evidence that M31's bulge has contributions from a classical bulge and a pseudobulge. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16311v1
FIREBall-2: flight preparation of a proven balloon payload to image the intermediate redshift circumgalactic medium by Vincent Picouet et al. on Wednesday 30 November FIREBall-2 is a stratospheric balloon-borne 1-m telescope coupled to a UV multi-object slit spectrograph designed to map the faint UV emission surrounding z~0.7 galaxies and quasars through their Lyman-alpha line emission. This spectro-imager had its first launch on September 22nd 2018 out of Ft. Sumner, NM, USA. Because the balloon was punctured, the flight was abruptly interrupted. Instead of the nominal 8 hours above 32 km altitude, the instrument could only perform science acquisition for 45 minutes at this altitude. In addition, the shape of the deflated balloon, combined with a full Moon, revealed a severe off-axis scattered light path, directly into the UV science detector and about 100 times larger than expected. In preparation for the next flight, and in addition to describing FIREBall-2's upgrade, this paper discusses the exposure time calculator (ETC) that has been designed to analyze the instrument's optimal performance (explore the instrument's limitations and subtle trade-offs). arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15491v1
Direct Statistical Constraints on the Natal Kick velocity of a Black Hole in an X-ray Quiet Binary by Sharan Banagiri et al. on Wednesday 30 November In recent years, a handful of ``dark" binaries have been discovered with a non-luminous compact object. Astrometry and radial velocity measurements of the bright companion allow us to measure the post-supernova orbital elements of such a binary. In this paper, we develop a statistical formalism to use such measurements to infer the pre-supernova orbital elements, and the natal kick imparted by the supernova (SN). We apply this formalism to the recent discovery of an X-ray quiet binary with a black hole, VFTS 243, in the Large Magellanic Cloud. Assuming an isotropic, Maxwellian distribution on natal kicks and using broad agnostic mass priors, we find that kick velocity can be constrained to less than $V_k < 56$ km/s and the dispersion of the kick distribution to $sigma_k < 66 $ km/s at 90 % confidence. Assuming that there was a natal kick, we find that at least about $0.6 M_{odot}$ was lost during the supernova with 90 % confidence. The pre-SN orbital separation is found to be robustly constrained to be around $0.43$ AU. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16361v1
Continuous gravitational wave emission from neutron stars with pinned superfluids in the core by Brynmor Haskell et al. on Wednesday 30 November We investigate the effect of a pinned superfluid component on the gravitational wave emission of a rotating neutron star. Pinning of superfluid vortices to the flux-tubes in the outer core (where the protons are likely to form a type-II superconductor) is a possible mechanism to sustain long-lived and non-axisymmetric neutron currents in the interior, that break the axial symmetry of the unperturbed hydrostatic configuration. We consider pinning-induced perturbations to a stationary corotating configuration, and determine upper limits on the strength of gravitational wave emission due to the pinning of vortices with a strong toroidal magnetic field of the kind predicted by recent magneto-hydrodynamic simulations of neutron star interiors. We estimate the contributions to gravitational wave emission from both the mass and current multipole generated by the pinned vorticity in the outer core, and find that the mass quadrupole can be large enough for gravitational waves to provide the dominant spindown torque in millisecond pulsars. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15507v1
Updating constraints on phantom crossing f T gravity by F. B. M. dos Santos. on Wednesday 30 November We establish constraints on $f(T)$ gravity by considering the possibility of a scenario that supports a phantom crossing of the equation of state parameter $w_{DE}$. After determining the viable parameter space of the model, while checking the impact on the background dynamics, we perform an analysis to obtain constraints on cosmological parameters and determine the viability of this scenario. To this end, we use combined data sets from cosmic chronometers (CC), baryonic acoustic oscillations (BAO), redshift space measurements (RSD), cosmic microwave background (CMB) and Big Bang nucleosynthesis (BBN) priors and Type Ia supernovae (SNe) measurements from the Pantheon set, in which the impact on the absolute magnitude is also considered. It is found that a state where a phantom crossing of $w_{DE}$ happens is favored by data, and while the $f(T)$ model is equivalent to the $Lambda$CDM one by statistical criteria when CC+SNe+BAO+BBN data is considered, a considerable tension is found when CMB data is added. We also find that the Hubble tension is alleviated in the $f(T)$ model, at the same time that it does not worsen the growth one, indicating a possibility of the scenario as an option to address the current cosmic tensions. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16370v1
Orientations of DM Halos in FIRE-2 Milky Way-mass Galaxies by Jay Baptista et al. on Wednesday 30 November The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo's orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30--90 degrees depending on the individual halo's formation history. In identical simulations using a model of self-interacting DM with $sigma = 1 , mathrm{cm}^2 , mathrm{g}^{-1}$, the halo remains aligned with the stellar disk out to $sim$200--400 kpc. Interactions with massive satellites ($M gtrsim 4 times 10^{10} , rm{M_odot}$ at pericenter; $M gtrsim 3.3 times 10^{10} , rm{M_odot}$ at infall) affect the orientation of the halo significantly, aligning the halo's major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16382v1
GECAM Localization of High Energy Transients and the Systematic Error by Yi Zhao et al. on Wednesday 30 November Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a pair of microsatellites (i.e. GECAM-A and GECAM-B) dedicated to monitoring gamma-ray transients including gravitational waves high-energy electromagnetic counterparts, Gamma-ray Bursts, Soft Gamma-ray Repeaters, Solar Flares and Terrestrial Gamma-ray Flashes. Since launch in December 2020, GECAM-B has detected hundreds of astronomical and terrestrial events. For these bursts, localization is the key for burst identification and classification as well as follow-up observations in multi-wavelength. Here, we propose a Bayesian localization method with Poisson data with Gaussian background profile likelihood to localize GECAM bursts based on the burst counts distribution in detectors with different orientations. We demonstrate that this method can work well for all kinds of bursts, especially for extremely short ones. In addition, we propose a new method to estimate the systematic error of localization based on a confidence level test, which can overcome some problems of the existing method in literature. We validate this method by Monte Carlo simulations, and then apply it to a burst sample with accurate location and find that the mean value of the systematic error of GECAM-B localization is $sim 2.5^{circ}$. By considering this systematic error, we can obtain a reliable localization probability map for GECAM bursts. Our methods can be applied to other gamma-ray monitors. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15570v1
On the sociology and hierarchy of voids: a study of seven CAVITY nearby galaxy voids and their dynamical CosmicFlows-3 environment by H. M. Courtois et al. on Wednesday 30 November Context. The present study addresses a key question for our understanding of the relation between void galaxies and their environment: the relationship between luminous and dark matter in and around voids. Aims. To explore how empty of matter local Universe voids are, we study the full (dark+luminous) matter content of seven nearby cosmic voids that are fully contained within the CosmicFlows-3 volume. Methods. The cosmic voids matter density profiles are independently obtained using two different methods. They are built on one hand from the galaxy redshift space 2 points-correlation function and, on the other hand, using peculiar velocity gradients from the CosmicFlows-3 dataset. Results. The results are noticeable since when using the redshift survey, all voids show a radial positive gradient of galaxies, while based on the dynamical analysis, only three of these voids display a clear underdensity of matter in their center. Conclusions. It is the first time such a detailed observational analysis of voids is conducted, showing that void emptiness should be derived from dynamical information. Yet, from this limited study, the Hercules void is the best candidate for a local Universe pure "pristine volume" expanding in 3 directions with no dark matter located in that void. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16388v1
Hubble Space Telescope transmission spectroscopy for the temperate sub-Neptune TOI-270d: a possible hydrogen-rich atmosphere containing water vapour by Thomas Mikal-Evans et al. on Wednesday 30 November TOI-270d is a temperate sub-Neptune discovered by the Transiting Exoplanet Survey Satellite (TESS) around a bright (J=9.1mag) M3V host star. With an approximate radius of 2RE and equilibrium temperature of 350K, TOI-270d is one of the most promising small exoplanets for atmospheric characterisation using transit spectroscopy. Here we present a primary transit observation of TOI-270d made with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) spectrograph across the 1.126-1.644 micron wavelength range, and a 95% credible upper limit of $8.2 times 10^{-14}$ erg s$^{-1}$ cm$^{-2}$ A$^{-1}$ arcsec$^{-2}$ for the stellar Ly-alpha emission obtained using the Space Telescope Imaging Spectrograph (STIS). The transmission spectrum derived from the TESS and WFC3 data provides evidence for molecular absorption by a hydrogen-rich atmosphere at 4-sigma significance relative to a featureless spectrum. The strongest evidence for any individual absorber is obtained for H2O, which is favoured at 3-sigma significance. When retrieving on the WFC3 data alone and allowing for the possibility of a heterogeneous stellar brightness profile, the detection significance of H2O is reduced to 2.8-sigma. Further observations are therefore required to robustly determine the atmospheric composition of TOI-270d and assess the impact of stellar heterogeneity. If confirmed, our findings would make TOI-270d one of the smallest and coolest exoplanets to date with detected atmospheric spectral features. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15576v1
Flares detected in ALMA single-dish images of the Sun by I. Skokić et al. on Wednesday 30 November The (sub)millimeter radiation of solar flares is poorly understood. Without spatial resolution, it cannot be compared easily to flare emissions in other wavelengths. The Atacama Large Millimeter-submillimeter Array (ALMA) offers sufficient resolution for the first time. However, used as an interferometer, its field of view is smaller than an active region and ALMA cannot observe on demand when a flare occurs. We use readily available large scale single-dish ALMA observations of solar millimeter flares and compare them to well-known features observed in other wavelengths. The properties of these other flare emissions, correlating in space and time, may then be used to interpret the millimeter brightenings and vice versa. The aim is to obtain reliable associations, limited by the time and space resolution of single-dish observations. We collected ALMA observations at 3 mm and 1 mm and searched for millimeter brightenings during times given in a flare catalog. We found five events with 9 or more images that can be used for comparison in time and space. The millimeter brightenings are associated with a variety of flare features in cool (H$alpha$, 30.4 nm), intermediate (17.1 nm), and hot (9.4 nm) lines. In several cases, the millimeter brightening peaked at the footpoint of a hot flare loop. In other cases the peak coincided with the top or footpoint of an active H{alpha} filament. We found correlations also with post-flare loops and tops of a hot loop, and in some cases to no features at all. The wide field of view provided by the single-dish observations allowed for completely overviewing the flare activity in millimeter waves for the first time. The associated phenomena often changed during the flare in type and location, and may explain the sometimes bewildering behavior of millimeter flare emissions observed previously without spatial resolution. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16935v1
Effects of the Central Mass Concentration on Bar Formation in Disk Galaxies by Dajeong Jang et al. on Wednesday 30 November While bars are common in disk galaxies, their formation conditions are not well understood. We use $N$-body simulations to study bar formation and evolution in isolated galaxies consisting of a stellar disk, a classical bulge, and a dark halo. We consider 24 galaxy models that are similar to the Milky Way but differ in the mass and compactness of the classical bulge and halo concentration. We find that the bar formation requires $(Q_{T,text{min}}/1.2)^2 + (text{CMC}/0.05)^2lesssim 1$, where $Q_{T,text{min}}$ and CMC refers to the minimum value of the Toomre stability parameter and the central mass concentration, respectively. Bars tend to be stronger, longer, and rotate slower in galaxies with a less massive and less compact bulge and halo. All bars formed in our models correspond to slow bars. A model with the bulge mass of $sim10$--$20$% of the disk under a concentrated halo produces a bar similar to the Milky-Way bar. We discuss our findings in relation to other bar formation criteria suggested by previous studies. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16816v1
Determination of 1929 Asteroid Rotation Periods from WISE Data by Adrian L. H. Lam et al. on Wednesday 30 November We used 22 $mu$m (W4) Wide-field Infrared Survey Explorer (WISE) observations of 4420 asteroids to analyze lightcurves and determined spin period estimates for 1929 asteroids. We fit second-order Fourier models at a large number of trial frequencies to the W4 data and analyzed the resulting periodograms. We initially excluded rotational frequencies exceeding 7.57 rotations per day (P < 3.17 hr), which are not sampled adequately by WISE, and periods that exceed twice the WISE observation interval, which is typically 36 hr. Three solutions accurately capture the vast majority of the rotational frequencies in our sample: the best-fit frequency and its mirrors around 3.78 and 7.57 rotations per day. By comparing our solutions to a high-quality control group of 752 asteroid spin periods, we found that one of our solutions is accurate (within 5%) in 88% of the cases. The best-fit, secondary, and tertiary solutions are accurate in 55%, 27%, and 6% of the cases, respectively. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16409v1
A Bayesian approach for torque modelling of BeXRB pulsars with application to super-Eddington accretors by A. S. Karaferias et al. on Wednesday 30 November In this study we present a method to estimate posterior distributions for standard accretion torque model parameters and binary orbital parameters for X-ray binaries using a nested sampling algorithm for Bayesian Parameter Estimation. We study the spin evolution of two Be X-ray binary systems in the Magellanic Clouds, RX J0520.5-6932 and RX J0209-7427, during major outbursts, in which they surpassed the Eddington-limit. Moreover, we apply our method to the recently discovered Swift J0243.6+6124; the only known Galactic pulsating ultra-luminous X-ray source. This is an excellent candidate for studying the disc evolution at super-Eddington accretion rates, for its luminosity span several orders of magnitude during its outburst, with a maximum $L_{rm X}$ that exceeded the Eddington limit by a factor of $sim 10$. Our method, when applied to RX J0520.5-6932 and RX J0209-7427, is able to identify the more favourable torque model for each system, while yielding meaningful ranges for the NS and orbital parameters. Our analysis for Swift J0243.6+6124 illustrates that, contrary to the standard torque model predictions, the magnetospheric radius and the Alfv'en radius are not proportional to each other when surpassing the Eddington limit. Reported distance estimates of this source range between 5 and 7 kpc. Smaller distances require non-typical neutron star properties (i.e. mass and radius) and possibly lower radiative efficiency of the accretion column. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16079v1
Gravity in the Local Universe : density and velocity fields using CosmicFlows-4 by H. M. Courtois et al. on Wednesday 30 November This article publicly releases three-dimensional reconstructions of the local Universe gravitational field below z=0.8 that were computed using the full catalogue CosmicFlows-4 of 56,000 galaxy distances and its sub-sample of 1,008 type Ia supernovae distances. The article also provides some first CF4 measurements of the growth rate of structure using the pairwise correlation of peculiar velocities fsigma8 = 0.44(+/-0.01) and of the bulk flow in the Local Universe of 200+/-88 kms-1 at distance 300 h-1Mpc. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16390v1
Exoplanet Detection by Machine Learning with Data Augmentation by Koray AydoÄŸan. on Wednesday 30 November It has recently been demonstrated that deep learning has significant potential to automate parts of the exoplanet detection pipeline using light curve data from satellites such as Kepler cite{borucki2010kepler} cite{koch2010kepler} and NASA's Transiting Exoplanet Survey Satellite (TESS) cite{ricker2010transiting}. Unfortunately, the smallness of the available datasets makes it difficult to realize the level of performance one expects from powerful network architectures. In this paper, we investigate the use of data augmentation techniques on light curve data from to train neural networks to identify exoplanets. The augmentation techniques used are of two classes: Simple (e.g. additive noise augmentation) and learning-based (e.g. first training a GAN cite{goodfellow2020generative} to generate new examples). We demonstrate that data augmentation has a potential to improve model performance for the exoplanet detection problem, and recommend the use of augmentation based on generative models as more data becomes available. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15577v1
Palatini formulation for gauge theory: implications for slow-roll inflation by Syksy Rasanen et al. on Wednesday 30 November We consider a formulation of gauge field theory where the gauge field $A_alpha$ and the field strength $F_{alphabeta}$ are independent variables, as in the Palatini formulation of gravity. For the simplest gauge field action, this is known to be equivalent to the usual formulation. We add non-minimal couplings between $F_{alphabeta}$ and a scalar field, solve for $F_{alphabeta}$ and insert it back into the action. This leads to modified gauge field and scalar field terms. We consider slow-roll inflation and show that because of the modifications to the scalar sector, adding higher order terms to the inflaton potential does not spoil its flatness, unlike in the usual case. Instead they make the effective potential closer to quadratic. The modifications also solve the problem that Higgs inflation in the Palatini formulation is sensitive to higher order terms. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15584v1
Observational studies of high-mass star formation by Igor I. Zinchenko. on Wednesday 30 November We present a review of observational studies of high-mass star formation, based mainly on our own research. It includes surveys of high-mass star-forming regions in various molecular lines and in continuum, investigations of filamentary infrared dark clouds, which represent the earliest phases of massive star formation, detailed studies of individual high-mass star-forming regions, dense cores and disks harboring massive (proto)stars, and associated outflows. Chemistry in these regions is discussed, too. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15586v1
Line-of-sight effects in strong gravitational lensing by Pierre Fleury et al. on Wednesday 30 November While most strong-gravitational-lensing systems may be roughly modelled by a single massive object between the source and the observer, in the details all the structures near the light path contribute to the observed images. These additional contributions, known as line-of-sight effects, are non-negligible in practice. This article proposes a new theoretical framework to model the line-of-sight effects, together with very promising applications at the interface of weak and strong lensing. Our approach relies on the dominant-lens approximation, where one deflector is treated as the main lens while the others are treated as perturbations. The resulting framework is technically simpler to handle than the multi-plane lensing formalism, while allowing one to consistently model any sub-critical perturbation. In particular, it is not limited to the usual external-convergence and external-shear parameterisation. As a first application, we identify a specific notion of line-of-sight shear that is not degenerate with the ellipticity of the main lens, and which could thus be extracted from strong-lensing images. This result supports and improves the recent proposal that Einstein rings might be powerful probes of cosmic shear. As a second application, we investigate the distortions of strong-lensing critical curves under line-of-sight effects, and more particularly their correlations across the sky. We find that such correlations may be used to probe, not only the large-scale structure of the Universe, but also the dark-matter halo profiles of strong lenses. This last possibility would be a key asset to improve the accuracy of the measurement of the Hubble-Lema^itre constant via time-delay cosmography. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2104.08883v5
Determination of 1929 Asteroid Rotation Periods from WISE Data by Adrian L. H. Lam et al. on Wednesday 30 November We used 22 $mu$m (W4) Wide-field Infrared Survey Explorer (WISE) observations of 4420 asteroids to analyze lightcurves and determined spin period estimates for 1929 asteroids. We fit second-order Fourier models at a large number of trial frequencies to the W4 data and analyzed the resulting periodograms. We initially excluded rotational frequencies exceeding 7.57 rotations per day (P < 3.17 hr), which are not sampled adequately by WISE, and periods that exceed twice the WISE observation interval, which is typically 36 hr. Three solutions accurately capture the vast majority of the rotational frequencies in our sample: the best-fit frequency and its mirrors around 3.78 and 7.57 rotations per day. By comparing our solutions to a high-quality control group of 752 asteroid spin periods, we found that one of our solutions is accurate (within 5%) in 88% of the cases. The best-fit, secondary, and tertiary solutions are accurate in 55%, 27%, and 6% of the cases, respectively. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16409v1
Estimating effective wind speed from Gemini Planet Imager's adaptive optics data using covariance maps by Daniel M. Levinstein et al. on Wednesday 30 November The Earth's turbulent atmosphere results in speckled and blurred images of astronomical objects when observed by ground based visible and near-infrared telescopes. Adaptive optics (AO) systems are employed to reduce these atmospheric effects by using wavefront sensors (WFS) and deformable mirrors. Some AO systems are not fast enough to correct for strong, fast, high turbulence wind layers leading to the wind butterfly effect, or wind-driven halo, reducing contrast capabilities in coronagraphic images. Estimating the effective wind speed of the atmosphere allows us to calculate the atmospheric coherence time. This is not only an important parameter to understand for site characterization but could be used to help remove the wind butterfly in post processing. Here we present a method for estimating the atmospheric effective wind speed from spatio-temporal covariance maps generated from pseudo open-loop (POL) WFS data. POL WFS data is used as it aims to reconstruct the full wavefront information when operating in closed-loop. The covariance maps show how different atmospheric turbulent layers traverse the telescope. Our method successfully recovered the effective wind speed from simulated WFS data generated with the soapy python library. The simulated atmospheric turbulence profiles consist of two turbulent layers of ranging strengths and velocities. The method has also been applied to Gemini Planet Imager (GPI) AO WFS data. This gives insight into how the effective wind speed can affect the wind-driven halo seen in the AO image point spread function. In this paper, we will present results from simulated and GPI WFS data. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16441v1
Bar formation in submaximal disks: a challenge for Λ CDM by Tahere Kashfi et al. on Wednesday 30 November Bar formation in cosmological simulations of galaxy formation remains challenging. It was previously shown that the fraction of barred galaxies at low stellar masses ($M_*
Accelerating BAO Scale Fitting Using Taylor Series by Matthew Hansen et al. on Wednesday 30 November The Universe is currently undergoing accelerated expansion driven by dark energy. Dark energy's essential nature remains mysterious: one means of revealing it is by measuring the Universe's size at different redshifts. This may be done using the Baryon Acoustic Oscillation (BAO) feature, a standard ruler in the galaxy 2-Point Correlation Function (2PCF). In order to measure the distance scale, one dilates and contracts a template for the 2PCF in a fiducial cosmology, using a scaling factor $alpha$. The standard method for finding the best-fit $alpha$ is to compute the likelihood over a grid of roughly 100 values of it. This approach is slow; in this work, we propose a significantly faster way. Our method writes the 2PCF as a polynomial in $alpha$ by Taylor-expanding it about $alpha = 1$, exploiting that we know the fiducial cosmology sufficiently well that $alpha$ is within a few percent of unity. The likelihood resulting from this expansion may then be analytically solved for the best-fit $alpha$. Our method is 48-85$times$ faster than a directly comparable approach in which we numerically minimize $alpha$, and $sim$$12,000 times$ faster than the standard iterative method. Our work will be highly enabling for upcoming large-scale structure redshift surveys such as that by Dark Energy Spectroscopic Instrument (DESI). arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2112.06438v2
New constraints on the dark matter density profiles of dwarf galaxies from proper motions of globular cluster streams by Khyati Malhan et al. on Wednesday 30 November The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter. Recently, in Malhan et al. (2021), we employed N-body simulations to show that the cuspy cold dark matter (CDM) subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripping within their parent halos prior to accretion onto the Milky Way. We previously found that clusters that are accreted within cuspy subhalos produce streams with larger physical widths and higher dispersions in line-of-sight velocity and angular momentum than streams that are accreted within cored subhalos. Here, we use the same suite of simulations to demonstrate that the dispersion in the tangential velocities of streams ($sigma_{v_mathrm{Tan}}$) is also sensitive to the central DM density profiles of their parent dwarfs and GCs that were accreted from: cuspy subhalos produce streams with larger $sigma_{v_mathrm{Tan}}$ than those accreted inside cored subhalos. Using Gaia EDR3 observations of multiple GC streams we compare their $sigma_{v_mathrm{Tan}}$ values with simulations. The measured $sigma_{v_mathrm{Tan}}$ values are consistent with both an ``in situ'' origin and with accretion inside cored subhalos of $Msim 10^{8-9}M_{odot}$ (or very low-mass cuspy subhalos of mass $sim 10^8M_{odot}$). Despite the large current uncertainties in $sigma_{v_mathrm{Tan}}$, we find a low probability that any of the progenitor GCs were accreted from cuspy subhalos of $M_{rm subhalo}buildrel > over sim$ $10^9 M_{odot}$. The uncertainties on Gaia tangential velocity measurements are expected to decrease in future and will allow for stronger constraints on subhalo DM density profiles. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2201.03571v2
First release of PLATO consortium stellar limb-darkening coefficients by Giuseppe Morello et al. on Wednesday 30 November We release the first grid of stellar limb-darkening coefficients (LDCs) and intensity profiles (IPs) computed by the consortium of the PLAnetary Transits and Oscillations of stars (PLATO), the next medium-class (M3) mission under development by the European Space Agency (ESA) to be launched in 2026. We have performed spectral synthesis with texttt{TurboSpectrum} on a grid of texttt{MARCS} model atmospheres. Finally, we adopted texttt{ExoTETHyS} to convolve the high-resolution spectra ($R=2times10^5$) with the state-of-the-art response functions for all the PLATO cameras, and computed the LDCs that best approximate the convolved IPs. In addition to the PLATO products, we provide new LDCs and IPs for the Kepler mission, based on the same grid of stellar atmospheric models and calculation procedures. The data can be downloaded from the following link: url{https://doi.org/10.5281/zenodo.7339706}. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16510v1
Diffuse Emission of Galactic High-Energy Neutrinos from a Global Fit of Cosmic Rays by Georg Schwefer et al. on Wednesday 30 November In the standard picture of galactic cosmic rays, a diffuse flux of high-energy gamma-rays and neutrinos is produced from inelastic collisions of cosmic ray nuclei with the interstellar gas. The neutrino flux is a guaranteed signal for high-energy neutrino observatories such as IceCube, but has not been found yet. Experimental searches for this flux constitute an important test of the standard picture of galactic cosmic rays. Both the observation and non-observation would allow important implications for the physics of cosmic ray acceleration and transport. We present DINECRAFT, a new model of galactic diffuse high-energy gamma-rays and neutrinos, fitted to recent cosmic ray data from AMS-02, DAMPE, IceTop as well as KASCADE. We quantify the uncertainties for the predicted emission from the cosmic ray model, but also from the choice of source distribution, gas maps and cross-sections. We consider the possibility of a contribution from unresolved sources. Our model predictions exhibit significant deviations from older models. Our fiducial model is available at https://doi.org/10.5281/zenodo.7373010 . arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15607v1
The diversity of rotation curves of simulated galaxies with cusps and cores by Finn A. Roper et al. on Wednesday 30 November We use $Lambda$CDM cosmological hydrodynamical simulations to explore the kinematics of gaseous discs in late-type dwarf galaxies. We create high-resolution 21-cm 'observations' of simulated dwarfs produced in two variations of the EAGLE galaxy formation model: one where supernova-driven gas flows redistribute dark matter and form constant-density central 'cores', and another where the central 'cusps' survive intact. We 'observe' each galaxy along multiple sight lines and derive a rotation curve for each observation using a conventional tilted-ring approach to model the gas kinematics. We find that the modelling process introduces systematic discrepancies between the recovered rotation curve and the actual circular velocity curve driven primarily by (i) non-circular gas orbits within the discs; (ii) the finite thickness of gaseous discs, which leads to overlap of different radii in projection; and (iii) departures from dynamical equilibrium. Dwarfs with dark matter cusps often appear to have a core, whilst the inverse error is less common. These effects naturally reproduce an observed trend which other models struggle to explain: late-type dwarfs with more steeply-rising rotation curves appear to be dark matter-dominated in the inner regions, whereas the opposite seems to hold in galaxies with core-like rotation curves. We conclude that if similar effects affect the rotation curves of observed dwarfs, a late-type dwarf population in which all galaxies have sizeable dark matter cores is most likely incompatible with current measurements. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2203.16652v3
Non-ideal magnetohydrodynamics of self-gravitating filaments by Nicol Gutiérrez-Vera et al. on Wednesday 30 November Filaments have been studied in detail through observations and simulations. A range of numerical works have separately investigated how chemistry and diffusion effects, as well as magnetic fields and their structure impact the gas dynamics of the filament. However, non-ideal effects have hardly been explored thus far. We investigate how non-ideal magnetohydrodynamic (MHD) effects, combined with a simplified chemical model affect the evolution and accretion of a star-forming filament. We modeled an accreting self-gravitating turbulent filament using lemongrab, a one-dimensional (1D) non-ideal MHD code that includes chemistry. We explore the influence of non-ideal MHD, the orientation and strength of the magnetic field, and the cosmic ray ionization rate, on the evolution of the filament, with particular focus on the width and accretion rate. We find that the filament width and the accretion rate are determined by the magnetic field properties, including the initial strength, the coupling with the gas controlled by the cosmic ray ionization rate, and the orientation of the magnetic field with respect to the accretion flow direction. Increasing the cosmic-ray ionization rate leads to a behavior closer to that of ideal MHD, reducing the magnetic pressure support and, hence, damping the accretion efficiency with a consequent broadening of the filament width. For the same reason, we obtained a narrower width and a larger accretion rate when we reduced the initial magnetic field strength. Overall, while these factors affect the final results by approximately a factor of~2, removing the non-ideal MHD effects results in a much greater variation (up to a factor of~7). The inclusion of non-ideal MHD effects and the cosmic-ray ionization is crucial for the study of self-gravitating filaments and in determining critical observable quantities, such as the filament width and accretion rate. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15615v1
Identifying the physical origin of gamma-ray bursts with supervised machine learning by Jia-Wei Luo et al. on Wednesday 30 November The empirical classification of gamma-ray bursts (GRBs) into long and short GRBs based on their durations is already firmly established. This empirical classification is generally linked to the physical classification of GRBs originating from compact binary mergers and GRBs originating from massive star collapses, or Type I and II GRBs, with the majority of short GRBs belonging to Type I and the majority of long GRBs belonging to Type II. However, there is a significant overlap in the duration distributions of long and short GRBs. Furthermore, some intermingled GRBs, i.e., short-duration Type II and long-duration Type I GRBs, have been reported. A multi-wavelength, multi-parameter classification scheme of GRBs is evidently needed. In this paper, we seek to build such a classification scheme with supervised machine learning methods, chiefly XGBoost. We utilize the GRB Big Table and Greiner's GRB catalog and divide the input features into three subgroups: prompt emission, afterglow, and host galaxy. We find that the prompt emission subgroup performs the best in distinguishing between Type I and II GRBs. We also find the most important distinguishing feature in prompt emission to be $T_{90}$, hardness ratio, and fluence. After building the machine learning model, we apply it to the currently unclassified GRBs to predict their probabilities of being either GRB class, and we assign the most probable class of each GRB to be its possible physical class. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16451v1
Parametrized Asymmetric Neutral hydrogen Disk Integrated Spectrum Characterization PANDISC I: Introduction to A Physically Motivated H I Model by Bo Peng et al. on Wednesday 30 November Modelling the integrated H I spectra of galaxies has been a difficult task due to their diverse shapes, but more dynamical information is waiting to be explored in Hi line profiles. Based on simple assumptions, we construct a physically motivated model for the integrated Hi spectra: Parametrized Asymmetric Neutral hydrogen Disk Integrated Spectrum Characterization (PANDISC). The model shows great flexibility in reproducing the diverse Hi profiles. We use Monte-Carlo Markov Chain (MCMC) for fitting the model to global H I profiles, producing statistically robust quantitative results. Comparing with several samples of H I data available in the literature , we find the model-fitted results agree with catalogued velocity widths (e.g., W50) down to the lowest S/N. The model is also shown to be useful for applications like the baryonic Tully-Fisher relation (BTFR) and profile-based sample control. By comparing v_r to v_flat , we uncover how the H I width is affected by the structure of the rotation curve, following a trend consistent with the difference in the BTFR slope. We also select a sample of spectra with broad wing-like features suggestive of a population of galaxies with unusual gas dynamics. The PANDISC model bears both promise and limitations for potential use beyond H I lines. Further application on the whole ALFALFA sample will enable us to perform large scale ensemble studies of the H I properties and dynamics in nearby galaxies. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16455v1
MOJAVE XX Decade-long linear polarization variability in AGN jets at parsec scales by D. I. Zobnina et al. on Wednesday 30 November We studied the variability properties of the linear polarization of active galactic nucleus (AGN) jets on parsec-scales using stacking. Our sample is drawn from the MOJAVE program and consists of 436 AGNs manifesting core-jet morphology and having at least five VLBA observing epochs at 15 GHz from January 1996 through August 2019, with some additional archival VLBA data reduced by us. We employed a stacking procedure and constructed maps of (i) standard deviation of fractional polarization and electric vector position angle (EVPA) over epochs as the measure of variability, (ii) median polarization degree to quantify typical values in time. The distributions of these values were analysed along and across the jet. We found that EVPA variability in the core is typically higher than in the jet, most likely due to changes in opacity and component blending in the core region. The EVPA becomes more stable down the outflow. Most of the sources showing this trend have a time coverage of more than 12 years and at least 15 epochs. The possible cause could be the increase of stability in the magnetic field direction. The majority of AGNs exhibit insignificant trends of the relative fractional polarization variability along the ridgeline or across the jet width. There are no significant optical-class-dependent or spectral-class-dependent relations in the properties of EVPA and relative fractional polarization variability. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15624v1
Accurate modelling of extragalactic microlensing by compact objects by VÃctor Boscá et al. on Wednesday 30 November Microlensing of extragalactic sources, in particular the probability of significant amplifications, is a potentially powerful probe of the abundance of compact objects outside the halo of the Milky Way. Accurate experimental constraints require an equally accurate theoretical model for the amplification statistics produced by such a population. In this article, we argue that the simplest (strongest-lens) model does not meet this demanding requirement. We thus propose an elaborate practical modelling scheme for extragalactic microlensing. We derive from first principles an expression for the amplification probability that consistently allows for: (i) the coupling between microlenses; (ii) realistic perturbations from the cosmic large-scale structure; (iii) extended-source corrections. An important conclusion is that the external shear applied on the dominant microlens, both by the other lenses and by the large-scale structure, is practically negligible. Yet, the predictions of our approach can still differ by a factor of a few with respect to existing models of the literature. Updated constraints on the abundance of compact objects accounting for such discrepancies may be required. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2205.00991v2
Separating the blue cloud and the red sequence using Otsu's method for image segmentation by Biswajit Pandey. on Wednesday 30 November The observed colour bimodality allows a classification of the galaxies into two distinct classes: the `blue cloud' and the `red sequence'. Such classification is often carried out using empirical cuts in colour and other galaxy properties that lack solid mathematical justifications. We propose a method for separating the galaxies in the `blue cloud' and the `red sequence' using Otsu's thresholding technique for image segmentation. We show that this technique provides a robust and parameter-free method for the classification of the red and blue galaxies based on the minimization of the inter-class variance and maximization of the intra-class variance. We also apply an iterative triclass thresholding technique based on Otsu's method to improve the classification. The same method can also be applied to classify the galaxies based on their physical properties, such as star formation rate, stellar mass function, bulge-to-disk mass ratio and age, all of which have bimodal distributions. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15642v1
Constraints on dark matter annihilation and decay from the large-scale structure of the nearby universe by Deaglan J. Bartlett et al. on Wednesday 30 November Decaying or annihilating dark matter particles could be detected through gamma-ray emission from the species they decay or annihilate into. This is usually done by modelling the flux from specific dark matter-rich objects such as the Milky Way halo, Local Group dwarfs, and nearby groups. However, these objects are expected to have significant emission from baryonic processes as well, and the analyses discard gamma-ray data over most of the sky. Here we construct full-sky templates for gamma-ray flux from the large-scale structure within $sim$200 Mpc by means of a suite of constrained $N$-body simulations (CSiBORG) produced using the Bayesian Origin Reconstruction from Galaxies algorithm. Marginalising over uncertainties in this reconstruction, small-scale structure, and parameters describing astrophysical contributions to the observed gamma-ray sky, we compare to observations from the Fermi Large Area Telescope to constrain dark matter annihilation cross sections and decay rates through a Markov Chain Monte Carlo analysis. We rule out the thermal relic cross section for $s$-wave annihilation for all $m_chi lesssim 7 {rm , GeV}/c^2$ at 95% confidence if the annihilation produces gluons or quarks less massive than the bottom quark. We infer a contribution to the gamma-ray sky with the same spatial distribution as dark matter decay at $3.3sigma$. Although this could be due to dark matter decay via these channels with a decay rate $Gamma approx 6 times 10^{-28} {rm , s^{-1}}$, we find that a power-law spectrum of index $p=-2.75^{+0.71}_{-0.46}$, likely of baryonic origin, is preferred by the data. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2205.12916v2
Calibrating cosmological simulations with implicit likelihood inference using galaxy growth observables by Yongseok Jo et al. on Wednesday 30 November In a novel approach employing implicit likelihood inference (ILI), also known as likelihood-free inference, we calibrate the parameters of cosmological hydrodynamic simulations against observations, which has previously been unfeasible due to the high computational cost of these simulations. For computational efficiency, we train neural networks as emulators on ~1000 cosmological simulations from the CAMELS project to estimate simulated observables, taking as input the cosmological and astrophysical parameters, and use these emulators as surrogates to the cosmological simulations. Using the cosmic star formation rate density (SFRD) and, separately, stellar mass functions (SMFs) at different redshifts, we perform ILI on selected cosmological and astrophysical parameters (Omega_m, sigma_8, stellar wind feedback, and kinetic black hole feedback) and obtain full 6-dimensional posterior distributions. In the performance test, the ILI from the emulated SFRD (SMFs) can recover the target observables with a relative error of 0.17% (0.4%). We find that degeneracies exist between the parameters inferred from the emulated SFRD, confirmed with new full cosmological simulations. We also find that the SMFs can break the degeneracy in the SFRD, which indicates that the SMFs provide complementary constraints for the parameters. Further, we find that the parameter combination inferred from an observationally-inferred SFRD reproduces the target observed SFRD very well, whereas, in the case of the SMFs, the inferred and observed SMFs show significant discrepancies that indicate potential limitations of the current galaxy formation modeling and calibration framework, and/or systematic differences and inconsistencies between observations of the stellar mass function. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16461v1
Statistical Study of the Star Formation Efficiency in Bars: Is Star Formation Suppressed in Gas-Rich Bars? by Fumiya Maeda et al. on Wednesday 30 November The dependence of star formation efficiency (SFE) on galactic structures, especially whether the SFE in the bar region is lower than those in the other regions, has recently been debated. We report the SFEs of 18 nearby gas-rich massive star-forming barred galaxies with a large apparent bar major axis ($geqq 75^{primeprime}$). We statistically measure the SFE by distinguishing the center, bar-end, and bar regions for the first time. The molecular gas surface density is derived from archival CO(1-0) and/or CO(2-1) data by assuming a constant CO-to-H$_2$ conversion factor ($alpha_{rm CO}$), and the star formation rate surface density is derived from a linear combination of far-ultraviolet and mid-infrared intensities. The angular resolution is $15^{primeprime}$, which corresponds to $0.3 - 1.8~rm kpc$. We find that the ratio of the SFE in the bar to that in the disk was systematically lower than unity (typically $0.6-0.8$), which means that the star formation in the bar is systematically suppressed. Our results are inconsistent with similar recent statistical studies that reported that SFE tends to be independent of galactic structures. This inconsistency can be attributed to the differences in the definition of the bar region, spatial resolution, $alpha_{rm CO}$, and sample galaxies. Furthermore, we find a negative correlation between SFE and velocity width of the CO spectrum, which is consistent with the idea that the large dynamical effects, such as strong shocks, large shear, and fast cloud-cloud collisions caused by the noncircular motion of the bar, result in a low SFE. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15681v1