Abstract visibility
Abstract text size

Papers for Wednesday, Jul 06 2022

Papers with local authors

Lachlan Lancaster, Sarah Pearson, Benjamin F. Williams, Kathryn V. Johnston, Tjitske K. Starkenburg, Erin Kado-Fong, Anil C. Seth, Eric F. Bell

16 pages, 6 figures, submitted to AJ. Comments welcome

0 votes
Paper 40 — arXiv:2207.02208
0 votes
Paper 40 — arXiv:2207.02208

Studies of resolved stellar populations in the Milky Way and nearby galaxies reveal an amazingly detailed and clear picture of galaxy evolution. Within the Local Group, the ability to probe the stellar populations of small and large galaxies opens up the possibility of exploring key questions such as the nature of dark matter, the detailed formation history of different galaxy components, and the role of accretion in galactic formation. Upcoming wide-field surveys promise to extend this ability to all galaxies within 10~Mpc, drastically increasing our capability to decipher galaxy evolution and enabling statistical studies of galaxies' stellar populations. To facilitate the optimum use of these upcoming capabilities we develop a simple formalism to predict the density of resolved stars for an observation of a stellar population at fixed surface brightness and population parameters. We provide an interface to calculate all quantities of interest to this formalism via a public release of the code: \texttt{walter}. This code enables calculation of (i) the expected number density of detected stars, (ii) the exposure time needed to reach certain population features, such as the horizontal branch, and (iii) an estimate of the crowding limit, among other features. These calculations will be very useful for planning surveys with NASA's upcoming Nancy Grace Roman Space Telescope (Roman, formerly WFIRST), which we use for example calculations throughout this work.

All other papers

Michael S. Calzadilla, Michael McDonald, Megan Donahue, Brian R. McNamara, Kevin Fogarty, Massimo Gaspari, Myriam Gitti, Helen R. Russell, Grant R. Tremblay, G. Mark Voit, Francesco Ubertosi

26 pages, 5 figures, 5 tables. Submitted to ApJ

We present new, deep, narrow- and broad-band Hubble Space Telescope observations of seven of the most star-forming brightest cluster galaxies (BCGs). Continuum-subtracted [O II] maps reveal the detailed, complex structure of warm ($T \sim 10^4$ K) ionized gas filaments in these BCGs, allowing us to measure spatially-resolved star formation rates (SFRs) of ~60-600 Msun/yr. We compare the SFRs in these systems and others from the literature to their intracluster medium (ICM) cooling rates (dM/dt), measured from archival Chandra X-ray data, finding a best-fit relation of log(SFR) = (1.67+/-0.17) log(dM/dt) + (-3.25+/-0.38) with an intrinsic scatter of 0.39+/-0.09 dex. This steeper-than-unity slope implies an increasingly efficient conversion of hot ($T \sim 10^7$ K) gas into young stars with increasing dM/dt, or conversely a gradual decrease in the effectiveness of AGN feedback in the strongest cool cores. We also seek to understand the physical extent of these multiphase filaments that we observe in cluster cores. We show, for the first time, that the average extent of the multiphase gas is always smaller than the radii at which the cooling time reaches 1 Gyr, the tcool/tff profile flattens, and that X-ray cavities are observed. This implies a close connection between the multiphase filaments, the thermodynamics of the cooling core, and the dynamics of X-ray bubbles. Interestingly, we find a one-to-one correlation between the average extent of cool multiphase filaments and the radius at which the cooling time reaches 0.5 Gyr, which may be indicative of a universal condensation timescale in cluster cores.

Matthew Alessi, Ralph E. Pudritz

30 pages (including 4 page Appendix), 10 figures in main text ( 3 in Appendix), 1 movie (clickable in caption of Figure 7). Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Recent surveys show that protoplanetary disks have lower levels of turbulence than expected based on their observed accretion rates. A viable solution to this is that magnetized disk winds dominate angular momentum transport. This has several important implications for planet formation processes. We compute the physical and chemical evolution of disks and the formation and migration of planets under the combined effects of angular momentum transport by turbulent viscosity and disk winds. We take into account the critical role of planet traps to limit Type I migration in all of these models and compute thousands of planet evolution tracks for single planets drawn from a distribution of initial disk properties and turbulence strengths. We do not consider multi-planet models nor include N-body planet-planet interactions. Within this physical framework we find that populations with a constant value disk turbulence and winds strength produce mass-semimajor axis distributions in the M-a diagram with insufficient scatter to compare reasonably with observations. However, populations produced as a consequence of sampling disks with a distribution of the relative strengths of disk turbulence and winds fit much better. Such models give rise to a substantial super Earth population at orbital radii 0.03-2 AU, as well as a clear separation between the produced hot Jupiter and warm Jupiter populations. Additionally, this model results in a good comparison with the exoplanetary mass-radius distribution in the M-R diagram after post-disk atmospheric photoevaporation is accounted for.

Silvan Fischbacher (1), Tomasz Kacprzak (1 and 2), Jonathan Blazek (3), Alexandre Refregier (1) ((1) ETH Zurich, (2) Swiss Data Science Center, (3) Northeastern University)

36 pages including 15 figures and 9 tables

Intrinsic alignment (IA) modelling and photometric redshift estimation are two of the main sources of systematic uncertainty in weak lensing surveys. We investigate the impact of redshift errors and their interplay with different IA models. Generally, errors on the mean $\delta_z$ and on the width $\sigma_z$ of the redshift bins can both lead to biases in cosmological constraints. We find that such biases can, however, only be partially resolved by marginalizing over $\delta_z$ and $\sigma_z$. For Stage-III surveys, $\delta_z$ and $\sigma_z$ cannot be well constrained due to limited statistics. The resulting biases are thus sensitive to prior volume effects. For Stage-IV surveys, we observe that marginalizing over the redshift parameters has an impact and reduces the bias. We derive requirements on the uncertainty of $\sigma_z$ and $\delta_z$ for both Stage-III and Stage-IV surveys. We assume that the redshift systematic errors on $S_8$ should be less than half of the statistical errors, and the median bias should be smaller than $0.25\sigma$. We find that the uncertainty on $\delta_z$ has to be $\lesssim0.025$ for the NLA IA model with a Stage-III survey. For $\sigma_z$, the requirement is met even for large uncertainties $\leq0.3$. For the TATT IA model, the uncertainty on $\delta_z$ has to be $\lesssim0.02$ and the uncertainty on $\sigma_z$ has to be $\lesssim0.2$. For Stage-IV surveys, the uncertainty on $\delta_z$ has to be $\lesssim0.005$ and the uncertainty on $\sigma_z$ should be $\lesssim0.1$, with no significant dependence on the IA model. This required high precision will be a challenge for the redshift calibration of these future surveys. Finally, we investigate whether the interplay between redshift systematics and IA modelling can explain the $S_8$-tension between cosmic shear results and CMB measurements. We find that this is unlikely to explain the current $S_8$-tension.

J. Casares, T. Munõz-Darias, M.A.P. Torres, D. Mata Sanchez, C.T. Britt, M. Armas Padilla, A. Alvarez-Hernandez, V.A. Cuneo, J.I. Gonzalez Hernandez, F. Jimenez-Ibarra, P.G. Jonker, G. Panizo-Espinar, J. Sanchez-Sierras, I.V. Yanes-Rizo

Accepted for publication in MNRAS, 21 pages, 15 figures, 6 Tables

We present a new method to derive binary inclinations in quiescent black hole (BH) X-ray transients (XRTs), based on the depth of the trough (T) from double-peaked Ha emission profiles arising in accretion discs. We find that the inclination angle (i) is linearly correlated with T in phase-averaged spectra with sufficient orbital coverage (>~50 per cent) and spectral resolution, following i (deg)=93.5 x T +23.7. The correlation is caused by a combination of line opacity and local broadening, where a leading (excess broadening) component scales with the de-projected velocity of the outer disc. Interestingly, such scaling allows to estimate the fundamental ratio M1/Porb by simply resolving the intrinsic width of the double-peak profile. We apply the T-i correlation to derive binary inclinations for GRO J0422+32 and Swift J1357-0933, two BH XRTs where strong flickering activity has hindered determining their values through ellipsoidal fits to photometric light curves. Remarkably, the inclination derived for GRO J0422+32 (i=55.6+-4.1) implies a BH mass of 2.7+0.7-0.5 Msun thus placing it within the gap that separates BHs from neutron stars. This result proves that low-mass BHs exist in nature and strongly suggests that the so-called "mass gap" is mainly produced by low number statistics and possibly observational biases. On the other hand, we find that Swift J1357-0933 contains a 10.9+1.7-1.6 Msun BH seen nearly edge on (i=87.4+2.6-5.6 deg). Such extreme inclination, however, should be treated with caution since it relies on extrapolating the T-i correlation beyond i>~75 deg, where it has not yet been tested.

Guido Roberts-Borsani, Tommaso Treu, Charlotte Mason, Richard S. Ellis, Nicolas Laporte, Thomas Schmidt, Maruša Bradač, Adriano Fontana, Takahiro Morishita, Paola Santini

14 pages, 4 figures, 2 tables. Submitted to ApJ

The 100% detection rate of Ly$\alpha$ in a sample of four luminous $z\sim8$ galaxies with red Spitzer/IRAC colors suggested objects with unusual ionizing capabilities and the presence of early (re)ionized bubbles in a neutral era. Whether such bubbles are a reflection of enhanced ionizing properties (nature) or an overdense environment (nurture), however, remains an open question. Here we aim to distinguish between these hypotheses with a spectroscopic search for Ly$\alpha$ in five fainter galaxies from the CANDELS-GOODS fields with Keck/MOSFIRE, selected on a similar IRAC excess and UV magnitudes that reflect reduced clustering effects. We detect $>4\sigma$ line emission in only two targets, placing them at redshifts of $z_{\rm Ly\alpha}=7.1081$ and $z_{\rm Ly\alpha}=7.9622$, with rest-frame EWs of 16-17 A that are $\sim1.5\times$ weaker than those of their bright counterparts. Comparing line detections we find a reduced rate of Ly$\alpha$ emission of $0.40^{+0.30}_{-0.25}$ in the faint sample compared to $1.00^{+0.00}_{-0.44}$ for the bright sample. The lower rate of the faint sample is well-matched by the predicted number of detections derived from simulations of a mostly neutral IGM and an intrinsic EW$_{0,Ly\alpha}$ distribution from $z\sim6$ galaxies. However, even with an extreme EW model, it is challenging to match the detection rate of the luminous objects. SED-fitting of the faint galaxies indicates generally young, metal- and dust-poor systems, albeit with reduced star-forming properties compared to their luminous counterparts. We argue the enhanced detection rate in UV bright galaxies is likely a byproduct of both extreme ionizing properties as well as environment, where nature and nurture both play an important role. Ultimately, UV emission lines and direct measurements of clustering with JWST are required to resolve the physical nature of this puzzling population.

As established in previous papers of this series, observables in highly distorted and magnified multiple images caused by the strong gravitational lensing effect can be used to constrain the distorting properties of the gravitational lens at the image positions. If the background source is extended and contains substructure, like star forming regions, which is resolved in multiple images, all substructure that can be matched across a minimum of three multiple images can be used to infer the local distorting properties of the lens. In this work, we replace the manual feature selection by an automated feature extraction based on SExtractor for Python and show its superior performance. Despite its aimed development to improve our lens reconstruction, it can be employed in any other approach, as well. Valuable insights on the definition of an `image position' in the presence of noise are gained from our calibration tests. Applying it to observations of a five-image configuration in galaxy cluster CL0024 and the triple-image configuration containing Hamilton's object, we determine local lens properties for multiple wavebands separately. Within current confidence bounds, all of them are consistent with each other, corroborating the wavelength-independence of strong lensing and offering a tool to detect deviations caused by micro-lensing and dust in further examples.

Roberto Soria, Ruican Ma, Lian Tao, Shuang-Nan Zhang

8 pages, 300 kB, accepted by MNRAS on 2022 July 4

We revisit various sets of published results from X-ray and optical studies of the Galactic black hole (BH) candidate MAXI J0637-430, which went into outburst in 2019. Combining the previously reported values of peak outburst luminosity, best-fitting radii of inner and outer accretion disk, viewing angle, exponential decay timescale and peak-to-peak separation of the He II 4686 disk emission line, we improve the constraints on the system parameters. We estimate a heliocentric distance d = (8.7 +/- 2.3) kpc, a projected Galactocentric distance R = (13.2 +/- 1.8) kpc and a height |z| = (3.1 +/- 0.8) kpc from the Galactic plane. It is the currently known Milky Way BH candidate located farthest from the Galactic Centre. We infer a BH mass M_1 = (5.1 +/- 1.6) M_{sun}, a spin parameter a* <~ 0.25, a donor star mass M_2 = (0.25 +/- 0.07) M_{sun}, a peak Eddington ratio lambda = 0.17 +/- 0.11 and a binary period P_{orb} = 2.2^{+0.8}_{-0.6} hr. This is the shortest period measured or estimated so far for any Galactic BH X-ray binary. If the donor star is a main-sequence dwarf, such a period corresponds to the evolutionary stage where orbital shrinking is driven by gravitational radiation and the star has regained contact with its Roche lobe (low end of the period gap). The three Galactic BHs with the shortest period (<~3 hr) are also those with the highest vertical distance from the Galactic plane (>~2 kpc). This is probably because binaries with higher binding energies can survive faster natal kicks.

Paolo Suin, Steven N. Shore, Václav Pavlík

9 pages, 2 appendices, 4 tables, 6 figures, accepted for publication to A&A

Context: Star clusters form within giant molecular clouds that are strongly altered by the feedback action of the massive stars, but the cluster still remains embedded in a dense, highly turbulent medium and interactions with ambient structures may modify its dynamical evolution from that expected if it were isolated. Aims: We aim to study coupling mechanisms between the dynamical evolution of the cluster, accelerated by the mass segregation process, with harassment effects caused by the gaseous environment. Methods: We simulated the cluster dynamical evolution combining $N$-body and hydrodynamic codes within the Astronomical Multipurpose Software Environment (AMUSE). Conclusions: Tidal harassment produces a sparser configuration more rapidly than the isolated reference simulations. The evolution of the asymptotic power-law density distribution exponent also shows substantially different behaviour in the two cases. The background is more effective on clusters in advanced stages of dynamical development.

Shreya Vardhan, Sašo Grozdanov, Samuel Leutheusser, Hong Liu

6+7 pages, 3 figures

We present a formulation of magnetohydrodynamics which can be used to describe the evolution of strong magnetic fields in neutron star interiors. Our approach is based on viewing magnetohydrodynamics as a theory with a one-form global symmetry and developing an effective field theory for the hydrodynamic modes associated with this symmetry. In the regime where the local velocity and temperature variations can be neglected, we derive the most general constitutive relation consistent with symmetry constraints for the electric field in the presence of a strong magnetic field. This constitutive relation not only reproduces the phenomena of Ohmic decay, ambipolar diffusion, and Hall drift derived in a phenomenological model by Goldreich and Reisenegger, but also reveals new terms in the evolution of the magnetic field which cannot easily be seen from such microscopic models. This formulation gives predictions for novel diffusion behaviors of small perturbations around a constant background magnetic field, and for the two-point correlation functions among various components of the electric and magnetic fields.

Jubee Sohn, Margaret J. Geller, Mark Vogelsberger, Josh Borrow

15 pages, 13 figures, 1 table, submitted to the ApJ

We explore the redshift evolution of dynamical scaling relations between massive clusters and their brightest cluster galaxies (BCGs) at $z < 2$ based on the IllustrisTNG-300 simulation. We select 270 massive clusters with $M_{200} > 10^{14}$ M$_{\odot}$ at $z = 0$ and trace their progenitors based on merger trees. From 67 redshift snapshots covering $z < 2$, we use the subhalo velocity dispersion to compute the cluster velocity dispersion ($\sigma_{cl}$). We also calculate the stellar velocity dispersion of the BCGs ($\sigma_{BCG}$). Both $\sigma_{cl}$ and $\sigma_{BCG}$ increase as universe ages. The BCG velocity dispersion grows more slowly than the cluster velocity dispersion. At $z > 1$, $\sigma_{BCG}$ is comparable with $\sigma_{cl}$ offering an interesting observational test. We also derive the $(\sigma_{BCG}/ \sigma_{cl}) - \sigma_{cl}$ scaling relation; the slope of the relation steepens at high redshift. The simulated redshift evolution of $\sigma_{cl}$ and $\sigma_{BCG}$ generally agrees with an observed cluster sample for $z < 0.3$ but with large scatter. Future large spectroscopic surveys will test the implications of the simulations for the mass evolution of both clusters and their BCGs.

Ryota Shimada, Hideyuki Hotta, Takaaki Yokoyama

17 pages, 20 figures, Accepted for publication in ApJ

Solar magnetic fields comprise an 11-year activity cycle, represented by the number of sunspots. The maintenance of such a solar magnetic field can be attributed to fluid motion in the convection zone, i.e. a dynamo. This study conducts the mean-field analyses of the global solar dynamo simulation presented by Hotta et al. (2016). Although the study succeeds in producing coherent large-scale magnetic fields at high Reynolds numbers, the detailed physics of the maintenance of this field have not been fully understood. This study extracts the alpha-tensor and the turbulent magnetic diffusivity tensor through mean-field analyses. The turbulent magnetic diffusivity exhibits a significant decrease towards high Reynolds numbers. The decrease in the turbulent magnetic diffusivity suppresses the energy conversion of large-scale field to small-scale field. This implies that the decrease in the turbulent magnetic diffusivity contributes to the maintenance of a large-scale magnetic field at high Reynolds numbers. A significant downward turbulent pumping is observed; it is enhanced in the weak phase of the large-scale field. This study proposes a cyclic reversal process of a large-scale field which is dominantly driven by the alpha-effect and is possibly triggered by downward pumping.

Lara Alegre, Jose Sabater, Philip Best, Rafaël I.J. Mostert, Wendy L. Williams, Gülay Gürkan, Martin J. Hardcastle, Rohit Kondapally, Tim W. Shimwell, Daniel J.B. Smith

Accepted for publication in MNRAS. 24 pages, 13 Figures. Submission includes a data table which is also available at this https URL

New-generation radio telescopes like LOFAR are conducting extensive sky surveys, detecting millions of sources. To maximise the scientific value of these surveys, radio source components must be properly associated into physical sources before being cross-matched with their optical/infrared counterparts. In this paper, we use machine learning to identify those radio sources for which either source association is required or statistical cross-matching to optical/infrared catalogues is unreliable. We train a binary classifier using manual annotations from the LOFAR Two-metre Sky Survey (LoTSS). We find that, compared to a classification model based on just the radio source parameters, the addition of features of the nearest-neighbour radio sources, the potential optical host galaxy, and the radio source composition in terms of Gaussian components, all improve model performance. Our best model, a gradient boosting classifier, achieves an accuracy of 95 per cent on a balanced dataset and 96 per cent on the whole (unbalanced) sample after optimising the classification threshold. Unsurprisingly, the classifier performs best on small, unresolved radio sources, reaching almost 99 per cent accuracy for sources smaller than 15 arcsec, but still achieves 70 per cent accuracy on resolved sources. It flags 68 per cent more sources than required as needing visual inspection, but this is still fewer than the manually-developed decision tree used in LoTSS, while also having a lower rate of wrongly accepted sources for statistical analysis. The results have an immediate practical application for cross-matching the next LoTSS data releases and can be generalised to other radio surveys.

Marco Simonte, Heinz Andernach, Marcus Brüggen, Dominik Schwarz, Isabella Prandoni, Anthony G. Willis

Accepted for publication in MNRAS. 23 pages, 16 figures, 3 tables

Giant radio galaxies (GRGs) are radio galaxies that have projected linear extents of more than 700 kpc or 1 Mpc, depending on definition. We have carried out a careful visual inspection in search of GRGs of the Bootes LOFAR Deep Field (BLDF) image at 150 MHz. We identified 74 GRGs with a projected size larger than 0.7 Mpc of which 38 are larger than 1 Mpc. The resulting GRG sky density is about 2.8 (1.43) GRGs per square degree for GRGs with linear size larger than 0.7 (1) Mpc. We studied their radio properties and the accretion state of the host galaxies using deep optical and infrared survey data and determined flux densities for these GRGs from available survey images at both 54 MHz and 1.4 GHz to obtain integrated radio spectral indices. We show the location of the GRGs in the P-D diagram. The accretion mode onto the central black holes of the GRG hosts is radiatively inefficient suggesting that the central engines are not undergoing massive accretion at the time of the emission. Interestingly, 14 out of 35 GRGs for which optical spectra are available show a moderate star formation rate. Based on the number density of optical galaxies taken from the DESI DR9 photometric redshift catalogue, we found no significant differences between the environments of GRGs and other radio galaxies, at least for redshift up to z = 0.7.

Giulia Valvano, Othon Cabo Winter, Rafael Sfair, Rai Machado Oliveira, Gabriel Borderes-Motta

12 pages, 13 figures

The first proposed Brazilian mission to deep space, the ASTER mission, has the triple asteroid system (153591) 2001 SN263 as a target. One of the mission's main goals is to analyze the physical and dynamical structures of the system to understand its origin and evolution. The present work aims to analyze how the asteroid's irregular shape interferes with the stability around the system. The results show that the irregular shape of the bodies plays an important role in the dynamics nearby the system. For instance, the perturbation due to the (153591) 2001 SN263 Alpha's shape affects the stability in the (153591) 2001 SN263 Gamma's vicinity. Similarly, the (153591) 2001 SN263 Beta's irregularity causes a significant instability in its nearby environment. As expected, the prograde case is the most unstable, while the retrograde scenario presents more stability. Additionally, we investigate how the solar radiation pressure perturbs particles of different sizes orbiting the triple system. We found that particles with a 10-50 cm radius could survive the radiation pressure for the retrograde case. Meanwhile, to resist solar radiation, the particles in prograde orbit must be larger than the particles in retrograde orbits, at least one order of magnitude.

Gianfranco Bino, Shantanu Basu, Masahiro N Machida, Aris Tritsis, Mahmoud Sharkawi, Kundan Kadam, Indrani Das

The canonical theory of star formation in a magnetized environment predicts the formation of hourglass-shaped magnetic fields during the prestellar collapse phase. In protostellar cores, recent observations reveal complex and strongly distorted magnetic fields in the inner regions that are sculpted by rotation and outflows. We conduct resistive, nonideal magnetohydrodynamic (MHD) simulations of a protostellar core and employ the radiative transfer code POLARIS to produce synthetic polarization segment maps. Comparison of our mock-polarization maps based on the toroidal-dominated magnetic field in the outflow zone with the observed polarization vectors of SiO lines in Orion Source I shows a reasonable agreement when the magnetic axis is tilted at an angle $\theta = 15^{\circ}$ with respect to the plane-of-sky and if the SiO lines have a net polarization parallel to the local magnetic field. Although the observed polarization is from SiO lines and our synthetic maps are due to polarized dust emission, a comparison is useful and allows us to resolve the ambiguity of whether the line polarization is parallel or perpendicular to the local magnetic field direction.

L. A. Rachmeler, J. Trujillo Bueno, D. E. McKenzie, R. Ishikawa, F. Auchere, K. Kobayashi, R. Kano, T. J. Okamoto, C. W. Bethge, D. Song, E. Alsina Ballester, L. Belluzzi, T. del Pino Aleman, M. Yoshida, T. Shimizu, A. Winebarger, A. R. Kobelski, G. D. Vigil, B. De Pontieu, Y. Suematsu, N. Narukage, M. Kubo, T. Sakao, H. Hara, J. Stepan, M. Carlsson, J. Leenaarts, A. Asensio Ramos

14 pages, 5 figures, to be published in the Astrophysical Journal (ApJ), initially submitted May 2022, revised submission July 2022

The CLASP2 (Chromospheric LAyer SpectroPolarimeter 2) sounding rocket mission was launched on 2019 April 11. CLASP2 measured the four Stokes parameters of the Mg II h & k spectral region around 2800 Angstroms along a 200 arcsecond slit at three locations on the solar disk, achieving the first spatially and spectrally resolved observations of the solar polarization in this near ultraviolet region. The focus of the work presented here is the center-to-limb variation of the linear polarization across these resonance lines, which is produced by the scattering of anisotropic radiation in the solar atmosphere. The linear polarization signals of the Mg II h & k lines are sensitive to the magnetic field from the low to the upper chromosphere through the Hanle and magneto-optical effects. We compare the observations to theoretical predictions from radiative transfer calculations in unmagnetized semi-empirical models, arguing that magnetic fields and horizontal inhomogeneities are needed to explain the observed polarization signals and spatial variations. This comparison is an important step in both validating and refining our understanding of the physical origin of these polarization signatures, and also in paving the way toward future space telescopes for probing the magnetic fields of the solar upper atmosphere via ultraviolet spectropolarimetry.

Raghunathan Srianand, Neeraj Gupta, Patrick Petitjean, Emmanuel Momjian, Sergei A. Balashev, Francoise Combes, Hsiao-Wen Chen, Jens-Kristian Krogager, Pasquier Noterdaeme, Hadi Rahmani, Andrew J. Baker, Kimberly L. Emig, Gyula I. G. Jozsa, Hans-Rainer Kloeckner, Kavilan Moodley

8 pages, 7 figures and accepted for publication in MNRAS

We report the emergence of a new HI 21-cm absorption at z_abs = 1.172635 in the damped Lyman-alpha absorber (DLA) towards the gamma-ray blazar PKS 2355-106 (z_em~1.639) using science verification observations (June 2020) from the MeerKAT Absorption Line Survey (MALS). Since 2006, this DLA is known to show a narrow HI 21-cm absorption at z_abs = 1.173019 coinciding with a distinct metal absorption line component. We do not detect significant HI 21-cm optical depth variations from this known HI component. A high resolution optical spectrum (August 2010) shows a distinct Mg I absorption at the redshift of the new HI 21-cm absorber. However, this component is not evident in the profiles of singly ionized species. We measure the metallicity ([Zn/H] = -(0.77\pm0.11) and [Si/H]= -(0.96\pm0.11)) and depletion ([Fe/Zn] = -(0.63\pm0.16)) for the full system. Using the apparent column density profiles of Si II, Fe II and Mg I we show that the depletion and the N(Mg I)/N(Si II) column density ratio systematically vary across the velocity range. The region with high depletion tends to have slightly larger N(Mg I)/N(Si II) ratio. The two HI 21-cm absorbers belong to this velocity range. The emergence of z_abs = 1.172635 can be understood if there is a large optical depth gradient over a length scale of ~0.35 pc. However, the gas producing the z_abs = 1.173019 component must be nearly uniform over the same scale. Systematic uncertainties introduced by the absorption line variability has to be accounted for in experiments measuring the variations of fundamental constants and cosmic acceleration even when the radio emission is apparently compact as in PKS 2355-106.

S. Miyake, T. Koi, Y. Muraki, Y. Matsubara, S. Masuda, P. Miranda, T. Naito, E. Ortiz, A. Oshima, T. Sakai, T. Sako, S. Shibata, H. Takamaru, M. Tokumaru, J. F. Valdes-Galicia

Paper presented in the 21st International Symposium on Very High Energy Cosmic Ray Intercations (ISVHE-CRI 2022) by online

In association with a large solar flare on November 7, 2004, the solar neutron detectors located at Mt. Chacaltaya (5,250m) in Bolivia and Mt. Sierra Negra (4,600m) in Mexico recorded very interesting events. In order to explain these events, we have performed a calculation solving the equation of motion of anti-protons inside the magnetosphere. Based on these results, the Mt. Chacaltaya event may be explained by the detection of solar neutrons, while the Mt. Sierra Negra event may be explained by the first detection of very high energy solar neutron decay protons (SNDPs) around 6 GeV.

Qiao Song, Xianyong Bai, Bo Chen, Xiuqing Hu, Yajie Chen, Zhenyong Hou, Xiaofan Zhang, Lingping He, Kefei Song, Peng Zhang, Jing-Song Wang, Xiaoxin Zhang, Weiguo Zong, Jinping Dun, Hui Tian, Yuanyong Deng

14 pages, 6 figures, Accepted for publication in RAA (Research in Astronomy and Astrophysics)

The extreme ultraviolet (EUV) observations are widely used in solar activity research and space weather forecasting since they can observe both the solar eruptions and the source regions of the solar wind. Flat field processing is indispensable to remove the instrumental non-uniformity of a solar EUV imager in producing high-quality scientific data from original observed data. Fengyun-3E (FY-3E) is a meteorological satellite operated in Sun-synchronous orbit, and the routine EUV imaging data from the Solar X-ray and Extreme Ultraviolet Imager (X-EUVI) onboard FY-3E has the characteristics of concentric rotation. Taking advantage of the concentric rotation, we propose a post-hoc flat field measurement method for its EUV 195 channel in this paper. This method removes small-scale and time-varying component of the coronal activities by taking the median value for each pixel along the time axis of a concentric rotation data cube, and then derives large-scale and invariable component of the quiet coronal radiation, and finally generates a flat field image. Analysis shows that our method is able to measure the instrumental spot-like non-uniformity possibly caused by contamination on the detector, which mostly disappears after the in-orbit self-cleaning process. It can also measure the quasi-periodic grid-like non-uniformity, possibly from the obscuration of the support mesh on the rear filter. After flat field correction, these instrumental non-uniformities from the original data are effectively removed. X-EUVI 195 data after dark and flat field corrections are consistent with the 193 channel data from SDO/AIA, verifying the suitability of the method. Our method is not only suitable for FY-3E/X-EUVI but also a candidate method for the flat field measurement of future solar EUV telescopes.

C. P. Folsom, R. Ignace, C. Erba, R. Casini, T. del Pino Alemán, K. Gayley, K. Hobbs, R. Manso Sainz, C. Neiner, V. Petit, M. E. Shultz, G. A. Wade

Submitted to Astrophysics and Space Science. arXiv admin note: text overlap with arXiv:2111.06434

Magnetic fields are important for stellar photospheres and magnetospheres, influencing photospheric physics and sculpting stellar winds. Observations of stellar magnetic fields are typically made in the visible, although infrared observations are becoming common. Here we consider the possibility of directly detecting magnetic fields at ultraviolet (UV) wavelengths using high resolution spectropolarimetry, specifically considering the capabilities of the proposed Polstar mission. UV observations are particularly advantageous for studying wind resonance lines not available in the visible, but they can also provide many photospheric lines in hot stars. Detecting photospheric magnetic fields using the Zeeman effect and Least Squares Deconvolution is potentially more effective in the UV due to the much higher density of strong lines. We investigate detecting magnetic fields in the magnetosphere of a star using the Zeeman effect in wind lines, and find that this could be detectable at high S/N in an O or B star with a strong magnetic field. We consider detecting magnetic fields using the Hanle effect in linear polarization, which is complementary to the Zeeman effect, and could be more sensitive in photospheric lines of rapid rotators. The Hanle effect can also be used to infer circumstellar magnetism in winds. Detecting the Hanle effect requires UV observations, and a multi-line approach is key for inferring magnetic field properties. This demonstrates that high resolution spectropolarimetry in the UV, and the proposed Polstar mission, has the potential to greatly expand our ability to detect and characterize magnetic fields in and around hot stars.

Isabelle Tanseri, Steffen Hagstotz, Sunny Vagnozzi, Elena Giusarma, Katherine Freese

31 pages, 7 sub-figures arranged in 6 figures, 2 tables

We revisit cosmological constraints on the sum of the neutrino masses $\Sigma m_\nu$ from a combination of full-shape BOSS galaxy clustering [$P(k)$] data and measurements of the cross-correlation between Planck Cosmic Microwave Background (CMB) lensing convergence and BOSS galaxy overdensity maps [$C^{\kappa \text{g}}_{\ell}$], using a simple but theoretically motivated model for the scale-dependent galaxy bias in auto- and cross-correlation measurements. We improve upon earlier related work in several respects, particularly through a more accurate treatment of the correlation and covariance between $P(k)$ and $C^{\kappa \text{g}}_{\ell}$ measurements. When combining these measurements with Planck CMB data, we find a 95% confidence level upper limit of $\Sigma m_\nu<0.14\,{\rm eV}$, while slightly weaker limits are obtained when including small-scale ACTPol CMB data, in agreement with our expectations. We confirm earlier findings that (once combined with CMB data) the full-shape information content is comparable to the geometrical information content in the reconstructed BAO peaks given the precision of current galaxy clustering data, discuss the physical significance of our inferred bias and shot noise parameters, and perform a number of robustness tests on our underlying model. While the inclusion of $C^{\kappa \text{g}}_{\ell}$ measurements does not currently appear to lead to substantial improvements in the resulting $\Sigma m_{\nu}$ constraints, we expect the converse to be true for near-future galaxy clustering measurements, whose shape information content will eventually supersede the geometrical one.

V. V. Bobylev, A. T. Bajkova, G. M. Karelin

19 pages, 5 figures, 4 tables

We have analyzed the kinematics of OB stars from the list by Xiang et al. (2021) that contains $\sim$13 000 single OB stars. For these stars there are photometric distance estimates and proper motions from the Gaia catalogue and line-of-sight velocities from the LAMOST catalogue. Based on a sample of single OB stars and using the photometric distances and proper motions of stars from the Gaia EDR3 catalogue, we have found the group velocity components $(U_\odot,V_\odot,W_\odot)=(9.63,9.93,7.45)\pm(0.27,0.34,0.10)$ km s$^{-1}$, and the following parameters of the angular velocity of Galactic rotation: $\Omega_0=29.20\pm0.18$ km s$^{-1}$ kpc$^{-1}$, $\Omega^{'}_0=-4.150\pm0.046$ km s$^{-1}$ kpc$^{-2}$ and $\Omega^{''}_0=0.795\pm0.018$ km s$^{-1}$ kpc$^{-3}$, where the error per unit weight $\sigma_0$ is 9.56 km s$^{-1}$ and $V_0=236.5\pm3.3$ km s$^{-1}$ (for the adopted $R_0=8.1\pm0.1$ kpc). Based on the same OB stars, we have found the residual velocity dispersions $(\sigma_1,\sigma_2,\sigma_3)=(15.13,9.69,7.98)\pm(0.07,0.05,0.04)$ km s$^{-1}$. We show that using the line-of-sight velocities increases significantly the space velocity dispersion and leads to a biased estimate of the velocity $U_\odot$. A comparison of the distances scales used has shown that the photometric distances from Xiang et al. (2021) should be lengthened approximately by 10%.

P. Gavras, L. Rimoldini, K. Nienartowicz, G. Jevardat de Fombelle, B. Holl, P. Ábrahám, M. Audard, M. Carnerero, G. Clementini, J. De Ridder, E. Distefano, P. Garcia-Lario, A. Garofalo, Á. Kóspál, K. Kruszyńska, M. Kun, I. Lecoeur-Taïbi, G. Marton, T. Mazeh, N. Mowlavi, C. Raiteri, V. Ripepi, L. Szabados, S. Zucker, L. Eyer

This paper is part of Gaia Data Release 3 (DR3). Submitted to A&A

Context. In the current ever increasing data volumes of astronomical surveys, automated methods are essential. Objects of known classes from the literature are necessary for training supervised machine learning algorithms, as well as for verification/validation of their results. Aims.The primary goal of this work is to provide a comprehensive data set of known variable objects from the literature cross-matched with \textit{Gaia}~DR3 sources, including a large number of both variability types and representatives, in order to cover as much as possible sky regions and magnitude ranges relevant to each class. In addition, non-variable objects from selected surveys are targeted to probe their variability in \textit{Gaia} and possible use as standards. This data set can be the base for a training set applicable in variability detection, classification, and validation. MethodsA statistical method that employed both astrometry (position and proper motion) and photometry (mean magnitude) was applied to selected literature catalogues in order to identify the correct counterparts of the known objects in the \textit{Gaia} data. The cross-match strategy was adapted to the properties of each catalogue and the verification of results excluded dubious matches. Results.Our catalogue gathers 7\,841\,723 \textit{Gaia} sources among which 1.2~million non-variable objects and 1.7~million galaxies, in addition to 4.9~million variable sources representing over 100~variability (sub)types. Conclusions.This data set served the requirements of \textit{Gaia}'s variability pipeline for its third data release (DR3), from classifier training to result validation, and it is expected to be a useful resource for the scientific community that is interested in the analysis of variability in the \textit{Gaia} data and other surveys.

Kastytis Zubovas, Jonas Bialopetravičius, Monika Kazlauskaitė

15 pages, 6 figures in main text. Accepted for publication in MNRAS

Large-scale outflows driven by active galactic nuclei (AGN) can have a profound influence on their host galaxies. The outflow properties themselves depend sensitively on the history of AGN energy injection during the lifetime of the outflow. Most observed outflows have dynamical timescales longer than the typical AGN episode duration, i.e. they have been inflated by multiple AGN episodes. Here, we present a neural-network based approach to inferring the most likely duty cycle and other properties of AGN based on the observable properties of their massive outflows. Our model recovers the AGN parameters of simulated outflows with typical errors $< 25\%$. We apply the method to a sample of 59 real molecular outflows and show that a large fraction of them have been inflated by AGN shining with a rather high duty cycle $\delta_{\rm AGN} > 0.2$. This result suggests that nuclear activity in galaxies is clustered hierarchically in time, with long phases of more frequent activity composed of many short activity episodes. We predict that $\sim \! 19\%$ of galaxies should have AGN-driven outflows, but half of them are fossils - this is consistent with currently available data. We discuss the possibilities to investigate AGN luminosity histories during outflow lifetimes and suggest ways to use our software to test other physical models of AGN outflows. The source code of all of the software used here is made public.

Núria Salvador-Rusiñol, Anna Ferré-Mateu, Alexandre Vazdekis, Michael A. Beasley

6 pages, 3 figures

Near-ultraviolet (NUV) spectroscopic studies have suggested that passively evolving massive, early-type galaxies host sub-one percent fractions of young stars in their innermost regions. We shed light on the origin of these stars by analysing NGC 1277, a widely studied nearby prototypical massive compact relic galaxy. These are rare galaxies that have survived without experiencing significant size evolution via accretion and mergers since their formation at high redshift. We obtain a spectrum in the UV range within the central 1 kpc region of NGC 1277. We compare a carefully selected set of optical and NUV line-strengths to model predictions with star formation histories characteristic of massive galaxies. We find a 0.8% mass fraction of young stars in the centre of NGC 1277, similar to that found in massive early-type galaxies. Given the limited accretion history of NGC 1277, these results favour an intrinsic, in-situ, process triggering star formation at later epochs. Our results suggest a general constraint on the amount of young stars in the cores of massive early-type galaxies. This amount should be assumed as an upper limit for the young stellar contribution in massive galaxies, as there might be present other contributions from evolved stars.

M. Sugizaki, T. Mihara, K. Kobayashi, H. Negoro, M. Shidatsu, S.N. Pike, W. Iwakiri, S. Urabe, M. Serino, N. Kawai, M. Nakajima, J. A. Kennea, Z. Liu

13 pages, 10 figures, accepted for publication in PASJ

We report on the discovery of a new supergiant fast X-ray transient (SFXT), MAXI J0709$-$159, and its identification with LY CMa (also known as HD 54786). On 2022 January 25, a new flaring X-ray object named MAXI J0709$-$159, was detected by Monitor of All-sky X-ray Image (MAXI). Two flaring activities were observed in the two scans of $\sim 3$ hours apart, where the 2-10 keV flux reached $5\times 10^{-9}$ erg cm$^{-2}$ s$^{-1}$. During the period, the source exhibited a large spectral change suggesting that the absorption column density $N_\mathrm{H}$ increased from $10^{22}$ cm$^{-2}$ to $10^{23}$ cm$^{-2}$. NuSTAR follow-up observation on January 29 identified a new X-ray source with a flux of $6\times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ at the position consistent with LY CMa, which has been identified as B supergiant as well as Be star, located at the 3 kpc distance. The observed X-ray activity characterized by the short ($\lesssim$ several hours) duration, the rapid ($\lesssim$ a few seconds) variabilities accompanied with spectral changes, and the large luminosity swing ($10^{32}$-$10^{37}$ erg s$^{-1}$) agree with those of SFXT. On the other hand, optical spectroscopic observations of LY CMa revealed a broad $H\alpha$ emission line, which may indicate the existence of a Be circumstellar disk. These obtained results suggest that the optical companion, LY CMa, certainly has a complex circumstellar medium including dense clumps.

Konstantinos Karampelas, James A. McLaughlin, Gert J. J. Botha, Stéphane Régnier

19 pages, 14 figures, accepted for publication in ApJ

Oscillatory reconnection can manifest through the interaction between the ubiquitous MHD waves and omnipresent null points in the solar atmosphere and is characterized by an inherent periodicity. In the current study, we focus on the relationship between the period of oscillatory reconnection and the strength of the wave pulse initially perturbing the null point, in a hot coronal plasma. We use the PLUTO code to solve the fully compressive, resistive MHD equations for a 2D magnetic X-point. Using wave pulses with a wide range of amplitudes, we perform a parameter study to obtain values for the period, considering the presence and absence of anisotropic thermal conduction separately. In both cases, we find that the resulting period is independent of the strength of the initial perturbation. The addition of anisotropic thermal conduction only leads to an increase in the mean value for the period, in agreement with our previous study. We also consider a different type of initial driver and we obtain an oscillation period matching the independent trend previously mentioned. Thus, we report for the first time on the independence between the type and strength of the initializing wave pulse and the resulting period of oscillatory reconnection in a hot coronal plasma. This makes oscillatory reconnection a promising mechanism to be used within the context of coronal seismology.

Amit Kashi, Amir Michaelis, Yarden Kaminetsky

Accepted for publication in MNRAS

We carry out a numerical experiment of ejecting winds in a massive colliding wind binary system, and quantifying the accretion onto the secondary star under different primary mass loss rates. We set a binary system comprising a Luminous Blue Variable (LBV) as the primary and a Wolf-Rayet (WR) star as the secondary, and vary the mass loss rate of the LBV to obtain different values of wind momentum ratio $\eta$. Our simulations include two sets of cases: one where the stars are stationary, and one that includes the orbital motion. As $\eta$ decreases the colliding wind structure moves closer to the secondary. We find that for $\eta \lesssim 0.05$ the accretion threshold is reached and clumps which originate by instabilities are accreted onto the secondary. For each value of $\eta$ we calculate the mass accretion rate and identify different regions in the $\dot{M}_{\rm acc}$ - $\eta$ diagram. For $0.001 \lesssim \eta \lesssim 0.05$ the accretion is sub- Bondi-Hoyle-Lyttleton (BHL) and the average accretion rate satisfies the power-law $\dot{M}_{\rm acc} \propto \eta^{-1.73}$ for static stars. The accretion is not continuous but rather changes from sporadic to a larger duty cycle as $\eta$ decreases. For $\eta\lesssim0.001$ the accretion becomes continuous in time and the accretion rate is BHL, up to a factor of 0.4--0.8. The simulations that include the orbital motion give qualitatively similar results, with the steeper power law $\dot{M}_{\rm acc} \propto \eta^{-1.86}$ for the sub-BHL region and lower $\eta$ as an accretion threshold.

About a third of the hot subdwarfs of spectral type B, which are mostly core-helium burning objects on the extreme horizontal branch, are found in close binaries with cool, low-mass stellar, substellar, or white dwarf companions. They can show light variations due to different phenomena. We used light curves from the Transiting Exoplanet Survey Satellite and the \textit{K2} space mission to look for more sdB binaries. Their light curves can be used to study the hot subdwarf primaries and their companions and get orbital, atmospheric, and absolute parameters for those systems. By classifying the light variations and combining this with the fit of the spectral energy distribution, the distance derived by the parallaxes obtained by \textit{Gaia} and the atmospheric parameters, we could derive the nature of the primary and secondary in 122 (75\%) of the known sdB binaries and 82 newly found reflection effect systems. We derive absolute masses, radii, and luminosities for a total of 39 hot subdwarfs with cool, low-mass companions, as well 29 known and newly found sdBs with white dwarf companions. The mass distribution of hot subdwarfs with cool, low-mass stellar and substellar companions differs from those with white dwarf companions, implying they come from different populations. By comparing the period and minimum companion mass distributions, we find that there are several different populations of hot subdwarfs with white dwarf binaries. We also derive the first orbital period distribution for hot subdwarfs with cool, low-mass stellar or substellar systems selected from light variations instead of radial velocity variations. It shows a period distribution from 1.5 hours to 35 hours compared to the distribution of hot subdwarfs with white dwarfs, which ranges from 1 hour to 30 days. These period distributions can be used to constrain the previous common envelope phase.

Morgan T. Chidester, Ebraheem Farag, F.X. Timmes

This manuscript has been accepted for publication to the Astrophysical Journal. 13 pages, 9 figures. Presented at the AAS 240th summer conference in Pasadena. Zenodo link for our figures, data, etc. found here: 10.5281/zenodo.6779983

We seek signatures of the current experimental $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate probability distribution function in the pulsation periods of carbon-oxygen white dwarf models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it forms during the evolution of low-mass stars under radiative helium burning conditions, mitigating the impact of convective mixing processes. We make direct quantitative connections between the pulsation periods of the identified trapped g-modes in variable WD models and the current experimental $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate probability distribution function. We find an average spread in relative period shifts of $\Delta P/P \simeq \pm$ 2\% for the identified trapped g-modes over the $\pm$ 3$\sigma$ uncertainty in the $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate probability distribution function -- across the effective temperature range of observed DAV and DBV white dwarfs and for different white dwarf masses, helium shell masses, and hydrogen shell masses. The g-mode pulsation periods of observed white dwarfs are typically given to 6-7 significant figures of precision. This suggests that an astrophysical constraint on the $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate could, in principle, be extractable from the period spectrum of observed variable white dwarfs.

Chris Hamilton, Roman R. Rafikov

30 pages, 12 figures. To be submitted to ApJ. Comments welcome

The dozens of compact object mergers detected by LIGO/Virgo raise a key theoretical question: how do initially wide binaries shrink sufficiently quickly that they are able to merge via gravitational wave (GW) radiation within a Hubble time? One promising class of answers involves secular driving of binary eccentricity by some external tidal perturbation. This perturbation can arise due to the presence of a tertiary point mass, in which case the system exhibits Lidov-Kozai (LK) dynamics, or it can stem from the tidal field of the stellar cluster in which the binary orbits. While these secular tide-driven mechanisms have been studied exhaustively in the case of no GW emission, when GWs are included the dynamical behavior is still incompletely understood. In this paper we consider compact object binaries driven to merger via high eccentricity excitation by (doubly-averaged, test-particle quadrupole level) cluster tides - which includes LK-driven mergers as a special case - and include the effects of both general relativistic precession and GW emission. We provide for the first time an analytical understanding of the different evolutionary stages of the binary's semimajor axis, secular oscillation timescale, and phase space structure all the way to merger. Our results will inform future population synthesis calculations of compact object binary mergers from hierarchical triples and stellar clusters.

Laura A. Hayes, Sophie Musset, Daniel M üller, S äm Krucker

18 pages, 11 figures, Book Chapter for Handbook of X-ray and Gamma-ray Astrophysics

The Spectrometer/Telescope for Imaging X-rays (STIX) is one of the 10 instruments on-board the scientific payload of ESA's Solar Orbiter mission. STIX provides hard X-ray imaging spectroscopy in the 4-150~keV energy range, observing hard X-ray bremsstrahlung emission from the Sun. These observations provide diagnostics of the hottest thermal plasmas ($>$10~MK) and information on the non-thermal energetic electrons accelerated above 10~keV during solar flares. STIX has a spectral resolution of 1~keV, and employs the use of in-direct bi-grid Fourier imaging to spatially locate hard X-ray emission. Given that STIX provides critical information about accelerated electrons at the Sun through hard X-ray diagnostics, it is a powerful contribution to the Solar Orbiter suite and has a significant role to explore the dynamics of solar inputs to the heliosphere. This chapter describes the STIX instrument, its design, objectives, first observations and outlines the new perspectives STIX provides over the mission lifetime of Solar Orbiter.

Susumu Inoue, Matteo Cerruti, Kohta Murase, Ruo-Yu Liu

13 pages including supplemental material, for submission to PRL

Powerful winds with wide opening angles, likely driven by accretion disks around black holes (BHs), are observed in the majority of active galactic nuclei (AGN) and can play a crucial role in AGN and galaxy evolution. If protons are accelerated in the wind near the BH via diffusive shock acceleration, p-gamma processes with AGN photons can generate neutrinos as well as pair cascade emission from the gamma-ray to radio bands. The TeV neutrinos tentatively detected by IceCube from the obscured Seyfert galaxy NGC 1068 can be interpreted consistently if the shock velocity is appreciably lower than the local escape velocity, which may correspond to a failed, line-driven wind that is physically well motivated. Although the p-gamma-induced cascade is gamma-gamma-attenuated above a few MeV, it can still contribute significantly to the sub-GeV gamma rays observed from NGC 1068. At higher energies, gamma rays can arise via $pp$ processes from a shock where an outgoing wind impacts the obscuring torus, along with some observable radio emission. Tests and implications of this model are discussed. Neutrinos and gamma rays may offer unique probes of AGN wind launching sites, particularly for objects obscured in other forms of radiation.

G. Rodrigues-da-Silva, L. G. Medeiros

30 pages plus appendices and references (single-column), 4 figures

Higher-order theories of gravity are extensions to general relativity (GR) motivated mainly by high-energy physics searching for GR ultraviolet completeness. They are characterized by the inclusion of correction terms in the Einstein-Hilbert action that leads to higher-order field equations. In this paper, we propose investigating inflation due to the GR extension built with all correction terms up to the second-order involving only the scalar curvature $R$, namely, $R^{2}$, $R^{3}$, $R\square R$. We investigate inflation within the Friedmann cosmological background, where we study the phase space of the model, as well as explore inflation in slow-roll leading-order. Furthermore, we describe the evolution of scalar perturbations and properly establish the curvature perturbation. Finally, we confront the proposed model with recent observations from Planck, BICEP3/Keck, and BAO data.

Francis Fortin, Federico Garcia, Sylvain Chaty

19 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

High-Mass X-ray Binaries (HMXBs) are produced after the first supernova event in a massive binary. These objects are intrinsically young, and can suffer from a significant natal kick. As such, the progenitors of HMXBs are likely to have formed away from the current location of the X-ray emitting systems. We aim to find the birthplace of the known HMXBs of our Milky Way. Specifically, we want to answer the question whether the formation of HMXBs can be associated to open stellar clusters and/or Galactic spiral structures, and infer from that the time elapsed since the first supernova event. We use astrometric data from the Gaia EDR3 to initialize the position and velocity of each known HMXBs from the Galaxy, and integrate their motion back in time. In parallel, we perform the same calculations on a sample of 1381 open clusters detected by Gaia as well as for four Galactic spiral arms which shape and motion have also been recently modelled using Gaia data. We report on all the encounter candidates between HMXBs and clusters or spiral arms in the past 100 Myr. In our sample of 26 HMXBs, we infer that 7 were born in clusters, 8 were born near a Galactic spiral arm, and conclude that 7 others could have formed isolated. The birthplaces of the remaining 4 HMXBs are still inconclusive due to a combination of great distance, poor astrometric data and lack of known open cluster in the vicinity. We provide the kinematical age since supernova of 15 HMXBs. The astrometry from Gaia and the orbit integration we employ are effective at finding the birthplaces of HMXBs in the Milky Way. By considering the biases in our data and method, we find it is likely that the progenitors of HMXBs preferentially formed alongside other massive stars in open clusters.

C. Baumgartner, A. C. Birch, H. Schunker, R.H. Cameron, L. Gizon

accepted in A&A

The twist of the magnetic field above a sunspot is an important quantity in solar physics. For example, magnetic twist plays a role in the initiation of flares and coronal mass ejections (CMEs). Various proxies for the twist above the photosphere have been found using models of uniformly twisted flux tubes, and are routinely computed from single photospheric vector magnetograms. One class of proxies is based on $\alpha_z$, the ratio of the vertical current to the vertical magnetic field. Another class of proxies is based on the so-called twist density, $q$, which depends on the ratio of the azimuthal field to the vertical field. However, the sensitivity of these proxies to temporal fluctuations of the magnetic field has not yet been well characterized. We aim to determine the sensitivity of twist proxies to temporal fluctuations in the magnetic field as estimated from time-series of SDO/HMI vector magnetic field maps. To this end, we introduce a model of a sunspot with a peak vertical field of 2370 Gauss at the photosphere and a uniform twist density $q= -0.024$ Mm$^{-1}$. We add realizations of the temporal fluctuations of the magnetic field that are consistent with SDO/HMI observations, including the spatial correlations. Using a Monte-Carlo approach, we determine the robustness of the different proxies to the temporal fluctuations. The temporal fluctuations of the three components of the magnetic field are correlated for spatial separations up to 1.4 Mm (more than expected from the point spread function alone). The Monte-Carlo approach enables us to demonstrate that several proxies for the twist of the magnetic field are not biased in each of the individual magnetograms. The associated random errors on the proxies have standard deviations in the range between $0.002$ and $0.006$ Mm$^{-1}$, which is smaller by approximately one order of magnitude than the mean value of $q$.

Hervé Le Coroller, Mathias Nowak, Kevin Wagner, Markus Kasper, Gael Chauvin, Celia Desgrange, Simon Conseil

9 pages, 11 figures, K-Stacker github link

Keplerian-Stacker is an algorithm able to combine multiple observations acquired at different epochs taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. In 2019, a total of 100 hours of observation were allocated to VLT VISIR-NEAR, a collaboration between ESO and Breakthrough Initiatives, to search for low mass planets in the habitable zone of the Alpha Cen AB binary system. A weak signal (S/N = 3) was reported around Alpha Cen A, at a separation of 1.1 a.u. which corresponds to the habitable zone. We have re-analysed the NEAR data using K-Stacker. This algorithm is a brute-force method able to find planets in time series of observations and to constrain their orbital parameters, even if they remain undetected in a single epoch. We scanned a total of about 3.5e+5 independent orbits, among which about 15 % correspond to fast moving orbits on which planets cannot be detected without taking into account the orbital motion. We find only a single planet candidate, which matches the C1 detection reported in Wagner et al. 2021. Despite the significant amount of time spent on this target, the orbit of this candidate remains poorly constrained due to these observations being closely distributed in 34 days. We argue that future single-target deep surveys would benefit from a K-Stacker based strategy, where the observations would be split over a significant part of the expected orbital period to better constrain the orbital parameters. This application of K-Stacker on high contrast imaging data in the mid-infrared demonstrates the capability of this algorithm to aid in the search for Earth-like planets in the habitable zone of the nearest stars with future instruments of the E-ELT such as METIS.

Florian Peißker, Andreas Eckart, Michal Zajaček, Silke Britzen

33 pages, 23 figures, published at ApJ

The ongoing monitoring of the Galactic center and Sgr A*, the central supermassive black hole, produces surprising and unexpected findings. This goes hand in hand with the technical evolution of ground- and space-based telescopes and instruments, but also with the progression of image filter techniques such as the Lucy Richardson algorithm. As we continue to trace the members of the S cluster close to Sgr A* on their expected trajectory around the supermassive black hole, we present the finding of a new stellar source, which we call S4716. The newly found star orbits SgrA* in about 4.0 yr and can be detected with NIRC2 (Keck), OSIRIS (Keck), SINFONI (VLT), NACO (VLT), and GRAVITY (VLTI). With a periapse distance of about 100 au, S4716 shows an equivalent distance toward Sgr A* as S4711. These fast-moving stars undergo a similar dynamical evolution, since S4711-S4716 share comparable orbital properties. We will furthermore draw a connection between the recent finding of a new faint star called S300 and the data presented here. Additionally, we observed a blend-star event with S4716 and another newly identified S star S148 in 2017.

Dark Energy is the largest fraction of the energy density of our Universe - yet it remains one of the enduring enigmas of our times. Here we show that Dark Energy can be used to solve 2 tantalizing mysteries of the observable universe. We build on existing models of Dark Energy linked to neutrino masses. In these models Dark Energy can undergo Phase Transitions and form Black Holes. Here we look at the implications of the family structure of neutrinos for the phase transitions in dark energy and associated peaks in black hole formation. It has been previously shown that one of these peaks in Black Hole formation is associated with the observed peak in Quasar formation. Here, we predict that there will also be an earlier peak in the Dark Energy Black Holes at high redshifts. These Dark Energy Black Holes formed at high redshifts are Intermediate Mass Black Holes (IMBHs).These Dark Energy Black Holes at large redshift can help explain both the EDGES observations and the observations of large Supermassive Black Holes (SMBHs) at redshifts $z \sim 7$. The existence of an earlier phase of Dark Energy Black Holes take care of some current challenges to theory implied by existing astronomical data and also helps us look for these Dark Energy Black Holes at high redshifts as predicted here through targeted searches for these Black Holes at the redshifts $z \sim 18$. There is a slight dependence of the location of the peak on the lightest neutrino mass - so the peak may be located at a slightly lower value of the redshift. This may actually enable a measurement of the lightest neutrino mass. Finding these Dark Energy Black Holes of Intermediate Mass should be within the reach of upcoming observations - particularly with the James Webb Space Telescope - but perhaps also through the use of other innovative techniques focusing specifically on the redshifts around $z \sim 18$.

Yann Gouttenoire

Submitted to Springer for publication

This book provides a thorough survey of the important questions at the interface between theoretical particle physics and cosmology. After discussing the theoretical and experimental physics revolution that led to the rise of the Standard Model in the past century, this volume reviews all major open puzzles, like the hierarchy problem, the small value of the cosmological constant, the matter-antimatter asymmetry, or the dark matter problem, and presents the state-of-the-art in the proposed solutions, with an extensive bibliography. This manuscript emphasises the fields of thermal dark matter, cosmological first-order phase transitions and gravitational-wave signatures. Comprehensive and encyclopedic, this book could be a rich resource for both researchers and students entering the field. In addition to the reviews composing two third of the material, one third presents the original PhD research work of the author.

Pilar Hernandez, Jacobo Lopez-Pavon, Nuria Rius, Stefan Sandner

54 pages, 18 figures

We revisit the generation of a matter-antimatter asymmetry in the minimal extension of the Standard Model with two singlet heavy neutral leptons (HNL) that can explain neutrino masses. We derive an accurate analytical approximation to the solution of the complete linearized set of kinetic equations, which exposes the non-trivial parameter dependencies in the form of parameterization-independent CP invariants. The identification of various washout regimes relevant in different regions of parameter space sheds light on the relevance of the mass corrections in the interaction rates and clarifies the correlations of baryogenesis with other observables. In particular, by requiring that the measured baryon asymmetry is reproduced, we derive robust upper or lower bounds on the HNL mixings depending on their masses, and constraints on their flavour structure, as well as on the CP-violating phases of the PMNS mixing matrix, and the amplitude of neutrinoless double-beta decay. We also find certain correlations between low and high scale CP phases. Especially emphasizing the testable part of the parameter space we demonstrate that our findings are in very good agreement with numerical results. The methods developed in this work can help in exploring more complex scenarios.

We re-examine sterile neutrino dark matter in gauged $U(1)_{B-L}$ model. Improvements have been made by tracing and careful evaluation of the evolution of the number densities of sterile neutrino $N$ and extra neutral gauge boson $Z'$. As a result, the cosmologically interesting gauge coupling of $U(1)_{B-L}$ for freeze-in sterile neutrino turns out to be smaller than the values reported in the literature, to avoid the overproduction of $Z'$ so that it is consistent with the Big Bang Nucleosynthesis and the Cosmic Microwave Background constraints on the effective number of neutrino species. Similarly, the free streaming length constraints exclude a large parameter space derived in previous studies. In addition to known freeze-in pair production of $N$ from the standard model fermion pairs, we find the case that $N$ is dominantly produced from a pair of $Z'$ at the temperature characterized by the $B-L$ breaking scalar mass. Thus, the naive truncation of the $U(1)_{B-L}$ scalar contribution made in the literature is not valid.

Since it was confirmed two decades ago that the expansion of the Universe is accelerating, it would be of theoretical interests to figure out what is the influence from cosmological constant on detection of stochastic gravitational wave background. This paper studies the overlap reduction functions in de-Sitter space-time for a pair of one-way tracking gravitational wave detectors. It is shown to be non-trivial in an expanding Universe, because the propagation of light along line of sight also has effect on the response of GW detectors. It is found that the expansion of the Universe can enhance the value of magnitude of the overlap reduction functions, when the detector pairs are close to each other. For nanohertz gravitational waves, this effect can dominate the values of overlap reduction functions when the galactic pulsar pairs are separated by milliarcsecond.

Ksh. Newton Singh, Megandhren Govender, Sudan Hansraj, Farook Rahaman

11 Pages, 16 Figures, 1 Table

We construct an exact anisotropic star model with a linear barotropic equation of state and with Finch-Skea potential within the framework of pure Lovelock gravity. A comparison with the corresponding Einstein model in a suitable limit is easily deduced. Evidently higher curvature effects induced by the Lovelock contributions generate lower densities, pressures, surface tensions and anisotropy factors when compared to its Einstein counterpart. The maximum moment of inertia is attained for the Einstein case and hence it may be inferred that Lovelock effects soften the equation of state. The model satisfies various stability tests.

Martina Muratore, Olaf Hartwig, Daniele Vetrugno, Stefano Vitale, William Joseph Weber

We present a study of the use and limits of the TDI null channels for in flight estimation of the LISA instrumental noise. The paper considers how the two main limiting noise sources, test-mass acceleration noise and interferometric phase measurement noise, propagate through different TDI channels: the Michelson combination X that is the most sensitive to gravitational waves, then the less-sensitive combinations $\alpha$, and finally the null channel $\zeta$ which includes the noise information carried by all existing null channels. We find that null channels such as $\zeta$ not only have a reduced sensitivity to the gravitational waves, but also feature a larger degree of cancellation of the test mass acceleration noise relative to the interferometry noise. This severely limits their use in quantifying the low frequency instrumental noise in the Michelson X combination, which is expected to be dominated by acceleration noise. However, we show that one can still use in-flight noise estimations from $\zeta$ to put an upper bound on the considered noises entering in the X channel, which allows to distinguish them from a strong stochastic gravitational wave background.

Vitor S. Barroso, August Geelmuyden, Zack Fifer, Sebastian Erne, Anastasios Avgoustidis, Richard J. A. Hill, Silke Weinfurtner

For the past two hundred years, parametric instabilities have been studied in various physical systems, such as fluids, mechanical devices and even inflationary cosmology. It was not until a few decades ago that this subharmonic unstable response arose as a central mechanism for the thermalisation of the Early Universe, in a theory known as preheating. Here we study a parametrically driven two-fluid interface to simulate the key aspects of inflationary preheating dynamics through the onset of nonlinear Faraday waves. We present a detailed analysis of the effective field theory description for interfacial waves through the factorization properties of higher-order correlations. Despite the intricacies of a damped and highly interacting hydrodynamical system, we show that the scattering of large amplitude Faraday waves is connected to a broadening of primary resonance bands and the subsequent appearance of secondary instabilities as predicted in preheating dynamics.