Analytical approach to subhalo population in dark matter haloes

In the standard model of cosmic structure formation, dark matter haloes form by gravitational instability. The process is hierarchical: smaller systems collapse earlier, and later merge to form larger haloes. The galaxy clusters, hosted by the largest dark matter haloes, are at the top of this hierarchy and representing the largest as well as the last structures formed in the Universe, while the smaller and first haloes are those Earth-sized dark subhaloes that have been both predicted by theoretical considerations and found in numerical simulations, though there do not exist any observational hints of their existence. The probability that a halo of mass m at redshift z will be part of a larger halo of mass M at the present time can be described in the frame of the extended Press & Schecter theory making use of the progenitor (conditional) mass function. Using the progenitor mass function, we calculate analytically, at redshift zero, the distribution of subhaloes in mass, formation epoch and rarity of the peak of the density field at the formation epoch. That is done for a Milky Way size system, assuming both a spherical and an ellipsoidal collapse model. Our calculation assumes that small progenitors do not lose mass due to dynamical processes after entering the parent halo, and that they do not interact with other subhaloes. For a Λ cold dark matter power spectrum, we obtain a subhalo mass function  d n /d m   proportional to   m ^−α  with a model-independent  α∼ 2  . Assuming that the dark matter is a weakly interacting massive particle, the inferred distributions are used to test the feasibility of an indirect detection in the γ-ray energy band of such a population of subhaloes with a Gamma-ray Large Area Space Telescope like satellite.     

A merger tree with microsolar mass resolution: application to γ-ray emission from subhalo population

Observations of the southern Cepheid ℓ Car to yield the mean angular diameter and angular pulsation amplitude have been made with the Sydney University Stellar Interferometer at a wavelength of 696 nm. The resulting mean limb-darkened angular diameter is 2.990 ± 0.017 mas (i.e. ± 0.6 per cent) with a maximum-to-minimum amplitude of 0.560 ± 0.018 mas corresponding to 18.7 ± 0.6 per cent in the mean stellar diameter. Careful attention has been paid to uncertainties, including those in measurements, in the adopted calibrator angular diameters, in the projected values of visibility squared at zero baseline, and to systematic effects. No evidence was found for a circumstellar envelope at 696 nm. The interferometric results have been combined with radial displacements of the stellar atmosphere derived from selected radial velocity data taken from the literature to determine the distance and mean diameter of ℓ Car. The distance is determined to be 525 ± 26 pc and the mean radius  169 ± 8 R_⊙  . Comparison with published values for the distance and mean radius shows excellent agreement, particularly when a common scaling factor from observed radial velocity to pulsation velocity of the stellar atmosphere (the p -factor) is used.     

High-Performance Parallel Solver for Integral Equations of Electromagnetics Based on Galerkin Method

A new parallel solver for the volumetric integral equations (IE) of electrodynamics is presented. The solver is based on the Galerkin method, which ensures convergent numerical solution. The main features include: (i) memory usage eight times lower compared with analogous IE-based algorithms, without additional restrictions on the background media; (ii) accurate and stable method to compute matrix coefficients corresponding to the IE; and (iii) high degree of parallelism. The solver’s computational efficiency is demonstrated on a problem of magnetotelluric sounding of media with large conductivity contrast, revealing good agreement with results obtained using the second-order finite-element method. Due to the effective approach to parallelization and distributed data storage, the program exhibits perfect scalability on different hardware platforms.     

Energy partitioning during seismic slip in pseudotachylyte-bearing faults (Gole Larghe Fault, Adamello, Italy)

The determination of the earthquake energy budget remains a challenging issue for Earth scientists, as understanding the partitioning of energy is a key towards the understanding the physics of earthquakes. Here we estimate the partition of the mechanical work density into heat and surface energy (energy required to create new fracture surface) during seismic slip on a location along a fault. Earthquake energy partitioning is determined from field and microstructural analyses of a fault segment decorated by pseudotachylyte (solidified friction-induced melt produced during seismic slip) exhumed from a depth of ~ 10 km—typical for earthquake hypocenters in the continental crust. Frictional heat per unit fault area estimated from the thickness of pseudotachylytes is ~ 27 MJ m^− 2. Surface energy, estimated from microcrack density inside clast (i.e., cracked grains) entrapped in the pseudotachylyte and in the fault wall rock, ranges between 0.10 and 0.85 MJ m^− 2. Our estimates for the studied fault segment suggest that ~ 97–99% of the energy was dissipated as heat during seismic slip. We conclude that at 10 km depth, less than 3% of the total mechanical work density is adsorbed as surface energy on the fault plane during earthquake rupture.     

基于冻融次数-物理时间比拟理论的冻土长期强度预测方法

冻融作用可以改变土的结构,引起其力学特征的变化,从而影响到冻土工程的稳定性。由于冻融周期的设定不尽相同,因此大量的试验结果不能有效地进行对比分析。另外,对于冻融作用下土体力学特征的预测也成为研究的难点。在冻土遗传蠕变理论的基础上,利用球型模板压入仪作为测试方法,提出了一种冻融次数-物理时间比拟的方法。把冻融周期(次数)转化成物理时间(min),将不同冻融周期土的长期强度曲线簇映射到同一坐标中,利用该归一化曲线进行长期抗剪强度的预测,并选取了两种土样进行了测试,分别得到了相关的长期抗剪强度的预测方程。该方法对冻融作用下土体力学行为的成果对比研究以及长期强度的预测有重要的理论意义,对寒区工程稳定性预测分析也具有工程实践价值。