Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stellar properties”

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.

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

TOI–757 b: an eccentric transiting mini–Neptune on a 17.5–d orbit

ABSTRACT We report the spectroscopic confirmation and fundamental properties of TOI$-$757 b, a mini$-$Neptune on a 17.5$-$d orbit transiting a bright star ($V\, =\, 9.7$ mag) discovered by the TESS mission. We acquired high$-$precision radial velocity measurements with the HARPS, ESPRESSO, and PFS spectrographs to confirm the planet detection and determine its mass. We also acquired space$-$borne transit photometry with the CHEOPS space telescope to place stronger constraints on the planet radius, supported with ground$-$based LCOGT photometry. WASP and KELT photometry were used to help constrain the stellar rotation period. We also determined the fundamental parameters of the host star. We find that TOI$-$757 b has a radius of $R_{\mathrm{p}} = 2.5 \pm 0.1 R_{\oplus }$ and a mass of $M_{\mathrm{p}} = 10.5^{+2.2}_{-2.1} M_{\oplus }$, implying a bulk density of $\rho _{\text{p}} = 3.6 \pm 0.8$ g cm$^{-3}$. Our internal composition modelling was unable to constrain the composition of TOI$-$757 b, highlighting the importance of atmospheric observations for the system. We also find the planet to be highly eccentric with e = 0.39$^{+0.08}_{-0.07}$, making it one of the very few highly eccentric planets among precisely characterized mini$-$Neptunes. Based on comparisons to other similar eccentric systems, we find a likely scenario for TOI$-$757 b’s formation to be high eccentricity migration due to a distant outer companion. We additionally propose the possibility of a more intrinsic explanation for the high eccentricity due to star$-$star interactions during the earlier epoch of the Galactic disc formation, given the low metallicity and older age of TOI$-$757.

Alqasim, A. (ORCID:0000000151025505)↗
DOE OSTI2024

Constraining neutrino-DM interactions with Milky Way dwarf spheroidals and supernova neutrinos

We constrain the neutrino-dark matter cross section using properties of the dark matter density profiles of Milky Way dwarf spheroidal galaxies. The constraint arises from core-collapse supernova neutrinos scattering on dark matter as a form of energy injection, allowing the transformation of the dark matter density profile from a cusped profile to a flatter profile. We assume a standard cosmology of dark energy and cold, collisionless, and non-self-interacting dark matter. By requiring that the dark matter cores do not lose too much mass or overshoot constraints from stellar kinematics, we place an upper limit on the cross section of σ ν − DM ( E ν = 15 MeV , m χ ≲ 130 GeV ) ≈ 3.4 × 10 − 23 cm 2 and σ ν − DM ( E ν = 15 MeV , m χ ≳ 130 GeV ) ≈ 3.2 × 10 − 27 ( m χ 1 GeV ) 2 cm 2 , which is stronger than previous bounds for these energies. Consideration of baryonic feedback or host galaxy effects on the dark matter profile can strengthen this constraint. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Multiple Explanations for the Single Transit of KIC 5951458 Based on Radial Velocity Measurements Extracted with a Novel Matched-template Technique

Planetary systems that show single-transit events are a critical pathway to increasing the yield of long-period exoplanets from transit surveys. From the primary Kepler mission, KIC 5951458 b (Kepler-456b) was thought to be a single-transit giant planet with an orbital period of 1310 days. However, radial velocity (RV) observations of KIC 5951458 from the HIRES instrument on the Keck telescope suggest that the system is far more complicated. To extract precise RVs for this V≈13 star, we develop a novel matched-template technique that takes advantage of a broad library of template spectra acquired with HIRES. We validate this technique and measure its noise floor to be 4–8 m s{sup −1} (in addition to internal RV error) for most stars that would be targeted for precision RVs. For KIC 5951458, we detect a long-term RV trend that suggests the existence of a stellar companion with an orbital period greater than a few thousand days. We also detect an additional signal in the RVs that is possibly caused by a planetary or brown dwarf companion with mass in the range of 0.6–82 M {sub Jup} and orbital period below a few thousand days. Curiously, from just the data on hand, it is not possible to determine which object caused the single “transit” event. We demonstrate how a modest set of RVs allows us to update the properties of this unusual system and predict the optimal timing for future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

On the Estimation of Circumbinary Orbital Properties

We describe a fast, approximate method to characterize the orbits of satellites around a central binary in numerical simulations. A goal is to distinguish the free eccentricity—the random motion of a satellite relative to a dynamically cool orbit—from oscillatory modes driven by the central binary’s time-varying gravitational potential. We assess the performance of the method using the Kepler-16, Kepler-47, and Pluto–Charon systems. We then apply the method to a simulation of orbital damping in a circumbinary environment, resolving relative speeds between small bodies that are slow enough to promote mergers and growth. These results illustrate how dynamical cooling can set the stage for the formation of Tatooine-like planets around stellar binaries and the small moons around the Pluto–Charon binary planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Rates and properties of Type Ia supernovae in galaxy clusters within the dark energy survey

We identify 66 photometrically classified Type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) that have occurred within red-sequence selected galaxy clusters. We compare light-curve and host galaxy properties of the cluster SNe to 1024 DES SNe Ia located in field galaxies, the largest comparison of two such samples at high redshift (z > 0.1). We find that cluster SN light curves decline faster than those in the field (97.7 per cent confidence). However, when limiting these samples to host galaxies of similar colour and mass, there is no significant difference in the SN light-curve properties. Motivated by previous detections of a higher-normalized SN Ia delay-time distribution in galaxy clusters, we measure the intrinsic rate of SNe Ia in cluster and field environments. We find the average ratio of the SN Ia rate per galaxy between high-mass (⁠|$10\le \log \mathrm{(\mathit{ M}_{*}/{\rm \mathit{ M}}_{\odot })} \le 11.25$|⁠) cluster and field galaxies to be 0.594 ± 0.068. This difference is mass-dependent, with the ratio declining with increasing mass, which suggests that the stellar populations in cluster hosts are older than those in field hosts. We show that the mass-normalized rate (or SNe per unit mass) in massive–passive galaxies is consistent between cluster and field environments. Additionally, both of these rates are consistent with rates previously measured in clusters at similar redshifts. We conclude that in massive–passive galaxies, which are the dominant hosts of cluster SNe, the cluster delay-time distribution is comparable to the field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Avalanche statistics of fluctuation-induced fluxes from the SLPM and the W7-AS stellarator

Measurements of fluctuating floating potentials and ion saturation currents at different radial locations in the Santander Linear Plasma Machine (Castellanos et al 2005 Plasma Phys. Control. Fusion47 2067) and at the edge of the W7-AS stellarator by means of radially movable Langmuir probes allow to infer the corresponding fluctuation-induced radial flux temporal series. Avalanche-like transport events are identified in the time series and statistically characterized in terms of avalanche size/duration/quiet-time distributions and size-duration scaling relations. Transport is diffusive in the inner and intermediate radial region of the SLPM r < r tr ≈ 2.6 cm, undergoing a transition at r tr , becoming non-diffusive in the outermost region of the device, r > r tr . Here, the results obtained at the edge of the W7-AS stellarator are similar to those found in SLPM for r > r tr , i.e. consistent with what would be expected for scale-free, self-similar plasma transport dynamics near a critical state.

Avalanches↗
DOE OSTI2022

Substructure in the stellar halo near the Sun: I. Data-driven clustering in integrals-of-motion space

Context. Merger debris is expected to populate the stellar haloes of galaxies. In the case of the Milky Way, this debris should be apparent as clumps in a space defined by the orbital integrals of motion of the stars. Aims. Our aim is to develop a data-driven and statistics-based method for finding these clumps in integrals-of-motion space for nearby halo stars and to evaluate their significance robustly. Methods. We used data from Gaia EDR3, extended with radial velocities from ground-based spectroscopic surveys, to construct a sample of halo stars within 2.5 kpc from the Sun. We applied a hierarchical clustering method that makes exhaustive use of the single linkage algorithm in three-dimensional space defined by the commonly used integrals of motion energy E, together with two components of the angular momentum, L z and L ⊥ . To evaluate the statistical significance of the clusters, we compared the density within an ellipsoidal region centred on the cluster to that of random sets with similar global dynamical properties. By selecting the signal at the location of their maximum statistical significance in the hierarchical tree, we extracted a set of significant unique clusters. By describing these clusters with ellipsoids, we estimated the proximity of a star to the cluster centre using the Mahalanobis distance. Additionally, we applied the HDBSCAN clustering algorithm in velocity space to each cluster to extract subgroups representing debris with different orbital phases. Results. Our procedure identifies 67 highly significant clusters (> 3σ), containing 12% of the sources in our halo set, and 232 subgroups or individual streams in velocity space. In total, 13.8% of the stars in our data set can be confidently associated with a significant cluster based on their Mahalanobis distance. Inspection of the hierarchical tree describing our data set reveals a complex web of relations between the significant clusters, suggesting that they can be tentatively grouped into at least six main large structures, many of which can be associated with previously identified halo substructures, and a number of independent substructures. This preliminary conclusion is further explored in a companion paper, in which we also characterise the substructures in terms of their stellar populations. Conclusions. Our method allows us to systematically detect kinematic substructures in the Galactic stellar halo with a data-driven and interpretable algorithm. The list of the clusters and the associated star catalogue are provided in two tables available at the CDS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Use of current-potential patches to obtain fundamental improvements to the coils of magnetic fusion devices

A central issue in the design of tokamaks or stellarators is the coils that produce the external magnetic fields. The freedom that remains unstudied in the design of coils is enormous. This freedom could be quickly studied computationally at low cost with high reliability. In particular, the space between toroidal field and modular coils that blocks access to the plasma chamber could be increased by a large factor. The concept of current-potential patches, which was developed in Todd Elder's thesis, provides a method for separating the study of the feasibility of coils with attractive physics properties from the engineering design of specific coils. In addition to enhanced accessibility, coils can be designed for increased plasma-coil separation, insensitivity to coil position errors, and plasma control.

Boozer, Allen H. [Columbia Univ., New York, NY (Un↗
DOE OSTI2021

Atacama Cosmology Telescope: Modeling the gas thermodynamics in BOSS CMASS galaxies from kinematic and thermal Sunyaev-Zel’dovich measurements

The thermal and kinematic Sunyaev-Zel'dovich effects (tSZ, kSZ) probe the thermodynamic properties of the circumgalactic and intracluster medium (CGM and ICM) of galaxies, groups, and clusters, since they are proportional, respectively, to the integrated electron pressure and momentum along the line-of-sight. We present constraints on the gas thermodynamics of CMASS galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) using new measurements of the kSZ and tSZ signals obtained in a companion paper. Combining kSZ and tSZ measurements, we measure within our model the amplitude of energy injection $\epsilon M_\star c^2$, where $M_\star$ is the stellar mass, to be $\epsilon=(40\pm9)\times10^{-6}$, and the amplitude of the non-thermal pressure profile to be $\alpha_{\rm Nth}<0.2$ (2$\sigma$), indicating that less than 20% of the total pressure within the virial radius is due to a non-thermal component. We estimate the effects of including baryons in the modeling of weak-lensing galaxy cross-correlation measurements using the best fit density profile from the kSZ measurement. Our estimate reduces the difference between the original theoretical model and the weak-lensing galaxy cross-correlation measurements in arXiv:1611.08606 by half, but does not fully reconcile it. Comparing the kSZ and tSZ measurements to cosmological simulations, we find that they under predict the CGM pressure and to a lesser extent the CGM density at larger radii. This suggests that the energy injected via feedback models in the simulations that we compared against does not sufficiently heat the gas at these radii. We do not find significant disagreement at smaller radii. These measurements provide novel tests of current and future simulations. Here, this work demonstrates the power of joint, high signal-to-noise kSZ and tSZ observations, upon which future cross-correlation studies will improve.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Velocity Map Asymmetry of Ionized Gas in MaNGA. I. The Catalog and General Properties

The SDSS-IV MaNGA survey has measured two-dimensional maps of emission-line velocities for a statistically powerful sample of nearby galaxies. The asymmetric features of these kinematics maps reflect the nonrotational component of a galaxy’s internal motion of ionized gas. In this study, we present a catalog of kinematic asymmetry measurements of the H α velocity map of a sample of 5353 MaNGA galaxies. Based on this catalog, we find that “special” galaxies (e.g., merging galaxies, barred galaxies, and active galactic nucleus host galaxies) contain more galaxies with highly asymmetric velocity maps. However, we notice that more than half of galaxies with high kinematic asymmetry in our sample are quite “regular.” For those “regular” galaxies, kinematic asymmetry shows a significant anticorrelation with stellar mass at log M * < 9.7 , while such a trend becomes very weak at log M * > 9.7 . Moreover, at a given stellar mass, the kinematic asymmetry shows weak correlations with photometric morphology, star formation rate, and environment, while it is independent of H i gas content. We also have quantified the observational effects in the kinematic asymmetry measurement. We find that both the signal-to-noise ratio of H α flux and disk inclination angle contribute to the measures of kinematic asymmetry, while the physical spatial resolution is an irrelevant factor inside the MaNGA redshift coverage.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The cosmic DANCe of Perseus

Context. Star-forming regions are excellent benchmarks for testing and validating theories of star formation and stellar evolution. The Perseus star-forming region, being one of the youngest (< 10 Myr), closest (280-320 pc), and most studied in the literature, is a fundamental benchmark. Aims. We aim to study the membership, phase-space structure, mass, and energy (kinetic plus potential) distribution of the Perseus star-forming region using public catalogues (Gaia, APOGEE, 2MASS, and Pan-STARRS). Methods. We used Bayesian methodologies that account for extinction to identify the Perseus physical groups in the phase-space, retrieve their candidate members, derive their properties (age, mass, 3D positions, 3D velocities, and energy), and attempt to reconstruct their origin. Results. We identify 1052 candidate members in seven physical groups (one of them new) with ages between 3 and 10 Myr, dynamical super-virial states, and large fractions of energetically unbounded stars. Their mass distributions are broadly compatible with that of Chabrier for masses ≳0.1 M ⊙ and do not show hints of over-abundance of low-mass stars in NGC 1333 with respect to IC 348. These groups’ ages, spatial structure, and kinematics are compatible with at least three generations of stars. Future work is still needed to clarify if the formation of the youngest was triggered by the oldest. Conclusions. The exquisite Gaia data complemented with public archives and mined with comprehensive Bayesian methodologies allow us to identify 31% more members than previous studies, discover a new physical group (Gorgophone: 7 Myr, 191 members, and 145 M ⊙ ), and confirm that the spatial, kinematic, and energy distributions of these groups support the hierarchical star formation scenario.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Cluster environment quenches the star formation of low-mass satellite galaxies from the inside-out

ABSTRACT Environment plays a critical role in the star formation history of galaxies. Tidal and hydrodynamical stripping, prominent in cluster environment, can remove the peripheral gas of galaxies and star formation may thus be environmentally suppressed from the outside-in. We revisit the environmental dependence of the radial gradient of specific star formation rate (sSFR) profile. We probe the radial gradient by using the archival spectral indices D4000n and HδA measured from SDSS fibre spectra, to indicate central sSFR and the total sSFR from fitting the spectral energy distribution. Despite the low spatial resolution, the wealth of SDSS data allows to disentangle the dependences on stellar mass, sSFR, and environment. We find that low-mass satellite galaxies in the mass range $9\lt \mathrm{log}\, \mathcal {M}_{\star }/\mathcal {M}_{\odot }\lt 9.8$ on average quench in more inside-out pattern compared to isolated galaxies matched in mass, sSFR, and fibre coverage. This environmental effect is particularly strong for galaxies below the star formation main sequence, and peaks for those in the core of massive clusters where the phase-space diagram reveals clear links between the inside-out quenching and orbital properties. Our results suggest that both tidal and hydrodynamical interactions in cluster environment suppress the star formation of satellites mainly from the inside-out. As accreted gas of low angular momentum from hot gas haloes is an important source for replenishing central gas reservoir, we discuss how gas stripping in clusters may lead to starvation and cause inside-out quenching when the outer star-forming discs are not significantly affected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Using DUNE to shed light on the electromagnetic properties of neutrinos

We study future DUNE sensitivity to various electromagnetic couplings of neutrinos, including magnetic moments, milli-charges, and charge radii. The DUNE PRISM capabilities play a crucial role in constraining the electron flavored couplings. We find that DUNE will be able to place the strongest beam based constraint on the muon-neutrino magnetic moment by improving on LSND’s bounds by roughly a factor of two, although Borexino’s constraint from solar neutrinos will be stronger. For the muon neutrino millicharge DUNE can place the leading beam based bound, with two orders of magnitude improvement compared to the existing COHERENT constraint, suggesting that DUNE can be useful for light mediators more generally. Despite this strength, the millicharge bounds are not competitive with strong bounds from stellar cooling, beta-decay, and matter stability. Finally, DUNE may be able to test the SM prediction for the muon neutrino charge radius, by placing a constraint two times better than CHARM-II and CCFR experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Counterpart identification and classification for eRASS1 and characterisation of the active galactic nuclei content

Context. Accurately accounting for the Active Galactic Nucleus (AGN) phase in galaxy evolution requires a large, clean AGN sample. This is now possible with SRG/eROSITA, which completed its first all-sky X-ray survey (eRASS1) on June 12, 2020. The public Data Release 1 (DR1, Jan 31, 2024) includes 930,203 sources from the western Galactic hemisphere. Aims. The data enable the selection of a large AGN sample and the discovery of rare sources. However, scientific return depends on accurate characterisation of the X-ray emitters, requiring high-quality multi-wavelength data. This paper presents the identification and classification of optical and infrared counterparts to eRASS1 sources. Methods. Counterparts to eRASS1 X-ray point sources were identified using Gaia DR3, CatWISE2020, and Legacy Survey DR10 (LS10) with the Bayesian NWAY algorithm and trained priors. Sources were classified as Galactic or extragalactic via a machine-learning model combining optical/IR and X-ray properties, trained on a reference sample. For extragalactic LS10 sources, photometric redshifts were computed using CIRCLEZ. Results. Within the LS10 footprint, all 656,614 eROSITA/DR1 sources have at least one possible optical counterpart; ∼570 000 are extragalactic and likely AGN. Half are new detections compared to AllWISE, Gaia, and Quaia AGN catalogues. Gaia and CatWISE2020 counterparts are less reliable, due to the survey’s shallowness and the limited amount of features available to assess the probability of being an X-ray emitter. In the Galactic plane, where the overdensity of stellar sources also increases the chance of associations, using conservative reliability cuts, we identified approximately 18 000 Gaia and 55 000 CatWISE2020 extragalactic sources. Conclusions. We have released three high-quality counterpart catalogues – plus the training and validation sets – as a benchmark for the field. These datasets have many applications, but in particular, they empower researchers to build AGN samples tailored for completeness and purity, accelerating the hunt for the Universe’s most energetic engines.

X-rays: general↗
DOE OSTI2024

Self-interacting Dark Matter Interpretation of Crater II

Abstract The satellite galaxy Crater II of the Milky Way is extremely cold and exceptionally diffuse. These unusual properties are challenging to understand in the standard model of cold dark matter. We use controlled N -body simulations to investigate the formation of Crater II in self-interacting dark matter (SIDM), where dark matter particles can scatter and thermalize. Taking the orbit motivated by the measurements from Gaia Early Data Release 3, we show a strong self-interacting cross section per particle mass of 60 cm 2 g −1 is favored for Crater II. The simulated SIDM halo, with a 1 kpc core, leads to both a low stellar velocity dispersion and a large half-light radius for Crater II. These characteristics remain robust regardless of the initial stellar distribution.

Astronomy & Astrophysics↗
DOE OSTI2023

How the dynamical properties of globular clusters impact their $γ$-ray and X-ray emission

The X-ray and $γ$-ray emission of globular clusters (GCs) is attributed to their large fraction of compact binary systems, especially those with millisecond pulsars (MSPs). We analyse a population of 124 Galactic GCs to investigate how their dynamical properties affect the formation and evolution of compact binary systems and how this can be translated into the clusters’ observed X-ray and γ-ray emission. We use mainly Chandra X-ray Observatory and Fermi Large Area Telescope observations to achieve our goals and start by detecting 39 GCs in γ rays, seven of which are not listed in previous Fermi -LAT catalogs. Additionally, we find that the total number of X-ray sources within a GC and its γ-ray luminosity are linearly correlated with the stellar encounter rate, indicating that compact binary systems are mainly formed via close stellar encounters. We also find an unexpected rise in the number of X-ray sources for GCs with low rates of stellar encounters, suggesting that there is a dynamical threshold where the formation of X-ray sources is dominated by stellar encounters. Furthermore, we use the Heggie-Hills law to find that subsequent stellar encounters in these compact binaries will, on average, make the binaries even harder, with basically no possibility of binary ionization. Finally, we find that all GCs are point-like sources in γ rays, indicating that the MSPs are concentrated in the clusters’ cores, likely due to dynamical friction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Constraining the Milky Way Mass Profile with Phase-space Distribution of Satellite Galaxies

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M={1.23}_{-0.18}^{+0.21}\times {10}^{12}{M}_{\odot }$ and a concentration of $c={9.4}_{-2.1}^{+2.8}$ with the prior based on the M–c relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of the kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. Finally, by varying the disruption rate of such satellites in this simulation, we estimate a ~5% scatter in the inferred MW halo mass among hydrodynamics-based simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Photometric Objects around Cosmic Webs (PAC) Delineated in a Spectroscopic Survey. I. Methods

We provide a method for estimating the projected density distribution ${\bar{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) 3 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 10 9.0 to 10 12.0 M ⊙ around massive galaxies at z s ≈ 0.6. Using the abundance-matching method, we model ${\bar{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.

79 ASTRONOMY AND ASTROPHYSICS↗