Papers for Wednesday, Jul 07 2021

Type IIn supernovae (SNe IIn) are a relatively infrequently observed subclass of SNe whose photometric and spectroscopic properties are varied. A common thread among SNe IIn are the complex multiple-component hydrogen Balmer lines. Owing to the heterogeneity of SNe IIn, online databases contain some outdated, erroneous, or even contradictory classifications. SN IIn classification is further complicated by SN impostors and contamination from underlying HII regions. We have compiled a catalogue of systematically classified nearby (redshift z < 0.02) SNe IIn using the Open Supernova Catalogue (OSC). We present spectral classifications for 115 objects previously classified as SNe IIn. Our classification is based upon results obtained by fitting multiple Gaussians to the H-alpha profiles. We compare classifications reported by the OSC and Transient Name Server (TNS) along with the best matched templates from SNID. We find that 28 objects have been misclassified as SNe IIn. TNS and OSC can be unreliable; they disagree on the classifications of 51 of the objects and contain a number of erroneous classifications. Furthermore, OSC and TNS hold misclassifications for 34 and twelve (respectively) of the transients we classify as SNe IIn. In total, we classify 87 SNe IIn. We highlight the importance of ensuring that online databases remain up to date when new or even contemporaneous data become available. Our work shows the great range of spectral properties and features that SNe IIn exhibit, which may be linked to multiple progenitor channels and environment diversity. We set out a classification sche me for SNe IIn based on the H-alpha profile which is not greatly affected by the inhomogeneity of SNe IIn.

(Aims) In this work, we predict the far-infrared polarisation signal emitted by non-spherical dust grains in nearby galaxies. We determine the angular resolution and sensitivity required to study the magnetic field configuration in these galaxies. (Methods) We post-process a set of Milky Way like galaxies from the Auriga project, assuming a dust mix consisting of spheroidal dust grains that are partially aligned with the model magnetic field. We constrain our dust model using Planck 353 GHz observations of the Milky Way. This model is then extrapolated to shorter wavelengths that cover the peak of interstellar dust emission and to observations of arbitrarily oriented nearby Milky Way like galaxies. (Results) Assuming an intrinsic linear polarisation fraction that does not vary significantly with wavelength for wavelengths longer than 50 micron, we predict a linear polarisation fraction with a maximum of $10-15\%$ and a median value of $\approx{}7\%$ for face-on galaxies and $\approx{}3\%$ for edge-on galaxies. The polarisation fraction anti-correlates with the line of sight density and with the angular dispersion function which expresses the large scale order of the magnetic field perpendicular to the line of sight. The maximum linear polarisation fraction agrees well with the intrinsic properties of the dust model. The true magnetic field orientation can be traced along low density lines of sight when it is coherent along the line of sight. These results also hold for nearby galaxies, where a coherent magnetic field structure is recovered over a range of different broad bands. (Conclusions) Polarised emission from non-spherical dust grains accurately traces the large scale structure of the galactic magnetic field in Milky Way like galaxies, with expected maximum linear polarisation fractions of $10-15\%$. To resolve this maximum, a spatial resolution of at least 1 kpc is required.

Primordial black holes (PBHs) are dark matter candidates that span broad mass ranges from $10^{-17}$ $M_\odot$ to $\sim 100$ $M_\odot$. We show that the stochastic gravitational wave background can be a powerful window for the detection of sub-solar mass PBHs and shed light on their formation channel via third-generation gravitational wave detectors such as Cosmic Explorer and the Einstein Telescope. By using the mass distribution of the compact objects and the redshift evolution of the merger rates, we can distinguish astrophysical sources from PBHs and will be able to constrain the fraction of sub-solar mass PBHs $\leq 1$ $M_\odot$ in the form of dark matter $f_{PBH}\leq 1\%$ at $68\%$ C.L. even for a pessimistic value of the suppression factor ($f_{sup} \sim 10^{-3}$). For $f_{sup} \sim 1$, the constraints on $f_{PBH}$ will be less than $0.001\%$. Furthermore, we will be able to measure the redshift evolution of the PBH merger rate with about $1\%$ accuracy, making it possible to uniquely distinguish between the Poisson and clustered PBH scenarios.

We combine the Santa-Cruz Semi-Analytic Model (SAM) for galaxy formation and evolution with the circumgalactic medium (CGM) model presented in Faerman et al. (2020) to explore the CGM properties of $L^{*}$ galaxies. We use the SAM to generate a sample of galaxies with halo masses similar to the Milky Way (MW) halo, $M_{\rm vir} \approx 10^{12}~{\rm M_{sun}}$, and find that the CGM mass and mean metallicity in the sample are correlated. We use the CGM masses and metallicities of the SAM galaxies as inputs for the FSM20 model, and vary the amount of non-thermal support. The density profiles in our models can be approximated by power-law functions with slopes in the range of $0.75 < a_n < 1.25$, with higher non-thermal pressure resulting in flatter distributions. We explore how the gas pressure, dispersion measure, OVI-OVIII column densities, and cooling rates behave with the gas distribution and total mass. We show that for CGM masses below $\sim 3 \times 10^{10}~{\rm M_{sun}}$, photoionization has a significant effect on the column densities of OVI and OVIII. The combination of different MW CGM observations favors models with similar fractions in thermal pressure, magnetic fields/cosmic rays, and turbulent support, and with $M_{\rm gas} \sim 3-10 \times 10^{10}~{\rm M_{sun}}$. The MW OVI column requires $t_{\rm cool}/t_{\rm dyn} \sim 4$, independent of the gas distribution. The AGN jet-driven heating rates in the SAM are enough to offset the CGM cooling, although exact balance is not required in star-forming galaxies. We provide predictions for the columns densities of additional metal ions - NV, NeVIII, and MgX.

The Humphreys-Davidson (HD) limit sets the boundary between evolutionary channels of massive stars that either end their lives as red supergiants (RSGs) or as the hotter blue supergiants (BSGs) and Wolf-Rayet stars. Mixing in the envelopes of massive stars close to their Eddington limit is crucial for investigating the upper luminosity limit of the coolest supergiants. We study the effects of excess mixing in superadiabatic layers that are dominated by radiation pressure, and we critically investigate the effects of mixing and mass loss on the evolution of RSGs with log (Teff/K) < 3.68 - as a function of metallicity. Using MESA, we produce grids of massive star models at three metallicities: Galactic (Zsol), LMC (1/2 Zsol) and SMC (1/5 Zsol), with both high and low amounts of overshooting to study the upper luminosity limit of RSGs. We systematically study the effects of excess mixing in the superadiabatic layers of post-main sequence massive stars, overshooting above the hydrogen core and yellow supergiant (YSG) mass-loss rates on the fraction of core helium burning time spent as a RSG. We find that the excess mixing in the superadiabatic layers is stronger at lower metallicities, as it depends on the opacities in the hydrogen bump at log (Teff/K) ~ 4, which become more pronounced at lower metallicity. This shifts the cutoff luminosities to lower values at lower metallicities, thus balancing the first-order effect of mass loss. The opposing effects of mass loss and excess envelope mixing during post-main sequence evolution of stars with higher overshooting potentially results in a metallicity-independent upper luminosity limit.

Supernova cosmology surveys are traditionally time consuming, especially for the critical spectroscopic data. However, a single spectrum at maximum light may provide accurate distance estimation if recent developments hold. This could open up a new type of supernova cosmology survey, with a useful interaction between the spectra and a focus on specific redshifts. We optimize the redshift selection and show that this condensed survey could efficiently deliver highly accurate dark energy constraints.

Using the global 21-cm signal measurement by the EDGES collaboration, we derive constraints on the fraction of the dark matter that is in the form of primordial black holes (PBHs) with masses in the range $10^{15}$-$10^{17}\,$g. Improving upon previous analyses, we consider the effect of the X-ray heating of the intergalactic medium on these constraints, and also use the full shape of the 21-cm absorption feature in our inference. In order to account for the anomalously deep absorption amplitude, we also consider an excess radio background motivated by LWA1 and ARCADE2 observations. Because the heating rate induced by PBH evaporation evolves slowly, the data favour a scenario in which PBH-induced heating is accompanied by X-ray heating. Also, for the same reason, using the full measurement across the EDGES observation band yields much stronger constraints on PBHs than just the redshift of absorption. We find that 21-cm observations exclude $f_{\mathrm{PBH}} \gtrsim 10^{-9.7}$ at 95% CL for $M_{\mathrm{PBH}}=10^{15}\,$g. This limit weakens approximately as $M_{\mathrm{PBH}}^4$ towards higher masses, thus providing the strongest constraints on ultralight evaporating PBHs as dark matter over the entire mass range $10^{15}$-$10^{17}\,$g. Under the assumption of a simple spherical gravitational collapse based on the Press-Schechter formalism, we also derive bounds on the curvature power spectrum at extremely small scales ($k\sim 10^{15}\,$Mpc$^{-1}$). This highlights the usefulness of global 21-cm measurements, including non-detections, across wide frequency bands for probing exotic physical processes.

Context. Multidimensional hydrodynamic simulations of convection in stellar interiors are numerically challenging, especially for flows at low Mach numbers. Methods. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M$_\odot$ star model. The results obtained with the low-Mach AUSM$^{+}$-up solver were compared to results when using its non low-Mach variant AUSM$_\mathrm{B}^{+}$-up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from $180 \times 90^2$ to $810 \times 540^2$ cells and the nuclear energy release was boosted by factors of $3 \times 10^3$, $1 \times 10^4$, and $3 \times 10^4$ to study the dependence of the results on these parameters. Results. The boosted energy input results in convection at Mach numbers in the range of $10^{-2}$ to $10^{-3}$. Standard mixing-length theory (MLT) predicts convective velocities of about $1.6 \times 10^{-4}$ if no boosting is applied. Simulations with AUSM$^{+}$-up show a Kolmogorov-like inertial range in the kinetic energy spectrum that extends further toward smaller scales compared with its non low-Mach variant. The kinetic energy dissipation of the AUSM$^{+}$-up solver already converges at a lower resolution compared to AUSM$^{+}_{\mathrm{B}}$ -up. The extracted entrainment rates at the boundaries of the convection zone are well represented by the bulk Richardson entrainment law and the corresponding fitting parameters are in agreement with published results for carbon shell burning.

Dark magnetic spots crossing the stellar disc lead to quasi-periodic brightness variations, which allow us to constrain stellar surface rotation and photometric activity. The current work is the second of this series (Santos et al. 2019; Paper I), where we analyze the Kepler long-cadence data of 132,921 main-sequence F and G stars and late subgiant stars. Rotation-period candidates are obtained by combining wavelet analysis with autocorrelation function. Reliable rotation periods are then selected via a machine learning (ML) algorithm (Breton et al. 2021), automatic selection, and complementary visual inspection. The ML training data set comprises 26,521 main-sequence K and M stars from Paper I. To supplement the training, we analyze in the same way as Paper I, i.e. automatic selection and visual inspection, 34,100 additional stars. We finally provide rotation periods Prot and associated photometric activity proxy Sph for 39,592 targets. Hotter stars are generally faster rotators than cooler stars. For main-sequence G stars, Sph spans a wider range of values with increasing effective temperature, while F stars tend to have smaller Sph values in comparison with cooler stars. Overall for G stars, fast rotators are photometrically more active than slow rotators, with Sph saturating at short periods. The combined outcome of the two papers accounts for average Prot and Sph values for 55,232 main-sequence and subgiant FGKM stars (out of 159,442 targets), with 24,182 new Prot detections in comparison with McQuillan et al. (2014). The upper edge of the Prot distribution is located at longer Prot than found previously.

The Pluto system is an archetype for the multitude of icy dwarf planets and accompanying satellite systems that populate the vast volume of the solar system beyond Neptune. New Horizons' exploration of Pluto and its five moons gave us a glimpse into the range of properties that their kin may host. Furthermore, the surfaces of Pluto and Charon record eons of bombardment by small trans-Neptunian objects, and by treating them as witness plates we can infer a few key properties of the trans-Neptunian population at sizes far below current direct-detection limits. This chapter summarizes what we have learned from the Pluto system about the origins and properties of the trans-Neptunian populations, the processes that have acted upon those members over the age of the solar system, and the processes likely to remain active today. Included in this summary is an inference of the properties of the size distribution of small trans-Neptunian objects and estimates on the fraction of binary systems present at small sizes. Further, this chapter compares the extant properties of the satellites of trans-Neptunian dwarf planets and their implications for the processes of satellite formation and the early evolution of planetesimals in the outer solar system. Finally, this chapter concludes with a discussion of near-term theoretical, observational, and laboratory efforts that can further ground our understanding of the Pluto system and how its properties can guide future exploration of trans-Neptunian space.

In recent years, integral-field spectroscopic surveys have revealed that the presence of kinematically-decoupled stellar components is not a rare phenomenon in nearby galaxies. However, complete statistics are still lacking because they depend on the detection limit of these objects. We aim at investigating the kinematic signatures of two large-scale counter-rotating stellar disks in mock integral-field spectroscopic data to address their detection limit as a function of the galaxy properties and instrumental setup. We build a set of mock data of two large-scale counter-rotating stellar disks as if they were observed with the Multi-Unit Spectroscopic Explorer (MUSE). We account for different photometric, kinematic, and stellar population properties of the two counter-rotating components as a function of galaxy inclination. We extract the stellar kinematics in the wavelength region of the Calcium triplet absorption lines by adopting a Gauss-Hermite (GH) parametrization of the line-of-sight velocity distribution (LOSVD). We confirm that the strongest signature of the presence of two counter-rotating stellar disks is the symmetric double peak in the velocity dispersion map, already known as $2\sigma$ feature. The size, shape, and slope of the 2$\sigma$ peak strongly depend on the velocity separation and relative light contribution of the two counter-rotating stellar disks. When the $2\sigma$ peak is difficult to detect due to the low signal-to-noise of the data, the large-scale structure in the $h_3$ map can be used as a diagnostic for strong and weak counter-rotation. The counter-rotating kinematic signatures become fainter at lower viewing angles as an effect of the smaller projected velocity separation between the two counter-rotating components. We confirm that the observed frequency of $2\sigma$ galaxies represents only a lower limit of the stellar counter-rotation phenomenon.

In this work we show that the cosmological lithium problem and the $H_0$ tension could be eased at the same time by allowing variations in the fundamental constants. We compute the primordial abundances of light elements resulting from Big Bang Nucleosynthesis considering the fine structure constant, the Higgs' vacuum expectation value and Newton's constant as free parameters. Using the observational data for abundances, we set constraints on the variations of the fundamental constants. An interpretation of the results in terms of the number of effective relativistic species gives $N_{\rm eff}=4.04\pm0.12$. If one extrapolates the fit of {\sffamily Planck} considering this $N_{\rm eff}$, the value of the inferred Hubble's constant shifts to $H_0=(71.85\pm0.77)\,\text{km}\,\text{s}^{-1} \text{Mpc}^{-1}$, compatible with current direct determinations.

Evidence of flare induced, large-amplitude, decay-less transverse oscillations is presented. A system of multi-thermal coronal loops as observed with the Atmospheric Imaging Assembly (AIA), exhibit decay-less transverse oscillations after a flare erupts nearby one of the loop footpoints. Measured oscillation periods lie between 4.2 min and 6.9 min wherein the displacement amplitudes range from 0.17 Mm to 1.16 Mm. A motion-magnification technique is employed to detect the pre-flare decay-less oscillations. These oscillations have similar periods (between 3.7 min and 5.0 min) like the previous ones but their amplitudes (0.04 Mm to 0.12 Mm) are found to be significantly smaller. No phase difference is found among oscillating threads of a loop when observed through a particular AIA channel or when their multi-channel signatures are compared. These features suggest that the occurrence of a flare in this case neither changed the nature of these oscillations (decaying vs decay-less) nor the oscillation periods. The only effect the flare has is to increase the oscillation amplitudes.

A radial velocity study by Donati et al. (2016) reported the detection of a close-in giant planet in a 4.93 d orbit around the ~2 Myr old weak-lined T Tauri star V830 Tau. Because of the stringent timescale constraints that a very young host star like V830 Tau would place on hot Jupiter formation models and inward migration mechanisms, independent confirmation of the planet's existence is needed but so far has not been obtained. We present new Chandra X-ray observations of V830 Tau. The Chandra observations in combination with previous XMM-Newton observations reveal strong variable X-ray emission with an X-ray luminosity spanning the range log Lx = 30.10 - 30.87 ergs/s. Chandra High Energy Transmission Grating (HETG) spectra show emission lines formed over a range of plasma temperatures from ~4 MK (Ne IX) to ~16 MK (S XV). At the separation of the reported planet (0.057 au) the X-ray flux is ~10$^{6}$ - 10$^{7}$ times greater than the Sun's X-ray flux at Jupiter. We provide estimates of the X-ray ionization and atmospheric heating rates at the planet's separation and identify areas of uncertainty that will need to be addressed in any future atmospheric models.

The high dynamic range between contaminating foreground emission and the fluctuating 21cm brightness temperature field is one of the most problematic characteristics of 21cm intensity mapping data. While these components would ordinarily have distinctive frequency spectra, making it relatively easy to separate them, instrumental effects and calibration errors further complicate matters by modulating and mixing them together. A popular class of foreground cleaning method are unsupervised techniques related to Principal Component Analysis (PCA), which exploit the different shapes and amplitudes of each component's contribution to the covariance of the data in order to segregate the signals. These methods have been shown to be effective at removing foregrounds, while also unavoidably filtering out some of the 21cm signal too. In this paper we examine, for the first time in the context of 21cm intensity mapping, a generalised method called Kernel PCA, which instead operates on the covariance of non-linear transformations of the data. This allows more flexible functional bases to be constructed, in principle allowing a cleaner separation between foregrounds and the 21cm signal to be found. We show that Kernel PCA is effective when applied to simulated single-dish (autocorrelation) 21cm data under a variety of assumptions about foregrounds models, instrumental effects etc. It presents a different set of behaviours to PCA, e.g. in terms of sensitivity to the data resolution and smoothing scale, outperforming it on intermediate to large scales in most scenarios.

The universe is composed of galaxies that have diverse shapes. Once the structure of a galaxy is determined, it is possible to obtain important information about its formation and evolution. Morphologically classifying galaxies means cataloging them according to their visual appearance and the classification is linked to the physical properties of the galaxy. A morphological classification made through visual inspection is subject to biases introduced by subjective observations made by human volunteers. For this reason, systematic, objective and easily reproducible classification of galaxies has been gaining importance since the astronomer Edwin Hubble created his famous classification method. In this work, we combine accurate visual classifications of the Galaxy Zoo project with \emph {Deep Learning} methods. The goal is to find an efficient technique at human performance level classification, but in a systematic and automatic way, for classification of elliptical and spiral galaxies. For this, a neural network model was created through an Ensemble of four other convolutional models, allowing a greater accuracy in the classification than what would be obtained with any one individual. Details of the individual models and improvements made are also described. The present work is entirely based on the analysis of images (not parameter tables) from DR1 (www.datalab.noao.edu) of the Southern Photometric Local Universe Survey (S-PLUS). In terms of classification, we achieved, with the Ensemble, an accuracy of $\approx 99 \%$ in the test sample (using pre-trained networks).

We pursue a program to confront observations with arbitrarily inhomogeneous cosmologies beyond the FLRW metric. The main idea is to test the Copernican principle rather than assuming it a priori. We consider the $\Lambda$CDM model endowed with a spherical $\Lambda$LTB inhomogeneity around us, that is, we assume isotropy and test the hypothesis of homogeneity. We confront the $\Lambda$LTB model with the latest available data from CMB, BAO, type Ia supernovae, local $H_0$, cosmic chronometers, Compton y-distortion and kinetic Sunyaev-Zeldovich effect. We find that these data can constrain tightly this extra inhomogeneity, almost to the cosmic variance level: on scales $\gtrsim 100$ Mpc structures can have a small non-Copernican effective contrast of just $\delta_L \sim 0.01$. Furthermore, the constraints on the standard $\Lambda$CDM parameters are not weakened after marginalizing over the parameters that model the local structure, to which we assign ignorance priors. In other words, dropping the FLRW metric assumption does not imply worse constraints on the cosmological parameters. This positive result confirms that the present and future data can be meaningfully analyzed within the framework of inhomogeneous cosmology.

ATLAS is the first Portuguese radar system that aims to detect space debris. The article introduces the system and provides a brief description of its capabilities. The system is capable of synthesizing arbitrary amplitude modulated pulse shapes with a resolution of 10 ns. Given that degree of freedom we decided to test an amplitude modulated chirp signal developed by us and a nested barker code. These waveforms are explained as well as their advantages and drawbacks for space debris detection. An experimental setup was developed to test the system receiver and waveforms are processed by digital matched filtering. The experiments test the system using different waveform shapes and noise levels. Experimental results are in agreement with simulation and show that the chirp signal is more resilient to Doppler shifts, has higher range resolution and lower peak-to-sidelobe ratio in comparison with the nested barker code. Future work in order to increase detection capabilities is discussed at the end.

As part of the Portuguese Space Surveillance and Tracking (SST) System, two new Wide Field of View (2.3deg x 2.3deg) small aperture (30cm) telescopes will be deployed in 2021, at the Pampilhosa da Serra Space Observatory (PASO), located in the center of the continental Portuguese territory, in the heart of a certified Dark Sky area. These optical systems will provide added value capabilities to the Portuguese SST network, complementing the optical telescopes currently in commissioning in Madeira and Azores. These telescopes are optimized for GEO and MEO survey operations and besides the required SST operational capability, they will also provide an important development component to the Portuguese SST network. The telescopes will be equipped with filter wheels, being able to perform observations in several optical bands including white light, BVRI bands and narrow band filters such as H(alpha) and O[III] to study potential different objects' albedos. This configuration enables us to conduct a study on space debris classification$/$characterization using combinations of different colors aiming the production of improved color index schemes to be incorporated in the automatic pipelines for classification of space debris. This optical sensor will also be used to conduct studies on image processing algorithms, including source extraction and classification solutions through the application of machine learning techniques. Since SST dedicated telescopes produce a large quantity of data per observation night, fast, efficient and automatic image processing techniques are mandatory. A platform like this one, dedicated to the development of Space Surveillance studies, will add a critical capability to keep the Portuguese SST network updated, and as a consequence it may provide useful developments to the European SST network as well.

The complexity of constraining the stellar initial mass function (IMF) in early-type galaxies cannot be overstated, given the necessity of both very high signal-to-noise (S/N) data and the difficulty of breaking the strong degeneracies that occur among several stellar population parameters including age, metallicity and elemental abundances. With this paper, the second in a series, we present a detailed analysis of the biases that can occur when retrieving the IMF shape by exploiting both optical and NIR IMF sensitive spectral indices. As a test case, here we analyze data for the nearby galaxy M89, for which we have high S/N spectroscopic data that cover the 3500-9000{\AA} spectral region and allow us to study the radial variation of the stellar population properties out to 1 R_e. Carrying out parallel simulations that mimic the retrieval of all the explored stellar parameters from a known input model, we quantify the amount of bias at each step of our analysis. From more general simulations we conclude that to accurately retrieve the IMF, it is necessary not only to retrieve accurate estimates of the age and metallicity, but also of all the elemental abundances that the spectral index fits are sensitive to. With our analysis technique applied to M89, we find consistency with a bottom-heavy IMF with a negative gradient from the center to half R_e when using the Conroy et al. 2018 as well as Vazdekis et al. 2016 EMILES stellar population models. We find agreement both with a parallel full spectral fitting of the same data and with literature results.

To study the quality of stellar spectra of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and the correctness of the corresponding stellar parameters derived by the LASP (LAMOST Stellar Parameter Pipeline), the outlier analysis method is applied to the archived AFGK stars in the fifth data release (DR5) of LAMOST. The outlier factor is defined in order to sort more than 3 million stellar spectra selected from the DR5 Stellar Parameter catalog. We propose an improved Local Outlier Factor (LOF) method based on Principal Component Analysis and Monte Carlo to enable the computation of the LOF rankings for randomly picked sub-samples that are computed in parallel by multiple computers, and finally to obtain the outlier ranking of each spectrum in the entire dataset. Totally 3,627 most outlier ranked spectra, around one-thousandth of all spectra, are selected and clustered into 10 groups, and the parameter density distribution of them conforms to the parameter distribution of LAMOST DR5, which suggests that in the whole parameter space the probability of bad spectra is uniformly distributed. By cross-matching the 3,627 spectra with APOGEE, we obtain 122 common ones. The published parameters calculated from LASP agree with APOGEE for the 122 spectra although there are bad pixels or bad flux calibrations in them. On the other hand, some outlier spectra show strong nebular contamination warning the corresponding parameters should be carefully used. A catalog and a spectral atlas of all the 3,627 outliers can be found at the link this http URL

Transitional millisecond pulsars are millisecond pulsars that switch between a rotation-powered millisecond pulsar state and an accretion-powered X-ray binary state, and are thought to be an evolutionary stage between neutron star low-mass X-ray binaries and millisecond pulsars. So far, only three confirmed systems have been identified in addition to a handful of candidates. We present the results of a multi-wavelength study of the low-mass X-ray binary NGC 6652B in the globular cluster NGC 6652, including simultaneous radio and X-ray observations taken by the Karl G. Jansky Very Large Array and the Chandra X-ray Observatory, and optical spectroscopy and photometry. This source is the second brightest X-ray source in NGC 6652 ($L_{\textrm{X}}\sim1.8\times10^{34}$ erg s$^{-1}$) and is known to be variable. We observe several X-ray flares over the duration of our X-ray observations, in addition to persistent radio emission and occasional radio flares. Simultaneous radio and X-ray data show no clear evidence of anti-correlated variability. Optical spectra of NGC 6652B indicate variable, broad H $\alpha$ emission which transitions from double-peaked emission to absorption over a time-scale of hours. We consider a variety of possible explanations for the source behaviour, and conclude that based on the radio and X-ray luminosities, short time-scale variability and X-ray flaring, and optical spectra, NGC 6652B is best explained as a transitional millisecond pulsar candidate that displays prolonged X-ray flaring behaviour. However, this could only be confirmed with observations of a change to the rotation-powered millisecond pulsar state.

We present a new catalog of 11209 O-rich AGB stars and 7172 C-rich AGB stars in our Galaxy identifying more AGB stars in the bulge component and considering more visual carbon stars. For each object, we cross-identify the IRAS, AKARI, MSX, WISE, 2MASS, and AAVSO counterparts. We present the new catalog in two parts: one is based on the IRAS PSC for brighter or more isolated objects, the other one is based on the ALLWISE source catalog for less bright or objects in crowded regions. We present various infrared two-color diagrams (2CDs) for the sample stars. We find that the theoretical dust shell models can roughly explain the observations of AGB stars on the various IR 2CDs. We investigate IR properties of SiO and OH maser emission sources in the catalog. For Mira variables in the sample stars, we find that the IR colors get redder for longer pulsation periods. We also study infrared variability of the sample stars using the WISE photometric data in the last 12 years: the ALLWISE multiepoch data and the Near-Earth Object WISE Reactivation (NEOWISE-R) 2021 data release. We generate light curves using the WISE data at W1 and W2 bands and compute the Lomb-Scargle periodograms for all of the sample stars. From the WISE light curves, we have found useful variation parameters for 3710 objects in the catalog, for which periods were either known or unknown in previous works.

Protoplanetary disks drive some of the formation process (e.g., accretion, gas dissipation, formation of structures, etc.) of stars and planets. Understanding such physical processes is one of the main astrophysical questions. HD 163296 is an interesting young stellar object for which infrared and sub-millimeter observations have shown a prominent circumstellar disk with gaps plausibly created by forming planets. This study aims at characterizing the morphology of the inner disk in HD 163296 with multi-epoch near-infrared interferometric observations performed with GRAVITY at the Very Large Telescope Interferometer (VLTI). Our goal is to depict the K-band (lambda_0 ~ 2.2 um) structure of the inner rim with milliarcsecond (sub-au) angular resolution. Our data is complemented with archival PIONIER (H-band; lambda_0 ~ 1.65 um) data of the source. We performed a Gradient Descent parametric model fitting to recover the sub-au morphology of our source. Our analysis shows the existence of an asymmetry in the disk surrounding the central star of HD 163296. We confirm variability of the disk structure in the inner ~2 mas (0.2 au). While variability of the inner disk structure in this source has been suggested by previous interferometric studies, this is the first time that it is confirmed in the H- and K-bands by using a complete analysis of the closure phases and squared visibilities over several epochs. Because of the separation from the star, position changes, and persistence of this asymmetric structure on timescales of several years, we argue that it is a dusty feature (e.g., a vortex or dust clouds), probably, made by a mixing of sillicate and carbon dust and/or refractory grains, inhomogeneously distributed above the mid-plane of the disk.

We aim to study the nature of the faint, polarised radio source population whose source composition and redshift dependence contain information about the strength, morphology, and evolution of magnetic fields over cosmic timescales. We use a 15 pointing radio continuum L-band mosaic of the Lockman Hole, observed in full polarisation, generated from archival data of the WSRT. The data were analysed using the RM-Synthesis technique. We achieved a noise of 7 {\mu}Jy/beam in polarised intensity, with a resolution of 15''. Using infrared and optical images and source catalogues, we were able to cross-identify and determine redshifts for one third of our detected polarised sources. We detected 150 polarised sources, most of which are weakly polarised with a mean fractional polarisation of 5.4 %. With a total area of 6.5 deg^2 and a detection threshold of 6.25 {\sigma} we find 23 polarised sources per deg^2. Based on our multi wavelength analysis, we find that our sample consists of AGN only. We find a discrepancy between archival number counts and those present in our data, which we attribute to sample variance. Considering the absolute radio luminosty, to distinguish weak and strong sources, we find a general trend of increased probability to detect weak sources at low redshift and strong sources at high redshift. Further, we find an anti-correlation between fractional polarisation and redshift for our strong sources sample at z{\geq}0.6. A decrease in the fractional polarisation of strong sources with increasing redshift cannot be explained by a constant magnetic field and electron density over cosmic scales, however the changing properties of cluster environments over the cosmic timemay play an important role. Disentangling these two effects requires deeper and wider polarisation observations, and better models of the morphology and strength of cosmic magnetic fields.

Theoretical prediction of surface stellar abundances of light elements -- lithium, beryllium, and boron -- represents one of the most interesting open problems in astrophysics. As well known, several measurements of 7-Li abundances in stellar atmospheres point out a disagreement between predictions and observations in different stellar evolutionary phases, rising doubts about the capability of present stellar models to precisely reproduce stellar envelope characteristics. Light elements are burned at relatively low temperatures (from 2 to 5 MK) and thus in the evolutionary stages of a star they are gradually destroyed at different depths of stellar interior, in dependence on the stellar mass. Their surface abundances are strongly influenced by the nuclear cross sections, by the extension of the convective envelope and by the temperature at its bottom, which depend on the characteristics of the star (mass and chemical composition) as well as on the energy transport in the convective stellar envelope. In recent years, a great effort has been made to improve the precision of light element burning cross sections. However, theoretical predictions surface light element abundance are challenging because they are also influenced by the uncertainties in the input physics adopted in the calculations as well as the efficiency of several standard and non-standard physical processes active in young stars (i.e. diffusion, radiative levitation, magnetic fields, rotation). Moreover, it is still not completely clear how much the previous protostellar evolution affects the characteristics of a stellar model and thus the light element depletion. This paper presents the state-of-the-art of theoretical predictions for protostars and PMS stars, discussing the role of several input physics on the stellar evolution.

Solar filaments are cold and dense materials situated in magnetic dips, which show distinct radiation characteristics compared to the surrounding coronal plasma. They are associated with coronal sheared and twisted magnetic field lines. However, the exact magnetic configuration supporting a filament material is not easy to be ascertained because of the absence of routine observations of the magnetic field inside filaments. Since many filaments lie above weak-field regions, it is nearly impossible to extrapolate their coronal magnetic structures by applying the traditional methods to noisy photospheric magnetograms, in particular the horizontal components. In this paper, we construct magnetic structures for some filaments with the regularized Biot--Savart laws and calculate their magnetic twists. Moreover, we make a parameter survey for the flux ropes of the Titov-Demoulin-modified model to explore the factors affecting the twist of a force-free magnetic flux rope. It is found that the twist of a force-free flux rope is proportional to its axial length to minor radius ratio, and is basically independent of the overlying background magnetic field strength. Thus, we infer that long quiescent filaments are likely to be supported by more twisted flux ropes than short active-region filaments, which is consistent with observations.

Starting in summer 2021, the Radio Neutrino Observatory in Greenland (RNO-G) will search for astrophysical neutrinos at energies >10 PeV by detecting the radio emission from particle showers in the ice around Summit Station, Greenland. We present an extensive simulation study that shows how RNO-G will be able to measure the energy of such particle cascades, which will in turn be used to estimate the energy of the incoming neutrino that caused them. The location of the neutrino interaction is determined using the differences in arrival times between channels and the electric field of the radio signal is reconstructed using a novel approach based on Information Field Theory. Based on these properties, the shower energy can be estimated. We show that this method can achieve an uncertainty of 13% on the logarithm of the shower energy after modest quality cuts and estimate how this can constrain the energy of the neutrino. The method presented in this paper is applicable to all similar radio neutrino detectors, such as the proposed radio array of IceCube-Gen2.

We present a method that isolates time-varying components from coronagraph and EUV images, allowing sub-streamer transients propagating within streamers to be tracked from the low to high corona. The method uses a temporal bandpass filter with a transmission bandwidth of ~2.5-10 hours that suppresses both high and low frequency variations in observations made by the STEREO/SECCHI suite. We demonstrate that this method proves crucial in linking the low corona where the magnetic field is highly non-radial, to their counterparts in the high corona where the magnetic field follows a radial path through the COR1 instrument. We also applied our method to observations by the COR2 and EUVI instruments onboard SECCHI and produced height-time profiles that revealed small density enhancements, associated with helmet streamers, propagating from ~1.2 Rs out to beyond 5 Rs. Our processing method reveals that these features are common during the period of solar minimum in this study. The features recur on timescales of hours, originate very close to the Sun, and remain coherent out into interplanetary space. We measure the speed of the features and classify them as: slow (a few to tens of km/s) and fast (~100 km/s). Both types of features serve as an observable tracer of a variable component of the slow solar wind to its source regions. Our methodology helps overcome the difficulties in tracking small-scale features through COR1. As a result, it proved successful in measuring the connectivity between the low and high corona and in measuring velocities of small-scale features.

Context: We present a new study of the eclipsing cataclysmic variable CzeV404 Her (Porb = 0.098 d) that is located in the period gap. Aims: This report determines the origin of the object and the system parameters and probes the accretion flow structure of the system. Methods: We conducted simultaneous time-resolved photometric and spectroscopic observations of CzeV404 Her. We applied our light-curve modelling techniques and the Doppler tomography method to determine the system parameters and analyse the structure of the accretion disk. Results: We found that the system has a massive white dwarf M_WD = 1.00(2) M_sun a mass ratio of q = 0.16, and a relatively hot secondary with an effective temperature T_2 = 4100(50) K. The system inclination is i = 78.8{\deg}. The accretion disk spreads out to the tidal limitation radius and has an extended hot spot or line region. The hot spot or line is hotter than the remaining outer part of the disk in quiescence or in intermediate state, but does not stand out completely from the disk flux in (super)outbursts. Conclusions: We claim that this object represents a link between two distinct classes of SU UMa-type and SW Sex-type cataclysmic variables. The accretion flow structure in the disk corresponds to the SW Sex systems, but the physical conditions inside the disk fit the behaviour of SU UMa-type objects.

We study the spherically averaged bispectrum of the 21-cm signal from the Epoch of Reionization (EoR). This metric provides a quantitative measurement of the level of non-Gaussianity of the signal which is expected to be high. We focus on the impact of the light-cone effect on the bispectrum and its detectability with the future SKA-Low telescope. Our investigation is based on semi-numerical light-cone simulation and an ensemble of 50 independent realisations of the 21-cm signal to estimate the cosmic variance errors. We calculate the bispectrum with a new, optimised direct estimation method, DviSukta which calculates the bispectrum for all possible unique triangles. We find that the light-cone effect becomes important on scales $k_1 \lesssim 0.1\,{\rm Mpc}^{-1}$ where for most triangle shapes the cosmic variance errors dominate. Only for the squeezed limit triangles, the impact of the light-cone effect exceeds the cosmic variance. Combining the effects of system noise and cosmic variance we find that $\sim 3\sigma$ detection of the bispectrum is possible for all unique triangle shapes around a scale of $k_1 \sim 0.2\,{\rm Mpc}^{-1}$, and cosmic variance errors dominate above and noise errors below this length scale. Only the squeezed limit triangles are able to achieve a more than $5\sigma$ significance over a wide range of scales, $k_1 \lesssim 0.8\,{\rm Mpc}^{-1}$. Our results suggest that among all the possible triangle combinations for the bispectrum, the squeezed limit one will be the most measurable and hence useful.

In the past, Kepler painstakingly derived laws of planetary motion using difficult to understand and hard to follow techniques. In 1843 William Hamilton created and described the quaternions, which extend the complex numbers and can easily describe rotations in three dimensional space. In this article, we will harness this system to provide a new and intuitive way to derive Kepler's laws. This will include using a quaternionic version of the spatial Kepler problem differential equation, and using the general solution to describe the motion of planets orbiting a central body. We use the standard method for regularizing celestial mechanics, but this article will be solely focused on showing the validity of Kepler's laws.

We consider the effects of uniform rapid stellar rotation on the minimum mass of stable hydrogen burning. To focus on the effects of rotation, we use a polytropic model of the star, and employ a suitable numerical scheme, relaxing the assumption of spherical symmetry. We obtain an analytical formula for the minimum mass of hydrogen burning as a function of the stellar rotation speed. Further, we show the possibility of a maximum mass of stable hydrogen burning in such stars, which is purely an artefact of rotation. The existence of this extremum in mass results in a maximum admissible value of the stellar rotation speed, beyond which no brown dwarfs can exist, within the ambits of our model. For a given angular speed, we predict a mass range beyond which no brown dwarf will evolve into a main sequence star.

Henize 2-10 (He 2-10) is a nearby (D = 9 Mpc) starbursting blue compact dwarf galaxy that boasts a high star formation rate and a low luminosity AGN. He 2-10 is also one of the first galaxies in which embedded superstar clusters (SSCs) were discovered. SSCs are massive, compact star clusters that will impact their host galaxies dramatically when their massive stars evolve. Here, we discuss radio, submillimeter, and infrared observations of He 2-10 from 1.87 microns to 6 cm in high angular resolution (~0.3 arcsec), which allows us to disentangle individual clusters from aggregate complexes as identified at lower resolution. These results indicate the importance of spatial resolution to characterize SSCs, as low resolution studies of SSCs average over aggregate complexes that may host SSCs at different stages of evolution. We explore the thermal, non-thermal, and dust emission associated with the clusters along with dense molecular tracers to construct a holistic review of the natal SSCs that have yet to dramatically disrupt their parent molecular clouds. We assess the production rate of ionizing photons, extinction, total mass, and the star formation efficiency associated with the clusters. Notably, we find that the star formation efficiency for the some of the natal clusters is high (>70%), which suggests that these clusters could remain bound even after the gas is dispersed from the system from stellar feedback mechanisms. If they remain bound, these SSCs could survive to become objects indistinguishable from globular clusters.

A detailed study of the BL Lacertae PKS 0903-57 has been done for the first time with 12 years of Fermi Large Area Telescope data. We have identified two bright gamma-ray flares in 2018 and 2020. Many sub-structures were observed during multiple time binning of these flares. We have performed detailed temporal and spectral study on all the sub-structures separately. A single-zone emission model is used for time-dependent leptonic modeling of the multi-wavelength spectral energy distributions. Our estimated values of variability time scale, magnetic field in the emission region, jet power obtained from leptonic modeling of PKS 0903-57 are presented in this work. Currently, we have a minimal number of observations in X-rays and other bands. Hence, more simultaneous multi-wavelength monitoring of this source is required to have a better understanding of the physical processes happening in the jet of the blazar PKS 0903-57.

The origin of the heavy elements in the Universe is not fully determined. Neutron star-black hole (NSBH) and neutron star-neutron star mergers may both produce heavy elements via rapid neutron-capture process (r-process). We use the recent detection of gravitational waves from NSBHs, improved measurements of neutron star equation-of-state, and the most modern numerical simulations of the ejected materials from binary collisions to investigate the production of heavy elements from binary mergers. As the amount of ejecta depends on the mass and spin distribution of compact objects, as well as on the equation-of-state of neutron stars, we consider various models for these quantities, informed by gravitational-wave and pulsar data. We find that even in the most favorable model, neutron star-black hole mergers are unlikely to account for more than 77% of the r-process elements in the local Universe. If black holes have preferentially small spins, this bound decreases to 35%. Finally, if black hole spins are small and there is a dearth of low mass ($<5M_{\odot}$) black holes in NSBH binaries, the NSBH contribution to r-process elements is negligible.

Cold, non-self-gravitating clumps occur in various astrophysical systems, ranging from the interstellar and circumgalactic medium (CGM), to AGN outflows and solar coronal loops. Cold gas has diverse origins such as turbulent mixing or precipitation from hotter phases. We obtain the analytic solution for a steady pressure-driven 1-D cooling flow around cold over-densities, irrespective of their origin. Our solutions describe the slow and steady radiative cooling-driven local gas inflow in the saturated regime of nonlinear thermal instability in clouds, sheets and filaments. We use a simple two-fluid treatment to include magnetic fields as an additional polytropic fluid. To test the limits of applicability of these analytic solutions, we compare with the gas structure found in and around small-scale cold clouds in the CGM of massive halos in the TNG50 cosmological MHD simulation from the IllustrisTNG suite. Despite qualitative resemblance of the gas structure, we find that deviations from steady state, complex geometries and turbulence all add complexity beyond our analytic solutions. We derive an exact relation between the mass cooling rate ($\dot{\rm M}_{\rm cool}$) and the radiative cooling rate ($\dot{\rm E}_{\rm cool}$) for a steady cooling flow. A comparison with the TNG50 clouds shows that this cooling flow relation applies in a narrow temperature range around $\rm \sim 10^{4.5}$ K where the isobaric cooling time is the shortest. In general, turbulence and mixing, instead of radiative cooling, may dominate the transition of gas between different temperature phases.

A debate is emerging regarding the recent inconsistent results of different studies for the Cosmic Star Formation Rate Density (CSFRD) at high-z. We employ UV and IR datasets to investigate the star formation rate function (SFRF) at ${\rm z \sim 0-9}$. We find that the SFRFs derived from the dust corrected ${\rm UV}$ (${\rm UV_{corr}}$) data contradict those from IR on some key issues since they are described by different distributions (Schechter vs double-power law), imply different physics for galaxy formation (${\rm UV_{corr}}$ data suggest a SFR limit/strong mechanism that diminish the number density of high star forming systems with respect IR) and compare differently with the stellar mass density evolution obtained from SED fitting (${\rm UV_{corr}}$ is in agreement, while IR in tension up to 0.5 dex). However, both tracers agree on a constant CSFRD evolution at ${\rm z \sim 1-4}$ and point to a plateau instead of a peak. In addition, using both indicators we demonstrate that the evolution of the {\it observed} CSFRD can be described by only {\bf 2} parameters and a function that has the form of a Gamma distribution (${\bf \Gamma(a,bt)}$). In contrast to previous parameterizations used in the literature our framework connects the parameters to physical properties like the star formation rate depletion time and cosmic baryonic gas density. The build up of stellar mass occurs in $\Gamma(a,bt)$ distributed steps and is the result of gas consumption up to the limit that there is no eligible gas for SF at t = ${\rm \infty}$, resulting to a final cosmic stellar mass density of $\sim 0.5 \times 10^9 \, {\rm \frac{M_{\odot}}{Mpc^3}}$.

We present the first short-duration candidate microlensing events from the Kepler K2 mission. From late April to early July 2016, Campaign 9 of K2 obtained high temporal cadence observations over a 3.7 square degree region of the Galactic bulge. Its primary objectives were to look for evidence of a free-floating planet (FFP) population using microlensing, and demonstrate the feasibility of space-based planetary microlensing surveys. Though Kepler K2 is far from optimal for microlensing, the recently developed MCPM photometric pipeline enables us to identify and model microlensing events. We describe our blind event-selection pipeline in detail and use it to recover 22 short-duration events with effective timescales of less than 10 days previously announced by the OGLE and KMTNet ground-based surveys. We also announce five new candidate events. One of these is a caustic-crossing binary event, consistent with a bound planet and modelled as such in a companion study. The other four have very short durations (effective timescales less than 0.1 days) typical of an Earth-mass FFP population. Whilst Kepler was not designed for crowded-field photometry, the K2C9 dataset clearly demonstrates the feasibility of conducting blind space-based microlensing surveys towards the Galactic bulge.

We present a pulsar candidate identification and confirmation procedure based on a position-switch mode during the pulsar search observations. This method enables the simultaneous search and confirmation of a pulsar in a single observation, by utilizing the different spatial features of a pulsar signal and a radio frequency interference (RFI). Based on this method, we performed test pulsar search observations in globular clusters M3, M15, and M92. We discovered and confirmed a new pulsar, M3F, and detected the known pulsars M3B, M15 A to G (except C), and M92A.

The magnetic fields of accretion disks play an important role in studying their evolution. We may assume that its generation is connected to the dynamo mechanism, which is similar with that in the galactic disks. Here, we propose a model of the magnetic field of the accretion disk that uses the same approaches that have been used for galaxies. It is necessary to obtain the field, which is expected to be less than the equipartition value, and without destroying the disk. To do so, it is necessary to formulate the basic properties of the ionized medium and to estimate the parameters governing the dynamo. We used the no-z approximation that has been developed for thin disks. We also take different boundary conditions that can change the value of the field significantly. We show that the magnetic field strictly depends on the boundary conditions. Taking zero conditions and the fixed magnetic field condition on the inner boundary, which are connected to the physical properties of the accretion disk, we can avoid solutions that are greater than the equipartition field.

In this work, we focused on observations of six trans-Neptunian objects (TNOs) whose apparent magnitudes are brighter than 20 mag. We present the results of astrometric and photometric observations of (134340) Pluto, (136108) Haumea, (136472) Makemake, (136199) Eris, (90482) Orcus and (20000) Varuna obtained at the Kyiv comet station (Code MPC 585) in 2017-2019. For observations, we used the 0.7-m reflector AZT-8 with the FLI PL4710 CCD camera and filters of the Johnson-Cousins photometric system. From our images, we measured the objects astrometric positions, calculated apparent magnitudes in the BVRI (mostly R) bands using the aperture photometry method, and found the absolute magnitudes together with the colour indices in several bands. Analysing our results, we investigate the limitation on the astrometry and photometry of faint objects with the 0.7-m telescope.

Everything we know about the environment around us is thanks to light, a kind of electromagnetic radiation. Astronomy takes advantage of it and all the electromagnetic spectrum with the help of many devices to record them and determine from which places in our Universe they come. These signals must be processed to obtain the images that are will be then exposed to the public. Astronomers know that these images will inspire and generate curiosity in each person who sees them. This Science is inclusive and wants to transmit the knowledge and the beautiful events that happen in the universe to all people. We try to do that by developing and showing other forms of teaching, taking advantage of new technologies that are available today to bring this knowledge closer to the minorities in our country like visually impaired people. Because of this problem, the AstroBVI project arrived in Peru for the first time, thanks to the distribution of tactile images of 3D galaxies that were delivered to us from the Centro de Astronom\'ia de la Universidad de Antofagasta of the and financed by the International Astronomical Union Office of Astronomy for Development (IAU-OAD). This allowed the holding of seven workshops during 2019, visiting various institutions and benefiting more than 160 participants with blindness and low vision, identifying in them a lot of interest, which shows the enormous potential presented by these 3D tactile materials.

The characteristics of the stellar populations in the Galactic Bulge inform and constrain the Milky Way's formation and evolution. The metal-poor population is particularly important in light of cosmological simulations, which predict that some of the oldest stars in the Galaxy now reside in its center. The metal-poor bulge appears to consist of multiple stellar populations that require dynamical analyses to disentangle. In this work, we undertake a detailed chemodynamical study of the metal-poor stars in the inner Galaxy. Using R$\sim$ 20,000 VLT/GIRAFFE spectra of 319 metal-poor (-2.55 dex$\leq$[Fe/H]$\leq$0.83 dex, with $\overline{\rm{[Fe/H]}}$=-0.84 dex) stars, we perform stellar parameter analysis and report 12 elemental abundances (C, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Zn, Ba, and Ce) with precisions of $\approx$0.10 dex. Based on kinematic and spatial properties, we categorise the stars into four groups, associated with the following Galactic structures: the inner bulge, the outer bulge, the halo, and the disk. We find evidence that the inner and outer bulge population is more chemically complex (i.e., higher chemical dimensionality and less correlated abundances) than the halo population. This result suggests that the older bulge population was enriched by a larger diversity of nucleosynthetic events. We also find one inner bulge star with a [Ca/Mg] ratio consistent with theoretical pair-instability supernova yields and two stars that have chemistry consistent with globular cluster stars.

We study quantized vortices in ${}^{3}P_{2}$ superfluids using a microscopic theory for the first time. The theory is based on the Eilenberger equation to determine the order parameters and the Bogoliubov-de Gennes (BdG) equation to obtain the eigenenergies and the core magnetization. Within axisymmetric vortex configurations, we find several stable and metastable vortex configurations which depend on the strength of a magnetic field, similar to a $v$ vortex and $o$ vortex in $^3$He superfluids. We demonstrate that the $o$ vortex is the most stable axisymmetric vortex in the presence of a strong magnetic field, and we find two zero-energy Majorana fermion bound states in the $o$-vortex core. We show that the profiles of the core magnetization calculated using the BdG equation are drastically different from those calculated using only the order parameter profiles known before.

We investigate the potential of type II supernovae (SNe) to constrain axion-like particles (ALPs) coupled simultaneously to nucleons and electrons. ALPs coupled to nucleons can be efficiently produced in the SN core via nucleon-nucleon bremsstrahlung and, for a wide range of parameters, leave the SN unhindered, producing a large ALP flux. For masses exceeding 1 MeV, these ALPs would decay into electron-positron pairs, generating a positron flux. In the case of Galactic SNe, the annihilation of the created positrons with the electrons present in the Galaxy would contribute to the 511 keV annihilation line. Using the SPI (SPectrometer on INTEGRAL) observation of this line, allows us to exclude a wide range of the axion-electron coupling, $10^{-19} \lesssim g_{ae} \lesssim 10^{-11}$, for $g_{ap}\sim 10^{-9}$. Additionally, ALPs from extra-galactic SNe decaying into electron-positron pairs would yield a contribution to the cosmic X-ray background. In this case, we constrain the ALP-electron coupling down to $g_{ae} \sim 10^{-20}$.

We investigate the spectrum of linearized excitations of global vortices in $2+1$ dimensions. After identifying the existence of two localized excitation modes, we compute their decay time scale and compare the results to the numerical evolution of the full non-linear equations. We show numerically how the interaction of vortices with an external source of radiation or other vortices can excite these modes dynamically. We then simulate the formation of vortices in a phase transition and their interaction with a thermal bath estimating the amplitudes of these modes in each case. These numerical experiments indicate that even though, in principle, vortices are capable of storing a large amount of energy in these internal excitations, this does not seem to happen dynamically. We then explore the evolution of a network of vortices in an expanding (2+1) dimensional background, in particular in a radiation dominated universe. We find that vortices are still excited after the course of the cosmological evolution but again the level of excitation is very small. The extra energy in the vortices in these cosmological simulations never exceeds the $1\%$ level of the total mass of the core of the vortex.

Using a extended version of the Quantum Hadrodynamics (QHD), I propose a new microscopic equation of state (EoS) able to correctly reproduce the main properties of symmetric nuclear matter at the saturation density, as well produce massive neutron stars and satisfactory results for the radius and the tidal parameter $\Lambda$. I show that even when hyperons are present, this EoS is able to reproduce at least 2.00 solar masses neutron star. The constraints about the radius of a $2.00M_\odot$ and the minimum mass that enable direct Urca effect are also checked.

Detection of a stochastic background of gravitational waves is likely to occur in the next few years. Beyond searches for the isotropic component of SGWBs, there have been various mapping methods proposed to target anisotropic backgrounds. Some of these methods have been applied to data taken by the Laser Interferometer Gravitational-wave Observatories (LIGO) and Virgo. Specifically, these directional searches have focused on mapping the intensity of the signal on the sky via maximum likelihood solutions. We compare this intensity mapping approach to a previously proposed, but never employed, amplitude-phase mapping method to understand whether this latter approach may be employed in future searches. We build up our understanding of the differences between these two approaches by analysing simple toy models of time-stream data, and run mock-data mapping tests for the two methods. We find that the amplitude-phase method is only applicable to the case of a background which is phase-coherent on large scales or, at the very least, has an intrinsic coherence scale that is larger than that of the detector. Otherwise, the amplitude-phase mapping method leads to a loss of overall information, with respect to both phase and amplitude. Since we do not expect these phase-coherent properties to hold for any of the gravitational-wave background signals we hope to detect in the near future, we conclude that intensity mapping is the preferred method for such backgrounds.

A double-phase argon detector is excellent in particle identification and position reconstruction. However, the properties of the electroluminescence (EL) process for secondary light emission in the gas phase are not fully understood. The EL process was thought to be explained using an ordinary EL mechanism because of an argon excimer, but there were no visible light (VL) emissions in this mechanism. However, recent measurements indicated there were visible components in the argon gas electroluminescence, which was proposed to explain the visible light components by a new mechanism called neutral bremsstrahlung (NBrS). In this article, we studied gaseous argon electroluminescence in the VL region from 300 to 600 nm at room temperature and normal pressure using a gaseous time projection chamber (TPC). The secondary emission light from the TPC luminescence region was dispersed using a spectrometer. Then, the interpretation of the observed spectrum using the ordinary EL model, NBrS model, and the effect of nitrogen impurity was discussed.

We investigate ion-scale kinetic plasma instabilities at the collisionless shock using linear theory and nonlinear Particle-in-Cell (PIC) simulations. We focus on the Alfv\'en-ion-cyclotron (AIC), mirror, and Weibel instabilities, which are all driven unstable by the effective temperature anisotropy induced by the shock-reflected ions within the transition layer of a strictly perpendicular shock. We conduct linear dispersion analysis with a homogeneous plasma model to mimic the shock transition layer by adopting a ring distribution with finite thermal spread to represent the velocity distribution of the reflected ions. We find that, for wave propagation parallel to the ambient magnetic field, the AIC instability at lower Alfv\'en Mach numbers tends to transition to the Weibel instability at higher Alfv\'en Mach numbers. The instability property is, however, also strongly affected by the sound Mach number. We conclude that the instability at a strong shock with Alfv\'en and sound Mach numbers both in excess of $\sim 20{\rm -}40$ may be considered as Weibel-like in the sense that the reflected ions behave essentially unmagnetized. Two-dimensional PIC simulations confirm the linear theory and find that, with typical parameters of young supernova remnant shocks, the ring distribution model produces magnetic fluctuations of the order of the background magnetic field, which is smaller than those observed in previous PIC simulations for Weibel-dominated shocks. This indicates that the assumption of the gyrotropic reflected ion distribution may not be adequate to quantitatively predict nonlinear behaviors of the dynamics in high Mach number shocks.

Gravitational waves from the distant sources are gravitationally lensed during their propagation through the intervening matter inhomogeneities before arriving at detectors. It has been proposed in the literature that the variance of the lensed waveform can be used to extract information of the matter power spectrum at very small scales and of low-mass dark halos. In this paper, we show that the variance of the amplitude fluctuation and that of the phase fluctuation of the lensed waveform obey a simple relation irrespective of the shape of the matter power spectrum. We study conditions under which this relation can be violated and discuss some potential applications of the relation. This relation may be used to confirm the robustness of claimed observations of gravitational lensing of gravitational waves and the subsequent reconstruction of the matter power spectrum.

The existence and stability of non-Abelian half-quantum vortices (HQVs) are established in ${}^{3}P_{2}$ superfluids in neutron stars with strong magnetic fields, the largest topological quantum matter in our Universe. Using a self-consistent microscopic framework, we find that one integer vortex is energetically destabilized into a pair of two non-Abelian HQVs due to the strong spin-orbit coupled gap functions. We find a topologically protected Majorana fermion on each HQV, thereby providing two-fold non-Abelian anyons characterized by both Majorana fermions and a non-Abelian first homotopy group.

The Transition-Edge Sensor (TES) is an extremely sensitive device which is used to measure the energy of individual X-ray photons. For astronomical spectrometry applications, SRON develops a Frequency Domain Multiplexing (FDM) read-out system for kilopixel arrays of such TESs. Each TES is voltage biased at a specific frequency in the range 1 to 5 MHz. Isolation between the individual pixels is obtained through very narrow-band (high-Q) lithographic LC resonators. To prevent energy resolution degradation due to intermodulation line noise, the bias frequencies are distributed on a regular grid. The requirements on the accuracy of the LC resonance frequency are very high. The deviation of the resonance frequencies due to production tolerances is significant with respect to the bandwidth, and a controller is necessary to compensate for the LC series impedance. We present two such controllers: a simple orthogonal proportional-integrating (PI) controller and a more complex impedance estimator. Both controllers operate in baseband and try to make the TES current in-phase with the bias voltage, effectively operating as phase-locked loops (PLL). They allow off-LC-resonance operation of the TES pixels, while preserving TES thermal response and energy resolution. Extensive experimental results -- published in a companion paper recently -- with the proposed methods, show that these controllers allow the preservation of single pixel energy resolution in multiplexed operation.

We describe the design of a gravitational wave timing array, a novel scheme that can be used to search for low-frequency gravitational waves by monitoring continuous gravitational waves at higher frequencies. We show that observations of gravitational waves produced by Galactic binaries using a space-based detector like LISA provide sensitivity in the nanohertz to microhertz band. While the expected sensitivity of this proposal is not competitive with other methods, it fills a gap in frequency space around the microhertz regime, which is above the range probed by current pulsar timing arrays and below the expected direct frequency coverage of LISA. The low-frequency extension of sensitivity does not require any experimental design change to space-based gravitational wave detectors, and can be achieved with the data products that would already be collected by them.