Showing 1–50 of 173 results for author: Müller, E

1D $Z_2$ lattice gauge theory in periodic Gauss law sectors

Authors: Vaibhav Sharma, Erich J Mueller

Abstract: We calculate the properties of a 1D $Z_2$ lattice gauge theory in different Gauss law sectors, corresponding to different configurations of static charges set by the orientations of the gauge spins. Importantly, in quantum simulator experiments these sectors can be accessed without adding any additional physical particles or changing the Hamiltonian: The Gauss law sectors are simply set by the ini… ▽ More We calculate the properties of a 1D $Z_2$ lattice gauge theory in different Gauss law sectors, corresponding to different configurations of static charges set by the orientations of the gauge spins. Importantly, in quantum simulator experiments these sectors can be accessed without adding any additional physical particles or changing the Hamiltonian: The Gauss law sectors are simply set by the initial conditions. We study the interplay between conservation laws and interactions when the static charges are chosen to form periodic patterns. We classify the different Gauss law sectors and use the density matrix renormalization group to calculate the ground state compressibility, density profiles, charge density wave order parameters, and single particle correlation functions as a function of matter density. We find confined and deconfined phases, charge density waves, correlated insulators, and supersolids. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 11 pages, 11 figures

Defects induce phase transition from dynamic to static rippling in graphene

Authors: Fabian L. Thiemann, Camille Scalliet, Erich A. Müller, Angelos Michaelides

Abstract: Many of graphene's remarkable properties are intrinsically linked to its inherent ripples. Defects, whether naturally present or artificially introduced, are known to have a strong impact on the rippling of graphene. However, how defects alter ripple dynamics in two-dimensional (2D) materials in general, and graphene in particular, remains largely unexplored. Here, using machine learning-driven mo… ▽ More Many of graphene's remarkable properties are intrinsically linked to its inherent ripples. Defects, whether naturally present or artificially introduced, are known to have a strong impact on the rippling of graphene. However, how defects alter ripple dynamics in two-dimensional (2D) materials in general, and graphene in particular, remains largely unexplored. Here, using machine learning-driven molecular dynamics simulations, we reveal a fundamental connection between defect concentration and ripple dynamics in freestanding graphene sheets. Specifically, we find that at a critical concentration of approximately 0.1%, dynamic rippling undergoes a transition from freely propagating to static ripples. This is in quantitative alignment with the experimentally observed turning point in the non-monotonic scaling of the Young's modulus and emphasises the critical interplay between defects and material dynamics. Our work not only unveils the significant impact of defects on rippling dynamics in graphene but also paves the way to design two-dimensional devices with tailored properties. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 8 pages, 20 pages supplementary information

Simulating X-ray absorption spectroscopy of battery materials on a quantum computer

Authors: Stepan Fomichev, Kasra Hejazi, Ignacio Loaiza, Modjtaba Shokrian Zini, Alain Delgado, Arne-Christian Voigt, Jonathan E. Mueller, Juan Miguel Arrazola

Abstract: X-ray absorption spectroscopy is a crucial experimental technique for elucidating the mechanisms of structural degradation in battery materials. However, extracting information from the measured spectrum is challenging without high-quality simulations. In this work, we propose simulating near-edge X-ray absorption spectra as a promising application for quantum computing. It is attractive due to th… ▽ More X-ray absorption spectroscopy is a crucial experimental technique for elucidating the mechanisms of structural degradation in battery materials. However, extracting information from the measured spectrum is challenging without high-quality simulations. In this work, we propose simulating near-edge X-ray absorption spectra as a promising application for quantum computing. It is attractive due to the ultralocal nature of X-ray absorption that significantly reduces the sizes of problems to be simulated, and because of the classical hardness of simulating spectra. We describe three quantum algorithms to compute the X-ray absorption spectrum and provide their asymptotic cost. One of these is a Monte-Carlo based time-domain algorithm, which is cost-friendly to early fault-tolerant quantum computers. We then apply the framework to an industrially relevant example, a CAS(22e,18o) active space for an O-Mn cluster in a Li-excess battery cathode, showing that practically useful simulations could be obtained with much fewer qubits and gates than ground-state energy estimation of the same material. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Dynamics of spin helices in the diluted one-dimensional $XX$ model

Authors: Darren Pereira, Erich J. Mueller

Abstract: Motivated by discrepancies between recent cold atom experiments and the associated theory, we explore the effect of immobile holes on the quantum dynamics of $x$-$z$ spin helices in the one-dimensional $XX$ model. We calculate the exact spin dynamics by mapping onto a system of non-interacting fermions, averaging over the distribution of holes. At small hole densities we find that the helical spin… ▽ More Motivated by discrepancies between recent cold atom experiments and the associated theory, we explore the effect of immobile holes on the quantum dynamics of $x$-$z$ spin helices in the one-dimensional $XX$ model. We calculate the exact spin dynamics by mapping onto a system of non-interacting fermions, averaging over the distribution of holes. At small hole densities we find that the helical spin pattern decays exponentially, with a pitch dependence that agrees with the experiments. At large hole densities we instead find persistent oscillations. While our analytic approach does not generalize to the $XXZ$ model with arbitrary anisotropies, we validate a matrix product state technique which might be used to model the experiments in those settings. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 9 pages, 8 figures

High temperature transport in the one dimensional mass-imbalanced Fermi-Hubbard model

Authors: Thomas G. Kiely, Erich J. Mueller

Abstract: We study transport in the one-dimensional mass-imbalanced Fermi-Hubbard model at infinite temperature, focusing on the case of strong interactions. Prior theoretical and experimental investigations have revealed unconventionally long transport timescales, with complications due to strong finite size effects. We compute the dynamical current-current correlation function directly in the thermodynami… ▽ More We study transport in the one-dimensional mass-imbalanced Fermi-Hubbard model at infinite temperature, focusing on the case of strong interactions. Prior theoretical and experimental investigations have revealed unconventionally long transport timescales, with complications due to strong finite size effects. We compute the dynamical current-current correlation function directly in the thermodynamic limit using infinite tensor network techniques. We show that transport in the strong-imbalance limit is dominated by AC resonances, which we compute with an analytic expansion. We study the dephasing of these resonances with mass imbalance, $η$. In the small-imbalance limit, the model is nearly integrable. We connect these unusual limits by computing the DC conductivity and transport decay time as a function of $η$ and the interaction strength $U/t$. We propose an experimental protocol to measure these correlation functions in cold atom experiments. △ Less

Submitted 31 May, 2024; v1 submitted 11 April, 2024; originally announced April 2024.

Comments: 13 pages, 6 figures

Coherent Acoustic Control of Defect Orbital States in the Strong-Driving Limit

Authors: B. A. McCullian, V. Sharma, H. Y. Chen, J. C. Crossman, E. J. Mueller, G. D. Fuchs

Abstract: We use a bulk acoustic wave resonator to demonstrate coherent control of the excited orbital states in a diamond nitrogen-vacancy (NV) center at cryogenic temperature. Coherent quantum control is an essential tool for understanding and mitigating decoherence. Moreover, characterizing and controlling orbital states is a central challenge for quantum networking, where optical coherence is tied to or… ▽ More We use a bulk acoustic wave resonator to demonstrate coherent control of the excited orbital states in a diamond nitrogen-vacancy (NV) center at cryogenic temperature. Coherent quantum control is an essential tool for understanding and mitigating decoherence. Moreover, characterizing and controlling orbital states is a central challenge for quantum networking, where optical coherence is tied to orbital coherence. We study resonant multi-phonon orbital Rabi oscillations in both the frequency and time domain, extracting the strength of the orbital-phonon interactions and the coherence of the acoustically driven orbital states. We reach the strong-driving limit, where the physics is dominated by the coupling induced by the acoustic waves. We find agreement between our measurements, quantum master equation simulations, and a Landau-Zener transition model in the strong-driving limit. Using perturbation theory, we derive an expression for the orbital Rabi frequency versus acoustic drive strength that is non-perturbative in the drive strength and agrees well with our measurements for all acoustic powers. Motivated by continuous wave spin resonance-based decoherence protection schemes, we model the orbital decoherence and find good agreement between our model and our measured few-to-several nanoseconds orbital decoherence times. We discuss the outlook for orbital decoherence protection. △ Less

Submitted 16 March, 2024; originally announced March 2024.

Comments: 19 pages, 11 figures including supplement

Accurate Crystal Structure Prediction of New 2D Hybrid Organic Inorganic Perovskites

Authors: Nima Karimitari, William J. Baldwin, Evan W. Muller, Zachary J. L. Bare, W. Joshua Kennedy, Gábor Csányi, Christopher Sutton

Abstract: Low dimensional hybrid organic-inorganic perovskites (HOIPs) represent a promising class of electronically active materials for both light absorption and emission. The design space of HOIPs is extremely large, since a diverse space of organic cations can be combined with different inorganic frameworks. This immense design space allows for tunable electronic and mechanical properties, but also nece… ▽ More Low dimensional hybrid organic-inorganic perovskites (HOIPs) represent a promising class of electronically active materials for both light absorption and emission. The design space of HOIPs is extremely large, since a diverse space of organic cations can be combined with different inorganic frameworks. This immense design space allows for tunable electronic and mechanical properties, but also necessitates the development of new tools for in silico high throughput analysis of candidate structures. In this work, we present an accurate, efficient, transferable and widely applicable machine learning interatomic potential (MLIP) for predicting the structure of new 2D HOIPs. Using the MACE architecture, an MLIP is trained on 86 diverse experimentally reported HOIP structures. The model is tested on 73 unseen perovskite compositions, and achieves chemical accuracy with respect to the reference electronic structure method. Our model is then combined with a simple random structure search algorithm to predict the structure of hypothetical HOIPs given only the proposed composition. Success is demonstrated by correctly and reliably recovering the crystal structure of a set of experimentally known 2D perovskites. Such a random structure search is impossible with ab initio methods due to the associated computational cost, but is relatively inexpensive with the MACE potential. Finally, the procedure is used to predict the structure formed by a new organic cation with no previously known corresponding perovskite. Laboratory synthesis of the new hybrid perovskite confirms the accuracy of our prediction. This capability, applied at scale, enables efficient screening of thousands of combinations of organic cations and inorganic layers. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 14 pages and 9 figures in the main text. Supplementary included in pdf

Emergence of interfacial magnetism in strongly-correlated nickelate-titanate superlattices

Authors: Teguh Citra Asmara, Robert J. Green, Andreas Suter, Yuan Wei, Wenliang Zhang, Grant Harris, Yi Tseng, Tianlun Yu, Davide Betto, Mirian Garcia-Fernandez, Stefano Agrestini, Yannick Maximilian Klein, Neeraj Kumar, Carlos William Galdino, Zaher Salman, Thomas Prokscha, Marisa Medarde, Elisabeth Müller, Yona Soh, Nicholas B. Brookes, Ke-Jin Zhou, Milan Radovic, Thorsten Schmitt

Abstract: Strongly-correlated transition-metal oxides are widely known for their various exotic phenomena. This is exemplified by rare-earth nickelates such as LaNiO$_{3}$, which possess intimate interconnections between their electronic, spin, and lattice degrees of freedom. Their properties can be further enhanced by pairing them in hybrid heterostructures, which can lead to hidden phases and emergent phe… ▽ More Strongly-correlated transition-metal oxides are widely known for their various exotic phenomena. This is exemplified by rare-earth nickelates such as LaNiO$_{3}$, which possess intimate interconnections between their electronic, spin, and lattice degrees of freedom. Their properties can be further enhanced by pairing them in hybrid heterostructures, which can lead to hidden phases and emergent phenomena. An important example is the LaNiO$_{3}$/LaTiO$_{3}$ superlattice, where an interlayer electron transfer has been observed from LaTiO$_{3}$ into LaNiO$_{3}$ and is predicted to result in a high-spin state. However, macroscopic emergence of magnetic order has so far not been observed. Here, by using muon spin rotation, x-ray absorption, and resonant inelastic x-ray scattering, we present direct evidence of an emergent antiferromagnetic order with high magnon energy and exchange interactions at the LaNiO$_{3}$/LaTiO$_{3}$ interface. As the magnetism is purely interfacial, a single LaNiO$_{3}$/LaTiO$_{3}$ interface can essentially behave as an atomically thin quasi-two-dimensional antiferromagnet, potentially allowing its technological utilisation in advanced spintronic devices. Furthermore, its strong quasi-two-dimensional magnetic correlations and orbitally-polarized planar ligand holes make its electronic and magnetic configurations resemble the precursor states of superconducting cuprates and nickelates, but with an S $\rightarrow$ 1 spin state instead. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 41 pages, 13 figures

Bose-enhanced relaxation of driven atom-molecule condensates

Authors: Dimitri Pimenov, Erich J. Mueller

Abstract: Motivated by recent experiments we study the interconversion between ultracold atomic and molecular condensates, quantifying the resulting oscillations and their slow decay. We find that near equilibrium the dominant damping source is the decay of condensed molecules into non-condensed pairs, with a pair kinetic energy that is resonant with the frequency of the oscillating atom-molecule interconve… ▽ More Motivated by recent experiments we study the interconversion between ultracold atomic and molecular condensates, quantifying the resulting oscillations and their slow decay. We find that near equilibrium the dominant damping source is the decay of condensed molecules into non-condensed pairs, with a pair kinetic energy that is resonant with the frequency of the oscillating atom-molecule interconversions. The decay, however, is non-exponential, as strong population of the resonant pairs leads to Bose enhancement. Introducing an oscillating magnetic field, which periodically modulates the molecular binding energy, enhances the oscillations at short times. However, the resulting enhancement of the pair-production process results in an accelerated decay which rapidly cuts off the initial oscillation growth. △ Less

Submitted 19 April, 2024; v1 submitted 5 January, 2024; originally announced January 2024.

Comments: journal version

Journal ref: Phys. Rev. A 109, 043311 (2024)

Superconducting Penetration Depth Through a Van Hove Singularity: Sr$_2$RuO$_4$ Under Uniaxial Stress

Authors: Eli Mueller, Yusuke Iguchi, Fabian Jerzembeck, Jorge O. Rodriguez, Marisa Romanelli, Edgar Abarca-Morales, Anastasios Markou, Naoki Kikugawa, Dmitry A. Sokolov, Gwansuk Oh, Clifford W. Hicks, Andrew P. Mackenzie, Yoshiteru Maeno, Vidya Madhavan, Kathryn A. Moler

Abstract: In the unconventional superconductor Sr$_2$RuO$_4$, uniaxial stress along the $[100]$ direction tunes the Fermi level through a Van Hove singularity (VHS) in the density of states, causing a strong enhancement of the superconducting critical temperature $T_\textrm{c}$. Here, we report measurements of the London penetration depth $λ$ as this tuning is performed. We find that the zero-temperature su… ▽ More In the unconventional superconductor Sr$_2$RuO$_4$, uniaxial stress along the $[100]$ direction tunes the Fermi level through a Van Hove singularity (VHS) in the density of states, causing a strong enhancement of the superconducting critical temperature $T_\textrm{c}$. Here, we report measurements of the London penetration depth $λ$ as this tuning is performed. We find that the zero-temperature superfluid density, here defined as $λ(0)^{-2}$, increases by $\sim$15%, with a peak that coincides with the peak in $T_\textrm{c}$. We also find that the low temperature form of $λ(T)$ is quadratic over the entire strain range. Using scanning tunneling microscopy, we find that the gap increases from $Δ_0 \approx 350~μ$eV in unstressed Sr$_2$RuO$_4$ to $Δ_0 \approx 600~μ$eV in a sample strained to near the peak in $T_c$. With a nodal order parameter, an increase in the superconducting gap could bring about an increase in the superfluid density through reduced sensitivity to defects and through reduced non-local effects in the Meissner screening. Our data indicate that tuning to the VHS increases the gap throughout the Brillouin zone, and that non-local effects are likely more important than reduced scattering. △ Less

Submitted 11 December, 2023; v1 submitted 8 December, 2023; originally announced December 2023.

Initial state preparation for quantum chemistry on quantum computers

Authors: Stepan Fomichev, Kasra Hejazi, Modjtaba Shokrian Zini, Matthew Kiser, Joana Fraxanet Morales, Pablo Antonio Moreno Casares, Alain Delgado, Joonsuk Huh, Arne-Christian Voigt, Jonathan E. Mueller, Juan Miguel Arrazola

Abstract: Quantum algorithms for ground-state energy estimation of chemical systems require a high-quality initial state. However, initial state preparation is commonly either neglected entirely, or assumed to be solved by a simple product state like Hartree-Fock. Even if a nontrivial state is prepared, strong correlations render ground state overlap inadequate for quality assessment. In this work, we addre… ▽ More Quantum algorithms for ground-state energy estimation of chemical systems require a high-quality initial state. However, initial state preparation is commonly either neglected entirely, or assumed to be solved by a simple product state like Hartree-Fock. Even if a nontrivial state is prepared, strong correlations render ground state overlap inadequate for quality assessment. In this work, we address the initial state preparation problem with an end-to-end algorithm that prepares and quantifies the quality of initial states, accomplishing the latter with a new metric -- the energy distribution. To be able to prepare more complicated initial states, we introduce an implementation technique for states in the form of a sum of Slater determinants that exhibits significantly better scaling than all prior approaches. We also propose low-precision quantum phase estimation (QPE) for further state quality refinement. The complete algorithm is capable of generating high-quality states for energy estimation, and is shown in select cases to lower the overall estimation cost by several orders of magnitude when compared with the best single product state ansatz. More broadly, the energy distribution picture suggests that the goal of QPE should be reinterpreted as generating improvements compared to the energy of the initial state and other classical estimates, which can still be achieved even if QPE does not project directly onto the ground state. Finally, we show how the energy distribution can help in identifying potential quantum advantage. △ Less

Submitted 8 February, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Resistivity of the two-dimensional Bose-Hubbard model at weak coupling

Authors: Eduardo O. Rizzatti, Erich J. Mueller

Abstract: We calculate the weak-coupling resistivity of the two-dimensional Bose Hubbard model, comparing with the more familiar fermionic case. At high temperature the resistivity is linear in $T$, while in the low temperature normal state it is exponentially suppressed. We explore the density dependence and calculate the momentum relaxation rate. We calculate the weak-coupling resistivity of the two-dimensional Bose Hubbard model, comparing with the more familiar fermionic case. At high temperature the resistivity is linear in $T$, while in the low temperature normal state it is exponentially suppressed. We explore the density dependence and calculate the momentum relaxation rate. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Constraints on a split superconducting transition under uniaxial strain in Sr$_2$RuO$_4$ from scanning SQUID microscopy

Authors: Eli Mueller, Yusuke Iguchi, Christopher Watson, Clifford Hicks, Yoshiteru Maeno, Kathryn Moler

Abstract: More than two decades after the discovery of superconductivity in Sr$_2$RuO$_4$, it is still unclear whether the order parameter has a single component or two degenerate components. For any two-component scenario, application of uniaxial strain is expected to lift the degeneracy, generating two distinct phase transitions. The presence of a second (lower-temperature) transition may be observable by… ▽ More More than two decades after the discovery of superconductivity in Sr$_2$RuO$_4$, it is still unclear whether the order parameter has a single component or two degenerate components. For any two-component scenario, application of uniaxial strain is expected to lift the degeneracy, generating two distinct phase transitions. The presence of a second (lower-temperature) transition may be observable by probes that are sensitive to changes in the London penetration depth, $λ$, as a function of temperature, $T$. Here, we use scanning SQUID microscopy combined with a uniaxial strain device to test for a second transition under strain. We only observe a single transition. Within the temperature range where a second transition has been suggested by $μ$SR measurements, we further place a tight upper bound of less than 1% on the change in the zero temperature superfluid density $n_s\proptoλ^{-2}(0)$ due to a second transition, suggesting that such a transition does not occur. These results constrain theories of the order parameter in Sr$_2$RuO$_4$. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Vortex dynamics induced by scanning SQUID susceptometry

Authors: Logan Bishop-Van Horn, Eli Mueller, Kathryn A. Moler

Abstract: We measured the local magnetic response of a niobium thin film by applying a millitesla-scale AC magnetic field using a micron-scale field coil and detecting the response with a micron-scale pickup loop in a scanning superconducting quantum interference device (SQUID) susceptometry measurement. Near the film's critical temperature, we observed a step-like nonlinear and dissipative magnetic respons… ▽ More We measured the local magnetic response of a niobium thin film by applying a millitesla-scale AC magnetic field using a micron-scale field coil and detecting the response with a micron-scale pickup loop in a scanning superconducting quantum interference device (SQUID) susceptometry measurement. Near the film's critical temperature, we observed a step-like nonlinear and dissipative magnetic response due to the dynamics of a small number of vortex-antivortex pairs induced in the film by the local applied AC field. We modeled the dynamics of the measurement using a combined two-dimensional London-Maxwell and time-dependent Ginzburg-Landau approach, allowing us to construct a detailed real-space picture of the vortex motion causing the observed dissipative response. This work pushes scanning SQUID susceptometry of two-dimensional superconductors beyond the regime of linear response and lays the foundation for microscopic studies of vortex dynamics and pinning in superconducting devices and more exotic materials systems. △ Less

Submitted 15 June, 2023; v1 submitted 25 April, 2023; originally announced April 2023.

Comments: 15 pages, 8 figures. The first two authors contributed equally. This is the accepted version

Journal ref: Phys. Rev. B 107, 224509 (2023)

Subsystem symmetry, spin glass order, and criticality from random measurements in a two-dimensional Bacon-Shor circuit

Authors: Vaibhav Sharma, Chao-Ming Jian, Erich J Mueller

Abstract: We study a 2D measurement-only random circuit motivated by the Bacon-Shor error correcting code. We find a rich phase diagram as one varies the relative probabilities of measuring nearest neighbor Pauli XX and ZZ check operators. In the Bacon-Shor code, these checks commute with a group of stabilizer and logical operators, which therefore represent conserved quantities. Described as a subsystem sy… ▽ More We study a 2D measurement-only random circuit motivated by the Bacon-Shor error correcting code. We find a rich phase diagram as one varies the relative probabilities of measuring nearest neighbor Pauli XX and ZZ check operators. In the Bacon-Shor code, these checks commute with a group of stabilizer and logical operators, which therefore represent conserved quantities. Described as a subsystem symmetry, these conservation laws lead to a continuous phase transition between an X-basis and Z-basis spin glass order. The two phases are separated by a critical point where the entanglement entropy between two halves of an L X L system scales as L ln L, a logarithmic violation of the area law. We generalize to a model where the check operators break the subsystem symmetries (and the Bacon-Shor code structure). In tension with established heuristics, we find that the phase transition is replaced by a smooth crossover, and the X- and Z-basis spin glass orders spatially coexist. Additionally, if we approach the line of subsystem symmetries away from the critical point in the phase diagram, some spin glass order parameters jump discontinuously △ Less

Submitted 3 August, 2023; v1 submitted 3 March, 2023; originally announced March 2023.

Comments: 11 pages, 7 figures ; Minor revisions in response to referee comments and corrected typos

Journal ref: Phys. Rev. B 108, 024205 (2023)

Quantum simulation of battery materials using ionic pseudopotentials

Authors: Modjtaba Shokrian Zini, Alain Delgado, Roberto dos Reis, Pablo A. M. Casares, Jonathan E. Mueller, Arne-Christian Voigt, Juan Miguel Arrazola

Abstract: Ionic pseudopotentials are widely used in classical simulations of materials to model the effective potential due to the nucleus and the core electrons. Modeling fewer electrons explicitly results in a reduction in the number of plane waves needed to accurately represent the states of a system. In this work, we introduce a quantum algorithm that uses pseudopotentials to reduce the cost of simulati… ▽ More Ionic pseudopotentials are widely used in classical simulations of materials to model the effective potential due to the nucleus and the core electrons. Modeling fewer electrons explicitly results in a reduction in the number of plane waves needed to accurately represent the states of a system. In this work, we introduce a quantum algorithm that uses pseudopotentials to reduce the cost of simulating periodic materials on a quantum computer. We use a qubitization-based quantum phase estimation algorithm that employs a first-quantization representation of the Hamiltonian in a plane-wave basis. We address the challenge of incorporating the complexity of pseudopotentials into quantum simulations by developing highly-optimized compilation strategies for the qubitization of the Hamiltonian. This includes a linear combination of unitaries decomposition that leverages the form of separable pseudopotentials. Our strategies make use of quantum read-only memory subroutines as a more efficient alternative to quantum arithmetic. We estimate the computational cost of applying our algorithm to simulating lithium-excess cathode materials for batteries, where more accurate simulations are needed to inform strategies for gaining reversible access to the excess capacity they offer. We estimate the number of qubits and Toffoli gates required to perform sufficiently accurate simulations with our algorithm for three materials: lithium manganese oxide, lithium nickel-manganese oxide, and lithium manganese oxyfluoride. Our optimized compilation strategies result in a pseudopotential-based quantum algorithm with a total Toffoli cost four orders of magnitude lower than the previous state of the art for a fixed target accuracy. △ Less

Submitted 4 July, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Journal ref: Quantum 7, 1049 (2023)

Observation of suppressed viscosity in the normal state of $^3$He due to superfluid fluctuations

Authors: Rakin N. Baten, Yefan Tian, Eric N. Smith, Erich Mueller, Jeevak M. Parpia

Abstract: By monitoring the quality factor of a quartz tuning fork oscillator we have observed a fluctuation-driven reduction in the viscosity of bulk $^3$He in the normal state near the superfluid transition temperature, $T_c$. These fluctuations, which are only found within $100 μ$K of $T_c$, play a vital role in the theoretical modeling of ordering; they encode details about the Fermi liquid parameters,… ▽ More By monitoring the quality factor of a quartz tuning fork oscillator we have observed a fluctuation-driven reduction in the viscosity of bulk $^3$He in the normal state near the superfluid transition temperature, $T_c$. These fluctuations, which are only found within $100 μ$K of $T_c$, play a vital role in the theoretical modeling of ordering; they encode details about the Fermi liquid parameters, pairing symmetry, and scattering phase shifts. They will be of crucial importance for transport probes of the topologically nontrivial features of superfluid $^3$He under strong confinement. Here we characterize the temperature and pressure dependence of the fluctuation signature, finding data collapse consistent with the predicted theoretical behavior. △ Less

Submitted 25 August, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

Comments: 5 figures

Supercooling of the A phase of $^3$He

Authors: Yefan Tian, Dmytro Lotnyk, Anna Eyal, Kuang Zhang, Nikolay Zhelev, T. S. Abhilash, Aldo Chavez, Eric Smith, Mark Hindmarsh, John Saunders, Erich Mueller, Jeevak Parpia

Abstract: Because of the extreme purity, lack of disorder, and complex order parameter, the first-order superfluid $^3$He A-B transition is the leading model system for first order transitions in the early universe. Here we report on the path dependence of the supercooling of the A phase over a wide range of pressures below 29.3 bar at nearly zero magnetic field. The A phase can be cooled significantly belo… ▽ More Because of the extreme purity, lack of disorder, and complex order parameter, the first-order superfluid $^3$He A-B transition is the leading model system for first order transitions in the early universe. Here we report on the path dependence of the supercooling of the A phase over a wide range of pressures below 29.3 bar at nearly zero magnetic field. The A phase can be cooled significantly below the thermodynamic A-B transition temperature. While the extent of supercooling is highly reproducible, it depends strongly upon the cooling trajectory: The metastability of the A phase is enhanced by transiting through regions where the A phase is more stable. We provide evidence that some of the additional supercooling is due to the elimination of B phase seeds formed upon passage through the superfluid transition. A greater understanding of the physics is essential before the $^3$He can be exploited to model transitions in the early universe. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: 6 Figures

Effect of intense x-ray free-electron laser transient gratings on the magnetic domain structure of Tm:YIG

Authors: Victor Ukleev, Max Burian, Sebastian Gliga, C. A. F. Vaz, Benedikt Rösner, Danny Fainozzi, Gediminas Seniutinas, Adam Kubec, Roman Mankowsky, Henrik T. Lemke, Ethan R. Rosenberg, Caroline A. Ross, Elisabeth Müller, Christian David, Cristian Svetina, Urs Staub

Abstract: Magnetic patterns can be controlled globally using fields or spin polarized currents. In contrast, the local control of the magnetization on the nanometer length scale remains challenging. Here, we demonstrate how magnetic domain patterns in a Tm-doped yttrium iron garnet (Tm:YIG) thin film with perpendicular magnetic anisotropy can be permanently and locally imprinted by high intensity photon pul… ▽ More Magnetic patterns can be controlled globally using fields or spin polarized currents. In contrast, the local control of the magnetization on the nanometer length scale remains challenging. Here, we demonstrate how magnetic domain patterns in a Tm-doped yttrium iron garnet (Tm:YIG) thin film with perpendicular magnetic anisotropy can be permanently and locally imprinted by high intensity photon pulses of a hard x-ray transient grating (XTG). Micromagnetic simulations provide a qualitative understanding of the observed changes in the orientation of magnetic domains in Tm:YIG and XTG-induced changes. The presented results offer a route for the local manipulation of the magnetic state using hard XTG. △ Less

Submitted 3 March, 2023; v1 submitted 8 August, 2022; originally announced August 2022.

Comments: 15 pages, 3 figures

Journal ref: J. Appl. Phys. 133, 123902 (2023)

Role of Dy on the magnetic properties of orthorhombic DyFeO3

Authors: Banani Biswas, Veronica F. Michel, Oystein S. Fjellvag, Gesara Bimashofer, Max Dobeli, Michal Jambor, Lukas Keller, Elisabeth Muller, Victor Ukleev, Ekaterina V. Pomjakushina, Deepak Singh, Uwe Stuhr, Carlos A. F. Vaz, Thomas Lippert, Christof W. Schneider

Abstract: Orthoferrites are a class of magnetic materials with a magnetic ordering temperature above 600 K, predominant G-type antiferromagnetic ordering of the Fe-spin system and, depending on the rare-earth ion, a spin reorientation of the Fe spin taking place at lower temperatures. DyFeO3 is of particular interest since the spin reorientation is classified as a Morin transition with the transition temper… ▽ More Orthoferrites are a class of magnetic materials with a magnetic ordering temperature above 600 K, predominant G-type antiferromagnetic ordering of the Fe-spin system and, depending on the rare-earth ion, a spin reorientation of the Fe spin taking place at lower temperatures. DyFeO3 is of particular interest since the spin reorientation is classified as a Morin transition with the transition temperature depending strongly on the Dy-Fe interaction. Here, we report a detailed study of the magnetic and structural properties of microcrystalline DyFeO3 powder and bulk single crystal using neutron diffraction and magnetometry between 1.5 and 450 K. We find that, while the magnetic properties of the single crystal are largely as expected, the powder shows strongly modified magnetic properties, including a modified spin reorientation and a smaller Dy-Fe interaction energy of the order of 10 μeV. Subtle structural differences between powder and single crystal show that they belong to distinct magnetic space groups. In addition, the Dy ordering at 2 K in the powder is incommensurate, with a modulation vector of 0.0173(5) c*, corresponding to a periodicity of ~58 unit cells. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: contains supplementary material

Journal ref: Published in Phys. Rev. Materials 6, 074401 (2022)

Role of conservation laws in the Density Matrix Renormalization Group

Authors: Thomas G. Kiely, Erich J. Mueller

Abstract: We explore matrix product state approximations to wavefunctions which have spontaneously broken symmetries or are critical. We are motivated by the fact that symmetries, and their associated conservation laws, lead to block-sparse matrix product states. Numerical calculations which take advantage of these symmetries run faster and require less memory. However, in symmetry-broken and critical phase… ▽ More We explore matrix product state approximations to wavefunctions which have spontaneously broken symmetries or are critical. We are motivated by the fact that symmetries, and their associated conservation laws, lead to block-sparse matrix product states. Numerical calculations which take advantage of these symmetries run faster and require less memory. However, in symmetry-broken and critical phases the block sparse ansatz yields less accurate energies. We characterize the role of conservation laws in matrix product states and determine when it is beneficial to make use of them. △ Less

Submitted 15 December, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 10 pages, 6 figures

Journal ref: Phys. Rev. B 106, 235126 (2022)

Resonant enhancement of particle emission from a parametrically driven condensate in a one-dimensional lattice

Authors: L. Q. Lai, Y. B. Yu, Erich J. Mueller

Abstract: Motivated by recent experiments, we investigate particle emission from a Bose-Einstein condensate in a one-dimensional lattice, where the interaction strength is periodically modulated. The modulated interactions parametrically excite a collective mode, leading to density oscillations. These collective oscillations in turn drive particle emission. This multistep process amplifies the drive, produc… ▽ More Motivated by recent experiments, we investigate particle emission from a Bose-Einstein condensate in a one-dimensional lattice, where the interaction strength is periodically modulated. The modulated interactions parametrically excite a collective mode, leading to density oscillations. These collective oscillations in turn drive particle emission. This multistep process amplifies the drive, producing larger particle jets. We find that the amplitude dependence of the emission rate has a characteristic threshold behavior, as seen in experiments. △ Less

Submitted 2 September, 2022; v1 submitted 4 May, 2022; originally announced May 2022.

Comments: 5 pages, 6 figures

Journal ref: Phys. Rev. A 106, 033302 (2022)

Simulating key properties of lithium-ion batteries with a fault-tolerant quantum computer

Authors: Alain Delgado, Pablo A. M. Casares, Roberto dos Reis, Modjtaba Shokrian Zini, Roberto Campos, Norge Cruz-Hernández, Arne-Christian Voigt, Angus Lowe, Soran Jahangiri, M. A. Martin-Delgado, Jonathan E. Mueller, Juan Miguel Arrazola

Abstract: There is a pressing need to develop new rechargeable battery technologies that can offer higher energy storage, faster charging, and lower costs. Despite the success of existing methods for the simulation of battery materials, they can sometimes fall short of delivering accurate and reliable results. Quantum computing has been discussed as an avenue to overcome these issues, but only limited work… ▽ More There is a pressing need to develop new rechargeable battery technologies that can offer higher energy storage, faster charging, and lower costs. Despite the success of existing methods for the simulation of battery materials, they can sometimes fall short of delivering accurate and reliable results. Quantum computing has been discussed as an avenue to overcome these issues, but only limited work has been done to outline how they may impact battery simulations. In this work, we provide a detailed answer to the following question: how can a quantum computer be used to simulate key properties of a lithium-ion battery? Based on recently-introduced first-quantization techniques, we lay out an end-to-end quantum algorithm for calculating equilibrium cell voltages, ionic mobility, and thermal stability. These can be obtained from ground-state energies of materials, which is the core calculation executed by the quantum computer using qubitization-based quantum phase estimation. The algorithm includes explicit methods for preparing approximate ground states of periodic materials in first quantization. We bring these insights together to perform the first estimation of the resources required to implement a quantum algorithm for simulating a realistic cathode material, dilithium iron silicate. △ Less

Submitted 6 February, 2023; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 31 pages, 14 figures

Journal ref: Phys. Rev. A 106, 032428 (2022)

Engineered Dissipation for Quantum Information Science

Authors: Patrick M. Harrington, Erich Mueller, Kater Murch

Abstract: Quantum information processing relies on precise control of non-classical states in the presence of many uncontrolled environmental degrees of freedom -- requiring careful orchestration of how the relevant degrees of freedom interact with that environment. These interactions are often viewed as detrimental, as they dissipate energy and decohere quantum states. Nonetheless, when controlled, dissipa… ▽ More Quantum information processing relies on precise control of non-classical states in the presence of many uncontrolled environmental degrees of freedom -- requiring careful orchestration of how the relevant degrees of freedom interact with that environment. These interactions are often viewed as detrimental, as they dissipate energy and decohere quantum states. Nonetheless, when controlled, dissipation is an essential tool for manipulating quantum information: Dissipation engineering enables quantum measurement, quantum state preparation, and quantum state stabilization. The progress of quantum device technology, marked by improvements of characteristic coherence times and extensible architectures for quantum control, has coincided with the development of such dissipation engineering tools which interface quantum and classical degrees of freedom. This Review presents dissipation as a fundamental aspect of the measurement and control of quantum devices and highlights the role of dissipation engineering for quantum error correction and quantum simulation that enables quantum information processing on a practical scale. △ Less

Submitted 30 May, 2022; v1 submitted 10 February, 2022; originally announced February 2022.

Comments: Review article: 32 pages, 4 figures, 2 boxes

Journal ref: Nat Rev Phys 4, 660-671 (2022)

Water flow in single-wall nanotubes: Oxygen makes it slip, hydrogen makes it stick

Authors: Fabian L Thiemann, Christoph Schran, Patrick Rowe, Erich A Müller, Angelos Michaelides

Abstract: Experimental measurements have reported ultra-fast and radius-dependent water transport in carbon nanotubes which are absent in boron nitride nanotubes. Despite considerable effort, the origin of this contrasting (and fascinating) behaviour is not understood. Here, with the aid of machine learning-based molecular dynamics simulations that deliver first-principles accuracy, we investigate water tra… ▽ More Experimental measurements have reported ultra-fast and radius-dependent water transport in carbon nanotubes which are absent in boron nitride nanotubes. Despite considerable effort, the origin of this contrasting (and fascinating) behaviour is not understood. Here, with the aid of machine learning-based molecular dynamics simulations that deliver first-principles accuracy, we investigate water transport in single-wall carbon and boron nitride nanotubes. Our simulations reveal a large, radius-dependent hydrodynamic slippage on both materials with water experiencing indeed a $\approx 5$ times lower friction on carbon surfaces compared to boron nitride. Analysis of the diffusion mechanisms across the two materials reveals that the fast water transport on carbon is governed by facile oxygen motion, whereas the higher friction on boron nitride arises from specific hydrogen-nitrogen interactions. This work not only delivers a clear reference of unprecedented accuracy for water flow in single-wall nanotubes, but also provides detailed mechanistic insight into its radius and material dependence for future technological application. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Superfluidity in the 1D Bose-Hubbard Model

Authors: Thomas G. Kiely, Erich J. Mueller

Abstract: We study superfluidity in the 1D Bose-Hubbard model using a variational matrix product state technique. We determine the superfluid density as a function of the Hubbard parameters by calculating the energy cost of phase twists in the thermodynamic limit. As the system is critical, correlation functions decay as power laws and the entanglement entropy grows with the bond dimension of our variationa… ▽ More We study superfluidity in the 1D Bose-Hubbard model using a variational matrix product state technique. We determine the superfluid density as a function of the Hubbard parameters by calculating the energy cost of phase twists in the thermodynamic limit. As the system is critical, correlation functions decay as power laws and the entanglement entropy grows with the bond dimension of our variational state. We relate the resulting scaling laws to the superfluid density. We compare two different algorithms for optimizing the infinite matrix product state and develop a physical explanation why one of them (VUMPS) is more efficient than the other (iDMRG). Finally, we comment on finite-temperature superfluidity in one dimension and how our results can be realized in cold atom experiments. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: 16 pages, 9 figures (some equations use graphicx)

The BrainScaleS-2 accelerated neuromorphic system with hybrid plasticity

Authors: Christian Pehle, Sebastian Billaudelle, Benjamin Cramer, Jakob Kaiser, Korbinian Schreiber, Yannik Stradmann, Johannes Weis, Aron Leibfried, Eric Müller, Johannes Schemmel

Abstract: Since the beginning of information processing by electronic components, the nervous system has served as a metaphor for the organization of computational primitives. Brain-inspired computing today encompasses a class of approaches ranging from using novel nano-devices for computation to research into large-scale neuromorphic architectures, such as TrueNorth, SpiNNaker, BrainScaleS, Tianjic, and Lo… ▽ More Since the beginning of information processing by electronic components, the nervous system has served as a metaphor for the organization of computational primitives. Brain-inspired computing today encompasses a class of approaches ranging from using novel nano-devices for computation to research into large-scale neuromorphic architectures, such as TrueNorth, SpiNNaker, BrainScaleS, Tianjic, and Loihi. While implementation details differ, spiking neural networks - sometimes referred to as the third generation of neural networks - are the common abstraction used to model computation with such systems. Here we describe the second generation of the BrainScaleS neuromorphic architecture, emphasizing applications enabled by this architecture. It combines a custom analog accelerator core supporting the accelerated physical emulation of bio-inspired spiking neural network primitives with a tightly coupled digital processor and a digital event-routing network. △ Less

Submitted 3 February, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 22 pages, 10 figures; amended funding acknowledgements, added one citation

Rotating Bose gas dynamically entering the lowest Landau level

Authors: Vaibhav Sharma, Erich J Mueller

Abstract: Motivated by recent experiments, we model the dynamics of a condensed Bose gas in a rotating anisotropic trap, where the equations of motion are analogous to those of charged particles in a magnetic field. As the rotation rate is ramped from zero to the trapping frequency, the condensate stretches along one direction and is squeezed along another, becoming long and thin. When the trap anisotropy i… ▽ More Motivated by recent experiments, we model the dynamics of a condensed Bose gas in a rotating anisotropic trap, where the equations of motion are analogous to those of charged particles in a magnetic field. As the rotation rate is ramped from zero to the trapping frequency, the condensate stretches along one direction and is squeezed along another, becoming long and thin. When the trap anisotropy is slowly switched off on a particular timescale, the condensate is left in the lowest Landau level. We use a time dependent variational approach to quantify these dynamics and give intuitive arguments about the structure of the condensate wavefunction. This preparation of a lowest Landau level condensate can be an important first step in realizing bosonic analogs of quantum Hall states. △ Less

Submitted 7 February, 2022; v1 submitted 19 November, 2021; originally announced November 2021.

Comments: 8 pages, 5 figures; Minor revisions in response to referee comments and corrected typos

Journal ref: Phys. Rev. A 105, 023310 (2022)

Semi-empirical Quantum Optics for Mid-Infrared Molecular Nanophotonics

Authors: Johan F. Triana, Mauricio Arias, Jun Nishida, Eric Muller, Roland Wilcken, Samuel C. Johnson, Aldo Delgado, Markus B. Raschke, Felipe Herrera

Abstract: Nanoscale infrared (IR) resonators with sub-diffraction limited mode volumes and open geometries have emerged as new platforms for implementing cavity quantum electrodynamics (QED) at room temperature. The use of infrared (IR) nano-antennas and tip nanoprobes to study strong light-matter coupling of molecular vibrations with the vacuum field can be exploited for IR quantum control with nanometer a… ▽ More Nanoscale infrared (IR) resonators with sub-diffraction limited mode volumes and open geometries have emerged as new platforms for implementing cavity quantum electrodynamics (QED) at room temperature. The use of infrared (IR) nano-antennas and tip nanoprobes to study strong light-matter coupling of molecular vibrations with the vacuum field can be exploited for IR quantum control with nanometer and femtosecond resolution. In order to accelerate the development of molecule-based quantum nano-photonic devices in the mid-IR, we develop a generally applicable semi-empirical quantum optics approach to describe light-matter interaction in systems driven by mid-IR femtosecond laser pulses. The theory is shown to reproduce recent experiments on the acceleration of the vibrational relaxation rate in infrared nanostructures, and also provide physical insights for the implementation of coherent phase rotations of the near-field using broadband nanotips. We then apply the quantum framework to develop general tip-design rules for the experimental manipulation of vibrational strong coupling and Fano interference effects in open infrared resonators. We finally propose the possibility of transferring the natural anharmonicity of molecular vibrational levels to the resonator near-field in the weak coupling regime, to implement intensity-dependent phase shifts of the coupled system response with strong pulses. Our semi-empirical quantum theory is equivalent to first-principles techniques based on Maxwell's equations, but its lower computational cost suggests its use a rapid design tool for the development of strongly-coupled infrared nanophotonic hardware for applications ranging from quantum control of materials to quantum information processing. △ Less

Submitted 14 October, 2021; originally announced October 2021.

Comments: 18 pages, 7 figures

Journal ref: J. Chem. Phys. 156, 124110, 2022

Dynamics of Spin Helices in the One-Dimensional $XX$ Model

Authors: Darren Pereira, Erich J. Mueller

Abstract: Motivated by cold-atom experiments and a desire to understand far-from-equilibrium quantum transport, we analytically study the dynamics of spin helices in the one-dimensional $XX$ model. We use a Jordan-Wigner transformation to map the spin chain onto a non-interacting Fermi gas with simple equations of motion. The resulting dynamics are nontrivial, however, as the spin-helix initial condition co… ▽ More Motivated by cold-atom experiments and a desire to understand far-from-equilibrium quantum transport, we analytically study the dynamics of spin helices in the one-dimensional $XX$ model. We use a Jordan-Wigner transformation to map the spin chain onto a non-interacting Fermi gas with simple equations of motion. The resulting dynamics are nontrivial, however, as the spin-helix initial condition corresponds to a highly nonequilibrium distribution of the fermions. We find a separation of timescales between the in-plane and out-of-plane spin dynamics. We gain insights from analyzing the case of a uniform spin chain and from a semiclassical model. One of our key findings is that the spin correlation functions decay as $t^{-1/2}$ at long time, in contrast to the experimentally observed exponential decay. △ Less

Submitted 20 October, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 10 pages, 6 figures

Journal ref: Phys. Rev. A 106, 043306 (2022)

Anomalous Resistivity at Weak Coupling

Authors: Thomas G. Kiely, Erich J. Mueller

Abstract: Recent cold atom experiments have observed bad and strange metal behaviors in strongly-interacting Fermi-Hubbard systems. Motivated by these results, we calculate the thermoelectric transport properties of a 2D Fermi-Hubbard system in the weak coupling limit using quantum kinetic theory. We find that many features attributed to strong correlations are also found at weak coupling. In particular, fo… ▽ More Recent cold atom experiments have observed bad and strange metal behaviors in strongly-interacting Fermi-Hubbard systems. Motivated by these results, we calculate the thermoelectric transport properties of a 2D Fermi-Hubbard system in the weak coupling limit using quantum kinetic theory. We find that many features attributed to strong correlations are also found at weak coupling. In particular, for temperatures $T\gtrsim t$ the electrical resistivity is nearly linear in temperature despite the fact that the quasiparticle scattering rate is non-linear and changes by nearly an order of magnitude. We argue that this asymptotic behavior is a general feature of systems with a finite spectral width, which implies that there is no MIR bound on the resistivity in single-band models. Due to nesting, the $T$-linear resistivity persists down to $T=0$ at half filling. Our work sheds light on the transport regime in ultracold atom experiments, which can differ substantially from that of condensed matter systems. Disentangling these band-structure effects from the physics of strong correlations is a major challenge for future experiments. △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: 5 pages, 2 figures

Density-Matrix Renormalization Group for Continuous Quantum Systems

Authors: Shovan Dutta, Anton Buyskikh, Andrew J. Daley, Erich J. Mueller

Abstract: We introduce a versatile and practical framework for applying matrix product state techniques to continuous quantum systems. We divide space into multiple segments and generate continuous basis functions for the many-body state in each segment. By combining this mapping with existing numerical Density-Matrix Renormalization Group routines, we show how one can accurately obtain the ground-state wav… ▽ More We introduce a versatile and practical framework for applying matrix product state techniques to continuous quantum systems. We divide space into multiple segments and generate continuous basis functions for the many-body state in each segment. By combining this mapping with existing numerical Density-Matrix Renormalization Group routines, we show how one can accurately obtain the ground-state wave function, spatial correlations, and spatial entanglement entropy directly in the continuum. For a prototypical mesoscopic system of strongly-interacting bosons we demonstrate faster convergence than standard grid-based discretization. We illustrate the power of our approach by studying a superfluid-insulator transition in an external potential. We outline how one can directly apply or generalize this technique to a wide variety of experimentally relevant problems across condensed matter physics and quantum field theory. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Main paper: 7 pages, 5 figures; Supplement: 13 pages, 5 figures

Journal ref: Phys. Rev. Lett. 128, 230401 (2022)

Driven dissipative preparation of few-body Laughlin states of Rydberg polaritons in twisted cavities

Authors: Kristina R. Colladay, Erich J. Mueller

Abstract: We present a driven dissipative protocol for creating an optical analog of the Laughlin state in a system of Rydberg polaritons in a twisted optical cavity. We envision resonantly driving the system into a 4-polariton state by injecting photons in carefully selected modes. The dissipative nature of the polariton-polariton interactions leads to a decay into a two-polariton analog of the Laughlin st… ▽ More We present a driven dissipative protocol for creating an optical analog of the Laughlin state in a system of Rydberg polaritons in a twisted optical cavity. We envision resonantly driving the system into a 4-polariton state by injecting photons in carefully selected modes. The dissipative nature of the polariton-polariton interactions leads to a decay into a two-polariton analog of the Laughlin state. Generalizations of this technique could be used to explore fractional statistics and anyon based quantum information processing. We also model recent experiments that attempt to coherently drive into this same state. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Comments: 12 pages, 10 figures

Defect-Dependent Corrugation in Graphene

Authors: Fabian L. Thiemann, Patrick Rowe, Andrea Zen, Erich A. Müller, Angelos Michaelides

Abstract: Graphene's intrinsically corrugated and wrinkled topology fundamentally influences its electronic, mechanical, and chemical properties. Experimental techniques allow the manipulation of pristine graphene and the controlled production of defects which allows to control the atomic out-of-plane fluctuations and, thus, tune graphene's properties. Here, we perform large scale machine learning-driven mo… ▽ More Graphene's intrinsically corrugated and wrinkled topology fundamentally influences its electronic, mechanical, and chemical properties. Experimental techniques allow the manipulation of pristine graphene and the controlled production of defects which allows to control the atomic out-of-plane fluctuations and, thus, tune graphene's properties. Here, we perform large scale machine learning-driven molecular dynamics simulations to understand the impact of defects on the structure of graphene. We find that defects cause significantly higher corrugation leading to a strongly wrinkled surface. The magnitude of this structural transformation strongly depends on the defect concentration and specific type of defect. Analysing the atomic neighborhood of the defects reveals that the extent of these morphological changes depends on the preferred geometrical orientation and the interactions between defects. While our work highlights that defects can strongly affect graphene's morphology, it also emphasises the differences between distinct types by linking the global structure to the local environment of the defects. △ Less

Submitted 14 September, 2021; v1 submitted 2 July, 2021; originally announced July 2021.

Emission of particles from a parametrically driven condensate in a one-dimensional lattice

Authors: L. Q. Lai, Y. B. Yu, Erich J. Mueller

Abstract: Motivated by recent experiments, we calculate particle emission from a Bose-Einstein condensate trapped in a single deep well of a one-dimensional lattice when the interaction strength is modulated. In addition to pair emission, which has been widely studied, we observe single-particle emission. Within linear response, we are able to write closed-form expressions for the single-particle emission r… ▽ More Motivated by recent experiments, we calculate particle emission from a Bose-Einstein condensate trapped in a single deep well of a one-dimensional lattice when the interaction strength is modulated. In addition to pair emission, which has been widely studied, we observe single-particle emission. Within linear response, we are able to write closed-form expressions for the single-particle emission rates and reduce the pair emission rates to one-dimensional integrals. The full nonlinear theory of single-particle emission is reduced to a single variable integrodifferential equation, which we numerically solve. △ Less

Submitted 7 September, 2021; v1 submitted 17 June, 2021; originally announced June 2021.

Comments: 8 pages, 6 figures

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

Transport in the 2D Fermi-Hubbard Model: Lessons from Weak Coupling

Authors: Thomas G. Kiely, Erich J. Mueller

Abstract: We use quantum kinetic theory to calculate the thermoelectric transport properties of the 2D single band Fermi-Hubbard model in the weak coupling limit. For generic filling, we find that the high-temperature limiting behaviors of the electrical ($\sim T$) and thermal ($\sim T^2$) resistivities persist down to temperatures of order the hopping matrix element $T\sim t$, almost an order of magnitude… ▽ More We use quantum kinetic theory to calculate the thermoelectric transport properties of the 2D single band Fermi-Hubbard model in the weak coupling limit. For generic filling, we find that the high-temperature limiting behaviors of the electrical ($\sim T$) and thermal ($\sim T^2$) resistivities persist down to temperatures of order the hopping matrix element $T\sim t$, almost an order of magnitude below the bandwidth. At half filling, perfect nesting leads to anomalous low temperature scattering and nearly $T$-linear electrical resistivity at all temperatures. We hypothesize that the $T$-linear resistivity observed in recent cold atom experiments is continuously connected to this weak coupling physics and suggest avenues for experimental verification. We find a number of other novel thermoelectric results, such as a low-temperature Wiedemann-Franz law with Lorenz coefficient $5π^2/36$. △ Less

Submitted 25 August, 2021; v1 submitted 8 June, 2021; originally announced June 2021.

Comments: 14 pages, 6 figures

Bath mediated decay of density waves in a disordered Bose lattice gas

Authors: Sridhar Prabhu, Erich J. Mueller

Abstract: Motivated by a recent experiment, we study the dynamics of bosons in a disordered optical lattice, interacting with a variably sized bath of disorder free atoms. As the number of particles in the bath is increased, there is a transition between "localized" and "ergodic" behavior, which are characterized by the long-time behavior of an initial density wave. We model the dynamics with a stochastic m… ▽ More Motivated by a recent experiment, we study the dynamics of bosons in a disordered optical lattice, interacting with a variably sized bath of disorder free atoms. As the number of particles in the bath is increased, there is a transition between "localized" and "ergodic" behavior, which are characterized by the long-time behavior of an initial density wave. We model the dynamics with a stochastic mean field theory, reproducing the central observations of the experiment. A key conclusion from our study is that particle loss plays an important role. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: 7 pages, 4 figures

ReaxFF Reactive Force Field Development for Cu/Si Systems and application to Copper Cluster Formation During Cu Diffusion Inside Silicon

Authors: Kamyar Akbari Roshan, Mahdi Khajeh Talkhoncheh, Jonathan E. Mueller, William A. Goddard III, Adri C. T. van Duin

Abstract: Transition metal impurities such as nickel, copper, and iron, in solid-state materials like silicon have a significant impact on the electrical performance of integrated circuits and solar cells. To study the impact of copper impurities inside bulk silicon on the electrical properties of the material, one needs to understand the configurational space of copper atoms incorporated inside the silicon… ▽ More Transition metal impurities such as nickel, copper, and iron, in solid-state materials like silicon have a significant impact on the electrical performance of integrated circuits and solar cells. To study the impact of copper impurities inside bulk silicon on the electrical properties of the material, one needs to understand the configurational space of copper atoms incorporated inside the silicon lattice. In this work, we performed ReaxFF reactive force field based molecular dynamics simulations, studying different configurations of individual and crystalline copper atoms inside bulk silicon by looking at the diffusional behavior of copper in silicon. The ReaxFF Cu/Si parameter set was developed by training against DFT data, including the energy barrier for an individual Cu-atom inside a silicon lattice. We found that the diffusion of copper atoms has a direct relationship with the temperature. Moreover, it is also shown that individual copper atoms start to clusterize inside bulk silicon at elevated temperatures. Our simulation results provide a comprehensive picture of the effects of temperature and copper concentration on the crystallization of individual copper inside silicon lattice. Finally, the stress-strain relationship of Cu/Si compounds under uniaxial tensile loading have been obtained. Our results indicate a decrease in the elastic modulus with increasing level of Cu-impurity concentration. We observe spontaneous microcracking of the Si during the stress-strain tests as a consequence of the formation of a small Cu clusters adjacent to the Si surface. △ Less

Submitted 1 February, 2021; originally announced February 2021.

Comments: 29 pages, 13 figures

MSC Class: H.1

Driven-dissipative control of cold atoms in tilted optical lattices

Authors: Vaibhav Sharma, Erich J Mueller

Abstract: We present a sequence of driven-dissipative protocols for controlling cold atoms in tilted optical lattices. These experimentally accessible examples are templates that demonstrate how dissipation can be used to manipulate quantum many-body systems. We consider bosonic atoms trapped in a tilted optical lattice, immersed in a superfluid bath, and excited by coherent Raman lasers. With these ingredi… ▽ More We present a sequence of driven-dissipative protocols for controlling cold atoms in tilted optical lattices. These experimentally accessible examples are templates that demonstrate how dissipation can be used to manipulate quantum many-body systems. We consider bosonic atoms trapped in a tilted optical lattice, immersed in a superfluid bath, and excited by coherent Raman lasers. With these ingredients, we are able to controllably transport atoms in the lattice and produce self-healing quantum states: a Mott insulator and the topologically ordered spin-1 AKLT state. △ Less

Submitted 19 April, 2021; v1 submitted 2 January, 2021; originally announced January 2021.

Comments: Published version after response to referee comments

Journal ref: Phys. Rev. A 103, 043322 (2021)

Path-Dependent Supercooling of the $^3$He Superfluid A-B transition

Authors: Dmytro Lotnyk, Anna Eyal, Nikolay Zhelev, Abhilash Sebastian, Yefan Tian, Aldo Chavez, Eric Smith, John Saunders, Erich Mueller, Jeevak Parpia

Abstract: We examine the discontinuous first-order superfluid $^3$He A to B transition in the vicinity of the polycritical point (2.232 mK and 21.22 bar). We find path-dependent transitions: cooling at fixed pressure yields a well defined transition line in the temperature-pressure plane, but this line can be reliably crossed by depressurizing at nearly constant temperature after transiting $T_{\rm c}$ at a… ▽ More We examine the discontinuous first-order superfluid $^3$He A to B transition in the vicinity of the polycritical point (2.232 mK and 21.22 bar). We find path-dependent transitions: cooling at fixed pressure yields a well defined transition line in the temperature-pressure plane, but this line can be reliably crossed by depressurizing at nearly constant temperature after transiting $T_{\rm c}$ at a higher pressure. This path dependence is not consistent with any of the standard B-phase nucleation mechanisms in the literature. This symmetry breaking transition is a potential simulator for first order transitions in the early universe. △ Less

Submitted 2 April, 2021; v1 submitted 27 December, 2020; originally announced December 2020.

Journal ref: Phys. Rev. Lett. 126, 215301 (2021)

Correlated Insulators in Twisted Bilayer Graphene

Authors: Ipsita Mandal, Jia Yao, Erich J. Mueller

Abstract: Experiments on graphene bilayers, where the top layer is rotated with respect to the one below, have displayed insulating behavior when the moiré bands are partially filled. We calculate the charge distributions in these phases, and estimate the excitation gaps. Experiments on graphene bilayers, where the top layer is rotated with respect to the one below, have displayed insulating behavior when the moiré bands are partially filled. We calculate the charge distributions in these phases, and estimate the excitation gaps. △ Less

Submitted 18 January, 2021; v1 submitted 31 August, 2020; originally announced August 2020.

Comments: version accepted in PRB

Journal ref: Phys. Rev. B 103, 125127 (2021)

Drag in Bose-Fermi Mixtures

Authors: Kai Yen Jee, Erich J. Mueller

Abstract: We use kinetic theory to model the dynamics of a small Bose condensed cloud of heavy particles moving through a larger degenerate Fermi gas of light particles. Varying the Bose-Fermi interaction, we find a crossover between bulk and surface dominated regimes -- where scattering occurs throughout the Bose cloud, or solely on the surface. We calculate the damping and frequency shift of the dipole mo… ▽ More We use kinetic theory to model the dynamics of a small Bose condensed cloud of heavy particles moving through a larger degenerate Fermi gas of light particles. Varying the Bose-Fermi interaction, we find a crossover between bulk and surface dominated regimes -- where scattering occurs throughout the Bose cloud, or solely on the surface. We calculate the damping and frequency shift of the dipole mode in a harmonic trap as a function of the magnetic field controlling an inter-species Feshbach resonance. We find excellent agreement between our stochastic model and the experimental studies of Cs-Li mixtures. △ Less

Submitted 14 January, 2021; v1 submitted 10 August, 2020; originally announced August 2020.

Comments: 10 pages, 7 figures

Journal ref: Phys. Rev. A 103, 033307 (2021)

Generating a Machine-learned Equation of State for Fluid Properties

Authors: Kezheng Zhu, Erich A. Müller

Abstract: Equations of State (EoS) for fluids have been a staple of engineering design and practice for over a century. Available EoS are based on the fitting of a closed-form analytical expression to suitable experimental data. The underlying mathematical structure and the underlying physical model significantly restrain the applicability and accuracy of the resulting EoS. This contribution explores the is… ▽ More Equations of State (EoS) for fluids have been a staple of engineering design and practice for over a century. Available EoS are based on the fitting of a closed-form analytical expression to suitable experimental data. The underlying mathematical structure and the underlying physical model significantly restrain the applicability and accuracy of the resulting EoS. This contribution explores the issues surrounding the substitution of analytical EoS for machine-learned models, in particular, we describe, as a proof of concept, the effectiveness of a machine-learned model to replicate statistical associating fluid theory (SAFT-VR-Mie) EoS for pure fluids. By utilizing Artificial Neural Network and Gaussian Process Regression, predictions of thermodynamic properties such as critical pressure and temperature, vapor pressures and densities of pure model fluids are performed based on molecular descriptors. To quantify the effectiveness of the Machine Learning techniques, a large data set is constructed using the comparisons between the Machine-Learned EoS and the surrogate data set suggest that the proposed approach shows promise as a viable technique for the correlation, extrapolation and prediction of thermophysical properties of fluids. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: Submitted to J. Phys. Chem. B

Machine Learning Potential for Hexagonal Boron Nitride Applied to Thermally and Mechanically Induced Rippling

Authors: Fabian L. Thiemann, Patrick Rowe, Erich A. Müller, Angelos Michaelides

Abstract: We introduce an interatomic potential for hexagonal boron nitride (hBN) based on the Gaussian approximation potential (GAP) machine learning methodology. The potential is based on a training set of configurations collected from density functional theory (DFT) simulations and is capable of treating bulk and multilayer hBN as well as nanotubes of arbitrary chirality. The developed force field faithf… ▽ More We introduce an interatomic potential for hexagonal boron nitride (hBN) based on the Gaussian approximation potential (GAP) machine learning methodology. The potential is based on a training set of configurations collected from density functional theory (DFT) simulations and is capable of treating bulk and multilayer hBN as well as nanotubes of arbitrary chirality. The developed force field faithfully reproduces the potential energy surface predicted by DFT while improving the efficiency by several orders of magnitude. We test our potential by comparing formation energies, geometrical properties, phonon dispersion spectra and mechanical properties with respect to benchmark DFT calculations and experiments. In addition, we use our model and a recently developed graphene-GAP to analyse and compare thermally and mechanically induced rippling in large scale two-dimensional (2D) hBN and graphene. Both materials show almost identical scaling behaviour with an exponent of $η\approx 0.85$ for the height fluctuations agreeing well with the theory of flexible membranes. Based on its lower resistance to bending, however, hBN experiences slightly larger out-of-plane deviations both at zero and finite applied external strain. Upon compression a phase transition from incoherent ripple motion to soliton-ripples is observed for both materials. Our potential is freely available online at [http://www.libatoms.org]. △ Less

Submitted 4 October, 2020; v1 submitted 22 July, 2020; originally announced July 2020.

Comments: This article has been accepted for publication in the Journal of Physical Chemistry C

Native point defects and low $p$-doping efficiency in $Mg_2 (Si,Sn)$ solid solutions: A hybrid-density functional study

Authors: Byungki Ryu, Eun-Ae Choi, Sungjin Park, Jaywan Chung, Johannes de Boor, Pawel Ziolkowski, Eckhard Müller, SuDong Park

Abstract: We perform hybrid-density functional calculations to investigate the charged defect formation energy of native point defects in $Mg_2 Si$, $Mg_2 Sn$, and their solid solutions. The band gap correction by hybrid-density functional is found to be critical to determine the charged defect density in these materials. For $Mg_2 Si$, $Mg$ interstitials are dominant and provide unintentional $n$-type cond… ▽ More We perform hybrid-density functional calculations to investigate the charged defect formation energy of native point defects in $Mg_2 Si$, $Mg_2 Sn$, and their solid solutions. The band gap correction by hybrid-density functional is found to be critical to determine the charged defect density in these materials. For $Mg_2 Si$, $Mg$ interstitials are dominant and provide unintentional $n$-type conductivity. Additionally, as the $Mg$ vacancies can dominate in $Mg$-poor $Mg_2 Sn$, $p$-type conductivity is possible for $Mg_2 Sn$. However, the existence of low formation energy defects such as $Mg_{Sn}^{1+}$ and $I_{Mg}^{2+}$ in $Mg_2 Sn$ and their diffusion can cause severe charge compensation of hole carriers resulting in low $p$-type doping efficiency and thermal degradation. Our results indicate that, in addition to the extrinsic doping strategy, alloying of $Mg_2 Si$ with $Mg_2 Sn$ under $Mg$-poor conditions would be necessary to enhance the $p$-type conductivity with less charge compensation. △ Less

Submitted 14 September, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: Main: 17 pages (including title, abstract, main, references, figure captions. 4 figures). This manuscript is accepted for publication in JALCOM. The article will be published as Gold Open Access

Journal ref: Journal of Alloys and Compounds 853, 157145 (2021)

Counterintuitive example on relation between ZT and thermoelectric efficiency

Authors: Byungki Ryu, Jaywan Chung, Eun-Ae Choi, Pawel Ziolkowski, Eckhard Müller, SuDong Park

Abstract: The thermoelectric figure of merit ZT, which is defined using electrical conductivity, Seebeck coefficient, thermal conductivity, and absolute temperature T, has been widely used as a simple estimator of the conversion efficiency of a thermoelectric heat engine. When material properties are constant or slowly varying with T, a higher ZT ensures a higher maximum conversion efficiency of thermoelect… ▽ More The thermoelectric figure of merit ZT, which is defined using electrical conductivity, Seebeck coefficient, thermal conductivity, and absolute temperature T, has been widely used as a simple estimator of the conversion efficiency of a thermoelectric heat engine. When material properties are constant or slowly varying with T, a higher ZT ensures a higher maximum conversion efficiency of thermoelectric materials. However, as material properties can vary strongly with T, efficiency predictions based on ZT can be inaccurate, especially for wide-temperature applications. Moreover, although ZT values continue to increase, there has been no investigation of the relationship between ZT and the efficiency in the higher ZT regime. In this paper, we report a counterintuitive situation by comparing two materials: although one material has a higher ZT value over the whole operational temperature range, its maximum conversion efficiency is smaller than that of the other. This indicates that, for material comparisons, the evaluation of exact efficiencies as opposed to a simple comparison of the ZTs is necessary in certain cases. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 12 pages, 2 tables, 2 figures

Journal ref: Appl. Phys. Lett. 116, 193903 (2020)

Dynamics of Bose-Einstein Recondensation in Higher Bands

Authors: Vaibhav Sharma, Sayan Choudhury, Erich J. Mueller

Abstract: Motivated by recent experiments, we explore the kinetics of Bose-Einstein condensation in the upper band of a double well optical lattice. These experiments engineer a non-equilibrium situation in which the highest energy state in the band is macroscopically occupied. The system subsequently relaxes and the condensate moves to the lowest energy state. We model this process, finding that the kineti… ▽ More Motivated by recent experiments, we explore the kinetics of Bose-Einstein condensation in the upper band of a double well optical lattice. These experiments engineer a non-equilibrium situation in which the highest energy state in the band is macroscopically occupied. The system subsequently relaxes and the condensate moves to the lowest energy state. We model this process, finding that the kinetics occurs in three phases: The condensate first evaporates, forming a highly non-equilibrium gas with no phase coherence. Energy is then redistributed among the noncondensed atoms. Finally the atoms recondense. We calculate the time-scales for each of these phases, and explain how this scenario can be verified through future experiments. △ Less

Submitted 20 February, 2020; v1 submitted 12 November, 2019; originally announced November 2019.

Comments: 12 Pages, 4 Figures; Revisions in response to referee comments and corrected typographical errors

Journal ref: Phys. Rev. A 101, 033609 (2020)

Exact Topological Flat Bands from Continuum Landau Levels

Authors: Junkai Dong, Erich Mueller

Abstract: We construct and characterize tight binding Hamiltonians which contain a completely flat topological band made of continuum lowest Landau level wavefunctions sampled on a lattice. We find an infinite family of such Hamiltonians, with simple analytic descriptions. These provide a valuable tool for constructing exactly solvable models. We also implement a numerical algorithm for finding the most loc… ▽ More We construct and characterize tight binding Hamiltonians which contain a completely flat topological band made of continuum lowest Landau level wavefunctions sampled on a lattice. We find an infinite family of such Hamiltonians, with simple analytic descriptions. These provide a valuable tool for constructing exactly solvable models. We also implement a numerical algorithm for finding the most local Hamiltonian with a flat Landau level. We find intriguing structures in the spatial dependence of the matrix elements for this optimized model. The models we construct serve as foundations for numerical and experimental studies of topological systems, both non-interacting and interacting. △ Less

Submitted 17 December, 2019; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. A 101, 013629 (2020)

Thermal transport of helium-3 in a strongly confining channel

Authors: D. Lotnyk, A. Eyal, N. Zhelev, T. S. Abhilash, E. N. Smith, M. Terilli, J. Wilson, E. Mueller, D. Einzel, J. Saunders, J. M. Parpia

Abstract: In a neutral system such as liquid helium-3, transport of mass, heat, and spin provide information analogous to electrical counterparts in metals, superconductors and topological materials. Of particular interest is transport in strongly confining channels of height approaching the superfluid coherence length, where new quantum states are found and excitations bound to surfaces and edges should be… ▽ More In a neutral system such as liquid helium-3, transport of mass, heat, and spin provide information analogous to electrical counterparts in metals, superconductors and topological materials. Of particular interest is transport in strongly confining channels of height approaching the superfluid coherence length, where new quantum states are found and excitations bound to surfaces and edges should be present. Here we report on the thermal conduction of helium-3 in a 1.1~$μ$m high microfabricated channel. In the normal state we observe a diffusive thermal conductivity that is approximately temperature independent, consistent with recent work on the interference of bulk and boundary scattering. In the superfluid state we measure diffusive thermal transport in the absence of thermal counterflow. An anomalous thermal response is also detected in the superfluid which we suggest may arise from a flux of surface excitations. △ Less

Submitted 18 October, 2019; originally announced October 2019.

Comments: A supplement is available. Please contact Jeevak Parpia (jmp9@cornell.edu) if you would like the supplement as well

Spreading of aqueous droplets with common and superspreading surfactants. A molecular dynamics study

Authors: Panagiotis E. Theodorakis, Edward R. Smith, Erich A. Müller

Abstract: The surfactant-driven spreading of droplets is an essential process in many applications ranging from coating flow technology to enhanced oil recovery. Despite the significant advancement in describing spreading processes in surfactant-laden droplets, including the exciting phenomena of superspreading, many features of the underlying mechanisms require further understanding. Here, we have carried… ▽ More The surfactant-driven spreading of droplets is an essential process in many applications ranging from coating flow technology to enhanced oil recovery. Despite the significant advancement in describing spreading processes in surfactant-laden droplets, including the exciting phenomena of superspreading, many features of the underlying mechanisms require further understanding. Here, we have carried out molecular dynamics simulations of a coarse-grained model with a force-field obtained from the statistical associating fluid theory to study droplets laden with common and superspreading surfactants. We have confirmed the important elements of the superspreading mechanism, i.e. the adsorption of surfactant at the contact line and the fast replenishment of surfactant from the bulk. Through a detailed comparison of a range of droplets with different surfactants, our analysis has indicated that the ability of surfactant to adsorb at the interfaces is the key feature of the superspreading mechanism. To this end, surfactants that tend to form aggregates and have a strong hydrophobic attraction in the aggregated cores prevent the fast replenishment of the interfaces, resulting in reduced spreading ability. We also show that various surfactant properties, such as diffusion and architecture, play a secondary role in the spreading process. Moreover, we discuss various drop properties such as the height, contact angle, and surfactant surface concentration, highlighting differences between superspreading and common surfactants. We anticipate that our study will provide further insight for applications requiring the control of wetting. △ Less

Submitted 2 September, 2019; originally announced September 2019.

Comments: 26 pages, 11 figures

Journal ref: Colloids and Surfaces A: Physicochemical and Engineering Aspects 581, 123810 (2019)