skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, July 11 until 2:00 AM ET on Saturday, July 12 due to maintenance. We apologize for the inconvenience.


Title: Fitting Methods for Probability Distribution Functions in Turbulent Star-forming Clouds
Abstract We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback and magnetic fields, to investigate the effectiveness of different fitting methods for volume and column density probability distribution functions (PDFs). The first method fits a piecewise lognormal and power-law (PL) function to recover PDF parameters such as the PL slope and transition density. The second method fits a polynomial spline function and examines the first and second derivatives of the spline to determine the PL slope and the functional transition density. The first PL (set by the transition between lognormal and PL function) can also be visualized in the derivatives directly. In general, the two methods produce fits that agree reasonably well for volume density but vary for column density, likely due to the increased statistical noise in the column density PDFs as compared to the volume density PDFs. We test a well-known conversion for estimating volume density PL slopes from column density slopes and find that the spline method produces a better match (χ2of 3.34 versusχ2of 5.92), albeit with a significant scatter. Ultimately, we recommend the use of both fitting methods on column density data to mitigate the effects of noise.  more » « less
Award ID(s):
2401740
PAR ID:
10567045
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
979
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 89
Size(s):
Article No. 89
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT The velocity gradients technique (VGT) and the probability density functions (PDFs) of mass density are tools to study turbulence, magnetic fields, and self-gravity in molecular clouds. However, self-absorption can significantly make the observed intensity different from the column density structures. In this work, we study the effects of self-absorption on the VGT and the intensity PDFs utilizing three synthetic emission lines of CO isotopologues 12CO (1–0), 13CO (1–0), and C18O (1–0). We confirm that the performance of VGT is insensitive to the radiative transfer effect. We numerically show the possibility of constructing 3D magnetic fields tomography through VGT. We find that the intensity PDFs change their shape from the pure lognormal to a distribution that exhibits a power-law tail depending on the optical depth for supersonic turbulence. We conclude the change of CO isotopologues’ intensity PDFs can be independent of self-gravity, which makes the intensity PDFs less reliable in identifying gravitational collapsing regions. We compute the intensity PDFs for a star-forming region NGC 1333 and find the change of intensity PDFs in observation agrees with our numerical results. The synergy of VGT and the column density PDFs confirms that the self-gravitating gas occupies a large volume in NGC 1333. 
    more » « less
  2. Abstract We perform X-ray spectral analyses to derive the characteristics (e.g., column density, X-ray luminosity) of ≈10,200 active galactic nuclei (AGNs) in the XMM-Spitzer Extragalactic Representative Volume Survey, which was designed to investigate the growth of supermassive black holes across a wide dynamic range of cosmic environments. Using physical torus models (e.g., Borus02) and a Bayesian approach, we uncover 22 representative Compton-thick (CT;NH> 1.5 × 1024cm−2) AGN candidates with good signal-to-noise ratios as well as a large sample of 136 heavily obscured AGNs. We also find an increasing CT fraction (fCT) from low (z< 0.75) to high (z> 0.75) redshift. Our CT candidates tend to show hard X-ray spectral shapes and dust extinction in their spectral energy distribution fits, which may shed light on the connection between AGN obscuration and host-galaxy evolution. 
    more » « less
  3. A<sc>bstract</sc> We study the interactions of systems of two and three nondegenerate mesons composed of pions and kaons at maximal isospin using lattice QCD, specificallyπ+K++π+K+andK+K+π+. Utilizing the stochastic LapH method, we determine the spectrum of these systems on two CLSNf= 2 + 1 ensembles with pion masses of 200 MeV and 340 MeV, and include many levels in different momentum frames. We constrain the K matrices describing two- and three-particle interactions by fitting the spectrum to the results predicted by the finite-volume formalism, including up topwaves. This requires also results for theπ+π+andK+K+spectrum, which have been obtained previously on the same configurations. We explore different fitting strategies, comparing fits to energy shifts with fits to energies boosted to the rest frame, and also comparing simultaneous global fits to all relevant two- and three-particle channels to those where we first fit two-particle channels and then add in the three-particle information. We provide the first determination of the three-particle K matrix inπ+π+K+andK+K+π+systems, finding statistically significant nonzero results in most cases. We includesandpwaves in the K matrix forπ+K+scattering, finding evidence for an attractivep-wave scattering length. We compare our results to Chiral Perturbation Theory, including an investigation of the impact of discretization errors, for which we provide the leading order predictions obtained using Wilson Chiral Perturbation Theory. 
    more » « less
  4. Abstract We use the known surface boulder-size distribution of the C-type rubble pile asteroid Ryugu (NEA 162173) to determine its macroporosity, assuming it is a homogeneous granular aggregate. We show that the volume-frequency distribution of its boulders, cobbles, and pebbles, is well-represented by a lognormal function withσ= 2.4 ± 0.1 andμ= 0.2 ± 0.05. Application of linear-mixture packing theory yields a value for the macroporosity ofϕ= 0.14 ± 0.04. Given its low bulk density of 1.19 gm cm−3, this implies an average density for Ryugu’s rocks of 1.38 ± 0.07 gm cm−3throughout its volume, consistent with a recent determination for surface boulders based on their thermal properties. This supports the spectrum-based argument that interplanetary dust particles may be the best analog material available on Earth, and it suggests that high-density, well-lithified objects such as chondrules and chondrule-bearing chondrites may be rare on Ryugu. Implications of this result for the origin of chondrules, a long-standing problem in cosmochemistry, are discussed. We propose that chondrules and most chondrites formed together in rare lithification events, which occurred during the accretion of chondritic envelopes to large, differentiated planetesimals at a time when they were still hot from26Al decay. 
    more » « less
  5. Abstract Reconstructing the history of polar temperature from ice core water isotope (δ18O) calibration has remained a challenge in paleoclimate research, because of our incomplete understanding of various temperature–δ18O relationships. This paper resolves this classical problem in a new framework called the unified slope equations (USE), which illustrates the general relations among spatial and temporalδ18O–surface temperature slopes. The USE is applied to the Antarctica temperature change during the last deglaciation in model simulations and observations. It is shown that the comparable Antarctica-mean spatial slope with deglacial temporal slope inδ18O–surface temperature reconstruction is caused, accidentally, by the compensation responses between theδ18O–inversion layer temperature relation and the inversion layer temperature itself. Furthermore, in light of the USE, we propose that the present seasonal slope ofδ18O–inversion layer temperature is an optimal paleothermometer that is more accurate and robust than the spatial slope. This optimal slope suggests the possibility of reconstructing past Antarctic temperature changes using present and future instrumental observations. Significance StatementThis paper develops a new framework called the unified slope equations (USE) to provide, for the first time, a general relation among various spatial and temporal water isotope–temperature slopes. The application of the USE to Antarctic deglacial temperature change shows that the optimal paleothermometer is the seasonal slope of the inversion layer temperature. 
    more » « less