3 resultados para Time-Fractional Multiterm Diffusion Equation
em Brock University, Canada
Resumo:
This investigation comprises a comparison of experimental and theoretical dechanneling of MeV protons in copper single crystals. Dechanneling results when an ion's transverse energy increases to the value where the ion can undergo small impact parameter collisions with individual atoms. Depth dependent dechanneling rates were determined as functions of lattice temperature, ion beam energy and crystal axis orientation. Ion beam energies were IMeV and 2MeV,temperatures ranged from 35 K to 280 K and the experiment was carried out along both the (lOa) and <110) axes. Experimental data took the form of aligned and random Rutherford backscattered energy spectra. Dechanneling rates were extracted from these spectra using a single scattering theory that took explicit account of the different stopping powers experienced by channeled and dechanneled ions and also included a correction factor to take into account multiple scattering effects along the ion's trajectory. The assumption of statistical equilibrium and small angle scattering of the channeled ions allows a description of dechanneling in terms of the solution of a diffusion like equation which contains a so called diffusion function. The diffusion function is shown to be related to the increase in average transverse energy. Theoretical treatments of increase in average transverse energy due to collisions of projectiles with channel electrons and thermal perturbations in the lattice potential are reviewed. Using the diffusion equation and the electron density in the channel centre as a fitting parameter dechanneling rates are extracted. Excellent agreement between theory and experiment has been demonstrated. Electron densities determined in the fitting procedure appear to be realistic. The surface parameters show themselves to be good indicators of the quality of the crystal.
Resumo:
We examined three different algorithms used in diffusion Monte Carlo (DMC) to study their precisions and accuracies in predicting properties of isolated atoms, which are H atom ground state, Be atom ground state and H atom first excited state. All three algorithms — basic DMC, minimal stochastic reconfiguration DMC, and pure DMC, each with future-walking, are successfully impletmented in ground state energy and simple moments calculations with satisfactory results. Pure diffusion Monte Carlo with future-walking algorithm is proven to be the simplest approach with the least variance. Polarizabilities for Be atom ground state and H atom first excited state are not satisfactorily estimated in the infinitesimal differentiation approach. Likewise, an approach using the finite field approximation with an unperturbed wavefunction for the latter system also fails. However, accurate estimations for the a-polarizabilities are obtained by using wavefunctions that come from the time-independent perturbation theory. This suggests the flaw in our approach to polarizability estimation for these difficult cases rests with our having assumed the trial function is unaffected by infinitesimal perturbations in the Hamiltonian.
Resumo:
Our objective is to develop a diffusion Monte Carlo (DMC) algorithm to estimate the exact expectation values, ($o|^|^o), of multiplicative operators, such as polarizabilities and high-order hyperpolarizabilities, for isolated atoms and molecules. The existing forward-walking pure diffusion Monte Carlo (FW-PDMC) algorithm which attempts this has a serious bias. On the other hand, the DMC algorithm with minimal stochastic reconfiguration provides unbiased estimates of the energies, but the expectation values ($o|^|^) are contaminated by ^, an user specified, approximate wave function, when A does not commute with the Hamiltonian. We modified the latter algorithm to obtain the exact expectation values for these operators, while at the same time eliminating the bias. To compare the efficiency of FW-PDMC and the modified DMC algorithms we calculated simple properties of the H atom, such as various functions of coordinates and polarizabilities. Using three non-exact wave functions, one of moderate quality and the others very crude, in each case the results are within statistical error of the exact values.