Papers for Thursday, Apr 08 2021

We present a new catalog of 40502 globular cluster (GC) candidates in NGC 5128 out to a projected radius of $\sim$150 kpc, based on data from the Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS), Gaia Data Release 2, and the NOAO Source Catalog. Ranking these candidates based on the likelihood that they are true GCs, we find that approximately 1900 belong to our top two ranking categories and should be the highest priority for spectroscopic follow-up for confirmation. Taking into account our new data and a vetting of previous GC catalogs, we estimate a total GC population of $1450 \pm 160$ GCs. We show that a substantial number of sources previously argued to be low-velocity GCs are instead foreground stars, reducing the inferred GC velocity dispersion. This work showcases the power of Gaia to identify slightly extended sources at the $\sim 4$ Mpc distance of NGC 5128, enabling accurate identification of GCs throughout the entire extended halo, not just the inner regions that have been the focus of most previous work.

Active region EUV loops are believed to trace a subset of magnetic field lines through the corona. Malanushenko et al. (2009) proposed a method, using loop images and line-of-sight photospheric magnetograms, to infer the three-dimensional shape and field strength along each loop. McCarthy et al. (2019) used this novel method to compute the total magnetic flux interconnecting a pair of active regions observed by SDO/AIA. They adopted the common assumption that each loop had a circular cross section. The accuracy of inferred shape and circularity of cross sections can both be tested using observations of the same loops from additional vantage points as provided by STEREO/EUVI. Here, we use multiple viewing angles to confirm the three-dimensional structure of loops. Of 151 viable cases, 105 (69.5%) matched some form of visible coronal structure when viewed approximately in quadrature. A loop with a circular cross-section should appear of a similar width in different perspectives. In contradiction to this, we find a puzzling lack of correlation between loop diameters seen from different perspectives, even an anti-correlation in some cases. Features identified as monolithic loops in AIA may, in fact, be more complex density enhancements. The 30.5% of reconstructions from AIA which did not match any feature in EUVI might be such enhancements. Others may be genuine loop structures, but with elliptical cross sections. We observe an anti-correlation between diameter and brightness, lending support to the latter hypothesis. Of 13 suitable for width analysis, four loops are consistent with non-circular cross sections, where we find anti-correlation in both comparisons.

Multi-messenger observations of binary neutron star mergers offer a promising path towards resolution of the Hubble constant ($H_0$) tension, provided their constraints are shown to be free from systematics such as the Malmquist bias. In the traditional Bayesian framework, accounting for selection effects in the likelihood requires calculation of the expected number (or fraction) of detections as a function of the parameters describing the population and cosmology; a potentially costly and/or inaccurate process. This calculation can, however, be bypassed completely by performing the inference in a framework in which the likelihood is never explicitly calculated, but instead fit using forward simulations of the data, which naturally include the selection. This is Likelihood-Free Inference (LFI). Here, we use density-estimation LFI, coupled to neural-network-based data compression, to infer $H_0$ from mock catalogues of binary neutron star mergers, given noisy redshift, distance and peculiar velocity estimates for each object. We demonstrate that LFI yields statistically unbiased estimates of $H_0$ in the presence of selection effects, with precision matching that of sampling the full Bayesian hierarchical model. Marginalizing over the bias increases the $H_0$ uncertainty by only $6\%$ for training sets consisting of $O(10^4)$ populations. The resulting LFI framework is applicable to population-level inference problems with selection effects across astrophysics.

We extend results first announced by Franz et al. (1998), that identified vA 351 = H346 in the Hyades as a multiple star system containing a white dwarf. With Hubble Space Telescope Fine Guidance Sensor fringe tracking and scanning, and more recent speckle observations, all spanning 20.7 years, we establish a parallax, relative orbit, and mass fraction for two components, with a period, $P=2.70$y and total mass 2.1 Msun. With ground-based radial velocities from the McDonald Observatory Otto Struve 2.1m telescope Sandiford Spectrograph, and Center for Astrophysics Digital Speedometers, spanning 37 years, we find that component B consists of BC, two M dwarf stars orbiting with a very short period (P_ BC=0.749 days), having a mass ratio M_C/M_B=0.95. We confirm that the total mass of the system can only be reconciled with the distance and component photometry by including a fainter, higher mass component. The quadruple system consists of three M dwarfs (A,B,C) and one white dwarf (D). We determine individual M dwarf masses M_A=0.53+/-0.10 Msun, M_B=0.43+/-0.04Msun, and M_C=0.41+/-0.04Msun. The WD mass, 0.54+/-0.04Msun, comes from cooling models, an assumed Hyades age of 670My, and consistency with all previous and derived astrometric, photometric, and RV results. Velocities from H-alpha and He I emission lines confirm the BC period derived from absorption lines, with similar (He I) and higher (H-alpha) velocity amplitudes. We ascribe the larger H-alpha amplitude to emission from a region each component shadows from the other, depending on the line of sight.

With the dramatic rise in high-quality galaxy data expected from Euclid and Vera C. Rubin Observatory, there will be increasing demand for fast high-precision methods for measuring galaxy fluxes. These will be essential for inferring the redshifts of the galaxies. In this paper, we introduce Lumos, a deep learning method to measure photometry from galaxy images. Lumos builds on BKGnet, an algorithm to predict the background and its associated error, and predicts the background-subtracted flux probability density function. We have developed Lumos for data from the Physics of the Accelerating Universe Survey (PAUS), an imaging survey using 40 narrow-band filter camera (PAUCam). PAUCam images are affected by scattered light, displaying a background noise pattern that can be predicted and corrected for. On average, Lumos increases the SNR of the observations by a factor of 2 compared to an aperture photometry algorithm. It also incorporates other advantages like robustness towards distorting artifacts, e.g. cosmic rays or scattered light, the ability of deblending and less sensitivity to uncertainties in the galaxy profile parameters used to infer the photometry. Indeed, the number of flagged photometry outlier observations is reduced from 10% to 2%, comparing to aperture photometry. Furthermore, with Lumos photometry, the photo-z scatter is reduced by ~10% with the Deepz machine learning photo-z code and the photo-z outlier rate by 20%. The photo-z improvement is lower than expected from the SNR increment, however currently the photometric calibration and outliers in the photometry seem to be its limiting factor.

Recent ALMA observations have found many protoplanetary discs with rings that can be explained by gap-opening planets less massive than Jupiter. Meanwhile, recent studies have suggested that protoplanetary discs should have low levels of turbulence. Past computational work on low-viscosity discs has hinted that these two developments might not be self-consistent because even low-mass planets can be accompanied by vortices instead of conventional double rings. We investigate this potential discrepancy by conducting hydrodynamic simulations of growing planetary cores in discs with various aspect ratios ($H/r=0.04$, 0.06, 0.08) and viscosities ($1.5 \times 10^{-5} \lesssim \alpha \lesssim 3 \times 10^{-4}$), having these cores accrete their gas mass directly from the disc. With $\alpha < 10^{-4}$, we find that sub-Saturn-mass planets in discs with $H/r \le 0.06$ are more likely to be accompanied by dust asymmetries compared to Jupiter-mass planets because they can trigger several generations of vortices in succession. We also find that vortices with $H/r = 0.08$ survive $>6000$ planet orbits regardless of the planet mass or disc mass because they are less affected by the planet's spiral waves. We connect our results to observations and find that the outward migration of vortices with $H/r \ge 0.08$ may be able to explain the cavity in Oph IRS 48 or the two clumps in MWC 758. Lastly, we show that the lack of observed asymmetries in the disc population in Taurus is unexpected given the long asymmetry lifetimes in our low viscosity simulations ($\alpha \sim 2 \times 10^{-5}$), a discrepancy we suggest is due to these discs having higher viscosities.

Swift J004929.5-733107 is an X-ray source in the Small Magellanic Cloud (SMC) that has been reported several times, but the optical counterpart has been unclear due to source confusion in a crowded region of the SMC. Previous works proposed [MA93] 302 as the counterpart, however we show here, using data obtained from the S-CUBED project, that the X-ray positio is inconsistent with that object. Instead we propose a previously unclassified object which has all the indications of being a newly identified Be star exhibiting strong HU emission. Evidence for the presence of significant I-band variability strongly suggest that this is, in fact, a Be type star with a large circumstellar disk. Over 18 years worth of optical monitoring by the OGLE project reveal a periodic modulation at a period of 413d, probably the binary period of the system. A SALT optical spectrum shows strong Balmer emission and supports a proposed spectral classification of B1-3 III-IVe. The X-ray data obtained from the S-CUBED project reveal a time-averaged spectrum well fitted by a photon index = 0.93 pm 0.16. Assuming the known distance to the SMC the flux corresponds to a luminosity 10E35 erg/s. All of these observational facts suggest that this is confirmed as a Be star-neutron star X-ray binary (BeXRB) in the SMC, albeit one with an unusually long binary period at the limits of the Corbet Diagram.

Solar twins are key in different areas of astrophysics, however only just over a hundred were identified and well-studied in the last two decades. In this work, we take advantage of the very precise \textit{Gaia} (DR2/EDR3), Tycho and 2MASS photometric systems to create the Inti survey of new solar twins in the Northern Hemisphere. The spectra of our targets were initially obtained with spectrographs of moderate resolution (ARCES and Goodman spectrographs with $R$ = 31500 and 11930, respectively) to find the best solar twin candidates and then observed at McDonald Observatory with higher resolving power (TS23, $R$ = 60000) and signal-to-noise ratio (SNR $\sim$ 300-500). The stellar parameters were estimated through the differential spectroscopic equilibrium relative to the Sun, which allow us to achieve a high internal precision ($\sigma(T_{\rm{eff}})$ = 15 K, $\sigma(\log g)$ = 0.03 dex, $\sigma$([Fe/H]) = 0.01 dex, and $\sigma(v_{t})$ = 0.03 km s$^{-1}$). We propose a new class of stars with evolution similar to the Sun: \textit{solar proxy}, which is useful to perform studies related to the evolution of the Sun, such as its rotational and magnetic evolution. Its definition is based on metallicity ($-$0.15 dex $\leq$ [Fe/H] $\leq$ +0.15 dex) and mass (0.95 M$_{\odot}$ $\leq$ M $\leq$ 1.05 M$_{\odot}$) constraints, thus assuring that the star follows a similar evolutionary path as the Sun along the main sequence. Based on this new definition, we report 70 newly identified solar proxies, 46 solar analogs and 13 solar-type stars. In addition, we identified 9 \textit{close solar twins} whose stellar parameters are the most similar to those of the Sun.

The discovery of large-amplitude narrowband whistler-mode waves at frequencies of tenths of the electron cyclotron frequency in large numbers both inside ~.3 AU and at ~1 AU provides an answer to longstanding questions about scattering and energization of solar wind electrons. The waves can have rapid nonlinear interactions with electrons over a broad energy range. Counter-propagation between electrons and waves is not required for resonance with the obliquely propagating waves in contrast to the case for parallel propagation. Using a full 3d particle tracing code, we have examined interactions of electrons with energies from 0 eV to 2 keV with whistler-mode waves with amplitudes of 20 mV/m and propagation angles from 0 to 180 degrees to the background magnetic field. Interactions with wave packets and single waves are both modeled based on observations at ~.3 AU and 1 AU. The simulations demonstrate the key role played by these waves in rapid scattering and energization of electrons. Results provide evidence for nonlinear effects, indicating that quasi-linear methods are not adequate for modeling the role of whistlers in the evolution of solar wind electrons. Strong scattering and energization for some initial energy and pitch angle ranges occurs for both counter-propagating and obliquely propagating waves. The strong scattering of strahl electrons counteracts the pitch angle narrowing due to conservation of the first adiabatic invariant as electrons propagate away from the sun into regions of smaller magnetic field. Scattering also produces the hotter isotropic halo. The concomitant limiting of the electron heat flux is also relevant in other astrophysical settings.

We present a brief description of the newly upgraded versions of the programs TLUSTY, version 208, and SYNSPEC, version 54. TLUSTY is used to generate model stellar atmospheres or accretion disks, and SYNSPEC produces detailed synthetic spectra and/or opacity tables. This paper complements published guides that provide a detailed description of earlier versions of the codes, TLUSTY205, and SYNSPEC51. The main upgrades include the flexible construction of opacity tables in SYNSPEC, and their use in producing hybrid models with TLUSTY}, with important species treated in NLTE, while the bulk of opacity of atomic and molecular lines and continua are considered in LTE using a pre-calculated opacity table. There is also a number of additional changes and upgrades that increase the versatility and flexibility of these codes.

Very recently, the Tibet-AS$\gamma$ collaboration reported the detection of $\gamma$ rays from the galactic disk in the energy range of 100 TeV -- 1 PeV. Remarkably, many of these $\gamma$ rays were observed apart from known very high energy (E$>$ 100 GeV) $\gamma$-ray sources. These results are best understood if these diffuse $\gamma$ rays: 1) were produced by a conventional rather than an exotic (i.e. dark matter decay or annihilation) process, 2) have a hadronic rather than a leptonic origin, 3) were produced in impulsive rather than stable sources or, alternatively, in optically thick sources. In addition to that, the detection of the sub-PeV diffuse $\gamma$ rays implies a limit on the flux of neutrinos from the Galactic disk and a lower limit on the rigidity of the cutoff in the Galactic cosmic ray spectrum.

Direct measurements of three-dimensional magnetic fields in the interstellar medium (ISM) are not achievable. However, the anisotropic nature of magnetohydrodynamic (MHD) turbulence provides a novel way of tracing the magnetic fields. Guided by the advanced understanding of turbulence's anisotropy in the Position-Position-Velocity (PPV) space, we extend the Structure-Function Analysis (SFA) to measure both the three-dimensional magnetic field orientation and Alfven Mach number $M_A$, which provides the information on magnetic field strength. Following the theoretical framework developed in Kandel et al. (2016), we find that the anisotropy in a given velocity channel is affected by the inclination angle between the 3D magnetic field direction and the line-of-sight as well as media magnetization. We analyze the synthetic PPV cubes generated by incompressible and compressible MHD simulations. We confirm that the PPV channel's intensity fluctuations measured in various position angles reveal plane-of-the-sky magnetic field orientation. We show that by varying the channel width, the anisotropies of the intensity fluctuations in PPV space can be used to simultaneously estimate both magnetic field inclination angle and strength of total magnetic fields.

We studied the spectral features of Si II $\lambda\lambda$4130, 5972, 6355 and S II W-trough for a large sample of Type Ia supernovae (SNe Ia). We find that in NV (Normal-Velocity) subclass of SNe Ia, these features tend to reach a maximum line strength near maximum light, except for Si II $\lambda$5972. Spectral features with higher excitation energy, such as S II W-trough, are relatively weak and have relatively low velocity. SNe Ia with larger $\Delta$m$_{15}$($B$) tend to have lower velocities especially at phases after maximum light. NV SNe show a trend of increasing line strength with increasing $\Delta$m$_{15}$($B$), while 91T/99aa-like SNe show an opposite trend. Near maximum light, the absorption depth of Si II $\lambda$5972 shows the strongest correlation with $\Delta$m$_{15}$($B$), while at early times the sum of the depths of Si II $\lambda\lambda$4130 and 5972 shows the strongest correlation with $\Delta$m$_{15}$($B$). The overall correlation between velocity and line strength is positive, but within NV SNe the correlation is negative or unrelated. In normal SNe Ia, the velocity-difference and depth-ratio of a longer-wavelength feature to a shorter-wavelength feature tend to increase with increasing $\Delta$m$_{15}$($B$). These results are mostly well explained with atomic physics, but some puzzles remain, possibly related to the effects of the saturation, line competition or other factors.

Context. We present our findings for the HCN/H13CN 1-0 line ratio in the molecular outflow of Arp 220 west based on high-resolution ALMA data. Aims. Molecular gas masses in the outflowing gas of galaxies driven by active galactic nuclei (AGNs) or starbursts are important parameters for understanding the feedback of these latter two phenomena and star-formation quenching. The conversion factor of line luminosities to masses is related to the optical depth of the molecular lines. Methods. Using H13CN 1-0, the isotopic line of HCN 1-0, to obtain the line ratio of HCN/H13CN 1-0 is an ideal way to derive the optical depth of HCN 1-0 in outflowing gas. Results. With the nondetection of H13CN 1-0 in the outflowing gas, a 3{\sigma} lower limit of HCN/H13CN 1-0 line ratio is obtained, which is at least three times higher than that found in the whole of the whole system of Arp 220. The high HCN/H13CN 1-0 line ratio indicates low opacity of HCN 1-0 in the outflowing gas, even though the upper limit of HCN 1-0 opacity obtained here is still not good enough to draw any robust conclusions if HCN 1-0 is optically thin. A lower conversion factor of HCN 1-0 luminosity to dense gas mass in the outflowing gas should be used than that used for the host galaxy of Arp 220.

This work proposes a thermophysical model for realistic surface layers on airless small bodies (RSTPM), for the use of interpreting their multi-epoch thermal lightcurves (e.g WISE/NEOWISE). RSTPM considers real orbital cycle, rotation cycle, rough surface, temperature dependent thermal parameters, as well as contributions of sunlight reflection to observations, hence being able to produce precise temperature distribution and thermal emission of airless small bodies regarding the variations in orbital time scales. Details of the physics, mathematics and numerical algorithms of RSTPM are presented. When used to interpret multi-epoch thermal lightcurves by WISE/NEOWISE, RSTPM can give constraints on the spin orientation and surface physical properties, like mean thermal inertia or mean size of dust grains, roughness fraction, albedo and so on via radiometric procedure. As an application example, we apply this model to the main-belt object (24) Themis, the largest object of the Themis family, which is believed to be the source region of many main-belt comets. We find multi-epoch (2010, 2014-2018) observations of Themis by WISE/NEOWISE, yielding 18 thermal lightcurves. By fitting these data with RSTPM, best-fit spin orientation of Themis is derived to be ($\lambda=137^\circ$, $\beta=59^\circ$) in ecliptic coordinates, the mean radius of dust grains on the surface is estimated to be $\tilde{b}=140^{+500}_{-114}(6\sim640)~\mu$m, indicating the surface thermal inertia to vary from $\sim3\rm~Jm^{-2}s^{-0.5}K^{-1}$ to $\sim60\rm~Jm^{-2}s^{-0.5}K^{-1}$ due to seasonal temperature variation. Further analysis found that thermal light curves of Themis show a weak rotation-phase dependent feature, indicative of heterogeneous thermal properties or imperfections of lightcurve inversion shape model.

To reveal the detail of the internal structure, the relationship between chromospheric activity and the Rossby number, N_R (= rotational period P / convective turnover time tau_c), has been extensively examined for main-sequence stars. The goal of our work is to apply the same methods to pre-main-sequence (PMS) stars and identify the appropriate model of tau_c for them. Yamashita et al. (2020) investigated the relationship between N_R and strengths of the Ca II infrared triplet (IRT; lambda 8498, 8542, 8662 A) emission lines of 60 PMS stars. Their equivalent widths are converted into the emission line to stellar bolometric luminosity ratio (R'). The 54 PMS stars have N_R < 10^{-1.0} and show R' \sim 10^{-4.2} as large as the maximum R' of the zero-age main-sequence (ZAMS) stars. However, because all R' was saturated against N_R, it was not possible to estimate the appropriate tau_c model for the PMS stars. We noticed that Mg I emission lines at 8808 A is an optically thin chromospheric line, appropriate for determination of the adequate tau_c for PMS stars. Using the archive data of the Anglo-Australian Telescope (AAT)/the University College London Echelle Spectrograph (UCLES), we investigated the Mg I line of 52 ZAMS stars. After subtracting photospheric absorption component, the Mg I line is detected as an emission line in 45 ZAMS stars, whose R' is between 10^{-5.9} and 10^{-4.1}. The Mg I line is not saturated yet in "the saturated regime for the Ca II emission lines", i.e. 10^{-1.6} < N_R < 10^{-0.8}. Therefore, the adequate tau_c for PMS stars can be determined by measuring of their R' values.

We present the first southern-hemisphere all-sky imager and radio-transient monitoring system implemented on two prototype stations of the low-frequency component of the Square Kilometre Array. Since its deployment the system has been used for real-time monitoring of the recorded commissioning data. Additionally, a transient searching algorithm has been executed on the resulting all-sky images. It uses a difference imaging technique, and has enabled identification of a wide variety of transient classes, ranging from human-made radio-frequency interference to genuine astrophysical events. Observations at the frequency 159.4 MHz and higher in a single coarse channel (0.926 MHz) were made with 2s time resolution, and multiple nights were analysed. Despite having modest sensitivity (~few Jy/beam), using a single coarse channel and 2-s imaging, the system detected bright transients from PSR B0950+08, proving that it can be used to detect bright transients of an astrophysical origin. The unusual, extreme activity of the pulsar PSR B0950+08 (up to ~155 Jy/beam) was initially detected in a "blind" search in the 2020-04-10/11 data and later assigned to this specific pulsar. The limitations of our data, however, prevent use from making firm conclusions of the effect being due to a combination of refractive and diffractive scintillation or intrinsic emission mechanisms. The system can routinely collect data over many days without interruptions; the large amount of recorded data at 159.4 and 229.7 MHz allowed us to determine a preliminary transient surface density upper limit of $1.32 \times 10^{-9} \text{deg}^{-2}$ for a timescale and limiting flux density of 2s and 42 Jy, respectively. We plan to extend the observing bandwidth to tens of MHz and improve time resolution to tens of milliseconds in order to increase the sensitivity and enable detections of Fast Radio Bursts below 300 MHz.

The disc instability model accounts well for most of the observed properties of dwarf novae and soft X-ray transients, but the rebrightenings, reflares, and echoes occurring at the end of outbursts or shortly after in WZ Sge stars or soft X-ray transients have not yet been convincingly explained by any model. We determine the additional ingredients that must be added to the DIM to account for the observed rebrightenings. We analyse in detail a recently discovered system, TCP J21040470+4631129, which has shown very peculiar rebrightenings, model its light curve using our numerical code including mass transfer variations from the secondary, inner-disc truncation, disc irradiation by a hot white dwarf and, in some cases, the mass-transfer stream over(under)flow. We show that the luminosity in quiescence is dominated by a hot white dwarf that cools down on time scales of months. The mass transfer rate from the secondary has to increase by several orders of magnitudes during the initial superoutburst for a reason that remains elusive, slowly returning to its secular average, causing the observed succession of outbursts with increasing quiescence durations, until the disc can be steady, cold, and neutral; its inner parts being truncated either by the white dwarf magnetic field or by evaporation. The very short, quiescence phases between reflares are reproduced when the mass-transfer stream overflows the disc. Using similar additions to the DIM, we have also produced light curves close to those observed in two WZ Sge stars, the prototype and EG Cnc. Our model successfully explains the reflares observed in WZ Sge systems. It requires, however, the inner disc truncation in dwarf novae to be due not (only) to the white dwarf magnetic field but, as in X-ray binaries, rather to evaporation of the inner disc. A similar model could also explain reflares observed in soft X-ray transients.

We develop a new method based on Gaussian process to reconstruct the mass distribution of binary black holes (BBHs). Instead of prespecifying the formalisms of mass distribution, we introduce a more flexible and nonparametric model with which the distribution can be mainly determined by the observed data. We first test our method with simulated data, and find that it can well recover the injected distribution. Then we apply this method to analyze the data of BBHs' observations from LIGO-Virgo Gravitational-Wave Transient Catalog 2. By reconstructing the chirp mass distribution, we find that there is a peak or a platform locating at $20-30\,M_{\odot}$ rather than a single-power-law-like decrease from low mass to high mass. Moreover, one or two peaks in the chirp mass range of $\mathcal{M}<20\,M_{\odot}$ may be favored by the data. Assuming a mass-independent mass ratio distribution of $p(q)\propto q^{1.4}$, we further obtain a distribution of primary mass, and find that there is a feature locating in the range of $(30, 40)\,M_{\odot}$, which can be related to \textsc{Broken Power Law} and \textsc{Power Law + Peak} distributions described in The LIGO Scientific Collaboration et al. (2020). Besides, the merger rate of BBHs is estimated to $\mathcal{R} = 25.53^{+13.87}_{-8.76}~{\rm Gpc^{-3}~yr^{-1}}$ supposing there is no redshift evolution.

Chemical bimodality of the Milky Way (MW) disc stars constitutes one of the most remarkable properties of MW. The cold accretion theory for the cosmological gas accretion provides one viable explanation to this phenomenon. In this scenario, the rapid cold-mode accretion in the early epoch creates the first generation stars relatively rich in $\alpha$-elements(O,Mg,Si,S,Ca,etc) and later cooling flow produces iron-rich second generation stars, creating the bimodality in the [$\alpha$/Fe] ratio. We employ a cosmologically motivated chemical evolution model for disc galaxies to elucidate the role played by type Ia supernovae (SNIa), which serve as the major source of iron, in the creation of the bimodality. To this end, we divide SNIa into two groups, those formed from the 1st generation stars (the first SNIa) and those formed from the 2nd generation stars (the second SNIa). The model with the first SNIa suppressed during the {\it second} star formation stage produces stars having high [$\alpha$/Fe] in the early phase of this stage, whereas the model which prohibits the second SNIa produces high [$\alpha$/Fe] stars in the late phase. Both models fail to create a well-defined bimodality. We thus conclude that the cooperation of the first and the second SNIa plays a crucial role in creating the bimodality by maintaining rich iron content in the interstellar gas throughout the second star formation stage.

There is a population of stars with velocities in excess of 500 km s$^{-1}$ relative to the Galactic center. Many, perhaps most, of these hyper-velocity stars (HVSs) are B stars, similar to the disk and S stars in a nuclear cluster around a super-massive black hole (SMBH) near $\rm Sgr~A^{\star}$. In the paper I of this series, we showed that the eccentricity of the stars emerged from a hypothetical disk around the SMBH can be rapidly excited by the secular perturbation of its intermediate-mass companion (IMC), and we suggested IRS 13E as a potential candidate for the IMC. Here, we show that this process leads to an influx of stars on parabolic orbits to the proximity of $\rm Sgr~A^{\star}$ on a secular timescale of a few Myr. This timescale is much shorter than the diffusion timescale into the lost cone through either the classical or the resonant relaxation. Precession of the highly-eccentric stars' longitude of periastron, relative to that of the IMC, brings them to its proximity within a few Myr. The IMC's gravitational perturbation scatters a fraction of the stars from nearly parabolic to hyperbolic orbits, with respect to the SMBH. Their follow-up close encounters with the SMBH induce them to escape with hyper-velocity. This scenario is a variant of the hypothesis proposed by Hills based on the anticipated breakup of some progenitor binary stars in the proximity of the SMBH, and its main objective is to account for the limited lifespan of the known HVSs. We generalize our previous numerical simulations of this process with a much wider range of orbital configurations. We demonstrate the robustness and evaluate the efficiency of this channel of HVS formation. From these numerical simulations, we infer observable kinematic properties for the HVSs.

What is habitability? Can we quantify it? What do we mean under the term habitable or potentially habitable planet? With estimates of the number of planets in our Galaxy alone running into billions, possibly a number greater than the number of stars, it is high time to start characterizing them, sorting them into classes/types just like stars, to better understand their formation paths, their properties and, ultimately, their ability to beget or sustain life. After all, we do have life thriving on one of these billions of planets, why not on others? Which planets are better suited for life and which ones are definitely not worth spending expensive telescope time on? We need to find sort of quick assessment score, a metric, using which we can make a list of promising planets and dedicate our efforts to them. Exoplanetary habitability is a transdisciplinary subject integrating astrophysics, astrobiology, planetary science, even terrestrial environmental sciences. We review the existing metrics of habitability and the new classification schemes of extrasolar planets and provide an exposition of the use of computational intelligence techniques to evaluate habitability scores and to automate the process of classification of exoplanets. We examine how solving convex optimization techniques, as in computing new metrics such as CDHS and CEESA, cross-validates ML-based classification of exoplanets. Despite the recent criticism of exoplanetary habitability ranking, this field has to continue and evolve to use all available machinery of astroinformatics, artificial intelligence and machine learning. It might actually develop into a sort of same scale as stellar types in astronomy, to be used as a quick tool of screening exoplanets in important characteristics in search for potentially habitable planets for detailed follow-up targets.

We present X-ray analysis of the ejecta of supernova remnant G350.1$-$0.3 observed with Chandra and Suzaku, and clarify the ejecta's kinematics over a decade and obtain a new observational clue to understanding the origin of the asymmetric explosion. Two images of Chandra X-ray Observatory taken in 2009 and 2018 are analyzed in several methods, and enable us to measure the velocities in the plane of the sky. A maximum velocity is 4640$\pm$290 km s$^{-1}$ (0.218$\pm$0.014 arcsec yr$^{-1}$) in the eastern region in the remnant. These findings trigger us to scrutinize the Doppler effects in the spectra of the thermal emission, and the velocities in the line-of-sight direction are estimated to be a thousand km s$^{-1}$. The results are confirmed by analyzing the spectra of Suzaku. Combining the proper motions and line-of-sight velocities, the ejecta's three-dimensional velocities are $\sim$3000-5000 km s$^{-1}$. The center of the explosion is more stringently constrained by finding the optimal time to reproduce the observed spatial expansion. Our findings that the age of the SNR is estimated at most to be 655 years, and the CCO is observed as a point source object against the SNR strengthen the 'hydrodynamical kick' hypothesis on the origin of the remnant.

High precision astrometry aims at source position determination to a very small fraction of the diffraction image size, in high SNR regime. One of the key limitations to such goal is the optical response variation of the telescope over a sizeable FOV, required to ensure bright reference objects to any selected target. The issue translates into severe calibration constraints, and/or the need for complex telescope and focal plane metrology. We propose an innovative system approach derived from the established TMA telescope concept, extended to achieve high filling factor of an annular field of view around the optical axis of the telescope. The proposed design is a very compact, 1 m class telescope compatible with modern CCD and CMOS detectors (EFL = 15 m). We describe the concept implementation guidelines and the optical performance of the current optical design. The diffraction limited FOV exceeds 1.25 square degrees, and the detector occupies the best 0.25 square degree with 66 devices.

We theoretically and observationally investigate different choices of initial conditions for the primordial mode function that are imposed during an epoch preceding inflation. By deriving predictions for the observables resulting from several alternate quantum vacuum prescriptions we show some choices of vacua are theoretically observationally distinguishable from others. Comparing these predictions to the Planck 2018 observations via a Bayesian analysis shows no significant evidence to favour any of the quantum vacuum prescriptions over the others. In addition we consider frozen initial conditions, representing a white-noise initial state at the big-bang singularity. Under certain assumptions the cosmological concordance model and frozen initial conditions are found to produce identical predictions for the cosmic microwave background anisotropies. Frozen initial conditions may thus provide an alternative theoretic paradigm to explain observations that were previously understood in terms of the inflation of a quantum vacuum.

We study the homologous collapse of stellar nuclear core, the virial theorem for hadron collisional relaxations, and photon productions from hadron collisions. We thus show the gravo-thermal dynamical process that transforms gravitational energy to photon energy. The process is energetically and entropically favourable. The total baryon number conservation, Euler equation for energy-momentum conservation and Poisson's equation for gravitational potential are adopted to describe homologous core collapses. The virial theorem determines the hadron collision energy gain from gravitational potential. The hadronic photon production rate determines the photon energy density. The time scales of macroscopic and microscopic processes are studied to verify approximations. As a result, we show the formation of opaque photon spheres, whose total energy, size, temperature and number density, accounting for the main energetic features of Gamma-Ray Burst progenitors. We obtain the intrinsic correlations of these quantities. They depend only on the averaged thermal index of the stellar core, and it is possible to confront them with observational data.

Many stars are in binaries or higher-order multiple stellar systems. Although in recent years a large number of binaries have been proven to host exoplanets, how planet formation proceeds in multiple stellar systems has not been studied much yet from the theoretical standpoint. In this paper we focus on the evolution of the dust grains in planet-forming discs in binaries. We take into account the dynamics of gas and dust in discs around each component of a binary system under the hypothesis that the evolution of the circumprimary and the circumsecondary discs is independent. It is known from previous studies that the secular evolution of the gas in binary discs is hastened due to the tidal interactions with their hosting stars. Here we prove that binarity affects dust dynamics too, possibly in a more dramatic way than the gas. In particular, the presence of a stellar companion significantly reduces the amount of solids retained in binary discs because of a faster, more efficient radial drift, ultimately shortening their lifetime. We prove that how rapidly discs disperse depends both on the binary separation, with discs in wider binaries living longer, and on the disc viscosity. Although the less-viscous discs lose high amounts of solids in the earliest stages of their evolution, they are dissipated slowly, while those with higher viscosities show an opposite behaviour. The faster radial migration of dust in binary discs has a striking impact on planet formation, which seems to be inhibited in this hostile environment, unless other disc substructures halt radial drift further in. We conclude that if planetesimal formation were viable in binary discs, this process would take place on very short time scales.

3C 84 (NGC 1275) is the bright radio core of the Perseus Cluster. Even in the absence of strong relativistic effects, the source has been detected at Gamma-rays up to TeV energies. Despite its intensive study, the physical processes responsible for the high-energy emission in the source remain unanswered. We present a detailed kinematics study of the source and its connection to Gamma-ray emission. The sub-parsec scale radio structure is dominated by slow-moving features in both the eastern and western lanes of the jet. The jet appears to have accelerated to its maximum speed within less than 125 000 gravitational radii. The fastest reliably detected speed in the jet was ~0.9 c. This leads to a minimum Lorentz factor of ~1.35. Our analysis suggests the presence of multiple high-energy sites in the source. If Gamma-rays are associated with kinematic changes in the jet, they are being produced in both eastern and western lanes in the jet. Three Gamma-ray flares are contemporaneous with epochs where the slowly moving emission region splits into two sub-regions. We estimate the significance of these events being associated as ~2-3 sigma. We tested our results against theoretical predictions for magnetic reconnection-induced mini-jets and turbulence and find them compatible.

Aims. We aim to infer information about the magnetic field in the low solar corona from coronal rain clumps using high-resolution spectropolarimetric observations in the Ca ii 8542 A line obtained with the Swedish 1-m Solar Telescope. Methods. The weak-field approximation (WFA) provides a simple tool to obtain the line-of-sight component of the magnetic field from spectropolarimetric observations. We adapted a method developed in a previous paper in order to assess the different conditions that must be satisfied in order to properly use the WFA for the data at hand. We also made use of velocity measurements in order to estimate the plane-of-the-sky magnetic field component, so that the magnetic field vector could be inferred. Results. We have inferred the magnetic field vector from a data set totalling 100 spectral scans in the Ca ii 8542 A line, containing an off-limb view of the lower portion of catastrophically cooled coronal loops in an active region. Our results, albeit limited by the cadence and signal-to-noise ratio of the data, suggest that magnetic field strengths of hundreds of Gauss, even reaching up to 1000 G, are omnipresent at coronal heights below 9 Mm from the visible limb. Our results are also compatible with the presence of larger magnetic field values such as those reported by previous works. However, for large magnetic fields, the Doppler width from coronal rain is not that much larger than the Zeeman width, thwarting the application of the WFA. Furthermore, we have determined the temperature, T , and microturbulent velocity, $\xi$, of coronal rain clumps and off-limb spicules present in the same data set, and we have found that the former ones have narrower T and $\xi$ distributions, their average temperature is similar, and coronal rain has microturbulent velocities smaller than those of spicules

Over the past decades, the kinematics of galaxies in the local Universe and at intermediate redshift (i.e., z~1-3) have been characterized in great detail, but only a handful of galaxies at high redshift (z>4) have been examined in such a way. The Atacama Large Millimeter/submillimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE) survey observed 118 star-forming main sequence galaxies at z=4.4-5.9 in [CII]158um emission, increasing the number of such observations by nearly an order of magnitude. To characterize the morpho-kinematics of this sample, we apply a well-tested tilted ring model fitting code (3DBarolo), a quantitative morphological classification (Gini-M20), and a set of disk identification criteria to the ALPINE data. By exploring the G-M20 of z>4 rest-frame FIR and [CII] data for the first time, we find that our 1"~6kpc resolution is insufficient to separate galaxy types based solely on these data. Of the 75 [CII]-detected ALPINE galaxies, 29 are detected at high enough significance and with sufficient spatial resolution to allow for tilted ring model fitting and the derivation of morpho-kinematic parameters. By combining these results with disk identification criteria, we are able to robustly classify 14 of the 29 fit sources (six rotators, five mergers, and three dispersion-dominated systems), with the remaining sources showing complex behaviour. We then compare the rotation curves and dynamical mass profiles of the six ALPINE rotators to the two previously detected z~4-6 unlensed main sequence rotators, finding high rotational velocities (~50-250km/s) and a range of rotation curve shapes.

We use the Gemini NIFS instrument to map the H$_2 2.1218\mu$m and Br$\gamma$ flux distributions in the inner 0.04-2 kpc of a sample of 36 nearby active galaxies ($0.001\lesssim z\lesssim0.056$) at spatial resolutions from 4 to 250 pc. We find extended emission in 34 galaxies. In $\sim$55% of them, the emission in both lines is most extended along the galaxy major axis, while in the other 45% the extent follows a distinct orientation. The emission of H$_2$ is less concentrated than that of Br$\gamma$, presenting a radius that contains half of the flux 60% greater, on average. The H$_2$ emission is driven by thermal processes - X-ray heating and shocks - at most locations for all galaxies, where $0.4<H_2/Br\gamma<6$. For regions where H$_2$/Br$\gamma>6$ (seen in 40% of the galaxies), shocks are the main H$_2$ excitation mechanism, while in regions with H$_2$/Br$\gamma<0.4$ (25% of the sample) the H$_2$ emission is produced by fluorescence. The only difference we found between type 1 and type 2 AGN was in the nuclear emission-line equivalent widths, that are smaller in type 1 than in type 2 due to a larger contribution to the continuum from the hot dusty torus in the former. The gas masses in the inner 125 pc radius are in the range $10^1-10^4$ M$_\odot$ for the hot H$_2$ and $10^3-10^6$ M$_\odot$ for the ionised gas and would be enough to power the AGN in our sample for $10^5-10^8$ yr at their current accretion rates.

It is thought that magnetic fields must be present in the interiors of stars to resolve certain discrepancies between theory and observation (e.g. angular momentum transport), but such fields are difficult to detect and characterise. Asteroseismology is a powerful technique for inferring the internal structures of stars by measuring their oscillation frequencies, and succeeds particularly with evolved stars, owing to their mixed modes, which are sensitive to the deep interior. The goal of this work is to present a phenomenological study of the combined effects of rotation and magnetism in evolved stars, where both are assumed weak enough that first-order perturbation theory applies, and we focus on the regime where Coriolis and Lorentz forces are comparable. Axisymmetric "twisted-torus" field configurations are used, which are confined to the core and allowed to be misaligned with respect to the rotation axis. Factors such as the field radius, topology and obliquity are examined. We observe that fields with finer-scale radial structure and/or smaller radial extent produce smaller contributions to the frequency shift. The interplay of rotation and magnetism is shown to be complex: we demonstrate that it is possible for nearly symmetric multiplets of apparently low multiplicity to arise even under a substantial field, which might falsely appear to rule out its presence. Our results suggest that proper modelling of rotation and magnetism, in a simultaneous fashion, may be required to draw robust conclusions about the existence/non-existence of a core magnetic field in any given object.

Recent high precision meteoritic data infers that Mars finished its accretion rapidly within 10 Myr of the beginning of the Solar system and had an accretion zone that did not entirely overlap with the Earth's. Here we present a detailed study of the accretion zone of planetary embryos from high resolution simulations of planetesimals in a disc. We found that all simulations with Jupiter and Saturn on their current eccentric orbits (EJS) result in a similar accretion zone between fast-forming Mars and Earth region embryos. Assuming more circular orbits for Jupiter and Saturn (CJS), on the other hand, has a significantly higher chance of forming Mars with an accretion zone not entirely dominated by Earth and Venus region embryos, however CJS in general forms Mars slower than in EJS. By further quantifying the degree of overlap between accretion zones of embryos in different regions with the average overlap coefficient (OVL), we found that the OVL of CJS shows a better match with the OVL from a chondritic isotopic mixing model of Earth and Mars, which indicates that the giant planets are likely to have resided on more circular orbits than today during gas disc dissipation, matching their suggested pre-instability orbits. More samples, including those from Mercury and Venus, could potentially confirm this hypothesis.

Inside-Out Planet Formation (IOPF) is a theory addressing the origin of Systems of Tightly-Packed Inner Planets (STIPs) via {\it in situ} formation and growth of the planets. It predicts that a pebble ring is established at the pressure maximum associated with the dead zone inner boundary (DZIB) with an inner disk magnetorotational instability (MRI)-active region. Using direct $N$-body simulations, we study the collisional evolution of planetesimals formed from such a pebble ring, in particular examining whether a single dominant planet emerges. We consider a variety of models, including some in which the planetesimals are continuing to grow via pebble accretion. We find that the planetesimal ring undergoes oligarchic evolution, and typically turns into 2 or 3 surviving oligarchs on nearly coplanar and circular orbits, independent of the explored initial conditions or form of pebble accretion. The most massive oligarchs typically consist of about $70\%$ of the total mass, with the building-up process typically finishing within $\sim 10^5$ years. However, a relatively massive secondary planet always remains with $\sim30-65\%$ of the mass of the primary. Such secondary planets have properties that are inconsistent with the observed properties of the innermost pairs of planets in STIPs. Thus, for IOPF to be a viable theory for STIP formation, it needs to be shown how oligarchic growth of a relatively massive secondary from the initial pebble ring can be avoided. We discuss some potential additional physical processes that should be included in the modeling and explored as next steps.

We revisit the constraints on evaporating primordial black holes (PBHs) from the isotropic X-ray and soft gamma-ray background in the mass range $10^{16}-10^{18}$ g. We find that they are stronger than usually inferred due to two neglected effects: i) The contribution of the annihilation radiation due to positrons emitted in the evaporation process. ii) The high-latitude, Galactic contribution to the measured isotropic flux. We study the dependence of the bounds from the datasets used, the positron annihilation conditions, and the inclusion of the astrophysical background. We derive competitive bounds excluding non-spinning PBH with monochromatic mass function as the totality of dark matter for masses below about 1.6$\times 10^{17}\,$g. We also show that the inclusion of spin and/or an extended, log-normal mass function lead to tighter bounds. Our study suggests that the isotropic flux is an extremely promising target for future missions in improving the sensitivity to PBHs as candidates for dark matter.

Astrometric Science and Technology Roadmap for Astrophysics (ASTRA) is a bilateral cooperation between China and Italy with the goal of consolidating astrometric measurement concepts and technologies. In particular, the objectives include critical analysis of the Gaia methodology and performance, as well as principle demonstration experiments aimed at future innovative astrometric applications requiring high precision over large angular separations (one to 180 degrees). Such measurement technologies will be the building blocks for future instrumentation focused on the "great questions" of modern cosmology, like General Relativity validity (including Dark Matter and Dark Energy behavior), formation and evolution of structure like proto-galaxies, and planetary systems formation in bio compatibles environments. We describe three principle demonstration tests designed to address some of the potential showstoppers for high astrometric precision experiments. The three tests are focused on the key concepts of multiple fields telescopes, astrometric metrology and very fine sub-pixel precision (goal: <1/2000 pixel) in white light.

HIP 41378 f is a temperate $9.2\pm0.1 R_{\oplus}$ planet with period of 542.08 days and an extremely low density of $0.09\pm0.02$ g cm$^{-3}$. It transits the bright star HIP 41378 (V=8.93), making it an exciting target for atmospheric characterization including transmission spectroscopy. HIP 41378 was monitored photometrically between the dates of 2019 November 19 and November 28. We detected a transit of HIP 41378 f with NGTS, just the third transit ever detected for this planet, which confirms the orbital period. This is also the first ground-based detection of a transit of HIP 41378 f. Additional ground-based photometry was also obtained and used to constrain the time of the transit. The transit was measured to occur 1.50 hours earlier than predicted. We use an analytic transit timing variation (TTV) model to show the observed TTV can be explained by interactions between HIP 41378 e and HIP 41378 f. Using our TTV model, we predict the epochs of future transits of HIP 41378 f, with derived transit centres of T$_{C,4} = 2459355.087^{+0.031}_{-0.022}$ (May 2021) and T$_{C,5} = 2459897.078^{+0.114}_{-0.060}$ (Nov 2022).

S stars are late-type giants that are transition objects between M-type stars and carbon stars on the asymptotic giant branch (AGB). They are classified into two types: intrinsic or extrinsic, based on the presence or absence of technetium (Tc). The Tc-rich or intrinsic S stars are thermally-pulsing (TP-)AGB stars internally producing s-process elements (including Tc) which are brought to their surface via the third dredge-up (TDU). Tc-poor or extrinsic S stars gained their s-process overabundances via accretion of s-process-rich material from an AGB companion which has since turned into a dim white dwarf. Our goal is to investigate the evolutionary status of Tc-rich S stars by locating them in a Hertzsprung-Russell (HR) diagram using the results of Gaia early Data Release 3 (EDR3). We combine the current sample of 13 Tc-rich stars with our previous studies of 10 Tc-rich stars to determine the observational onset of the TDU in the metallicity range [-0.7; 0]. We also compare our abundance determinations with dedicated AGB nucleosynthesis predictions. The stellar parameters are derived using an iterative tool which combines HERMES high-resolution spectra, accurate Gaia EDR3 parallaxes, stellar evolution models and tailored MARCS model atmospheres for S-type stars. Using these stellar parameters we determine the heavy-element abundances by line synthesis. In the HR diagram, the intrinsic S stars are located at higher luminosities than the predicted onset of the TDU. These findings are consistent with Tc-rich S stars being genuinely TP-AGB stars. The comparison of the derived s-process abundance profiles of our intrinsic S stars with the nucleosynthesis predictions provide an overall good agreement. Stars with highest [s/Fe] tend to have the highest C/O ratios.

The radio sky at frequencies below $\sim10$ MHz is still largely unknown, this remains the last unexplored part of the electromagnetic spectrum in astronomy. The upcoming space experiments aiming at such low frequencies (ultra-long wavelength or ultra-low frequency) would benefit from reasonable expectations of the sky brightness distribution at relevant frequencies. In this work, we develop a radio sky model that is valid down to $\sim1$ MHz. In addition to the discrete HII objects, we take into account the free-free absorption by thermal electrons in the Milky Way's warm ionized medium (WIM). This absorption effect becomes obvious at $\lesssim10$ MHz, and could make the global radio spectrum turn over at $\sim3$ MHz. Our sky map shows unique features at the ultra-long wavelengths, including a darker Galactic plane in contrast to the sky at higher frequencies, and the huge shadows of the spiral arms on the sky map. It would be a useful guidance for designing the future ultra-long wavelength observations. Our model could be downloaded at https://doi.org/10.5281/zenodo.4454153.

We report the first survey of molecular emission from cometary volatiles using standalone Atacama Compact Array (ACA) observations of the Atacama Large Millimeter/Submillimeter Array (ALMA) toward comet C/2015 ER61 (PanSTARRS) carried out on UT 2017 April 11 and 15, shortly after its April 4 outburst. These measurements of HCN, CS, CH$_3$OH, H$_2$CO, and HNC (along with continuum emission from dust) probed the inner coma of C/2015 ER61, revealing asymmetric outgassing and discerning parent from daughter/distributed source species. This work presents spectrally integrated flux maps, autocorrelation spectra, production rates, and parent scale lengths for each molecule, and a stringent upper limit for CO. HCN is consistent with direct nucleus release in C/2015 ER61, whereas CS, H$_2$CO, HNC, and potentially CH$_3$OH are associated with distributed sources in the coma. Adopting a Haser model, parent scale lengths determined for H$_2$CO (L$_p$ $\sim$ 2200 km) and HNC (L$_p$ $\sim$ 3300 km) are consistent with previous work in comets, whereas significant extended source production (L$_p$ $\sim$ 2000 km) is indicated for CS, suggesting production from an unknown parent in the coma. The continuum presents a point-source distribution, with a flux density implying an excessively large nucleus, inconsistent with other estimates of the nucleus size. It is best explained by the thermal emission of slowly-moving outburst ejectas, with total mass 5--8 $\times$ 10$^{10}$ kg. These results demonstrate the power of the ACA for revealing the abundances, spatial distributions, and locations of molecular production for volatiles in moderately bright comets such as C/2015 ER61.

The modern accuracy of measurements allows the residual/peculiar (Galactocentric) velocity of the supermassive black hole (SMBH) in our Galaxy, Sgr A*, on the order of several kilometers per second. We integrate possible orbits of the SMBH along with the surrounding nuclear star cluster (NSC) for a barred model of the Galaxy using modern constraints on the components of the SMBH Galactocentric velocity. Is is shown that the range of oscillations of the SMBH + NSC in a regular Galactic field in the plane of the Galaxy allowed by these constraints strongly depends on the set of central components of the Galactic potential. If the central components are represented only by a bulge/bar, for a point estimate of the SMBH Galactocentric velocity, the oscillation amplitude does not exceed 7 pc in the case that a classical bulge is present and reaches 25 pc if there is no bulge; with SMBH velocity components within the $2\sigma$ significance level, the amplitude can reach 15 and 50 pc, respectively. However, when taking into account the nuclear stellar disk (NSD), even in the absence of a bulge, the oscillation amplitude is only 5 pc for the point estimate of the SMBH velocity, and 10 pc for the $2\sigma$ significance level. Thus, the possible oscillations of the SMBH + NSC complex from the confirmed components of the Galaxy's potential are mostly limited by the NSD, and even taking into account the uncertainty of the mass of the latter, the oscillation amplitude can hardly exceed $13\,\text{pc}=6'$.

Accurate mass-loss rates and terminal velocities from massive stars winds are essential to obtain synthetic spectra from radiative transfer calculations and to determine the evolutionary path of massive stars. From a theoretical point of view, analytical expressions for the wind parameters and velocity profile would have many advantages over numerical calculations that solve the complex non-linear set of hydrodynamic equations. In a previous work, we obtained an analytical description for the fast wind regime. Now, we propose an approximate expression for the line-force in terms of new parameters and obtain a velocity profile closed-form solution (in terms of the Lambert $W$ function) for the $\delta$-slow regime. Using this analytical velocity profile, we were able to obtain the mass-loss rates based on the m-CAK theory. Moreover, we established a relation between this new set of line-force parameters with the known stellar and m-CAK line-force parameters. To this purpose, we calculated a grid of numerical hydrodynamical models and performed a multivariate multiple regression. The numerical and our descriptions lead to good agreement between their values.

Because many of our X-ray telescopes are optimized towards observing faint sources, observations of bright sources like X-ray binaries in outburst are often affected by instrumental biases. These effects include dead time and photon pile-up, which can dramatically change the statistical inference of physical parameters from these observations. While dead time is difficult to take into account in a statistically consistent manner, simulating dead time-affected data is often straightforward. This structure makes the issue of inferring physical properties from dead time-affected observations fall into a class of problems common across many scientific disciplines. There is a growing number of methods to address them under the names of Approximate Bayesian Computation (ABC) or Simulation-Based Inference (SBI), aided by new developments in density estimation and statistical machine learning. In this paper, we introduce SBI as a principled way to infer variability properties from dead time-affected light curves. We use Sequential Neural Posterior Estimation to estimate the posterior probability for variability properties. We show that this method can recover variability parameters on simulated data even when dead time is variable, and present results of an application of this approach to NuSTAR observations of the galactic black hole X-ray binary GRS 1915+105.

We developed a parametric Galactic model toward the Galactic bulge by fitting to spatial distributions of the Gaia DR2 disk velocity, VVV proper motion, BRAVA radial velocity, OGLE-III red clump star count, and OGLE-IV star count and microlens rate, optimized for use in microlensing studies. We include the asymmetric drift of Galactic disk stars and the dependence of velocity dispersion on Galactic location in kinematic model, which has been ignored in most previous models used for microlensing studies. We show that our model predicts a microlensing parameter distribution that is significantly different from the one with a typically used model in the previous studies. Through our modeling, we estimate various fundamental model parameters for our Galaxy, including the initial mass function (IMF) in the inner Galaxy. Combined constraints from star counts and the microlensing event timescale distribution from the OGLE-IV survey, in addition to a prior on the bulge stellar mass, enable us to successfully measure IMF slopes using a broken power law form over a broad mass range, $\alpha_{\rm bd} = 0.22^{+0.20}_{-0.55}$ for $M < 0.08 \, M_{\odot}$, $\alpha_{\rm ms} = 1.16^{+0.08}_{-0.15}$ for $0.08 \, M_{\odot} \leq M < M_{\rm br}$, and $\alpha_{\rm hm} = 2.32^{+0.14}_{-0.10}$ for $M \geq M_{\rm br}$, as well as a break mass at $M_{\rm br} = 0.90^{+0.05}_{-0.14} \, M_{\odot}$. This is significantly different from the Kroupa IMF for local stars, but similar to the Zoccali IMF measured from a bulge luminosity function. We also estimate the dark matter mass fraction in the bulge region of $28 \pm 7$\% which could be larger than a previous estimate. Because our model is purely parametric, it can be universally applied using the parameters provided in this paper.

There is an accepted approach to calculation of the neutrino flavor density-matrix in the halo of a supernova, in which neutrino amplitudes, not cross-sections, need to be followed carefully in the region above the region of frequent scatterings. The same reasoning and techniques, applied to the evolution of neutrino flavors and energy distributions in the early universe in the era of neutrino decoupling, leads to radical changes in the predictions of the effects of the neutrino-neutrino interaction. Predictions for the production of sterile neutrinos, should they exist, will also be changed.

Using the fact that we only observe those modes which exit the Hubble horizon during inflation, one can calculate the entanglement entropy of such long-wavelength perturbations by tracing out unobservable sub-Hubble fluctuations they are coupled to. On requiring that this perturbative entanglement entropy, which increases with time, obey the covariant entropy bound for an accelerating background, we find an upper bound on the duration of inflation. This presents a new perspective on the (meta-)stability of de Sitter spacetime and an associated lifetime for it.

In this work, the behaviour of magneto-hydrodynamic waves in optically thin plasmas considering dissipative processes, thermal and magnetic diffusion, a given ionization and the heating and cooling functions are investigated for several particular cases. A numerical eigenvalues analysis of the dimensionless secular equations according to various cases is performed for the entire set of MHD equations.

Neutrino interaction with photons induced from a strong magnetic field approximation in the electroweak theory is presented. Contrary to the weak filed approximation, the interaction in the magnetized plasma has non-negligible strength and a topological form expressed by the Chern-Simons form of the neutrino current and the electromagnetic vector potential. The derivation of the interaction Lagrangian and its properties are presented.

In this paper, we discuss interesting potential implications for the supersymmetric (SUSY) universe in light of cosmological problems on (1) the number of the satellite galaxies of the Milky Way (missing satellite problem) and (2) a value of the matter density fluctuation at the scale around 8$h^{-1}$Mpc ($S_{8}$ tension). The implications are extracted by assuming that the gravitino of a particular mass can be of help to alleviate the cosmological tension. We consider two gravitino mass regimes vastly separated, that is, $m_{3/2}\simeq100{\rm eV}$ and $m_{3/2}\simeq100{\rm GeV}$. We discuss non-trivial features of each supersymmetric universe associated with a specific gravitino mass by projecting potential resolutions of the cosmological problems on each of associated SUSY models.

Electric field measurements of the Time Domain Sampler (TDS) receiver, a part of the Radio and Plasma Waves (RPW) instrument on-board Solar Orbiter, often exhibit very intense broadband wave emissions at frequencies below 20~kHz in the spacecraft frame. During the first year of the mission, the RPW/TDS instrument has been operating from the first perihelion in mid-June 2020 and through the first flyby of Venus in late December 2020. In this paper, a year-long study of electrostatic fluctuations observed in the solar wind at an interval of heliocentric distances from 0.5 to 1~AU is shown. The on-board processed properties of waveform snapshots that are continuously collected allow mapping plasma waves at frequencies between 200~Hz and 20~kHz. For a detailed spectral and polarization analysis, the triggered waveform snapshots and a Doppler-shifted solution of the dispersion relation for wave mode identification were used. The occurrence rate of low-frequency waves peaks around perihelion at distances of 0.5~AU and decreases with increasing distances, with only a few waves detected per day at 0.9~AU. A more detailed analysis of more than ten thousand triggered waveform snapshots shows the median wave frequency at about 2.3~kHz and wave amplitude about 1.1~mV/m. The relative phase distribution between two components of E-field projected in the Y-Z Spacecraft Reference Frame (SRF) plane shows a mostly linear wave polarization. Electric field fluctuations are closely aligned with the directions of the ambient field lines. The observed waves are interpreted as the strongly Doppler-shifted electrostatic ion-acoustic mode generated by the resonant interaction with ion beams or by the current-driven instability.

We consider a generic dispersive massive gravity theory and numerically study its resulting modified energy and strain spectra of gravitational waves (GWs) sourced by (i) fully developed turbulence during the electroweak phase transition (EWPT); and (ii) forced hydromagnetic turbulence during the QCD phase transition (QCDPT). The GW spectra are then computed in both the spatial and temporal Fourier domains. We find, from the spatial spectra, that the slope modifications are independent of the eddy size at QCDPT, and, from the temporal spectra, that the modifications are pronounced in the nHz-10nHz range -- the sensitivity range of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) -- for a graviton mass $m_{\rm g}$ in the range $2.2\times10^{-23}{\rm eV}/c^2\lesssim m_{\rm g}\lesssim7.4\times10^{-22}{\rm eV}/c^2$.

The Fermilab E989 experiment has recently reported their result of muon anomalous magnetic moment ($g$-2). Combined with the E821 result, the discrepancy with the Standard Model (SM) reaches $4.2\sigma$, which may indicate the light new physics related with the electroweak interactions. On the other hand, the observed Galactic Center GeV Excess (GCE) and anti-proton excess can also be explained by the light weakly interacting massive particle dark matter. In this work, we attempt to pin down a common origin of these anomalies in the Next-to-Minimal Supersymmetric Model. By considering various constraints, we find that the eletroweakinos and sleptons have to be lighter than about 1 TeV, and the geometric mean of their masses need to be less than about 375 GeV to interpret the muon $g-2$ within $2\sigma$ range. In order to accommodate both the muon $g-2$ and GCE, the bino-like neutralino DM is needed and has to resonantly annihilate through the $Z$ boson or Higgs boson to give the correct relic density. Besides, the DM annihilating cross section for GCE can be achieved by a singlet-like Higgs boson via $s$-channel. If the anti-proton excess is explained together, only the Higgs funnel is feasible. We point out that the favored parameter space to explain all these anomalies can be probed by the future direct detection experiments.

Magnetic monopoles have been a subject of study for more than a century since the first ideas by A. Vaschy and P. Curie, circa 1890. In 1974, Y. Nambu proposed a model for magnetic monopoles exploring a parallelism between the broken symmetry Higgs and the superconductivity Ginzburg-Landau theories in order to describe the pions quark-antiquark confinement states. There, Nambu describes an energetic string where its end points behave like two magnetic monopoles with opposite magnetic charges -- quark and antiquark. Consequently, not only the interaction among monopole and antimonopole, mediated by a massive vector boson (Yukawa potential), but also the energetic string (linear potential) contributes to the effective interaction potential. We propose here a monopole-antimonopole non confining attractive interaction of the Nambu-type, and then investigate the formation of bound states, the monopolium. Some necessary conditions for the existence of bound states to be fulfilled by the proposed Nambu-type potential, Kato weakness, Set\^o and Bargmann conditions, are verified. In the following, ground state energies are estimated for a variety of monopolium reduced mass, from $10^2$MeV to $10^2$TeV, and Compton interaction lengths, from $10^{-2}$am to $10^{-1}$pm, where discussion about non relativistic and relativistic limits validation is carried out.

The Lepton Portal Dark Matter model, in which dark matter states only coupling to the charged leptons, can explain the excess of the muon anomalous magnetic moment measured by the Muon $g - 2$ experiment. In this paper, we demonstrate that real, charge-neutral scalar dark matter with a large number of internal degrees of freedom and a mass approximately degenerate with the charged fermionic mediator state can accommodate the $(g - 2)_\mu$ excess. The model remains consistent with the dark matter relic abundance, direct detection, and indirect detection constraints. The dark matter and its charged fermion partner masses are constrained to be below around 200 GeV. The high-luminosity LHC and future lepton colliders, as well as indirect searches at CTA and GAMMA-400, can test this scenario.