Computational study of tetrahedral Al–Si and octahedral Al–Mg ordering in phengite

 As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K₂ ^[6](Al₃Mg)^[4](Si₇Al)O₂₀(OH)₄. Values of the atomic exchange interaction parameters J _n used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range order was observed. Received: 8 August 2002 / Accepted: 14 February 2003 Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility.     

A TWO-LAYER PRIMITIVE EQUATION OCEAN CURRENT MODEL AND ITS SIMULATION

A two-layer primitive equation model is developed in this paper. The capabilities of this model aretested by the use of multiyearly averaged January and July sea surface level pressure fields and windfields which can be diagnosed from the pressure fields. The results show that the ocean surface currentsand undercurrents in the second layer driven by the sea surface wind and the sea surface pressure areclose to the observation. The results are also compared with that of the IAP OGCM and the OSUOGCM.     

巨分子云聚合形成N体模拟中粒子数的效应

通过N体数值模拟以不同粒子数对星系中巨分子云的形成的影响进行了研究。结果表明：在聚合形成机制下，当对不同数目的分子云取相同的平均密度时，基本分子云的有效半径与其相应质量的立方根成正比，巨分子云的碎裂率与模拟基本分子云的数目无关。     

The evolution of binary fractions in globular clusters

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code ( startrack ) and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing three- and four-body interactions ( fewbody ). We find that the combination of stellar evolution and dynamical interactions (binary–single and binary–binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster such as 47 Tuc, assuming an initial binary fraction of 100 per cent, is only ∼ 5–10 per cent. We show that this is in good agreement with recent Hubble Space Telescope observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modelling of globular clusters, suggesting that 'realistic models' should incorporate much larger initial binary fractions than has usually been the case in the past.     

Elucidating Dominant Factors Affecting Land Surface Hydrological Simulations of the Community Land Model over China

In order to compare the impacts of the choice of land surface model(LSM) parameterization schemes, meteorological forcing, and land surface parameters on land surface hydrological simulations, and explore to what extent the quality can be improved, a series of experiments with different LSMs, forcing datasets, and parameter datasets concerning soil texture and land cover were conducted. Six simulations are run for the Chinese mainland on 0.1° × 0.1° grids from 1979 to 2008, and the simulated mon...     

Pyroclastic Flow Hazard at Volcán Citlaltépetl

Volcán Citlaltépetl (Pico de Orizaba) with an elevation of 5,675 m is the highest volcano in North America. Its most recent catastrophic events involved the production of pyroclastic flows that erupted approximately 4,000, 8,500, and 13,000 years ago. The distribution of mapped deposits from these eruptions gives an approximate guide to the extent of products from potential future eruptions. Because the topography of this volcano is constantly changing computer simulations were made on the present topography using three computer algorithms: energy cone, FLOW2D, and FLOW3D. The Heim Coefficient (), used as a code parameter for frictional sliding in all our algorithms, is the ratio of the assumed drop in elevation (H) divided by the lateral extent of the mapped deposits (L). The viscosity parameter for the FLOW2D and FLOW3D codes was adjusted so that the paths of the flows mimicked those inferred from the mapped deposits. We modeled two categories of pyroclastic flows modeled for the level I and level II events. Level I pyroclastic flows correspond to small but more frequent block-and-ash flows that remain on the main cone. Level II flows correspond to more widespread flows from catastrophic eruptions with an approximate 4,000-year repose period. We developed hazard maps from simulations based on a National Imagery and Mapping Agency (NIMA) DTED-1 DEM with a 90 m grid and a vertical accuracy of ±30 m. Because realistic visualization is an important aid to understanding the risks related to volcanic hazards we present the DEM as modeled by FLOW3D. The model shows that the pyroclastic flows extend for much greater distances to the east of the volcano summit where the topographic relief is nearly 4,300 m. This study was used to plot hazard zones for pyroclastic flows in the official hazard map that was published recently.     

蒙特卡洛模型在现代黄河三角洲浅层沉积物压实速率模拟中的应用

在进行黄河三角洲浅层沉积物固结压实研究中,过去传统的方法通常基于钻孔和剖面,进行局部计算分析。为克服钻孔分布不均的局限性需研究整个现代黄河三角洲地区浅层沉积物的压实沉降特征,基于1987年黄河三角洲地区综合工程地质勘察获取的钻孔资料,采用了数学地质中的蒙特卡洛方法,结合土力学的分层总和法对整个现代黄河三角洲浅层沉积物的压实速率进行模拟,并分析了压实模拟的具体影响因素。结果表明当前浅层沉积物的90%累积概率(P90)的压实速率变化范围为0.18~9.07 mm/a,影响模拟沉积层的压实速率主要因素为沉积物的初始孔隙比、压缩系数和平均沉积速率,上述要素均与压实速率呈显著正相关,软土层应是浅沉积地层压实沉降的主要贡献层。提出按照时间序列分别对不同流路时期叶瓣进行模拟,更符合现代黄河三角洲地区的沉积环境和历史,并对比23个地面沉降水准监测数据,定量分离了浅层沉积物固结压实分量,其约占沉降总量的13%。     

HD 82943行星系统的2:1共振与稳定性

用直接数值积分方法通过模拟不同的行星构型探讨了HD82943行星系统(由两颗共振的巨行星组成)的长期动力演化,同时,还研究了在相空间的轨道运动.在对系统长达107年数值积分中,发现所有的稳定轨道均与2:1共振相关联.典型地,在相同的时标内,两个共振幅角θ1和θ2同时(或其中之一)存在秤动.由于共振幅角在一定范围内的秤动,因而使两颗行星轨道半长径被约束而表现为规则运动模式.另外,利用分析模型(包含了外行星偏心率e2的因素),还讨论了对于不同取值的e2和相对近星点经度θ时,内行星在相空间的运动,并发现2:1轨道共振对于相对较小的e2以及当θ=0°时易于保持.此外,适中的e2将导致系统的两颗行星进入深度共振状态.再者,分析模型和数值计算的结果吻合得很好,两者都揭示了行星系统的2:1共振结构.     

Numerical simulations of the Hyades dynamics and the nature of the moving Hyades cluster

We present our numerical simulations of the dynamical evolution of the Hyades open cluster. The simulations were performed usinga modified NBODY6 algorithm that included tidal forces and a realistic orbit of the cluster in a gravitational field described by the Miyamoto-Nagai potential. Our goal was to study the nature of movingclu sters. We show that the stars that were earlier cluster members could be later identified within a sphere of 50 pc in diameter around the Sun. The number of such stars for the chosen initial mass and virial radius of the cluster does not exceed ten. The maximum space velocity of these stars relative to the core of the current cluster does not exceed 3 km s^?1. Our numerical simulations confirm the assumption that some of the moving clusters near the Sun could consist of stars that have escaped from open clusters in the course of their dynamical evolution.     

Simulating the formation of a protocluster at z∼ 2

We present results from two high-resolution hydrodynamical simulations of protocluster regions at   z ≃ 2.1  . The simulations have been compared to observational results for the so-called Spiderweb galaxy system, the core of a putative protocluster region at   z = 2.16  , found around a radio galaxy. The simulated regions have been chosen so as to form a poor cluster with   M ₂₀₀≃ 10¹⁴  h ⁻¹ M_⊙  (C1) and a rich cluster with   M ₂₀₀≃ 2 × 10¹⁵  h ⁻¹ M_⊙  (C2) at   z = 0  . The simulated protoclusters show evidence of ongoing assembly of a dominating central galaxy. The stellar mass of the brightest cluster galaxy of the C2 system is in excess with respect to observational estimates for the Spiderweb galaxy, with a total star formation rate which is also larger than indicated by observations. We find that the projected velocities of galaxies in the C2 cluster are consistent with observations, while those measured for the poorer cluster C1 are too low compared with the observed velocities. We argue that the Spiderweb complex resembles the high-redshift progenitor of a rich galaxy cluster. Our results indicate that the included supernovae feedback is not enough to suppress star formation in these systems, supporting the need of introducing active galactic nuclei feedback. According to our simulations, a diffuse atmosphere of hot gas in hydrostatic equilibrium should already be present at this redshift, and enriched at a level comparable to that of nearby galaxy clusters. The presence of this gas should be detectable with future deep X-ray observations.