A fast parallel clustering algorithm for molecular simulation trajectories

We implemented a GPU‐powered parallel k‐centers algorithm to perform clustering on the conformations of molecular dynamics (MD) simulations. The algorithm is up to two orders of magnitude faster than the CPU implementation. We tested our algorithm on four protein MD simulation datasets ranging from the small Alanine Dipeptide to a 370‐residue Maltose Binding Protein (MBP). It is capable of grouping 250,000 conformations of the MBP into 4000 clusters within 40 seconds. To achieve this, we effectively parallelized the code on the GPU and utilize the triangle inequality of metric spaces. Furthermore, the algorithm's running time is linear with respect to the number of cluster centers. In addition, we found the triangle inequality to be less effective in higher dimensions and provide a mathematical rationale. Finally, using Alanine Dipeptide as an example, we show a strong correlation between cluster populations resulting from the k‐centers algorithm and the underlying density. © 2012 Wiley Periodicals, Inc.     

Construction of conformal mappings by generalized polarization tensors

We present a new systematic method to compute the Riemann mapping from the outside of the unit disc to the outside of a simply connected domain. We derive explicit relations between the coefficients of the Riemann mapping and the generalized polarization tensors associated with the domain. Because the generalized polarization tensors can be computed numerically, we are able to compute the coefficients of the Riemann mapping using these relations. Effectiveness of the method is validated by numerical examples. Copyright © 2014 John Wiley & Sons, Ltd.     

Fast enclosure for the minimum norm least squares solution of the matrix equation AXB = C

Fast algorithms for enclosing the minimum norm least squares solution of the matrix equation AXB = C are proposed. To develop these algorithms, theories for obtaining error bounds of numerical solutions are established. The error bounds obtained by these algorithms are verified in the sense that all the possible rounding errors have been taken into account. Techniques for accelerating the enclosure and obtaining smaller error bounds are introduced. Numerical results show the properties of the proposed algorithms. Copyright © 2015 John Wiley & Sons, Ltd.     

Are accurate computations of the 13C′ shielding feasible at the DFT level of theory?

The goal of this study is twofold. First, to investigate the relative influence of the main structural factors affecting the computation of the ¹³C′ shielding, namely, the conformation of the residue itself and the next nearest‐neighbor effects. Second, to determine whether calculation of the ¹³C′ shielding at the density functional level of theory (DFT), with an accuracy similar to that of the ¹³C^α shielding, is feasible with the existing computational resources. The DFT calculations, carried out for a large number of possible conformations of the tripeptide Ac‐G XY ‐NMe, with different combinations of X and Y residues, enable us to conclude that the accurate computation of the ¹³C′ shielding for a given residue X depends on the: (i) (?,ψ) backbone torsional angles of X ; (ii) side‐chain conformation of X ; (iii) (?,ψ) torsional angles of Y ; and (iv) identity of residue Y . Consequently, DFT‐based quantum mechanical calculations of the ¹³C′ shielding, with all these factors taken into account, are two orders of magnitude more CPU demanding than the computation, with similar accuracy, of the ¹³C^α shielding. Despite not considering the effect of the possible hydrogen bond interaction of the carbonyl oxygen, this work contributes to our general understanding of the main structural factors affecting the accurate computation of the ¹³C′ shielding in proteins and may spur significant progress in effort to develop new validation methods for protein structures. © 2013 Wiley Periodicals, Inc.     

基于遗传编程方法的微孔磷酸铝的定向合成研究

根据微孔磷酸铝合成数据的特点,针对合成数据库中合成参数过多和交叉描述等问题,利用改进的遗传编程算法对具有(6,8)元环的微孔磷酸铝合成参数进行特征提取,优化出新的复合特征来更好地描述磷酸铝合成中溶剂和模板剂的特征,并通过参数进化过程的研究,考察了模板剂和溶剂对产物生成的具体影响,从而指导具有特定结构磷酸铝的定向合成.     

Basis Sets Dependency in Constructing Spectroscopy-Accuracy Ab Initio Global Electric Dipole Moment Functions

Recently, more attention have been paid on the construction of dipole moment functions (DMF) using theoretical methods. However, the computational methods to construct DMFs are not validated as much as those for potential energy surfaces do. In this letter, using Ar...He as an example, we tested how spectroscopyaccuracy DMFs can be constructed using ab initio methods. We especially focused on the basis set dependency in this scenario, i.e., the convergence of DMF with the sizes of basis sets, basis set superposition error, and mid-bond functions. We also tested the explicitly correlated method, which converges with smaller basis sets than the conventional methods do. This work can serve as a pictorial sample of all these computational technologies behaving in the context of constructing DMFs.     

A survey on genomic data by privacy-preserving techniques perspective

Nowadays, the purpose of human genomics is widely emerging in health-related problems and also to achieve time and cost-efficient healthcare. Due to advancement in genomics and its research, development in privacy concerns is needed regarding querying, accessing and, storage and computation of the genomic data. While the genomic data is widely accessible, the privacy issues may emerge due to the untrusted third party (adversaries/researchers), they may reveal the information or strategy plans regarding the genome data of an individual when it is requested for research purposes. To mitigate this problem many privacy-preserving techniques are used along with cryptographic methods are briefly discussed. Furthermore, efficiency and accuracy in a secure and private genomic data computation are needed to be researched in future.     

Infrared/laser multi-sensor fusion and tracking based on the multi-scale model

The state estimation problem of targets detected by infrared/laser composite detection system with different sampling rates was studied in this paper. An effective state estimation algorithm based on data fusion is presented. Because sampling rate of infrared detection system is much higher than that of the laser detection system, the theory of multi-scale analysis is used to establish multi-scale model in this algorithm. At the fine scale, angle information provided by infrared detection system is used to estimate the target state through the unscented Kalman filter. It makes full use of the high frequency characteristic of infrared detection system to improve target state estimation accuracy. At the coarse scale, due to the sampling ratio of infrared and laser detection systems is an integer multiple, the angle information can be fused directly with the distance information of laser detection system to determine the target location. The fused information is served as observation, while the converted measurement Kalman filter (CMKF) is used to estimate the target state, which greatly reduces the complexity of filtering process and gets the optimal fusion estimation. The simulation results of tracking a target in 3-D space by infrared and laser detection systems demonstrate that the proposed algorithm in this paper is efficient and can obtain better performance than traditional algorithm.     

Data-fusion of high resolution X-ray CT,SEM and EDS for 3D and pseudo-3D chemical and structural characterization of sandstone

When dealing with the characterization of the structure and composition of natural stones, problems of representativeness and choice of analysis technique almost always occur. Since feature-sizes are typically spread over the nanometer to centimeter range, there is never one single technique that allows a rapid and complete characterization. Over the last few decades, high resolution X-ray CT (μ-CT) has become an invaluable tool for the 3D characterization of many materials, including natural stones. This technique has many important advantages, but there are also some limitations, including a tradeoff between resolution and sample size and a lack of chemical information. For geologists, this chemical information is of importance for the determination of minerals inside samples. We suggest a workflow for the complete chemical and structural characterization of a representative volume of a heterogeneous geological material. This workflow consists of combining information derived from CT scans at different spatial resolutions with information from scanning electron microscopy and energy-dispersive X-ray spectroscopy.