Differential diffusion quantum Monte Carlo method: determination of potential energy surfaces of molecules

A differential approach for self-optimizing diffusion Monte Carlo calculation was proposed in this paper, which is a new algorithm combining three techniques such as optimizing, diffusion and correlation sampling. This method can be used to directly compute the energy differential between two systems in the diffusion process, making the statistical error of calculation be reduced to order of 10-5 hartree, and recover about more than 80% of the correlation energy. We employed this approach to set up a potential energy surface of a molecule, used a "rigid move" model, and utilized Jacobi transformation to make energy calculation for two configurations of a molecule having good positive correlation. So, an accurate energy differential could be obtained, and the potential energy surface with good quality can be depicted. In calculation, a technique called "post-equilibrium remaining sample was set up firstly, which can save about 50% of computation expense. This novel algorithm was used to study the potenti     

Reliable radical stabilization energies from diffusion Monte Carlo calculations

We assess the performance of variational (VMC) and diffusion (DMC) quantum Monte Carlo methods for calculating the radical stabilization energies of a set of 43 carbon-centered radical species. Even using simple single-determinant trial wavefunctions, both methods perform exceptionally well, with mean absolute deviations from reference values well under the chemical accuracy standard of 1 kcal/mol. In addition, the use of DMC results in a highly concentrated spread of errors, with all 43 results within chemical accuracy at the 95% confidence level. These results indicate that DMC is an extremely reliable method for calculating radical stabilization energies and could be used as a benchmark method for larger systems in future.     

Monte Carlo simulation of diffusion effects on chain‐extension reactions

The diffusion effects on chain‐extension reactions using carboxyl‐terminated polyamide‐12 as a model reactant with bisoxazolines were investigated by the stochastic Monte Carlo method. Thus, complicated direct modeling and numerical calculations were avoided. The chain‐length dependence and detailed diffusive behavior were discussed in depth. The diffusion effects retarded the progress of chain‐extension reactions and led to lower coupling efficiency. The simulated results indicated that the diffusion effects could make the final molecular weight distributions wider. In the presence of diffusion and with the progress of the coupling efficiency, peaks in the evolution curves of the weight‐average molecular weight and valleys in the evolution curves of the polydispersity index were observed, respectively, when the coupling efficiency was low enough. These phenomena were different from those without diffusion effects and were analyzed in detail. The critical entanglement chain length had strong effects on the simulated results of the diffusion effects, especially when its value was near the average chain length. The results also showed that the effects of the reactant degradation made the molecular weight distribution of the reaction system wider and weakened the diffusion effects on the coupling reaction. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2902–2911, 2006     

Monte‐Carlo simulation of x‐ray photoelectron emission from silver

Silver 3d x‐ray photoelectron spectroscopy (XPS) spectra were simulated with the Monte‐Carlo method using an effective energy‐loss function that was derived from a reflected electron energy‐loss spectroscopy (REELS) analysis based on an extended Landau approach. After confirming that Monte‐Carlo simulation based on the use of the effective energy‐loss function can successfully describe the experimental REELS spectrum and Ag 3d XPS spectrum, we applied Monte‐Carlo simulation to predict the angular distribution of Ag 3d x‐ray photoelectrons for different x‐ray incidence angles and different photoelectron take‐off angles. The experimental photoelectron emission microscope that we are constructing was confirmed as being close to the optimum configuration, in which the x‐ray incident angle as measured from the surface normal direction is 74° and the photoelectron take‐off angle is set normal to the surface. The depth distribution functions of the Ag 3d X‐ray photoelectrons for different energy windows suggest that the photoelectron emission microscope will exhibit greater surface sensitivity for narrower photoelectron energy windows. Copyright © 2003 John Wiley & Sons, Ltd.     

Self‐diffusion coefficient of ring polymers in semidilute solution

In a topologically constraining environment the size of a flexible nonconcatenated ring polymer (macrocycles) and its dynamics are known to differ from that of linear polymers. Hence, the diffusion coefficient of ring polymers can be expected to be different from linear chains. We present here scaling arguments for the concentration and molecular weight dependence of self‐diffusion coefficient of ring polymers in semidilute solutions, and show that contrary to expectations these scaling relations are identical to what is known for linear polymers. At higher concentrations excluded volume interactions arising from possibilities of segmental overlap can become effective for large ring polymers. In this regime the diffusion coefficient of large ring polymers shows a relatively weaker dependence on concentration and molecular weight. ©2008 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 46: 2370–2379, 2008     

Quantum Monte Carlo calculations for ground and excited states

A brief overview of the diffusion quantum Monte Carlo method is given. We illustrate the application to ground‐state calculations by a study of the relative stability of carbon clusters near the crossover to fullerene stability, thereby determining the smallest stable fullerene. The application to excited states is illustrated via a study of excitonic states in small hydrogenated silicon clusters. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

精确固定节面量子Monte：Carlo差值法

本文提出了精确固定节面量子Monte Carlo差值法，这个新算法能够在精确固定节面量子Monte Carlo方法的基础上直接计算两个体系之间的能量差，且使计算结果的统计误差达到10-2 kJ/mol 数量级，获得电子相关能90%以上。我们把这个新算法应用于分子势能面的研究中，使用一个“刚性移动”模型，利用Jacobi变换使分子两个几何构型的能量计算具有很好的正相关性，因而能得到准确的能量差值，于是精确的分子势能面就可以得到。这个新算法已经被使用到BH分子基态势能曲线和H₃分子势能面的研究。这个算法还可应用于分子光谱、化学反应能量变化值等领域的研究。     

Variational Monte Carlo calculations for some cations and anions of the first‐row atoms using explicitly correlated wave functions

The variational Monte Carlo method is applied to calculate ground‐state energies of some cations and anions of the first‐row atoms. Accurate values providing between 80 and 90% of the correlation energy are obtained. Explicitly correlated wave functions including up to 42 variational parameters are used. The nondynamic correlation due to the 2s ? 2p near degeneracy effect is included by using a multideterminant wave function. The variational free parameters have been fixed by minimizing the energy that has shown to be a more convenient functional than the variance of the local energy, which is the most commonly employed method in variational Monte Carlo calculations. The energies obtained improve previous works using similar wave functions. © 2002 Wiley Periodicals, Inc.; DOI 10.1002/qua.10125     

Monte‐Carlo simulation of secondary electron emission by x‐ray irradiation—an application of x‐ray absorption near‐edge structure (XANES)

A Monte‐Carlo simulation program has been developed for describing x‐ray absorption near‐edge structure (XANES) observed by synchrotron radiation. The Monte‐Carlo simulation was applied for interpreting XANES spectroscopy on a polycrystalline Ag specimen under synchrotron irradiation with photon energy 3340–3390 eV around the absorption edge of the Ag L_α line at 3352 eV. The results clearly indicate that Monte‐Carlo simulation describes the experimental results with considerable success. Dependence of secondary electron yield on the incident angle of synchrotron radiation was also studied. Copyright © 2004 John Wiley & Sons, Ltd.     

A Hundred‐Year‐Old Experiment Re‐evaluated: Accurate Ab Initio Monte Carlo Simulations of the Melting of Radon

State‐of‐the‐art relativistic coupled‐cluster theory is used to construct many‐body potentials for the noble‐gas element radon to determine its bulk properties including the solid‐to‐liquid phase transition from parallel tempering Monte Carlo simulations through either direct sampling of the bulk or from a finite cluster approach. The calculated melting temperature are 200(3) K and 200(6) K from bulk simulations and from extrapolation of finite cluster values, respectively. This is in excellent agreement with the often debated (but widely cited) and only available value of 202 K, dating back to measurements by Gray and Ramsay in 1909.     

Segregation in zirconia: equilibrium versus non‐equilibrium segregation

The present work considers the differences between the following two phenomena: equilibrium segregation, which is a thermodynamic phenomenon where the driving force is excess interfacial energy; and non‐equilibrium segregation, which can be produced by phenomenological shocks applied to the solid surface. Fully stabilized cubic yttria‐stabilized zirconia (YSZ) is used as the exemplar for these considerations. Equilibrium segregation in YSZ results in enrichment of the surface and near‐surface layers in constituent elements, typically yttria and impurities. This segregation is an intrinsic material property and has an impact on the performance of zirconia at elevated temperatures. On the other hand, non‐equilibrium segregation leads to a complex distribution of properties (structure and concentration gradients) that are determined by the experimental procedures used rather than being a material property. However, such non‐equilibrium segregation can result in localized structural changes. The present paper also considers the effect of the gas phase on surface properties of metal oxides, including YSZ, and the surface dynamics of YSZ at temperatures below that required to reach equilibrium. Copyright © 2005 John Wiley & Sons, Ltd.     

MOCASSIN: Metropolis and kinetic Monte Carlo for solid electrolytes

The combination of density functional theory and Monte Carlo simulations is a powerful approach for the atomistic modeling of defect transport in solid electrolytes. The present contribution introduces the MOCASSIN software (Monte Carlo for Solid State Ionics) for kinetic and Metropolis Monte Carlo simulations of crystalline materials. MOCASSIN combines model building, visualization, and simulation, aiming to provide accessible MC for end users. Developed for the investigation of solid electrolytes, MOCASSIN is ideal for screening common variation parameters, such as temperature and doping fraction. The input effort is minimized using space groups for processing symmetry. The graphical interface for model building allows complex model input, including multiple mobile species, multiple migration paths, small polaron hopping, vehicle movements, multiple complex migration mechanisms, and custom interaction clusters. The software is provided free of charge for noncommercial usage.     

Free energy calculations,enhanced by a gaussian ansatz,for the “chemical work” distribution

The evaluation of the free energy is essential in molecular simulation because it is intimately related with the existence of multiphase equilibrium. Recently, it was demonstrated that it is possible to evaluate the Helmholtz free energy using a single statistical ensemble along an entire isotherm by accounting for the “chemical work” of transforming each molecule, from an interacting one, to an ideal gas. In this work, we show that it is possible to perform such a free energy perturbation over a liquid vapor phase transition. Furthermore, we investigate the link between a general free energy perturbation scheme and the novel nonequilibrium theories of Crook's and Jarzinsky. We find that for finite systems away from the thermodynamic limit the second law of thermodynamics will always be an inequality for isothermal free energy perturbations, resulting always to a dissipated work that may tend to zero only in the thermodynamic limit. The work, the heat, and the entropy produced during a thermodynamic free energy perturbation can be viewed in the context of the Crooks and Jarzinsky formalism, revealing that for a given value of the ensemble average of the “irreversible” work, the minimum entropy production corresponded to a Gaussian distribution for the histogram of the work. We propose the evaluation of the free energy difference in any free energy perturbation based scheme on the average irreversible “chemical work” minus the dissipated work that can be calculated from the variance of the distribution of the logarithm of the work histogram, within the Gaussian approximation. As a consequence, using the Gaussian ansatz for the distribution of the “chemical work,” accurate estimates for the chemical potential and the free energy of the system can be performed using much shorter simulations and avoiding the necessity of sampling the computational costly tails of the “chemical work.” For a more general free energy perturbation scheme that the Gaussian ansatz may not be valid, the free energy calculation can be expressed in terms of the moment generating function of the “chemical work” distribution. © 2014 Wiley Periodicals, Inc.     

Self-optimizing Diffusion Quantum Monte Carlo Calculation: the Potential Energy Curve of C_2

11NTRODUCTIONDiffusionquantumMonteCarlo(DMC)isoneofthesimplestofthevariousMonteCarlotechniquesavailabletosolvetheSchrodingerequation,forarecentre-viewofDMC(seeRef[1i).Foravarietyofsmallatomsandmolecules,DMChasbeenshowntobecapableofprovidinganestimateoftheground-state(nonrelativistic)energywithanaccuracycomparabletogoodqualityClcalculations,evenwhenarelativelysimpletrialwavefunctionisemployed.However,todate,thereareseveralobstaclesinaDMCcalculation:(1)BeforetheDMCcomputation,thepar…     

Diffusion Monte Carlo Calculations of Zero-Point Energies of Methanol and Deuterated Methanol

Diffusion Monte Carlo (DMC) simulations have been used to obtain quantum zero-point energies of methanol and all its isotopologs and isotopomers, using a new, accurate semi-global potential energy surface. This potential energy surface is a precise, permutationally invariant fit to 6676 ab initio energies, obtained at the CCSD(T)-F12b/aug-cc-pVDZ level of theory. Quantum zero-point energies of deuterated methanol isotopomers are very close to each other and so a simple statistical argument can be used to estimate the populations of each isotopomer at very low-temperatures. The DMC simulations also indicate that there is virtually zero probability for H/D exchange in the zero-point state. © 2019 Wiley Periodicals, Inc.     

A sparse algorithm for the evaluation of the local energy in quantum Monte Carlo

A new algorithm is presented for the sparse representation and evaluation of Slater determinants in the quantum Monte Carlo (QMC) method. The approach, combined with the use of localized orbitals in a Slater-type orbital basis set, significantly extends the size molecule that can be treated with the QMC method. Application of the algorithm to systems containing up to 390 electrons confirms that the cost of evaluating the Slater determinant scales linearly with system size.     

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 …     

Quantum Monte Carlo calculations of bond dissociation energies for some nitro and amino molecules

Bond dissociation energies (BDEs) for some nitro or amino contained prototypical molecules in energetic materials are computed by fixed‐node diffusion quantum Monte Carlo method. The nodes are determined from a Slater determinant calculated within density functional theory at the B3LYP/6‐311G** level. The possible errors, the nodal error, and the cancellation of nodal errors in calculating BDE are discussed, and the accuracy is compared with other available ab initio computations and experimental results. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2010