Search

We present the completed KMOS3D survey, an integral field spectroscopic survey of 739 galaxies at 0.6 < z < 2.7 using the K-band Multi Object Spectrograph (KMOS) at the Very Large Telescope. The KMOS3D survey provides a population-wide census of kinematics, star formation, outflows, and nebular gas conditions both on and off the star-forming galaxy main sequence through the spatially resolved and integrated properties of Hα, [N ii], and [S ii] emission lines. We detect Hα emission for 91% of galaxies on the main sequence of star formation and 79% overall. The depth of the survey has allowed us to detect galaxies with star formation rates below 1 M⊙ yr−1, as well as to resolve 81% of detected galaxies with ≥3 resolution elements along the kinematic major axis. The detection fraction of Hα is a strong function of both color and offset from the main sequence, with the detected and nondetected samples exhibiting different spectral energy distribution shapes. Comparison of Hα and UV+IR star formation rates reveal that dust attenuation corrections may be underestimated by 0.5 dex at the highest masses (). We confirm our first year results of a high rotation-dominated fraction (monotonic velocity gradient and vrot/) of 77% for the full KMOS3D sample. The rotation-dominated fraction is a function of both stellar mass and redshift, with the strongest evolution measured over the redshift range of the survey for galaxies with . With this paper, we include a final data release of all 739 observed objects (http://www.mpe.mpg.de/ir/KMOS3D). [ABSTRACT FROM AUTHOR]

Although the study of X-ray binaries has led to major breakthroughs in high-energy astrophysics, their circumbinary environment at scales of ∼100–10,000 au has not been thoroughly investigated. In this paper, we undertake a novel and exploratory study by employing direct and high-contrast imaging techniques on a sample of X-ray binaries, using adaptive optics and the vortex coronagraph on Keck/NIRC2. High-contrast imaging opens up the possibility to search for exoplanets, brown dwarfs, circumbinary companion stars, and protoplanetary disks in these extreme systems. Here we present the first near-infrared high-contrast images of 13 high-mass X-ray binaries located within ∼2–3 kpc. The key results of this campaign involve the discovery of several candidate circumbinary companions ranging from substellar (brown dwarf) to stellar masses. By conducting an analysis based on Galactic population models, we discriminate sources that are likely background/foreground stars and isolate those that have a high probability (≳60%–99%) of being gravitationally bound to the X-ray binary. This paper seeks to establish a preliminary catalog for future analyses of proper motion and subsequent observations. With our preliminary results, we calculate the first estimate of the companion frequency and the multiplicity frequency for X-ray binaries: ≈0.6 and 1.8 ± 0.9, respectively, considering only the sources that are most likely bound to the X-ray binary. In addition to extending our comprehension of how brown dwarfs and stars can form and survive in such extreme systems, our study opens a new window to our understanding of the formation of X-ray binaries. [ABSTRACT FROM AUTHOR]

In this paper, we revisit metal-enriched rotating pair-instability supernovae (PISNe) models for metallicities consistent with the Small Magellanic Cloud, the Large Magellanic Cloud (LMC), and 0.1 Z ⊙. By calculating multiple models, we intend to clarify mass ranges and the ejected 56Ni masses from the PISNe, and mass-loss histories for progenitors. We find that the choice of the Wolf–Rayet (WR) mass-loss rates is important, and we adopt the recently proposed rate of Sander & Vink, which covers the mass ranges for PISNe progenitors. We show that slow rotation lowers the PISN range because the core mass increases by rotational mixing. On the other hand, when we assume a typical rotation speed for observed OB stars, the mass-loss increase becomes more significant, and the final stellar masses are lower than in nonrotating models. As a result, the typical mass range for bright supernovae (SNe), with a 56Ni mass higher than 10 M ⊙ for these fast-rotating models is more than 400 and 350 M ⊙ for LMC and 0.1 Z ⊙ metallicities, respectively. It is interesting that unlike in previous works, we find oxygen-rich progenitors for most cases. This O-rich progenitor may be consistent with the recently identified PISN candidate SN2018ibb. He-rich progenitors are seen only for relatively dim and metal-poor (Z ≲ 0.1 Z ⊙) PISNe. We also discuss the black hole mass gap for metal-enriched PISNe, and we show that the upper bound for the gap is lower than in the Population III case. [ABSTRACT FROM AUTHOR]

We present JWST and Atacama Large Millimeter/submillimeter Array (ALMA) imaging for the lensing system SPT0418−47, which includes a strongly lensed, dusty, star-forming galaxy at redshift z = 4.225 and an associated multiply imaged companion. The JWST NIRCam and MIRI imaging observations presented in this paper were acquired as part of the Early Release Science program Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star formation (TEMPLATES). This data set provides robust mutiwavelength detections of stellar light in both the main (SPT0418A) and companion (SPT0418B) galaxies, while the ALMA detection of [C ii ] emission confirms that SPT0418B lies at the same redshift as SPT0418A. We infer that the projected physical separation of the two galaxies is 4.42 ± 0.05 kpc. We derive total magnifications of μ = 29 ± 1 and μ = 4.1 ± 0.7 for SPT0418A and SPT0418B, respectively. We use both prospector and cigale to derive stellar masses. We find that SPT0418A has a stellar mass of M * = 3.4 − 0.6 + 1.1 × 10 10 M ⊙ from prospector or M * = 1.5 ± 0.3 × 1010 M ⊙ from cigale. The stellar mass ratio of SPT0418A and SPT0418B is roughly between 4 and 7 ( 4.2 − 1.6 + 1.9 for prospector and 7.5 ± 3.7 for cigale). We see evidence of extended structure associated with SPT0418A that is suggestive of a tidal feature. These features, along with the close projected proximity, imply that the system is interacting. Interestingly, the star formation rates and stellar masses of both galaxies are consistent with the main sequence of star-forming galaxies at this epoch, indicating that this ongoing interaction has not noticeably elevated the star formation levels. [ABSTRACT FROM AUTHOR]

Observed breakBRD ("break bulges in red disks") galaxies are a nearby sample of face-on disk galaxies with particularly centrally concentrated star formation: they have red disks but recent star formation in their centers as measured by the D n 4000 spectral index. In Kopenhafer et al., a comparable population of breakBRD analogs was identified in the TNG simulation, in which the central concentration of star formation was found to reflect a central concentration of dense, star-forming gas caused by a lack of dense gas in the galaxy outskirts. In this paper, we examine the circumgalactic medium of the central breakBRD analogs to determine if the extended halo gas also shows differences from that around comparison galaxies with comparable stellar mass. We examine the circumgalactic medium gas mass, specific angular momentum, and metallicity in these galaxy populations. We find less gas in the circumgalactic medium of breakBRD galaxies, and that the breakBRD circumgalactic medium is slightly more concentrated than that of comparable M * galaxies. In addition, we find that the angular momentum in the circumgalactic medium of breakBRD galaxies tends to be low for their stellar mass, and shows more misalignment to the angular momentum vector of the stellar disk. Finally, we find that the circumgalactic medium metallicity of breakBRD galaxies tends to be high for their stellar mass. Together with their low star formation rate, we argue that these circumgalactic medium properties indicate a small amount of disk feeding concentrated in the central regions and a lack of low-metallicity gas accretion from the intergalactic medium. [ABSTRACT FROM AUTHOR]

This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. In the first paper, we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method, which includes a new method for calculating fully relativistic stellar self-gravity, and provide details about the outcomes of full disruptions, focusing on the stellar mass dependence of the outcomes for a black hole of mass 106M⊙. We consider eight different stellar masses, from 0.15 M⊙ to 10 M⊙. We find that, relative to the traditional order-of-magnitude estimate rt, the physical tidal radius of low-mass stars (M⋆ ≲ 0.7 M⊙) is larger by tens of percent, while for high-mass stars (M⋆ ≳ 1 M⊙) it is smaller by a factor of 2–2.5. The traditional estimate of the range of energies found in the debris is ≈1.4× too large for low-mass stars, but is a factor of ∼2 too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass loss continues for many stellar vibration times because the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface until the star has reached a distance from the black hole ∼O(10)rt. [ABSTRACT FROM AUTHOR]

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses () and six black hole masses (). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of M⋆ for M⋆/M⊙ ≲ 3. Independent of M⋆, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of ∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of M⋆; at MBH/M⊙ ≃ 5 × 106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of 0.35 − 2, depending on M⋆ and MBH, implying a corresponding change (∝(ΔE)−3/2) in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star's angular momentum and the fraction of M⋆ remaining. Within the "empty loss-cone" regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions. [ABSTRACT FROM AUTHOR]

We present Hubble Space Telescope (HST) imaging data and Canada-France-Hawaii Telescope (CFHT) near-infrared ground-based images for the final sample of 56 candidate galaxy-scale lenses uncovered in the CFHT Legacy Survey as part of the Strong Lensing in the Legacy Survey project. The new images are used to perform lens modeling, measure surface photometry, and estimate stellar masses of the deflector early-type galaxies (ETGs). Lens modeling is performed on the HST images (or CFHT when HST is not available) by fitting the spatially extended light distribution of the lensed features assuming a singular isothermal ellipsoid mass profile and by reconstructing the intrinsic source light distribution on a pixelized grid. Based on the analysis of systematic uncertainties and comparison with inference based on different methods, we estimate that our Einstein radii are accurate to ∼3%. HST imaging provides a much higher success rate in confirming gravitational lenses and measuring their Einstein radii than CFHT imaging does. Lens modeling with ground-based images, however, when successful, yields Einstein radius measurements that are competitive with space-based images. Information from the lens models is used together with spectroscopic information from companion Paper IV to classify the systems, resulting in a final sample of 39 confirmed (grade A) lenses and 17 promising candidates (grade B,C). This represents an increase of half an order of magnitude in sample size with respect to the sample of confirmed lenses studied in Papers I and II. The Einstein radii of the confirmed lenses in our sample span the range 5-15 kpc and are typically larger than those of other surveys, probing the mass in regions where the dark matter contribution is more important. Stellar masses are in the range 1011-1012M☼, covering the range of massive ETGs. The redshifts of the main deflector span a range 0.3 ⩽ zd ⩽ 0.8, which nicely complements low-redshift samples like the Sloan Lens ACS survey and thus provides an excellent sample for the study of the cosmic evolution of the mass distribution of ETGs over the second half of the history of the universe. [ABSTRACT FROM AUTHOR]

Recent work has pointed out the potential existence of a tight relation between the cosmological parameter Ωm, at fixed Ωb, and the properties of individual galaxies in state-of-the-art cosmological hydrodynamic simulations. In this paper, we investigate whether such a relation also holds for galaxies from simulations run with a different code that makes use of a distinct subgrid physics: Astrid. We also find that in this case, neural networks are able to infer the value of Ωm with a ∼10% precision from the properties of individual galaxies, while accounting for astrophysics uncertainties, as modeled in Cosmology and Astrophysics with MachinE Learning (CAMELS). This tight relationship is present at all considered redshifts, z ≤ 3, and the stellar mass, the stellar metallicity, and the maximum circular velocity are among the most important galaxy properties behind the relation. In order to use this method with real galaxies, one needs to quantify its robustness: the accuracy of the model when tested on galaxies generated by codes different from the one used for training. We quantify the robustness of the models by testing them on galaxies from four different codes: IllustrisTNG, SIMBA, Astrid, and Magneticum. We show that the models perform well on a large fraction of the galaxies, but fail dramatically on a small fraction of them. Removing these outliers significantly improves the accuracy of the models across simulation codes. [ABSTRACT FROM AUTHOR]

Liu and collaborators recently proposed an elliptical accretion disk model for tidal disruption events (TDEs). They showed that the accretion disks of optical/UV TDEs are large and highly eccentric and suggested that the broad optical emission lines with complex and diverse profiles originate in a cool eccentric accretion disk of random inclination and orientation. In this paper, we calculate the radiation efficiency of the elliptical accretion disk and investigate the implications for observations of TDEs. We compile observational data for the peak bolometric luminosity and total radiation energy after peak brightness of 18 TDE sources and compare these data to the predictions from the elliptical accretion disk model. Our results show that the observations are consistent with the theoretical predictions and that the majority of the orbital energy of the stellar debris is advected into the black hole (BH) without being converted into radiation. Furthermore, we derive the masses of the disrupted stars and the masses of the BHs of the TDEs. The BH masses obtained in this paper are also consistent with those calculated with the MBH–σ* relation. Our results provide an effective method for measuring the masses of BHs in large numbers of TDEs to be discovered in ongoing and next-generation sky surveys, regardless of whether the BHs are located at the centers of galactic nuclei or wander in disks and halos. [ABSTRACT FROM AUTHOR]

Nuclear star clusters (NSCs) are dense stellar clusters that are found at the centers of a majority of galaxies. In this paper, we study the density profiles for 29 galaxies in a volume-limited survey within 10 Mpc to characterize their NSCs. These galaxies span a 3 × 108–8 × 1010M⊙ and a wide range of Hubble types. We use high-resolution Hubble Space Telescope archival data to create luminosity models for the galaxies using Sérsic profiles to parameterize the NSCs. We also provide estimates for photometric masses of NSCs and their host galaxies using color–M/L relationships and examine their correlation. We use the multi-Gaussian expansion to derive the NSC densities and their 3D mass-density profiles. The 3D density profiles characterize the NSC densities on scales as small as ∼1 pc, approaching the likely spheres of influence for BHs in these objects. We find that these densities correlate with galaxy mass, with NSC density profiles becoming both denser and flatter at higher galaxy masses. Most galaxy NSCs are denser than typical globular clusters. We parameterize the 3D NSC density profiles and their scatter and slope as a function of galaxy stellar mass to enable the construction of realistic nuclear mass profiles. Our fitted profiles and the derived relations are useful in predicting the rate of tidal disruption events in galaxies. We will verify the results of this paper in a follow-up paper that presents the dynamical modeling of the same sample of NSCs. [ABSTRACT FROM AUTHOR]

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius () for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little correlation between SFR and at fixed stellar mass. The Kennicutt–Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed), and this is supported by observations of gas in local star-forming galaxies. We investigate the implications by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time, and the net SFR is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a two-parameter family whose properties are set mainly by halo mass and concentration. These determine the radius and gas accretion rate, which in turn predict how wind strength needs to vary with to keep the SFR constant. [ABSTRACT FROM AUTHOR]

There has been a long-standing factor-of-two tension between the observed star formation rate density and the observed stellar mass buildup after z ∼ 2. Recently, we have proposed that sophisticated panchromatic SED models can resolve this tension, as these methods infer systematically higher masses and lower star formation rates than standard approaches. In a series of papers, we now extend this analysis and present a complete, self-consistent census of galaxy formation over 0.2 < z < 3 inferred with the Prospector galaxy SED-fitting code. In this work, Paper I, we present the evolution of the galaxy stellar mass function using new mass measurements of ∼105 galaxies in the 3D-HST and COSMOS-2015 surveys. We employ a new methodology to infer the mass function from the observed stellar masses: instead of fitting independent mass functions in a series of fixed redshift intervals, we construct a continuity model that directly fits for the redshift evolution of the mass function. This approach ensures a smooth picture of galaxy assembly and makes use of the full, non-Gaussian uncertainty contours in our stellar mass inferences. The resulting mass function has higher number densities at a fixed stellar mass than almost any other measurement in the literature, largely owing to the older stellar ages inferred by Prospector. The stellar mass density is ∼50% higher than previous measurements, with the offset peaking at z ∼ 1. The next two papers in this series will present the new measurements of the star-forming main sequence and the cosmic star formation rate density, respectively. [ABSTRACT FROM AUTHOR]

We continue our comprehensive theoretical investigation of the spectral evolution of white dwarfs based on sophisticated simulations of element transport. In this paper, we focus on the transformation of PG 1159 stars into DO/DB white dwarfs due to the gravitational settling of heavy elements and then into DQ white dwarfs through the convective dredge-up of carbon. We study the impact of several physical parameters on the evolution of the surface carbon abundance over a wide range of effective temperatures. In the hot PG 1159 and DO phases, our calculations confirm that the temperature of the PG 1159-to-DO transition depends sensitively on the stellar mass and the wind mass-loss rate. We show that measured carbon abundances of DOZ white dwarfs are mostly accounted for by our models, with the notable exception of the coolest DOZ stars. In the cooler DB and DQ phases, the predicted atmospheric composition is strongly influenced by the stellar mass, the thickness of the envelope, the initial carbon content, the efficiency of convective overshoot, and the presence of residual hydrogen. We demonstrate that, under reasonable assumptions, our simulations reproduce very well the observed carbon abundance pattern of DQ stars, which thus allows us to constrain the extent of the overshoot region in cool helium-rich white dwarfs. We also argue that our calculations naturally explain a number of recent empirical results, such as the relative excess of low-mass DQ stars and the presence of trace hydrogen and/or carbon at the surface of most DC and DZ stars. [ABSTRACT FROM AUTHOR]

We present kinematics of six local extremely metal-poor galaxies (EMPGs) with low metallicities (0.016–0.098 Z ⊙) and low stellar masses (104.7–107.6 M ⊙). Taking deep medium/high-resolution (R ∼ 7500) integral-field spectra with 8.2 m Subaru, we resolve the small inner velocity gradients and dispersions of the EMPGs with H α emission. Carefully masking out substructures originating by inflow and/or outflow, we fit three-dimensional disk models to the observed H α flux, velocity, and velocity dispersion maps. All the EMPGs show rotational velocities (v rot) of 5–23 km s−1 smaller than the velocity dispersions (σ 0) of 17–31 km s−1, indicating dispersion-dominated (v rot/ σ 0 = 0.29–0.80 < 1) systems affected by inflow and/or outflow. Except for two EMPGs with large uncertainties, we find that the EMPGs have very large gas-mass fractions of f gas ≃ 0.9–1.0. Comparing our results with other H α kinematics studies, we find that v rot/ σ 0 decreases and f gas increases with decreasing metallicity, decreasing stellar mass, and increasing specific star formation rate. We also find that simulated high- z (z ∼ 7) forming galaxies have gas fractions and dynamics similar to the observed EMPGs. Our EMPG observations and the simulations suggest that primordial galaxies are gas-rich dispersion-dominated systems, which would be identified by the forthcoming James Webb Space Telescope observations at z ∼ 7. [ABSTRACT FROM AUTHOR]

We present Keck/LRIS follow-up spectroscopy for 13 photometric candidates of extremely metal-poor galaxies (EMPGs) selected by a machine-learning technique applied to the deep (∼26 AB mag) optical and wide-area (∼500 deg2) Subaru imaging data in the EMPRESS survey. Nine out of the 13 candidates are EMPGs with an oxygen abundance (O/H) less than ∼10% solar value (O/H)⊙, and four sources are contaminants of moderately metal-rich galaxies or no emission-line objects. Notably, two out of the nine EMPGs have extremely low stellar masses and oxygen abundances of 5 Ă— 104â€"7 Ă— 105 M ⊙ and 2%â€"3% (O/H)⊙, respectively. With a sample of five EMPGs with (Fe/O) measurements, two (three) of which are taken from this study (the literature), we confirm that two EMPGs with the lowest (O/H) ratios of ∼2% (O/H)⊙ show high (Fe/O) ratios of ∼0.1, close to the solar abundance ratio. Comparing galaxy chemical enrichment models, we find that the two EMPGs cannot be explained by a scenario of metal-poor gas accretion/episodic star formation history due to their low (N/O) ratios. We conclude that the two EMPGs can be reproduced by the inclusion of bright hypernovae and/or hypothetical pair-instability supernovae (SNe) preferentially produced in a metal-poor environment. This conclusion implies that primordial galaxies at z ∼ 10 could have a high abundance of Fe that did not originate from Type Ia SNe with delays and that Fe may not serve as a cosmic clock for primordial galaxies. [ABSTRACT FROM AUTHOR]

Observations indicate that the star formation rate (SFR) of nuclear rings varies considerably with time and is sometimes asymmetric rather than being uniform across a ring. To understand what controls temporal and spatial distributions of ring star formation, we run semiglobal, hydrodynamic simulations of nuclear rings subject to time-varying and/or asymmetric mass inflow rates. These controlled variations in the inflow lead to variations in the star formation, while the ring orbital period (18 Myr) and radius (600 pc) remain approximately constant. We find that both the mass inflow rate and supernova feedback affect the ring SFR. An oscillating inflow rate with period Î" Ď„ in and amplitude 20 causes large-amplitude (a factor of ≳5), quasi-periodic variations of the SFR when Î" Ď„ in ≳ 50 Myr. We find that the time-varying interstellar medium (ISM) weight and midplane pressure track each other closely, establishing an instantaneous vertical equilibrium. The measured time-varying depletion time is consistent with the prediction from self-regulation theory provided the time delay between star formation and supernova feedback is taken into account. The supernova feedback is responsible only for small-amplitude (a factor of ∼2) fluctuations of the SFR with a timescale ≲40 Myr. Asymmetry in the inflow rate does not necessarily lead to asymmetric star formation in nuclear rings. Only when the inflow rate from one dust lane is suddenly increased by a large factor do the rings undergo a transient period of lopsided star formation. [ABSTRACT FROM AUTHOR]

Fuzzy dark matter (FDM), a scalar particle coupled to the gravitational field without self-interaction, whose mass range is m ∼ 10âˆ'24â€"10âˆ'20 eV, is one of the promising alternative dark matter candidates to cold dark matter. The quantum interference pattern, which is a unique structure of FDM, can be seen in halos in cosmological FDM simulations. In this paper, we first provide an analytic model of the subgalactic matter power spectrum originating from quantum clumps in FDM halos, in which the density distribution of the FDM is expressed by a superposition of quantum clumps whose size corresponds to the de Broglie wavelength of the FDM. These clumps are assumed to be distributed randomly, such that the ensemble average density follows a halo profile such as the Navarroâ€"Frenkâ€"White profile. We then compare the convergence power spectrum projected along the line of sight around the Einstein radius, which is converted from the subgalactic matter power spectrum, to that measured in the strong lens system SDSS J0252 + 0039. While we find that the current observation provides no useful constraint on the FDM mass, we show that future deep, high spatial resolution observations of strong lens systems can tightly constrain FDM with a mass around 10âˆ'22 eV. [ABSTRACT FROM AUTHOR]

This is the first paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the intermediate to the (super)massive regime. In the present work, we focus on stellar BHs and provide an ab initio computation of their mass function across cosmic times; we mainly consider the standard, and likely dominant, production channel of stellar-mass BHs constituted by isolated single/binary star evolution. Specifically, we exploit the state-of-the-art stellar and binary evolutionary code SEVN , and couple its outputs with redshift-dependent galaxy statistics and empirical scaling relations involving galaxy metallicity, star formation rate and stellar mass. The resulting relic mass function as a function of the BH mass m• features a rather flat shape up to m• ≈ 50 M⊙ and then a log-normal decline for larger masses, while its overall normalization at a given mass increases with decreasing redshift. We highlight the contribution to the local mass function from isolated stars evolving into BHs and from binary stellar systems ending up in single or binary BHs. We also include the distortion on the mass function induced by binary BH mergers, finding that it has a minor effect at the high-mass end. We estimate a local stellar BH relic mass density of ρ• ≈ 5 × 107M⊙ Mpc−3, which exceeds by more than two orders of magnitude that in supermassive BHs; this translates into an energy density parameter Ω• ≈ 4 × 10−4, implying that the total mass in stellar BHs amounts to ≲1% of the local baryonic matter. We show how our mass function for merging BH binaries compares with the recent estimates from gravitational wave observations by LIGO/Virgo, and discuss the possible implications for dynamical formation of BH binaries in dense environments like star clusters. We address the impact of adopting different binary stellar evolution codes (SEVN and COSMIC) on the mass function, and find the main differences to occur at the high-mass end, in connection with the numerical treatment of stellar binary evolution effects. We highlight that our results can provide a firm theoretical basis for a physically motivated light seed distribution at high redshift, to be implemented in semi-analytic and numerical models of BH formation and evolution. Finally, we stress that the present work can constitute a starting point to investigate the origin of heavy seeds and the growth of (super)massive BHs in high-redshift star-forming galaxies, that we will pursue in forthcoming papers. [ABSTRACT FROM AUTHOR]

AllS single stars that are born with masses up to 8.5–10 M ⊙ will end their lives as white dwarf (WD) stars. In this evolutionary stage, WDs enter the cooling sequence, where the stars radiate away their thermal energy and are basically cooling. As these stars cool, they reach temperatures and conditions that cause the stars to pulsate. Using differential photometry to produce light curves, we can determine the observed periods of pulsation from the WD. We used the White Dwarf Evolution Code (WDEC) to calculate a grid of over one million models with various temperature, stellar mass, and mass of helium and hydrogen layers and calculated their theoretical pulsation periods. In this paper, we describe our approach to WD asteroseismology using WDEC models, and we present seismological studies for 29 observed DAVs in the Kepler and K2 data sets, 25 of which have never been analyzed using these observations and 19 of which have never been seismically analyzed in any capacity before. Learning about the internal structure of WDs places important constraints on the WD cooling sequence and our overall understanding of stellar evolution for low-mass stars. [ABSTRACT FROM AUTHOR]

We present a suite of high-resolution cosmological zoom-in simulations to z = 4 of a 1012M⊙ halo at z = 0, obtained using seven contemporary astrophysical simulation codes (A rt-I , E nzo , R amses , C hanga , G adget-3 , G ear , and G izmo) widely used in the numerical galaxy formation community. The physics prescriptions for gas cooling and heating and star formation are the same as the ones used in our previous Assembling Galaxies of Resolved Anatomy (AGORA) disk comparison but now account for the effects of cosmological processes such as the expansion of the universe, intergalactic gas inflow, and the cosmic ultraviolet background radiation emitted by massive stars and quasars. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. In the first two steps, we methodically calibrate the gas physics, such as cooling and heating, in simulations without star formation. In the third step, we seek agreement on the total stellar mass produced with the common star formation prescription used in the AGORA disk comparison, in stellar-feedback-free simulations. In the last calibration step, we activate stellar feedback, where each code group is asked to set the feedback prescription to as close to the most widely used one in its code community as possible, while aiming for convergence in the stellar mass at z = 4 to the values predicted by semiempirical models. After all the participating code groups successfully complete the calibration steps, we achieve a suite of cosmological simulations with similar mass assembly histories down to z = 4. With numerical accuracy that resolves the internal structure of a target halo (≲100 physical pc at z = 4), we find that the codes overall agree well with one another, e.g., in gas and stellar properties, but also show differences, e.g., in circumgalactic medium (CGM) properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, high-resolution cosmological zoom-in simulations can have robust and reproducible results. New code groups are invited to join and enrich this comparison by generating equivalent models or to test the code's compatibility on their own, by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the zoom-in simulations presented here will be presented in forthcoming reports from the AGORA Collaboration, including studies of the CGM, simulations by additional codes, and results at lower redshift. [ABSTRACT FROM AUTHOR]

We develop a physical framework for interpreting complex circumstellar patterns whorled around asymptotic giant branch (AGB) stars by investigating stable, coplanar triple systems using hydrodynamic and particle simulations. The introduction of a close tertiary body causes an additional periodic variation in the orbital velocity and trajectory of the AGB star. As a result, the circumstellar outflow builds a fine non-Archimedean spiral pattern superimposed upon the Archimedean spiral produced by the outer binary alone. This fine spiral can be approximated by off-centered circular rings that become tangent to each other at the location of the Archimedean spiral. The superimposed fine pattern fades out relatively quickly as a function of distance from the center of the system, in contrast to the dominant Archimedean spiral pattern, which presents a much slower fractional density decrease with radius. The different rates of radial decrease of the density contrast in the two superimposed patterns, coupled with their different time and spatial scales, lead to an apparent, but illusory radial change in the observed pattern interval, as has been reported, for example, in CW Leo. The function describing the detailed radial dependence of the expansion velocity is different in the two patterns, which may be used to distinguish them. The shape of the circumstellar whorled pattern is further explored as a function of the orbital eccentricity and the inner companion's mass. Although this study is confined to stable, coplanar triple systems, the results are likely applicable to moderately noncoplanar systems and open interesting avenues for studying noncoplanar systems. [ABSTRACT FROM AUTHOR]

We measure chemical abundances for over 20 elements of 15 N-rich field stars with high-resolution (R ∼ 30,000) optical spectra. We find that Na, Mg, Al, Si, and Ca abundances of our N-rich field stars are mostly consistent with those of stars from globular clusters (GCs). Seven stars are estimated to have [Al/Fe ] > 0.5, which is not found in most GC "first generation" stars. On the other hand, α element abundances (especially Ti) could show distinguishable differences between in situ stars and accreted stars. We discover that one interesting star, with consistently low [Mg/Fe], [Si/Fe], [Ca/Fe], [Ti/Fe], [Sc/Fe], [V/Fe], and [Co/Fe], show similar kinematics and [Ba/Eu] as other stars from the dissolved dwarf galaxy "Gaia–Sausage–Enceladus." The α-element abundances and the iron-peak element abundances of the N-rich field stars with metallicities − 1.25 ≤ [Fe/H] ≤ − 0.95 show consistent values with Milky Way field stars (we refer to Milky Way field stars as Milky Way halo field stars unless otherwise specified in this paper) rather than stars from dwarf galaxies, indicating that they were formed in situ. In addition, the neutron-capture elements of N-rich field stars show that most of them could be enriched by asymptotic giant branch stars with masses around 3–5 M⊙. [ABSTRACT FROM AUTHOR]

This paper constructs a theoretical framework for calculating the distribution of masses for gas giant planets forming via the core accretion paradigm. Starting with known properties of circumstellar disks, we present models for the planetary mass distribution over the range 0.1MJ ≤ Mp < 10MJ. If the circumstellar disk lifetime is solely responsible for the end of planetary mass accretion, the observed (nearly) exponential distribution of disk lifetime would imprint an exponential falloff in the planetary mass function. This result is in apparent conflict with observations, which suggest that the mass distribution has a (nearly) power-law form of , with an index of p ≈ 1.3, over the relevant planetary mass range (and for stellar masses ∼0.5–2M⊙). The mass accretion rate onto the planet depends on the fraction of the (circumstellar) disk accretion flow that enters the Hill sphere, and on the efficiency with which the planet captures the incoming material. Models for the planetary mass function that include distributions for these efficiencies, with uninformed priors, can produce nearly power-law behavior, consistent with current observations. The disk lifetimes, accretion rates, and other input parameters depend on the mass of the host star. We show how these variations lead to different forms for the planetary mass function for different stellar masses. Compared to stars with masses M* = 0.5–2M⊙, stars with smaller masses are predicted to have a steeper planetary mass function (fewer large planets). [ABSTRACT FROM AUTHOR]

In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a nonspinning 106M⊙ supermassive black hole. Here we focus on the stellar mass dependence of the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a t−5/3 power law, becoming more like a single peak with a tail declining as t−p with p ≃ 2–5. Surviving remnants are spun-up in the prograde direction and are hotter than main-sequence stars of the same mass. Their specific orbital energy is ≃10−3× that of the debris, but of either sign with respect to the black hole potential, while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential ≲300 km s−1, so they are unable to travel far out into the galactic bulge. The remnants most deeply bound to the black hole go through a second tidal disruption event upon their first return to pericenter; if they have not thermally relaxed, they will be completely disrupted. [ABSTRACT FROM AUTHOR]

A recent comparison by Merritt of simulated and observed Milky Wayâ€"mass galaxies has identified a significant tension between the outskirts (r > 20 kpc) of the stellar halos in simulated and observed galaxies. Using observations from the Dragonfly telescope and simulated galaxies from the Illustris-TNG100 project, Merritt found that the outskirts of stellar halos in simulated galaxies have surface densities 1â€"2 dex higher than those of observed galaxies. In this paper, we compare two suites of 15 simulated Milky Wayâ€"like galaxies, each drawn from the same initial conditions, simulated with the same hydrodynamical code, but with two different models for feedback from supernovae. We find that the McMaster Unbiased Galaxy Simulations (MUGS), which use an older “delayed-cooling” model for feedback, also produce too much stellar mass in the outskirts of the halo, with median surface densities well above observational constraints. The MUGS2 simulations, which instead use a new, physically motivated “superbubble” model for stellar feedback, have 1â€"2 dex lower outer stellar halo masses and surface densities. The MUGS2 simulations generally match both the median surface density profile as well as the scatter in stellar halo surface density profiles seen in observed stellar halos. We conclude that there is no “missing outskirts” problem in cosmological simulations, provided that supernova feedback is modeled in a way that allows it to efficiently regulate star formation in the low-mass progenitor environments of stellar halo outskirts. [ABSTRACT FROM AUTHOR]

X-ray studies of jellyfish galaxies opened a window into the physics of the interplay between the intracluster medium (ICM) and interstellar medium (ISM). In this paper, we present the study of an archival Chandra observation of the GASP jellyfish galaxy JO194. We observe X-ray emission extending from the stellar disk to the unwinding spiral arms with an average temperature of kT = 0.79 ± 0.03 keV. To investigate the origin of the X-ray emission, we compare the observed X-ray luminosities with those expected from the star formation rates (SFRs) obtained from H α emission. We estimate an X-ray luminosity excess of a factor ∼2â€"4 with respect to the SF; therefore, we conclude that SF is not the main event responsible for the extended X-ray emission of JO194. The metallicity in the spiral arms ( Z = 0.24 âˆ' 0.12 + 0.19 Z ⊙ ) is consistent with that of the ICM around JO194 (Z = 0.35 ± 0.07); thus, we suggest that ICM radiative cooling dominates the X-ray emission of the arms. We speculate that the X-ray plasma results from the ISMâ€'ICM interplay, although the nature of this interplay is still mostly unknown. Finally, we observe that the X-ray properties of JO194 are consistent with those of two other GASP galaxies with different stellar mass, phase-space conditions in their hosting clusters, and local ICM conditions. We suggest that the conditions required to induce extended X-ray emission in jellyfish galaxies are established at the beginning of the stripping, and they can persist on long timescales so that galaxies in different clusters and evolutionary stages can present a similar extended X-ray emission. [ABSTRACT FROM AUTHOR]

The semianalytical model a-sloth (Ancient Stars and Local Observables by Tracing Halos) is the first public code that connects the formation of the first stars and galaxies to observables. After several successful projects with this model, we publish the source code (https://gitlab.com/thartwig/asloth) and describe the public version in this paper. The model is based on dark matter merger trees that can either be generated based on Extended Pressâ€"Schechter theory or be imported from dark matter simulations. On top of these merger trees, a-sloth applies analytical recipes for baryonic physics to model the formation of both metal-free and metal-poor stars and the transition between them with unprecedented precision and fidelity. a-sloth samples individual stars and includes radiative, chemical, and mechanical feedback. It is calibrated based on six observables, such as the optical depth to Thomson scattering, the stellar mass of the Milky Way and its satellite galaxies, the number of extremely metal-poor stars, and the cosmic star formation rate density at high redshift. a-sloth has versatile applications with moderate computational requirements. It can be used to constrain the properties of the first stars and high- z galaxies based on local observables, predicts properties of the oldest and most metal-poor stars in the Milky Way, can serve as a subgrid model for larger cosmological simulations, and predicts next-generation observables of the early universe, such as supernova rates or gravitational wave events. [ABSTRACT FROM AUTHOR]

The radial spatial distribution of low-mass satellites around a Milky Way (MW)-like host is an important benchmark for simulations of small-scale structure. The distribution is sensitive to the disruption of subhalos by the central disk and can indicate whether the disruption observed in simulations of MW analogs is artificial (i.e., numerical) or physical in origin. We consider a sample of 12 well-surveyed satellite systems of MW-like hosts in the Local Volume (D < 12 Mpc) that are complete to MV < −9 and within 150 projected kpc. We investigate the radial distribution of satellites and compare with ΛCDM cosmological simulations, including big-box cosmological simulations and high-resolution zoom-in simulations of MW-sized halos. We find that the observed satellites are significantly more centrally concentrated than the simulated systems. Several of the observed hosts, including the MW, are ∼2σ outliers relative to the simulated hosts in being too concentrated, while none of the observed hosts are less centrally concentrated than the simulations. This result is robust to different ways of measuring the radial concentration. We find that this discrepancy is more significant for bright, MV < −12 satellites, suggestive that this is not the result of observational incompleteness. We argue that the discrepancy is possibly due to artificial disruption in the simulations, but, if so, this has important ramifications for what relation between stellar mass and halo mass is allowed in the low-mass regime by the observed abundance of satellites. [ABSTRACT FROM AUTHOR]

We present an analysis of the quenching of star formation in massive galaxies (M * > 109.5 M ⊙) within the first 0.5–3 Gyr of the Universe's history utilizing JWST-CEERS data. We utilize a combination of advanced statistical methods to accurately constrain the intrinsic dependence of quenching in a multidimensional and intercorrelated parameter space. Specifically, we apply random forest classification, area statistics, and a partial correlation analysis to the JWST-CEERS data. First, we identify the key testable predictions from two state-of-the-art cosmological simulations (IllustrisTNG and EAGLE). Both simulations predict that quenching should be regulated by supermassive black hole mass in the early Universe. Furthermore, both simulations identify the stellar potential (ϕ *) as the optimal proxy for black hole mass in photometric data. In photometric observations, where we have no direct constraints on black hole masses, we find that the stellar potential is the most predictive parameter of massive galaxy quenching at all epochs from z = 0–8, exactly as predicted by simulations for this sample. The stellar potential outperforms stellar mass, galaxy size, galaxy density, and Sérsic index as a predictor of quiescence at all epochs probed in JWST-CEERS. Collectively, these results strongly imply a stable quenching mechanism operating throughout cosmic history, which is closely connected to the central gravitational potential in galaxies. This connection is explained in cosmological models via massive black holes forming and growing in deep potential wells, and subsequently quenching galaxies through a mix of ejective and preventative active galactic nucleus feedback. [ABSTRACT FROM AUTHOR]

Recently, a new class of white dwarfs (“slowly cooling WDs”) has been identified in the globular cluster M13. The cooling time of these stars is increased by stable thermonuclear hydrogen burning in their residual envelope. These WDs are thought to be originated by horizontal branch (HB) stars populating the HB blue tail that skipped the asymptotic giant branch phase. To further explore this phenomenon, we took advantage of deep photometric data acquired with the Hubble Space Telescope in the near-ultraviolet and investigate the bright portion of the WD cooling sequence in NGC 6752, another Galactic globular cluster with a metallicity, age, and HB morphology similar to M13. The normalized WD luminosity function derived in NGC 6752 turns out to be impressively similar to that observed in M13, in agreement with the fact that the stellar mass distribution along the HB of these two systems is almost identical. As in the case of M13, the comparison with theoretical predictions is consistent with ∼70% of the investigated WDs evolving at slower rates than standard, purely cooling WDs. Thanks to its relatively short distance from Earth, NGC 6752 photometry reaches a luminosity 1 order of a magnitude fainter than the case of M13, allowing us to sample a regime where the cooling time delay, with respect to standard WD models, reaches ∼300 Myr. The results presented in this paper provide new evidence for the existence of slowly cooling WDs and further support to the scenario proposing a direct causal connection between this phenomenon and the HB morphology of the host stellar cluster. [ABSTRACT FROM AUTHOR]

In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in H i, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2â€"1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA H i survey. The sample spans a stellar mass range of 9 ≤ log M ⋆ / M ⊙ ≤ 11. We study molecular gas radial profiles, isodensity radii, and surface densities as a function of galaxy H i deficiency and morphology. There is a weak correlation between global H i and H2 deficiencies, and resolved properties of molecular gas correlate with H i deficiency: galaxies that have large H i deficiencies have relatively steep and truncated molecular gas radial profiles, which is due to the removal of low-surface-density molecular gas on the outskirts. Therefore, while the environmental mechanisms observed in H i also affect molecular gas reservoirs, there is only a moderate reduction of the total amount of molecular gas. [ABSTRACT FROM AUTHOR]

Accurate stellar ages are essential for our understanding of the star formation history of the Milky Way and Galactic chemical evolution, as well as to constrain exoplanet formation models. Gyrochronology, a relationship between stellar rotation and age, appears to offer a reliable age indicator for main-sequence (MS) stars over the mass range of approximately 0.6â€"1.3 M ⊙. Those stars lose their angular momentum due to magnetic braking and as a result their rotation speeds decrease with age. Although current gyrochronology relations have been fairly well tested for young MS stars with masses greater than 1 M ⊙, primarily in young open clusters, insufficient tests exist for older and lower mass MS stars. Binary stars offer the potential to expand and fill in the range of ages and metallicity over which gyrochronology can be empirically tested. In this paper, we demonstrate a Monte Carlo approach to evaluate gyrochronology models using binary stars. As examples, we used five previously published wide binary pairs. We also demonstrate a Monte Carlo approach to assess the precision and accuracy of ages derived from each gyrochronology model. For the traditional Skumanich models, the age uncertainties are σ age/age = 15%â€"20% for stars with B âˆ' V = 0.65 and σ age/age = 5%â€"10% for stars with B âˆ' V = 1.5 and rotation period P ≤ 20 days. [ABSTRACT FROM AUTHOR]

A rare population of massive disk galaxies have been found to invade the red sequence dominated by early-type galaxies. The formation and origins of these red/quenched massive disk galaxies have recently gained great interest. The quenching mechanisms that are usually proposed, such as bar quenching and environment quenching, do not seem to be suitable for those bulgeless quenched disks in a low-density environment. In this paper, we use the TNG300 simulation to investigate the formation of massive quenched central disk galaxies. It is found that these galaxies contain less gas than their star-forming counterparts and harbor giant super massive black holes (SMBHs; above 108 M ⊙). By tracing their formation history, we found that quenched disk galaxies formed early and preserved disk morphology for cosmological timescales. They have experienced less than one major merger on average, and mini-mergers (mass ratio <1/10) have mainly contributed to the growth of their SMBHs. In the IllustrisTNG simulation, the black hole feedback mode switches from thermal to kinetic feedback when the black hole mass is more massive than ∼108 M ⊙, which is more efficient to eject gas outside of the galaxy and to suppress further cooling of the hot gaseous halo. We conclude that the dominant quenching mechanism in massive red/quenched disk galaxies is kinetic active galactic nuclei feedback. [ABSTRACT FROM AUTHOR]

Physical and chemical properties of the interstellar medium (ISM) at subgalactic (∼kiloparsec) scales play an indispensable role in controlling the ability of gas to form stars. In this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of eight resolved ISM properties in star-forming regions to constrain the areas of this hyperspace where most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M⋆ = 107–1011M⊙) and redshift (0 ≤ z ≤ 3) into kiloparsec-sized regions and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency, and dark-matter volumetric density representative of each region in the context of their star formation activity and environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two distinct environmental regimes, which are brought about by an underlying bicomponent radial SFR profile in galaxies. We examine how the relative prominence of these regimes depends on galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and find similarities as well as departures. Further, using dimensionality reduction, we characterize the aforementioned hyperspace to reveal a high degree of multicollinearity in relationships among ISM properties that drive the distribution of star formation at kiloparsec scales. Based on this, we show that a reduced 3D representation underpinned by a multivariate radius relationship is sufficient to capture most of the variance in the original 8D space. [ABSTRACT FROM AUTHOR]

Magnetospheric clouds have been proposed as explanations for depth-varying dips in the phased light curves of young, magnetically active stars such as σ Ori E and RIK-210. However, the stellar theory that first predicted magnetospheric clouds also anticipated an associated mass-balancing mechanism known as centrifugal breakout for which there has been limited empirical evidence. In this paper, we present data from the Transiting Exoplanet Survey Satellite, Las Cumbres Observatory, All-Sky Automated Survey for Supernovae, and Veloce on the 45 Myr M3.5 star TIC 234284556 , and propose that it is a candidate for the direct detection of centrifugal breakout. In assessing this hypothesis, we examine the sudden (∼1 day timescale) disappearance of a previously stable (∼1 month timescale) transit-like event. We also interpret the presence of an anomalous brightening event that precedes the disappearance of the signal, analyze rotational amplitudes and optical flaring as a proxy for magnetic activity, and estimate the mass of gas and dust present immediately prior to the potential breakout event. After demonstrating that our spectral and photometric data support a magnetospheric cloud and centrifugal breakout model and disfavor alternate scenarios, we discuss the possibility of a coronal mass ejection or stellar wind origin of the corotating material and we introduce a reionization mechanism as a potential explanation for more gradual variations in eclipse parameters. Finally, after comparing TIC 234284556 with previously identified “flux-dip” stars, we argue that TIC 234284556 may be an archetypal representative of a whole class of young, magnetically active stars. [ABSTRACT FROM AUTHOR]

We provide a method for estimating the projected density distribution n ÂŻ 2 w p (r p) of photometric objects around spectroscopic objects in a spectroscopic survey. This quantity describes the distribution of photometric sources with certain physical properties (e.g., luminosity, mass, and color) around cosmic webs (PAC) traced by the spectroscopic objects. The method can make full use of current and future deep and wide photometric surveys to explore the formation of galaxies up to medium redshift (z s < 2) 33 Throughout the paper, we use z s for redshift, z for the z -band magnitude, z p for the photometric redshift. with the aid of cosmological spectroscopic surveys that sample only a fairly limited species of objects (e.g., emission line galaxies). As an example, we apply the PAC method to the CMASS spectroscopic and HSC-SSP PDR2 photometric samples to explore the distribution of galaxies for a wide range of stellar masses from 109.0 to 1012.0 M ⊙ around massive stars at z s ≈ 0.6. Using the abundance-matching method, we model n ÂŻ 2 w p (r p) in N -body simulation using Markov chain Monte Carlo sampling, and we accurately measure the stellarâ€"halo mass relation and stellar mass function for the whole mass range. We can also measure the conditional stellar mass function of satellites for central galaxies of different mass. The PAC method has many potential applications for studying the evolution of galaxies. [ABSTRACT FROM AUTHOR]

We present the lifetime star formation histories (SFHs) for six ultrafaint dwarf (UFD; M V > − 7.0, 4.9 < log 10 (M * (z = 0) / M ⊙) < 5.5) satellite galaxies of M31 based on deep color–magnitude diagrams constructed from Hubble Space Telescope imaging. These are the first SFHs obtained from the oldest main-sequence turnoff of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50% of their stellar mass by z = 5 (12.6 Gyr ago), similar to known UFDs around the MW, but that 10%–40% of their stellar mass formed at later times. We uncover one remarkable UFD, And xiii, which formed only 10% of its stellar mass by z = 5, and 75% in a rapid burst at z ∼ 2–3, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This "young" UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least-massive MW UFDs (M *(z = 5) ≲ 5 × 104 M ⊙) are likely quenched by reionization, whereas more-massive M31 UFDs (M *(z = 5) ≳ 105 M ⊙) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs. [ABSTRACT FROM AUTHOR]

Cold, low-mass, field brown dwarfs are important for constraining the terminus of the stellar mass function, and also for optimizing atmospheric studies of exoplanets. In 2020 new model grids for such objects were made available: Sonora-Bobcat and ATMO 2020. Also, new candidate cold brown dwarfs were announced, and new spectroscopic observations at λ ≈ 4.8 μm were published. In this paper we present new infrared photometry for some of the coldest brown dwarfs, and put the new data and models together to explore the properties of these objects. We reconfirm the importance of mixing in these atmospheres, which leads to CO and NH3 abundances that differ by orders of magnitude from chemical equilibrium values. We also demonstrate that the new models retain the known factor ≳3 discrepancy with observations at 2 ≲ λμm ≲ 4, for brown dwarfs cooler than 600 K. We show that the entire 1 ≲ λμm ≲ 20 energy distribution of six brown dwarfs with 260 ≤ Teff K ≤ 475 can be well reproduced, for the first time, by model atmospheres which include disequilibrium chemistry as well as a photospheric temperature gradient which deviates from the standard radiative/convective equilibrium value. This change to the pressure–temperature profile is not unexpected for rotating and turbulent atmospheres that are subject to diabatic processes. A limited grid of modified-adiabat model colors is generated, and used to estimate temperatures and metallicities for the currently known Y dwarfs. A compilation of the photometric data used here is given in Appendix C. [ABSTRACT FROM AUTHOR]

We developed a parametric Galactic model toward the Galactic bulge by fitting to spatial distributions of the Gaia DR2 disk velocity, VVV proper motion, BRAVA radial velocity, OGLE-III red clump star count, and OGLE-IV star count and microlens rate, optimized for use in microlensing studies. We include the asymmetric drift of Galactic disk stars and the dependence of velocity dispersion on Galactic location in the kinematic model, which has been ignored in most previous models used for microlensing studies. We show that our model predicts a microlensing parameter distribution significantly different from those typically used in previous studies. We estimate various fundamental model parameters for our Galaxy through our modeling, including the initial mass function (IMF) in the inner Galaxy. Combined constraints from star counts and the microlensing event timescale distribution from the OGLE-IV survey, in addition to a prior on the bulge stellar mass, enable us to successfully measure IMF slopes using a broken power-law form over a broad mass range, for M < 0.08 M⊙, for 0.08 M⊙ ≤ M < Mbr, and for M ≥ Mbr, as well as a break mass at. This is significantly different from the Kroupa IMF for local stars, but similar to the Zoccali IMF measured from a bulge luminosity function. We also estimate the dark matter mass fraction in the bulge region of 28% ± 7% which could be larger than a previous estimate. Because our model is purely parametric, it can be universally applied using the parameters provided in this paper. (A tool for microlensing simulation using our Galactic model has been published, Koshimoto & Ranc 2021, and can be downloaded at https://github.com/nkoshimoto/genulens.) [ABSTRACT FROM AUTHOR]

Solar-type stars exhibit their highest levels of magnetic activity during their early convective pre-main-sequence (PMS) phase of evolution. The most powerful PMS flares, superflares and megaflares, have peak X-ray luminosities of erg s−1 and total energies erg. Among >24,000 X-ray-selected young (t ≲ 5 Myr) members of 40 nearby star-forming regions from our earlier Chandra MYStIX and SFiNCs surveys, we identify and analyze a well-defined sample of 1086 X-ray superflares and megaflares, the largest sample ever studied. Most are considerably more powerful than optical/X-ray superflares detected on main-sequence stars. This study presents energy estimates of these X-ray flares and the properties of their host stars. These events are produced by young stars of all masses over evolutionary stages ranging from protostars to diskless stars, with the occurrence rate positively correlated with stellar mass. Flare properties are indistinguishable for disk-bearing and diskless stars indicating star–disk magnetic fields are not involved. A slope of α ≃ 2 in the flare energy distributions is consistent with that of optical/X-ray flaring from older stars and the Sun. Megaflares (erg) from solar-mass stars have an occurrence rate of flares/star/year and contribute at least 10%–20% to the total PMS X-ray energetics. These explosive events may have important astrophysical effects on protoplanetary disk photoevaporation, ionization of disk gas, production of spallogenic radionuclides in disk solids, and hydrodynamic escape of young planetary atmospheres. Our following paper details plasma and magnetic loop modeling of the >50 brightest X-ray megaflares. [ABSTRACT FROM AUTHOR]

Interstellar thermal pressures can be measured using C i absorption lines that probe the pressure-sensitive populations of the fine-structure levels of its ground state. In a survey of C i absorption toward Galactic hot stars, Jenkins & Tripp found evidence of small amounts (∼0.05%) of gas at high pressures () mixed with a more general presence of lower pressure material exhibiting a log-normal distribution that spanned the range. In this paper, we study Milky Way C i lines in the spectra of extragalactic sources instead of Galactic stars and thus measure the pressures without being influenced by regions where stellar mass loss and H ii region expansions could create localized pressure elevations. We find that the distribution of low pressures in the current sample favors slightly higher pressures than the earlier survey, and the fraction of gaseous material at extremely high pressures is about the same as that found earlier. Thus we conclude that the earlier survey was not appreciably influenced by the stellar environments, and the small amounts of high pressure gas indeed exist within the general interstellar medium. [ABSTRACT FROM AUTHOR]

The structure of protostellar cores can often be approximated by isothermal Bonnor–Ebert spheres (BES), which are stabilized by an external pressure. For the typical pressure of 104kB K cm−3 to 105kB K cm−3 found in molecular clouds, cores with masses below 1.5 M⊙ are stable against gravitational collapse. In this paper, we analyze the efficiency of triggering gravitational collapse with a nearby stellar wind, which represents an interesting scenario for triggering low-mass star formation. We analytically derive a new stability criterion for a BES compressed by a stellar wind, which depends on its initial nondimensional radius. If the stability limit is violated the wind triggers a core collapse. Otherwise, the core is destroyed by the wind. We estimate its validity range to and confirm this in simulations with the SPH-Code GADGET-3. The efficiency of triggering a gravitational collapse strongly decreases for since in this case destruction and acceleration of the whole sphere begin to dominate. We were unable to trigger a collapse for , which leads to the conclusion that a stellar wind can move the smallest unstable stellar mass to 0.5 M⊙ and that destabilizing even smaller cores would require external pressure larger than 105kB K cm−3. For the expected wind strength according to our criterion is small enough that the compression is slower than the sound speed of the BES and sound waves can be triggered. In this case our criterion somewhat underestimates the onset of collapse and detailed numerical analyses are required. [ABSTRACT FROM AUTHOR]

We present emission-line measurements and physical interpretations for a sample of 117 [O iii ] emitting galaxies at z = 5.33–6.93, using the first deep JWST/NIRCam wide-field slitless spectroscopic observations. Our 9.7 hr integration is centered upon the z = 6.3 quasar J0100+2802—the first of six fields targeted by the EIGER survey—and covers λ = 3–4 μ m. We detect 133 [O iii ] doublets, but close pairs motivated by their small scale clustering excess. The galaxies are characterized by a UV luminosity M UV ∼ −19.6 (−17.7 to −22.3), stellar mass ∼108 (106.8−10.1) M ⊙, H β and [O iii ]4960+5008 EWs ≈ 850 Å (up to 3000 Å), young ages, a highly excited interstellar medium, and low dust attenuations. These high EWs are very rare in the local universe, but we show they are ubiquitous at z ∼ 6 based on the measured number densities. The stacked spectrum reveals H γ and [O iii ]4364, which shows that the galaxies are typically dust- and metal-poor (E (B − V) = 0.1, 12 + log (O / H) = 7.4) with a high electron temperature (2 × 104 K) and a production efficiency of ionizing photons (ξ ion = 1025.3 Hz erg−1). We further show the existence of a strong mass–metallicity relation. The properties of the stars and gas in z ∼ 6 galaxies conspire to maximize the [O iii ] output from galaxies, yielding an [O iii ] luminosity density at z ≈ 6 that is significantly higher than that at z ≈ 2. Thus, [O iii ] emission-line surveys with JWST prove a highly efficient method to trace the galaxy density in the Epoch of Reionization. [ABSTRACT FROM AUTHOR]

We present element abundance ratios and ionizing radiation of local young low-mass (∼106 M⊙) extremely metal-poor galaxies (EMPGs) with a 2% solar oxygen abundance (O/H)⊙ and a high specific star formation rate (sSFR ∼ 300 Gyr−1) and other (extremely) metal-poor galaxies, which are compiled from Extremely Metal-Poor Representatives Explored by the Subaru Survey (EMPRESS) and the literature. Weak emission lines such as [Fe iii ] λ4658 and He ii λ4686 are detected in very deep optical spectra of the EMPGs taken with 8 m class telescopes, including Keck and Subaru, enabling us to derive element abundance ratios with photoionization models. We find that neon-to-oxygen and argon-to-oxygen ratios are comparable to those of known local dwarf galaxies and that the nitrogen-to-oxygen abundance ratios (N/O) are lower than 20% (N/O)⊙, consistent with the low oxygen abundance. However, the iron-to-oxygen abundance ratios (Fe/O) of the EMPGs are generally high; the EMPGs with the 2%-solar oxygen abundance show high Fe/O ratios of ∼90%–140% (Fe/O)⊙, which are unlikely to be explained by suggested scenarios of Type Ia supernova iron productions, iron's dust depletion, and metal-poor gas inflow onto previously metal-riched galaxies with solar abundances. Moreover, the EMPGs with the 2%-solar oxygen abundance have very high He ii λ4686/Hβ ratios of ∼1/40, which are not reproduced by existing models of high-mass X-ray binaries with progenitor stellar masses <120 M⊙. Comparing stellar-nucleosynthesis and photoionization models with a comprehensive sample of EMPGs identified by this and previous EMPG studies, we propose that both the high Fe/O ratios and the high He ii λ4686/Hβ ratios are explained by the past existence of supermassive (>300 M⊙) stars, which may evolve into intermediate-mass black holes (≳100 M⊙). [ABSTRACT FROM AUTHOR]

The See Change survey was designed to make z > 1 cosmological measurements by efficiently discovering high-redshift Type Ia supernovae (SNe Ia) and improving cluster mass measurements through weak lensing. This survey observed twelve galaxy clusters with the Hubble Space Telescope (HST) spanning the redshift range z = 1.13–1.75, discovering 57 likely transients and 27 likely SNe Ia at z ∼ 0.8–2.3. As in similar previous surveys, this proved to be a highly efficient use of HST for supernova observations; the See Change survey additionally tested the feasibility of maintaining, or further increasing, the efficiency at yet higher redshifts, where we have less detailed information on the expected cluster masses and star formation rates. We find that the resulting number of SNe Ia per orbit is a factor of ∼8 higher than for a field search, and 45% of our orbits contained an active SN Ia within 22 rest-frame days of peak, with one of the clusters by itself yielding 6 of the SNe Ia. We present the survey design, pipeline, and supernova discoveries. Novel features include fully blinded supernova searches, the first random forest candidate classifier for undersampled IR data (with a 50% detection threshold within 0.05 mag of human searchers), real-time forward-modeling photometry of candidates, and semi-automated photometric classifications and follow-up forecasts. We also describe the spectroscopic follow-up, instrumental in measuring host galaxy redshifts. The cosmology analysis of our sample will be presented in a companion paper. [ABSTRACT FROM AUTHOR]

This paper systematically investigates the comoving megaparsec-scale intergalactic medium (IGM) environment around galaxies traced by the Lyα forest. Using our cosmological hydrodynamic simulations, we investigate the IGM–galaxy connection at z = 2 by two methods: (i) cross-correlation analysis between galaxies and the fluctuation of Lyα forest transmission (δF) and (ii) comparison of the overdensity of neutral hydrogen (H i) and galaxies. Our simulations reproduce observed cross-correlation functions (CCFs) between Lyα forest and Lyman-break galaxies. We further investigate the variation of the CCF using subsamples divided by dark matter halo mass (MDH), galaxy stellar mass (M⋆), and star formation rate (SFR) and find that the CCF signal becomes stronger with increasing MDH, M⋆, and SFR. The CCFs between galaxies and gas density fluctuation are also found to have similar trends. Therefore, the variation of δF–CCF depending on MDH, M⋆, and SFR is due to varying gas densities around galaxies. We find that the correlation between galaxies and the IGM H i distribution strongly depends on MDH as expected from linear theory. Our results support the ΛCDM paradigm, confirming a spatial correlation between galaxies and IGM H i , with more massive galaxies being clustered in higher-density regions. [ABSTRACT FROM AUTHOR]

We present detailed observations of photoionization conditions and galaxy kinematics in 11 z = 1.39–2.59 radio-loud quasar host galaxies. Data were taken with the OSIRIS integral field spectrograph and the adaptive optics system at the W. M. Keck Observatory that targeted nebular emission lines (Hβ, [O iii ], Hα, [N ii ]) redshifted into the near-infrared (1–2.4 μm). We detect extended ionized emission on scales ranging from 1 to 30 kpc photoionized by stars, shocks, and active galactic nuclei (AGN). Spatially resolved emission-line ratios indicate that our systems reside off the star formation and AGN-mixing sequence on the Baldwin, Phillips, & Terlevich diagram at low redshift. The dominant cause of the difference between line ratios of low-redshift galaxies and our sample is due to lower gas-phase metallicities, which are 2–5× less compared to galaxies with AGN in the nearby universe. Using gas velocity dispersion as a proxy to stellar velocity dispersion and dynamical mass measurement through inclined disk modeling, we find that the quasar host galaxies are undermassive relative to their central supermassive black hole mass, with all systems residing off the local scaling (M•–σ, M•–M*) relationship. These quasar host galaxies require substantial growth, up to an order of magnitude in stellar mass, to grow into present-day massive elliptical galaxies. Combining these results with part I of our sample paper, we find evidence for winds capable of causing feedback before the AGN host galaxies land on the local scaling relation between black hole and galaxy stellar mass, and before the enrichment of the interstellar medium to a level observed in local galaxies with AGN. [ABSTRACT FROM AUTHOR]

Elliptical accretion disk models for tidal disruption events (TDEs) have been recently proposed and independently developed by two groups. Although these two models are characterized by a similar geometry, their physical properties differ considerably. In this paper, we further investigate the properties of the elliptical accretion disk of the nearly uniform distribution of eccentricity within the disk plane. Our results show that the elliptical accretion disks have distinctive hydrodynamic structures and spectral energy distributions, associated with TDEs. The soft X-ray photons generated at pericenter and nearby are trapped in the disk and advected around the ellipse because of large electron scattering opacity. They are absorbed and reprocessed into emission lines and low-frequency continuum via recombination and bremsstrahlung emission. Because of the rapid increase of bound–free and free–free opacities with radius, the low-frequency continuum photons become trapped in the disk at large radius and are advected through apocenter and back to the photon-trapping radius. Elliptical accretion disks predict sub-Eddington luminosities and emit mainly at the photon-trapping radius of thousands of Schwarzschild radii with a blackbody spectrum of nearly single temperature of typically about 3 × 104 K. Because of the self-regulation, the photon-trapping radius expands and contracts following the rise and fall of accretion rate. The radiation temperature is nearly independent of BH mass and accretion rate and varies weakly with the stellar mass and the viscosity parameter. Our results are well consistent with the observations of optical/UV TDEs. [ABSTRACT FROM AUTHOR]

The Event Horizon Telescope mass estimate for M87* is consistent with the stellar dynamics mass estimate and inconsistent with the gas-dynamics mass estimates by up to 2σ. We have previously explored a new gas-dynamics model that incorporated sub-Keplerian gas velocities and could, in principle, explain the discrepancy in the stellar and gas-dynamics mass estimate. In this paper, we extend this gas-dynamical model to also include non-trivial disk heights, which may also resolve the mass discrepancy independent of sub-Keplerian velocity components. By combining the existing velocity measurements and the Event Horizon Telescope mass estimate, we place constraints on the gas disk inclination and sub-Keplerian fraction. These constraints require the parsec-scale ionized gas disk to be misaligned with the milliarcsecond radio jet by at least 11°, and more typically 27°. Modifications to the gas-dynamics model either by introducing sub-Keplerian velocities or thick disks produce further misalignment with the radio jet. If the jet is produced in a Blandford–Znajek-type process, the angular momentum of the black hole is decoupled with the angular momentum of the large-scale gas feeding M87*. [ABSTRACT FROM AUTHOR]