Several generalizations of the well-known fluid model of Braginskii (1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multispecies plasmas with arbitrary masses and temperatures, where all of the fluid moments are described by their evolution equations. The 21-moment model contains two "heat flux vectors" (third- and fifth-order moments) and two "viscosity tensors" (second- and fourth-order moments). The Braginskii model is then obtained as a particular case of a one ion–electron plasma with similar temperatures, with decoupled heat fluxes and viscosity tensors expressed in a quasistatic approximation. We provide all of the numerical values of the Braginskii model in a fully analytic form (together with the fourth- and fifth-order moments). For multispecies plasmas, the model makes the calculation of the transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasistatic approximation that makes some Braginskii coefficients divergent in a weakly collisional regime. Importantly, we show that the heat fluxes and viscosity tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the fourth-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad.

The Lomb–Scargle periodogram is a well-known algorithm for detecting and characterizing periodic signals in unevenly sampled data. This paper presents a conceptual introduction to the Lomb–Scargle periodogram and important practical considerations for its use. Rather than a rigorous mathematical treatment, the goal of this paper is to build intuition about what assumptions are implicit in the use of the Lomb–Scargle periodogram and related estimators of periodicity, so as to motivate important practical considerations required in its proper application and interpretation.

To date, 241 individual molecular species, composed of 19 different elements, have been detected in the interstellar and circumstellar medium by astronomical observations. These molecules range in size from two atoms to 70 and have been detected across the electromagnetic spectrum from centimeter wavelengths to the ultraviolet. This census presents a summary of the first detection of each molecular species, including the observational facility, wavelength range, transitions, and enabling laboratory spectroscopic work, as well as listing tentative and disputed detections. Tables of molecules detected in interstellar ices, external galaxies, protoplanetary disks, and exoplanetary atmospheres are provided. A number of visual representations of these aggregate data are presented and briefly discussed in context.

We determine the low-redshift X-ray luminosity function, active black hole mass function (BHMF), and Eddington ratio distribution function (ERDF) for both unobscured (Type 1) and obscured (Type 2) active galactic nuclei (AGNs), using the unprecedented spectroscopic completeness of the BAT AGN Spectroscopic Survey (BASS) data release 2. In addition to a straightforward 1/V_max approach, we also compute the intrinsic distributions, accounting for sample truncation by employing a forward-modeling approach to recover the observed BHMF and ERDF. As previous BHMFs and ERDFs have been robustly determined only for samples of bright, broad-line (Type 1) AGNs and/or quasars, ours are the first directly observationally constrained BHMF and ERDF of Type 2 AGNs. We find that after accounting for all observational biases, the intrinsic ERDF of Type 2 AGNs is significantly more skewed toward lower Eddington ratios than the intrinsic ERDF of Type 1 AGNs. This result supports the radiation-regulated unification scenario, in which radiation pressure dictates the geometry of the dusty obscuring structure around an AGN. Calculating the ERDFs in two separate mass bins, we verify that the derived shape is consistent, validating the assumption that the ERDF (shape) is mass-independent. We report the local AGN duty cycle as a function of mass and Eddington ratio, by comparing the BASS active BHMF with the local mass function for all supermassive black holes. We also present the $\mathrm{log}N-\mathrm{log}S$ of the Swift/BAT 70 month sources.

A complete accounting of nearby objects—from the highest-mass white dwarf progenitors down to low-mass brown dwarfs—is now possible, thanks to an almost complete set of trigonometric parallax determinations from Gaia, ground-based surveys, and Spitzer follow-up. We create a census of objects within a Sun-centered sphere of 20 pc radius and check published literature to decompose each binary or higher-order system into its separate components. The result is a volume-limited census of ∼3600 individual star formation products useful in measuring the initial mass function across the stellar (<8M_⊙) and substellar (≳5M_Jup) regimes. Comparing our resulting initial mass function to previous measurements shows good agreement above 0.8M_⊙ and a divergence at lower masses. Our 20 pc space densities are best fit with a quadripartite power law, $\xi (M)={dN}/{dM}\propto {M}^{-\alpha }$, with long-established values of α = 2.3 at high masses (0.55 < M < 8.00M_⊙), and α = 1.3 at intermediate masses (0.22 < M < 0.55M_⊙), but at lower masses, we find α = 0.25 for 0.05 < M < 0.22M_⊙, and α = 0.6 for 0.01 < M < 0.05M_⊙. This implies that the rate of production as a function of decreasing mass diminishes in the low-mass star/high-mass brown dwarf regime before increasing again in the low-mass brown dwarf regime. Correcting for completeness, we find a star to brown dwarf number ratio of, currently, 4:1, and an average mass per object of 0.41 M_⊙.

The value of the eclipse solar radius during the 2017 August 21 total solar eclipse was estimated to be S_⊙ = (95995 ± 005) at 1 au with no significant dependence on wavelength. The measurement was obtained from analysis of a video of the eclipse flash spectrum recorded at the southern limit of the umbral shadow path. Our analysis was conducted by extracting light curves from the flash spectrum and comparing them to simulated light curves. The simulations were performed by integrating the limb darkening function over the exposed area of the photosphere. These numerical integrations relied on very precise computations of the relative movement of the lunar and solar limbs.

Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultrawide millimeter spectral scan survey from the Atacama Large Millimeter/submillimeter Array Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 16 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH⁺ and the first interferometric images of C₃H⁺, NO, and HCS⁺. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our data set is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as (i) young starburst tracers characterized by high-excitation transitions of HC₃N and complex organic molecules versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing photodissociation regions, (ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH₃OH, HNCO, HOCO⁺, and OCS) versus shocks from star formation-induced outflows (high-excitation transitions of SiO), as well as (iii) outflows showing emission from HOC⁺, CCH, H₃O⁺, CO isotopologues, HCN, HCO⁺, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback.

With the advent of more sensitive all-sky instruments, the transient universe is being probed in greater depth than ever before. Taking advantage of available resources, we have established a comprehensive database of black hole (and black hole candidate) X-ray binary (BHXB) activity between 1996 and 2015 as revealed by all-sky instruments, scanning surveys, and select narrow-field X-ray instruments on board the INTErnational Gamma-Ray Astrophysics Laboratory, Monitor of All-Sky X-ray Image, Rossi X-ray Timing Explorer, and Swift telescopes; the Whole-sky Alberta Time-resolved Comprehensive black-Hole Database Of the Galaxy or WATCHDOG. Over the past two decades, we have detected 132 transient outbursts, tracked and classified behavior occurring in 47 transient and 10 persistently accreting BHs, and performed a statistical study on a number of outburst properties across the Galactic population. We find that outbursts undergone by BHXBs that do not reach the thermally dominant accretion state make up a substantial fraction (∼40%) of the Galactic transient BHXB outburst sample over the past ∼20 years. Our findings suggest that this "hard-only" behavior, observed in transient and persistently accreting BHXBs, is neither a rare nor recent phenomenon and may be indicative of an underlying physical process, relatively common among binary BHs, involving the mass-transfer rate onto the BH remaining at a low level rather than increasing as the outburst evolves. We discuss how the larger number of these "hard-only" outbursts and detected outbursts in general have significant implications for both the luminosity function and mass-transfer history of the Galactic BHXB population.

We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to streamline the extraction, preservation, integration and cross-linking of valuable data from astrophysics literature into major astronomical databases, and consequently facilitate new modes of science discovery that will better exploit the vast quantities of panchromatic and multidimensional data associated with the literature. We encourage authors, journal editors, referees, and publishers to implement the best practices reviewed here, as well as related recommendations from international astronomical organizations such as the International Astronomical Union for publication of nomenclature, data, and metadata. A convenient Checklist of Recommendations for Publishing Data in the Literature (Appendix A) is included for authors to consult before the submission of the final version of their journal articles and associated data files. We recommend that publishers of journals in astronomy and astrophysics incorporate a link to this document in their Instructions to Authors.

This paper presents a table with estimates of the absolute magnitude of the Sun and the conversions from vegamag to the AB and ST systems for several wide-band filters used in ground-based and space-based observatories. These estimates use the dustless spectral energy distribution (SED) of Vega, calibrated absolutely using the SED of Sirius, to set the vegamag zero-points and a composite spectrum of the Sun that coadds space-based observations from the ultraviolet to the near-infrared with models of the Solar atmosphere. The uncertainty of the absolute magnitudes is estimated by comparing the synthetic colors with photometric measurements of solar analogs and is found to be ∼0.02 mag. Combined with the uncertainty of ∼2% in the calibration of the Vega SED, the errors of these absolute magnitudes are ∼3%–4%. Using these SEDs, for three of the most utilized filters in extragalactic work the estimated absolute magnitudes of the Sun are M_B = 5.44, M_V = 4.81, and M_K = 3.27 mag in the vegamag system and M_B = 5.31, M_V = 4.80, and M_K = 5.08 mag in AB.

We use the spectroscopic data collected by the Magellanic Quasars Survey (MQS) and the photometric V- and I-band data from the Optical Gravitational Lensing Experiment (OGLE) to measure the physical parameters for active galactic nuclei (AGNs) located behind the Magellanic Clouds. The flux-uncalibrated MQS spectra were obtained with the 4 m Anglo-Australian Telescope and the AAOmega spectroscope (R = 1300) in a typical ∼1.5 hr visit. They span a spectral range of 3700–8500 Å and have signal-to-noise ratios in a range of 3–300. We report the discovery and observational properties of 161 AGNs in this footprint, which expands the total number of spectroscopically confirmed AGNs by MQS to 919. After the conversion of the OGLE mean magnitudes to the monochromatic luminosities at 5100, 3000, and 1350 Å, we were able to reliably measure the black hole masses for 165 out of 919 AGNs. The remaining physical parameters we provide are the bolometric luminosities and the Eddington ratios. A fraction of these AGNs have been observed by the OGLE survey since 1997 (all of them since 2001), enabling studies of correlations between the variability and physical parameters of these AGNs.

We implement an algorithm based on the weighted stacking of astronomical images that can combine different observations of the same region of the sky removing the interfering signals. We develop a C++ code that takes as input a set of spectral cubes and computes the local weights of the intensity for each pixel of every channel. The weights are calculated as the inverse variance of the nearby pixels and are used to compute the weighted merge of the input files. Astronomical sources, present in all cubes, are preserved by the weighted average. However, interfering signals, present in specific cubes and a certain frequency range, are down-weighted in the average and removed from the output spectral cube. We present the results obtained by analyzing simulated spectral cubes containing astronomical sources, noise, and a random set of interferences of different intensities and spectral occupations. The performance of the algorithm is evaluated by comparing the output result with the input sky-model image. Finally, we present the results obtained by applying the method to a set of real data consisting of observations of the Andromeda galaxy, Messier 31 (M31), at 6.6 GHz obtained with the Sardinia Radio Telescope.

We introduce Weakhub, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino–matter, neutrino–neutrino interactions and plasma process, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations of β-processes, incorporating various weak corrections and medium modifications due to the nuclear equation of state. Calculations of plasma processes, electron neutrino–antineutrino annihilation, and nuclear de-excitation are also included. We also present the detailed derivations of weak interactions and the coupling to the two-moment based general-relativistic multigroup radiation transport in the general-relativistic multigrid numerical (Gmunu) code. We compare the neutrino opacity spectra for all interactions and estimate their contributions at hydrodynamical points in core-collapse supernovae and binary neutron star (BNS) postmerger remnants, and predict the effects of improved opacities in comparison to conventional ones for a BNS postmerger at a specific hydrodynamical point. We test the implementation of the conventional set of interactions by comparing it to an open-source neutrino library NuLib in a core-collapse supernova simulation. We demonstrate good agreement with discrepancies of less than ∼10% in luminosity for all neutrino species, while also highlighting the reasons contributing to the differences. To compare the advanced interactions to the conventional set in core-collapse supernova modeling, we perform simulations to analyze their impacts on neutrino signatures, hydrodynamical behaviors, and shock dynamics, showing significant deviations.

A detailed measurement is made of the metallicity distributions, kinematics, and dynamics of the thin and thick disks across a large disk volume (5.0 ≤ R ≤ 15.0 kpc and ∣Z∣ ≤ 3.0 kpc) by using the LAMOST–APOGEE red clump stars. The metallicity distribution results show that the radial metallicity gradient Δ[Fe/H]/ΔR of the thin disk weakens with ∣Z∣ from −0.06 dex kpc⁻¹ at around ∣Z∣ < 0.25 kpc to −0.02 dex kpc⁻¹ at around ∣Z∣ > 2.75 kpc, while the thick disk displays a global weak positive Δ[Fe/H]/ΔR that is generally weaker than 0.01 dex kpc⁻¹. The vertical metallicity gradient Δ[Fe/H]/Δ∣Z∣ steadily weakened from −0.36 dex kpc⁻¹ at R ∼ 5.5 kpc to −0.05 dex kpc⁻¹ at around R > 11.5 kpc for the thin disk, while the thick disk presents an almost constant value (nearly −0.06∼−0.08 dex kpc⁻¹) for all the R bins. These results indicate the contribution of the radial migration to the disk evolution, and the obvious north–south asymmetry in [Fe/H] may be linked to disk warp and/or disk perturbation events. The oscillations in the corrected Δ[Fe/H]/Δ∣Z∣ with R likely arise from the resonances with the Galactic bar. Our detailed measurements of ΔV_ϕ/Δ[Fe/H] indicate an inside-out and upside-down star formation scenario for the thick disk. The results of eccentricity distributions and [α/Fe]–velocity dispersion relations are likely to suggest that thick-disk stars require an obvious contribution from other heating mechanisms, such as mergers and accretion, or are born in the chaotic mergers of gas-rich systems and/or the turbulent interstellar medium.

We present the updated galaxy cluster catalog of the second Planck catalog of Sunyaev–Zel'dovich sources (PSZ2) through the compilation of the data for clusters and galaxies with spectroscopically measured redshifts in the literature. The original version of PSZ2 comprises 1653 Sunyaev–Zel'dovich (SZ) sources, of which 1203 have been validated as genuine galaxy clusters, while the remaining 450 sources are yet to be validated. To increase the number of genuine clusters in PSZ2, we first update the validations of the cluster candidates and their redshift information using the data compiled for the confirmed clusters and the member galaxies in the literature. We then use the galaxy redshift data in the fields of the remaining cluster candidates by searching for possible member galaxies with measured spectroscopic redshifts around the SZ centroids. In this search process, we classify clusters as strong candidates if they contain more than nine galaxies within a 4500 km s⁻¹ velocity range and within 15' around the SZ centroids. This process results in the validation of 139 new genuine clusters, the update of redshift information on 399 clusters, and the identification of 10 strong candidates, which increases the number of validated clusters up to 1334 among the 1653 SZ sources. Our updated galaxy cluster catalog will be very useful for studies of galaxy formation and cosmology through a combination with other all-sky surveys including the Wide-field Infrared Survey Explorer and SPHEREx.