43 research outputs found
Calculation of the energies of the multideterminant states of the nitrogen vacancy center in diamond with quantum Monte Carlo
Certain point defects in solids can efficiently be used as qubits for
applications in quantum technology. They have spin states that are
initializable, readable, robust, and can be manipulated optically. New
theoretical methods are needed to find the best host materials and defect
configurations. Most methods proposed so far rely either on cluster models or
restrict the many-body treatment of the defects to a subspace of
single-particle orbitals. We explore best practices and theory for the use of
quantum Monte Carlo to predict the excitation spectra for spin defects, by
using the negatively charged nitrogen vacancy (NV) center in diamond as a
test system. Quantum Monte Carlo can be used to explicitly simulate electronic
correlations with larger systems and sets of orbitals than previous methods due
to favourable scaling with respect to system size and computing power. We
consider different trial wave functions for variational and diffusion Monte
Carlo methods, explore the nodal surface errors of the ground and excited state
wave functions and study whether the variational principle holds for the
excited states. We compute the vertical excitation energies in different
simulation cell sizes and extrapolate to infinite system size, and include
backflow corrections to the extrapolated energies. The final results for
vertical excitation energies are found to overestimate the experimental
estimates, but the triplet-to-triplet and singlet-to-singlet transitions are
accurate against experiment. Finally, we list further developments for QMC
needed to address the problem of accurately predicting structural and spin
properties of the solid-state defects
Split Ga vacancies : Abundant defects in β-Ga2O3
We have applied positron annihilation spectroscopy to study a wide range of β-Ga2O3bulk crystals and thin films with various doping levels. The Doppler broadening of the 511 keV positron-electron annihilation line exhibits colossal anisotropy compared to other three-dimensional crystalline semiconductors. State-of-the-art theoretical calculations of the positron characteristics in the β-Ga2O3lattice reveal that the positron state is effectively 1-dimensional, giving rise to strong anisotropy. Strongly relaxed split Ga vacancies are found to exhibit even stronger anisotropy and to dominate the positron annihilation signals in almost all experiments. The evidence leads to the conclusion that split Ga vacancies are abundant, with concentration of 1018 cm-3 or more, in β-Ga2O3samples irrespective of conductivity.Peer reviewe
Identifying vacancy complexes in compound semiconductors with positron annihilation spectroscopy: a case study of InN
We present a comprehensive study of vacancy and vacancy-impurity complexes in
InN combining positron annihilation spectroscopy and ab-initio calculations.
Positron densities and annihilation characteristics of common vacancy-type
defects are calculated using density functional theory and the feasibility of
their experimental detection and distinction with positron annihilation methods
is discussed. The computational results are compared to positron lifetime and
conventional as well as coincidence Doppler broadening measurements of several
representative InN samples. The particular dominant vacancy-type positron traps
are identified and their characteristic positron lifetimes, Doppler ratio
curves and lineshape parameters determined. We find that In vacancies and their
complexes with N vacancies or impurities act as efficient positron traps,
inducing distinct changes in the annihilation parameters compared to the InN
lattice. Neutral or positively charged N vacancies and pure N vacancy complexes
on the other hand do not trap positrons. The predominantly introduced positron
trap in irradiated InN is identified as the isolated In vacancy, while in
as-grown InN layers In vacancies do not occur isolated but complexed with one
or more N vacancies. The number of N vacancies per In vacancy in these
complexes is found to increase from the near surface region towards the
layer-substrate interface.Comment: 10 pages, 6 figure
Exchange-correlation potentials for inhomogeneous electron systems in two dimensions from exact diagonalization: comparison with the local-spin-density approximation
We consider electronic exchange and correlation effects in density-functional
calculations of two-dimensional systems. Starting from wave function
calculations of total energies and electron densities of inhomogeneous model
systems, we derive corresponding exchange-correlation potentials and energies.
We compare these with predictions of the local-spin-density approximation and
discuss its accuracy. Our data will be useful as reference data in testing,
comparing and parametrizing exchange and correlation functionals for
two-dimensional electronic systems.Comment: Submitted to Physical Review B on January 3, 2012. Second revised
version submitted on April 13, 201
Defect identification in semiconductors with positron annihilation: Experiment and theory
Peer reviewe
PGET Monte Carlo simulations using Serpent
Since 2017, over 100 spent nuclear fuel assemblies at the Finnish nuclear power plants have been imaged with the Passive Gamma Emission Tomography (PGET) device in preparation of the implementation of PGET in the safeguards infrastructure of the Finnish geological repository. In order to increase understanding of the PGET method and guide its further development, we have recently implemented PGET in Serpent, a widely-used neutron and photon transport Monte Carlo simulation code. We will discuss the major aspects of this implementation and illustrate the usefulness of the simulations with a few examples. The PGET device as used in the measurements (which was developed under the guidance of IAEA and is approved for safeguards inspections) was implemented in a very realistic way based on its technical drawings. The simulation produces sinograms in user-defined energy windows as well as the uncertainty on these sinograms. Tomographic images are then reconstructed using the exact same algorithm as used for the measured data. A dedicated variance reduction scheme was implemented, increasing the computational efficiency by about a factor of 30. The simulation of the PGET response at one angular measurement position for 1 billion primary photons takes a few hours on a single 40-core node. The 1-sigma uncertainty in the highest intensity sinogram pixels is about a few percent. Aiming at improving the imaging of VVER-440 assemblies, we have simulated assemblies containing one or a few missing fuel rods or having only one emitting rod (the other rods being present but not emitting) in various well-chosen places, configurations that are not accessible in practice. The single-emitting rod results show in great detail those parts of the sinogram that contain most of the information for the particular rod position. How this information might be used for obtaining better images, especially of the central region of a fuel assembly, will be discussed
PGET Monte Carlo simulations using Serpent
Since 2017, over 100 spent nuclear fuel assemblies at the Finnish nuclear power plants have been imaged with the Passive Gamma Emission Tomography (PGET) device in preparation of the implementation of PGET in the safeguards infrastructure of the Finnish geological repository. In order to increase understanding of the PGET method and guide its further development, we have recently implemented PGET in Serpent, a widely-used neutron and photon transport Monte Carlo simulation code. We will discuss the major aspects of this implementation and illustrate the usefulness of the simulations with a few examples. The PGET device as used in the measurements (which was developed under the guidance of IAEA and is approved for safeguards inspections) was implemented in a very realistic way based on its technical drawings. The simulation produces sinograms in user-defined energy windows as well as the uncertainty on these sinograms. Tomographic images are then reconstructed using the exact same algorithm as used for the measured data. A dedicated variance reduction scheme was implemented, increasing the computational efficiency by about a factor of 30. The simulation of the PGET response at one angular measurement position for 1 billion primary photons takes a few hours on a single 40-core node. The 1-sigma uncertainty in the highest intensity sinogram pixels is about a few percent. Aiming at improving the imaging of VVER-440 assemblies, we have simulated assemblies containing one or a few missing fuel rods or having only one emitting rod (the other rods being present but not emitting) in various well-chosen places, configurations that are not accessible in practice. The single-emitting rod results show in great detail those parts of the sinogram that contain most of the information for the particular rod position. How this information might be used for obtaining better images, especially of the central region of a fuel assembly, will be discussed
Effect of interstitial carbon on the evolution of early-stage irradiation damage in equi-atomic FeMnNiCoCr high-entropy alloys
Owing to their excellent radiation tolerance, some of the high-entropy alloys (HEAs) are considered as potential candidates for structural materials in extreme conditions. In order to shed light on the early-stage irradiation damage in HEAs, we performed positron annihilation spectroscopy on hydrogen implanted equiatomic FeMnNiCoCr and interstitial carbon-containing FeMnNiCoCr HEAs. We reveal primary damage as monovacancies in low dose irradiated HEAs. The enhancement of Frenkel pair recombination by C addition is observed in C-containing HEAs. In addition, the C interstitials suppress the vacancy cluster formation in high dose irradiated HEAs.Peer reviewe
Quantum Monte Carlo study of Doppler broadening of positron annihilation radiation in semiconductors and insulators
The measurement of the momentum distribution of positron annihilation
radiation is a powerful method to detect and identify open-volume defects in
crystalline solids. The Doppler broadening of the 511 keV line of the
electron-positron annihilation event reflects the momentum density of
annihilating pairs and local electron momenta at positron annihilation sites.
It can provide information on the chemical surroundings of vacancies, such as
the impurity atoms around them. Accurate methods based on first-principles
calculations are crucial for interpreting measured Doppler spectra. In this
work we will validate such a method based on variational quantum Monte Carlo by
benchmarking results in aluminium nitride and silicon against experimental data
measured from defect-free reference samples. The method directly models
electron-positron correlations using variational wave functions. We achieve
better agreement with experiments for our test set than conventional
state-of-the-art methods. We show that normalized Doppler broadening spectra
calculated with quantum Monte Carlo converge rapidly as a function of
simulation cell size, and backflow transformations have only a minor effect.
This makes the method robust and practical to support positron-based
spectroscopies.Comment: 10 pages, 4 figure
Interfacial N Vacancies in GaN/(Al,Ga)N/GaN Heterostructures
We show that N-polar GaN/(Al, Ga)N/GaN heterostructures exhibit significant N deficiency at the bottom (Al, Ga)N/GaN interface, and that these N vacancies are responsible for the trapping of holes observed in unoptimized N-polar GaN/(Al, Ga)N/GaN high electron mobility transistors. We arrive at this conclusion by performing positron annihilation experiments on GaN/(Al, Ga)N/GaN heterostructures of both N and Ga polarity, as well as state-of-the-art theoretical calculations of the positron states and positron-electron annihilation signals. We suggest that the occurrence of high interfacial N vacancy concentrations is a universal property of nitride semiconductor heterostructures at net negative polarization interfaces.Peer reviewe