Optimization of the valence energy variance of the CuH molecule

We developed the concept of split-'t to deal with the large molecules (in terms of the number of electrons and nuclear charge Z). This naturally leads to partitioning the local energy into components due to each electron shell. The minimization of the variation of the valence shell local energy is used to optimize a simple two parameter CuH wave function. Molecular properties (spectroscopic constants and the dipole moment) are calculated for the optimized and nearly optimized wave functions using the Variational Quantum Monte Carlo method. Our best results are comparable to those from the single and double configuration interaction (SDCI) method.

Variational Monte Carlo study of core-valence separation schemes for first-row atoms and positive ions /

All-electron partitioning of wave functions into products ^core^vai of core and valence parts in orbital space results in the loss of core-valence antisymmetry, uncorrelation of motion of core and valence electrons, and core-valence overlap. These effects are studied with the variational Monte Carlo method using appropriately designed wave functions for the first-row atoms and positive ions. It is shown that the loss of antisymmetry with respect to interchange of core and valence electrons is a dominant effect which increases rapidly through the row, while the effect of core-valence uncorrelation is generally smaller. Orthogonality of the core and valence parts partially substitutes the exclusion principle and is absolutely necessary for meaningful calculations with partitioned wave functions. Core-valence overlap may lead to nonsensical values of the total energy. It has been found that even relatively crude core-valence partitioned wave functions generally can estimate ionization potentials with better accuracy than that of the traditional, non-partitioned ones, provided that they achieve maximum separation (independence) of core and valence shells accompanied by high internal flexibility of ^core and Wvai- Our best core-valence partitioned wave function of that kind estimates the IP's with an accuracy comparable to the most accurate theoretical determinations in the literature.

Variational Monte Carlo estimation of the dissociation energy of CuH using correlated sampling

A new approach to treating large Z systems by quantum Monte Carlo has been developed. It naturally leads to notion of the 'valence energy'. Possibilities of the new approach has been explored by optimizing the wave function for CuH and Cu and computing dissociation energy and dipole moment of CuH using variational Monte Carlo. The dissociation energy obtained is about 40% smaller than the experimental value; the method is comparable with SCF and simple pseudopotential calculations. The dipole moment differs from the best theoretical estimate by about 50% what is again comparable with other methods (Complete Active Space SCF and pseudopotential methods).

Exact property estimation from diffusion Monte Carlo with minimal stochastic reconfiguration /

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.

Significance of DNA base excision repair in the maintenance of mitochondrial function in Saccharomyces cerevisiae /

Mitochondria have an important role in cell metabolism, being the major site of ATP production via oxidative phosphorylation (OXPHOS). Accumulation of mtDNA mutations have been linked to the development of respiratory dysfunction, apoptosis, and aging. Base excision repair (BER) is the major and the only certain repair pathway existing in mitochondria that is in responsible for removing and repairing various base modifications as well as abasic sites (AP sites). In this research, Saccharomyces cerevisiae (S. cerevisiae) BER gene knockout strains, including 3 single DNA glycosylase gene knockout strains and Ap endonuclease (Apn 1 p) knockout strain were used to examine the importance of this DNA repair pathway to the maintenance of respiratory function. Here, I show that individual DNA glycosylases are nonessential in maintenance of normal function in yeast mitochondria, corroborating with previous research in mammalian experimental models. The yeast strain lacking Apn 1 p activity exhibits respiratory deficits, including inefficient and significantly low intracellular ATP level, which maybe due to partial uncoupling of OXPHOS. Growth of this yeast strain on respiratory medium is inhibited, but no evidence was found for increased ROS level in Apn 1 p mitochondria. This strain also shows an increased cell size, and this observation combined with an uncoupled OXPHOS may indicate a premature aging in the Apnlp knockout strain, but more evidence is needed to support this hypothesis. However, the BER is necessary for maintenance of mitochondrial function in respiring S.cerevisiae.

An experimental investigation of dechanneling in copper single crystals

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.

On the formulation and the calculation of the harmonic contributions to the Debye-Waller factor in metals (sodium)

The algebraic expressions for the anharmonic contributions to the Debye-Waller factor up to 0(A ) and 0 L% ) £ where ^ is the scattering wave-vector] have been derived in a form suitable for cubic metals with small ion cores where the interatomic potential extends to many neighbours. This has been achieved in terms of various wave-vector dependent tensors, following the work of Shukla and Taylor (1974) on the cubic anharmonic Helmholtz free energy. The contribution to the various wave-vector dependent tensors from the coulomb and the electron-ion terms in the interatomic metallic potential has been obtained by the Ewald procedure. All the restricted multiple whole B r i l l o u i n zone (B.Z.) sums are reduced to single whole B.Z. sums by using the plane wave representation of the delta function. These single whole B.Z. sums are further reduced to the •%?? portion of the B.Z. following Shukla and Wilk (1974) and Shukla and Taylor (1974). Numerical calculations have been performed for sodium where the Born-Mayer term in the interatomic potential has been neglected because i t is small £ Vosko (1964)3 • *n o^er to compare our calculated results with the experimental results of Dawton (1937), we have also calculated the r a t io of the intensities at different temperatures for the lowest five reflections (110), (200), (220), (310) and (400) . Our calculated quasi-harmonic results agree reasonably well with the experimental results at temperatures (T) of the order of the Debye temperature ( 0 ). For T » © ^ 9 our calculated anharmonic results are found to be in good agreement with the experimental results.The anomalous terms in the Debye-Waller factor are found not to be negligible for certain reflections even for T ^ ©^ . At temperature T yy Op 9 where the temperature is of the order of the melting temperature (Xm) » "the anomalous terms are found to be important almost for all the f i ve reflections.

Optimization of trial wave functions for use in quantum Monte Carlo with application to LiH

Methods for both partial and full optimization of wavefunction parameters are explored, and these are applied to the LiH molecule. A partial optimization can be easily performed with little difficulty. But to perform a full optimization we must avoid a wrong minimum, and deal with linear-dependency, time step-dependency and ensemble-dependency problems. Five basis sets are examined. The optimized wavefunction with a 3-function set gives a variational energy of -7.998 + 0.005 a.u., which is comparable to that (-7.990 + 0.003) 1 of Reynold's unoptimized \fin ( a double-~ set of eight functions). The optimized wavefunction with a double~ plus 3dz2 set gives ari energy of -8.052 + 0.003 a.u., which is comparable with the fixed-node energy (-8.059 + 0.004)1 of the \fin. The optimized double-~ function itself gives an energy of -8.049 + 0.002 a.u. Each number above was obtained on a Bourrghs 7900 mainframe computer with 14 -15 hrs CPU time.

THE INFLUENCE OF LENGTH CHANGE SPEED AND DIRECTION ON DYNAMIC FUNCTION POTENTIATION IN FAST MOUSE MUSCLE

The purpose of this study was to test the hypothesis that the potentiation of dynamic function was dependent upon both length change speed and direction. Mouse EDL was cycled in vitro (25º C) about optimal length (Lo) with constant peak strain (± 2.5% Lo) at 1.5, 3.3 and 6.9 Hz before and after a conditioning stimulus. A single pulse was applied during shortening or lengthening and peak dynamic (concentric or eccentric) forces were assessed at Lo. Stimulation increased peak concentric force at all frequencies (range: 19 ± 1 to 30 ± 2%) but this increase was proportional to shortening speed, as were the related changes to concentric work/power (range: -15 ± 1 to 39 ± 1 %). In contrast, stimulation did not increase eccentric force, work or power at any frequency. Thus, results reveal a unique hysteresis like effect for the potentiation of dynamic output wherein concentric and eccentric forces increase and decrease, respectively, with work cycle frequency.

THE INFLUENCE OF LENGTH CHANGE SPEED AND DIRECTION ON DYNAMIC FUNCTION POTENTIATION IN FAST MOUSE MUSCLE

The purpose of this study was to test the hypothesis that the potentiation of dynamic function was dependent upon both length change speed and direction. Mouse EDL was cycled in vitro (250 C) about optimal length (Lo) with constant peak strain (± 2.5% Lo) at 1.5,3.3 and 6.9 Hz before and after a conditioning stimulus. A single pulse was applied during shortening or lengthening and peak dynamic (concentric or eccentric) forces were assessed at Lo. Stimulation increased peak concentric force at all frequencies (range: 19±1 to 30 ± 2%) but this increase was proportional to shortening speed, as were the related changes to concentric work/power (range: -15 ± 1 to 39 ± 1 %). In contrast, stimulation did not increase eccentric force, work or power at any frequency. Thus, results reveal a unique hysteresis like effect for the potentiation of dynamic output wherein concentric and eccentric forces increase and decrease, respectively, with work cycle frequency.