High-density Core Research Articles

Kernel Density Estimation for the Interpretation of Seismic Big Data in Tectonics Using QGIS: The Türkiye–Syria Earthquakes (2023)

Numerous studies have utilized remote sensing techniques to analyze seismic data in active areas. Point density techniques, widely used in remote sensing, examine the spatial distribution of point clouds related to specific variables. Applying these techniques to complex tectonic settings, such as the East Anatolian Fault Zone, helps identify major active fractures using both surface and deep information. This study employed kernel density estimation (KDE) to compare two distinct point-cloud populations from the seismic event along the Türkiye–Syria border on 6 February 2023, providing insights into the main active orientations supporting the Global Tectonics framework. This study considered two populations of seismic foci point clouds containing over 40,000 events, recorded by the Turkish Disaster and Emergency Management Authority (AFAD) and Kandilli Observatory and Earthquake Research Institute (KOERI). These populations were divided into two datasets: crude and relocated-filtered. Kernel density analysis demonstrated that both datasets yielded similar geological interpretations. The high-density cores of both datasets perfectly matched, exhibiting identical structures consistent with geological knowledge. Areas with a minimal concentration of earthquakes at depth were also identified, separating different crustal strength levels.

Spatial distribution characteristics and influencing factors of youth soccer specialty schools in China

Youth soccer specialty schools are the primary project to strengthen the foundation of soccer talents, improve the level of soccer development, and achieve the dream of Chinese soccer. Taking 32,775 soccer specialty schools in China as the research sample, ArcGIS spatial analysis technology is used to explore the spatial distribution characteristics of soccer specialty schools and their influencing factors. The analysis shows that the overall distribution of the country forms five quantity concentration area and eight high-density core areas, with a high degree of spatial concentration and significant differences between the East and West. The factors influencing soccer specialty schools include nine factors in four dimensions: historical development, geographical environment, socio-economic, and educational resources. It was found that geographic environment and socio-economy provide the foundation and background for developing educational resources. However, the direct role of educational resources is more significant in a specific educational program, such as soccer specialty schools, followed by socio-economy and geographic environment. This study comprehensively analyzes and interprets soccer specialty schools’ spatial distribution characteristics and influencing factors. These results enrich the research on the spatial layout of sports, providing necessary theoretical support and practical references for the layout and high-quality development of related industries in China and other countries or regions.

Full-wave simulations on helicon and parasitic excitation of slow waves near the edge plasma

Helicon waves are thought to be promising in various tokamaks, such as DIII-D, because they can penetrate reactor-grade high-density cores and drive the off-axis current with higher efficiency. In the frequency regime ∼ 476 MHz, both slow electrostatic and fast electromagnetic helicon waves can coexist in DIII-D. If the antenna parasitically excites the slow mode, these waves can propagate along the magnetic field line into the scrape-off layer (SOL). Although the importance of the misalignment of the Faraday screen and the electron density in the SOL on the excitation and propagation of slow modes is well known, the conditions for minimizing slow mode excitation have yet to be optimized. Using the Petra-M simulation code in the 2D domain, we analyze the effects of the misalignment of the antenna in the poloidal direction, the misalignment of the Faraday screen in the toroidal direction, and the density in front of the antenna on slow mode generation. Our results suggest that the misalignment of the Faraday screen is a critical factor in reducing the slow mode and that the misalignment angle should be below ∼ 5° to minimize the slow wave excitation. When the electron density is higher than 3.5×1018 m−3 in the SOL, the generation of the slow mode from the antenna is minimized and unaffected by the misalignment of the Faraday screen.

Formation of high areal density core using an efficient and robust implosion method for fast ignition

Abstract A new fuel compression method for a fast ignition scheme is discussed. To form a high areal density fuel plasma for the ignition condition, homogenous isentropic compression (HIC) with solid spherical target is effective. We improve a multi-step pulse shape method that uses progressive shockwaves and reflected shockwaves for the compression, where a precisely controlled step-pulse laser drives the shockwaves to compress the fuel and suppress entropy increase. Another advantage of this approach is the relatively smooth high dense fuel is distributed at maximum compression time, compared to our previous design based on Kidder’s HIC method. In addition, we insert a power dip as a preconditioning before the last pulse step to reduce the electron and ion temperature near critical density. As a result, an optimum implosion is designed using 245 kJ of implosion laser energy to meet the ignition condition.

From expansion to efficiency: Machine learning-based forecasting of Japan's building material stocks under demographic declines

Japan's unique demographic trajectory, marked by population decline and aging, coupled with continued urbanization, presents distinct challenges for aligning built environment capacity with resource efficiency. This study aims to investigate the historical evolution and project future scenarios of building material stock (MS) and their spatial distribution across Japan's three major metropolitan areas. Through a comprehensive material flow and stock analysis, the historical accumulation of building materials from 2009 to 2020 was quantified, revealing a dominance of concrete and an increasing overall stock. The contributions of various driving factors to changes in construction areas were explored, identifying population dynamics as the predominant influence. Leveraging shared socioeconomic pathway scenarios (SSPs), this research forecasted building floor area and MS until 2050 under five distinct SSPs. The results indicated an overall reduction across all scenarios, yet with a continued concentration in high-density urban cores. The substantial gap between the highest and lowest projected MS scenarios highlighted opportunities for material conservation and emission reductions through sustainable practices. Sustainable urban development in densely populated regions necessitates a balance between infrastructure provision and environmental conservation, while in sparsely populated areas, the focus shifts to the efficient management and utilization of vacant properties and materials to cope with the impacts of significant population declines. By offering insights into the building floor area and MS implications of Japan's demographic changes, this study underscores the necessity of sustainable urban planning and resource management strategies to navigate the challenges posed by demographic shifts, ultimately contributing to sustainable development and environmental conservation goals.

The effect of flow-derived mechanical cues on the growth and morphology of platelet aggregates under low, medium, and high shear rates

Platelet aggregation is a dynamic process that can obstruct blood flow, leading to cardiovascular diseases. While many studies have demonstrated clear connections between shear rate and platelet aggregation, the impact of flow-derived mechanical signals on this process is not fully understood. The objective of this work is to investigate the role of flow conditions on platelet aggregation dynamics, including effects on growth, shape, density composition, and their potential correlation with binding processes that are characterised by longer (e.g., via αIIbβ3 integrin) and shorter (e.g., via VWF) initial binding times. In vitro blood perfusion experiments were conducted at wall shear rates of 800, 1600 and 4000 s-1. Detailed analysis of two modalities of experimental images was performed to offer insights into the morphology of platelet aggregates. A consistent structural pattern was observed across all samples: a high-density core enveloped by a low-density outer shell. An image-based 3D computational blood flow model was subsequently employed to study the local flow conditions, including binding availability time and flow-derived mechanical signals via shear rate and rate of elongation. The results show substantial dependence of the aggregation dynamics on these flow parameters. We found that the different binding mechanisms that prefer different flow regimes do not have a monotonic cross-over in efficiency as the flow increases. There is a significant dip in the cumulative aggregation potential in-between the preferred regimes. The results suggest that treatments targeting the biomechanical pathways could benefit from creating conditions that exploit these low-efficiency zones of aggregation.

К оценке геодинамической устойчивости геосфер Земли

The Earth, as a system, is characterized by certain average state parameters (density, temperature). These parameters cannot but depend on similar states of internal geospheres (subsystems), which follow the unified law of increasing entropy of the cooling Earth from the center in the direction of the Earth's radius. As a model for assessing the dynamic stability of the Earth's internal geospheres relative to each other, the principle of the golden proportion (golden section) was used with the given parameters of density, temperature, power of the geospheres to the extreme ratios of their states at points (0.618; 1.618; 2.618), which determine the conditions for the stable dynamic equilibrium of the compared parameters of geospheres. And if these relations satisfied the known theoretical, model, geophysical data, then the conclusion was made that the geospheres in the compared parameter relations were in a state close to a stable dynamic equilibrium of the exchange of matter and energy between them. Otherwise, it was recognized that the compared geospheres of the Earth were in a state of either unstable dynamic equilibrium with an assessment of the level of this deviation from the golden numbers, or the known values of the density and temperature of the geospheres, as well as their structure, should be subject to clarification. Average model estimates of the parameters of the Earth's internal geospheres reach a stable dynamic equilibrium at a level of deviation from the golden numbers of up to 2.72 %. This allows estimating the model duration of the modern geological activity of the Earth for another 124 million years until the formation of a new state of the existing continents. The instability of the dynamic equilibrium of density and temperature of matter is established at the boundaries: continental and oceanic crust; oceanic crust and Golitsin layer; at the boundary between the outer part of the Earth's core and the lower mantle; inside the Earth's core at the level of identified transition zones in the four-layer core model. The presence of spheres with a high core density from 16.05 to 26.45 g/cm3 suggests the stratification of its central part into essentially gold and platinoid. The presence of instability boundaries in the dynamic state of geospheres is in good agreement with the concept of plate tectonics and with the depths of plume origin. It can be assumed that the relative equilibrium of the internal structure of the Earth will continue for about 124 million years until the beginning of the formation of a new Pangea Proxima, which, according to Cristofor Scotese, will arise in 200 million years.

Spatial Distribution Characteristics and Influencing Factors of Intangible Cultural Heritage in the Yunnan, Guangxi, and Guizhou Rocky Desertification Area

The Yunnan, Guangxi, and Guizhou rocky desertification area is characterized by a fragile ecological environment, along with the inhabitation of many ethnic minorities, making it the region with the most poverty-alleviation targets among the 14 concentrated poverty-stricken areas in China. However, this area has obvious cultural resources, but it lacks systematic study. Therefore, integral investigation and research which can promote the sustainable development with abundant culture resource should pay more attention to it. This study, with 1113 national and provincial intangible cultural heritage projects in this area as research objects, used ArcGIS to visualize their geographic data for research. Geodetector was used to quantify the influence of relevant factors and analyze the distribution characteristics and influencing factors of intangible cultural heritage in this region. By providing a correlative strategy, the aim is to promote the coordinated development of the regional cultural economy. The results are as follows: (1) Guizhou exhibits the highest distribution of intangible cultural heritage projects, followed by Guangxi and Yunnan, indicating an uneven number of intangible cultural heritage projects. (2) A spatial analysis shows the clustered spatial distribution of intangible cultural heritage projects in this area, forming a high-density core area and multiple sub high-density areas with high spatial autocorrelation. (3) The spatial distribution of intangible cultural heritage in the rocky desertification areas of Yunnan, Guangxi, and Guizhou is generally correlated with survival orientation and natural environmental factors and positively correlated with development orientation and social, economic, ethnic, and cultural factors. Among many factors, cultural environment quantity exhibits the strongest correlation with the spatial distribution of intangible cultural heritage. (4) Based on the degree of correlation with different influencing factors, four development strategies are proposed: digitization of intangible cultural heritage inheritance and development, multi-media intangible cultural heritage creation and dissemination, innovative construction of intangible cultural heritage exhibitions, and tourism-driven dynamic inheritance of intangible cultural heritage.

Evaluation of Satisfaction with Spatial Reuse of Industrial Heritage in High-Density Urban Areas: A Case Study of the Core Area of Beijing’s Central City

Industrial heritage is regarded as an important stock of spatial resources in cities, which highlights its utilization value for urban regeneration in high-density urban areas. With the dramatic increase in the number of industrial heritage reuse projects, how to scientifically evaluate the satisfaction with their spatial reuse is a key part of the solution for the mutual balance between heritage preservation and urban renewal. This paper takes eight industrial heritage conversion and utilization projects in the high-density core area of Beijing’s central city as examples; establishes an evaluation system for the satisfaction with the spatial reuse through six dimensions, namely, functional replacement, transportation accessibility, carrying capacity, public space, boundary form, and recognition of value; and uses the IPA method to evaluate the cases. This method is used to determine the degree of user satisfaction with the spatial reuse of industrial heritage in the core area of Beijing’s central city and to summarize the advantages and problems of its reuse. The results of this study reveal a trend toward the “community-oriented” re-generation of industrial heritage in the core area of Beijing’s central city, and this paper proposes recommendations for adaptive use to support high-quality urban regeneration work.

Seismic lithology and geological implications of pre-salt jurassic carbonate rocks in the Amu Darya Right Bank

The formation and evolution of oil and gas reservoirs are typically related to the environment where sedimentary facies occur and regional sedimentary facies are distributed. The accurate characterization of reef–shoal complexes is one of the important tasks for exploring oil and gas in carbonate reservoirs. However, most traditional methods for characterizing sedimentary facies are based on high-density drilling data, core data, and thin-section data, often making them less suitable in areas with few drilled wells or sparse well networks. Hence, they often cannot characterize sedimentary facies accurately, severely hindering the oil and gas exploration success rate. This study seeks to characterize sedimentary facies using limited amount of data from drilled wells. An innovative model representing the relationship between lithofacies and seismic attributes was established by comprehensively analyzing the relationship among lithologic associations, principal component analysis, and regression analysis based on seismic lithology and data from various sources. The sedimentary facies of the Upper Jurassic (Oxfordian) formation in the eastern region of the Amu Darya Right Bank were accurately characterized using Red-Green-Blue (RGB) fusion. These analyses indicate that the relationship among micrite, bioclastic micrite, and bioclastic limestone can effectively represent the developmental pattern of sedimentary facies in the study area. The planar distribution of sedimentary facies based on the relationship between the lithologic associations identified through seismic interpretation revealed the relationship between sedimentary facies and the paleomorphology of the region. Results show that the paleo-landform in the study area's western zone is low-lying, where shoals are poorly developed; the paleo-landform in the central zone is high, with high-energy water and well-developed and widely distributed thick shoals; and the eastern zone is far from the open sea and affected by the central-zone landform, with low-energy water and low-quality reef–shoal complexes. Controlled by the unique landforms of the study area, the western reservoirs are thin and have low quality, central reservoirs are thick with good petrophysical properties, and eastern reservoirs have differing thickness and petrophysical properties.

The principle of magnetic flux switch

This paper introduces an innovative Magnetic Switch (MFS) designed to control and alter magnetic flux within energy system components, offering an alternative to conventional power electronic devices. The MFS comprises a low-current control coil, a control core, and a high-density magnetic flux-carrying main core combined with a main coil energy system. In this novel magnetic configuration, when a low-power current excites the control coil, the magnetic flux in the main core (supplied by the main coil) decreases to nearly zero. Conversely, when the control coil disconnects from the power source, the magnetic flux within the main core attains its maximum value. This operation positions the MFS as a groundbreaking concept within magnetic-based energy systems, akin to transistors in power electronics. The main outcomes of the MFS concept are that it can vary the magnetic flux of the main core in the large range, and it is a fast magnetic switch with a simple and low power loss control circuit and an independent control coil from the main coil. Analytical studies thoroughly elucidate the performance and advantages of this proposed magnetic switch, substantiated by Finite Element Method simulations and experimental prototype outcomes.

Resolving Phase Transition Properties of Dense Matter through Tidal-excited g-mode from Inspiralling Neutron Stars

The investigation of the phase state of dense matter is hindered by complications of first-principle nonperturbative quantum chromodynamics. By performing the first consistent general-relativistic calculations of tidal-excited g-mode of neutron stars with a first-order strong interaction phase transition in the high-density core, we demonstrate that the gravitational-wave signal during binary neutron star inspiral probes their innermost hadron-quark transition and provides potent constraints from present and future gravitational-wave detectors.

Chapter 3: The Origins and Evolution of Planetary Systems.

The materials that form the diverse chemicals and structures on Earth-from mountains to oceans and biological organisms-all originated in a universe dominated by hydrogen and helium. Over billions of years, the composition and structure of the galaxies and stars evolved, and the elements of life, CHONPS, were formed through nucleosynthesis in stellar cores. Climactic events such as supernovae and stellar collisions produced heavier elements and spread them throughout the cosmos, often to be incorporated into new, more metal-rich stars. Stars typically form in molecular clouds containing small amounts of dust through the collapse of a high-density core. The surrounding nebular material is then pulled into a protoplanetary disk, from which planets, moons, asteroids, and comets eventually accrete. During the accretion of planetary systems, turbulent mixing can expose matter to a variety of different thermal and radiative environments. Chemical and physical changes in planetary system materials occur before and throughout the process of accretion, though many factors such as distance from the star, impact history, and level of heating experienced combine to ultimately determine the final geophysical characteristics. In Earth's planetary system, called the Solar System, after the orbits of the planets had settled into their current configuration, large impacts became rare, and the composition of and relative positions of objects became largely fixed. Further evolution of the respective chemical and physical environments of the planets-geosphere, hydrosphere, and atmosphere-then became dependent on their local geochemistry, their atmospheric interactions with solar radiation, and smaller asteroid impacts. On Earth, the presence of land, air, and water, along with an abundance of important geophysical and geochemical phenomena, led to a habitable planet where conditions were right for life to thrive.

High-density cluster core-based [formula omitted]-means clustering with an unknown number of clusters

The k-means algorithm, known for its simplicity and adaptability, faces challenges related to manual cluster number selection and sensitivity to initial centroid placement. This paper introduces an innovative framework aimed at overcoming these challenges. By proposing a data-driven cluster number estimation method and a robust initialization strategy based on high-density cluster cores, our approach revolutionizes k-means, unlocking its full unsupervised potential and ensuring superior performance, even in scenarios involving overlapping clusters. The method employs a novel density-based technique to accurately identify cluster cores, resulting in substantial improvements over existing methods. Rigorous experimentation on synthetic and real-world datasets demonstrates an average performance enhancement of 15% in terms of the Adjusted Rand Index for datasets with overlapping clusters, surpassing the capabilities of state-of-the-art density-based clustering methods and traditional k-means. Moreover, our method autonomously determines the optimal number of clusters, facilitating true unsupervised learning and eliminating the impact of initial centroid placement on clustering outcomes. This leads to stable and consistent results, addressing key limitations of the conventional k-means algorithm. The practical applicability of our approach is exemplified in image segmentation tasks, showcasing its versatility and reliability in real-world scenarios.

A New Kilohertz Gravitational-wave Feature from Rapidly Rotating Core-collapse Supernovae

We present self-consistent three-dimensional core-collapse supernova simulations of a rotating 20M ⊙ progenitor model with various initial angular velocities from 0.0 to 4.0 rad s−1 using the smoothed particle hydrodynamics code SPHYNX and the grid-based hydrodynamics code FLASH. We identify two strong gravitational-wave features with peak frequencies of ∼300 Hz and ∼1.3 kHz in the first 100 ms postbounce. We demonstrate that these two features are associated with the m = 1 deformation from the proto-neutron star (PNS) modulation induced by the low-T/∣W∣ instability, regardless of the simulation code. The 300 Hz feature is present in models with an initial angular velocity between 1.0 and 4.0 rad s−1, while the 1.3 kHz feature is only present in a narrower range, from 1.5 to 3.5 rad s−1. We show that the 1.3 kHz signal originates from the high-density inner core of the PNS, and the m = 1 deformation triggers a strong asymmetric distribution of electron antineutrinos. In addition to the 300 Hz and 1.3 kHz features, we also observe one weaker but noticeable gravitational-wave feature from higher-order modes in the range between 1.5 and 3.5 rad s−1. Its initial peak frequency is around 800 Hz, and it gradually increases to 900–1000 Hz. Therefore, in addition to the gravitational bounce signal, the detection of the 300 Hz, 1.3 kHz, the higher-order mode, and even the related asymmetric emission of neutrinos could provide additional diagnostics for estimating the initial angular velocity of a collapsing core.

Development of High-Density and Low-Loss Powder Magnetic Core for Reactor in Hybrid Vehicles

Development of High-Density and Low-Loss Powder Magnetic Core for Reactor in Hybrid Vehicles

Illuminating evaporating protostellar outflows: ERIS/SPIFFIER reveals the dissociation and ionization of HH 900

ABSTRACT Protostellar jets and outflows are signposts of active star formation. In H ii regions, molecular tracers like CO only reveal embedded portions of the outflow. Outside the natal cloud, outflows are dissociated, ionized, and eventually completely ablated, leaving behind only the high-density jet core. Before this process is complete, there should be a phase where the outflow is partially molecular and partially ionized. In this paper, we capture the HH 900 outflow while this process is in action. New observations from the Enhanced Resolution Imager and Spectrograph/SPIFFIER near-infrared (IR) integral field unit spectrograph using the K-middle filter (λ = 2.06–2.34 μm) reveal H2 emission from the dissociating outflow and Br-γ tracing its ionized skin. Both lines trace the wide-angle outflow morphology but H2 only extends ∼5000 au into the H ii region while Br-γ extends the full length of the outflow (∼12 650 au), indicating rapid dissociation of the molecules. H2 has higher velocities further from the driving source, consistent with a jet-driven outflow. Diagnostic line ratios indicate that photoexcitation, not just shocks, contributes to the excitation in the outflow. We argue that HH 900 is the first clear example of an evaporating molecular outflow and predict that a large column of neutral material that may be detectable with Atacama Large Millimeter Array accompanies the dissociating molecules. Results from this study will help guide the interpretation of near-IR images of externally irradiated jets and outflows such as those obtained with the JWST in high-mass star-forming regions where these conditions may be common.

Resonance in radio frequency sheath admittance and enhanced impurity emission near the ion cyclotron frequency

Ion cyclotron resonance heating (ICRH) is of considerable interest among all auxiliary heating techniques, because it transfers power directly to ions, targets the high-density core, and involves the cheapest radio-frequency (RF) components. During ICRH operation, RF sheaths form on the ICRH antenna itself, nearby hardware, and far-field surfaces. These sheaths are associated with large hot-spot formation and impurity emissions. This work presents high-resolution numerical modelling of RF sheaths in nuclear fusion scenarios using hPIC2, a Debye-scale particle-in-cell code. The modelling reveals a new RF sheath phenomenon which occurs when the RF is resonant with a species ion cyclotron frequency or harmonic in an oblique magnetic field. This resonance of the RF sheath causes modifications of the energy-angle distributions of the ions impacting on the walls, with consequent increase in wall impurity emission. Due to the 1/R scaling of the magnetic field in a tokamak, such cases are possible on divertor or vacuum chamber surfaces, depending on the geometry of the tokamak and the RF heating scenario. A simple physics interpretation using a driven damped harmonic oscillator is proposed, where the ion sheath transit time plays the role of damping the RF sheath admittance. Critically, the resonance leads to increased ion heat flux at the wall as well as increased physical sputtering, despite a lack of increase in RF rectified sheath potential.

Influence of Aspect Ratio on the Flexural and Buckling Behavior of an Aluminium Sandwich Composite: A Numerical and Experimental Approach.

In the field of engineering materials, lightweight and ultra-lightweight composites are used in real time to a greater extent, with high-performance targeting for tailor-made systems in aerospace, automotive, and biomedical applications. Sandwich composites are among the most popular lightweight materials used in structural and vehicle-building applications. In the present investigation, one such sandwich composite laminate composed of aluminum face sheets and a high-density polyethylene core was considered to analyze sandwich composites' flexural and buckling behavior experimentally and numerically. The influence of aspect ratios, such as length to thickness and width to thickness, on the flexural and buckling performance of sandwich composite laminates was explored in the study. Laminates with different widths, namely, 10, 12, and 15 mm, and a uniform thickness and length of 3 mm and 150 mm, respectively, were used for flexural analysis, whereas laminates with widths of 10, 12, and 15 mm and a uniform thickness and length of 3 mm and 350 mm, respectively, were used for buckling analysis. The geometrical influence of the laminates on mechanical performance was studied through performance measures such as critical bending load, flexural stiffness, inter-laminar shear stress, and critical buckling load. A significant influence of aspect ratio on the mechanical behavior of the laminates was observed using both experimental and numerical approaches. Flexural behavior was observed to be better at greater widths, namely, 15 mm, and with a minimum support span of 90 mm due to reduced spring back effects and increased bending resistance. A maximum width of 15 mm allowed for a higher buckling load capacity similar to that of bending resistance. A critical buckling load of 655.8 N seemed to be the maximum and was obtained for the highest aspect ratio, b/t = 5. The soft core and ductile metal face sheets offered combined resistance to both bending and buckling. A lower aspect ratio (span to thickness) rendered these sandwich laminates better in terms of both bending and buckling.

Dosimetric Validation for Prospective Clinical Trial of GRID Collimator-Based Spatially Fractionated Radiation Therapy: Dose Metrics Consistency and Heterogeneous Pattern Reproducibility

Dose heterogeneity within a tumor target is likely responsible for the biologic effects and local tumor control from spatially fractionated radiation therapy (SFRT). This study used a commercially available GRID-pattern dose mudulated nonuniform radiation therapy (GRID) collimator to assess the interplan variability of heterogeneity dose metrics in patients with various bulky tumor sizes and depths. The 3-dimensional heterogeneity metrics of 14 bulky tumors, ranging from 155 to 2161 cm3 in volume, 6 to 23 cm in equivalent diameter, and 3 to 13 cm in depth, and treated with GRID collimator-based SFRT were studied. A prescription dose of 15 Gy was given at the tumor center with 6 MV photons. The dose-volume histogram indices, dose heterogeneity parameters, and peak/valley dose ratios were derived; the equivalent uniform doses of cancer cells with various radiosensitivities in each plan were estimated. To account for the spatial fractionation, high dose core number density of the tumor target was defined and calculated. Among 14 plans, the dose-volume histogram indices D5, D10, D50, D90, and D95 (doses covering 5%, 10%, 50%, 90%, and 95% of the target volume) were found within 10% variation. The dose ratio of D10/D90 also showed a moderate consistency (range, 3.9-5.0; mean, 4.4). The equivalent uniform doses were consistent, ranging from 4.3 to 5.5 Gy, mean 4.6 Gy, for radiosensitive cancer cells and from 5.8 to 6.9 Gy, mean 6.2 Gy, for radioresistant cancer cells. The high dose core number density was within 20% among all plans. GRID collimator-based SFRT delivers a consistent heterogeneity dose distribution and high dose core density across bulky tumor plans. The interplan reproducibility and simplicity of GRID therapy may be useful for certain clinical indications and interinstitutional clinical trial design, and its heterogeneity metrics may help guide multileaf-collimator-based SFRT planning to achieve similar or further optimized dose distributions.