RASPA: molecular simulation software for adsorption and diffusion in flexible nanoporous materials

A new software package, RASPA, for simulating adsorption and diffusion of molecules in flexible nanoporous materials is presented. The code implements the latest state-of-the-art algorithms for molecular dynamics and Monte Carlo (MC) in various ensembles including symplectic/measure-preserving integrators, Ewald summation, configurational-bias MC, continuous fractional component MC, reactive MC and Baker's minimisation. We show example applications of RASPA in computing coexistence properties, adsorption isotherms for single and multiple components, self- and collective diffusivities, reaction systems and visualisation. The software is released under the GNU General Public License.     

Force and torque modelling of drilling simulation for orthopaedic surgery

The advent of haptic simulation systems for orthopaedic surgery procedures has provided surgeons with an excellent tool for training and preoperative planning purposes. This is especially true for procedures involving the drilling of bone, which require a great amount of adroitness and experience due to difficulties arising from vibration and drill bit breakage. One of the potential difficulties with the drilling of bone is the lack of consistent material evacuation from the drill's flutes as the material tends to clog. This clogging leads to significant increases in force and torque experienced by the surgeon. Clogging was observed for feed rates greater than 0.5 mm/s and spindle speeds less than 2500 rpm. The drilling simulation systems that have been created to date do not address the issue of drill flute clogging. This paper presents force and torque prediction models that account for this phenomenon. The two coefficients of friction required by these models were determined via a set of calibration experiments. The accuracy of both models was evaluated by an additional set of validation experiments resulting in average R² regression correlation values of 0.9546 and 0.9209 for the force and torque prediction models, respectively. The resulting models can be adopted by haptic simulation systems to provide a more realistic tactile output.     

Predicting the effects of climate change on population connectivity and genetic diversity of an imperiled freshwater mussel,Cumberlandia monodonta (Bivalvia: Margaritiferidae), in riverine systems

In the face of global climate change, organisms may respond to temperature increases by shifting their ranges poleward or to higher altitudes. However, the direction of range shifts in riverine systems is less clear. Because rivers are dendritic networks, there is only one dispersal route from any given location to another. Thus, range shifts are only possible if branches are connected by suitable habitat, and stream‐dwelling organisms can disperse through these branches. We used Cumberlandia monodonta (Bivalvia: Unionoida: Margaritiferidae) as a model species to investigate the effects of climate change on population connectivity because a majority of contemporary populations are panmictic. We combined ecological niche models (ENMs) with population genetic simulations to investigate the effects of climate change on population connectivity and genetic diversity of C. monodonta. The ENMs were constructed using bioclimatic and landscape data to project shifts in suitable habitat under future climate scenarios. We then used forward‐time simulations to project potential changes in genetic diversity and population connectivity based on these range shifts. ENM results under current conditions indicated long stretches of highly suitable habitat in rivers where C. monodonta persists; populations in the upper Mississippi River remain connected by suitable habitat that does not impede gene flow. Future climate scenarios projected northward and headwater‐ward range contraction and drastic declines in habitat suitability for most extant populations throughout the Mississippi River Basin. Simulations indicated that climate change would greatly reduce genetic diversity and connectivity across populations. Results suggest that a single, large population of C. monodonta will become further fragmented into smaller populations, each of which will be isolated and begin to differentiate genetically. Because C. monodonta is a widely distributed species and purely aquatic, our results suggest that persistence and connectivity of stream‐dwelling organisms will be significantly altered in response to future climate change.     

动态适应性生态经济区划模型及其应用

提高城市生态经济区划的精确性和动态性,对科学指导城市化发展具有重要的理论意义和应用价值。利用夜间灯光数据与人口密度建立线性模型,探索了以往用行政区域为最小统计单元数据的模拟细化问题;然后通过引入可变参数构建了动态适应性生态经济区划模型,在增强模型动态适应性的同时,将一级区划结果统一划分为生态管控区域、生态优先区域、优化开发区域和重点开发区域4个区域。以广州市增城区为典型案例,通过改进的动态适应性生态经济区划模型,运用GIS将增城区在两种情景下进行了模拟和对比,并提出了政策建议。区划结果符合当地发展特征,也为其他城市与区域的生态经济区划研究提供了科学方法。     

Wind and fairness in ski jumping: A computer modelling analysis

Wind is closely associated with the discussion of fairness in ski jumping. To counter-act its influence on the jump length, the International Ski Federation (FIS) has introduced a wind compensation approach. We applied three differently accurate computer models of the flight phase with wind (M1, M2, and M3) to study the jump length effects of various wind scenarios. The previously used model M1 is accurate for wind blowing in direction of the flight path, but inaccuracies are to be expected for wind directions deviating from the tangent to the flight path. M2 considers the change of airflow direction, but it does not consider the associated change in the angle of attack of the skis which additionally modifies drag and lift area time functions. M3 predicts the length effect for all wind directions within the plane of the flight trajectory without any mathematical simplification. Prediction errors of M3 are determined only by the quality of the input data: wind velocity, drag and lift area functions, take-off velocity, and weight. For comparing the three models, drag and lift area functions of an optimized reference jump were used. Results obtained with M2, which is much easier to handle than M3, did not deviate noticeably when compared to predictions of the reference model M3. Therefore, we suggest to use M2 in future applications. A comparison of M2 predictions with the FIS wind compensation system showed substantial discrepancies, for instance: in the first flight phase, tailwind can increase jump length, and headwind can decrease it; this is opposite of what had been anticipated before and is not considered in the current wind compensation system in ski jumping.     

Can altered muscle synergies control unimpaired gait?

Recent studies have postulated that the human motor control system recruits groups of muscles through low-dimensional motor commands, or muscle synergies. This scheme simplifies the neural control problem associated with the high-dimensional structure of the neuromuscular system. Several lines of evidence have suggested that neurological injuries, such as stroke or cerebral palsy, may reduce the dimensions that are available to the motor control system, and these altered dimensions or synergies are thought to contribute to impaired walking patterns. However, no study has investigated whether impaired low-dimensional control spaces necessarily lead to impaired walking patterns. In this study, using a two-dimensional model of walking, we developed a synergy-based control framework that can simulate the dynamics of walking. The simulation analysis showed that a synergy-based control scheme can produce well-coordinated movements of walking matching unimpaired gait. However, when the dimensions available to the controller were reduced, the simplified emergent pattern deviated from unimpaired gait. A system with two synergies, similar to those seen after neurological injury, could not produce an unimpaired walking pattern. These findings provide further evidence that altered muscle synergies can contribute to impaired gait patterns and may need to be directly addressed to improve gait after neurological injury.     

Extended foot-ankle musculoskeletal models for application in movement analysis

Multibody simulations of human motion require representative models of the anatomical structures. A model that captures the complexity of the foot is still lacking. In the present work, two detailed 3D multibody foot-ankle models generated based on CT scans using a semi-automatic tool are described. The proposed models consists of five rigid segments (talus, calcaneus, midfoot, forefoot and toes), connected by five joints (ankle, subtalar, midtarsal, tarsometatarsal and metatarsophalangeal), one with 15DOF and the other with 8DOF. The calculated kinematics of both models were evaluated using gait trials and compared against literature, both presenting realistic results. An inverse dynamic analysis was performed for the 8DOF model, again presenting feasible dynamic results.     

Molecular dynamics simulation of isothermal crystallisation of polymer chains around single polymer lamella

The isothermal crystallisation of polyethylene (PE) chains around single PE lamella in vacuum is investigated by molecular dynamic simulation. The crystallisation process is analysed in terms of the orientational order parameters, principal moments of inertia for the simulated systems. The effects of charge interactions between the polymer chains and lamella are discussed. It is found that the crystallisation process for uncharged systems can be divided into three stages: (1) adsorption, (2) orientation and (3) arrangement. The single polymer lamella changes a little during the three stages. PE chains are arranged parallel to the chain direction of the stems in the crystalline state. When considering the effect of charge interactions between the polymer chains and lamella, a different crystallisation process appears. The single polymer lamella is affected by the charged polymer chains.     

Computational Oral Absorption Simulation for Low‐Solubility Compounds

Bile micelles play an important role in oral absorption of low‐solubility compounds. Bile micelles can affect solubility, dissolution rate, and permeability. For the pH–solubility profile in bile micelles, the Henderson–Hasselbalch equation should be modified to take bile‐micelle partition into account. For the dissolution rate, in the Nernst–Brunner equation, the effective diffusion coefficient in bile‐micelle media should be used instead of the monomer diffusion coefficient. The diffusion coefficient of bile micelles is 8‐ to 18‐fold smaller than that of monomer molecules. For permeability, the effective diffusion coefficient in the unstirred water layer adjacent to the epithelial membrane, and the free fraction at the epithelial membrane surface should be taken into account. The importance of these aspects is demonstrated here using several in vivo and clinical oral‐absorption data of low‐solubility model compounds. Using the theoretical equations, the food effect on oral absorption is further discussed.