Papers for Thursday, Apr 15 2021

Nuclear star clusters (NSCs) are the densest stellar systems in the Universe and are found in the centres of all types of galaxies. They are thought to form via mergers of star clusters such as ancient globular clusters (GCs) that spiral to the centre as a result of dynamical friction or through in-situ star formation directly at the galaxy centre. There is evidence that both paths occur, but the relative contribution of either channel and their correlation with galaxy properties are not yet constrained observationally. We aim to derive the dominant NSC formation channel for a sample of 25 nucleated galaxies, mostly in the Fornax galaxy cluster, with stellar masses between $M_\rm{gal} \sim 10^8$ and $10^{10.5} M_\odot$ and NSC masses between $M_\rm{NSC} \sim 10^5$ and $10^{8.5} M_\odot$. Using Multi-Unit Spectroscopic Explorer (MUSE) data from the Fornax 3D survey and the ESO archive, we derive star formation histories, mean ages and metallicities of NSCs, and compare them to the host galaxies. In many low-mass galaxies, the NSCs are significantly more metal-poor than the hosts with properties similar to GCs. In contrast, in the massive galaxies, we find diverse star formation histories and cases of ongoing or recent in-situ star formation. Massive NSCs ($> 10^7 M_\odot$) occupy a different region in the mass-metallicity diagram than lower mass NSCs and GCs, indicating a different enrichment history. We find a clear transition of the dominant NSC formation channel with both galaxy and NSC mass. We hypothesise that while GC-accretion forms the NSCs of the dwarf galaxies, central star formation is responsible for the efficient mass build up in the most massive NSCs in our sample. At intermediate masses, both channels can contribute. The transition between these formation channels seems to occur at galaxy masses $M_\rm{gal} \sim 10^9 M_\odot$ and NSC masses $M_\rm{NSC} \sim 10^7 M_\odot$.

We present the open source few-body gravity integration toolkit {\tt SpaceHub}. {\tt SpaceHub} offers a variety of algorithmic methods, including the unique algorithms AR-Radau, AR-Sym6, AR-ABITS and AR-chain$^+$ which we show out-perform other methods in the literature and allow for fast, precise and accurate computations to deal with few-body problems ranging from interacting black holes to planetary dynamics. We show that AR-Sym6 and AR-chain$^+$, with algorithmic regularization, chain algorithm, active round-off error compensation and a symplectic kernel implementation, are the fastest and most accurate algorithms to treat black hole dynamics with extreme mass ratios, extreme eccentricities and very close encounters. AR-Radau, the first regularized Radau integrator with round off error control down to 64 bits floating point machine precision, has the ability to handle extremely eccentric orbits and close approaches in long-term integrations. AR-ABITS, a bit efficient arbitrary precision method, achieves any precision with the least CPU cost compared to other open source arbitrary precision few-body codes. With the implementation of deep numerical and code optimization, these new algorithms in {\tt SpaceHub} prove superior to other popular high precision few-body codes in terms of performance, accuracy and speed.

ASASSN-14ko is a recently discovered periodically flaring transient at the center of the AGN ESO 253-G003 with a slowly decreasing period. Here we show that the flares originate from the northern, brighter nucleus in this dual-AGN, post-merger system. The light curves for the two flares that occurred in May 2020 and September 2020 are nearly identical over all wavelengths. For both events, Swift observations showed that the UV and optical wavelengths brightened in unison. The effective temperature of the UV/optical emission rises and falls with the increase and subsequent decline in the luminosity. The X-ray flux, in contrast, first rapidly drops over $\sim$2.6 days, rises for $\sim$5.8 days, drops again over $\sim$4.3 days and then recovers. The X-ray spectral evolution of the two flares differ, however. During the May 2020 peak the spectrum softened with increases in the X-ray luminosity, while we observed the reverse for the September 2020 peak.

This is the second of a series of papers that focuses on searching large sets of photometric light curves for evidence of close binaries with a dormant black hole, and, in some cases, a dormant neutron star. The detection of such a binary is based on identifying a star that displays a large ellipsoidal periodic modulation, induced by tidal interaction with its companion. Based on the observed ellipsoidal amplitude and the primary mass and radius, one can derive a minimum mass ratio of the binary. A binary with a minimum mass ratio significantly larger than unity might be a candidate for having a dormant compact-object companion. Unfortunately, the photometric search is hampered by the fact that in many cases the primary mass and radius are not well known. In this paper we present a simple approach that circumvents this problem by suggesting a robust modified minimum mass ratio, assuming the primary fills its Roche lobe. The newly defined modified minimum mass ratio is always smaller than the minimum mass ratio, which is, in its turn, smaller than the actual mass ratio. Therefore, binaries with a modified minimum mass ratio larger than unity are candidates for having a compact-object secondary.

Hydrogen deuteride (HD) is prevalent in a wide variety of astrophysical environments, and measuring its large-scale distribution at different epochs can in principle provide information about the properties of these environments. In this paper, we explore the prospects for accessing this distribution using line intensity mapping of emission from the lowest rotational transition in HD, focusing on observations of the epoch of reionization ($z\sim6-10$) and earlier. We find the signal from the epoch of reionization to be strongest most promising, through cross-correlations within existing [CII] intensity mapping surveys. While the signal we predict is out of reach for current-generation projects, planned future improvements should be able to detect reionization-era HD without any additional observations, and would help to constrain the properties of the star-forming galaxies thought to play a key role in reionization. We also investigate several avenues for measuring HD during "cosmic dawn" ($z\sim10-30$), a period in which HD could provide one of the only complementary observables to 21$\,$cm intensity maps. We conclude that existing and planned facilities are poorly matched to the specifications desirable for a significant detection, though such a measurement may be achievable with sustained future effort. Finally, we explain why HD intensity mapping of the intergalactic medium during the cosmic dark ages ($z\gtrsim 30$) appears to be out of reach of any conceivable experiment.

Context. The different theoretical models concerning the formation of high-mass stars make distinct predictions regarding their progenitors, i.e. the high-mass prestellar cores. However, so far no conclusive observation of such objects has been made. Aims. We aim to study the very early stages of high-mass star formation in two infrared-dark, massive clumps, to identify the core population that they harbour. Methods. We obtained ALMA observations of continuum emission at 0.8mm and of the ortho-$\rm H_2D^+$ transition at 372GHz towards the two clumps. We use the SCIMES algorithm to identify cores in the position-position-velocity space, finding 16 cores. We model their observed spectra in the LTE approximation, deriving the centroid velocity, linewidth, and column density maps. We also study the correlation between the continuum and molecular data, which in general do not present the same structure. Results. We report for the first time the detection of ortho-$\rm H_2D^+$ in high-mass star-forming regions performed with an interferometer. The molecular emission shows narrow and subsonic lines, suggesting that locally the temperature of the gas is less than 10K. From the continuum emission we estimate the cores' total masses, and compare them with the respective virial masses. We also compute the volume density values, which are found to be higher than $10^{6}\, \rm cm^{-3}$. Conclusions. Our data confirm that ortho-$\rm H_2D^+$ is an ideal tracer of cold and dense gas. Interestingly, almost all the $\rm H_2D^+$-identified cores are less massive than 13M_sun , with the exception of one core in AG354. Furthermore, most of them are subvirial and larger than their Jeans masses. These results are difficult to explain in the context of the turbulent accretion models, which predict massive and virialised prestellar cores.

Stars form by accreting material from their surrounding disks. There is a consensus that matter flowing through the disk is channelled onto the stellar surface by the stellar magnetic field. This is thought to be strong enough to truncate the disk close to the so-called corotation radius where the disk rotates at the same rate as the star. Spectro-interferometric studies in young stellar objects show that Hydrogen is mostly emitted in a region of a few milliarcseconds across, usually located within the dust sublimation radius. Its origin is still a matter of debate and it can be interpreted as coming from the stellar magnetosphere, a rotating wind or a disk. In the case of intermediate-mass Herbig AeBe stars, the fact that the Br gamma emission is spatially resolved rules out that most of the emission comes from the magnetosphere. This is due to the weak magnetic fields (some tenths of G) detected in these sources, resulting in very compact magnetospheres. In the case of T Tauri sources, their larger magnetospheres should make them easier to resolve. However, the small angular size of the magnetosphere (a few tenths of milliarcseconds), along with the presence of winds emitting in Hydrogen make the observations interpretation challenging. Here, we present direct evidence of magnetospheric accretion by spatially resolving the inner disk of the 60 pc T Tauri star TW Hydrae through optical long baseline interferometry. We find that the hydrogen near-infrared emission comes from a region approximately 3.5 stellar radii (R*) across. This region is within the continuum dusty disk emitting region (Rcont = 7 R*) and smaller than the corotation radius which is twice as big. This indicates that the hydrogen emission originates at the accretion columns, as expected in magnetospheric accretion models, rather than in a wind emitted at much larger distance (>1au).

The modern search for extraterrestrial intelligence (SETI) began with the seminal publications of Cocconi & Morrison (1959) and Schwartz & Townes (1961), who proposed to search for narrow-band signals in the radio spectrum, and for optical laser pulses. Over the last six decades, more than one hundred dedicated search programs have targeted these wavelengths; all with null results. All of these campaigns searched for classical communications, that is, for a significant number of photons above a noise threshold; with the assumption of a pattern encoded in time and/or frequency space. I argue that future searches should also target quantum communications. They are preferred over classical communications with regards to security and information efficiency, and they would have escaped detection in all previous searches. The measurement of Fock state photons or squeezed light would indicate the artificiality of a signal. I show that quantum coherence is feasible over interstellar distances, and explain for the first time how astronomers can search for quantum transmissions sent by ETI to Earth, using commercially available telescopes and receiver equipment.

We present observations of three protoplanetary disks in visible scattered light around M-type stars in the Upper Scorpius OB association using the STIS instrument on the Hubble Space Telescope. The disks around stars 2MASS J16090075-1908526, 2MASS J16142029-1906481 and 2MASS J16123916-1859284 have all been previously detected with ALMA, and 2MASS J16123916-1859284 has never previously been imaged at scattered light wavelengths. We process our images using Reference Differential Imaging, comparing and contrasting three reduction techniques - classical subtraction, Karhunen-Loeve Image Projection and Non-Negative Matrix Factorisation, selecting the classical method as the most reliable of the three for our observations. Of the three disks, two are tentatively detected (2MASS J16142029-1906481 and 2MASS J16123916-1859284), with the third going undetected. Our two detections are shown to be consistent when varying the reference star or reduction method used, and both detections exhibit structure out to projected distances of > 200 au. Structures at these distances from the host star have never been previously detected at any wavelength for either disk, illustrating the utility of visible-wavelength observations in probing the distribution of small dust grains at large angular separations.

We present detailed characterization of the extremely dusty main sequence star TYC 8830 410 1. This system hosts inner planetary system dust (Tdust~300 K) with a fractional infrared luminosity of ~1%. Mid-infrared spectroscopy reveals a strong, mildy-crystalline solid-state emission feature. TYC 8830 410 1 (spectral type G9V) has a 49.5" separation M4-type companion co-moving and co-distant with it, and we estimate a system age of ~600 Myr. TYC 8830 410 1 also experiences "dipper"-like dimming events as detected by ASAS-SN, TESS, and characterized in more detail with the LCOGT. These recurring eclipses suggest at least one roughly star-sized cloud of dust orbits the star in addition to assorted smaller dust structures. The extreme properties of the material orbiting TYC 8830 410 1 point to dramatic dust-production mechanisms that likely included something similar to the giant-impact event thought to have formed the Earth-Moon system, although hundreds of millions of years after such processes are thought to have concluded in the solar system. TYC 8830 410 1 holds promise to deliver significant advances in our understanding of the origin, structure, and evolution of extremely dusty inner planetary systems.

The evolution of a relativistic blastwave is usually delineated under the assumption of pressure balance between forward- and reverse-shocked regions. However, such a treatment usually violates the energy conservation law. A mechanical model of non-magnetized blastwaves was proposed in previous works to solve the problem. In this paper, we generalize the mechanical model to the case of a blastwave driven by an ejecta with an arbitrary magnetization parameter $\sigma_{\rm ej}$. We test our modified mechanical model by considering a long-lasting magnetized ejecta and found that it is much better in terms of energy conservation than the pressure-balance treatment. For a constant luminosity $L_{\rm ej} = 10^{47}{\rm erg~s^{-1}}$ and $\sigma_{\rm ej} < 10$, the deviation from energy conservation is negligibly small at small radii, but only reaches less than $25\%$ even at $10^{19}{\rm cm}$ from the central engine. Assuming a finite lifetime of the central engine, the reverse shock crosses the magnetized ejecta earlier for the ejecta wiith a higher $\sigma_{\rm ej}$. After shock crossing, all the ejecta energy is injected into the blastwave.

We study the properties of the faint X-ray activity of Galactic transient black hole candidate XTE~J1908+094 during its 2019 outburst. Here, we report the results of detailed spectral and temporal analysis during this outburst using observations from {\it Nuclear Spectroscopic Telescope Array (NuSTAR)}. We have not observed any quasi-periodic-oscillations (QPOs) in the power density spectrum (PDS). The spectral study suggests that the source remained in the softer (more precisely in soft-intermediate) spectral state during this short period of the X-ray activity. We notice a faint but broad Fe K$\alpha$ emission line at around 6.5 keV. We also estimate the probable mass of the black hole to be $6.5^{+0.5}_{-0.7}~M_\odot$ with 90\% confidence.

We performed a detailed analysis of the detectability of a wide range of gravitational waves derived from core-collapse supernova simulations using gravitational-wave detector noise scaled to the sensitivity of the upcoming fourth and fifth observing runs of the Advanced LIGO, Advanced Virgo, and KAGRA. We use the coherent WaveBurst algorithm, which was used in the previous observing runs to search for gravitational waves from core-collapse supernovae. As coherent WaveBurst makes minimal assumptions on the morphology of a gravitational-wave signal, it can play an important role in the first detection of gravitational waves from an event in the Milky Way. We predict that signals from neutrino-driven explosions could be detected up to an average distance of 10 kpc, and distances of over 100 kpc can be reached for explosions of rapidly rotating progenitor stars. An estimated minimum signal-to-noise ratio of 10-25 is needed for the signals to be detected. We quantify the accuracy of the waveforms reconstructed with coherent WaveBurst and we determine that the most challenging signals to reconstruct are those produced in long-duration neutrino-driven explosions and models that form black holes a few seconds after the core bounce.

Oxygen is a promising exoplanet biosignature due to the evolutionary advantage conferred by harnessing starlight for photosynthesis, and the apparent low likelihood of maintaining oxygen-rich atmospheres without life. Hypothetical scenarios have been proposed for non-biological oxygen accumulation on planets around late M-dwarfs, where the extended pre-main sequence may favor abiotic O$_2$ accumulation. In contrast, abiotic oxygen accumulation on planets around F, G, and K-type stars is seemingly less likely, provided they possess substantial non-condensable gas inventories. The comparative robustness of oxygen biosignatures around larger stars has motivated plans for next-generation telescopes capable of oxygen detection on planets around sun-like stars. However, the general tendency of terrestrial planets to develop oxygen-rich atmospheres across a broad range of initial conditions and evolutionary scenarios has not been explored. Here, we use a coupled thermal-geochemical-climate model of terrestrial planet evolution to illustrate three scenarios whereby significant abiotic oxygen can accumulate around sun-like stars, even when significant non-condensable gas inventories are present. For Earth-mass planets, we find abiotic oxygen can accumulate to modern levels if (1) the CO$_2$:H$_2$O ratio of the initial volatile inventory is high, (2) the initial water inventory exceeds ~50 Earth oceans, or (3) the initial water inventory is very low. Fortunately, these three abiotic oxygen scenarios could be distinguished from biological oxygen with observations of other atmospheric constituents or characterizing the planetary surface. This highlights the need for broadly capable next-generation telescopes that are equipped to constrain surface water inventories via time-resolved photometry and search for temporal biosignatures or disequilibrium biosignatures to assess whether oxygen is biogenic.

Detached circumplanetary disks are unstable to tilting as a result of the stellar tidal potential. We examine how a tilted circumplanetary disk affects the evolution of the spin axis of an oblate planet. The disk is evolved using time-dependent equations for linear wave-like warp evolution, including terms representing the effect of the tidal potential and planetary oblateness. For a disk with a sufficiently large mass, we find that the planet spin quickly aligns to the misaligned disk. The tilt of the planetary spin axis then increases on the same timescale as the disk. This can be an efficient mechanism for generating primordial obliquity in giant planets. We suggest that directly imaged exoplanets at large orbital radii, where the disk mass criterion is more likely to be satisfied, could have significant obliquities due to the tilt instability of their circumplanetary disks.

We use the known surface boulder-size distribution of the C-type rubble pile asteroid Ryugu (NEA 162173) to determine its macroporosity, assuming it is a homogeneous granular aggregate. We show that the volume-frequency distribution of its boulders, cobbles and pebbles, is well represented by a lognormal function with $\sigma = 2.4 \pm 0.1$ and $\mu = 0.2 \pm 0.05$. Application of linear-mixture packing theory yields a value for the macroporosity of $\phi = 0.14 \pm 0.04$. Given its low bulk density of 1.19 gm cm$^{-3}$, this implies an average density for Ryugu's rocks of $1.38 \pm 0.07$ gm cm$^{-3}$ throughout its volume, consistent with a recent determination for surface boulders based on their thermal properties. This supports the spectrum-based argument that IDP's may be the best analog material available on Earth and suggests that high-density, well-lithified objects such as chondrules and chondrule-bearing chondrites may be rare on Ryugu. Implications of this result for the origin of chondrules, a long-standing problem in cosmochemistry, are discussed. We propose that chondrules and most chondrites formed together in rare lithification events, which occurred during the accretion of chondritic envelopes to large, differentiated planetesimals at a time when they were still hot from $^{26}$Al decay.

Star Formation Histories (SFHs) reveal physical processes that influence how galaxies form their stellar mass. We compare the SFHs of a sample of 36 nearby (D $\leq$ 4 Mpc) dwarf galaxies from the ACS Nearby Galaxy Survey Treasury (ANGST), inferred from the Color Magnitude Diagrams (CMDs) of individually resolved stars in these galaxies, with those reconstructed by broad-band Spectral Energy Distribution (SED) fitting using the Dense Basis SED fitting code. When comparing individual SFHs, we introduce metrics for evaluating SFH reconstruction techniques. For both the SED and CMD methods, the median normalized SFH of galaxies in the sample shows a period of quiescence at lookback times of 3-6 Gyr followed by rejuvenated star formation over the past 3 Gyr that remains active until the present day. To determine if these represent special epochs of star formation in the D $\leq$ 4 Mpc portion of the Local Volume, we break this ANGST dwarf galaxy sample into subsets based on specific star formation rate and spatial location. Modulo offsets between the methods of about 1 Gyr, all subsets show significant decreases and increases in their median normalized SFHs at the same epochs, and the majority of the individual galaxy SFHs are consistent with these trends. These results motivate further study of potential synchronized star formation quiescence and rejuvenation in the Local Volume as well as development of a hybrid method of SFH reconstruction that combines CMDs and SEDs, which have complementary systematics.

The advent of increasingly large and complex datasets has fundamentally altered the way that scientists conduct astronomy research. The need to work closely to the data has motivated the creation of online science platforms, which include a suite of software tools and services, therefore going beyond data storage and data access. We present two example applications of Jupyter as a part of astrophysical science platforms for professional researchers and students. First, the Astro Data Lab is developed and operated by NOIRLab with a mission to serve the astronomy community with now over 1500 registered users. Second, the Dark Energy Spectroscopic Instrument science platform serves its geographically distributed team comprising about 900 collaborators from over 90 institutions. We describe the main uses of Jupyter and the interfaces that needed to be created to embed it within science platform ecosystems. We use these examples to illustrate the broader concept of empowering researchers and providing them with access to not only large datasets but also cutting-edge software, tools, and data services without requiring any local installation, which can be relevant for a wide range of disciplines. Future advances may involve science platform networks, and tools for simultaneously developing Jupyter notebooks to facilitate collaborations.

Light curves produced by the Kepler mission demonstrate stochastic brightness fluctuations (or "flicker") of stellar origin which contribute to the noise floor, limiting the sensitivity of exoplanet detection and characterization methods. In stars with surface convection, the primary driver of these variations on short (sub-eight-hour) timescales is believed to be convective granulation. In this work, we improve existing models of this granular flicker amplitude, or $F_8$, by including the effect of the Kepler bandpass on measured flicker, by incorporating metallicity in determining convective Mach numbers, and by using scaling relations from a wider set of numerical simulations. To motivate and validate these changes, we use a recent database of convective flicker measurements in Kepler stars, which allows us to more fully detail the remaining model--prediction error. Our model improvements reduce the typical misprediction of flicker amplitude from a factor of 2.5 to 2. We rule out rotation period and strong magnetic activity as possible explanations for the remaining model error, and we show that binary companions may affect convective flicker. We also introduce an "envelope" model which predicts a range of flicker amplitudes for any one star to account for some of the spread in numerical simulations, and we find that this range covers 78% of observed stars. We note that the solar granular flicker amplitude is lower than most Sun-like stars. This improved model of convective flicker amplitude can better characterize this source of noise in exoplanet studies as well as better inform models and simulations of stellar granulation.

We present the serendipitous discovery of an extremely broad ($\Delta V_{LSR} \sim 150$ km/s), faint ($T_{mb} < 10 \textrm{mK}$), and ubiquitous 1667 and 1665 MHz ground-state thermal OH emission towards the 2nd quadrant of the outer Galaxy ($R_{gal}$ > 8 kpc) with the Green Bank Telescope. Originally discovered in 2015, we describe the redundant experimental, observational, and data quality tests of this result over the last five years. The longitude-velocity distribution of the emission unambiguously suggests large-scale Galactic structure. We observe a smooth distribution of OH in radial velocity that is morphologically similar to the HI radial velocity distribution in the outer Galaxy, showing that molecular gas is significantly more extended in the outer Galaxy than previously expected. Our results imply the existence of a thick ($-200< z < 200$ pc) disk of diffuse ($n_{H_{2}}$ $\sim$ 5 $\times$ 10$^{-3}$ cm$^{-3}$) molecular gas in the Outer Galaxy previously undetected in all-sky CO surveys.

We search for microlensing planets with signals exhibiting no caustic-crossing features, considering the possibility that such signals may be missed due to their weak and featureless nature. For this purpose, we reexamine the lensing events found by the KMTNet survey before the 2019 season. From this investigation, we find two new planetary lensing events, KMT-2018-BLG-1976 and KMT-2018-BLG-1996. We also present the analysis of the planetary event OGLE-2019-BLG-0954, for which the planetary signal was known, but no detailed analysis has been presented before. We identify the genuineness of the planetary signals by checking various interpretations that can generate short-term anomalies in lensing light curves. From Bayesian analyses conducted with the constraint from available observables, we find that the host and planet masses are $(M_1, M_2)\sim (0.65~M_\odot, 2~M_{\rm J})$ for KMT-2018-BLG-1976L, $\sim (0.69~M_\odot, 1~M_{\rm J})$ for KMT-2018-BLG-1996L, and $\sim (0.80~M_\odot, 14~M_{\rm J})$ for OGLE-2019-BLG-0954L. The estimated distance to OGLE-2019-BLG-0954L, $3.63^{+1.22}_{-1.64}$~kpc, indicates that it is located in the disk, and the brightness expected from the mass and distance matches well the brightness of the blend, indicating that the lens accounts for most of the blended flux. The lens of OGLE-2019-BLG-0954 could be resolved from the source by conducting high-resolution follow-up observations in and after 2024.

We present spatially resolved ($0.1'' - 1.0''$) radio maps of Neptune taken from the Very Large Array and Atacama Large Submillimeter/Millimeter Array between $2015-2017$. Combined, these observations probe from just below the main methane cloud deck at $\sim 1$ bar down to the NH$_4$SH cloud at $\sim50$ bar. Prominent latitudinal variations in the brightness temperature are seen across the disk. Depending on wavelength, the south polar region is $5-40$ K brighter than the mid-latitudes and northern equatorial region. We use radiative transfer modeling coupled to Markov Chain Monte Carlo methods to retrieve H$_2$S, NH$_3$, and CH$_4$ abundance profiles across the disk, though only strong constraints can be made for H$_2$S. Below all cloud formation, the data are well fit by $53.8^{+18.9}_{-13.4}\times$ and $3.9^{+2.1}_{-3.1}\times$ protosolar enrichment in the H$_2$S and NH$_3$ abundances, respectively, assuming a dry adiabat. Models in which the radio-cold mid-latitudes and northern equatorial region are supersaturated in H$_2$S are statistically favored over models following strict thermochemical equilibrium. H$_2$S is more abundant at the equatorial region than at the poles, indicative of strong, persistent global circulation. Our results imply that Neptune's sulfur-to-nitrogen ratio exceeds unity as H$_2$S is more abundant than NH$_3$ in every retrieval. The absence of NH$_3$ above 50 bar can be explained either by partial dissolution of NH$_3$ in an ionic ocean at GPa pressures or by a planet formation scenario in which hydrated clathrates preferentially delivered sulfur rather than nitrogen onto planetesimals, or a combination of these hypotheses.

Cosmological probes based on galaxy clusters rely on cluster number counts and large-scale structure information. X-ray cluster surveys are well suited for this purpose, since they are far less affected than optical surveys by projection effects, and cluster properties can be predicted with good accuracy. The XMM Cluster Archive Super Survey, X-CLASS, is a serendipitous search of X-ray-detected galaxy clusters in 4176 XMM-Newton archival observations until August 2015. All observations are clipped to exposure times of 10 and 20 ks to obtain uniformity and they span ~269 deg$^2$ across the high-Galactic latitude sky ($|b|> 20^o$). The main goal of the survey is the compilation of a well-selected cluster sample suitable for cosmological analyses. We describe the detection algorithm, the visual inspection, the verification process and the redshift validation of the cluster sample, as well as the cluster selection function computed by simulations. We also present the various metadata that are released with the catalogue, along with the redshifts of 124 clusters obtained with a dedicated multi-object spectroscopic follow-up programme. With this publication we release the new X-CLASS catalogue of 1646 well-selected X-ray-detected clusters over a wide sky area, along with their selection function. The sample spans a wide redshift range, from the local Universe up to z~1.5, with 982 spectroscopically confirmed clusters, and over 70 clusters above z=0.8. Because of its homogeneous selection and thorough verification, the cluster sample can be used for cosmological analyses, but also as a test-bed for the upcoming eROSITA observations and other current and future large-area cluster surveys. It is the first time that such a catalogue is made available to the community via an interactive database which gives access to a wealth of supplementary information, images, and data.

Because of its distinctive compositional properties and variability, low-speed ($\lesssim 450$ km s$^{-1}$) solar wind is widely believed to originate from coronal streamers, unlike high-speed wind, which comes from coronal holes. An alternative scenario is that the bulk of the slow wind (excluding that in the immediate vicinity of the heliospheric current sheet) originates from rapidly diverging flux tubes rooted inside small coronal holes or just within the boundaries of large holes. This viewpoint is based largely on photospheric field extrapolations, which are subject to considerable uncertainties and do not include dynamical effects, making it difficult to be certain whether a source is located just inside or outside a hole boundary, or whether a high-latitude hole will be connected to Earth. To minimize the dependence on field-line extrapolations, we have searched for cases where equatorial coronal holes at central meridian are followed by low-speed streams at Earth. We describe 14 examples from the period 2014--2017, involving Fe XIV 21.1 nm coronal holes located near active regions and having equatorial widths of $\sim$3$^\circ$--10$^\circ$. The associated in situ wind was characterized by speeds $v\sim 300$--450 km s$^{-1}$ and by O$^{7+}$/O$^{6+}$ ratios of $\sim$0.05--0.15, with $v$ showing the usual correlation with proton temperature. In addition, consistent with other recent studies, this slow wind had remarkably high Alfv\'enicity, similar to that in high-speed streams. We conclude that small coronal holes are a major contributor to the slow solar wind during the maximum and early post-maximum phases of the solar cycle.

Earlier studies using extreme-ultraviolet images and line-of-sight magnetograms from the Solar Dynamics Observatory (SDO) have suggested that active region (AR) plages and strong network concentrations often have small, looplike features embedded within them, even though no minority-polarity flux is visible in the corresponding magnetograms. Because of the unexpected nature of these findings, we have searched the SDO database for examples of inverted-Y structures rooted inside "unipolar" plages, with such jetlike structures being interpreted as evidence for magnetic reconnection between small bipoles and the dominant-polarity field. Several illustrative cases are presented from the period 2013--2015, all of which are associated with transient outflows from AR "moss." The triangular or dome-shaped bases have horizontal dimensions of $\sim$2--4 Mm, corresponding to $\sim$1--3 granular diameters. We also note that the spongy-textured Fe IX 17.1 nm moss is not confined to plages, but may extend into regions where the photospheric field is relatively weak or even has mixed polarity. We again find a tendency for bright coronal loops seen in the 17.1, 19.3, and 21.1 nm passbands to show looplike fine structure and compact brightenings at their footpoints. These observations provide further confirmation that present-day magnetograms are significantly underrepresenting the amount of minority-polarity flux inside AR plages and again suggest that footpoint reconnection and small-scale flux cancellation may play a major role in coronal heating, both inside and outside ARs.

Solar active region 12673 produced two successive X-class flares (X2.2 and X9.3) approximately 3 hours apart in September 2017. The X9.3 was the recorded largest solar flare in Solar Cycle 24. In this study we perform a data-constrained magnetohydrodynamic simulation taking into account the observed photospheric magnetic field to reveal the initiation and dynamics of the X2.2 and X9.3 flares. According to our simulation, the X2.2 flare is first triggered by magnetic reconnection at a local site where at the photosphere the negative polarity intrudes into the opposite-polarity region. This magnetic reconnection expels the innermost field lines upward beneath which the magnetic flux rope is formed through continuous reconnection with external twisted field lines. Continuous magnetic reconnection after the X2.2 flare enhances the magnetic flux rope, which is lifted up and eventually erupts via the torus instability. This gives rise to the X9.3 flare.

Pollution of white dwarf atmospheres may be caused by asteroids that originate from the locations of secular and mean-motion resonances in planetary systems. Asteroids in these locations experience increased eccentricity, leading to tidal disruption by the white dwarf. We examine how the $\nu_6$ secular resonance shifts outwards into a previously stable region of the asteroid belt, as the star evolves to a white dwarf. Analytic secular models require a planet to be engulfed in order to shift the resonance. We show with numerical simulations that as a planet gets engulfed by the evolving star, the secular resonance shifts and the rate of tidal disruption events increases with the engulfed planet's mass and its orbital separation. We also investigate the behaviour of mean-motion resonances. The width of a mean-motion resonance increases as the star loses mass and becomes a white dwarf. The $\nu_6$ secular resonance is more efficient at driving tidal disruptions than mean-motion resonances with Jupiter. By examining 230 observed exoplanetary systems whose central star will evolve into a white dwarf, we find that along with an Earth mass planet at $1\,\rm au$, hot Jupiters at a semi--major axis $a\gtrsim 0.05\,\rm au$ and super--Earths of mass $10\,\rm M_\oplus$ at $a\gtrsim 0.3\,\rm au$ represent planet types whose engulfment shifts resonances enough to cause pollution of the white dwarfs to a degree in agreement with observations.

We investigated fundamental processes of collisional sticking and fragmentation of dust aggregates by carrying out N-body simulations of submicron-sized icy dust monomers. We examined the condition for collisional growth of two colliding dust aggregates in a wide range of the mass ratio, 1-64. We found that the mass transfer from a larger dust aggregate to a smaller one is a dominant process in collisions with a mass ratio of 2-30 and impact velocity of \approx 30-170 m s^-1. As a result, the critical velocity, v_fra, for fragmentation of the largest body is considerably reduced for such unequal-mass collisions; v_fra of collisions with a mass ratio of 3 is about half of that obtained from equal-mass collisions. The impact velocity is generally higher for collisions between dust aggregates with higher mass ratios because of the difference between the radial drift velocities in the typical condition of protoplanetary disks. Therefore, the reduced v_fra for unequal-mass collisions would delay growth of dust grains in the inner region of protoplanetary disks.

We use the data of modern digital sky surveys (PanSTARRS-1, SDSS) combined with HI-line and far ultraviolet (GALEX) surveys to reclassify 165 early-type galaxies from the Catalog of Isolated Galaxies (KIG). As a result, the number of E- and S0-type galaxies reduced to 91. Our search for companions of early-type KIG galaxies revealed 90 companions around 45 host galaxies with line-of-sight velocity differences $|dV| < 500$ km s$^{-1}$ and linear projected separations $R_{p} < 750$ kpc. We found no appreciable differences in either integrated luminosity or color of galaxies associated with the presence or absence of close neighbors. We found a characteristic orbital mass-to-luminosity ratio for 26 systems "KIG galaxy--companion" to be $M_{\odot}/L_{K} = (74\pm26) M_{\odot}/L_{\odot}$, which is consistent with the $M_{\rm orb}/L_{K}$ estimates for early-type isolated galaxies in the 2MIG catalog ($63 M_{\odot}/L_{\odot}$), and also with the $M_{\rm orb}/L_{K}$ estimates for E- and S0-type galaxies in the Local Volume: $38\pm22$ (NGC 3115), $82\pm26$ (NGC 5128), $65\pm20$ (NGC 4594). The high halo-to-stellar mass ratio for E- and S0-type galaxies compared to the average $(20\pm3) M_{\odot}/L_{\odot}$ ratio for bulgeless spiral galaxies is indicative of a significant difference between the dynamic evolution of early- and late-type galaxies.

The thermodynamic structure of protoplanetary discs is determined by the dust opacities which depend on the size of the dust grains and their chemical composition. In the inner regions, the grain sizes are regulated by the level of turbulence (e.g., $\alpha$-viscosity) and by the dust fragmentation velocity that represents the maximal velocity grains can have at a collision before they fragment. Here, we perform self-consistently calculated 2D hydrodynamical simulations considering a full grain size distribution of dust grains with a transition in the dust fragmentation velocity at the water-iceline, reflecting the results of previous laboratory experiments of particle collisions, where silicate particles normally have a lower dust fragmentation velocity than water-ice particles. Furthermore, we probe the effects of variations in the water abundance, the dust-to-gas ratio, and the turbulence parameter on the disc structure. For the discs with a transition in the dust fragmentation velocity at the water-iceline, we find a narrow but striking zone of planetary outward migration, also for low viscosity. In addition, we find a bump in the radial pressure gradient profile that tends to be located slightly interior to the iceline. Both of these features are present for all tested disc parameters. Thus, we conclude that the iceline can function both as a migration trap which can extend the growth times of planets before they migrate to the inner edge of the protoplanetary disc, and as a pressure trap where planetesimal formation might be initiated or enhanced.

The current role of data-driven science is constantly increasing its importance within Astrophysics, due to the huge amount of multi-wavelength data collected every day, characterized by complex and high-volume information requiring efficient and as much as possible automated exploration tools. Furthermore, to accomplish main and legacy science objectives of future or incoming large and deep survey projects, such as JWST, LSST and Euclid, a crucial role is played by an accurate estimation of photometric redshifts, whose knowledge would permit the detection and analysis of extended and peculiar sources by disentangling low-z from high-z sources and would contribute to solve the modern cosmological discrepancies. The recent photometric redshift data challenges, organized within several survey projects, like LSST and Euclid, pushed the exploitation of multi-wavelength and multi-dimensional data observed or ad hoc simulated to improve and optimize the photometric redshifts prediction and statistical characterization based on both SED template fitting and machine learning methodologies. But they also provided a new impetus in the investigation on hybrid and deep learning techniques, aimed at conjugating the positive peculiarities of different methodologies, thus optimizing the estimation accuracy and maximizing the photometric range coverage, particularly important in the high-z regime, where the spectroscopic ground truth is poorly available. In such a context we summarize what learned and proposed in more than a decade of research.

In this work we use seeing-limited SINFONI H- and K-band spectroscopy to analyse the spatially resolved kinematics of the stellar component in the inner r<1-2 kpc of ten nearby (mean z=0.014) LIRGs. We compare the stellar kinematics with those for different gas phases, and analyse the relative effects of using different tracers when estimating dynamical masses. The stellar velocity and velocity dispersion maps are extracted in both near-IR bands by fitting the continuum emission using pPXF while we use the gas kinematics from previous works. We find that the different gas phases have similar kinematics, whereas the stellar component is rotating with slightly lower velocities (V$_*$~0.8V$_g$) but in significantly warmer orbits ($\sigma_*$~2$\sigma_g$) than the gas phases. These values indicate that stars are rotating in thick discs while the gas phases are confined in dynamically cooler rotating discs. However, these differences do not lead to significant discrepancies between the dynamical mass estimations based on the stellar and gas kinematics. This result suggests that the gas kinematics can be used to estimate M$_{dyn}$ also in z~2 SFGs, a galaxy population that shares many structural and kinematic properties with local LIRGs.

Recent observations of X-ray pulsars at low luminosities allow, for the first time, to compare theoretical models for the emission from highly magnetized neutron star atmospheres at low mass accretion rates ($\dot{M} \lesssim 10^{15}$ g s$^{-1}$) with the broadband X-ray data. The purpose of this paper is to investigate the spectral formation in the neutron star atmosphere at low $\dot{M}$ and to conduct a parameter study of physical properties of the emitting region. We obtain the structure of the static atmosphere, assuming that Coulomb collisions are the dominant deceleration process. The upper part of the atmosphere is strongly heated by the braking plasma, reaching temperatures of 30-40 keV, while its denser isothermal interior is much cooler (~2 keV). We numerically solve the polarized radiative transfer in the atmosphere with magnetic Compton scattering, free-free processes, and non-thermal cyclotron emission due to possible collisional excitations of electrons. The strongly polarized emitted spectrum has a double-hump shape that is observed in low-luminosity X-ray pulsars. A low-energy "thermal" component is dominated by extraordinary photons that can leave the atmosphere from deeper layers due to their long mean free path at soft energies. We find that a high-energy component is formed due to resonant Comptonization in the heated non-isothermal part of the atmosphere even in the absence of collisional excitations. The latter, however, affect the ratio of the two components. A strong cyclotron line originates from the optically thin, uppermost zone. A fit of the model to NuSTAR and Swift/XRT observations of GX 304-1 provides an accurate description of the data with reasonable parameters. The model can thus reproduce the characteristic double-hump spectrum observed in low-luminosity X-ray pulsars and provides insights into spectral formation.

In protoplanetary discs and planetary atmospheres, dust grains coagulate to form fractal dust aggregates. The geometric cross-section of these aggregates is a crucial parameter characterizing aerodynamical friction, collision rates, and opacities. However, numerical measurements of the cross-section are often time-consuming as aggregates exhibit complex shapes. In this study, we derive a novel analytic expression for geometric cross-sections of fractal aggregates. If an aggregate consists of $N$ monomers of radius $R_0$, its geometric cross-section $G$ is expressed as \begin{equation} \frac{G}{N\pi R_0^2}=\frac{A}{1+(N-1)\tilde{\sigma}}, \nonumber \end{equation} where $\tilde{\sigma}$ is an overlapping efficiency, and $A$ is a numerical factor connecting the analytic expression to the small non-fractal cluster limit. The overlapping efficiency depends on the fractal dimension, fractal prefactor, and $N$ of the aggregate, and its analytic expression is derived as well. The analytic expressions successfully reproduce numerically measured cross-sections of aggregates. We also find that our expressions are compatible with the mean-field light scattering theory of aggregates in the geometrical optics limit. The analytic expressions greatly simplify an otherwise tedious calculation and will be useful in model calculations of fractal grain growth in protoplanetary discs and planetary atmospheres.

We compare spectra of the zonal harmonics of the large-scale magnetic field of the Sun using observation results and solar dynamo models. The main solar activity cycle as recorded in these tracers is a much more complicated phenomenon than the eigen solution of solar dynamo equations with the growth saturated by a back reaction of the dynamo-driven magnetic field on solar hydrodynamics. The nominal 11(22)-year cycle as recorded in each mode has a specific phase shift varying from cycle to cycle; the actual length of the cycle varies from one cycle to another and from tracer to tracer. Both the observation and the dynamo model show an exceptional role of the axisymmetric $\ell_{5}$ mode. Its origin seems to be readily connected with the formation and evolution of sunspots on the solar surface. The results of observations and dynamo models show a good agreement for the low $\ell_{1}$ and $\ell_{3}$ modes. The results for these modes do not differ significantly for the axisymmetric and nonaxisymmetric models. Our findings support the idea that the sources of the solar dynamo arise as a result of both the distributed dynamo processes in the bulk of the convection zone and the surface magnetic activity.

To investigate the effects of massive star evolution on surrounding molecules, we select 9 massive clumps previously observed with the Atacama Pathfinder Experiment (APEX) telescope and the Submillimeter Array (SMA) telescope. Based on the observations of APEX, we obtain luminosity to mass ratio L$_{\rm clump}$/M$_{\rm clump}$ that range from 10 to 154 L$_{\sun}$/M$_{\sun}$, where some of them embedded Ultra Compact (UC) H\,{\footnotesize II} region. Using the SMA, CH$_3$CN (12$_{\rm K}$--11$_{\rm K}$) transitions were observed toward 9 massive star-forming regions. We derive the CH$_3$CN rotational temperature and column density using XCLASS program, and calculate its fractional abundance. We find that CH$_3$CN temperature seems to increase with the increase of L$_{\rm clump}$/M$_{\rm clump}$ when the ratio is between 10 to 40 L$_{\sun}$/M$_{\sun}$, then decrease when L$_{\rm clump}$/M$_{\rm clump}$ $\ge$ 40 L$_{\sun}$/M$_{\sun}$. Assuming the CH$_3$CN gas is heated by radiation from the central star, the effective distance of CH$_3$CN relative to the central star is estimated. The distance range from $\sim$ 0.003 to $\sim$ 0.083 pc, which accounts for $\sim$ 1/100 to $\sim$ 1/1000 of clump size. The effective distance increases slightly as L$_{\rm clump}$/M$_{\rm clump}$ increases (R$_{\rm eff}$ $\sim$ (L$_{\rm clump}$/M$_{\rm clump}$)$^{0.5\pm0.2}$). Overall, the CH$_3$CN abundance is found to decrease as the clumps evolve, e.g., X$_{\rm CH_3CN}$ $\sim$ (L$_{\rm clump}$/M$_{\rm clump}$)$^{-1.0\pm0.7}$. The steady decline of CH$_3$CN abundance as the clumps evolution can be interpreted as a result of photodissociation.

The mergers of double helium white dwarfs are believed to form isolated helium-rich hot subdwarfs. Observation shows that the helium-rich hot subdwarfs can be divided into two subgroups based on whether the surface is carbon-rich or carbon-normal. But it is not clear whether this distribution directly comes from binary evolution. We adopt the binary population synthesis (BPS) to obtain the population of single helium-rich hot subdwarfs according to the channel of double helium white dwarfs merger. We find that the merger channel can represent the two subgroups in the $T_{\rm{eff}}-\log g$ plane related to different masses of progenitor helium white dwarfs. For $Z$ = 0.02, the birth rates and local density of helium-rich hot subdwarf stars by the mergers of two helium white dwarfs is $\sim 4.82 \times 10^{-3}$ $\rm yr^{-1}$ and $\sim$ 290.0 $\rm kpc^{-3}$ at 13.7 Gyr in our Galaxy, respectively. The proportion of carbon-rich and carbon-normal helium-rich hot subdwarfs are 32$\%$ and 68$\%$, respectively.

Gaps in protoplanetary disks have long been hailed as signposts of planet formation. However, a direct link between exoplanets and disks remains hard to identify. We present a large sample study of ALMA disk surveys of nearby star-forming regions to disentangle this connection. All disks are classified as either structured (transition, ring, extended) or non-structured (compact) disks. Although low-resolution observations may not identify large scale substructure, we assume that an extended disk must contain substructure from a dust evolution argument. A comparison across ages reveals that structured disks retain high dust masses up to at least 10 Myr, whereas the dust mass of compact, non-structured disks decreases over time. This can be understood if the dust mass evolves primarily by radial drift, unless drift is prevented by pressure bumps. We identify a stellar mass dependence of the fraction of structured disks. We propose a scenario linking this dependence with that of giant exoplanet occurrence rates. We show that there are enough exoplanets to account for the observed disk structures if transitional disks are created by exoplanets more massive than Jupiter, and ring disks by exoplanets more massive than Neptune, under the assumption that most of those planets eventually migrate inwards. On the other hand, the known anti-correlation between transiting super-Earths and stellar mass implies those planets must form in the disks without observed structure, consistent with formation through pebble accretion in drift-dominated disks. These findings support an evolutionary scenario where the early formation of giant planets determines the disk's dust evolution and its observational appearance.

Within just two years, two interstellar objects (ISOs) - Oumuamuas and Borisov - have been discovered. Large quantities of planetesimals form as a by-product of planet formation. Therefore, it seems likely that ISOs are former planetesimals that became unbound from their parent star. The discoveries raise the question of the dominant ISO formation process. Here, we concentrate on planetesimals released during another star's close flybys. We quantify the amount of planetesimals released during close stellar flybys, their ejection velocity and likely composition. We study the dependence of the effect of parabolic flybys on the mass ratio between the perturber and parent star, the periastron distance, inclination, and angle of periastron. Whenever ISOs are produced, they leave their parent system typically with velocities of 0.5-2 km/s. This ejection velocity is distinctly different to that of ISOs produced by planet scattering (4-8 km/s) and those shed during the stellar post-main-sequence phase 0.1-0.2 km/s). Using the typical disc truncation radius in various cluster environments, we find that clusters like the Orion nebula cluster are likely to produce the equivalent of 0.85 Earth-masses of ISOs per star. In contrast, clusters like NGC 3603 could produce up to 50 Earth-masses of ISOs per star. Our solar system probably produced the equivalent of 2-3 Earth masses of ISOs, which left our solar system at a mean ejection velocity of 0.7 km/s. Most ISOs produced by flybys should be comet-like, similar to Borisov. ISOs originating from compact long-lived clusters would often show a deficiency in CO. As soon as a statistically significant sample of ISOs is discovered, the combined information of their observed velocities and composition might help in constraining the dominant production process (abridged).

In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass approximately 6.5 x 10^9 M_solar. The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87's spectrum. We can exclude that the simultaneous gamma-ray emission is produced via inverse Compton emission in the same region producing the EHT mm-band emission, and further conclude that the gamma-rays can only be produced in the inner jets (inward of HST-1) if there are strongly particle-dominated regions. Direct synchrotron emission from accelerated protons and secondaries cannot yet be excluded.

Recent water line observations toward several low-mass protostars suggest low water gas fractional abundances in the inner warm envelopes. Water destruction by X-rays has been proposed to influence the water abundances in these regions, but the detailed chemistry, including the nature of alternative oxygen carriers, is not yet understood. In this study, we aim to understand the impact of X-rays on the composition of low-mass protostellar envelopes, focusing specifically on water and related oxygen bearing species. We compute the chemical composition of two low-mass protostellar envelopes using a 1D gas-grain chemical reaction network, under various X-ray field strengths. According to our calculations, outside the water snowline, the water gas abundance increases with $L_{\mathrm{X}}$. Inside the water snowline, water maintains a high abundance of $\sim 10^{-4}$ for small $L_{\mathrm{X}}$, with water and CO being the dominant oxygen carriers. For large $L_{\mathrm{X}}$, the water gas abundances significantly decrease just inside the water snowline (down to $\sim10^{-8}-10^{-7}$) and in the innermost regions ($\sim10^{-6}$). For these cases, the O$_{2}$ and O gas abundances reach $\sim 10^{-4}$ within the water snowline, and they become the dominant oxygen carriers. The HCO$^{+}$ and CH$_{3}$OH abundances, which have been used as tracers of the water snowline, significantly increase/decrease within the water snowline, respectively, as the X-ray fluxes become larger. The abundances of some other dominant molecules, such as CO$_{2}$, OH, CH$_{4}$, HCN, and NH$_{3}$, are also affected by strong X-ray fields, especially within their own snowlines. These X-ray effects are larger in lower density envelope models. Future observations of water and related molecules (using e.g., ALMA and ngVLA) will access the regions around protostars where such X-ray induced chemistry is effective.

We present spectroscopy of individual stars in 26 Magellanic Cloud (MC) star clusters with the aim of estimating dynamical masses and $V$-band mass-to-light ($M/L_V$) ratios over a wide range in age and metallicity. We obtained 3137 high-resolution stellar spectra with M2FS on the \textit{Magellan}/Clay Telescope. Combined with 239 published spectroscopic results of comparable quality, we produced a final sample of 2787 stars with good quality spectra for kinematic analysis in the target clusters. Line-of-sight velocities measured from these spectra and stellar positions within each cluster were used in a customized expectation-maximization (EM) technique to estimate cluster membership probabilities. Using appropriate cluster structural parameters and corresponding single-mass dynamical models, this technique ultimately provides self-consistent total mass and $M/L_V$ estimates for each cluster. Mean metallicities for the clusters were also obtained and tied to a scale based on calcium IR triplet metallicites. We present trends of the cluster $M/L_V$ values with cluster age, mass and metallicity, and find that our results run about 40 per cent on average lower than the predictions of a set of simple stellar population (SSP) models. Modified SSP models that account for internal and external dynamical effects greatly improve agreement with our results, as can models that adopt a strongly bottom-light IMF. To the extent that dynamical evolution must occur, a modified IMF is not required to match data and models. In contrast, a bottom-heavy IMF is ruled out for our cluster sample as this would lead to higher predicted $M/L_V$ values, significantly increasing the discrepancy with our observations.

Galaxy clusters are a promising probe of late-time structure growth, but constraints on cosmology from cluster abundances are currently limited by systematics in their inferred masses. One unmitigated systematic effect in weak-lensing mass inference is ignoring the presence of baryons and treating the entire cluster as a dark matter halo. In this work we present a new flexible model for cluster densities that captures both the baryonic and dark matter profiles, a new general technique for calculating the lensing signal of an arbitrary density profile, and a methodology for stacking those lensing signal to appropriately model stacked weak-lensing measurements of galaxy cluster catalogues. We test this model on 1400 simulated clusters. Similarly to previous studies, we find that a dark matter-only model overestimates the average mass by $7.5\%$, but including our baryonic term reduces that to $0.7\%$. Additionally, to mitigate the computational complexity of our model, we construct an emulator (surrogate model) which accurately interpolates our model for parameter inference, while being much faster to use than the raw model. We also provide an open-source software framework for our model and emulator, called maszcal, which will serve as a platform for continued efforts to improve these mass-calibration techniques. In this work, we detail our model, the construction of the emulator, and the tests which we used to validate that our model does mitigate bias. Lastly, we describe tests of the emulator's accuracy

We describe the X-ray pulse profile models we use, and how we use them, to analyze Neutron Star Interior Composition Explorer (NICER) observations of rotation-powered millisecond pulsars to obtain information about the mass-radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we analyzed in Miller et al. (2019) and Riley et al. (2019) and describe a cross-verification of computations of the pulse profiles of a star with R/M 3, in case stars this compact need to be considered in future analyses. We also present our early cross-verification efforts of the parameter estimation procedures used by Miller et al. (2019) and Riley et al. (2019) by analyzing two distinct synthetic data sets. Both codes yielded credible regions in the mass-radius plane that are statistically consistent with one another and both gave posterior distributions for model parameter values consistent with the values that were used to generate the data. We also summarize the additional tests of the parameter estimation procedure of Miller et al. (2019) that used synthetic pulse profiles and the NICER pulse profile of PSR J0030+0451. We then illustrate how the precision of mass and radius estimates depends on the pulsar's spin rate and the size of its hot spot by analyzing four different synthetic pulse profiles. Finally, we assess possible sources of systematic error in these estimates made using this technique, some of which may warrant further investigation.

Using the first long-term photometry from the Transiting Exoplanet Survey Satellite, we have carried out a detailed time-resolved timing analysis of an intermediate polar TX Col. The power spectra of almost 52 days continuous time-series data reveal the orbital period of $5.691 \pm 0.006$ hr, spin period of $1909.5 \pm 0.2$ s, and beat period of $2105.76 \pm 0.25$ s, which is consistent with the earlier results. We have also found the presence of quasi-periodic oscillations (QPOs) for a few days with a period of 5850-5950 s, which appears to be due to the beating of the Keplerian period with the spin period of the white dwarf. The continuous data allowed us to look thoroughly at the day-wise evolution of the system's accretion geometry. We report here that the TX Col changes its accretion mechanism even on a time-scale of one day, confirming its variable disk-overflow accretion nature. For the majority of the time, it was found to be disk-overflow system with stream-fed dominance, however, pure disk-fed and pure stream-fed accretions cannot be ruled out during the observations.

Of profound astrobiological interest is that not only does Enceladus have a water ocean, but it also appears to be salty, important for its likely habitability. Here, we investigate how salinity affects ocean dynamics and equilibrium ice shell geometry and use knowledge of ice shell geometry and tidal heating rates to help constrain ocean salinity. We show that the vertical overturning circulation of the ocean, driven from above by melting and freezing and the temperature dependence of the freezing point of water on pressure, has opposing signs at very low and very high salinities. In both cases, heat and freshwater converges toward the equator, where the ice is thick, acting to homogenize thickness variations. In order to maintain observed ice thickness variations, ocean heat convergence should not overwhelm heat loss rates through the equatorial ice sheet. This can only happen when the ocean's salinity has intermediate values, order $20$~psu. In this case polar-sinking driven by meridional temperature variations is largely canceled by equatorial-sinking circulation driven by salinity variations and a consistent ocean circulation, ice shell geometry and tidal heating rate can be achieved.

A typical galaxy survey geometry results in galaxy pairs of different separation and angle to the line-of-sight having different distributions in redshift and consequently a different effective redshift. However, clustering measurements are analysed assuming that the clustering is representative of that at a single effective redshift. We investigate the impact of variations in the galaxy-pair effective redshift on the large-scale clustering measured in galaxy surveys. We find that galaxy surveys spanning a large redshift range have different effective redshifts as a function of both pair separation and angle. Furthermore, when considering tracers whose clustering amplitude evolves strongly with redshift, this combination can result in an additional scale-dependent clustering anisotropy. We demonstrate the size of this effect on the eBOSS DR16 Quasar sample and show that, while the impact on monopole is negligible, neglecting this effect can result in a large-scale tilt of $\sim 4\%$ and $\sim40\%$ in quadrupole and hexadecapole, respectively. We discuss strategies to mitigate this effect when making measurements.

Single gap glass Resistive Plate Chambers with areas of 17x17 cm2 and 30x30 cm2 are constructed in the Particle Detector Laboratory in the Sahand University of Technology. Simulation of the gas flow through the chambers is performed using commercial ANSYS-Fluent package. We have shown that flow rate has a linear relation with the pressure of the gas inside the chamber. We have also investigated the dependence of the gas pressure on the length of the hose connected to the outlet of the chamber and shown that it varies linearly. Simulation results are compared with experimentally measured values.

We consider a specific dark energy model, which only includes the Lagrangian up to the cubic order in terms of the vector field self-interactions in the generalized Proca theory. We examine the cosmological parameters in the model by using the data sets of CMB and CMB+HST, respectively. In particular, the Hubble constant is found to be $H_0=71.80^{+1.07}_{-0.72}$ ($72.48^{+0.72}_{-0.60}$) $\rm kms^{-1}Mpc^{-1}$ at $68\%$~C.L. with CMB (CMB+HST), which would alleviate the Hubble constant tension. We also obtain that the reduced $\chi^2$ values in our model are close to unity when fitting with CMB and CMB+HST, illustrating that our model is a good candidate to describe the cosmological evolutions of the universe.

In this paper, we propose a generalized natural inflation (GNI) model to study axion-like particle (ALP) inflation and dark matter (DM). GNI contains two additional parameters $(n_1, n_2)$ in comparison with the natural inflation, that make GNI more general. The $n_1$ build the connection between GNI and other ALP inflation model, $n_2$ controls the inflaton mass. After considering the cosmic microwave background and other cosmological observation limits, the model can realize small-field inflation with a wide mass range, and the ALP inflaton considering here can serve as the DM candidate for certain parameter spaces.

Decays of mesons produced in cosmic ray induced air showers in Earth's atmosphere can lead to a flux of light exotic particles which can be detected in underground experiments. We evaluate the energy spectra of the light neutral mesons $\pi^0$, $\eta$, $\rho^0$, $\omega$, $\phi$ and $J/\psi$ produced in interactions of cosmic ray protons and helium nuclei with air using QCD inspired event generators. Summing up the mesons produced in the individual hadronic interactions of air showers, we obtain the resulting fluxes of undecayed mesons. As an application, we re-consider the case of millicharged particles created in the electromagnetic decay channels of neutral mesons.

We derive two new forms of the K\'arm\'an-Howarth-Monin equation for decaying compressible Hall magnetohydrodynamic (MHD) turbulence. We test them on results of a weakly-compressible, two-dimensional, moderate-Reynolds-number Hall MHD simulation and compare them with an isotropic spectral transfer (ST) equation. The KHM and ST equations are automatically satisfied during the whole simulation owing to the periodic boundary conditions and have complementary cumulative behavior. They are used here to analyze the onset of turbulence and its properties when it is fully developed. These approaches give equivalent results characterizing: the decay of the kinetic + magnetic energy at large scales, the MHD and Hall cross-scale energy transfer/cascade, the pressure dilatation, and the dissipation. The Hall cascade appears when the MHD one brings the energy close to the ion inertial range and is related to the formation of reconnecting current sheets. At later times, the pressure-dilation energy-exchange rate oscillates around zero with no net effect on the cross-scale energy transfer when averaged over a period of its oscillations. A reduced one-dimensional analysis suggests that all three methods may be useful to estimate the energy cascade rate from in situ observations.

For the first time the scientific community in Latin America working at the forefront of research in high energy, cosmology and astroparticle physics (HECAP) have come together to discuss and provide scientific input towards the development of a regional strategy. The present document, the Latin American HECAP Physics Briefing Book, is the result of this ambitious bottom-up effort. This report contains the work performed by the Preparatory Group to synthesize the main contributions and discussions for each of the topical working groups. This briefing book discusses the relevant emerging projects developing in the region and considers potentially impactful future initiatives and participation of the Latin American HECAP community in international flagship projects to provide the essential input for the creation of a long-term HECAP strategy in the region.

Recent astronomical data have provided the primordial deuterium abundance with percent precision. As a result, Big Bang nucleosynthesis may provide a constraint on the universal baryon to photon ratio that is as precise as, but independent from, analyses of the cosmic microwave background. However, such a constraint requires that the nuclear reaction rates governing the production and destruction of primordial deuterium are sufficiently well known. Here, a new measurement of the $^2$H($p,\gamma$)$^3$He cross section is reported. This nuclear reaction dominates the error on the predicted Big Bang deuterium abundance. A proton beam of 400-1650keV beam energy was incident on solid titanium deuteride targets, and the emitted $\gamma$-rays were detected in two high-purity germanium detectors at angles of 55$^\circ$ and 90$^\circ$, respectively. The deuterium content of the targets has been obtained in situ by the $^2$H($^3$He,$p$)$^4$He reaction and offline using the Elastic Recoil Detection method. The astrophysical S-factor has been determined at center of mass energies between 265 and 1094 keV, addressing the uppermost part of the relevant energy range for Big Bang nucleosynthesis and complementary to ongoing work at lower energies. The new data support a higher S-factor at Big Bang temperatures than previously assumed, reducing the predicted deuterium abundance.

In the absence of an initial seed, the Biermann battery term of a non-ideal induction equation acts as a source that generates weak magnetic fields. These fields are then amplified via a dynamo mechanism. The Kelvin-Helmholtz instability is a fluid phenomenon that takes place in many astrophysical scenarios and can trigger the action of the Biermann battery and dynamo processes. We aim to investigate the effect that the ionisation degree of the plasma and the interaction between the charged and neutral species have on the generation and amplification of magnetic fields during the different stages of the instability. We use the two-fluid model implemented in the numerical code Mancha-2F. We perform 2D simulations starting from a configuration with no initial magnetic field and which is unstable due to a velocity shear. We vary the ionisation degree of the plasma and we analyse the role that the different collisional terms included in the equations of the model play on the evolution of the instability and the generation of magnetic field. We find that when no collisional coupling is considered between the two fluids, the effect of the Biermann battery mechanism does not depend on the ionisation degree. However, when elastic collisions are taken into account, the generation of magnetic field is increased as the ionisation degree is reduced. This behaviour is slightly enhanced if the process of charge-exchange is also considered. We also find a dependence on the total density of the plasma related to the dependence on the coupling degree between the two fluids. As the total density is increased, the results from the two-fluid model converge to the predictions of single-fluid models.

We present fully general relativistic simulations of the quasi-circular inspiral and merger of charged, non-spinning, binary black holes with charge-to-mass ratio $\lambda \le 0.3$. We discuss the key features that enabled long term and stable evolutions of these binaries. We also present a formalism for computing the angular momentum carried away by electromagnetic waves, and the electromagnetic contribution to black-hole horizon properties. We implement our formalism and present the results for the first time in numerical-relativity simulations. In addition, we compare our full non-linear solutions with existing approximate models for the inspiral and ringdown phases. We show that Newtonian models based on the quadrupole approximation have errors of 20 % - 100 % in key gauge-invariant quantities. On the other hand, for the systems considered, we find that estimates of the remnant black hole spin based on the motion of test particles in Kerr-Newman spacetimes agree with our non-linear calculations to within a few percent. Finally, we discuss the prospects for detecting black hole charge by future gravitational-wave detectors using either the inspiral-merger-ringdown signal or the ringdown signal alone.

Earth's modern climate is characterized by wet, rainy deep tropics, however paleoclimate and planetary science have revealed a wide range of hydrological cycle regimes connected to different external parameters. Here we investigate how surface wetness affects the tropical hydrological cycle. When surface wetness is decreased in an Earth-like general circulation model, the tropics remain wet but transition from a rainy to rain-free regime. The rain-free regime occurs when surface precipitation is suppressed as negative evaporation (surface condensation) balances moisture flux convergence. The regime transition is dominated by near-surface relative humidity changes in contrast to the hypothesis that relative humidity changes are small. We show near-surface relative humidity changes responsible for the regime transition are controlled by re-evaporation of stratiform precipitation near the lifting condensation level. Re-evaporation impacts the near-surface through vertical mixing. Our results reveal a new rain-free tropical hydrological cycle regime that goes beyond the wet/dry paradigm.

We apply the analogy between gravitational fields and optical media in the general relativistic geometric optics framework to describe how light can acquire orbital angular momentum (OAM) when it traverses the gravitational field of a massive rotating compact object and the interplay between OAM and polarization. Kerr spacetimes are known not only to impose a gravitational Faraday rotation on the polarization of a light beam, but also to set a characteristic fingerprint in the orbital angular momentum distribution of the radiation passing nearby a rotating black hole (BH). Kerr spacetime behaves like an inhomogeneous and anisotropic medium, in which light can acquire orbital angular momentum and spin-to-orbital angular momentum conversion can occur, acting as a polarization and phase changing medium for the gravitationally lensed light, as confirmed by the data analysis of M87* black hole.

In this paper, we proposed a new Monte Carlo radiative transport (MCRT) scheme, which is based completely on the Neumann series solution of Fredholm integral equation. This scheme indicates that the essence of MCRT is the calculation of infinite terms of multiple integrals in Neumann solution simultaneously. Under this perspective we redescribed MCRT procedure systematically, in which the main work amounts to choose an associated probability distribution function (PDF) for a set of random variables and the corresponding unbiased estimation functions. We can select a relatively optimal estimation procedure that has a lower variance from an infinite possible choices, such as the term by term estimation. In this scheme, MCRT can be regarded as a pure problem of integral evaluation, rather than as the tracing of random walking photons. Keeping this in mind, one can avert some subtle intuitive mistakes. In addition the $\delta$-functions in these integrals can be eliminated in advance by integrating them out directly. This fact together with the optimal chosen random variables can remarkably improve the Monte Carlo (MC) computational efficiency and accuracy, especially in systems with axial or spherical symmetry. An MCRT code, Lemon (Linear Integral Equations' Monte Carlo Solver Based on the Neumann solution), has been developed completely based on this scheme. Finally, we intend to verify the validation of Lemon, a suite of test problems mainly restricted to flat spacetime have been reproduced and the corresponding results are illustrated in detail.