Showing 1–50 of 205 results for author: Smith, J

Tunable Second-Order Structural Transition in As-Deficient MnAs

Authors: B. D. White, K. Huang, I. L. Fipps, J. J. Hamlin, S. Jang, G. J. Smith, B. Xia, J. W. Simonson, C. S. Nelson, M. C. Aronson, M. B. Maple

Abstract: We report measurements of magnetization, specific heat, and thermal expansion performed on As-deficient MnAs single crystals (MnAs$_{0.968}$). Ferromagnetic order is observed near $T_C \simeq$ 306 K on warming and $T_C \simeq$ 302 K on cooling, which is consistent with previously-reported values for stoichiometric MnAs samples. In contrast, the second-order structural phase transition is observed… ▽ More We report measurements of magnetization, specific heat, and thermal expansion performed on As-deficient MnAs single crystals (MnAs$_{0.968}$). Ferromagnetic order is observed near $T_C \simeq$ 306 K on warming and $T_C \simeq$ 302 K on cooling, which is consistent with previously-reported values for stoichiometric MnAs samples. In contrast, the second-order structural phase transition is observed at $T_S \simeq$ 353 K, which is nearly 50 K lower than in the stoichiometric compound. We observe differences in the thermal expansion of our samples when compared to reports of stoichiometric MnAs including: (1) the $\sim$1.5% volume decrease at $T_C$ is smaller than the expected value of 1.9%, (2) the lattice parameters perpendicular to the basal plane exhibit a discontinuous jump of $\sim$1.1% at $T_C$ instead of being continuous across $T_C$, and (3) thermal expansion perpendicular to the basal plane for $T_C \le T \le$ 315 K is negative rather than positive. We also observe a correlation between the ratio of hexagonal lattice parameters, $c/a$, and $T_S$, strongly suggesting that the degree of structural anisotropy in MnAs could play an important role in tuning $T_S$. △ Less

Submitted 14 August, 2024; originally announced August 2024.

Predicting quantum materials properties using novel faithful machine learning embeddings

Authors: Gavin Nop, Micah Mundy, Durga Paudyal, Jonathan Smith

Abstract: Machine Learning (ML) is accelerating the progress of materials prediction and classification, with particular success in CGNN designs. While classical ML methods remain accessible, advanced deep networks are still challenging to build and train. We introduce two new adaptations and refine two existing ML networks for generic crystalline quantum materials properties prediction and optimization. Th… ▽ More Machine Learning (ML) is accelerating the progress of materials prediction and classification, with particular success in CGNN designs. While classical ML methods remain accessible, advanced deep networks are still challenging to build and train. We introduce two new adaptations and refine two existing ML networks for generic crystalline quantum materials properties prediction and optimization. These new models achieve state-of-the-art performance in predicting TQC classification and strong performance in predicting band gaps, magnetic classifications, formation energies, and symmetry group. All networks easily generalize to all quantum crystalline materials property predictions. To support this, full implementations and automated methods for data handling and materials predictions are provided, facilitating the use of deep ML methods in quantum materials science. Finally, dataset error rates are analyzed using an ensemble model to identify and highlight highly atypical materials for further investigations. △ Less

Submitted 25 July, 2024; originally announced July 2024.

A green solvent system for precursor phase-engineered sequential deposition of stable formamidinium lead triiodide for perovskite solar cells

Authors: Benjamin M. Gallant, Philippe Holzhey, Joel A. Smith, Saqlain Choudhary, Karim A. Elmestekawy, Pietro Caprioglio, Igal Levine, Alex Sheader, Fengning Yang, Daniel T. W. Toolan, Rachel C. Kilbride, Augustin K. A. Zaininger, James M. Ball, M. Greyson Christoforo, Nakita Noel, Laura M. Herz, Dominik J. Kubicki, Henry J. Snaith

Abstract: Perovskite solar cells (PSCs) offer an efficient, inexpensive alternative to current photovoltaic technologies, with the potential for manufacture via high-throughput coating methods. However, challenges for commercial-scale solution-processing of metal-halide perovskites include the use of harmful solvents, the expense of maintaining controlled atmospheric conditions, and the inherent instabiliti… ▽ More Perovskite solar cells (PSCs) offer an efficient, inexpensive alternative to current photovoltaic technologies, with the potential for manufacture via high-throughput coating methods. However, challenges for commercial-scale solution-processing of metal-halide perovskites include the use of harmful solvents, the expense of maintaining controlled atmospheric conditions, and the inherent instabilities of PSCs under operation. Here, we address these challenges by introducing a high volatility, low toxicity, biorenewable solvent system to fabricate a range of 2D perovskites, which highly effective precursor phases for subsequent transformation to alpha-formamidinium lead triiodide (FAPbI3), fully processed under ambient conditions. PSCs utilising our FAPbI3 reproducibly show remarkable stability under illumination and elevated temperature (ISOS-L-2) and "damp heat" (ISOS-D-3) stressing, surpassing other state-of-the-art perovskite compositions. We determine that this enhancement is a consequence of the 2D precursor phase crystallisation route, which simultaneously avoids retention of residual low-volatility solvents (such as DMF and DMSO) and reduces the rate of degradation of FA+ in the material. Our findings highlight both the critical role of the initial crystallisation process in determining the operational stability of perovskite materials, and that neat FA+-based perovskites can be competitively stable despite the inherent metastability of the alpha-phase. △ Less

Submitted 14 June, 2024; v1 submitted 12 June, 2024; originally announced June 2024.

Materials for Quantum Technologies: a Roadmap for Spin and Topology

Authors: N. Banerjee, C. Bell, C. Ciccarelli, T. Hesjedal, F. Johnson, H. Kurebayashi, T. A. Moore, C. Moutafis, H. L. Stern, I. J. Vera-Marun, J. Wade, C. Barton, M. R. Connolly, N. J. Curson, K. Fallon, A. J. Fisher, D. A. Gangloff, W. Griggs, E. Linfield, C. H. Marrows, A. Rossi, F. Schindler, J. Smith, T. Thomson, O. Kazakova

Abstract: This Roadmap provides an overview of the critical role of materials in exploiting spin and topology for next-generation quantum technologies including computing, sensing, information storage and networking devices. We explore the key materials systems that support spin and topological phenomena and discuss their figures of merit. Spin and topology-based quantum technologies have several advantages… ▽ More This Roadmap provides an overview of the critical role of materials in exploiting spin and topology for next-generation quantum technologies including computing, sensing, information storage and networking devices. We explore the key materials systems that support spin and topological phenomena and discuss their figures of merit. Spin and topology-based quantum technologies have several advantages over their classical, charged-based counterparts, including non-volatility, faster data processing speeds, higher integration densities and lower power consumption. We discuss the main challenges facing the field, identify strategies to overcome them, and provide a realistic outlook on future possibilities of spin-based and topological materials in quantum technology applications. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Roadmap of the UKRI EPSRC Materials for Quantum Network (M4QN) Spin & Topology group. 25 pages, 5 figures, includes supplement. Comments welcome!

Formation of Beta-Indium Selenide Layers Grown via Selenium Passivation of InP(111)B Substrate

Authors: Kaushini S. Wickramasinghe, Candice Forrester, Martha R. McCartney, David J. Smith, Maria C. Tamargo

Abstract: Indium selenide, In2Se3, has recently attracted growing interest due to its novel properties, including room temperature ferroelectricity, outstanding photoresponsivity, and exotic in-plane ferroelectricity, which open up new regimes for next generation electronics. In2Se3 also provides the important advantage of tuning the electrical properties of ultra-thin layers with an external electrical and… ▽ More Indium selenide, In2Se3, has recently attracted growing interest due to its novel properties, including room temperature ferroelectricity, outstanding photoresponsivity, and exotic in-plane ferroelectricity, which open up new regimes for next generation electronics. In2Se3 also provides the important advantage of tuning the electrical properties of ultra-thin layers with an external electrical and magnetic field, making it a potential platform to study novel two-dimensional physics. Yet, In2Se3 has many different polymorphs, and it has been challenging to synthesize single-phase material, especially using scalable growth methods, as needed for technological applications. In this paper, we use aberration-corrected scanning transmission electron microscopy to characterize the microstructure of twin-free single-phase ultra-thin layers of beta-In2Se3, prepared by a unique molecular beam epitaxy approach. We emphasize features of the In2Se3 layer and In2Se3/InP interface which provide evidence for understanding the growth mechanism of the single-phase In2Se3. This novel approach for forming high-quality twin-free single phase two-dimensional crystals on InP substrates is likely to be applicable to other technologically important substrates. △ Less

Submitted 15 May, 2024; v1 submitted 15 May, 2024; originally announced May 2024.

Comments: 19 pages, 4 figures

MatterSim: A Deep Learning Atomistic Model Across Elements, Temperatures and Pressures

Authors: Han Yang, Chenxi Hu, Yichi Zhou, Xixian Liu, Yu Shi, Jielan Li, Guanzhi Li, Zekun Chen, Shuizhou Chen, Claudio Zeni, Matthew Horton, Robert Pinsler, Andrew Fowler, Daniel Zügner, Tian Xie, Jake Smith, Lixin Sun, Qian Wang, Lingyu Kong, Chang Liu, Hongxia Hao, Ziheng Lu

Abstract: Accurate and fast prediction of materials properties is central to the digital transformation of materials design. However, the vast design space and diverse operating conditions pose significant challenges for accurately modeling arbitrary material candidates and forecasting their properties. We present MatterSim, a deep learning model actively learned from large-scale first-principles computatio… ▽ More Accurate and fast prediction of materials properties is central to the digital transformation of materials design. However, the vast design space and diverse operating conditions pose significant challenges for accurately modeling arbitrary material candidates and forecasting their properties. We present MatterSim, a deep learning model actively learned from large-scale first-principles computations, for efficient atomistic simulations at first-principles level and accurate prediction of broad material properties across the periodic table, spanning temperatures from 0 to 5000 K and pressures up to 1000 GPa. Out-of-the-box, the model serves as a machine learning force field, and shows remarkable capabilities not only in predicting ground-state material structures and energetics, but also in simulating their behavior under realistic temperatures and pressures, signifying an up to ten-fold enhancement in precision compared to the prior best-in-class. This enables MatterSim to compute materials' lattice dynamics, mechanical and thermodynamic properties, and beyond, to an accuracy comparable with first-principles methods. Specifically, MatterSim predicts Gibbs free energies for a wide range of inorganic solids with near-first-principles accuracy and achieves a 15 meV/atom resolution for temperatures up to 1000K compared with experiments. This opens an opportunity to predict experimental phase diagrams of materials at minimal computational cost. Moreover, MatterSim also serves as a platform for continuous learning and customization by integrating domain-specific data. The model can be fine-tuned for atomistic simulations at a desired level of theory or for direct structure-to-property predictions, achieving high data efficiency with a reduction in data requirements by up to 97%. △ Less

Submitted 10 May, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Low-symmetry polymorph of GaP upends bonding paradigms of metallic high-pressure III-V compounds

Authors: Barbara Lavina, Enrique Zanardi, Andrés Mujica, Hyunchae Cynn, Yue Meng, Vitali Prakapenka, Jesse S. Smith

Abstract: The pressure-induced polymorphism of binary octect compounds has long been considered a settled problem although the possible atomic disordering of some phases remains a puzzling observation. Taking GaP as a case study, we conclude, through x-ray microdiffraction and first-principles calculations, that its high-pressure phase II (previously reported as being disordered) adopts in fact an ordered b… ▽ More The pressure-induced polymorphism of binary octect compounds has long been considered a settled problem although the possible atomic disordering of some phases remains a puzzling observation. Taking GaP as a case study, we conclude, through x-ray microdiffraction and first-principles calculations, that its high-pressure phase II (previously reported as being disordered) adopts in fact an ordered base-centered monoclinic structure previously unknown in this class of compounds. The formation of layered patterns with variable degrees of interlayer dimerization, as observed in GaP, marks a paradigm shift of our understanding of ordering in octect high-pressure phases which calls for a more extensive re-examination. A rich polymorphism with fine tuning of chemical and physical properties can be envisioned. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Phase Boundaries, Isotope Effect and Superconductivity of Lithium Under Hydrostatic Conditions

Authors: Stefano Racioppi, Iren Saffarian-Deemyad, William Holle, Francesco Belli, Richard Ferry, Curtis Kenney-Benson, Jesse S. Smith, Eva Zurek, Shanti Deemyad

Abstract: We present theoretical and experimental studies of superconductivity and low temperature structural phase boundaries in lithium. We mapped the structural phase diagram of 6Li and 7Li under hydrostatic conditions between 5 top 55GPa and within the temperature range of 15 to 75K, observing the FCC-hR1-cI16 phase transitions. 6Li and 7Li show some differences at the structural boundaries, with a pote… ▽ More We present theoretical and experimental studies of superconductivity and low temperature structural phase boundaries in lithium. We mapped the structural phase diagram of 6Li and 7Li under hydrostatic conditions between 5 top 55GPa and within the temperature range of 15 to 75K, observing the FCC-hR1-cI16 phase transitions. 6Li and 7Li show some differences at the structural boundaries, with a potential shift of the phase boundaries of 6Li to lower pressures. Density functional theory calculations and topological analysis of the electron density elucidates the superconducting properties and interatomic interactions within these phases of lithium. △ Less

Submitted 4 January, 2024; v1 submitted 29 December, 2023; originally announced December 2023.

Comments: (24 pages, 10 figures, 2 tables, 2 parts including main manuscript and the supplementary information section)

40 Years of SCES at Los Alamos

Authors: Z. Fisk, J. L. Smith, J. D. Thompson

Abstract: Reports of unconventional superconductivity in UBe13 in 1983 and soon thereafter of the possible coexistence of bulk superconductivity and spin fluctuations in UPt3 marked the beginning of a 40-year adventure in the study of strongly correlated quantum materials and phenomena at Los Alamos. The subsequent discovery and exploration of heavy-fermion magnetism, cuprates, Kondo insulators, Ce- and Pu-… ▽ More Reports of unconventional superconductivity in UBe13 in 1983 and soon thereafter of the possible coexistence of bulk superconductivity and spin fluctuations in UPt3 marked the beginning of a 40-year adventure in the study of strongly correlated quantum materials and phenomena at Los Alamos. The subsequent discovery and exploration of heavy-fermion magnetism, cuprates, Kondo insulators, Ce- and Pu-115 superconductors and, more broadly, quantum states of narrow-band systems provided challenges for the next 30 years. Progress was not made in a vacuum but benefitted from significant advances in the Americas, Asia and Europe as well as from essential collaborations, visitors and Los Alamos students and postdocs, many subsequently setting their own course in SCES. As often the case, serendipity played a role in shaping this history. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: Special session at the 2023 International Conference on Strongly Correlated Electron Systems

Journal ref: New Physics: Sae Mulli, vol. 73, No. 12 December 2023, pp.1023-1036

AI-guided inverse design and discovery of recyclable vitrimeric polymers

Authors: Yiwen Zheng, Prakash Thakolkaran, Agni K. Biswal, Jake A. Smith, Ziheng Lu, Shuxin Zheng, Bichlien H. Nguyen, Siddhant Kumar, Aniruddh Vashisth

Abstract: Vitrimer is a new, exciting class of sustainable polymers with the ability to heal due to their dynamic covalent adaptive network that can go through associative rearrangement reactions. However, a limited choice of constituent molecules restricts their property space, prohibiting full realization of their potential applications. To overcome this challenge, we couple molecular dynamics (MD) simula… ▽ More Vitrimer is a new, exciting class of sustainable polymers with the ability to heal due to their dynamic covalent adaptive network that can go through associative rearrangement reactions. However, a limited choice of constituent molecules restricts their property space, prohibiting full realization of their potential applications. To overcome this challenge, we couple molecular dynamics (MD) simulations and a novel graph variational autoencoder (VAE) machine learning model for inverse design of vitrimer chemistries with desired glass transition temperature (Tg) and synthesize a novel vitrimer polymer. We build the first vitrimer dataset of one million chemistries and calculate Tg on 8,424 of them by high-throughput MD simulations calibrated by a Gaussian process model. The proposed novel VAE employs dual graph encoders and a latent dimension overlapping scheme which allows for individual representation of multi-component vitrimers. By constructing a continuous latent space containing necessary information of vitrimers, we demonstrate high accuracy and efficiency of our framework in discovering novel vitrimers with desirable Tg beyond the training regime. To validate the effectiveness of our framework in experiments, we generate novel vitrimer chemistries with a target Tg = 323 K. By incorporating chemical intuition, we synthesize a vitrimer with Tg of 311-317 K, and experimentally demonstrate healability and flowability. The proposed framework offers an exciting tool for polymer chemists to design and synthesize novel, sustainable vitrimer polymers for a facet of applications. △ Less

Submitted 13 August, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

MatterGen: a generative model for inorganic materials design

Authors: Claudio Zeni, Robert Pinsler, Daniel Zügner, Andrew Fowler, Matthew Horton, Xiang Fu, Sasha Shysheya, Jonathan Crabbé, Lixin Sun, Jake Smith, Bichlien Nguyen, Hannes Schulz, Sarah Lewis, Chin-Wei Huang, Ziheng Lu, Yichi Zhou, Han Yang, Hongxia Hao, Jielan Li, Ryota Tomioka, Tian Xie

Abstract: The design of functional materials with desired properties is essential in driving technological advances in areas like energy storage, catalysis, and carbon capture. Generative models provide a new paradigm for materials design by directly generating entirely novel materials given desired property constraints. Despite recent progress, current generative models have low success rate in proposing s… ▽ More The design of functional materials with desired properties is essential in driving technological advances in areas like energy storage, catalysis, and carbon capture. Generative models provide a new paradigm for materials design by directly generating entirely novel materials given desired property constraints. Despite recent progress, current generative models have low success rate in proposing stable crystals, or can only satisfy a very limited set of property constraints. Here, we present MatterGen, a model that generates stable, diverse inorganic materials across the periodic table and can further be fine-tuned to steer the generation towards a broad range of property constraints. To enable this, we introduce a new diffusion-based generative process that produces crystalline structures by gradually refining atom types, coordinates, and the periodic lattice. We further introduce adapter modules to enable fine-tuning towards any given property constraints with a labeled dataset. Compared to prior generative models, structures produced by MatterGen are more than twice as likely to be novel and stable, and more than 15 times closer to the local energy minimum. After fine-tuning, MatterGen successfully generates stable, novel materials with desired chemistry, symmetry, as well as mechanical, electronic and magnetic properties. Finally, we demonstrate multi-property materials design capabilities by proposing structures that have both high magnetic density and a chemical composition with low supply-chain risk. We believe that the quality of generated materials and the breadth of MatterGen's capabilities represent a major advancement towards creating a universal generative model for materials design. △ Less

Submitted 29 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 13 pages main text, 35 pages supplementary information

Accurate Machine Learning Predictions of Coercivity in High-Performance Permanent Magnets

Authors: Churna Bhandari, Gavin N. Nop, Jonathan D. H. Smith, Durga Paudyal

Abstract: Increased demand for high-performance permanent magnets in the electric vehicle and wind turbine industries has prompted the search for cost-effective alternatives.Discovering new magnetic materials with the desired intrinsic and extrinsic permanent magnet properties presents a significant challenge to researchers because of issues with the global supply of rare-earth elements, material stability,… ▽ More Increased demand for high-performance permanent magnets in the electric vehicle and wind turbine industries has prompted the search for cost-effective alternatives.Discovering new magnetic materials with the desired intrinsic and extrinsic permanent magnet properties presents a significant challenge to researchers because of issues with the global supply of rare-earth elements, material stability, and a low maximum magnetic energy product BH$_{max}$.While first-principle density functional theory (DFT) predicts materials' magnetic moments, magneto-crystalline anisotropy constants, and exchange interactions, it cannot compute coercivity ($H_c$).Although it is possible to calculate $H_c$ theoretically with micromagnetic simulations, the predicted value is larger than the experiment by almost an order of magnitude, due to the Brown paradox.To circumvent these, we employ machine learning (ML) methods on an extensive database obtained from experiments, DFT calculations, and micromagnetic modeling.The use of a large dataset enables realistic $H_c$ predictions for materials such as Ce-doped Nd$_2$Fe$_{14}$B, comparing favorably against micromagnetically simulated coercivities.Remarkably, our ML model accurately identifies uniaxial magneto-crystalline anisotropy as the primary contributor to $H_c$. With DFT calculations, we predict the Nd-site dependent magnetic anisotropy behavior in Nd$_2$Fe$_{14}$B, confirming that Nd $4g$-sites mainly contribute to uniaxial magneto-crystalline anisotropy, and also calculate Curie temperature (T$_{C}$).Both calculated results are in good agreement with experiment.The coupled experimental dataset and ML modeling with DFT input predict $H_c$ with far greater accuracy and speed than was previously possible using micromagnetic modeling.Further, we reverse-engineer the inter-grain exchange coupling with micromagnetic simulations by employing the ML predictions. △ Less

Submitted 25 July, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 18 pages, 12 figures

Magnetic Tunnel Junction Random Number Generators Applied to Dynamically Tuned Probability Trees Driven by Spin Orbit Torque

Authors: Andrew Maicke, Jared Arzate, Samuel Liu, Jaesuk Kwon, J. Darby Smith, James B. Aimone, Shashank Misra, Catherine Schuman, Suma G. Cardwell, Jean Anne C. Incorvia

Abstract: Perpendicular magnetic tunnel junction (pMTJ)-based true-random number generators (RNG) can consume orders of magnitude less energy per bit than CMOS pseudo-RNG. Here, we numerically investigate with a macrospin Landau-Lifshitz-Gilbert equation solver the use of pMTJs driven by spin-orbit torque to directly sample numbers from arbitrary probability distributions with the help of a tunable probabil… ▽ More Perpendicular magnetic tunnel junction (pMTJ)-based true-random number generators (RNG) can consume orders of magnitude less energy per bit than CMOS pseudo-RNG. Here, we numerically investigate with a macrospin Landau-Lifshitz-Gilbert equation solver the use of pMTJs driven by spin-orbit torque to directly sample numbers from arbitrary probability distributions with the help of a tunable probability tree. The tree operates by dynamically biasing sequences of pMTJ relaxation events, called 'coinflips', via an additional applied spin-transfer-torque current. Specifically, using a single, ideal pMTJ device we successfully draw integer samples on the interval 0,255 from an exponential distribution based on p-value distribution analysis. In order to investigate device-to-device variations, the thermal stability of the pMTJs are varied based on manufactured device data. It is found that while repeatedly using a varied device inhibits ability to recover the probability distribution, the device variations average out when considering the entire set of devices as a 'bucket' to agnostically draw random numbers from. Further, it is noted that the device variations most significantly impact the highest level of the probability tree, iwth diminishing errors at lower levels. The devices are then used to draw both uniformly and exponentially distributed numbers for the Monte Carlo computation of a problem from particle transport, showing excellent data fit with the analytical solution. Finally, the devices are benchmarked against CMOS and memristor RNG, showing faster bit generation and significantly lower energy use. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 10 pages, 8 figures, 2 tables

The influence of pinholes and weak-points in aluminium-oxide Josephson junctions

Authors: K. Bayros, M. J. Cyster, J. S. Smith, J. H. Cole

Abstract: Josephson junctions are the key components used in superconducting qubits for quantum computing. The advancement of quantum computing is limited by a lack of stability and reproducibility of qubits which ultimately originates in the amorphous tunnel barrier of the Josephson junctions and other material imperfections. Pinholes in the junction have been suggested as one of the possible contributors… ▽ More Josephson junctions are the key components used in superconducting qubits for quantum computing. The advancement of quantum computing is limited by a lack of stability and reproducibility of qubits which ultimately originates in the amorphous tunnel barrier of the Josephson junctions and other material imperfections. Pinholes in the junction have been suggested as one of the possible contributors to these instabilities, but evidence of their existence and the effect they might have on transport is unclear. We use molecular dynamics to create three-dimensional atomistic models to describe Al-AlOx-Al tunnel junctions, showing that pinholes form when oxidation of the barrier is incomplete. Following this we use the atomistic model and simulate the electronic transport properties for tunnel junctions with different barrier thicknesses using the non-equilibrium Green's function formalism. We observe that pinholes may contribute to excess quasiparticle current flow in Al-AlOx-Al tunnel junctions with thinner barriers, and in thicker barriers we observe weak-points which facilitate leakage currents even when the oxide is continuous. We find that the disordered nature of the amorphous barrier results in significant variations in the transport properties. Additionally, we determine the current-phase relationship for our atomistic structures, confirming that devices with pinholes and weak-points cause a deviation from the ideal sinusoidal Josephson relationship. △ Less

Submitted 29 February, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 17 pages, 15 figures

Journal ref: Phys. Rev. Materials 8, 046202 (2024)

MOFDiff: Coarse-grained Diffusion for Metal-Organic Framework Design

Authors: Xiang Fu, Tian Xie, Andrew S. Rosen, Tommi Jaakkola, Jake Smith

Abstract: Metal-organic frameworks (MOFs) are of immense interest in applications such as gas storage and carbon capture due to their exceptional porosity and tunable chemistry. Their modular nature has enabled the use of template-based methods to generate hypothetical MOFs by combining molecular building blocks in accordance with known network topologies. However, the ability of these methods to identify t… ▽ More Metal-organic frameworks (MOFs) are of immense interest in applications such as gas storage and carbon capture due to their exceptional porosity and tunable chemistry. Their modular nature has enabled the use of template-based methods to generate hypothetical MOFs by combining molecular building blocks in accordance with known network topologies. However, the ability of these methods to identify top-performing MOFs is often hindered by the limited diversity of the resulting chemical space. In this work, we propose MOFDiff: a coarse-grained (CG) diffusion model that generates CG MOF structures through a denoising diffusion process over the coordinates and identities of the building blocks. The all-atom MOF structure is then determined through a novel assembly algorithm. Equivariant graph neural networks are used for the diffusion model to respect the permutational and roto-translational symmetries. We comprehensively evaluate our model's capability to generate valid and novel MOF structures and its effectiveness in designing outstanding MOF materials for carbon capture applications with molecular simulations. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 19 pages, 12 figures

Optical and Acoustic Phonons in Turbostratic and Cubic Boron Nitride Thin Films on Diamond Substrates

Authors: Erick Guzman, Fariborz Kargar, Avani Patel, Saurabh Vishwakarma, Dylan Wright, Richard B. Wilson, David J. Smith, Robert J. Nemanich, Alexander A. Balandin

Abstract: We report an investigation of the bulk optical, bulk acoustic, and surface acoustic phonons in thin films of turbostratic boron nitride (t-BN) and cubic boron nitride (c-BN) grown on B-doped polycrystalline and single-crystalline diamond (001) and (111) substrates. The characteristics of different types of phonons were determined using Raman and Brillouin-Mandelstam light scattering spectroscopies… ▽ More We report an investigation of the bulk optical, bulk acoustic, and surface acoustic phonons in thin films of turbostratic boron nitride (t-BN) and cubic boron nitride (c-BN) grown on B-doped polycrystalline and single-crystalline diamond (001) and (111) substrates. The characteristics of different types of phonons were determined using Raman and Brillouin-Mandelstam light scattering spectroscopies. The atomic structure of the films was determined using high-resolution transmission electron microscopy (HRTEM) and correlated with the Raman and Brillouin-Mandelstam spectroscopy data. The HRTEM analysis revealed that the cubic boron nitride thin films consisted of a mixture of c-BN and t-BN phases, with c-BN being the dominant phase. It was found that while visible Raman spectroscopy provided information for characterizing the t-BN phase, it faced challenges in differentiating the c-BN phase either due to the presence of high-density defects or the overlapping of the Raman features with those from the B-doped diamond substrates. In contrast, Brillouin-Mandelstam spectroscopy clearly distinguishes the bulk longitudinal and surface acoustic phonons of the c-BN thin films grown on diamond substrates. Additionally, the angle-dependent surface Brillouin-Mandelstam scattering data show the peaks associated with the Rayleigh surface acoustic waves, which have higher phase velocities in c-BN films on diamond (111) substrates. These findings provide valuable insights into the phonon characteristics of the c-BN and diamond interfaces and have important implications for the thermal management of electronic devices based on ultra-wide-band-gap materials. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 27 pages; 4 figures

Improved Flexible Coaxial Ribbon Cable for High-Density Superconducting Arrays

Authors: Jennifer Pearl Smith, Benjamin A. Mazin, Alirio Boaventura, Kyle J. Thompson, Miguel Daal

Abstract: Superconducting arrays often require specialized, high-density cryogenic cabling capable of transporting electrical signals across temperature stages with minimal loss, crosstalk, and thermal conductivity. We report improvements to the design and fabrication of previously published superconducting 53 wt% Nb-47 wt% Ti (Nb47Ti) FLexible coAXial ribbon cables (FLAX). We used 3D electromagnetic simula… ▽ More Superconducting arrays often require specialized, high-density cryogenic cabling capable of transporting electrical signals across temperature stages with minimal loss, crosstalk, and thermal conductivity. We report improvements to the design and fabrication of previously published superconducting 53 wt% Nb-47 wt% Ti (Nb47Ti) FLexible coAXial ribbon cables (FLAX). We used 3D electromagnetic simulations to inform design changes to improve the characteristic impedance of the cable and the connector transition. We increased the center conductor diameter from 0.003 inches to 0.005 inches which lowered the cable characteristic impedance from $\sim$60 $Ω$ to $\sim$53 $Ω$. This change had a negligible impact on the computed heat load which we estimate to be 5 nW per trace from 1 K to 100 mK with a 1-ft cable. This is approximately half the heat load calculated for the smallest commercially available superconducting coax. We also modified the transition board to include a capacitive coupling between the upper ground plane and signal traces that mitigates the inductive transition. We tested these changes in a 5-trace, 1-ft long cable at 4 K and found the microwave transmission improved from 6 dB to 1.5 dB of attenuation at 8 GHz. This loss is comparable to commercial superconducting coax and 3$\times$ lower than commercial NbTi-on-polyimide flex cables at 8 GHz. The nearest-neighbor forward crosstalk remained less than -40 dB at 8 GHz. We compare key performance metrics with commercially available superconducting coax and NbTi-on-polyimide flex cables and we share initial progress on commercialization of this technology by Maybell Quantum Industries. △ Less

Submitted 29 June, 2023; v1 submitted 23 June, 2023; originally announced June 2023.

Effect of a hybrid transition moment on Stark-modulated photon echoes in Er$^{3+}$:Y$_2$SiO$_5$

Authors: Rose L. Ahlefeldt, Alexey Lyasota, Jodie Smith, Jinliang Ren, Matthew J. Sellars

Abstract: The 1538 nm ${}^4$I$_{15/2}$ - ${}^4$I$_{13/2}$ transition of Er$^{3+}$ has an unusual hybrid electric-magnetic dipole character, and signatures of the hybrid moment can be expected in coherent transient measurements. Here, we investigate the effect of the hybrid moment in both sites of Er$^{3+}$:Y$_2$SiO$_5$ on Stark-modulated photon echo measurements, showing that it results in a reduction of vi… ▽ More The 1538 nm ${}^4$I$_{15/2}$ - ${}^4$I$_{13/2}$ transition of Er$^{3+}$ has an unusual hybrid electric-magnetic dipole character, and signatures of the hybrid moment can be expected in coherent transient measurements. Here, we investigate the effect of the hybrid moment in both sites of Er$^{3+}$:Y$_2$SiO$_5$ on Stark-modulated photon echo measurements, showing that it results in a reduction of visibility of the modulated signal as well as phase and polarization shifts. We interpret these effects using a simple optical Bloch equation model, showing that site 1 has a strongly mixed moment and site 2 is predominantly magnetic dipole. We discuss the implications of the hybrid moment for quantum information applications of quantum memories. We also use the echo measurements to extract the Stark shift of the optical transition along three orthogonal directions, finding values between 10.50 and 11.93 kHz/(V/cm) for site 1 and 1.61 and 15.35 kHz/(V/cm) for site 2. We observe a modification of the Zeeman shift by the electric field in both sites and discuss how this may be used to electrically control Er$^{3+}$ spin qubits. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 10 pages, 6 figures, supplemental information included

Signatures of pressure-enhanced helimagnetic order in van der Waals multiferroic NiI$_2$

Authors: Connor A. Occhialini, Luiz G. P. Martins, Qian Song, Jesse S. Smith, Jesse Kapeghian, Danila Amoroso, Joshua J. Sanchez, Paolo Barone, Bertrand Dupé, Matthieu j. Verstraete, Jing Kong, Antia S. Botana, Riccardo Comin

Abstract: The van der Waals (vdW) type-II multiferroic NiI$_2$ has emerged as a candidate for exploring non-collinear magnetism and magnetoelectric effects in the 2D limit. Frustrated intralayer exchange interactions on a triangular lattice result in a helimagnetic ground state, with spin-induced improper ferroelectricity stabilized by the interlayer interactions. Here we investigate the magnetic and struct… ▽ More The van der Waals (vdW) type-II multiferroic NiI$_2$ has emerged as a candidate for exploring non-collinear magnetism and magnetoelectric effects in the 2D limit. Frustrated intralayer exchange interactions on a triangular lattice result in a helimagnetic ground state, with spin-induced improper ferroelectricity stabilized by the interlayer interactions. Here we investigate the magnetic and structural phase transitions in bulk NiI$_2$, using high-pressure Raman spectroscopy, optical linear dichroism, and x-ray diffraction. We obtain evidence for a significant pressure enhancement of the antiferromagnetic and helimagnetic transition temperatures, at rates of $\sim15.3/14.4$ K/GPa, respectively. These enhancements are attributed to a cooperative effect of pressure-enhanced interlayer and third-nearest-neighbor intralayer exchange. These results reveal a general path for obtaining high-temperature type-II multiferroicity via high pressures in vdW materials. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Effects of Pressure on the Electronic and Magnetic Properties of Bulk NiI$_{2}$

Authors: Jesse Kapeghian, Danila Amoroso, Connor A. Occhialini, Luiz G. P. Martins, Qian Song, Jesse S. Smith, Joshua J. Sanchez, Jing Kong, Riccardo Comin, Paolo Barone, Bertrand Dupé, Matthieu J. Verstraete, Antia S. Botana

Abstract: Transition metal dihalides have recently garnered interest in the context of two-dimensional van der Waals magnets as their underlying geometrically frustrated triangular lattice leads to interesting competing exchange interactions. In particular, NiI$_{2}$ is a magnetic semiconductor that has been long known for its exotic helimagnetism in the bulk. Recent experiments have shown that the helimagn… ▽ More Transition metal dihalides have recently garnered interest in the context of two-dimensional van der Waals magnets as their underlying geometrically frustrated triangular lattice leads to interesting competing exchange interactions. In particular, NiI$_{2}$ is a magnetic semiconductor that has been long known for its exotic helimagnetism in the bulk. Recent experiments have shown that the helimagnetic state survives down to the monolayer limit with a layer-dependent magnetic transition temperature that suggests a relevant role of the interlayer coupling. Here, we explore the effects of hydrostatic pressure as a means to enhance this interlayer exchange and ultimately tune the electronic and magnetic response of NiI$_{2}$. We study first the evolution of the structural parameters as a function of external pressure using first-principles calculations combined with x-ray diffraction measurements. We then examine the evolution of the electronic structure and magnetic exchange interactions via first-principles calculations and Monte Carlo simulations. We find that the leading interlayer coupling is an antiferromagnetic second-nearest neighbor interaction that increases monotonically with pressure. The ratio between isotropic third- and first-nearest neighbor intralayer exchanges, which controls the magnetic frustration and determines the magnetic propagation vector $\mathbf{q}$ of the helimagnetic ground state, is also enhanced by pressure. As a consequence, our Monte Carlo simulations show a monotonic increase in the magnetic transition temperature, indicating that pressure is an effective means to tune the magnetic response of NiI$_{2}$. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Assessing the stability fields of molecular and polymeric CO2

Authors: Alexander F. Goncharov, Elena Bykova, Maxim Bykov, Eric Edmund, Jesse S. Smith, Stella Chariton, Vitali B. Prakapenka

Abstract: We investigated the stability of polymeric CO2 over a wide range of pressures, temperatures, and chemical environments. We find that the I-42d polymeric structure, consisting of a three-dimensional network of corner sharing CO4 tetrahedra, forms at 40-140 GPa and from a CO-N2 mixture at 39 GPa. An exceptional stability field of 0 to 286 GPa and 100 to 2500 K is documented for this structure, makin… ▽ More We investigated the stability of polymeric CO2 over a wide range of pressures, temperatures, and chemical environments. We find that the I-42d polymeric structure, consisting of a three-dimensional network of corner sharing CO4 tetrahedra, forms at 40-140 GPa and from a CO-N2 mixture at 39 GPa. An exceptional stability field of 0 to 286 GPa and 100 to 2500 K is documented for this structure, making it a viable candidate for planetary interiors. The stability of the tetrahedral polymeric motif of CO2-V is a consequence of the rigidity of sp3 hybridized orbitals of carbon in a closed-packed oxygen sublattice. △ Less

Submitted 21 April, 2023; originally announced April 2023.

Comments: 15 pages, 9 figures Supplementary Materials 6 pages, 1 table

Plastic strain-induced phase transformations in silicon: drastic reduction of transformation pressures, change in transformation sequence, and particle size effect

Authors: Sorb Yesudhas, Valery I. Levitas, Feng Lin, K. K. Pandey, Jesse Smith

Abstract: Pressure-induced phase transformations (PTs) between numerous phases of Si, the most important electronic material, have been studied for decades. This is not the case for plastic strain-induced PTs. Here, we revealed in-situ various unexpected plastic strain-induced PT phenomena. Thus, for 100 nm Si, strain-induced PT Si-I to Si-II (and Si-I to Si-III) initiates at 0.4 GPa (0.6 GPa) versus 16.2 G… ▽ More Pressure-induced phase transformations (PTs) between numerous phases of Si, the most important electronic material, have been studied for decades. This is not the case for plastic strain-induced PTs. Here, we revealed in-situ various unexpected plastic strain-induced PT phenomena. Thus, for 100 nm Si, strain-induced PT Si-I to Si-II (and Si-I to Si-III) initiates at 0.4 GPa (0.6 GPa) versus 16.2 GPa ($\infty$, since it does not occur) under hydrostatic conditions; for 30 nm Si, it is 6.1 GPa versus $\infty$. The predicted theoretical correlation between the direct and inverse Hall-Petch effect of the grain size on the yield strength and the minimum pressure for strain-induced PT is confirmed for the appearance of Si-II. Retaining Si-II at ambient pressure and obtaining reverse Si-II to Si-I PT are achieved, demonstrating the possibilities of manipulating different synthetic paths. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 26 pages, 10 figures

Transonic Dislocation Propagation in Diamond

Authors: Kento Katagiri, Tatiana Pikuz, Lichao Fang, Bruno Albertazzi, Shunsuke Egashira, Yuichi Inubushi, Genki Kamimura, Ryosuke Kodama, Michel Koenig, Bernard Kozioziemski, Gooru Masaoka, Kohei Miyanishi, Hirotaka Nakamura, Masato Ota, Gabriel Rigon, Youichi Sakawa, Takayoshi Sano, Frank Schoofs, Zoe J. Smith, Keiichi Sueda, Tadashi Togashi, Tommaso Vinci, Yifan Wang, Makina Yabashi, Toshinori Yabuuchi , et al. (2 additional authors not shown)

Abstract: The motion of line defects (dislocations) has been studied for over 60 years but the maximum speed at which they can move is unresolved. Recent models and atomistic simulations predict the existence of a limiting velocity of dislocation motions between the transonic and subsonic ranges at which the self-energy of dislocation diverges, though they do not deny the possibility of the transonic disloc… ▽ More The motion of line defects (dislocations) has been studied for over 60 years but the maximum speed at which they can move is unresolved. Recent models and atomistic simulations predict the existence of a limiting velocity of dislocation motions between the transonic and subsonic ranges at which the self-energy of dislocation diverges, though they do not deny the possibility of the transonic dislocations. We use femtosecond x-ray radiography to track ultrafast dislocation motion in shock-compressed single-crystal diamond. By visualizing stacking faults extending faster than the slowest sound wave speed of diamond, we show the evidence of partial dislocations at their leading edge moving transonically. Understanding the upper limit of dislocation mobility in crystals is essential to accurately model, predict, and control the mechanical properties of materials under extreme conditions. △ Less

Submitted 6 October, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

The human factor: results of a small-angle scattering data analysis Round Robin

Authors: Brian R. Pauw, Glen J. Smales, Andy S. Anker, Daniel M. Balazs, Frederick L. Beyer, Ralf Bienert, Wim G. Bouwman, Ingo Breßler, Joachim Breternitz, Erik S Brok, Gary Bryant, Andrew J. Clulow, Erin R. Crater, Frédéric De Geuser, Alessandra Del Giudice, Jérôme Deumer, Sabrina Disch, Shankar Dutt, Kilian Frank, Emiliano Fratini, Elliot P. Gilbert, Marc Benjamin Hahn, James Hallett, Max Hohenschutz, Martin Hollamby , et al. (24 additional authors not shown)

Abstract: A Round Robin study has been carried out to estimate the impact of the human element in small-angle scattering data analysis. Four corrected datasets were provided to participants ready for analysis. All datasets were measured on samples containing spherical scatterers, with two datasets in dilute dispersions, and two from powders. Most of the 46 participants correctly identified the number of pop… ▽ More A Round Robin study has been carried out to estimate the impact of the human element in small-angle scattering data analysis. Four corrected datasets were provided to participants ready for analysis. All datasets were measured on samples containing spherical scatterers, with two datasets in dilute dispersions, and two from powders. Most of the 46 participants correctly identified the number of populations in the dilute dispersions, with half of the population mean entries within 1.5% and half of the population width entries within 40%, respectively. Due to the added complexity of the structure factor, much fewer people submitted answers on the powder datasets. For those that did, half of the entries for the means and widths were within 44% and 86% respectively. This Round Robin experiment highlights several causes for the discrepancies, for which solutions are proposed. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 23 pages, 10 figures. For the original information sent to RR participants, see https://zenodo.org/record/7506365 . For the anonymized results and Jupyter notebook for analysis, see https://zenodo.org/record/7509710

Supersolidity and crystallization of a dipolar Bose gas in an infinite tube

Authors: Joseph C. Smith, D. Baillie, P. B. Blakie

Abstract: We calculate the ground states of a dipolar Bose gas confined in an infinite tube potential. We use the extended Gross-Pitaevskii equation theory and present a novel numerical method to efficiently obtain solutions. A key feature of this method is an analytic result for a truncated dipole-dipole interaction potential that enables the long-ranged interactions to be accurately evaluated within a uni… ▽ More We calculate the ground states of a dipolar Bose gas confined in an infinite tube potential. We use the extended Gross-Pitaevskii equation theory and present a novel numerical method to efficiently obtain solutions. A key feature of this method is an analytic result for a truncated dipole-dipole interaction potential that enables the long-ranged interactions to be accurately evaluated within a unit cell. Our focus is on the transition of the ground state to a crystal driven by dipole-dipole interactions as the short ranged interaction strength is varied. We find that the transition is continuous or discontinuous depending upon average system density. These results give deeper insight into the supersolid phase transition observed in recent experiments, and validate the utility of the reduced three-dimensional theory developed in [Phys. Rev. Res. 2, 043318 (2020)] for making qualitatively accurate predictions. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: 10 pages, 7 figures

Journal ref: Phys. Rev. A 107, 033301 (2023)

Random Bitstream Generation using Voltage-Controlled Magnetic Anisotropy and Spin Orbit Torque Magnetic Tunnel Junctions

Authors: Samuel Liu, Jaesuk Kwon, Paul W. Bessler, Suma Cardwell, Catherine Schuman, J. Darby Smith, James B. Aimone, Shashank Misra, Jean Anne C. Incorvia

Abstract: Probabilistic computing using random number generators (RNGs) can leverage the inherent stochasticity of nanodevices for system-level benefits. The magnetic tunnel junction (MTJ) has been studied as an RNG due to its thermally-driven magnetization dynamics, often using spin transfer torque (STT) current amplitude to control the random switching of the MTJ free layer magnetization, here called the… ▽ More Probabilistic computing using random number generators (RNGs) can leverage the inherent stochasticity of nanodevices for system-level benefits. The magnetic tunnel junction (MTJ) has been studied as an RNG due to its thermally-driven magnetization dynamics, often using spin transfer torque (STT) current amplitude to control the random switching of the MTJ free layer magnetization, here called the stochastic write method. There are additional knobs to control the MTJ-RNG, including voltage-controlled magnetic anisotropy (VCMA) and spin orbit torque (SOT), and there is need to systematically study and compared these methods. We build an analytical model of the MTJ to characterize using VCMA and SOT to generate random bit streams. The results show that both methods produce high quality, uniformly distributed bitstreams. Biasing the bitstreams using either STT current or an applied magnetic field shows an approximately sigmoidal distribution vs. bias amplitude for both VCMA and SOT, compared to less sigmoidal for stochastic write for the same conditions. The energy consumption per sample is calculated to be 0.1 pJ (SOT), 1 pJ (stochastic write), and 20 pJ (VCMA), revealing the potential energy benefit of using SOT and showing using VCMA may require higher damping materials. The generated bitstreams are then applied to two example tasks: generating an arbitrary probability distribution using the Hidden Correlation Bernoulli Coins method and using the MTJ-RNGs as stochastic neurons to perform simulated annealing in a Boltzmann machine, where both VCMA and SOT methods show the ability to effectively minimize the system energy with small delay and low energy. These results show the flexibility of the MTJ as a true RNG and elucidate design parameters for optimizing the device operation for applications. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: 9 pages, 6 figures

Ethylenediamine Addition Improves Performance and Suppresses Phase Instabilities in Mixed-Halide Perovskites

Authors: Margherita Taddei, Joel A. Smith, Benjamin M. Gallant, Suer Zhou, Robert J. E. Westbrook, Yangwei Shi, Jian Wang, James N. Drysdale, Declan P. McCarthy, Stephen Barlow, Seth R. Marder, Henry J. Snaith, David S. Ginger

Abstract: We show that adding ethylenediamine (EDA) to perovskite precursor solution improves the photovoltaic device performance and material stability of high-bromide-content, methylammonium-free, formamidinium cesium lead halide perovskites FA1-xCsxPb(I1-yBry)3 which are currently of interest for perovskite-on-Si tandem solar cells. Using spectroscopy and hyperspectral microscopy, we show that the additi… ▽ More We show that adding ethylenediamine (EDA) to perovskite precursor solution improves the photovoltaic device performance and material stability of high-bromide-content, methylammonium-free, formamidinium cesium lead halide perovskites FA1-xCsxPb(I1-yBry)3 which are currently of interest for perovskite-on-Si tandem solar cells. Using spectroscopy and hyperspectral microscopy, we show that the additive improves film homogeneity and suppresses the phase instability that is ubiquitous in high-Br perovskite formulations, producing films that remain stable for over 100 days in ambient conditions. With the addition of 1 mol% EDA we demonstrate 1.69 eV-gap perovskite single-junction p-i-n devices with a VOC of 1.22 V, and a champion maximum power point tracked power conversion efficiency of 18.8%, comparable to the best reported methylammonium-free perovskites. Using nuclear magnetic resonance (NMR) spectroscopy and X-ray diffraction techniques, we show that EDA reacts with FA+ in solution, rapidly and quantitatively forming imidazolinium cations. It is the presence of imidazolinium during crystallization which drives the improved perovskite thin-film properties. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Stochastic magnetic actuated random transducer devices based on perpendicular magnetic tunnel junctions

Authors: Laura Rehm, Corrado Carlo Maria Capriata, Misra Shashank, J. Darby Smith, Mustafa Pinarbasi, B. Gunnar Malm, Andrew D. Kent

Abstract: True random number generators are of great interest in many computing applications such as cryptography, neuromorphic systems and Monte Carlo simulations. Here we investigate perpendicular magnetic tunnel junction nanopillars (pMTJs) activated by short duration (ns) pulses in the ballistic limit for such applications. In this limit, a pulse can transform the Boltzmann distribution of initial free… ▽ More True random number generators are of great interest in many computing applications such as cryptography, neuromorphic systems and Monte Carlo simulations. Here we investigate perpendicular magnetic tunnel junction nanopillars (pMTJs) activated by short duration (ns) pulses in the ballistic limit for such applications. In this limit, a pulse can transform the Boltzmann distribution of initial free layer magnetization states into randomly magnetized down or up states, i.e. a bit that is 0 or 1, easily determined by measurement of the junction's tunnel resistance. It is demonstrated that bitstreams with millions of events: 1) are very well described by the binomial distribution; 2) can be used to create a uniform distribution of 8-bit random numbers; 3) pass multiple statistical tests for true randomness, including all the National Institute of Standards tests for random number generators with only one XOR operation; and 4) can have no drift in the bit probability with time. The results presented here show that pMTJs operated in the ballistic regime can generate true random numbers at GHz bitrates, while being more robust to environmental changes, such as their operating temperature, compared to other stochastic nanomagnetic devices. △ Less

Submitted 15 September, 2022; v1 submitted 3 September, 2022; originally announced September 2022.

Comments: 8 pages, 3 figures, 2 tables, will be submitted to peer-reviewed journal

Journal ref: Phys. Rev. Applied 19, 024035 (2023)

Growth of Ultrathin Bi$_2$Se$_3$ Films by Molecular Beam Epitaxy

Authors: Saadia Nasir, Walter J. Smith, Thomas E. Beechem, Stephanie Law

Abstract: Bi$_2$Se$_3$ is a widely studied 3D topological insulator having potential applications in optics, electronics, and spintronics. When the thickness of these films decrease to less than approximately 6 nm, the top and bottom surface states couple, resulting in the opening of a small gap at the Dirac point. In the 2D limit, Bi$_2$Se$_3$ may exhibit quantum spin Hall states. However, growing coalesce… ▽ More Bi$_2$Se$_3$ is a widely studied 3D topological insulator having potential applications in optics, electronics, and spintronics. When the thickness of these films decrease to less than approximately 6 nm, the top and bottom surface states couple, resulting in the opening of a small gap at the Dirac point. In the 2D limit, Bi$_2$Se$_3$ may exhibit quantum spin Hall states. However, growing coalesced ultra-thin Bi$_2$Se$_3$ films with a controllable thickness and typical triangular domain morphology in the few nanometer range is challenging. Here, we explore the growth of Bi$_2$Se$_3$ films having thickness down to 4 nm on sapphire substrates using molecular beam epitaxy that were then characterized with Hall measurements, atomic force microscopy, and Raman imaging. We find that substrate pre-treatment -- growing and decomposing a few layers of \BiSe before the actual deposition -- is critical to obtaining a completely coalesced film. In addition, higher growth rates and lower substrate temperatures led to improvement in surface roughness, in contrast to what is observed for conventional epitaxy. Overall, coalesced ultra-thin Bi$_2$Se$_3$ films with lower surface roughness enables thickness-dependent studies across the transition from a 3D-topological insulator to one with gapped surface states in the 2D regime. △ Less

Submitted 31 August, 2022; v1 submitted 30 August, 2022; originally announced August 2022.

Ambient temperature pressure driven alkane dehydrogenation by palladium metal

Authors: Mungo Frost, Emma E. McBride, Dean Smith, Jesse S. Smith, Siegfried Glenzer

Abstract: Dehydrogenation of alkanes is of increasing importance in fulfilling global demand for olefins and offers a potential source of carbon-neutral hydrogen as a co-product. Currently commercial dehydrogenation processes occur at high-temperatures (500-900$^\circ$C) which is energy intensive and results in side reactions and rapid coking of the catalysts. In addition the hydrogen produced is often burn… ▽ More Dehydrogenation of alkanes is of increasing importance in fulfilling global demand for olefins and offers a potential source of carbon-neutral hydrogen as a co-product. Currently commercial dehydrogenation processes occur at high-temperatures (500-900$^\circ$C) which is energy intensive and results in side reactions and rapid coking of the catalysts. In addition the hydrogen produced is often burned to maintain temperature and to inhibit the back reaction. Here we demonstrate pressure as a parameter to enable novel chemical catalytic processes and demonstrate ambient-temperature dehydrogenation of alkanes by palladium at 50-100 MPa pressures, with both hydrogen gas and olefins recovered on decompression. This reaction follows a fundamentally different path to current commercial high-temperature low-pressure dehydrogenation processes with the palladium catalyst reversibly forming a hydride intermediate. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Main article: 12 pages and 8 figures, Supplemental follows main article with 4 pages and 3 figures

The Validity Window of Space-Charge-Limited Current Measurements of Metal Halide Perovskite Devices

Authors: William R. Fisher, Philip Calado, Jason A. Röhr, Joel A. Smith, Xingyuan Shi, Onkar Game, Jenny Nelson, Piers R. F. Barnes

Abstract: Space-charge-limited current (SCLC) measurements are used to estimate charge carrier mobilities and electronic trap densities of semiconductors by analysing the current density-voltage (JV) relationship for unipolar devices predicted by the Mott-Gurney (MG) law. However, the interpretation of SCLC measurements for metal-halide perovskites is problematic due to mobile ionic defects which redistribu… ▽ More Space-charge-limited current (SCLC) measurements are used to estimate charge carrier mobilities and electronic trap densities of semiconductors by analysing the current density-voltage (JV) relationship for unipolar devices predicted by the Mott-Gurney (MG) law. However, the interpretation of SCLC measurements for metal-halide perovskites is problematic due to mobile ionic defects which redistribute to screen electrostatic fields in devices during measurements. To overcome this, an SCLC measurement protocol was recently suggested that minimises ionic charge redistribution by probing the current during millisecond voltage pulses superimposed on a background bias. Here, we use drift-diffusion simulations with mobile ions to assess the validity of the MG law for analysing both the standard and new protocol JV measurements. We simulated idealised perovskite devices with differing mobile ion densities and compared them with simulations and measurements of devices with typical contact materials. We found the validity region for the MG law is limited to perovskites with mobile ion densities lower than the device's equilibrium charge carrier density (<10^17 cm-3 for 400 nm thick methylammonium lead iodide films) and contacts with injection/extraction barriers <=0.1 eV. The latter limitation can be partially overcome by increasing the device thickness, whereas the former limitation cannot. This restricts the range of perovskite layer compositions and viable contact materials that can be reliably analysed with the MG law. Approaches such as estimating trap densities from the apparent voltage onset to trap-free SCLC regime should also be critically reviewed since they rely on the same potentially invalid assumptions as the MG law. Our results demonstrate that extracting meaningful and accurate values for metal halide perovskite material properties from SCLC maybe challenging, or often not possible. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: Main text 22 pages, supplementary materials 9 pages

Full Configuration Interaction Excited-State Energies in Large Active Spaces from Subspace Iteration with Repeated Random Sparsification

Authors: Samuel M. Greene, Robert J. Webber, James E. T. Smith, Jonathan Weare, Timothy C. Berkelbach

Abstract: We present a stable and systematically improvable quantum Monte Carlo (QMC) approach to calculating excited-state energies, which we implement using our fast randomized iteration method for the full configuration interaction problem (FCI-FRI). Unlike previous excited-state quantum Monte Carlo methods, our approach, which is an asymmetric variant of subspace iteration, avoids the use of dot product… ▽ More We present a stable and systematically improvable quantum Monte Carlo (QMC) approach to calculating excited-state energies, which we implement using our fast randomized iteration method for the full configuration interaction problem (FCI-FRI). Unlike previous excited-state quantum Monte Carlo methods, our approach, which is an asymmetric variant of subspace iteration, avoids the use of dot products of random vectors and instead relies upon trial vectors to maintain orthogonality and estimate eigenvalues. By leveraging recent advances, we apply our method to calculate ground- and excited-state energies of strongly correlated molecular systems in large active spaces, including the carbon dimer with 8 electrons in 108 orbitals (8e,108o), an oxo-Mn(salen) transition metal complex (28e,28o), ozone (18e,87o), and butadiene (22e,82o). In the majority of these test cases, our approach yields total excited-state energies that agree with those from state-of-the-art methods -- including heat-bath CI, the density matrix renormalization group approach, and FCIQMC -- to within sub-milliHartree accuracy. In all cases, estimated excitation energies agree to within about 0.1 eV. △ Less

Submitted 12 October, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: 16 pages, 5 figures, 3 tables

AutoDisk: Automated Diffraction Processing and Strain Mapping in 4D-STEM

Authors: Sihan Wang, Tim Eldred, Jacob Smith, Wenpei Gao

Abstract: Development in lattice strain mapping using four-dimensional scanning transmission electron microscopy (4D-STEM) method now offers improved precision and feasibility. However, automatic and accurate diffraction analysis is still challenging due to noise and the complexity of intensity in diffraction patterns. In this work, we demonstrate an approach, employing the blob detection on cross-correlate… ▽ More Development in lattice strain mapping using four-dimensional scanning transmission electron microscopy (4D-STEM) method now offers improved precision and feasibility. However, automatic and accurate diffraction analysis is still challenging due to noise and the complexity of intensity in diffraction patterns. In this work, we demonstrate an approach, employing the blob detection on cross-correlated diffraction patterns followed by lattice fitting algorithm, to automate the processing of four-dimensional data, including identifying and locating disks, and extracting local lattice parameters without prior knowledge about the material. The approach is both tested using simulated diffraction patterns and applied on experimental data acquired from a Pd@Pt core-shell nanoparticle. Our method shows robustness against various sample thicknesses and high noise, capability to handle complex patterns, and picometer-scale accuracy in strain measurement, making it a promising tool for high-throughput 4D-STEM data processing. △ Less

Submitted 2 January, 2022; originally announced January 2022.

Compositional Control and Optimization of Molecular Beam Epitaxial Growth of (Sb2Te3)x(MnSb2Te4)y Magnetic Topological Insulators

Authors: Ido Levy, Candice Forrester, Haiming Deng, Manuel Goldan, Martha R. McCartney, David J. Smith, Christophe Testelin, Lia Krusin-Elbaum, Maria C. Tamargo

Abstract: Magnetic topological insulators such as MnBi2Te4 and MnSb2Te4 are promising hosts of novel physical phenomena such as quantum anomalous Hall effect and intrinsic axion insulator state, both potentially important for the implementation in topological spintronics and error-free quantum computing. In the bulk, the materials are antiferromagnetic but appropriate stacking with non-magnetic layers or ex… ▽ More Magnetic topological insulators such as MnBi2Te4 and MnSb2Te4 are promising hosts of novel physical phenomena such as quantum anomalous Hall effect and intrinsic axion insulator state, both potentially important for the implementation in topological spintronics and error-free quantum computing. In the bulk, the materials are antiferromagnetic but appropriate stacking with non-magnetic layers or excess Mn in the crystal lattice can induce a net ferromagnetic alignment. In this work we report the growth of (Sb2Te3)x(MnSb2Te4)y layers with varying Mn content by molecular beam epitaxy. The Mn flux fraction provided during growth controls the percent of MnSb2Te4 that is formed in the resulting layers by a self-assembly process. Highly crystalline layers with compositions varying between Sb2Te3 (y=0) and MnSb2Te4 (x=0) were obtained. The results show that Mn incorporates as a structural component to form MnSb2Te4, and as an impurity element both in Sb2Te3 and in MnSb2Te4. Two modifications of the growth conditions were implemented to enhance the incorporation of Mn as a structural element to form MnSb2Te4. Annealing of a thin portion of the layer at the beginning of growth (pre-anneal step), and increasing the growth temperature, both result in a larger percent of MnSb2Te4 for similar Mn flux fractions during growth. Samples having at least a few percent of MnSb2Te4 layers exhibit ferromagnetic behavior likely due to the excess Mn in the system which stabilizes on Sb sites as MnSb antisite defects. △ Less

Submitted 12 December, 2021; originally announced December 2021.

Lower pressure phases and metastable states of superconducting photo-induced carbonaceous sulfur hydride

Authors: G. Alexander Smith, Ines E. Collings, Elliot Snider, Dean Smith, Sylvain Petitgirard, Jesse Smith, Melanie White, Elyse Jones, Paul Ellison, Keith V. Lawler, Ranga P. Dias, Ashkan Salamat

Abstract: Room-temperature superconductivity was recently discovered in carbonaceous sulfur hydride (C-S-H) close to 3\,Mbar. We report significant differences in the superconducting response of C-S-H, with a maximum $T_{C}$ of 191(1)\,K, below a 1\,Mbar. Variations in intensity of the C-H Raman modes reveal carbon content can vary between crystals synthesized with the same photo-induced method. Synchrotron… ▽ More Room-temperature superconductivity was recently discovered in carbonaceous sulfur hydride (C-S-H) close to 3\,Mbar. We report significant differences in the superconducting response of C-S-H, with a maximum $T_{C}$ of 191(1)\,K, below a 1\,Mbar. Variations in intensity of the C-H Raman modes reveal carbon content can vary between crystals synthesized with the same photo-induced method. Synchrotron single crystal x-ray diffraction identifies polymorphism with increasing degrees of covalency. These unique metastable states are highly sensitive to thermodynamic pathways. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Microscopic quantum point contact formation as the electromigration mechanism in granular superconductor nanowires

Authors: T. C. Bartolo, J. S. Smith, Yannick Schön, Jan Nicolas Voss, M. J. Cyster, A. V. Ustinov, H. Rotzinger, J. H. Cole

Abstract: Granular aluminium is a high kinetic inductance thin film superconductor which, when formed into nanowires can undergo an intrinsic electromigration process. We use a combination of experimental and computational approaches to investigate the role of grain morphology and distribution in granular aluminium thin films, when formed into nanowire constrictions. Treating the granular aluminium film as… ▽ More Granular aluminium is a high kinetic inductance thin film superconductor which, when formed into nanowires can undergo an intrinsic electromigration process. We use a combination of experimental and computational approaches to investigate the role of grain morphology and distribution in granular aluminium thin films, when formed into nanowire constrictions. Treating the granular aluminium film as a network of randomly distributed resistors with parameters motivated by the film microstructure allows us to model the electrical characteristics of the nanowires. This model provides estimates of the dependence of sheet resistance on grain size and distribution, and the resulting device to device variation for superconducting nanowires. By fabricating a series of different length nanowires, we study the electromigration process as a function of applied current, and then compare directly to the results of our computational model. In doing so we show that the electromigration is driven by the formation of quantum point contacts between metallic aluminium grains. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 11 pages, 11 figures

Structure and vibrational properties of methane up to 71 GPa

Authors: Maxim Bykov, Elena Bykova, Chris J. Pickard, Miguel Martinez-Canales, Konstantin Glazyrin, Jesse S. Smith, Alexander F. Goncharov

Abstract: Single-crystal synchrotron X-ray diffraction, Raman spectroscopy, and first principles calculations have been used to identify the structure of the high-pressure (HP) phase of molecular methane above 20 GPa up to 71 GPa. The structure of HP phase is trigonal R3, which can be represented as a distortion of the cubic phase B, previously documented at 7-15 GPa and confirmed here. The positions of hyd… ▽ More Single-crystal synchrotron X-ray diffraction, Raman spectroscopy, and first principles calculations have been used to identify the structure of the high-pressure (HP) phase of molecular methane above 20 GPa up to 71 GPa. The structure of HP phase is trigonal R3, which can be represented as a distortion of the cubic phase B, previously documented at 7-15 GPa and confirmed here. The positions of hydrogen atoms in HP phase have been 19 obtained from first principles calculations. The molecules occupy four different crystallographic sites in phases B and eleven sites in the HP phase, which result in splitting of molecular stretching modes detected in Raman spectroscopy and assigned here based on 22 a good agreement with the Raman spectra calculated from the first principles. △ Less

Submitted 14 October, 2021; originally announced October 2021.

Comments: 29 pages, 7 figures, 3 tables, supplemental material

Squirmer hydrodynamics near a periodic surface topography

Authors: Kenta Ishimoto, Eamonn A. Gaffney, David J. Smith

Abstract: The behaviour of microscopic swimmers has previously been explored near large scale confining geometries and in the presence of very small-scale surface roughness. Here we consider an intermediate case of how a simple microswimmer, the tangential spherical squirmer, behaves adjacent to singly- and doubly-periodic sinusoidal surface topographies that spatially oscillate with an amplitude that is an… ▽ More The behaviour of microscopic swimmers has previously been explored near large scale confining geometries and in the presence of very small-scale surface roughness. Here we consider an intermediate case of how a simple microswimmer, the tangential spherical squirmer, behaves adjacent to singly- and doubly-periodic sinusoidal surface topographies that spatially oscillate with an amplitude that is an order of magnitude less than the swimmer size and wavelengths within an order of magnitude of this scale. The nearest neighbour regularised Stokeslet method is used for numerical explorations after validating its accuracy for a spherical tangential squirmer that swims stably near a flat surface. The same squirmer is then introduced to the different surface topographies. The key governing factor in the resulting swimming behaviour is the size of the squirmer relative to the surface topography wavelength. When the squirmer is much larger, or much smaller, than the surface topography wavelength then directional guidance is not observed. Once the squirmer size is on the scale of the topography wavelength limited guidance is possible, often with local capture in the topography troughs. However, complex dynamics can also emerge, especially when the initial configuration is not close to alignment along topography troughs or above topography crests. In contrast to sensitivity in alignment and the topography wavelength, reductions in the amplitude of the surface topography or variations in the shape of the periodic surface topography do not have extensive impacts on the squirmer behaviour. Our findings more generally highlight that detailed numerical explorations are necessary to determine how surface topographies may interact with a given swimmer that swims stably near a flat surface and whether surface topographies may be designed to provide such swimmers with directional guidance. △ Less

Submitted 12 September, 2021; originally announced September 2021.

Comments: 13 pages, 9 figures

Data Mining for Terahertz Generation Crystals

Authors: Gabriel A. Valdivia-Berroeta, Zachary B. Zaccardi, Sydney K. F. Pettit, Sin-Hang Ho, Bruce Wayne Palmer, Matthew J. Lutz, Claire Rader, Brittan P. Hunter, Natalie K. Green, Connor Barlow, Coriantumr Z. Wayment, Daisy J. Harmon, Paige Petersen, Stacey J. Smith, David J. Michaelis, Jeremy A. Johnson

Abstract: We demonstrate a data mining approach to discover and develop new organic nonlinear optical crystals that produce intense pulses of terahertz radiation. We mine the Cambridge Structural Database for non-centrosymmetric materials and use this structural data in tandem with density functional theory calculations to predict new materials that efficiently generate terahertz radiation. This enables us… ▽ More We demonstrate a data mining approach to discover and develop new organic nonlinear optical crystals that produce intense pulses of terahertz radiation. We mine the Cambridge Structural Database for non-centrosymmetric materials and use this structural data in tandem with density functional theory calculations to predict new materials that efficiently generate terahertz radiation. This enables us to (in a relatively short time) discover, synthesize, and grow large, high-quality crystals of four promising materials and characterize them for intense terahertz generation. In a direct comparison to the current state-of-the-art organic terahertz generation crystals, these new materials excel. The discovery and characterization of these novel terahertz generators validates the approach of combining data mining with density functional theory calculations to predict properties of high-performance organic materials, potentially for a host of exciting applications. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 16 pages, 5 figures

Pressure-induced insulator-metal transition in two-dimensional Mott insulator NiPS3

Authors: Takahiro Matsuoka, Amanda Haglund, Rui Xue, Jesse S. Smith, Maik Lang, Antonio M. dos Santos, David Mandrus

Abstract: The pressure-induced insulator to metal transition (IMT) of layered magnetic nickel phosphorous tri-sulfide NiPS3 was studied in-situ under quasi-uniaxial conditions by means of electrical resistance (R) and X-ray diffraction (XRD) measurements. This sluggish transition is shown to occur at 35 GPa. Transport measurements show no evidence of superconductivity to the lowest measured temperature (~ 2… ▽ More The pressure-induced insulator to metal transition (IMT) of layered magnetic nickel phosphorous tri-sulfide NiPS3 was studied in-situ under quasi-uniaxial conditions by means of electrical resistance (R) and X-ray diffraction (XRD) measurements. This sluggish transition is shown to occur at 35 GPa. Transport measurements show no evidence of superconductivity to the lowest measured temperature (~ 2 K). The structure results presented here differ from earlier in-situ work that subjected the sample to a different pressure state, suggesting that in NiPS3 the phase stability fields are highly dependent on strain. It is suggested that careful control of the strain is essential when studying the electronic and magnetic properties of layered van der Waals solids. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Sub-Unit Cell Spatial Resolution of Phase and Polarization Mapping in Scanning Electron Diffraction

Authors: Tim B. Eldred, Jacob G. Smith, Wenpei Gao

Abstract: Diffraction analysis in four dimensional scanning transmission electron microscopy now enables the mapping of local structures including symmetry, strain, and polarization of materials. However, measuring the distribution of these configurations at the unit cell level remains a challenge because most analysis methods require the diffraction disks to be separated, limiting the electron probe sizes… ▽ More Diffraction analysis in four dimensional scanning transmission electron microscopy now enables the mapping of local structures including symmetry, strain, and polarization of materials. However, measuring the distribution of these configurations at the unit cell level remains a challenge because most analysis methods require the diffraction disks to be separated, limiting the electron probe sizes to be larger than a unit cell. Here we show improved spatial resolution in mapping the polarization displacement and phases of BaTiO3 sampled at a rate equivalent to the size of the projected unit cells using 4DSTEM. This improvement in spatial resolution is accomplished by masking out the overlapping regions in partially overlapped CBED patterns. By reducing the probe size to the order of single projected unit cells in size, the fluctuations of symmetry and structure inside the rhombohedral symmetry nanoregions inside of tetragonal domains can be observed. With the mapping performed at the unit cell resolution the nano-sized rhombohedral domains break down with local fluctuation, showing different character at different length scales. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Approximate theories for binary magnetic quantum droplets

Authors: Joseph C. Smith, P. B. Blakie, D. Baillie

Abstract: We develop two approximate theories to describe the miscible and immiscible droplets that can occur in a binary mixture of highly magnetic bosonic atoms. In addition to allowing simpler calculations, the approximate theories provide insight into the role of quantum fluctuations in the two regimes. Results are validated by comparison to those from the extended Gross-Pitaevskii equation. As an appli… ▽ More We develop two approximate theories to describe the miscible and immiscible droplets that can occur in a binary mixture of highly magnetic bosonic atoms. In addition to allowing simpler calculations, the approximate theories provide insight into the role of quantum fluctuations in the two regimes. Results are validated by comparison to those from the extended Gross-Pitaevskii equation. As an application we solve for the ground state droplets crossing the miscible-immiscible transition as function of the short-ranged interspecies interaction parameter. We consider regimes where the transition occurs suddenly or as a smooth cross-over. Using dynamical calculations we show that the character of the transition is revealed in the number of domains produced when ramping the droplet into the immiscible regime. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: 10 pages, 6 figures

Journal ref: Phys. Rev. A 104, 053316 (2021)

Imaging the breaking of electrostatic dams in graphene for ballistic and viscous fluids

Authors: Zachary J. Krebs, Wyatt A. Behn, Songci Li, Keenan J. Smith, Kenji Watanabe, Takashi Taniguchi, Alex Levchenko, Victor W. Brar

Abstract: The flow of charge carriers in materials can, under some circumstances, mimic the flow of viscous fluids. In order to visualize the consequences of such effects, new methodologies must be developed that can probe the quasiparticle flow profile with nm-scale resolution as the geometric parameters of the system are continuously evolved. In this work, scanning tunneling potentiometry (STP) is used to… ▽ More The flow of charge carriers in materials can, under some circumstances, mimic the flow of viscous fluids. In order to visualize the consequences of such effects, new methodologies must be developed that can probe the quasiparticle flow profile with nm-scale resolution as the geometric parameters of the system are continuously evolved. In this work, scanning tunneling potentiometry (STP) is used to image quasiparticle flow around engineered electrostatic barriers in graphene/hBN heterostructures. Measurements are performed as electrostatic dams - defined by lateral pn-junction barriers - are broken within the graphene sheet, and carriers move through conduction channels with physical widths that vary continuously from pinch-off to um-scale. Local, STP measurements of the electrochemical potential allow for direct characterization of the evolving flow profile, which we compare to finite-element simulations of a Stokesian fluid with varying parameters. Our results reveal distinctly non-Ohmic flow profiles, with charge dipoles forming across barriers due to carrier scattering and accumulation on the upstream side, and depletion downstream. Conductance measurements of individual channels, meanwhile, reveal that at low temperatures the quasiparticle flow is ballistic, but as the temperature is raised there is a Knudsen-to-Gurzhi regime crossover where the fluid becomes viscous and the channel conductance exceeds the ballistic limit set by Sharvin conductance. These results provide a clear illustration of how carrier flow in a Fermi fluid evolves as a function of carrier density, channel width, and temperature. They also demonstrate how STP can be used to extract key parameters of quasiparticle transport, with a spatial resolution that exceeds that of other methods by orders of magnitude. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Journal ref: Science 379, 671-676 (2023)

Microscopic processes during ultra-fast laser generation of Frenkel defects in diamond

Authors: Benjamin Griffiths, Andrew Kirkpatrick, Shannon S. Nicley, Rajesh L. Patel, Joanna M. Zajac, Gavin W. Morley, Martin J. Booth, Patrick S. Salter, Jason M. Smith

Abstract: Engineering single atomic defects into wide bandgap materials has become an attractive field in recent years due to emerging applications such as solid-state quantum bits and sensors. The simplest atomic-scale defect is the lattice vacancy which is often a constituent part of more complex defects such as the nitrogen-vacancy (NV) centre in diamond, therefore an understanding of the formation mecha… ▽ More Engineering single atomic defects into wide bandgap materials has become an attractive field in recent years due to emerging applications such as solid-state quantum bits and sensors. The simplest atomic-scale defect is the lattice vacancy which is often a constituent part of more complex defects such as the nitrogen-vacancy (NV) centre in diamond, therefore an understanding of the formation mechanisms and precision engineering of vacancies is desirable. We present a theoretical and experimental study into the ultra-fast laser generation of vacancy-interstitial pairs (Frenkel defects) in diamond. The process is described by a set of coupled rate equations of the pulsed laser interaction with the material and of the non-equilibrium dynamics of charge carriers during and in the wake of the pulse. We find that a model for Frenkel defect generation via the recombination of a bound biexciton as the electron plasma cools provides good agreement with experimental data, reproducing an effective non-linearity of $\sim$ 40 for Frenkel defect generation with respect to laser pulse energy. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: 17 pages, 11 figures

Influence of device geometry and imperfections on the interpretation of transverse magnetic focusing experiments

Authors: Yik K. Lee, Jackson S. Smith, Jared H. Cole

Abstract: Spatially separating electrons of different spins and efficiently generating spin currents are crucial steps towards building practical spintronics devices. Transverse magnetic focusing is a potential technique to accomplish both those tasks. In a material where there is significant Rashba spin--orbit interaction, electrons of different spins will traverse different paths in the presence of an ext… ▽ More Spatially separating electrons of different spins and efficiently generating spin currents are crucial steps towards building practical spintronics devices. Transverse magnetic focusing is a potential technique to accomplish both those tasks. In a material where there is significant Rashba spin--orbit interaction, electrons of different spins will traverse different paths in the presence of an external magnetic field. Experiments have demonstrated the viability of this technique by measuring conductance spectra that indicate the separation of spin-up and spin-down electrons. However the effect that the geometry of the leads has on these measurements is not well understood. We show that the resolution of features in the conductance spectra is affected by the shape, separation and width of the leads. Furthermore, the number of subbands occupied by the electrons in the leads affects the ratio between the amplitudes of the spin-split peaks in the spectra. We simulated devices with random onsite potentials and observed that transverse magnetic focusing devices are sensitive to disorder. Ultimately we show that careful choice and characterisation of device geometry is crucial for correctly interpreting the results of transverse magnetic focusing experiments. △ Less

Submitted 18 March, 2022; v1 submitted 26 April, 2021; originally announced April 2021.

Comments: 13 pages, 13 figures. Final version with minor revisions from previous version

Long Phonon Mean Free Paths Observed in Cross-plane Thermal-Conductivity Measurements of Exfoliated Hexagonal Boron Nitride

Authors: Gabriel R. Jaffe, Keenan J. Smith, Kenji Watanabe, Takashi Taniguchi, Max G. Lagally, Mark A. Eriksson, Victor W. Brar

Abstract: Sub-micron-thick layers of hexagonal boron nitride (hBN) exhibit high in-plane thermal conductivity, useful optical properties, and serve as dielectric encapsulation layers with low electrostatic inhomogeneity for graphene devices. Despite the promising applications of hBN as a heat spreader, the thickness dependence of the cross-plane thermal conductivity is not known, and the cross-plane phonon… ▽ More Sub-micron-thick layers of hexagonal boron nitride (hBN) exhibit high in-plane thermal conductivity, useful optical properties, and serve as dielectric encapsulation layers with low electrostatic inhomogeneity for graphene devices. Despite the promising applications of hBN as a heat spreader, the thickness dependence of the cross-plane thermal conductivity is not known, and the cross-plane phonon mean free paths in hBN have not been measured. We measure the cross-plane thermal conductivity of hBN flakes exfoliated from bulk crystals. We find that the thermal conductivity is extremely sensitive to film thickness. We measure a forty-fold increase in the cross-plane thermal conductivity between 7 nm and 585 nm flakes at 295 K. We attribute the large increase in thermal conductivity with increasing thickness to contributions from phonons with long mean free paths (MFPs), spanning many hundreds of nanometers in the thickest flakes. When planar twist interfaces are introduced into the crystal by mechanically stacking multiple thin flakes, the cross-plane thermal conductivity of the stack is found to be a factor of seven below that of individual flakes with similar total thickness, thus providing strong evidence that phonon scattering at twist boundaries limits the maximum phonon MFPs. These results have important implications for hBN integration in nanoelectronics and improve our understanding of thermal transport in two-dimensional materials. △ Less

Submitted 8 October, 2021; v1 submitted 12 March, 2021; originally announced March 2021.

Comments: 4 pages, 3 figures

Journal ref: ACS Appl. Mater. Interfaces 2023, 15, 9, 12545-12550

Localized Wannier function based tight-binding models for two-dimensional allotropes of bismuth

Authors: Qile Li, Jackson S. Smith, Yuefeng Yin, Chutian Wang, Mykhailo V. Klymenko, Jared H. Cole, Nikhil V. Medhekar

Abstract: With its monoelemental composition, various crystalline forms and an inherently strong spin-orbit coupling, bismuth has been regarded as an ideal prototype material to expand our understanding of topological electronic structures. In particular, two-dimensional bismuth thin films have attracted a growing interest due to potential applications in topological transistors and spintronics. This calls… ▽ More With its monoelemental composition, various crystalline forms and an inherently strong spin-orbit coupling, bismuth has been regarded as an ideal prototype material to expand our understanding of topological electronic structures. In particular, two-dimensional bismuth thin films have attracted a growing interest due to potential applications in topological transistors and spintronics. This calls for an effective physical model to give an accurate interpretation of the novel topological phenomena shown by two-dimensional bismuth. However, the conventional semi-empirical approach of adapting bulk bismuth hoppings fails to capture the topological features of two-dimensional bismuth allotropes because the electronic band topology is heavily influenced by crystalline symmetries as well as atom spacings. Here we provide a new parameterization using localized Wannier functions derived from the Bloch states in first-principles calculations. We construct new tight-binding models for three types of two-dimensional bismuth allotropes: a Bi (111) bilayer, bismuthene and a Bi(110) bilayer. We demonstrate that our tight-binding models can successfully reproduce the band structures, symmetries and topological features of these two-dimensional allotropes. We anticipate that these models can be extended to other similar two-dimensional topological structures such as antimonene and arsenene. Moreover, these models can serve as a starting point for investigating the electron/spin transport and electromagnetic response in low-dimensional topological devices. △ Less

Submitted 22 February, 2021; originally announced February 2021.

Comments: 20 pages, 6 figures

Journal ref: New J. Phys. 23, 063042 (2021)

Ubiquitous Dynamical Time Asymmetry in Measurements on Materials and Biological Systems

Authors: Alessio Lapolla, Jeremy C. Smith, Aljaž Godec

Abstract: Many measurements on soft condensed matter (e.g., biological and materials) systems track low-dimensional observables projected from the full system phase space as a function of time. Examples are dynamic structure factors, spectroscopic and rheological response functions, and time series of distances derived from optical tweezers, single-molecule spectroscopy and molecular dynamics simulations. I… ▽ More Many measurements on soft condensed matter (e.g., biological and materials) systems track low-dimensional observables projected from the full system phase space as a function of time. Examples are dynamic structure factors, spectroscopic and rheological response functions, and time series of distances derived from optical tweezers, single-molecule spectroscopy and molecular dynamics simulations. In many such systems the projection renders the reduced dynamics non-Markovian and the observable is not prepared in, or initially sampled from and averaged over, a stationary distribution. We prove that such systems always exhibit non-equilibrium, time asymmetric dynamics. That is, they evolve in time with a broken time-translation invariance in a manner closely resembling aging dynamics. We identify the entropy associated with the breaking of time-translation symmetry that is a measure of the instantaneous thermodynamic displacement of latent, hidden degrees of freedom from their stationary state. Dynamical time asymmetry is a general phenomenon, independent of the underlying energy surface, and is frequently even visible in measurements on systems that have fully reached equilibrium. This finding has fundamental implications for the interpretation of many experiments on, and simulations of, biological and materials systems. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Mixed Precision Fermi-Operator Expansion on Tensor Cores From a Machine Learning Perspective

Authors: Joshua Finkelstein, Justin Smith, Susan M. Mniszewski, Kipton Barros, Christian F. A. Negre, Emanuel H. Rubensson, Anders M. N. Niklasson

Abstract: We present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia's A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generali… ▽ More We present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia's A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent fractional occupation numbers of the electronic states at finite temperatures. △ Less

Submitted 16 January, 2021; originally announced January 2021.

Comments: 10 pages, 6 figures

Report number: LA-UR-21-20350

Measurements of Optical Scatter Versus Annealing Temperature for Amorphous Ta2O5 and TiO2:Ta2O5 Thin Films

Authors: Elenna M. Capote, Amy Gleckl, Jazlyn Guerrero, Michael Rezac, Robert Wright, Joshua R. Smith

Abstract: Optical coatings formed from amorphous oxide thin films have many applications in precision measurements. The Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) and Advanced Virgo use coatings of SiO$_2$ (silica) and TiO$_2$:Ta$_2$O$_5$ (titania doped tantala) and post-deposition annealing to 500$^\circ$C to achieve low thermal noise and low optical absorption. Optical scattering… ▽ More Optical coatings formed from amorphous oxide thin films have many applications in precision measurements. The Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) and Advanced Virgo use coatings of SiO$_2$ (silica) and TiO$_2$:Ta$_2$O$_5$ (titania doped tantala) and post-deposition annealing to 500$^\circ$C to achieve low thermal noise and low optical absorption. Optical scattering by these coatings is a key limit to the detectors' sensitivity. This paper describes optical scattering measurements for single-layer ion-beam-sputtered thin films on fused silica substrates: two samples of Ta$_2$O$_5$ and two of TiO$_2$:Ta$_2$O$_5$. Using an imaging scatterometer at a fixed scattering angle of 12.8$^\circ$, in-situ changes in the optical scatter of each sample were assessed during post-deposition annealing to 500$^\circ$C in vacuum. The scatter of three of the four coated optics was observed to decrease during the annealing process, by 25-30$\%$ for tantala and up to 74$\%$ for titania-doped tantala, while scatter from the fourth sample held constant. Angle-resolved scatter measurements performed before and after vacuum annealing suggests some improvement in three of four samples. These results demonstrate that post-deposition high-temperature annealing of single-layer tantala and titania-doped tantala thin films in vacuum does not lead to an increase in scatter, and may actually improve their scatter. △ Less

Submitted 28 April, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: 12 pages, 7 figures, research paper

Report number: LIGO Document P2000483