Fractional occupation numbers and spin density functional calculations of degenerate systems

We implemented ab initio self‐consistent field (SCF) fractional occupation numbers (FON) calculation with Dunlap's interpolation scheme for the twisted ethylene, which is a prototype molecule of a σ–π biradical system. The calculational results are compared with those of complete‐active‐space (CAS) SCF and spin‐unrestricted Kohn–Sham (UKS) calculations on potential surfaces, occupation numbers of natural orbitals, and correlation entropies. It was found that the UKS methods gave similar results to CASSCF, while the FON solutions appeared in only the nearly complete degenerate region. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 93: 317–323, 2003     

Dynamics of the formation of two-dimensional ordered structures

The growth of ordered domains in lattice gas models, which occurs after the system is quenched from infinite temperature to a state below the critical temperatureT _c, is studied by Monte Carlo simulation. For a square lattice with repulsion between nearest and next-nearest neighbors, which in equilibrium exhibits fourfold degenerate (2×1) superstructures, the time-dependent energy E(t), domain size L(t), and structure functionS(q, t) are obtained, both for Glauber dynamics (no conservation law) and the case with conserved density (Kawasaki dynamics). At late times the energy excess and halfwidth of the structure factor decrease proportional tot ^–x, whileL(t) t ^x, where the exponent x=1/2 for Glauber dynamics and x1/3 for Kawasaki dynamics. In addition, the structure factor satisfies a scaling lawS(k,t)=t ^2xS(kt^x). The smaller exponent for the conserved density case is traced back to the excess density contained in the walls between ordered domains which must be redistributed during growth. Quenches toT>T _c, T=T_c (where we estimate dynamic critical exponents) andT=0 are also considered. In the latter case, the system becomes frozen in a glasslike domain pattern far from equilibrium when using Kawasaki dynamics. The generalization of our results to other lattices and structures also is briefly discussed.     

Succession of stacking disorder in hard-sphere crystal under gravity by Monte Carlo simulation

We have performed Monte Carlo simulations of hard spheres under gravity. Vertical boundaries are hard walls, which are well separated with each other. On the other hand, the periodic boundary condition is imposed in the horizontal direction. While we previously reported enhancement of crystallinity as well as crystallization due to gravity, we present here the results that demonstrate the succession of a defect. In case that the crystal formed at the bottom of the system includes kinds of stacking disorders for the (0 0 1) growth, twin band structure develops as mediated by a stacking disorder succeeded in the crystal formed in the fluid region which lies on the bottom crystal. In case that the stacking structure along horizontal direction changes from the (1 1 1) stacking to the (0 0 1) stacking, twin band structure in the (0 0 1) stacking region develops as succeeded in the crystal transformed. The twin band structure also becomes large with its upward growth.     

Quasi‐random integration in quantum chemistry: Efficiency,stability, and application to the study of small atoms and molecules constrained in spherical boxes

This project consists of two parts. In the first part, a series of test calculations is performed to verify that the integrals involved in the determination of atomic and molecular properties by standard self‐consistent field (SCF) methods can be obtained through Halton, Korobov, or Hammersley quasi‐random integration procedures. Through these calculations, we confirm that all three methods lead to results that meet the levels of precision required for their use in the calculation of properties of small atoms or molecules at least at a Hartree–Fock level. Moreover, we have ensured that the efficiency of quasi‐random integration methods that we have tested is Halton=Korobov>Hammersley?pseudo‐random. We also find that these results are comparable to those yielded by ordinary Monte Carlo (pseudo‐random) integration, with a calculation effort of two orders of smaller magnitude. The second part, which would not have been possible without the integration method previously analyzed, contains a first study of atoms constrained in spherical boxes through SCF calculations with basis functions adapted to the features of the problem: Slater‐type orbitals (STOs) trimmed by multiplying them by a function that yields 1 for 0 < r < (R‐δ), polynomial values for (R‐δ) < r < R and null for r > R, R being the radius of the box and δ a variationally determined interval. As a result, we obtain a equation of state for electrons of small systems, valid just in the limit of low temperatures, but fairly simple. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

The index of ideality of correlation: improvement of models for toxicity to algae

Toxicity to algae is important characteristic of substances from ecologic point of view. The CORAL software (http://www.insilico.eu/coral) gives possibility to build up model of toxicity to algae using data on the molecular architecture and experimental toxicity, without additional data on physicochemical and/or biochemical parameters. Considerable improvement of the model is observed in the case of using the index of ideality of correlation (IIC) in the role of additional criterion of predictive potential. The IIC is calculated with using of the correlation coefficient between experimental and calculated values of endpoint for the calibration set, with taking into account the positive and negative dispersions between experimental and calculated values. The best model calculated with use the IIC is characterized (the validation set) by n?=?50, r²?=?0.947, RMSE?=?0.401 whereas, model calculated without use the IIC is characterized by n?=?50, r²?=?0.805, and RMSE?=?0.539. The suggested models are built up in accordance to five OECD principles.

     

Effect of chemical potential on the computer simulation of hydrogen storage in single walled carbon nanotubes

Hydrogen is a kind of clean, sustainable and renewable energy carrier. Of the problems to be solved for the utilization of hydrogen energy, how to store and transport hydrogen has been given high priority on the research agenda. Recently, carbon nanotubes (CNTs) were reported to be very promising candidates for hydrogen uptake[1], which may have possibility to satisfy the benchmark set by the US Department of Energy (DOE) Hydrogen Plan for fuel cell powered vehicles: a gravimetric density …     

Evaluation of the use of consensus values in proficiency testing programmes

Proficiency testing (PT) is an essential tool used by laboratory accreditation bodies to assess the competency of laboratories. Because of limited resources of PT providers or for other reasons, the assigned reference value used in the calculation of z-score values has usually been derived from some sort of consensus value obtained by central tendency estimators such as the arithmetic mean or robust mean. However, if the assigned reference value deviates significantly from the ‘true value’ of the analyte in the test material, laboratories’ performance will be evaluated incorrectly. This paper evaluates the use of consensus values in proficiency testing programmes using the Monte Carlo simulation technique. The results indicated that the deviation of the assigned value from the true value could be as large as 40%, depending on the parameters of the proficiency testing programmes under investigation such as sample homogeneity, number of participant laboratories, concentration level, method precision and laboratory bias. To study how these parameters affect the degree of discrepancy between the consensus value and the true value, a fractional factorial design was also applied. The findings indicate that the number of participating laboratories and the distribution of laboratory bias were the prime two factors affecting the deviation of the consensus value from the true value.     

Effect of conformational rigidity on the glass transition and initial decomposition temperatures of polyimides

The dependence of the glass transition and initial decomposition temperatures of polyperyleneimide and polynaphthylimide on the conformational rigidity was studied by the Monte Carlo, Kuhn segment, and quantum chemical AM1 methods. The corresponding linear plots can be used for estimation of the glass transition and initial decomposition temperatures when experimental determination is difficult.     

Rigid-body dynamics in the isothermal-isobaric ensemble: a test on the accuracy and computational efficiency

We have developed a time-reversible rigid-body (rRB) molecular dynamics algorithm in the isothermal-isobaric (NPT) ensemble. The algorithm is an extension of rigid-body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non-Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state-of-the-art constraint-dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time-step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS-MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL.