Papers for Monday, May 23 2022

We measure the anisotropy of the Milky Way stellar halo traced by a dense sample of 18<r<21 mag F-type main sequence turnoff stars using Gaia eDR3 proper motions and new radial velocity measurements published here.

The vast datasets associated with extrasolar systems promise to offer sensitive probes of new physics in the near future. We consider the possibility that such systems may capture primordial black holes (PBHs) or other exotic compact objects, giving rise to unique observational signatures. We estimate the rate of captures by extrasolar systems, accounting for several distinct mechanisms. We find that the capture rate is negligible unless PBHs account for the entirety of dark matter in a narrow mass range just above the threshold of existing constraints from evaporation. In this scenario, luminous evaporating PBHs may be detectable by exoplanet searches.

The cross-correlation between the 21 cm field and the galaxy distribution is a potential probe of the Epoch of Reionization (EoR). The 21 cm signal traces neutral gas in the intergalactic medium and, on large spatial scales, this should be anti-correlated with the high-redshift galaxy distribution which partly sources and tracks the ionized gas. In the near future, interferometers such as the Hydrogen Epoch of Reionization Array (HERA) are projected to provide extremely sensitive measurements of the 21 cm power spectrum. At the same time, the Nancy Grace Roman Space Telescope (Roman) will produce the most extensive catalog to date of bright galaxies from the EoR. Using semi-numeric simulations of reionization, we explore the prospects for measuring the cross-power spectrum between the 21 cm and galaxy fields during the EoR. We forecast a 14$\sigma$ detection between HERA and Roman, assuming an overlapping survey area of 500 deg$^2$, redshift uncertainties of $\sigma_z = 0.01$ (as expected for the high-latitude spectroscopic survey of Ly$\alpha$-emitting galaxies), and an effective Ly$\alpha$ emitter duty cycle of $f_\mathrm{LAE} = 0.1$. Thus the HERA-Roman cross-power spectrum may be used to help verify 21 cm detections from HERA. We find that the shot-noise in the galaxy distribution is a limiting factor for detection, and so supplemental observations using Roman should prioritize deeper observations, rather than covering a wider field of view.

Pre-stellar cores are the basic unit for the formation of stars and stellar systems. The anatomy of the physical and chemical structures of pre-stellar cores is critical for understanding the star formation process. L1544 is a prototypical pre-stellar core, which shows significant chemical differentiation surrounding the dust peak. We aim to constrain the physical conditions at the different molecular emission peaks. This study allows us to compare the abundance profiles predicted from chemical models together with the classical density structure of Bonnor-Ebert (BE) sphere. We conducted multi-transition pointed observations of CH$_{3}$OH, c-C$_{3}$H$_{2}$ and HNCO with the IRAM 30m telescope, towards the dust peak and the respective molecular peaks of L1544. With non-LTE radiative transfer calculations and a 1-dimensional model, we revisit the physical structure of L1544, and benchmark with the abundance profiles from current chemical models. We find that the HNCO, c-C$_{3}$H$_{2}$ and CH$_{3}$OH lines in L1544 are tracing progressively higher density gas, from $\sim$10$^{4}$ to several times 10$^{5}$ cm$^{-3}$. Particularly, we find that to produce the observed intensities and ratios of the CH$_{3}$OH lines, a local gas density enhancement upon the BE sphere is required. This suggests that the physical structure of an early-stage core may not necessarily follow a smooth decrease of gas density profile locally, but can be intercepted by clumpy substructures surrounding the gravitational center. Multiple transitions of molecular lines from different molecular species can provide a tomographic view of the density structure of pre-stellar cores. The local gas density enhancement deviating from the BE sphere may reflect the impact of accretion flows that appear asymmetric and are enhanced at the meeting point of large-scale cloud structures.

Thomson scattered photospheric light is the dominant constituent of the lower solar corona's spectral continuum viewed off-limb at optical wavelengths. Known as the K-corona, it is also linearly polarized. We investigate the possibility of using the a priori polarized characteristics of the K-corona, together with polarized emission lines, to measure and correct instrument-induced polarized crosstalk. First, we derive the Stokes parameters of Thomson scattering of unpolarized light in an irreducible spherical tensor formalism. This allows forward synthesis of the Thomson scattered signal for the more complex scenario of symmetry-breaking features in the incident radiation field, which could limit the accuracy of our proposed technique. For this, we make use of an advanced 3D radiative magnetohydrodynamic coronal model. Together with synthesized polarized signals in the Fe XIII 10746 Angstrom emission line, we find that an ad hoc correction of telescope and instrument-induced polarization crosstalk is possible under the assumption of a non-depolarizing optical system.

Strong magnetic fields can modify the microscopic composition of matter with consequences on stellar macroscopic properties. Within this context, we study, for the first time, the possibility of the appearance of spin-3/2 $\Delta$ baryons in magnetars. We make use of two different relativistic models for the equation of state of dense matter under the influence of strong magnetic fields considering the effects of Landau levels and the anomalous magnetic moment (AMM) proportional to the spin of all baryons and leptons. In particular, we analyze the effects of the AMM of $\Delta$ baryons in dense matter for the first time. {We also obtain global properties corresponding to the EoS models numerically and study the corresponding role of the $\Delta$ baryons.} We find that they are favored over hyperons, which causes an increase in isopin asymmetry and a decrease in spin polarization. We also find that, contrary to what generally occurs when new degrees of freedom are introduced, the $\Delta$s do not make the EoS significantly softer and magnetars less massive. Finally, the magnetic field distribution inside a given star is not affected by the presence of $\Delta$s.

The analysis of the photospheric velocity field is essential for understanding plasma turbulence in the solar surface, which may be responsible for driving processes such as magnetic reconnection, flares, wave propagation, particle acceleration, and coronal heating. Currently, the only available methods to estimate velocities at the solar photosphere transverse to an observer's line of sight infer flows from differences in image structure in successive observations. Due to data noise, algorithms such as local correlation tracking (LCT) may lead to a vector field with wide gaps where no velocity vectors are provided. In this letter, a novel method for image inpainting of highly corrupted data is proposed and applied to the restoration of horizontal velocity fields in the solar photosphere. The restored velocity field preserves all the vector field components present in the original field. The method shows robustness when applied to both simulated and observational data.

We present models designed to quantify the effects of stellar activity on exoplanet transit spectroscopy and atmospheric characterization at low (R = 100) and high (R = 100,000) spectral resolution. We study three model classes mirroring planetary system archetypes: a hot Jupiter around an early-K star (HD 189733 b); a mini-Neptune around an early-M dwarf (K2-18 b); and terrestrial planets around a late M dwarf (TRAPPIST-1). We map photospheres with temperatures and radial velocities (RV) and integrate specific intensity stellar models. We obtain transit spectra affected by stellar contamination, the Rossiter--McLaughlin effect (RME), and center-to-limb variations (CLV). We find that, at low resolution, for later-type stars, planetary water features become difficult to distinguish from contamination. Many distributions of unocculted active regions can induce planetary-like features of similar amplitudes in the case of a late M dwarf. Atmospheric characterization of planets around late-type stars will likely continue to suffer from degeneracy with stellar activity unless active regions' parameters can be constrained using additional information. For the early-K star, stellar contamination mostly manifests itself through a slope at optical wavelengths similar to Rayleigh scattering. In all cases, contamination induces offsets in measured planet radii. At high resolution, we show that we can determine the origin of $\text{H}_2$O and CO detection signals and lift the degeneracy observed at low resolution, provided sufficient planet RV variation during transit and adequate correction for the RME and CLV when required. High-resolution spectroscopy may therefore help resolve issues arising from stellar contamination for favorable systems.

The abundances of mixing--sensitive elements including lithium, [C/N], and 12C/13C are known to change near the red giant branch bump. The explanation most often offered for these alterations is double diffusive thermohaline mixing in the stellar interior. In this analysis, we investigate the ability of thermohaline mixing to explain the observed timing of these chemical depletion events. Recent observational measurements of lithium and [C/N] show that the abundance of lithium decreases before the abundance of [C/N], whereas numerical simulations of the propagation of the thermohaline mixing region computed with MESA show that the synthetic abundances drop simultaneously. We therefore conclude that thermohaline mixing alone cannot explain the distinct events of lithium depletion and [C/N] depletion, as the simultaneity predicted by simulations is not consistent with the observation of separate drops. We thus invite more sophisticated theoretical explanations for the observed temporal separation of these chemical depletion episodes as well as more extensive observational explorations across a range of masses and metallicities.

In the single-degenerate scenario of Type Ia supernovae (SNe Ia), the interaction between high-speed ejected material and the donor star in a binary system is expected to lead to mass being stripped from the donor. A series of multi-dimensional hydrodynamical simulations of ejecta-donor interaction have been performed in previous studies most of which adopt either a simplified analytical model or the W7 model to represent a normal SN Ia explosion. Whether different explosion mechanisms can significantly affect the results of ejecta-donor interaction is still unclear. In this work, we simulate hydrodynamical ejecta interactions with a main-sequence (MS) donor star in two dimensions for two near-Chandrasekhar-mass explosion models of SNe Ia, the W7 and N100 models. We find that about 0.30 and 0.37 Msun of hydrogen-rich material are stripped from a 2.5 Msun donor star in a 2 day orbit by the SN Ia explosion in simulations with the W7 deflagration and N100 delayed-detonation explosion model, respectively. The donor star receives a kick of about 74 and 86 km/s, respectively, in each case. The modal velocity, about 500 km/s, of stripped hydrogen-rich material in the N100 model is faster than the W7 model, with modal velocity of about 350km/s, by a factor 1.4. Based on our results, we conclude that the choice of near-Chandrasekhar-mass explosion model for normal SNe Ia seems to not significantly alter the ejecta-donor interaction for a given main-sequence donor model, at least in 2D.

We study evolution of single subhaloes with their masses of $\sim10^9 M_\odot$ in a Milky-Way-sized host halo for self-interacting dark matter (SIDM) models. We perform ideal dark-matter-only N-body simulations of halo-subhalo mergers by varying self-scattering cross sections (including a velocity-dependent scenario), subhalo orbits, and internal properties of the subhalo. We calibrate a gravothermal fluid model to predict time evolution in spherical mass density profiles of isolated SIDM haloes with the simulations. We find that tidal effects of SIDM subhaloes can be described with a framework developed for the case of collision-less dark matter, but minor revisions are necessary to explain our SIDM simulation results. As long as the cross section is less than $\sim10\, \mathrm{cm}^2/\mathrm{g}$ and a plausible range of subhalo density profiles at redshifts of $\sim2$ is assumed for initial states, our simulations do not exhibit a prominent feature of gravothermal collapse in the subhalo central density for 10 Gyr. We develop a semi-analytic model of SIDM subhaloes in a time-evolving density core of the host with tidal stripping and self-scattering ram pressure effects. Our semi-analytic approach provides a simple, efficient and physically-intuitive prediction of SIDM subhaloes, but further improvements are needed to account for baryonic effects in the host and the gravothermal instability accelerated by tidal stripping effects.

The observed large variation in the abundance of deuterium (D) in the interstellar medium (ISM) suggests that a significant fraction of D may be depleted into polycyclic aromatic hydrocarbons (PAHs). Signatures of deuteration of PAHs are expected to appear most clearly through C-D stretching modes at 4.4--4.7 micron, whose strengths in emission spectra relative to those of C-H stretching modes at 3.3--3.5 micron provide the relative abundance of D to hydrogen (H) in PAHs once we have accurate relative band strengths of both stretching modes. We report experimental results of the band strength of C-D stretching modes relative to C-H. We employ a laboratory analogue of interstellar carbonaceous dust, Quenched Carbonaceous Composite (QCC), and synthesize deuterated QCC (D-QCC) by replacing the starting gas of CH$_4$ of QCC by mixtures of CH$_4$ and CD$_4$ with various ratios. Infrared spectra of D-QCC are taken to estimate the relative band strengths of the stretching modes, while the D/H ratios in the D-QCC samples are measured with a nano-scale secondary ion mass spectrometer (NanoSIMS). We obtain that the relative strength of aromatic and aliphatic C-D to C-H stretches is 0.56 +/- 0.04 and 0.38 +/- 0.01 per D/H, respectively. The ratio for the aromatic stretches is in good agreement with the results of theoretical calculations, while that of aliphatic stretches is smaller than that of aromatic. The present results do not significantly change the D/H ratios in the interstellar PAHs previously estimated from observed spectra.

Tidal disruption of stars in dense nuclear star clusters containing supermassive central black holes (SMBH) is modeled by high-accuracy direct N-body simulation. Stars getting too close to the SMBH are tidally disrupted and a tidal disruption event (TDE) happens. TDEs probe properties of SMBH, their accretion disks, and the surrounding nuclear stellar cluster. In this paper we compare rates of full tidal disruption events (FTDE) with partial tidal disruption events (PTDE). Since a PTDE does not destroy the star, a leftover object emerges; we use the term 'leftover star' for it; two novel effects occur in the simulation: (1) variation of the leftover star's mass and radius, (2) variation of the leftover star's orbital energy. After switching on these two effects in our simulation, the number of FTDEs is reduced by roughly 28%, and the reduction is mostly due to the ejection of the leftover stars from PTDEs coming originally from relatively large distance. The number of PTDEs is about 75% higher than the simple estimation given by Stone et al. (2020), and the enhancement is mainly due to the multiple PTDEs produced by the leftover stars residing in the diffusive regime. We compute the peak mass fallback rate for the PTDEs and FTDEs recorded in the simulation, and find 58% of the PTDEs have peak mass fallback rate exceeding the Eddington limit, and the number of super-Eddington PTDEs is 2.3 times the number of super-Eddington FTDEs.

This paper proposes the idea that the observed dependence of stellar activity cycles on rotation rate can be a manifestation of a stronger dependence on the effective temperature. Observational evidence is recalled and theoretical arguments are given for the presence of cyclic activity in the case of sufficiently slow rotation only. Slow rotation means proximity to the observed upper bound on the rotation period of solar-type stars. This maximum rotation period depends on temperature and shortens for hotter stars. The maximum rotation period is interpreted as the minimum rotation rate for operation of a large-scale dynamo. A combined model for differential rotation and the dynamo is applied to stars of different mass rotating with a rate slightly above the threshold rate for the dynamo. Computations show shorter dynamo cycles for hotter stars. The faster rotation of hotter stars is explained by the larger amplitude of the $\alpha$-effect required for their dynamo. The amplitude of the (cycling) magnetic energy in the computations is proportional to the difference between the rotation period and its upper bound for the dynamo. Stars with moderately different rotation rates can differ significantly in super-criticality of their dynamos and therefore in their magnetic activity, as observed.

Preparing for the first detection of the cosmic 21-cm signal from large-scale interferometer experiments requires rigorous testing of the data analysis and reduction pipelines. To validate that these pipelines do not erroneously remove or add features that can mimic the cosmic signal (e.g. from side-lobes or large-scale power leakage), we require reionisation simulations larger than the experiments primary field of view. For an experiment such as the MWA, with a field of view of $\sim25^{2}$ deg.$^{2}$, this would require a simulation of several Gpcs, which is currently infeasible. To overcome this, we developed a simplified version of the semi-numerical reionisation simulation code 21CMFAST preferencing large volumes over some physical accuracy by assuming linear theory for structure formation. With this, we constructed a 7.5 Gpc comoving volume with voxel resolution of $\sim1.17$ cMpc tailored specifically to the binned spectral resolution of the MWA. This simulation was used for validating the pipelines for the 2020 MWA 21-cm power spectrum (PS) upper limits (Trott et al.). We then use this large-volume simulation to explore: (i) whether smaller volume simulations are biased by the missing large-scale modes, (ii) non-Gaussianity in estimates of the cosmic variance, (iii) biases in the recovered 21-cm PS following foreground wedge removal and (iv) the impact of tiling smaller volume simulations to achieve extremely large volumes. In summary, we find: (i) no biases from missing large-scale power, (ii) significant contribution from non-Gaussianity in the cosmic variance as expected following Mondal et al. (iii) an over-estimate of the 21-cm PS of 10-20 per cent following wedge mode excision for our particular model and (iv) tiling smaller volume simulations under-estimates the large-scale power and also the estimated cosmic variance.

Context. Stellar coronal mass ejections (CMEs) are the primary driver of the exoplanetary space weather and they could affect the habitability of exoplanets. However, detections of possible stellar CME signatures are extremely rare. Aims. This work aims to detect stellar CMEs from time-domain spectra observed through the LAMOST Medium-Resolution Spectroscopic Survey (LAMOST-MRS). Our sample includes 1,379,408 LAMOST-MRS spectra of 226,194 late-type main-sequence stars ($\rm T_{eff} < 6000$ K, $\rm log [g/(cm\ s^{-2})] > 4.0$). Methods. We first identified stellar CME candidates by examining the asymmetries of H$\alpha$ line profiles, and then performed double Gaussian fitting for H$\alpha$ contrast profiles (differences between the CME spectra and reference spectra) of the CME candidates to analyze the temporal variation of the asymmetric components. Results. Three stellar CME candidates were detected on three M dwarfs. The H$\alpha$ and Mg I triplet lines (at 5168.94 {\AA}, 5174.13 {\AA}, 5185.10 {\AA}) of candidate 1 all exhibit a blue-wing enhancement, and the corresponding Doppler shift of this enhancement shows a gradually increasing trend. The H$\alpha$ line also shows an obvious blue-wing enhancement in candidate 2. In candidate 3, the H$\alpha$ line shows an obvious red-wing enhancement, and the corresponding projected maximum velocity exceeds the surface escape velocity of the host star. The lower limit of the CME mass was estimated to be $\sim$$8 \times 10^{17}$ g to $4 \times 10^{18}$ g for these three candidates.

HD49798/RX J0648.0-4418 is a peculiar binary including a hot subdwarf of O spectral type and a compact companion in a 1.55 day orbit. According to the steady spin period derivative $\dot{P}=(-2.17\pm0.01)\times10^{-15} ~\rm s\,s^{-1}$ , the compact object was thought to be a contracting young white dwarf (WD). However, the X-ray luminosity producing by the wind accretion of massive WD is one order of magnitude smaller than the observed value. In this work, we propose an alternative model to account for the observed X-ray luminosity. If the WD was surrounded by a debris disk, the accretion from the debris disk can produce the observed X-ray luminosity and X-ray pulses. Based on the time-varying accretion rate model, the current mass of the debris disk is constrained to be $3.9\times10^{-6}~\rm M_{\odot}$. Comparing with the contraction of the WD, the accretion torque exerting by such a debris disk can only influence the spin evolution of the WD in the early stage. According to the accretion theory, the magnetic field of the WD is constrained to be $\sim (0.7-7)\times10^{4}$ G. The calculated conventional polar cap radius of the WD is larger than the observed emitting-zone radius, which probably originate from the existence of strong and small-scale local magnetic field in the polar cap surface. We expect that further multiband observations on this source can help us to confirm or rule out the existence of a debris disk.

We derive exact analytic solutions for density and velocity fields to all orders in Eulerian perturbation theory for $\Lambda$CDM cosmology. In particular, we show that density and velocity field kernels can be written in a separable form in time and momenta at each perturbative order. The kernel solutions are built from an analytic basis of momentum operators and their time-dependent coefficients, which solve a set of recursive differential equations. We also provide an exact closed perturbative solution for such coefficients, expanding around the (quasi-)EdS approximation. We find that the perturbative solution rapidly converges towards the numerically obtained solutions and its leading order result suffices for any practical requirements. To illustrate our findings, we compute the exact two-loop dark matter density and velocity power spectra in $\Lambda$CDM cosmology. We show that the difference between the exact $\Lambda$CDM and the (quasi-)EdS approximated result can reach the level of several percent. This deviation can be partially mitigated by exploiting the degeneracy with the EFT counterterms. As an additional benefit of our algorithm for the solutions of time-dependent coefficients, the computational complexity of power spectra loops in $\Lambda$CDM is reduced to the same level of the EdS case. In performing the two-loop computation, we devise an explicit method to implement the so-called IR cancellations, as well as the cancellations arising as a consequence of mass and momentum conservation.

Recent high-resolution X-ray spectroscopy revealed possible presence of charge exchange (CX) X-ray emission in supernova remnants (SNRs). Although CX is expected to take place at outermost edges of SNR shells, no significant measurement has been reported so far due to the lack of nearby SNR samples. Here we present an X-ray study of SNR G296.1$-$0.5, which has a complicated multiple-shell structure, with the Reflection Grating Spectrometer (RGS) onboard XMM-Newton. We select two shells in different regions and find that in both regions O VII line shows a high forbidden-to-resonance ($f/r$) ratio that cannot be reproduced by a simple thermal model. Our spectral analysis suggests a presence of CX and the result is also supported by our new radio observation, where we discover evidence of molecular clouds associated with these shells. Assuming G296.1$-$0.5 has a spherical shock, we estimate that CX is dominant in a thin layer with a thickness of 0.2--0.3\% of the shock radius. The result is consistent with a previous theoretical expectation and we therefore conclude that CX occurs in G296.1$-$0.5.

Streaming instability is hypothesized to be triggered at particular protoplanetary disk locations where the volume density of the solid particles is enriched comparable to that of the gas. A ring of planetesimals thus forms when this condition is fulfilled locally. These planetesimals collide with each other and accrete inward drifting pebbles from the outer disk to further increase masses. We investigate the growth of the planetesimals that form in a ring-belt at various disk radii. Their initial mass distributions are calculated based on the formula summarized from the streaming instability simulations. We simulate the subsequent dynamical evolution of the planetesimals with a protoplanetary disk model based either on the minimum mass solar nebula (MMSN) or on the Toomre stability criterion. For the MMSN model, both pebble accretion and planetesimal accretion are efficient at a close-in orbit of $0.3$ AU, resulting in the emergence of several super-Earth mass planets after $1$ Myr. For comparison, only the most massive planetesimals undergo substantial mass growth when they are born at $r{=}3$ AU, while the planetesimals at $r{=}30$ AU experience little or no growth. On the other hand, in the denser Toomre disk, the most massive forming planets can reach Earth mass at $t{=}1$ Myr and reach a mass between that of Neptune and that of Saturn within $3$ Myr at $30$ AU and $100$ AU. Both the pebble and planetesimal accretion rate decrease with disk radial distance. Nevertheless, planetesimal accretion is less pronounced than pebble accretion at more distant disk regions. Taken together, the planets acquire higher masses when the disk has a higher gas density, a higher pebble flux, and/or a lower Stokes number of pebbles.

We study a primordial black hole (PBH) formation scenario based on the Affleck-Dine (AD) mechanism and investigate two PBH mass regions: $M \sim 30 M_\odot$ motivated by the LIGO-Virgo observations of the binary black hole mergers and $M \gtrsim 10^4 M_\odot$ motivated by the observations of supermassive black holes at the center of galaxies. In the previous studies, it has been considered that the inhomogeneous AD baryogenesis generates regions with a large baryon asymmetry, some of which collapse into PBHs. In this paper, we show that this scenario is severely constrained due to the baryon asymmetry remaining outside PBHs, which would spoil the success of the big bang nucleosynthesis. Then, we propose an alternative scenario where the AD leptogenesis results in the inhomogeneous formation of Q-balls with lepton charges, which collapse into PBHs. As a result, we find that our scenario can explain the favorable PBH abundance without conflicting with the observational constraints.

We report XMM-Newton and TESS observations of V496 UMa, an AM Herculis-type cataclysmic variable. The XMM-Newton observation reveals that at times, two poles on the white dwarf accrete simultaneously, but accretion onto the secondary magnetic pole is erratic and can nearly cease in less than one binary orbit (1.5 h). Modelling of the X-ray spectrum during the primary maximum reveals no change in the accretion structures onto the primary pole when accretion onto the secondary pole is disrupted, suggesting that the disruption of accretion onto the secondary pole may be caused by mass-transfer variations from the donor star. The TESS observation, which spanned eight weeks at a two-minute cadence, shows a stable, double-humped orbital modulation due to cyclotron emission from the post-shock region, while the observed times of maximum light show a slow systematic drift that does not correlate with the system's overall brightness.

Catalogs of the Zurich Observatory contain positional information on sunspots, prominences and faculae in late 19th and early 20th centuries. This database is given in handwritten tabular form and was not systematically analysed earlier. It is different from the sunspot number time series made in Zurich and was obtained with a larger telescope. We trained a neural-network model for handwritten text recognition and present the database of reconstructed coordinates. The database obtained connects the earlier observations by Sp\"orer with later programs of the 20th century and supplements the sunspot-group catalogs of the Royal Greenwich Observatory. We also expect that the presented machine-learning approach and its deep capabilities will motivate the processing of a wide bulk of astronomical data, which is still given in non-digitized form or as plain scanned images.

The disk of FU Orionis is marginally resolved with MATISSE, suggesting that the region emitting in the thermal infrared is rather compact. An upper limit of $\sim1.3\pm0.1$ mas (in $L$) can be given for the diameter of the disk region probed in the $L$ band, corresponding to 0.5 au at the adopted Gaia EDR3 distance. This represents the hot, gaseous region of the accretion disk. The $N$-band data indicate that the dusty passive disk is silicate-rich. Only the innermost region of said dusty disk is found to emit strongly in the $N$ band, and it is resolved at an angular size of $\sim5$ mas, which translates to a diameter of about 2 au. The observations therefore place stringent constraints for the outer radius of the inner accretion disk. Dust radiative transfer simulations with RADMC-3D provide adequate fits to the spectral energy distribution from the optical to the submillimeter and to the interferometric observables when opting for an accretion rate $\dot{M}\sim 2\times 10^{-5}\, M_\odot$ yr$^{-1}$ and assuming $M_*=0.6\, M_\odot$. Most importantly, the hot inner accretion disk's outer radius can be fixed at 0.3 au. The outer radius of the dusty disk is placed at 100 au, based on constraints from scattered-light images in the literature. The dust mass contained in the disk is $2.4\times10^{-4}\, M_\odot$, and for a typical gas-to-dust ratio of 100, the total mass in the disk is approximately 0.02 $M_\odot$. We did not find any evidence for a nearby companion in the current interferometric data, and we tentatively explored the case of disk misalignment. For the latter, our modeling results suggest that the disk orientation is similar to that found in previous imaging studies by ALMA. Should there be an asymmetry in the very compact, inner accretion disk, this might be resolved at even smaller spatial scales ($\leq1$ mas).

Context: The space telescope Gaia is dedicated mainly to performing high-precision astrometry, but also spectroscopy and epoch photometry which can be used to study various types of photometric variability. One such variability type is exoplanetary transits. The photometric data accumulated so far have finally matured enough to allow the detection of some exoplanets. Aims: In order to fully exploit the scientific potential of Gaia, we search its photometric data for the signatures of exoplanetary transits. Methods: The search relies on a version of the Box-Least-Square (BLS) method, applied to a set of stars prioritized by machine-learning classification methods. An independent photometric validation was obtained using the public full-frame images of TESS. In order to validate the first two candidates, radial-velocity follow-up observations were performed using the spectrograph PEPSI of the Large Binocular Telescope (LBT). Results: The radial-velocity measurements confirm that two of the candidates are indeed hot Jupiters. Thus, they are the first exoplanets detected by Gaia - Gaia-1b and Gaia-2b. Conclusions: Gaia-1b and Gaia-2b demonstrate that the approach presented in this paper is indeed effective. This approach will be used to assemble a set of additional exoplanet candidates, to be released in Gaia third data release, ensuring better fulfillment of the exoplanet detection potential of Gaia.

In this paper the standard theory of collisional stellar systems is improved by considering the presence of a continuous mass distribution. The calculus of the diffusion coefficients is generalized and a new expression of the Fokker-Planck equation is found for multi-mass systems. A King-like distribution function, which validates the basic assumptions of most multi-mass models for Globular Clusters existing in literature, is obtained.

Galaxy rotation curves are considered to be convincing evidence for dark matter or some dynamically equivalent alternative mechanism. Starting only from the rotation curve data, we present a model independent approach of testing a very general hypothesis that dark matter has the properties of a barotropic fluid. It is shown how the speed of sound squared can be expressed in terms of rotation curve data and their radial derivatives and how model independent constraints can be obtained from the requirements that it is confined between 0 and $c^2$. Using the Milky Way rotation curve data available in the literature, we obtain the constraints on the barotropic fluid speed of sound and illustrate the potential of this approach. Technical challenges, limitations and possible future extensions and improvements of the proposed approach are discussed.

We aim to study the temporal and spectral behaviour of \vum from the optical to the X-ray regimes. We used archival \xmmn and \tes observations obtained in 2017 and 2019 to perform a spectral and timing analysis of the highly variable polar. The light curves of both satellites, TESS and XMM-Newton, reveal a double-humped pattern modulated with the periodicity of $91.058467 \pm 0.00001$ minutes. V496 UMa displays a two-pole accretion geometry in the high accretion state. X-ray spectra from both regions are composed of thermal plasma radiation and soft blackbody components with almost identical temperatures and a total accretion rate of $\dot{M}=1.4(8)\times10^{-11}$ M$_{\odot}$ yr$^{-1}$. The X-ray centers of the humps show longitudinal shifts of -18 and 4 degrees and -172 and -186 degrees size around photometric phase zero, for the main hump and second hump, respectively. The long-term ZTF light curves reveal high and low accretions states. Low-state ZTF and SDSS photometric data are consistent with an $0.8$ M$_{\odot}$ white dwarf at 10000 K and a main-sequence donor star with a spectral type of M5.0 at a \gai-determined distance of 758 pc. V496 UMa is a very bright polar in X-rays when it is in the high state. Due to its unusual geometric structure, mass accretion onto the second accretion pole is interrupted occasionally. This discontinuous behavior does not follow a certain pattern in time and is observed so far only in the high state. The X-ray light curves display clear evidence for an accretion stream at the photometric phase of $\phi=$0.81, which does not show up in optical light curves. An accurate period was derived using the combined TESS and XMM-Newton data, which differs by 3.8 $\sigma$ from published results.

(Abridged) Exoplanetary research has provided us with exciting discoveries of planets around very low mass (VLM) stars (e.g., TRAPPIST-1 and Proxima Centauri). However, current theoretical models strive to explain planet formation in these conditions and do not predict the development of giant planets. Recent high-resolution observations from ALMA of the disk around CIDA 1, a VLM star in Taurus, show substructures hinting at the presence of a massive planet. We aim to reproduce the dust ring of CIDA 1, observed in the dust continuum emission in ALMA Band 7 (0.9 mm) and Band 4 (2.1 mm), along with its $^{12}$CO (J=3-2) and $^{13}$CO (J=3-2) channel maps, assuming the structures are shaped by the interaction of the disk with a massive planet. We seek to retrieve the mass and position of the putative planet. We model the protoplanetary disk with a set of hydrodynamical simulations, varying the mass and locations of the embedded planet. We compute the dust and gas emission using radiative transfer simulations, and, finally, we obtain the synthetic observations treating the images as the actual ALMA observations. Our models indicate that a planet with a minimum mass of $\sim1.4\,\text{M}_\text{Jup}$ orbiting at a distance of $\sim 9-10$ au can explain the morphology and location of the observed dust ring at Band 7 and Band 4. We can reproduce the low spectral index ($\sim 2$) observed where the dust ring is detected. Our synthetic images reproduce the morphology of the $^{12}$CO and $^{13}$CO observed channel maps where the cloud absorption allowed a detection. Applying an empirical relation between planet mass and gap width in the dust, we predict a maximum planet mass of $\sim4 - 8\,\text{M}_\text{Jup}$. Our results suggest the presence of a massive planet orbiting CIDA 1, thus challenging our understanding of planet formation around VLM stars.

This paper explicates the concept of an Intermediate Point (IP), its incorporation as a node along an interplanetary trajectory, and how it permits the determination and optimization of trajectories to interstellar objects (ISOs). IPs can be used to model Solar Oberth Manoeuvres (SOM) as well as Vinfnity Leveraging Manoeuvres (VLM). The SOM has been established theoretically as an important mechanism by which rapid exit from the solar system and intercept of ISOs can both be achieved; the VLM has been demonstrated practically as a method of reducing overall mission Delta-V as well as the Characteristic Energy, C3, at launch. Thus via these two applications, the feasibility of missions to interstellar objects (ISOs) such as 1I/Oumuamua can be analysed. The interplanetary trajectory optimizer tool exploited for this analysis, OITS, permits IP selection as an encounter option along the interplanetary trajectory in question. OITS adopts the assumption of impulsive thrust at discrete points along the trajectory, an assumption which is generally valid for high thrust propulsion scenarios, like chemical or nuclear thermal for instance.

We study the vertical distribution of the highly inclined galaxies from the Continuum Halos in Nearby Galaxies - an EVLA Survey (CHANG-ES). We explore the feasibility of photometrically deriving the HI disk scale-heights from the moment-0 images of the relatively edge-on galaxies with inclination >80 deg, by quantifying the systematic broadening effects and thus deriving correction equations for direct measurements. The corrected HI disk scale-heights of the relatively edge-on galaxies from the CHANG-ES sample show trends consistent with the quasi-equilibrium model of the vertical structure of gas disks. The procedure provide a convenient way to derive the scale-heights and can easily be applied to statistical samples in the future.

This paper presents an interdisciplinary effort aiming to develop and share sustainable knowledge necessary to analyze, understand, and use published scientific results to advance reproducibility in multi-messenger astrophysics. Specifically, we target the breakthrough work associated with the generation of the first image of a black hole, called M87. The image was computed by the Event Horizon Telescope Collaboration. Based on the artifacts made available by EHT, we deliver documentation, code, and a computational environment to reproduce the first image of a black hole. Our deliverables support new discovery in multi-messenger astrophysics by providing all the necessary tools for generalizing methods and findings from the EHT use case. Challenges encountered during the reproducibility of EHT results are reported. The result of our effort is an open-source, containerized software package that enables the public to reproduce the first image of a black hole in the galaxy M87.

3C84 (NGC1275) is one of the brightest radio sources in the mm radio-bands, which led to a plethora of VLBI observations at numerous frequencies over the years. They reveal a two-sided jet structure, with an expanding but not well-collimated parsec-scale jet, pointing southward. High resolution mm-VLBI observations allow the study and imaging of the jet base on the sub-parsec scale. This could facilitate the investigation of the nature of the jet origin, also in view of the previously detected two-railed jet structure and east-west oriented core region seen with RadioAstron at 22 GHz. We produce VLBI images of this core and inner jet region, observed during the past twenty years at 15, 43, and 86 GHz. We determine the kinematics of the inner jet and ejected features at 43 and 86 GHz and compare their ejection times with radio and $\gamma$-ray variability. For the moving jet features, we find an average velocity of $\beta^\textrm{avg}_\textrm{app} = 0.055-0.22$c ($\mu^\textrm{avg} = 0.04-0.18\,$mas/yr). From the time-averaged VLBI images at the three frequencies, we measure the transverse jet width along the bulk flow. On the $\leq 1.5$ parsec scale, we find a clear trend of the jet width being frequency dependent, with the jet being narrower at higher frequencies. This stratification is discussed in the context of a spine-sheath scenario, and is compared to other possible interpretations. From quasi-simultaneous observations at 43 and 86\,GHz we obtain spectral index maps, revealing a time-variable orientation of the spectral index gradient, due to structural variability of the inner jet.

Soft gamma repeaters (SGR) are identified as single neutron stars (NS) inside the Galaxy, or nearby galaxies, with sporadic transient gamma radiation. A total number of discovered SGR, including relative Anomalous X-ray pulsars (AXP), is few tens of objects. Many of them show periodic radiation, connected with NS rotation, with periods 2-12 s. The slow rotation is accompanied by small rate of loss of rotational energy, which is considerably smaller than the observed sporadic gamma ray luminosity, and is many orders less that the luminosity during giant bursts, observed in 4 SGR. Therefore the energy source is usually connected with annihilation of very strong NS magnetic field. Another model is based on release of a nuclear energy stored in the NS nonequilibrium layer. We suggest here an observational test with could distinguish between these two models.

Robust position reconstruction is paramount for enabling discoveries in astroparticle physics as backgrounds are significantly reduced by only considering interactions within the fiducial volume. In this work, we present for the first time a method for position reconstruction using a Bayesian network which provides per interaction uncertainties. We demonstrate the utility of this method with simulated data based on the XENONnT detector design, a dual-phase xenon time-projection chamber, as a proof-of-concept. The network structure includes variables representing the 2D position of the interaction within the detector, the number of electrons entering the gaseous phase, and the hits measured by each sensor in the top array of the detector. The precision of the position reconstruction (difference between the true and expectation value of position) is comparable to the state-of-the-art methods -- an RMS of 0.69 cm, ~0.09 of the sensor spacing, for the inner part of the detector (<60 cm) and 0.98 cm, ~0.12 of the sensor spacing, near the wall of the detector (>60 cm). More importantly, the uncertainty of each interaction position was directly computed, which is not possible with other reconstruction methods. The method found a median 3-$\sigma$ confidence region of 11 cm$^2$ for the inner part of the detector and 21 cm$^2$ near the wall of the detector. We found the Bayesian network framework to be well suited to the problem of position reconstruction. The performance of this proof-of-concept, even with several simplifying assumptions, shows that this is a promising method for providing per interaction uncertainty, which can be extended to energy reconstruction and signal classification.

Accurate wavelength calibration is critical for qualitative and quantitative spectroscopic measurements. Many spectrometers for planetary exploration are equipped with onboard calibration sources. However, such calibration sources are not always available because planetary lander missions often have strong limitations in size and mass. In this study, we propose and validate a wavelength calibration method using solar Fraunhofer lines observed in reflective spectra. As a result, for a visible Raman spectrometer, the accuracy is better than 0.6 cm-1 in 0-4000 cm-1 range, and the magnesium abundance of olivine is estimated more accurately than 2%.

Non-linear growth of structure causes the gravitational potentials to grow with time, and this leaves an imprint on the small-scale temperature fluctuations of the Cosmic Microwave Background (CMB), a signal known as the Rees-Sciama (RS) effect. Building on previous studies, here we investigate the detectability of the RS effect by cross-correlating upcoming CMB and Large-Scale Structure surveys. We include tracers with realistic number density and bias, realistic noise for upcoming and future CMB experiments, and importantly, the contribution from CMB foregrounds to the noise budget. We also derive optimal redshift weights, which are crucial to the detection due to the mismatch between the redshift kernel of the RS effect and the typical redshift distribution of current and upcoming galaxy surveys. In agreement with previous work, we confirm that the signal would in principle be detectable at high significance by "white noise" versions of future CMB experiments, when foregrounds are not included as part of the noise. However, we show that inclusion of foregrounds limits the statistical detectability of the signal: an optimally-weighted high-redshift sample from Rubin Observatory LSST, together with CMB maps from CMB-S4 or CMB-HD, can yield a detection with signal-to-noise 6 - 8, when taking $\ell_{\rm max} = 6000$, provided that foreground-induced biases can be successfully controlled. Improvements are possible if the total power from foregrounds is further reduced, for example by more aggressive masking, or if the signal can be modeled down to smaller scales.

Ultralight axions (ULAs) are promising dark matter candidates that can have a distinct impact on the formation and evolution of structure on nonlinear scales relative to the cold, collisionless dark matter (CDM) paradigm. However, most studies of structure formation in ULA models do not include the effects of self-interactions, which are expected to arise generically. Here, we study how the tidal evolution of solitons is affected by ULA self-interaction strength and sign. Specifically, using the pseudospectral solver UltraDark.jl, we simulate the tidal disruption of self-interacting solitonic cores as they orbit a $10^{11}~M_{\mathrm{\odot}}$ Navarro-Frenk-White CDM host halo potential for a range of orbital parameters, assuming a fiducial ULA particle mass of $10^{-22}\mathrm{eV}$. We find that repulsive (attractive) self-interactions significantly accelerate (decelerate) soliton tidal disruption. We also identify a degeneracy between the self-interaction strength and soliton mass that determines the efficiency of tidal disruption, such that disruption timescales are affected at the $\sim 50\%$ level for variations in the dimensionless ULA self-coupling from $\lambda=-10^{-92}$ to $\lambda=10^{-92}$.

The spatial curvature ($\Omega_K$) of the Universe is one of the most fundamental quantities that could give a link to the early universe physics. In this paper we develop an approximate method to compute the nonlinear matter power spectrum, $P(k)$, for "non-flat" $\Lambda$CDM models using the separate universe (SU) ansatz which states that the effect of the curvature on structure formation is equivalent to that of long-wavelength density fluctuation ($\delta_{\rm b}$) in a local volume in the "flat" $\Lambda$CDM model, via the specific mapping between the background cosmological parameters and redshifts in the non-flat and flat models. By utilizing the fact that the normalized response of $P(k)$ to $\delta_{\rm b}$ (equivalently $\Omega_K$), which describes how the non-zero $\Omega_K$ alters $P(k)$ as a function of $k$, is well approximated by the response to the Hubble parameter $h$ within the flat model, our method allows one to generalize the prediction of $P(k)$ for flat cosmologies via fitting formulae or emulators to that for non-flat cosmologies. We use $N$-body simulations for the non-flat $\Lambda$CDM models with $|\Omega_K|\leq 0.1$ to show that our method can predict $P(k)$ for non-flat models up to $k \simeq 6\,h{\rm Mpc}^{-1}$ in the redshift range $z\simeq [0,1.5]$, to the fractional accuracy within $\sim 1$% that roughly corresponds to requirements for weak lensing cosmology with upcoming surveys. We find that the emulators, those built for flat cosmologies such as EuclidEmulator, can predict the non-flat $P(k)$ with least degradation.

The nature of dark matter is one of the most important unsolved questions in science. Some dark matter candidates do not have sufficient nongravitational interactions to be probed in laboratory or accelerator experiments. It is thus important to develop astrophysical probes which can constrain or lead to a discovery of such candidates. We illustrate this using state-of-the-art measurements of strong gravitationally-lensed quasars to constrain four of the most popular sterile neutrino models, and also report the constraints for other independent methods that are comparable in procedure. First, we derive effective relations to describe the correspondence between the mass of a thermal relic warm dark matter particle and the mass of sterile neutrinos produced via Higgs decay and GUT-scale scenarios, in terms of large-scale structure and galaxy formation astrophysical effects. Second, we show that sterile neutrinos produced through the Higgs decay mechanism are allowed only for mass 26 keV, and GUT-scale scenario 5.3 keV. Third, we show that the single sterile neutrino model produced through active neutrino oscillations is allowed for mass 92 keV, and the 3 sterile neutrino minimal standard model ($\nu$MSM) for mass 16 keV. These are the most stringent experimental limits on these models.

Using idealized models of the accretion disk, we investigate the effects induced by the modified theories of gravity on the annihilation of the neutrino pair annihilation into electron-positron pairs ($\nu{\bar \nu}\to e^-e^+$), occurring near the rotational axis. For the accretion disk, we have considered the models with temperature $T=constant$ and $T\propto r^{-1}$. In both cases, we find that the modified theories of gravity lead to an enhancement, up to more than one order of magnitude with respect to General Relativity, of the rate of the energy deposition rate of neutrino pair annihilation.

Next-generation detectors are expected to be sensitive to postmerger signals from binary neutron star coalescences and thus to directly probe the remnant dynamics. We investigate the scientific potential of postmerger detections with the Einstein Telescope using full Bayesian analyses with the state-of-the-art waveform model ${\tt NRPMw}$. We find that: (i) Postmerger signals with signal-to-noise ratio (SNR) ${\sim}7$ can be confidently detected with a Bayes' factor of $\log{\cal B}\simeq 5$ ($\rm e$-folded) and the posterior distributions report informative measurements already at SNR ${\sim}6$ for some noise realizations. (ii) The postmerger peak frequency $f_2$ can be confidently identified at SNR $7$ with errors of $O(1~{\rm kHz})$, that decrease below $O(100~{\rm Hz})$ for SNR 10. (iii) The remnant's time of collapse to black hole can be constrained to $O(20~{\rm ms})$ at SNR 10. However, the inference can be biased by noise fluctuationsif the latter exceed the signal's amplitude before collapse. (iv) Violations of the EOS-insentive relations for $f_2$ can be detected at SNR $\gtrsim 8$ if the frequency shifts are $\gtrsim 500~{\rm Hz}$; they can be smoking guns for EOS softening effects at extreme densities. However, the $f_2$ measurement can be significantly biased by subdominant frequency components for short-lived remnants. In these cases, an EOS softening might be better inferred from the remnant's earlier collapse.

Dyson spheres are hypothetical megastructures encircling stars in order to harvest most of their energy output. During the 11th edition of the GTOC challenge, participants were tasked with a complex trajectory planning related to the construction of a precursor Dyson structure, a heliocentric ring made of twelve stations. To this purpose, we developed several new approaches that synthesize techniques from machine learning, combinatorial optimization, planning and scheduling, and evolutionary optimization effectively integrated into a fully automated pipeline. These include a machine learned transfer time estimator, improving the established Edelbaum approximation and thus better informing a Lazy Race Tree Search to identify and collect asteroids with high arrival mass for the stations; a series of optimally-phased low-thrust transfers to all stations computed by indirect optimization techniques, exploiting the synodic periodicity of the system; and a modified Hungarian scheduling algorithm, which utilizes evolutionary techniques to arrange a mass-balanced arrival schedule out of all transfer possibilities. We describe the steps of our pipeline in detail with a special focus on how our approaches mutually benefit from each other. Lastly, we outline and analyze the final solution of our team, ACT&Friends, which ranked second at the GTOC 11 challenge.

The Standard Model Higgs potential becomes unstable at large Higgs field values where its quartic coupling becomes negative. While the tunneling lifetime of our current electroweak vacuum is comfortably longer than the age of the universe, quantum fluctuations during inflation might push the Higgs over the barrier, forming patches which might be lethal for our universe. We study the cosmological evolution of such regions and find that, at least in the thin wall approximation, they may be harmless as they collapse due to the backreaction of the Higgs itself. The presence of the Standard Model Higgs instability can provide a novel mechanism to end inflation and to reheat the universe through the evaporation of the black holes left over by the collapse of the Higgs bubbles. The bound on the Hubble rate during inflation may be therefore relaxed.

We describe and implement a procedure for determining the couplings of a Relativistic Mean Field Theory (RMFT) that is optimized for application to neutron star phenomenology. In the standard RMFT approach, the couplings are constrained by comparing the theory's predictions for symmetric matter at saturation density with measured nuclear properties. The theory is then applied to neutron stars which consist of neutron-rich matter at densities ranging up to several times saturation density, which allows for additional astrophysical constraints. In our approach, rather than using the RMFT to extrapolate from symmetric to neutron-rich matter and from finite-sized nuclei to uniform matter we fit the RMFT to properties of uniform pure neutron matter obtained from chiral effective field theory. Chiral effective field theory incorporates the experimental data for nuclei in the framework of a controlled expansion for nuclear forces valid at nuclear densities and enables us to account for theoretical uncertainties when fitting the RMFT. We construct four simple RMFTs that span the uncertainties provided by chiral effective field theory for neutron matter, and are consistent with current astrophysical constraints on the equation of state. Our new RMFTs can be used to model the properties of neutron-rich matter across the vast range of densities and temperatures encountered in simulations of neutron stars and their mergers.

We analyze whether a black hole can exist and survive in a universe that goes through a cosmological bounce. To this end, we investigate a central inhomogeneity embedded in a bouncing cosmological background modeled by the comoving generalized McVittie metric. Contrary to other dynamical metrics available in the literature, this solution allows for the interaction of the central object with the cosmological fluid. We show that the horizons associated with this metric change with cosmic time because they are coupled to the cosmic evolution as the mass of the central object is always proportional to the scale factor: it decreases during contraction and increases during expansion phases. After a full analysis of the causal structure of this spacetime, we determine that a dynamical black hole persists during the contraction, bounce, and expansion of the universe. This result implies that there is a class of bouncing models that admits black holes at all cosmological epochs. If these models are correct approximation to the real universe, then black holes surviving a cosmic collapse could play some role in the subsequent expanding phase.