Quantum Monte Carlo study of the Ne atom and the Ne+ ion

We report all-electron and pseudopotential calculations of the ground-stateenergies of the neutral Ne atom and the Ne+ ion using the variational and diffusion quantum Monte Carlo (DMC) methods. We investigate different levels of Slater-Jastrow trial wave function: (i) using Hartree-Fock orbitals, (ii) using orbitals optimized within a Monte Carlo procedure in the presence of a Jastrow factor, and (iii) including backflow correlations in the wave function. Small reductions in the total energy are obtained by optimizing the orbitals, while more significant reductions are obtained by incorporating backflow correlations. We study the finite-time-step and fixed-node biases in the DMC energy and show that there is a strong tendency for these errors to cancel when the first ionization potential (IP) is calculated. DMC gives highly accurate values for the IP of Ne at all the levels of trial wave function that we have considered

Preliminary Investigation of Novel Bone Graft Substitutes based on Strontium-Calcium-Zinc-Silicate Glasses

Bone graft procedures typically require surgeons to harvest bone from a second site on a given patient (Autograft) before repairing a bone defect. However, this results in increased surgical time, excessive blood loss and a significant increase in pain. In this context a synthetic bone graft with excellent histocompatibility, built in antibacterial efficacy and the ability to regenerate healthy tissue in place of diseased tissue would be a significant step forward relative to current state of the art philosophies. We developed a range of calcium-strontium-zinc-silicate glass-based bone grafts and characterized their structure and physical properties, then evaluated their in vitro cytotoxicity and in vivo biocompatibility using standardized models from the literature. A graft (designated BT109) of composition 0.28SrO/0.32ZnO/0.40 SiO2 (mol fraction) was the best performing formulation in vitro shown to induce extremely mild cytopathic effects (cell viability up to 95%) in comparison with the commercially available bone graft Novabone® (cell viability of up to 72%). Supplementary to this, the grafts were examined using the standard rat femur healing model on healthy Wister rats. All grafts were shown to be equally well tolerated in bone tissue and new bone was seen in close apposition to implanted particles with no evidence of an inflammatory response within bone. Complimentary to this BT109 was implanted into the femurs of ovariectomized rats to monitor the response of osteoporotic tissue to the bone grafts. The results from this experiment indicate that the novel grafts perform equally well in osteoporotic tissue as in healthy tissue, which is encouraging given that bone response to implants is usually diminished in ovariectomized rats. In conclusion these materials exhibit significant potential as synthetic bone grafts to warrant further investigation and optimisation. © 2008 Springer Science+Business Media, LLC

Effects of Adding Sodium and Fluoride Ions to Glass Ionomer on its Interactions with Sodium Fluoride Solution

This investigates the effects of the addition of Na and F ions to a glass ionomer cement in which those ions are not inherently present on its interactions with dilute (0.2%) NaF solution. Both the effect of the solution on the cement\u27s surface morphology and the effect of the cement on the solution in terms of take up of Na+ and F- and of change in pH are to be investigated. These results are to be compared to previous results obtained with glasses which contained both, one, or neither of the ions as components of their glasses. NaF (1.3% by weight in the mixed cement) was added to the powder components of a glass ionomer based on LG30 glass (which contains Al, Si, Ca, P, and O only). Discs of cement were set in moulds at 37°C for 1h then stored in water at 37°C for 3 days. Each test disc was then immersed in 10ml 0.2% NaF solution whereas controls remained immersed in water (N=3 for test and control). Test and control disc surfaces were assessed both qualitatively by electron microscopy and quantitatively by linear profilometry (R(a) values). Potentiometry was used to measure solution pH and Na and F concentrations using a pH electrode and suitable ion selective electrodes both before and after cement immersion. The surface of test specimens was subject considerable disruption with the polysalt cement matrix being removed and residual glass particles being disclosed. The controls showed no such disruption. This effect was reflected in a significant difference of R(a). Such an effect was not shown by test and control surfaces of LG30 but a similar effect was to that shown by LG26 (which contains F as a glass component). Solution pH changed by 1 unit which was much more than the change shown by LG30 or LG26 but is similar to that of AH2 and MP4 cements which both contain Na. The Na and F uptake was much lower than for LG30 whereas that of LG26 was higher than LG30. The Na:F ratio was 0.29:1 compared to 1.26:1 for LG30 (LG26=1.01:1, AH2=1.02:1, MP4=1.04:1). Fluoride addition to a F-free glass ionomer renders it vulnerable to surface disruption by NaF solution showing that fluoride complexes produced in glass dissolution are not necessarily involved in this process. Sodium addition to a Na-free glass ionomer confirms the role of this cement in enhancing pH change in NaF solution. The level of uptake of F- from a NaF solution in much lower than that for the F-free glass ionomer which shows there is no direct relationship between F- uptake and surface disruption. The ratio of Na:F uptake is below 0.3:1, but the pH change is similar to cements where the ratio is close to unity which indicates that F-/OH- interchange is not a significant mechanism even when anion/cation uptake is not balanced. Copyright (C) 2000

The impact of population-based faecal occult blood test screening on colorectal cancer mortality:a matched cohort study

BACKGROUND: Randomised trials show reduced colorectal cancer (CRC) mortality with faecal occult blood testing (FOBT). This outcome is now examined in a routine, population-based, screening programme. METHODS: Three biennial rounds of the UK CRC screening pilot were completed in Scotland (2000–2007) before the roll out of a national programme. All residents (50–69 years) in the three pilot Health Boards were invited for screening. They received a FOBT test by post to complete at home and return for analysis. Positive tests were followed up with colonoscopy. Controls, selected from non-pilot Health Boards, were matched by age, gender, and deprivation and assigned the invitation date of matched invitee. Follow-up was from invitation date to 31 December 2009 or date of death if earlier. RESULTS: There were 379 655 people in each group (median age 55.6 years, 51.6% male). Participation was 60.6%. There were 961 (0.25%) CRC deaths in invitees, 1056 (0.28%) in controls, rate ratio (RR) 0.90 (95% confidence interval (CI) 0.83–0.99) overall and 0.73 (95% CI 0.65–0.82) for participants. Non-participants had increased CRC mortality compared with controls, RR 1.21 (95% CI 1.06–1.38). CONCLUSION: There was a 10% relative reduction in CRC mortality in a routine screening programme, rising to 27% in participants

Muonium as a hydrogen analogue in silicon and germanium; quantum effects and hyperfine parameters

We report a first-principles theoretical study of hyperfine interactions, zero-point effects and defect energetics of muonium and hydrogen impurities in silicon and germanium. The spin-polarized density functional method is used, with the crystalline orbitals expanded in all-electron Gaussian basis sets. The behaviour of hydrogen and muonium impurities at both the tetrahedral and bond-centred sites is investigated within a supercell approximation. To describe the zero-point motion of the impurities, a double adiabatic approximation is employed in which the electron, muon/proton and host lattice degrees of freedom are decoupled. Within this approximation the relaxation of the atoms of the host lattice may differ for the muon and proton, although in practice the difference is found to be slight. With the inclusion of zero-point motion the tetrahedral site is energetically preferred over the bond-centred site in both silicon and germanium. The hyperfine and superhyperfine parameters, calculated as averages over the motion of the muon, agree reasonably well with the available data from muon spin resonance experiments.Comment: 20 pages, including 9 figures. To appear in Phys. Rev.

Unrestricted Hartree-Fock theory of Wigner crystals

We demonstrate that unrestricted Hartree-Fock theory applied to electrons in a uniform potential has stable Wigner crystal solutions for $r_s \geq 1.44$ in two dimensions and $r_s \geq 4.5$ in three dimensions. The correlation energies of the Wigner crystal phases are considerably smaller than those of the fluid phases at the same density.Comment: 4 pages, 5 figure

Quantum Monte Carlo calculations of the one-body density matrix and excitation energies of silicon

Quantum Monte Carlo (QMC) techniques are used to calculate the one-body density matrix and excitation energies for the valence electrons of bulk silicon. The one-body density matrix and energies are obtained from a Slater-Jastrow wave function with a determinant of local density approximation (LDA) orbitals. The QMC density matrix evaluated in a basis of LDA orbitals is strongly diagonally dominant. The natural orbitals obtained by diagonalizing the QMC density matrix resemble the LDA orbitals very closely. Replacing the determinant of LDA orbitals in the wave function by a determinant of natural orbitals makes no significant difference to the quality of the wave function's nodal surface, leaving the diffusion Monte Carlo energy unchanged. The Extended Koopmans' Theorem for correlated wave functions is used to calculate excitation energies for silicon, which are in reasonable agreement with the available experimental data. A diagonal approximation to the theorem, evaluated in the basis of LDA orbitals, works quite well for both the quasihole and quasielectron states. We have found that this approximation has an advantageous scaling with system size, allowing more efficient studies of larger systems.Comment: 13 pages, 4 figures. To appear in Phys. Rev.

Ground-state properties of rutile: electron-correlation effects

Electron-correlation effects on cohesive energy, lattice constant and bulk compressibility of rutile are calculated using an ab-initio scheme. A competition between the two groups of partially covalent Ti-O bonds is the reason that the correlation energy does not change linearly with deviations from the equilibrium geometry, but is dominated by quadratic terms instead. As a consequence, the Hartree-Fock lattice constants are close to the experimental ones, while the compressibility is strongly renormalized by electronic correlations.Comment: 1 figure to appear in Phys. Rev.

Structural, Electronic, and Magnetic Properties of MnO

We calculate the structural, electronic, and magnetic properties of MnO from first principles, using the full-potential linearized augmented planewave method, with both local-density and generalized-gradient approximations to exchange and correlation. We find the ground state to be of rhombohedrally distorted B1 structure with compression along the [111] direction, antiferromagnetic with type-II ordering, and insulating, consistent with experiment. We show that the distortion can be understood in terms of a Heisenberg model with distance dependent nearest-neighbor and next-nearest-neighbor couplings determined from first principles. Finally, we show that magnetic ordering can induce significant charge anisotropy, and give predictions for electric field gradients in the ground-state rhombohedrally distorted structure.Comment: Submitted to Physical Review B. Replaced: regenerated figures to resolve font problems in automatically generated pd

Carbon clusters near the crossover to fullerene stability

The thermodynamic stability of structural isomers of $\mathrm{C}_{24}$, $\mathrm{C}_{26}$, $\mathrm{C}_{28}$ and $\mathrm{C}_{32}$, including fullerenes, is studied using density functional and quantum Monte Carlo methods. The energetic ordering of the different isomers depends sensitively on the treatment of electron correlation. Fixed-node diffusion quantum Monte Carlo calculations predict that a $\mathrm{C}_{24}$ isomer is the smallest stable graphitic fragment and that the smallest stable fullerenes are the $\mathrm{C}_{26}$ and $\mathrm{C}_{28}$ clusters with $\mathrm{C}_{2v}$ and $\mathrm{T}_{d}$ symmetry, respectively. These results support proposals that a $\mathrm{C}_{28}$ solid could be synthesized by cluster deposition.Comment: 4 pages, includes 4 figures. For additional graphics, online paper and related information see http://www.tcm.phy.cam.ac.uk/~prck