Exact Fixed-node Quantum Monte Carlo： Self-optimizing Procedure

In this paper, a novel exact fixed-node quantum Monte Carlo (EFNQMC) algorithm was proposed, which is a self-optimizing and self-improving procedure. In contrast to the previous EFN-QMC method, the importance function of this method is optimized synchronistically in the diffusion procedure, but not be-fore beginning the EFNQMC computation. In order to optimize the importance function, the improved steepest descent tech-nique is used, in which the step size is automatically adjustable.The procedure is quasi-Newton type and converges super linear-ly. The present method also uses a novel trial function, which has correct electron-electron and electron-nucleus cusp condi-tious. The novel EFNQMC algorithm and the novel trial func-tion are employed to calculate the energies of 1 ^1A1 state of CH2, ^1Ag state of Cs and the ground-states of H2, LiH, Li2 and H2O.     

A new algorithm for the fixed-node quantum Monte Carlo method

A novel algorithm is proposed for the fixed-node quantum Monte Carlo (FNQMC) method.In contrast to previous procedures,its "guiding function" is not optimized prior to diffusion quantum Monte Carlo (DMC) computation but synchronistically in the diffusion process The new algorithm can not only save CPU time,but also make both of the optimization and diffusion carried out according to the same sampling fashion,reaching the goal to improve each other This new optimizing procedure converges super-linearly,and thus can accelerate the particle diffusion During the diffusion process,the node of the "guiding function" changes incessantly,which is conducible to reducing the "fixed-node error" The new algorithm has been used to calculate the total energies of states X3B1 and a1A1 of CH2 as well as π-X2B1 and λ-2A1 of NH2 The singlet-triplet energy splitting (λEsT) in CH2 and π energy splitting in NH2 obtained with this present method are (45 542±1.840) and (141.644±1.589) kJ/mol,respectively The calculated     

Monte Carlo Simulation of the Assembly of bis‐Biotinylated DNA and Streptavidin

We present Monte Carlo simulations of the self‐assembly of bivalent bis‐biotinylated DNA molecules with the tetravalent biotin‐binding protein streptavidin (STV). By fitting the STV binding probabilities for the four possible valencies, the modelling correctly reproduces the dependencies of various network parameters experimentally observed in an earlier study. The combined results from the experimental and theoretical studies suggest that the binding probability for divalent STV formation is about 50 times larger than for the formation of trivalent and about 200 times larger than for tetravalent STV. In accordance with the experimental results, the modelling also indicates that the mixture of an equimolar ratio of DNA and STV leads to a maximum in size of the oligomeric DNA–STV clusters formed. Furthermore, we found a percolation transition in which the DNA cluster size increases rapidly with increasing DNA concentration resulting in the formation of a single supercluster at elevated concentrations. This behaviour coincides with the occurrence of an immobile band previously observed in electrophoretic experiments, indicating the formation of extremely large DNA–STV aggregate networks.     

Monte Carlo simulation of polarizable systems: Early rejection scheme for improving the performance of adiabatic nuclear and electronic sampling Monte Carlo simulations

An early rejection scheme for trial moves in adiabatic nuclear and electronic sampling Monte Carlo simulation (ANES-MC) of polarizable intermolecular potential models is presented. The proposed algorithm is based on Swendsen–Wang filter functions for prediction of success or failure of trial moves in Monte Carlo simulations. The goal was to reduce the amount of calculations involved in ANES-MC electronic moves, by foreseeing the success of an attempt before making those moves. The new method was employed in Gibbs ensemble Monte Carlo (GEMC) simulations of the polarizable simple point charge-fluctuating charge (SPC-FQ) model of water. The overall improvement in GEMC depends on the number of swap attempts (transfer molecules between phases) in one Monte Carlo cycle. The proposed method allows this number to increase, enhancing the chemical potential equalization. For a system with 300 SPC-FQ water molecules, for example, the fractions of early rejected transfers were about 0.9998 and 0.9994 at 373 and 423 K, respectively. This means that the transfer moves consume only a very small part of the overall computing effort, making GEMC almost equivalent to a simulation in the canonical ensemble.     

Computation of density of perfluoroalkyl methacrylates: a molecular modeling approach

Molecular dynamics simulations have been carried out on different perfluoroalkyl methacrylates to predict their densities. Density calculations on selected perfluoroalkyl methacrylates have been performed using molecular dynamics in the NPT ensemble by employing COMPASS force field. The calculated density values compared quite well with the experimental data reported in the literature.     

Structure and hydration of polycytidylic acid from the data of infrared spectroscopy, EHF dielectrometry and computer modeling

The aim of this paper is to summarize the original results concerned with the elucidation of the role of water environment in the formation of different structures of polyribocytidylic acid (poly(rC)), depending on the pH, temperature and ion content. To solve this problem, we studied the hydration of poly(rC)-K⁺ in films differing in water content, by methods of infrared spectroscopy and piezogravimetry, and in solution, using the dielectric measurements at a wavelength of 7.6 mm (extremely high frequencies — EHF). The experimental results were confirmed by a Monte Carlo simulation of the interaction between water molecules and single-strand and double-strand poly(rC) fragments in clusters of 800 water molecules. A model of hydration of double-strand complex of poly(rC)⁺-poly(rC) has been proposed based on our results and the known X-ray parameters of the complex. The obtained results and proposed structure of poly(rC)⁺-poly(rC) suggest that, the stabilization of this complex occurs due to intra- and inter-chain water bridges, together with the hydrogen bonds between neutral and protonated cytosines in pairs.     

Orbital Energy Order and the Through Bond Interaction in Homonuclear Diatomic Molecules

The homonuclear diatomic molecules are the simplest systems having both the σ framework and the lone pair orbitals n_a and _b for investigating their through space and through bond interaction. The striking orbital energy order n_g～ n_a+ n_b > n_n ～ n_a - n_b has been accounted for by the through bond interaction. However, when the p-content in the lone pair orbitals n_a and n_b decreases, one may have the reverse orbital energy order: n_g < n_g. A reverse orbital energy order has been found in F₂ and Cl₂, whose n_a and n_b are almost pure s-type atomic orbitals. The reverse order also occurs in molecule N₂ when the internuclear distance is larger man 1.5 Å. It is also found that the detail through space and through bond interaction and the eventual orbital energy order for n_g and n_u can be accounted for by the Fock operator within the localized molecular orbital space.     

Monte‐Carlo simulation of spatial distribution of characteristic x‐rays in multi‐film targets: source size of Al Kα x‐rays in W/Al film targets

A Monte‐Carlo simulation approach has been applied to describe the spatial distribution of characteristic x‐rays in W/Al film targets of different combinations of film thicknesses for the optimal design of a small‐sized x‐ray source having a high x‐ray intensity. The result has led to optimal combinations of W and Al film targets for 100 kV electrons, e.g. W(1 µm)/Al(20 µm), W(3 µm)/Al(15 µm) and W(5 µm)/Al(8 µm). These Al/W targets could be used as x‐ray sources for a medical instrument currently under development. Copyright © 2004 John Wiley & Sons, Ltd.     

DFT study and Monte Carlo simulation on proton transfers of 2-amino-2-oxazoline, 2-amino-2-thiazoline, and 2-amino-2-imidazoline in the gas phase and in water

Density functional theory (DFT) and Monte Carlo (MC) simulation with free energy perturbation (FEP) techniques have been used to study the tautomeric proton transfer reaction of 2-amino-2-oxazoline, 2-amino-2-thiazoline, and 2-amino-2-imidazoline in the gas phase and in water. Two reaction pathways were considered: the direct and water-assisted transfers. The optimized structures and thermodynamic properties of stationary points for the title reaction system in the gas phase were calculated at the B3LYP/6-311+G(d, p) level of theory. The potential energy profiles along the minimum energy path in the gas phase and in water were obtained. The study of the solvent effect of water on the proton transfer of 2-amino-2-oxozoline, 2-amino-2-thiazoline, and 2-amino-2-imidazoline indicates that water as a solvent is favorable for the water-assisted process and slows down the rate of the direct transfer pathway.