SEARCH · Engineering Papers

Results for “Star formation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Fundamental problems in astroparticle physics and cosmology (Final Technical Report)

In the standard cosmological model, structure forms when gas cools and condenses to the center of dark matter halos. The onset of structure is governed by rapid star formation and subsequent evolution of the first massive stars. One of the endpoints of stellar evolution is the formation of a black hole. We showed that all black hole merger events must occur after the universe is about 85 million years old. This prediction is important because it means that if a future experiment makes a detection of a black hole merger event at an earlier time then our understanding of cosmology and the onset of structure formation is not fully developed. In addition, we explored dark matter models that have either self-interactions or experience decays. In the first case, self-interactions may lead to the formation of a thick disk that can be co-planar to the Galactic disc due to correlations in angular momentum between dark matter and baryons that are predicted by tidal torque theory. We showed that for heavy dark matter particles, either the self interaction may not be strong enough to solve the small scale structure motivation or a dark disc cannot be present in the Milky Way. In the second-case, we showed that if dark matter experiences decays into two particles, one of which behaves like radiation, then it is possible to relieve the tension that exists between measurements of the local expansion rate. Finally, we developed methods on how to distinguish gamma-ray excesses from promising dark matter targets as compared to conventional astrophysical sources. This is an important toolset that is applicable to any suspected dark matter source.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Searching for Spectroscopic Signatures of Ongoing Quenching in SDSS Galaxies

In this paper, we estimate the “star formation change parameter,” SFR 79 , which characterizes the current SFR relative to the average during the last 800 Myr, for ~300,000 galaxies selected from the Sloan Digital Sky Survey (SDSS). The goals are to examine, in a much larger and independent sample, the trends previously reported in a sample of star-forming (SF) MaNGA galaxies and also to search for spectroscopic signatures of ongoing quenching in the so-called “Green Valley,” which is generally believed to contain galaxies that are migrating from the SF population to the quenched population of galaxies. Measuring SFR 79 for our large sample of SDSS galaxies, we first confirm the basic results of SF galaxies published by Wang & Lilly. We then discuss in detail the calibration and meaning of SFR 79 for galaxies that are well below the SFMS and establish the expected statistical signature of systematic ongoing quenching from modeling the z ~ 0 quenching rate of the SF population. We conclude that it is not possible at present to establish unambiguous observational evidence for systematic ongoing quenching processes, due to limitations both in the noise of the observational data, in particular in the measurements of H δ absorption, and in the calibration of SFR 79 , as well as biases introduced by the necessity of selecting objects with significant H α emission. We do however see plausible indications of ongoing quenching, which are quantitatively consistent with expectations from “growth+quenching” models of galaxy evolution and a typical e-folding timescale for quenching of ~500 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

SDSS-IV MaNGA: stellar rotational support in disc galaxies versus central surface density and stellar population age

ABSTRACT We investigate how the stellar rotational support changes as a function of spatially resolved stellar population age ($D_n4000$) and relative central stellar surface density ($\Delta \Sigma _1$) for MaNGA isolated/central disc galaxies. We find that the galaxy rotational support indicator $\lambda _{R_\mathrm{e}}$ varies smoothly as a function of $\Delta \Sigma _1$ and $D_n4000$. $D_n4000$ versus $\Delta \Sigma _1$ follows a ‘J-shape’, with $\lambda _{R_\mathrm{e}}$ contributing to the scatters. In this ‘J-shaped’ pattern rotational support increases with central $D_n4000$ when $\Delta \Sigma _1$ is low but decreases with $\Delta \Sigma _1$ when $\Delta \Sigma _1$ is high. Restricting attention to low-$\Delta \Sigma _1$ (i.e. large-radius) galaxies, we suggest that the trend of increasing rotational support with $D_n4000$ for these objects is produced by a mix of two different processes, a primary trend characterized by growth in $\lambda _{R_\mathrm{e}}$ along with mass through gas accretion, on top of which disturbance episodes are overlaid, which reduce rotational support and trigger increased star formation. An additional finding is that star-forming galaxies with low $\Delta \Sigma _1$ have relatively larger radii than galaxies with higher $\Delta \Sigma _1$ at fixed stellar mass. Assuming that these relative radii rankings are preserved while galaxies are star forming then implies clear evolutionary paths in central $D_n4000$ versus $\Delta \Sigma _1$. The paper closes with comments on the implications that these paths have for the evolution of pseudo-bulges versus classical bulges. The utility of using $\rm D_n4000$–$\Delta \Sigma _1$ to study $\lambda _{R_\mathrm{e}}$ reinforces the notion that galaxy kinematics correlate both with structure and with stellar-population state, and indicates the importance of a multidimensional description for understanding bulge and galaxy evolution.

Wang, Xiaohan (ORCID:0009000283121002)↗

DOE OSTI2021

Measuring Young Stars in Space and Time. II. The Pre-main-sequence Stellar Content of N44

The Hubble Space Telescope survey Measuring Young Stars in Space and Time (MYSST) entails some of the deepest photometric observations of extragalactic star formation, capturing even the lowest-mass stars of the active star-forming complex N44 in the Large Magellanic Cloud. We employ the new MYSST stellar catalog to identify and characterize the content of young pre-main-sequence (PMS) stars across N44 and analyze the PMS clustering structure. To distinguish PMS stars from more evolved line of sight contaminants, a non-trivial task due to several effects that alter photometry, we utilize a machine-learning classification approach. This consists of training a support vector machine (SVM) and a random forest (RF) on a carefully selected subset of the MYSST data and categorize all observed stars as PMS or non-PMS. Combining SVM and RF predictions to retrieve the most robust set of PMS sources, we find ∼26,700 candidates with a PMS probability above 95% across N44. Employing a clustering approach based on a nearest neighbor surface density estimate, we identify 16 prominent PMS structures at 1σ significance above the mean density with sub-clusters persisting up to and beyond 3σ significance. The most active star-forming center, located at the western edge of N44's bubble, is a subcluster with an effective radius of ∼5.6 pc entailing more than 1100 PMS candidates. Furthermore, we confirm that almost all identified clusters coincide with known H ii regions and are close to or harbor massive young O stars or YSOs previously discovered by MUSE and Spitzer observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

HD 143811 AB b: A Directly Imaged Planet Orbiting a Spectroscopic Binary in Sco-Cen

We present confirmation of HD 143811 AB b, a substellar companion to spectroscopic binary HD 143811 AB through direct imaging with the Gemini Planet Imager (GPI) and Keck NIRC2. HD 143811 AB was observed as a part of the GPI Exoplanet Survey in 2016 and 2019 and is a member of the Sco-Cen star formation region. The exoplanet is detected ∼430 mas from the host star by GPI. With two GPI epochs and one from Keck/NIRC2 in 2022, we confirm through common proper motion analysis that the object is bound to its host star. We derive an orbit with a semimajor axis of $64^{+32}_{-14}$ au and eccentricity $0.23^{+0.24}_{-0.16}$. Spectral analysis of the GPI H-band spectrum and NIRC2 L′ photometry provides additional proof that this object is a substellar companion. We compare the spectrum of HD 143811 AB b to PHOENIX stellar models and Exo-Radioactive-Convective Equilibrium Model (REM) exoplanet atmosphere models and find that Exo-REM models provide the best fits to the data. From the Exo-REM models, we derive an effective temperature of $1042^{+178}_{-132}$ K for the planet and translate the derived luminosity of the planet to a mass of 5.6 ± 1.1 M Jup assuming hot-start evolutionary models. HD 143811 AB b is the first directly imaged planet around a binary that is not on an ultrawide orbit. Future characterization of this object will shed light on the formation of planets around binary star systems.

Astronomy and AstroPhysics↗

DOE OSTI2022

MaNGA 8313-1901: Gas Accretion Observed in a Blue Compact Dwarf Galaxy?

Gas accretion is an important process in the evolution of galaxies, but it has limited direct observational evidences. In this paper, we report the detection of a possible ongoing gas accretion event in a blue compact dwarf (BCD) galaxy, MaNGA 8313-1901, observed by the Mapping Nearby Galaxies and Apache Point Observatory (MaNGA) program. This galaxy has a distinct off-centered blue clump to the northeast (the NE clump) that shows low metallicity and enhanced star formation. The kinematics of the gas in the NE clump also seems to be detached from the host BCD galaxy. Together with the metallicity drop of the NE clump, it suggests that the NE clump likely has an external origin, such as gas accretion or galaxy interaction, rather than an internal origin, such as an H II complex in the disk. After removing the underlying host component, we find that the spectrum of the "pure" clump can match very well with a modeled spectrum containing a stellar population of the young stars (≤7 Myr) only. This may imply that the galaxy is experiencing an accretion of cold gas, instead of a merger event involving galaxies with significant preexisting old stars. We also find signs of another clump (the SW clump) at the southwest corner of the host galaxy, and the two clumps may share the same origin of gas accretion.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

3D turbulent reconnection: Theory, tests, and astrophysical implications

Magnetic reconnection, topological changes in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical fluid flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence, magnetic field lines constantly reconnect everywhere and on all scales, thus making magnetic reconnection an intrinsic part of the turbulent cascade. We note in particular that this is inconsistent with the usual practice of magnetic field lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies starting with the paper done by Lazarian and Vishniac [Astrophys. J. 517, 700–718 (1999)] proposed that 3D turbulence makes magnetic reconnection fast and that magnetic reconnection and turbulence are intrinsically connected. In particular, we discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence. We demonstrate that in the presence of turbulence, the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. The latter fact justifies a magnetohydrodynamiclike treatment of magnetic reconnection on all scales much larger than the relevant plasma scales. We discuss the numerical and observational evidence supporting the turbulent reconnection model. In particular, we demonstrate that the tearing reconnection is suppressed in 3D, and unlike the 2D settings, 3D reconnection induces turbulence that makes magnetic reconnection independent of resistivity. We show that turbulent reconnection dramatically affects key astrophysical processes, e.g., star formation, turbulent dynamo, and acceleration of cosmic rays. We provide criticism of the concept of “reconnection-mediated turbulence” and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity and why the latter have fatal conceptual flaws.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Quenching in the Right Place at the Right Time: Tracing the Shared History of Starbursts, Active Galactic Nuclei, and Poststarburst Galaxies Using Their Structures and Multiscale Environments

This work uses multiscale environments and the structures of galaxies in the Sloan Digital Sky Survey as consistency checks on the evolution from starburst to quiescence at redshift z < 0.2. The environmental indicators include fixed-aperture mass overdensities (δ xMpc , x ϵ {0.5, 1, 2, 4, 8} h –1 Mpc), k-nearest neighbor distances, the tidal parameter, halo mass (M h ), and satellite/central classification. The residuals of the specific star formation rates (Δ SSFRs) are used to select starbursts (ΔSSFR > 0.6 dex, N ≈ 8600). Quenched poststarbursts (QPSBs) are selected using Hα < 3 Å emission and Hδ A > 4 Å absorption (N ≈ 750). The environments of starbursts and QPSBs are compared with those of active galactic nuclei (AGNs) and inactive galaxies of varying ΔSSFR. The environments of starbursts, AGNs, and QPSBs are unlike the environments of most quiescent galaxies (QGs). About 70%–90% of starbursts, AGNs with Hδ A > 4 Å, and QPSBs are centrals; ~80%–90% have M h < 10 13 M ⊙ , and only ~2%–4% have M h > 10 14 M ⊙ or live in clusters. Their M h and satellite fractions are also different from those of QGs. All QPSBs are matched to some SFGs, starbursts, AGNs, and QGs of similar M $\star$ , environments, concentration indices, and velocity dispersions. A significant fraction (~20%–30%) of starbursts cannot be matched to QPSBs or QGs. The implications are that: (1) some starbursts do not quench rapidly; (2) the satellite-quenching mechanisms operating in high-density environments cannot account for most QPSBs; (3) the evolution from starbursts to QPSBs to QGs is not the dominant path at z < 0.2; and (4) starbursts are not mainly triggered by tidal interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Galaxy zoo: stronger bars facilitate quenching in star-forming galaxies

ABSTRACT We have used Galaxy Zoo DECaLS (GZD) to study strong and weak bars in disc galaxies. Out of the 314 000 galaxies in GZD, we created a volume-limited sample (0.01 < z < 0.05, Mr < − 18.96) which contains 1867 galaxies with reliable volunteer bar classifications in the ALFALFA footprint. In keeping with previous Galaxy Zoo surveys (such as GZ2), the morphological classifications from GZD agree well with previous morphological surveys. GZD considers galaxies to either have a strong bar (15.5 per cent), a weak bar (28.1 per cent) or no bar (56.4 per cent), based on volunteer classifications on images obtained from the DECaLS survey. This places GZD in a unique position to assess differences between strong and weak bars. We find that the strong bar fraction is typically higher in quiescent galaxies than in star-forming galaxies, while the weak bar fraction is similar. Moreover, we have found that strong bars facilitate the quenching process in star-forming galaxies, finding higher fibre star formation rates (SFRs), lower gas masses, and shorter depletion time-scales in these galaxies compared to unbarred galaxies. However, we also found that any differences between strong and weak bars disappear when controlling for bar length. Based on this, we conclude that weak and strong bars are not fundamentally different phenomena. Instead, we propose that there is a continuum of bar types, which varies from ‘weakest’ to ‘strongest’.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Observational Evidence for Cosmological Coupling of Black Holes and its Implications for an Astrophysical Source of Dark Energy

Observations have found black holes spanning 10 orders of magnitude in mass across most of cosmic history. The Kerr black hole solution is, however, provisional as its behavior at infinity is incompatible with an expanding universe. Black hole models with realistic behavior at infinity predict that the gravitating mass of a black hole can increase with the expansion of the universe independently of accretion or mergers, in a manner that depends on the black hole's interior solution. We test this prediction by considering the growth of supermassive black holes in elliptical galaxies over 0 < z ≲ 2.5. We find evidence for cosmologically coupled mass growth among these black holes, with zero cosmological coupling excluded at 99.98% confidence. The redshift dependence of the mass growth implies that, at z ≲ 7, black holes contribute an effectively constant cosmological energy density to Friedmann's equations. The continuity equation then requires that black holes contribute cosmologically as vacuum energy. We further show that black hole production from the cosmic star formation history gives the value of Ω Λ measured by Planck while being consistent with constraints from massive compact halo objects. We thus propose that stellar remnant black holes are the astrophysical origin of dark energy, explaining the onset of accelerating expansion at z ~ 0.7.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Correlation between the gas-phase metallicity and ionization parameter in extragalactic H II regions

The variations of the metallicity and ionization parameter in H II regions are usually thought to be the dominant factors that produce the variations we see in the observed emission line spectra. There is an increasing amount of evidence that these two quantities are physically correlated, although the exact form of this correlation is debatable in the literature. Simulated emission line spectra from photoionized clouds provide important clues about the physical conditions of H II regions and are frequently used for deriving metallicities and ionization parameters. Through a systematic investigation on the assumptions and methodology used in applying photoionization models, we find that the derived correlation has a strong dependence on the choice of model parameters. On the one hand, models that give consistent predictions over multiple emission-line ratios yield a positive correlation between the metallicity and ionization parameter for the general population of H II regions or star-forming galaxies. On the other hand, models that are inconsistent with the data locus in high-dimensional line ratio space yield discrepant correlations when different subsets of line ratios are used in the derivation. The correlation between the metallicity and ionization parameter has a secondary dependence on the surface density of the star formation rate (SFR), with the higher SFR regions showing a higher ionization parameter but weaker correlations. The existence of the positive correlation contradicts the analytical wind-driven bubble model for H II regions. We explore assumptions in both dynamical models and photoionization models, and conclude that there is a potential bias associated with the geometry. However, this is still insufficient to explain the correlation. Mechanisms that suppress the dynamical influence of stellar winds in realistic H II regions might be the key to solving this puzzle, though more sophisticated combinations of dynamical models and photoionization models to test are required.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

What drives galaxy quenching? A deep connection between galaxy kinematics and quenching in the local Universe

ABSTRACT We develop a 2D inclined rotating disc model, which we apply to the stellar velocity maps of 1862 galaxies taken from the MaNGA survey (SDSS public Data Release 15) . We use a random forest classifier to identify the kinematic parameters that are most connected to galaxy quenching. We find that kinematic parameters that relate predominantly to the disc (such as the mean rotational velocity) and parameters that characterize whether a galaxy is rotation- or dispersion-dominated (such as the ratio of rotational velocity to velocity dispersion) are not fundamentally linked to the quenching of star formation. Instead, we find overwhelmingly that it is the absolute level of velocity dispersion (a property that relates primarily to a galaxy’s bulge/spheroidal component) that is most important for separating star-forming and quenched galaxies. Furthermore, a partial correlation analysis shows that many commonly discussed correlations between galaxy properties and quenching are spurious, and that the fundamental correlation is between quenching and velocity dispersion. In particular, we find that at fixed velocity dispersion, there is only a very weak dependence of quenching on the disc properties, whereby more discy galaxies are slightly more likely to be forming stars. By invoking the tight relationship between black hole mass and velocity dispersion, and noting that black hole mass traces the total energy released by active galactic nuclei (AGNs), we argue that these data support a scenario in which quenching occurs by preventive feedback from AGNs. The kinematic measurements from this work are publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Warped Disk Galaxies. I. Linking U-type Warps in Groups/Clusters to Jellyfish Galaxies

Warped disk galaxies are classified into two morphologies: S and U types. Conventional theories routinely attribute both types to galactic tidal interaction and/or gas accretion, but reproducing U types in simulations is extremely challenging. Here we investigate whether both types are governed by the same mechanisms using the most extensive sample of ~8000 nearby (0.02 < z < 0.06) massive (M*/M ⊙ > 10 9 ) edge-on disks from the Sloan Digital Sky Survey. We find that U types show on average bluer optical colors and a higher specific star formation rate (sSFR) than S types, with more strongly warped U types having a higher sSFR. We also find that while the S-type warp properties correlate with the tidal force by the nearest neighbor regardless of the environment, there is no such correlation for U types in groups/clusters, suggesting a nontidal mechanism could be at play for U types, such as ram pressure stripping (RPS). Indeed, U types are more common in groups/clusters than in fields and they have stellar mass, gas fraction, sSFR enhancement, and phase-space distribution closely analogous to RPS-induced jellyfish galaxies in clusters. We furthermore show that the stellar disks of most RPS galaxies in the IllustrisTNG simulation are warped in a U shape and bent in the opposite direction of stripped gas tails, satisfying theoretical expectations for stellar warps embedded in jellyfishes. We therefore suggest that despite the majority of U types that live in fields being still less explained, RPS can be an alternative origin for those in groups/clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Spatially resolved mass–metallicity relation at $z$ ~ 0.26 from the MUSE-Wide Survey

Galaxies in the local universe have a spatially resolved star-forming main sequence (rSFMS) and mass–metallicity relation (rMZR). We know that the global mass–metallicity relation (MZR) results from the integral of the rMZR, and it evolves with redshift. However, the evolution of the rMZR with redshift is still unclear because the spatial resolution and signal-to-noise ratio are low. Currently, too few observations beyond the local universe are available, and only simulations can reproduce the evolution of the rMZR with redshift. We selected ten emission-line galaxies with an average redshift of $z$ ~ 0.26 from the MUSE-Wide DR1. We obtained the spatially resolved star formation rate (SFR) and metallicity from integral field spectroscopy (IFS), as well as the stellar mass surface density from 3D-HST photometry. We derived the rSFMS and rMZR at $z$ ~ 0.26 and compared them with those of local galaxies. We find that the rSFMS of galaxies at $z$ ~ 0.26 has a slope of ~0.771. The rMZR exists at z ~ 0.26, showing a similar shape to that of the local universe, but a lower average metallicity that is about ~0.11 dex lower than the local metallicity. In addition, we also study the spatially resolved fundamental metallicity relation (rFMR) of these galaxies. However, there is no obvious evidence that an rFMR exists at $z$ ~ 0.26, and it is not an extension of rMZR at a high SFR. Similar to their global versions, the rSFMS and rMZR of galaxies also evolve with redshift. At fixed stellar mass, galaxies at higher redshift show a higher SFR and lower metallicity. In conclusion, these suggest that the evolution of the global galaxy properties with redshift may result from integrating the evolution of the spatially resolved galaxy properties.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

The most luminous, merger-free AGNs show only marginal correlation with bar presence

The role of large-scale bars in the fuelling of active galactic nuclei (AGNs) is still debated, even as evidence mounts that black hole growth in the absence of galaxy mergers cumulatively dominates and may substantially influence disc (i.e. merger-free) galaxy evolution. We investigate whether large-scale galactic bars are a good candidate for merger-free AGN fuelling. Specifically, we combine slit spectroscopy and Hubble Space Telescope imagery to characterize star formation rates (SFRs) and stellar masses of the unambiguously disc-dominated host galaxies of a sample of luminous, Type-1 AGN with 0.02 < z < 0.24. After carefully correcting for AGN signal, we find no clear difference in SFRs between AGN hosts and a stellar mass-matched sample of galaxies lacking an AGN (0.013 < z < 0.19), although this could be due to small sample size (nAGN = 34). We correct for SFR and stellar mass to minimize selection biases, and compare the bar fraction in the two samples. We find that AGNs are marginally (~1.7σ) more likely to host a bar than inactive galaxies, with AGN hosts having a bar fraction, $f_{\mathrm{bar}}=0.59^{+0.08}_{-0.09}$ and inactive galaxies having a bar fraction, $f_{\mathrm{bar}}=0.44^{+0.08}_{-0.09}$. However, we find no further differences between SFR- and mass-matched AGNs and inactive samples. While bars could potentially trigger AGN activity, they appear to have no further, unique effect on a galaxy’s stellar mass or SFR.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

FORECASTOR – II. Simulating galaxy surveys with the Cosmological Advanced Survey Telescope for Optical and UV Research

The Cosmological Advanced Survey Telescope for Optical and UV Research (CASTOR) is a planned flagship space telescope, covering the blue-optical and UV part of the spectrum. Here, we introduce the CASTOR image simulator, a python GalSim package-based script capable of generating mock CASTOR images from an input catalogue. We generate example images from the CASTOR Wide, Deep, and Ultra-Deep surveys using simulated lightcones from the Santa Cruz semi-analytic model. We make predictions for the performance of these surveys by comparing galaxies that are extracted from each image using Source Extractor to the input catalogue. We find that the Wide, Deep, and Ultra-Deep surveys will be 75 per cent complete for point sources down to $\sim 27$, 29, and 30 mag, respectively, in the UV, u, and g filters, with the UV-split and u-split filters reaching a shallower depth. With a large area of $\sim 2200$ deg$^2$, the Wide survey will detect hundreds of millions of galaxies out to $z\sim 4$, mostly with $M_\ast \gtrsim 10^{9}\,{\rm M}_{\odot }$. The Ultra-Deep survey will probe to $z\sim 5$, detecting galaxies with $M_\ast \gtrsim 10^{7}{\rm M}_{\odot }$. These galaxy samples will enable precision measurements of the distribution of star formation in the cosmic web, connecting the growth of stellar mass to the assembly of dark matter haloes over two thirds of the history of the Universe, and other core goals of CASTOR’s legacy surveys. These image simulations and the tools developed to generate them will be a vital planning tool to estimate CASTOR’s performance and iterate the telescope and survey designs prior to launch.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Can we observe the ion-neutral drift velocity in prestellar cores?

Given the low-ionization fraction of molecular clouds, ambipolar diffusion is thought to be an integral process in star formation. However, chemical and radiative-transfer effects, observational challenges, and the fact that the ion-neutral drift velocity is inherently very small render a definite detection of ambipolar diffusion extremely non-trivial. Here, we study the ion-neutral drift velocity in a suite of chemodynamical, non-ideal magnetohydrodynamic (MHD), two-dimensional axisymmetric simulations of prestellar cores where we alter the temperature, cosmic-ray ionization rate, visual extinction, mass-to-flux ratio, and chemical evolution. Subsequently, we perform a number of non-local thermodynamic equilibrium (non-LTE) radiative-transfer calculations considering various idealized and non-idealized scenarios in order to assess which factor (chemistry, radiative transfer, and/or observational difficulties) is the most challenging to overcome in our efforts to detect the ion-neutral drift velocity. We find that temperature has a significant effect in the amplitude of the drift velocity with the coldest modelled cores (T = 6 K) exhibiting drift velocities comparable to the sound speed. Against expectations, we find that in idealized scenarios (where two species are perfectly chemically co-evolving) the drift velocity ‘survives’ radiative-transfer effects and can in principle be observed. However, we find that observational challenges and chemical effects can significantly hinder our view of the ion-neutral drift velocity. Finally, we propose that HCN and HCNH + ⁠, being chemically co-evolving, could be used in future observational studies aiming to measure the ion-neutral drift velocity.

79 ASTRONOMY AND ASTROPHYSICS↗