Engineering topics

Johnston, Steven

Publications and source records attributed to Johnston, Steven.

At least 19 records

DOE OSTI2025

Sunny.jl

Sunny is a Julia package designed to serve the needs of the quantum magnetism community. It supports the specification of a very broad class of spin models and a diverse suite of numerical solvers. These include powerful methods for simulating spin dynamics both in and out of equilibrium. Uniquely, it features a broad generalization of classical and semiclassical approaches to SU(N) coherent states, which is useful for studying systems exhibiting strong spin-orbit coupling or local entanglement effects. Sunny also offers a well-developed framework for calculating the dynamical spin structure factor, enabling direct comparison with scattering experiments. Ease of use is a priority, with tools for symmetry-guided modeling and interactive visualization.

Dahlbom, David↗

DOE OSTI2024

SmoQyDEAC.jl: A differential evolution package for the analytic continuation of imaginary time correlation functions

We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl package. It also provides detailed benchmarks of the package's output against the popular maximum entropy and stochastic analytic continuation methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Kekulé valence bond order in the honeycomb lattice optical Su-Schrieffer-Heeger model and its relevance to graphene

We perform sign-problem-free determinant quantum Monte Carlo simulations of the optical Su- Schrieffer-Heeger model on a half-filled honeycomb lattice. In particular, we investigate the model’s semi-metal (SM) to Kekulé Valence Bond Solid (KVBS) phase transition at zero and finite temper- atures as a function of phonon energy and interaction strength. Using hybrid Monte Carlo sampling methods we can simulate the model near the adiabatic regime, allowing us to access regions of parameter space relevant to graphene. Our simulations suggest that the SM-KVBS transition is weakly first-order at all temperatures, with graphene situated close to the phase boundary in the SM region of the phase diagram. Furthermore, our results highlight the important role bond-stretching phonon modes play in the formation of KVBS order in strained graphene-derived systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Data Files for P. Mai et al., "Fluctuating charge-density-wave correlations in the three-band Hubbard model" (2024)

The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model’s fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.

cuprates↗

DOE OSTI2024

Classical dynamics of the antiferromagnetic Heisenberg spin ladder

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the S = 1/2 ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule. Here, we anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly-coupled dimers, situations that remain challenging for alternative numerical methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

BEST Field Degradation Research

Data for control and 5-year fielded bifacial modules of PERC and SHJ technologies, including IV Curves, EQE, EL and PL, for degradation and performance studies.

14 SOLAR ENERGY↗

DOE OSTI2024

Magnetic, charge, and bond order in the two-dimensional Su-Schrieffer-Heeger-Holstein model

Most nonperturbative numerical studies of electron-phonon interactions focus on model Hamiltonians where the electrons interact with a phonon branch via a single type of microscopic mechanism. Two commonly explored couplings in this context are the Holstein and Su-Schrieffer-Heeger (SSH) interactions, which describe phonons modulating the on-site energy and intersite electron hopping, respectively. Many materials, however, have multiple phonon branches that can each interact with electronic degrees of freedom in different ways. We present here a determinant quantum Monte Carlo study of the half-filled two-dimensional (bond) SSH-Holstein Hamiltonian, where electrons couple to different phonon branches via either the Holstein or SSH mechanism. As a result, we map the model's phase diagram and determine the nature of the transitions between charge-density wave, bond-order wave, and antiferromagnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Suppressed superexchange interactions in the cuprates by bond-stretching oxygen phonons

We study a multiorbital Hubbard–Su-Schrieffer-Heeger model for the one-dimensional (1D) corner-shared cuprates in the adiabatic and nonadiabatic limits using exact diagonalization and determinant quantum Monte Carlo. At half filling and in the adiabatic limit, lattice dimerization can be achieved only over a narrow range of couplings slightly below a critical coupling g c . Beyond this value, the sign of the effective hopping changes, and the lattice becomes unstable. Strong lattice fluctuations replace the dimerization state in the nonadiabatic case. We also examine the model's temperature-dependent uniform and dynamical magnetic susceptibilities and compare them to the results of an effective spin-1/2 Heisenberg model. In doing so, we demonstrate that lattice fluctuations induced by the e-ph interaction suppress the effective superexchange interaction when g < g c . Furthermore, our results elucidate the effect of bond-stretching phonons in the parent cuprate compounds in general and are particularly relevant to 1D cuprates, where strong e-ph interactions have recently been inferred.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Comparative study of the superconductivity in the Holstein and optical Su-Schrieffer-Heeger models

Theoretical studies suggest that Su-Schrieffer-Heeger-like electron-phonon (e-ph) interactions can mediate high-temperature bipolaronic superconductivity that is robust against repulsive electron-electron interactions. Here we present a comparative analysis of the pairing and competing charge/bond correlations in the two-dimensional Holstein and optical Su-Schrieffer-Heeger (SSH) models using numerically exact determinant quantum Monte Carlo. We find that the SSH interactions support light bipolarons and strong superconducting correlations out to relatively large values of the e-ph coupling λ and densities near half-filling, while the Holstein interaction does not due to the formation of heavy bipolarons and competing charge-density-wave order. We further find that the Holstein and SSH models have comparable pairing correlations in the weak coupling limit for carrier concentrations < n > <<1, where competing orders and polaronic effects are absent. These results support the proposal that SSH (bi)polarons can support superconductivity to larger values of λ in comparison to the Holstein polaron, but that the resulting T c gains are small in the weak coupling limit. We also find that the SSH model's pairing correlations are suppressed after including a weak on-site Hubbard repulsion. Furthermore, these results have important implications for identifying and engineering bipolaronic superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Comparative determinant quantum Monte Carlo study of the acoustic and optical variants of the Su-Schrieffer-Heeger model

We compare the acoustic Su–Schrieffer–Heeger (SSH) model with two of its optical variants where the phonons are defined on either on the sites or bonds of the system. First, we discuss how to make fair comparisons between these models in any dimension by ensuring their dimensionless coupling λ and relevant phonon energies are the same. We then use determinant quantum Monte Carlo to perform nonperturbative and sign-problem-free simulations of all three models on one-dimensional chains at and away from half-filling. By comparing the results obtained from each model, we demonstrate that the optical and acoustic models produce near identical results within error bars for suitably chosen phonon energies and λ at half-filling. In contrast, the bond model has quantitatively different behavior due to its coupling to the q = 0 phonon mode. These differences also manifest in the total length of the chain, which shrinks for the bond model but not for the acoustic and optical models when λ ≠ 0. Furthermore, our results have important implications for quantum Monte Carlo modeling of SSH-like interactions, where these models are sometimes regarded as being interchangeable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Optoelectronic and structural characterization of trapezoidal defects in 4H-SiC epilayers and the effect on MOSFET reliability

To this day, trapezoidal defects are found in clusters and high counts in wafers representing the industry standard in terms of material quality being produced. This study sheds light on the nature, origin, behavior, and impact of this defect on device yield and reliability. Trapezoidal defects in 4H-SiC epitaxial layers were investigated by photoluminescence (PL) imaging, scanning electron microscopy (SEM), cathodoluminescence spectrum imaging (CLSI), SEM electron beam induced current (EBIC) imaging, and by transmission electron microscopy (TEM) observation. The bar-shaped stacking faults were identified by the PL and CL measurements with a peak emission wavelength of 420 and 450 nm. An optoelectronic behavioral study based on the recombination enhanced dislocation glide mechanism revealed how expanding dislocations and stacking faults interact with each other. Combining the luminescence and microscopy results, the nature of the stacking faults was identified as being a combination of Shockley-type and Frank-type stacking faults. The TEM analysis showed that these defects originate from the substrate and the stacking sequences of some of the faults were determined as (…2, 4, 2…) and (…2, 3, 2…) in the Zhdanov's notation by high-resolution TEM. The origin of this defect is speculated based on our results and previous reports. The EBIC imaging showed that the high density of SFs in these towers is a strong site of carrier recombination, which presumably has an impact on the transfer characteristics of SiC devices. Furthermore, these defects have shown to impact metal oxide semiconductor field effect transistors electrical performance via an increase in the on-state resistance depending on the coverage percentage of the tower of defects in the active area of the device.

4H-SiC↗

DOE OSTI2023

Stacking Faults Originating from Star-Defects in 4H-SiC

Intense efforts are currently in progress to study various sources of basal plane dislocations (BPDs) in SiC epitaxial layers. BPDs can generate Shockley-type stacking faults (SSFs) in SiC epitaxial layers, which have been shown to be associated with the degradation of power devices. This study shows that the star-shaped defect can be a source of several BPDs in the epitaxial layer. We investigate the complex microstructure of the star defect, the generation of BPDs, and expansion of SSFs using various complementary microscopy and optical techniques. We show direct evidence that star-defects can be a nucleation point of single-SSFs that can expand at the core of the defect. Newly found secondary dislocation arrays extending over a few centimeters away are found to be emanating from the primary arms of the star defect. The presence of such dislocation walls and the expansion of single-SSFs will affect the yield of numerous die on a wafer. Further understanding of the formation mechanism of stacking faults generated from star-defects as provided in this study helps understand their effect on SiC-based devices, which is crucial to assess device reliability.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Robust charge-density-wave correlations in the electron-doped single-band Hubbard model

There is growing evidence that the hole-doped single-band Hubbard and t–J models do not have a superconducting ground state reflective of the high-temperature cuprate superconductors but instead have striped spin- and charge-ordered ground states. Nevertheless, it is proposed that these models may still provide an effective low-energy model for electron-doped materials. Here we study the finite temperature spin and charge correlations in the electron-doped Hubbard model using quantum Monte Carlo dynamical cluster approximation calculations and contrast their behavior with those found on the hole-doped side of the phase diagram. We find evidence for a charge modulation with both checkerboard and unidirectional components decoupled from any spin-density modulations. These correlations are inconsistent with a weak-coupling description based on Fermi surface nesting, and their doping dependence agrees qualitatively with resonant inelastic x-ray scattering measurements. Our results provide evidence that the single-band Hubbard model describes the electron-doped cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Charge order in the kagome lattice Holstein model: a hybrid Monte Carlo study

The Holstein model is a paradigmatic description of the electron-phonon interaction, in which electrons couple to local dispersionless phonon modes, independent of momentum. The model has been shown to host a variety of ordered ground states such as charge density wave (CDW) order and superconductivity on several geometries, including the square, honeycomb, and Lieb lattices. In this work, we study CDW formation in the Holstein model on the kagome lattice, using a recently developed hybrid Monte Carlo simulation method. We present evidence for √3 × √3 CDW order at an average electron filling of < n > = 2/3 per site, with an ordering wavevector at the K-points of the Brillouin zone. We estimate a phase transition occurring at T c ≈ t/18, where t is the nearest-neighbor hopping parameter. Our simulations find no signature of CDW order at other electron fillings or ordering momenta for temperatures T ≥ t/20.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Group-equivariant autoencoder for identifying spontaneously broken symmetries

We introduce the group-equivariant autoencoder (GE autoencoder), a deep neural network (DNN) method that locates phase boundaries by determining which symmetries of the Hamiltonian have spontaneously broken at each temperature. We use group theory to deduce which symmetries of the system remain intact in all phases, and then use this information to constrain the parameters of the GE autoencoder such that the encoder learns an order parameter invariant to these “never-broken” symmetries. This procedure produces a dramatic reduction in the number of free parameters such that the GE-autoencoder size is independent of the system size. We include symmetry regularization terms in the loss function of the GE autoencoder so that the learned order parameter is also equivariant to the remaining symmetries of the system. By examining the group representation by which the learned order parameter transforms, we are then able to extract information about the associated spontaneous symmetry breaking. We test the GE autoencoder on the 2D classical ferromagnetic and antiferromagnetic Ising models, finding that the GE autoencoder (1) accurately determines which symmetries have spontaneously broken at each temperature; (2) estimates the critical temperature in the thermodynamic limit with greater accuracy, robustness, and time efficiency than a symmetry-agnostic baseline autoencoder; and (3) detects the presence of an external symmetry-breaking magnetic field with greater sensitivity than the baseline method. Lastly, we describe various key implementation details, including a quadratic-programming-based method for extracting the critical temperature estimate from trained autoencoders and calculations of the DNN initialization and learning rate settings required for fair model comparisons.

42 ENGINEERING↗

DOE OSTI2023

A hybrid Monte Carlo study of bond-stretching electron–phonon interactions and charge order in BaBiO3

Abstract The relationship between electron–phonon ( e -ph) interactions and charge-density-wave (CDW) order in the bismuthate family of high-temperature superconductors remains unresolved. We address this question using nonperturbative hybrid Monte Carlo calculations for the parent compound BaBiO 3 . Our model includes the Bi 6 s and O 2 p σ orbitals and coupling to the Bi-O bond-stretching branch of optical phonons via modulations of the Bi-O hopping integral. We simulate three-dimensional clusters of up to 4000 orbitals, with input model parameters taken from ab initio electronic structure calculations and a phonon energy ℏΩ 0 = 60 meV. Our results demonstrate that the coupling to the bond-stretching modes is sufficient to reproduce the CDW transition in this system, despite a relatively small dimensionless coupling. We also find that the transition deviates from the weak-coupling Peierls’ picture. This work demonstrates that off-diagonal e -ph interactions in orbital space are vital in establishing the bismuthate phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Charge correlations suppress unconventional pairing in the Holstein model

In a recent work by Schrodi et al. [Phys. Rev. B. 104, L140506 (2021)], the authors find an unconventional superconducting state with a sign-changing order parameter using the Migdal-Eliashberg theory, including the first vertex correction. This unconventional solution arises despite using an isotropic bare electron-phonon coupling in the Hamiltonian. We examine this claim using hybrid quantum Monte Carlo for a single-band Holstein model with a cuprate-like noninteracting band structure and identical parameters to Schrodi et al.. Our Monte Carlo results for these parameters suggest that unconventional pairing correlations do not exceed their noninteracting values at any carrier concentration we have checked. Instead, strong charge-density-wave correlations persist at the lowest accessible temperatures for dilute and nearly half-filled bands. As a result, we present arguments for how vertex-corrected Migdal-Eliashberg calculation schemes can lead to uncontrolled results in the presence of Fermi surface nesting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Fluctuating intertwined stripes in the strange metal regime of the Hubbard model

Strongly correlated electron systems host a variety of poorly understood correlations in their high-temperature normal state. Unlike ordered phases defined by order parameters, regions of the normal state are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground-state orders. Here, in this paper, we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short ranged; yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗