Abstract:
Quantum Monte Carlo and quantum chemistry techniques are used to investigate pseudopotential models of the lithium hydride (LiH) molecule. Interatomic potentials are calculated and tested by comparing with the experimental spectroscopic constants and well depth. Two recently-developed pseudopotentials are tested, and the effects of introducing a Li core polarization potential are investigated. The calculations are sufficiently accurate to isolate the errors from the pseudopotentials and core polarization potential. Core-valence correlation and core relaxation are found to be important in determining the interatomic potential.

Abstract:
We study the potential energy surface of the ozone molecule by means of Quantum Monte Carlo simulations based on the resonating valence bond concept. The trial wave function consists of an antisymmetrized geminal power arranged in a single-determinant that is multiplied by a Jastrow correlation factor. Whereas the determinantal part incorporates static correlation effects, the augmented real-space correlation factor accounts for the dynamics electron correlation. The accuracy of this approach is demonstrated by computing the potential energy surface for the ozone molecule in three vibrational states: symmetric, asymmetric and scissoring. We find that the employed wave function provides a detailed description of rather strongly-correlated multi-reference systems, which is in quantitative agreement with experiment.

Abstract:
Path integral Monte Carlo approach is used to study the coupled quantum dynamics of the electron and nuclei in hydrogen molecule ion. The coupling effects are demonstrated by comparing differences in adiabatic Born--Oppenheimer and non-adiabatic simulations, and inspecting projections of the full three-body dynamics onto adiabatic Born--Oppenheimer approximation. Coupling of electron and nuclear quantum dynamics is clearly seen. Nuclear pair correlation function is found to broaden by 0.040 a_0 and average bond length is larger by 0.056 a_0. Also, non-adiabatic correction to the binding energy is found. Electronic distribution is affected less, and therefore, we could say that the adiabatic approximation is better for the electron than for the nuclei.

Abstract:
In this article we present a caracterization of the vibrational spectrum of the H5+ molecule using the correlation function quantum Monte Carlo (CFQMC) method and a genetic algorithm study of the topology of the potential energy surface used in this work. The vibrational modes associated with the H3+ - H2 torsion and stretching posses very flat minima. As a consequence the fundamental frequencies corresponding to these modes are poorly described in the harmonic approximation. The vibrational frequencies obtained in this work are in good agreement with the available experimental data as well as other computational methods found in literature. In our genetic algorithm study of the potential energy surface using cartesian coordinates we have found some unexpected minima. A careful analysis shows that some of these minima are described by the same curviliniar coordinates in which the potential is described. However, they represent nonequivalent molecular geometries.

Abstract:
Kinesin is a stepping molecular motor travelling along the microtubule. It moves primarily in the plus end direction of the microtubule and occasionally in the minus-end, backward, direction. Recently, the backward steps of kinesin under different loads and temperatures start to attract interests, and the relations among them are revealed. This paper aims to theoretically understand these relations observed in experiments. After introducing a backward pathway into the previous model of the ATPase cycle of kinesin movement, the dependence of the backward movement on the load and the temperature is explored through Monte Carlo simulation. Our results agree well with previous experiments.

Abstract:
The variational Monte Carlo method is used to evaluate the ground-state energy of the confined hydrogen molecule, H_2. Accordingly, we considered the case of hydrogen molecule confined by a hard prolate spheroidal cavity when the nuclear positions are clamped at the foci (on-focus case). Also, the case of off-focus nuclei in which the two nuclei are not clamped to the foci is studied. This case provides flexibility for the treatment of the molecular properties by selecting an arbitrary size and shape of the confining spheroidal box. An accurate trial wave function depending on many variational parameters is used for this purpose. The obtained results are in good agreement with the most recent results.

Abstract:
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely: the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudo potential, and the basis set for QMC calculations. We also introduce a new strategy for the definition of the atomic orbitals involved in the Jastrow - Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.

Abstract:
By Using the variational Monte Carlo (VMC) method, we calculate the 1s{\sigma}_g state energies, the dissociation energies and the binding energies of the hydrogen molecule and its molecular ion in the presence of an aligned magnetic field regime between 0 a.u. and 10 a.u. The present calculations are based on using two types of compact and accurate trial wave functions, which are put forward for consideration in calculating energies in the absence of magnetic field. The obtained results are compared with the most recent accurate values. We conclude that the applications of VMC method can be extended successfully to cover the case of molecules under the effect of the magnetic field.

Abstract:
Spin crossover molecules have recently emerged as a family of compounds potentially useful for implementing molecular spintronics devices. The calculations of the electronic properties of such molecules is a formidable theoretical challenge as one has to describe the spin ground state of a transition metal as the legand field changes. The problem is dominated by the interplay between strong electron correlation at the transition metal site and charge delocalization over the ligands, and thus it fits into a class of problems where density functional theory may be inadequate. Furthermore, the crossover activity is extremely sensitive to environmental conditions, which are difficult to fully characterize. Here we discuss the phase transition of a prototypical spin crossover molecule as obtained with diffusion Monte Carlo simulations. We demonstrate that the ground state changes depending on whether the molecule is in the gas or in the solid phase. As our calculation provides a solid benchmark for the theory we then assess the performances of density functional theory. We find that the low spin state is always over-stabilized, not only by the (semi-)local functionals, but even by the most commonly used hybrids (such as B3LYP and PBE0). We then propose that reliable results can be obtained by using hybrid functionals containing about 50% of exact-exchange.

Abstract:
In this paper, we combine thermal effects with Landau-Zener (LZ) quantum tunneling effects in a dynamical Monte Carlo (DMC) framework to produce satisfactory magnetization curves of single-molecule magnet (SMM) systems. We use the giant spin approximation for SMM spins and consider regular lattices of SMMs with magnetic dipolar interactions (MDI). We calculate spin reversal probabilities from thermal-activated barrier hurdling, direct LZ tunneling, and thermal-assisted LZ tunnelings in the presence of sweeping magnetic fields. We do systematical DMC simulations for Mn$_{12}$ systems with various temperatures and sweeping rates. Our simulations produce clear step structures in low-temperature magnetization curves, and our results show that the thermally activated barrier hurdling becomes dominating at high temperature near 3K and the thermal-assisted tunnelings play important roles at intermediate temperature. These are consistent with corresponding experimental results on good Mn$_{12}$ samples (with less disorders) in the presence of little misalignments between the easy axis and applied magnetic fields, and therefore our magnetization curves are satisfactory. Furthermore, our DMC results show that the MDI, with the thermal effects, have important effects on the LZ tunneling processes, but both the MDI and the LZ tunneling give place to the thermal-activated barrier hurdling effect in determining the magnetization curves when the temperature is near 3K. This DMC approach can be applicable to other SMM systems, and could be used to study other properties of SMM systems.