Source author record

Andrew D. Baczewski

Andrew D. Baczewski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

First-principles simulation of light-ion microscopy of graphene

The extreme sensitivity of 2D materials to defects and nanostructure requires precise imaging techniques to verify presence of desirable and absence of undesirable features in the atomic geometry. Helium-ion beams have emerged as a promising materials imaging tool, achieving up to 20 times higher resolution and 10 times larger depth-of-field than conventional or environmental scanning electron microscopes. Here, we offer first-principles theoretical insights to advance ion-beam imaging of atomically thin materials by performing real-time time-dependent density functional theory simulations of single impacts of 10-200 keV light ions in free-standing graphene. We predict that detecting electrons emitted from the back of the material (the side from which the ion exits) would result in up to 3 times higher signal and up to 5 times higher contrast images, making 2D materials especially compelling targets for ion-beam microscopy. We also find that the charge induced in the graphene equilibrates on a sub-fs time scale, leading to only slight disturbances in the carbon lattice that are unlikely to damage the atomic structure for any of the beam parameters investigated here.

preprint2022arXiv

Precision tomography of a three-qubit donor quantum processor in silicon

Nuclear spins were among the first physical platforms to be considered for quantum information processing, because of their exceptional quantum coherence and atomic-scale footprint. However, their full potential for quantum computing has not yet been realized, due to the lack of methods to link nuclear qubits within a scalable device combined with multi-qubit operations with sufficient fidelity to sustain fault-tolerant quantum computation. Here we demonstrate universal quantum logic operations using a pair of ion-implanted 31P donor nuclei in a silicon nanoelectronic device. A nuclear two-qubit controlled-Z gate is obtained by imparting a geometric phase to a shared electron spin, and used to prepare entangled Bell states with fidelities up to 94.2(2.7)%. The quantum operations are precisely characterised using gate set tomography (GST), yielding one-qubit average gate fidelities up to 99.95(2)%, two-qubit average gate fidelity of 99.37(11)% and two-qubit preparation/measurement fidelities of 98.95(4)%. These three metrics indicate that nuclear spins in silicon are approaching the performance demanded in fault-tolerant quantum processors. We then demonstrate entanglement between the two nuclei and the shared electron by producing a Greenberger-Horne-Zeilinger three-qubit state with 92.5(1.0)% fidelity. Since electron spin qubits in semiconductors can be further coupled to other electrons or physically shuttled across different locations, these results establish a viable route for scalable quantum information processing using donor nuclear and electron spins.

preprint2022arXiv

Reaction pathways of BCl$_3$ for acceptor delta-doping of silicon

BCl$_3$ is a promising candidate for atomic-precision acceptor doping in Si, but optimizing the electrical properties of structures created with this technique requires a detailed understanding of adsorption and dissociation pathways for this precursor. Here, we use density functional theory and scanning tunneling microscopy (STM) to identify and explore these pathways for BCl$_3$ on Si(100) at different annealing temperatures. We demonstrate that BCl$_3$ adsorbs selectively without a reaction barrier, and subsequently dissociates relatively easily with reaction barriers $\approx$1 eV. Using this dissociation pathway, we parameterize a Kinetic Monte Carlo model to predict B incorporation rates as a function of dosing conditions. STM is used to image BCl$_{3}$ adsorbates, identifying several surface configurations and tracking the change in their distribution as a function of the annealing temperature, matching predictions of the kinetic model well. This straightforward pathway for atomic-precision acceptor doping helps enable a wide range of applications including bipolar nanoelectronics, acceptor-based qubits, and superconducting Si.

preprint2022arXiv

Robust incorporation in multi-donor patches created using atomic-precision advanced manufacturing

Atomic-precision advanced manufacturing enables the placement of dopant atoms within $\pm$1 lattice site in crystalline Si. However, it has recently been shown that reaction kinetics can introduce uncertainty in whether a single donor will incorporate at all in a minimal 3-dimer lithographic window. In this work, we explore the combined impact of lithographic variation and stochastic kinetics on P incorporation as the size of such a window is increased. We augment a kinetic model for PH$_3$ dissociation leading to P incorporation on Si(100)-2$\times$1 to include barriers for reactions across distinct dimer rows. Using this model, we demonstrate that even for a window consisting of 2$\times$3 silicon dimers, the probability that at least one donor incorporates is nearly unity. We also examine the impact of size of the lithographic window, finding that the incorporation fraction saturates to $δ$-layer like coverage as the circumference-to-area ratio approaches zero. We predict that this incorporation fraction depends strongly on the dosage of the precursor, and that the standard deviation of the number of incorporations scales as $\sim \sqrt{n}$, as would be expected for a series of largely independent incorporation events. Finally, we characterize an array of experimentally prepared multi-donor lithographic windows and use our kinetic model to study variability due to the observed lithographic roughness, predicting a negligible impact on incorporation statistics. We find good agreement between our model and the inferred incorporation in these windows from scanning tunneling microscope measurements, indicating the robustness of atomic-precision advanced manufacturing to errors in patterning for multi-donor patches.

preprint2021arXiv

Hole in one: Pathways to deterministic single-acceptor incorporation in Si(100)-2$\times$1

Stochastic incorporation kinetics can be a limiting factor in the scalability of semiconductor fabrication technologies using atomic-precision techniques. While these technologies have recently been extended from donors to acceptors, the extent to which kinetics will impact single-acceptor incorporation has yet to be assessed. We develop and apply an atomistic model for the single-acceptor incorporation rates of several recently demonstrated precursor molecules: diborane (B$_2$H$_6$), boron trichloride (BCl$_3$), and aluminum trichloride in both monomer (AlCl$_3$) and dimer forms (Al$_2$Cl$_6$), to identify the acceptor precursor and dosing conditions most likely to yield deterministic incorporation. While all three precursors can achieve single-acceptor incorporation, we predict that diborane is unlikely to achieve deterministic incorporation, boron trichloride can achieve deterministic incorporation with modest heating (50 $^{\circ}$C), and aluminum trichloride can achieve deterministic incorporation at room temperature. We conclude that both boron and aluminum trichloride are promising precursors for atomic-precision single-acceptor applications, with the potential to enable the reliable production of large arrays of single-atom quantum devices.

preprint2021arXiv

The impact of stochastic incorporation on atomic-precision Si:P arrays

Scanning tunneling microscope lithography can be used to create nanoelectronic devices in which dopant atoms are precisely positioned in a Si lattice within $\sim$1 nm of a target position. This exquisite precision is promising for realizing various quantum technologies. However, a potentially impactful form of disorder is due to incorporation kinetics, in which the number of P atoms that incorporate into a single lithographic window is manifestly uncertain. We present experimental results indicating that the likelihood of incorporating into an ideally written three-dimer single-donor window is $63 \pm 10\%$ for room-temperature dosing, and corroborate these results with a model for the incorporation kinetics. Nevertheless, further analysis of this model suggests conditions that might raise the incorporation rate to near-deterministic levels. We simulate bias spectroscopy on a chain of comparable dimensions to the array in our yield study, indicating that such an experiment may help confirm the inferred incorporation rate.

preprint2019arXiv

Coherent electrical control of a single high-spin nucleus in silicon

Nuclear spins are highly coherent quantum objects. In large ensembles, their control and detection via magnetic resonance is widely exploited, e.g. in chemistry, medicine, materials science and mining. Nuclear spins also featured in early ideas and demonstrations of quantum information processing. Scaling up these ideas requires controlling individual nuclei, which can be detected when coupled to an electron. However, the need to address the nuclei via oscillating magnetic fields complicates their integration in multi-spin nanoscale devices, because the field cannot be localized or screened. Control via electric fields would resolve this problem, but previous methods relied upon transducing electric signals into magnetic fields via the electron-nuclear hyperfine interaction, which severely affects the nuclear coherence. Here we demonstrate the coherent quantum control of a single antimony (spin-7/2) nucleus, using localized electric fields produced within a silicon nanoelectronic device. The method exploits an idea first proposed in 1961 but never realized experimentally with a single nucleus. Our results are quantitatively supported by a microscopic theoretical model that reveals how the purely electrical modulation of the nuclear electric quadrupole interaction, in the presence of lattice strain, results in coherent nuclear spin transitions. The spin dephasing time, 0.1 seconds, surpasses by orders of magnitude those obtained via methods that require a coupled electron spin for electrical drive. These results show that high-spin quadrupolar nuclei could be deployed as chaotic models, strain sensors and hybrid spin-mechanical quantum systems using all-electrical controls. Integrating electrically controllable nuclei with quantum dots could pave the way to scalable nuclear- and electron-spin-based quantum computers in silicon that operate without the need for oscillating magnetic fields.

preprint2016arXiv

X-Ray Thomson scattering without the Chihara decomposition

X-Ray Thomson Scattering (XRTS) is an important experimental technique used to measure the temperature, ionization state, structure, and density of warm dense matter (WDM). The fundamental property probed in these experiments is the electronic dynamic structure factor (DSF). In most models, this is decomposed into three terms [Chihara, J. Phys. F: Metal Phys. {\bf 17}, 295 (1987)] representing the response of tightly bound, loosely bound, and free electrons. Accompanying this decomposition is the classification of electrons as either bound or free, which is useful for gapped and cold systems but becomes increasingly questionable as temperatures and pressures increase into the WDM regime. In this work we provide unambiguous first principles calculations of the dynamic structure factor of warm dense beryllium, independent of the Chihara form, by treating bound and free states under a single formalism. The computational approach is real-time finite-temperature time-dependent density functional theory (TDDFT) being applied here for the first time to WDM. We compare results from TDDFT to Chihara-based calculations for experimentally relevant conditions in shock-compressed beryllium.

preprint2015arXiv

Multivalley effective mass theory simulation of donors in silicon

Last year, Salfi et al. made the first direct measurements of a donor wave function and found extremely good theoretical agreement with atomistic tight-binding [Salfi et al., Nat. Mater. 13, 605 (2014)]. Here, we show that multi-valley effective mass theory, applied properly, does achieve close agreement with tight-binding and hence gives reliable predictions. To demonstrate this, we variationally solve the coupled six-valley Shindo-Nara equations, including silicon's full Bloch functions. Surprisingly, we find that including the full Bloch functions necessitates a tetrahedral, rather than spherical, donor central cell correction to accurately reproduce the experimental energy spectrum of a phosphorus impurity in silicon. We cross-validate this method against atomistic tight-binding calculations, showing that the two theories agree well for the calculation of donor-donor tunnel coupling. Further, we benchmark our results by performing a statistical uncertainty analysis, confirming that derived quantities such as the wave function profile and tunnel couplings are robust with respect to variational energy fluctuations. Finally, we apply this method to exhaustively enumerate the tunnel coupling for all donor-donor configurations within a large search volume, demonstrating conclusively that the tunnel coupling has no spatially stable regions. Though this instability is problematic for reliably coupling donor pairs for two-qubit operations, we identify specific target locations where donor qubits can be placed with scanning tunneling microscopy technology to achieve reliably large tunnel couplings.

preprint2015arXiv

The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus

An outstanding challenge of theoretical electronic structure is the description of van der Waals (vdW) interactions in molecules and solids. Renewed interest in resolving this is in part motivated by the technological promise of layered systems including graphite, transition metal dichalcogenides, and more recently, black phosphorus, in which the interlayer interaction is widely believed to be dominated by these types of forces. We report a series of quantum Monte Carlo (QMC) calculations for bulk black phosphorus and related few-layer phosphorene, which elucidate the nature of the forces that bind these systems and provide benchmark data for the energetics of these systems. We find a significant charge redistribution due to the interaction between electrons on adjacent layers. Comparison to density functional theory (DFT) calculations indicate not only wide variability even among different vdW corrected functionals, but the failure of these functionals to capture the trend of reorganization predicted by QMC. The delicate interplay of steric and dispersive forces between layers indicate that few-layer phosphorene presents an unexpected challenge for the development of vdW corrected DFT.

preprint2014arXiv

A Discontinuous Galerkin Time Domain Framework for Periodic Structures Subject To Oblique Excitation

A nodal Discontinuous Galerkin (DG) method is derived for the analysis of time-domain (TD) scattering from doubly periodic PEC/dielectric structures under oblique interrogation. Field transformations are employed to elaborate a formalism that is free from any issues with causality that are common when applying spatial periodic boundary conditions simultaneously with incident fields at arbitrary angles of incidence. An upwind numerical flux is derived for the transformed variables, which retains the same form as it does in the original Maxwell problem for domains without explicitly imposed periodicity. This, in conjunction with the amenability of the DG framework to non-conformal meshes, provides a natural means of accurately solving the first order TD Maxwell equations for a number of periodic systems of engineering interest. Results are presented that substantiate the accuracy and utility of our method.

preprint2014arXiv

Comment on "Self-Averaging Stochastic Kohn-Sham Density-Functional Theory"

In a recent Letter, Baer et al. present a stochastic method for Kohn-Sham density functional theory calculations. Their convergence criterion is the self-averaging total energy per electron, which requires a number of statistical samples that decreases with system size and enables a sublinear-scaling computational cost. However, the electron density, atomic forces, orbital energies, and many other physical quantities do not self-average, and this procedure causes their statistical errors to grow with system size. Convergence of non-self-averaging quantities requires that the number of samples be maintained with system size and is incompatible with a sublinear-scaling cost.

preprint2013arXiv

Numerical Integration of the Extended Variable Generalized Langevin Equation with a Positive Prony Representable Memory Kernel

Generalized Langevin dynamics (GLD) arise in the modeling of a number of systems, ranging from structured fluids that exhibit a viscoelastic mechanical response, to biological systems, and other media that exhibit anomalous diffusive phenomena. Molecular dynamics (MD) simulations that include GLD in conjunction with external and/or pairwise forces require the development of numerical integrators that are efficient, stable, and have known convergence properties. In this article, we derive a family of extended variable integrators for the Generalized Langevin equation (GLE) with a positive Prony series memory kernel. Using stability and error analysis, we identify a superlative choice of parameters and implement the corresponding numerical algorithm in the LAMMPS MD software package. Salient features of the algorithm include exact conservation of the first and second moments of the equilibrium velocity distribution in some important cases, stable behavior in the limit of conventional Langevin dynamics, and the use of a convolution-free formalism that obviates the need for explicit storage of the time history of particle velocities. Capability is demonstrated with respect to accuracy in numerous canonical examples, stability in certain limits, and an exemplary application in which the effect of a harmonic confining potential is mapped onto a memory kernel.

preprint2011arXiv

An O(N) Method for Rapidly Computing Periodic Potentials Using Accelerated Cartesian Expansions

The evaluation of long-range potentials in periodic, many-body systems arises as a necessary step in the numerical modeling of a multitude of interesting physical problems. Direct evaluation of these potentials requires O(N^2) operations and O(N^2) storage, where N is the number of interacting bodies. In this work, we present a method, which requires O(N) operations and O(N) storage, for the evaluation of periodic Helmholtz, Coulomb, and Yukawa potentials with periodicity in 1-, 2-, and 3-dimensions, using the method of Accelerated Cartesian Expansions (ACE). We present all aspects necessary to effect this acceleration within the framework of ACE including the necessary translation operators, and appropriately modifying the hierarchical computational algorithm. We also present several results that validate the efficacy of this method with respect to both error convergence and cost scaling, and derive error bounds for one exemplary potential.

preprint2011arXiv

The Rapid Analysis of Scattering from Periodic Dielectric Structures Using Accelerated Cartesian Expansions (ACE)

The analysis of fields in periodic dielectric structures arise in numerous applications of recent interest, ranging from photonic bandgap (PBG) structures and plasmonically active nanostructures to metamaterials. To achieve an accurate representation of the fields in these structures using numerical methods, dense spatial discretization is required. This, in turn, affects the cost of analysis, particularly for integral equation based methods, for which traditional iterative methods require O(N^2) operations, N being the number of spatial degrees of freedom. In this paper, we introduce a method for the rapid solution of volumetric electric field integral equations used in the analysis of doubly periodic dielectric structures. The crux of our method is the ACE algorithm, which is used to evaluate the requisite potentials in O(N) cost. Results are provided that corroborate our claims of acceleration without compromising accuracy, as well as the application of our method to a number of compelling photonics applications.

preprint2007arXiv

Electronic and Structural Characteristics of Zinc-Blende Wurtzite Biphasic Homostructure GaN Nanowires

We report a new biphasic crystalline wurtzite/zinc-blende homostructure in gallium nitride nanowires. Cathodoluminescence was used to quantitatively measure the wurtzite and zinc-blende band gaps. High resolution transmission electron microscopy was used to identify distinct wurtzite and zinc-blende crystalline phases within single nanowires through the use of selected area electron diffraction, electron dispersive spectroscopy, electron energy loss spectroscopy, and fast Fourier transform techniques. A mechanism for growth is identified.