-
Distribution of solids in the rings of the HD 163296 disk: a multiwavelength study
Authors:
G. Guidi,
A. Isella,
L. Testi,
C. J. Chandler,
H. B. Liu,
H. M. Schmid,
G. Rosotti,
C. Meng,
J. Jennings,
J. P. Williams,
J. M. Carpenter,
I. de Gregorio-Monsalvo,
H. Li,
S. F. Liu,
S. Ortolani,
S. P. Quanz,
L. Ricci,
M. Tazzari
Abstract:
In this paper we analyze new observations from ALMA and VLA, at a high angular resolution corresponding to 5 - 8 au, of the protoplanetary disk around HD 163296 to determine the dust spatial distribution and grain properties. We fit the spectral energy distribution as a function of the radius at five wavelengths from 0.9 to 9\,mm, using a simple power law and a physical model based on an analytic…
▽ More
In this paper we analyze new observations from ALMA and VLA, at a high angular resolution corresponding to 5 - 8 au, of the protoplanetary disk around HD 163296 to determine the dust spatial distribution and grain properties. We fit the spectral energy distribution as a function of the radius at five wavelengths from 0.9 to 9\,mm, using a simple power law and a physical model based on an analytic description of radiative transfer that includes isothermal scattering. We considered eight dust populations and compared the models' performance using Bayesian evidence. Our analysis shows that the moderately high optical depth ($τ$>1) at $λ\leq$ 1.3 mm in the dust rings artificially lower the millimeter spectral index, which should therefore not be considered as a reliable direct proxy of the dust properties and especially the grain size. We find that the outer disk is composed of small grains on the order of 200 $μ$m with no significant difference between rings at 66 and 100 au and the adjacent gaps, while in the innermost 30 au, larger grains ($\geq$mm) could be present. We show that the assumptions on the dust composition have a strong impact on the derived surface densities and grain size. In particular, increasing the porosity of the grains to 80\% results in a total dust mass about five times higher with respect to grains with 25\% porosity. Finally, we find that the derived opacities as a function of frequency deviate from a simple power law and that grains with a lower porosity seem to better reproduce the observations of HD163296. While we do not find evidence of differential trapping in the rings of HD163296, our overall results are consistent with the postulated presence of giant planets affecting the dust temperature structure and surface density, and possibly originating a second-generation dust population of small grains.
△ Less
Submitted 4 July, 2022;
originally announced July 2022.
-
The Ultramarine Simulation: properties of dark matter haloes before redshift 5.5
Authors:
Qiao Wang,
Liang Gao,
Cheng Meng
Abstract:
We introduce the Ultramarine simulation, an extremely large $N$-body simulation of the structure formation and evolution to redshift 5.5 at which cosmic reionization was just completed. The simulation evolves 2.1 trillion particles within a $512~h^{-1}$Mpc cube and has an unprecedented mass and force resolution for large volume simulations of this kind, $5.6\times 10^6 h^{-1}$M$_\odot$ and…
▽ More
We introduce the Ultramarine simulation, an extremely large $N$-body simulation of the structure formation and evolution to redshift 5.5 at which cosmic reionization was just completed. The simulation evolves 2.1 trillion particles within a $512~h^{-1}$Mpc cube and has an unprecedented mass and force resolution for large volume simulations of this kind, $5.6\times 10^6 h^{-1}$M$_\odot$ and $1.2~h^{-1}$kpc, respectively. We present some basic statistical results of the simulation, including the halo mass function, halo bias parameter as well as halo mass-concentration relation at high redshifts, and compare them with some existing representative models. We find excellent agreement with some models on the high redshift halo mass functions, but neither the halo bias factor nor halo mass-concentration relation. All halo bias models for comparison over-predicate high redshift halo bias by large factors, an accurate fit to our simulation is given. High redshift dark matter haloes still can be reasonably described with NFW model, the halo mass-concentration relations are monotonic, with more massive haloes having lower concentration, in disfavor of the upturn feature reported by some studies. The mass concentration relation has little evolution between $z=5.5$ to $z=10$, in contrast to strong evolution predicted by most existing models. In addition, concentration parameters of high redshift dark matter haloes are much lower than most model predictions.
△ Less
Submitted 23 October, 2022; v1 submitted 13 June, 2022;
originally announced June 2022.
-
PhotoNs-GPU:A GPU accelerated cosmological simulation code
Authors:
Qiao Wang,
Chen Meng
Abstract:
We present a GPU-accelerated cosmological simulation code, PhotoNs-GPU, based on algorithm of Particle Mesh Fast Multipole Method (PM-FMM), and focus on the GPU utilization and optimization. A proper interpolated method for truncated gravity is introduced to speed up the special functions in kernels. We verify the GPU code in mixed precision and different levels of interpolated method on GPU. A ru…
▽ More
We present a GPU-accelerated cosmological simulation code, PhotoNs-GPU, based on algorithm of Particle Mesh Fast Multipole Method (PM-FMM), and focus on the GPU utilization and optimization. A proper interpolated method for truncated gravity is introduced to speed up the special functions in kernels. We verify the GPU code in mixed precision and different levels of interpolated method on GPU. A run with single precision is roughly two times faster that double precision for current practical cosmological simulations. But it could induce a unbiased small noise in power spectrum. Comparing with the CPU version of PhotoNs and Gadget-2, the efficiency of new code is significantly improved. Activated all the optimizations on the memory access, kernel functions and concurrency management, the peak performance of our test runs achieves 48% of the theoretical speed and the average performance approaches to 35% on GPU.
△ Less
Submitted 28 July, 2021;
originally announced July 2021.
-
PHoToNs--A Parallel Heterogeneous & Threads oriented code for cosmological N-body simulation
Authors:
Qiao Wang,
Zongyan Cao,
Liang Gao,
Xuebin Chi,
Chen Meng,
Jie Wang,
Long Wang
Abstract:
We introduce a new code for cosmological simulations, PHoToNs, which has features on performing massive cosmological simulations on heterogeneous high performance Computer (HPC) and threads oriented programming. PHoToNs adopts a hybrid scheme to compute gravity force, with the conventional PM to compute the long-range force, the Tree algorithm to compute the short range force, and the direct summa…
▽ More
We introduce a new code for cosmological simulations, PHoToNs, which has features on performing massive cosmological simulations on heterogeneous high performance Computer (HPC) and threads oriented programming. PHoToNs adopts a hybrid scheme to compute gravity force, with the conventional PM to compute the long-range force, the Tree algorithm to compute the short range force, and the direct summation PP to compute the gravity from very close particles. A self-similar space filling Peano-Hilbert curve is used to decompose computing domain. Threads programming is highly used to more flexibly manage the domain communication, PM calculation and synchronization, as well as Dual Tree Traversal on the CPU+MIC platform. The scalability of the PHoToNs performs well and the efficiency of PP kernel achieves 68.6% of peak performance on MIC and 74.4% on CPU platforms. We also test the accuracy of the code against the much used Gadget-2 in the community and found excellent agreement.
△ Less
Submitted 27 February, 2018;
originally announced February 2018.