Showing 1–50 of 135 results for author: Alvarez, G

Magnetic phase diagram of a two-orbital model for bilayer nickelates varying doping

Authors: Ling-Fang Lin, Yang Zhang, Nitin Kaushal, Gonzalo Alvarez, Thomas A. Maier, Adriana Moreo, Elbio Dagotto

Abstract: Motivated by the recently discovered high-$T_c$ bilayer nickelate superconductor La$_3$Ni$_2$O$_7$, we comprehensively research a bilayer $2\times2\times2$ cluster for different electronic densities $n$ by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying $n$. Based on the spin… ▽ More Motivated by the recently discovered high-$T_c$ bilayer nickelate superconductor La$_3$Ni$_2$O$_7$, we comprehensively research a bilayer $2\times2\times2$ cluster for different electronic densities $n$ by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying $n$. Based on the spin structure factor $S(q)$, we have obtained a rich magnetic phase diagram in the plane defined by $n$ and $U/W$, at fixed Hund coupling. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. For half-filling $n=2$ (two electrons per Ni site, corresponding to $N$ = 16 electrons), the canonical superexchange interaction leads to a robust G-AFM state $(π,π,π)$ with antiferromagnetic couplings in plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the ``half-empty'' and ``half-full'' mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At $n = 1.5$ (or $N=12$), density corresponding to La$_3$Ni$_2$O$_7$, we obtained the ``Stripe 2'' ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the $z$-axis) in the $2\times2\times2$ cluster. In addition, we obtained a much stronger AFM coupling along the $z$-axis than the magnetic coupling in the $xy$ plane. The random-phase approximation calculations with varying $n$ give very similar results as Lanczos. Meanwhile, a state with $q/π= (0.6, 0.6, 1)$ close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to $n=1.25$, possibly responsible for the E-phase SDW recently observed in experiments. △ Less

Submitted 10 August, 2024; originally announced August 2024.

Comments: 12 pages, 9 figures

Thermal transport of flexural phonons in a rectangular wire

Authors: G. Rivas Álvarez, E. Benítez Rodríguez, M. A. Bastarrachea-Magnani, M. Martínez-Mares, R. A. Méndez-Sánchez

Abstract: We study the quantum thermal transport of elastic excitations through a two-dimensional elastic waveguide between two thermal reservoirs. We solve the classical Kirchhoff-Love equation for rectangular wires and explore the dispersion relation for both the symmetric and antisymmetric solutions. Then, we study the phonon transport of these modes within the second quantization framework by analyzing… ▽ More We study the quantum thermal transport of elastic excitations through a two-dimensional elastic waveguide between two thermal reservoirs. We solve the classical Kirchhoff-Love equation for rectangular wires and explore the dispersion relation for both the symmetric and antisymmetric solutions. Then, we study the phonon transport of these modes within the second quantization framework by analyzing the mean quadratic displacement, from which the energy density current, the temperature field, and conductance are determined. We identify the relevant modes contributing to thermal transport and explore the average temperature difference to reach the high-temperature limit. We expect our results to pave the way for understanding phonon-mediated thermal transport in two-dimensional mesoscopic quantum devices. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: Submitted to Journal of Applied Physics

Crystalline-Symmetry-Protected Majorana Modes in Coupled Quantum Dots

Authors: Bradraj Pandey, Gonzalo Alvarez, Elbio Dagotto, Rui-Xing Zhang

Abstract: We propose a minimalist architecture for achieving various crystalline-symmetry-protected Majorana modes in an array of coupled quantum dots. Our framework is motivated by the recent experimental demonstrations of two-site and three-site artificial Kitaev chains in a similar setup. We find that introducing a $π$-phase domain wall in the Kitaev chain leads to a pair of mirror-protected Majorana zer… ▽ More We propose a minimalist architecture for achieving various crystalline-symmetry-protected Majorana modes in an array of coupled quantum dots. Our framework is motivated by the recent experimental demonstrations of two-site and three-site artificial Kitaev chains in a similar setup. We find that introducing a $π$-phase domain wall in the Kitaev chain leads to a pair of mirror-protected Majorana zero modes located at or near the junction. Joining two $π$-junctions into a closed loop, we can simulate two distinct classes of two-dimensional higher-order topological superconducting phases, both carrying symmetry-protected Majorana modes around the sample corners. As an extension of the $π$-junction, we further consider a general vertex structure where $n$ Kitaev chains meet, i.e., a Kitaev $n$-vertex. We prove that such an $n$-vertex, if respecting a dihedral symmetry group $D_n$, necessarily carries $n$ vertex-bound Majorana modes protected by the $D_n$ symmetry. Resilience of the junction and vertex Majorana bound states against disorder and correlation effects is also discussed. Our architecture paves the way for designing, constructing, and exploring a wide variety of artificial topological crystalline phases in experiments. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 8 pages, 5 figures

Anomalous properties of spark plasma sintered boron nitride solids

Authors: Abhijit Biswas, Peter Serles, Gustavo A. Alvarez, Jesse Schimpf, Michel Hache, Jonathan Kong, Pedro Guerra Demingos, Bo Yuan, Tymofii S. Pieshkov, Chenxi Li, Anand B. Puthirath, Bin Gao, Tia Gray, Xiang Zhang, Jishnu Murukeshan, Robert Vajtai, Pengcheng Dai, Chandra Veer Singh, Jane Howe, Yu Zou, Lane W. Martin, James Patrick Clancy, Zhiting Tian, Tobin Filleter, Pulickel M. Ajayan

Abstract: Hexagonal boron nitride (h-BN) is brittle, however, its atomic-scale structural engineering can lead to unprecedented physical properties. Here we report the bulk synthesis of high-density crystalline h-BN solids by using high-temperature spark plasma sintering (SPS) of micron size h-BN powders. In addition to the high mechanical strength and ductile response of such materials, we have obtained an… ▽ More Hexagonal boron nitride (h-BN) is brittle, however, its atomic-scale structural engineering can lead to unprecedented physical properties. Here we report the bulk synthesis of high-density crystalline h-BN solids by using high-temperature spark plasma sintering (SPS) of micron size h-BN powders. In addition to the high mechanical strength and ductile response of such materials, we have obtained anomalous values of dielectric constant beyond theoretical limits, high thermal conductivity, and exceptional neutron radiation shielding capability. Through exhaustive characterizations we reveal that SPS induces non-basal plane crystallinity, twisting of layers, and facilitates inter-grain fusion with a high degree of in-plane alignment across macroscale dimensions, resulting in near-theoretical density and anomalous properties. Our findings highlight the importance of material design, via new approaches such as twisting and interconnections between atomically thin layers, to create novel ceramics with properties that could go beyond their intrinsic theoretical predictions. △ Less

Submitted 10 July, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: Authors revised version, 46 pages, 4 figures

Sensing Out-of-Equilibrium and Quantum Non-Gaussian environments via induced Time-Reversal Symmetry Breaking on the quantum-probe dynamics

Authors: Martin Kuffer, Analía Zwick, Gonzalo A. Álvarez

Abstract: Advancing quantum sensing tools for investigating systems at atomic and nanoscales is crucial for the progress of quantum technologies. While numerous protocols employ quantum probes to extract information from stationary or weakly coupled environments, the challenges intensify at atomic- and nano-scales where the environment is inherently out-of-equilibrium or strongly coupled with the sensor. We… ▽ More Advancing quantum sensing tools for investigating systems at atomic and nanoscales is crucial for the progress of quantum technologies. While numerous protocols employ quantum probes to extract information from stationary or weakly coupled environments, the challenges intensify at atomic- and nano-scales where the environment is inherently out-of-equilibrium or strongly coupled with the sensor. We here prove that the time-reversal symmetry in the quantum-sensor control dynamics is broken, when partial information is probed from an environment that is out-of-equilibrium with non stationary fluctuations or is described by quantum non-Gaussian, strongly coupled environmental correlations. We exploit this phenomenon as a quantum sensing paradigm with proof-of principle experimental quantum simulations using solid-state nuclear magnetic resonance (NMR). This introduces a signal contrast on a qubit-probe that quantifies how far the sensed environment is from equilibrium or its quantum non-Gaussian nature. Protocols are also presented to discern and filter a variety of environmental properties including stationary, non-stationary and non-Gaussian quantum noise fluctuations as a step toward sensing the ubiquitous environments of a quantum-sensor at atomic and nanoscales. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Maximizing information obtainable by quantum sensors through the Quantum Zeno Effect

Authors: Bruno Ronchi, Analia Zwick, Gonzalo A. Alvarez

Abstract: Efficient quantum sensing technologies rely on precise control of quantum sensors, particularly two-level systems or qubits, to optimize estimation processes. We here exploit the Quantum Zeno Effect (QZE) as a tool for maximizing information obtainable by quantum sensors, with a specific focus on the level avoided crossing (LAC) phenomenon in qubit systems. While the estimation of the LAC energy s… ▽ More Efficient quantum sensing technologies rely on precise control of quantum sensors, particularly two-level systems or qubits, to optimize estimation processes. We here exploit the Quantum Zeno Effect (QZE) as a tool for maximizing information obtainable by quantum sensors, with a specific focus on the level avoided crossing (LAC) phenomenon in qubit systems. While the estimation of the LAC energy splitting has been extensively studied, we emphasize the crucial role that the QZE can play in estimating the coupling strength. We introduce the concept of information amplification by the QZE for a LAC system under off-resonant conditions. The proposed approach has implications for AC magnetic field sensing and the caracterization of complex systems, including many-spin systems requiring the estimation of spin-spin couplings. Overall, our findings contribute to the advancement of quantum sensing by leveraging the QZE for improved control and information extraction. △ Less

Submitted 17 March, 2024; originally announced March 2024.

Comments: 14 pages, 6 figures

Computational supremacy in quantum simulation

Authors: Andrew D. King, Alberto Nocera, Marek M. Rams, Jacek Dziarmaga, Roeland Wiersema, William Bernoudy, Jack Raymond, Nitin Kaushal, Niclas Heinsdorf, Richard Harris, Kelly Boothby, Fabio Altomare, Andrew J. Berkley, Martin Boschnak, Kevin Chern, Holly Christiani, Samantha Cibere, Jake Connor, Martin H. Dehn, Rahul Deshpande, Sara Ejtemaee, Pau Farré, Kelsey Hamer, Emile Hoskinson, Shuiyuan Huang , et al. (37 additional authors not shown)

Abstract: Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. Establishing this capability, especially for impactful and meaningful problems, remains a central challenge. One such problem is the simulation of nonequilibrium dynamics of a magnetic spin system quenched through a quantum phase transition. State-of-the-art classical simulations dem… ▽ More Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. Establishing this capability, especially for impactful and meaningful problems, remains a central challenge. One such problem is the simulation of nonequilibrium dynamics of a magnetic spin system quenched through a quantum phase transition. State-of-the-art classical simulations demand resources that grow exponentially with system size. Here we show that superconducting quantum annealing processors can rapidly generate samples in close agreement with solutions of the Schrödinger equation. We demonstrate area-law scaling of entanglement in the model quench in two-, three- and infinite-dimensional spin glasses, supporting the observed stretched-exponential scaling of effort for classical approaches. We assess approximate methods based on tensor networks and neural networks and conclude that no known approach can achieve the same accuracy as the quantum annealer within a reasonable timeframe. Thus quantum annealers can answer questions of practical importance that classical computers cannot. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Block Mott insulating state induced by next-nearest neighbor hopping in the S = 3/2 zigzag chain BaCoTe2O7

Authors: Ling-Fang Lin, Yang Zhang, Gonzalo Alvarez, Adriana Moreo, Elbio Dagotto

Abstract: Quasi-one-dimensional correlated electronic multi-orbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-fields splitting, the electronic correlations, and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe$_2$O$_7$, wit… ▽ More Quasi-one-dimensional correlated electronic multi-orbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-fields splitting, the electronic correlations, and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe$_2$O$_7$, with electronic density $n = 7$, was prepared experimentally. Here, we systematically study the electronic and magnetic properties of this quasi-one-dimensional compound from the theory perspective. Based on first-principles density functional theory calculations, strongly anisotropic one-dimensional electronic Co $3d$ bands were found near the Fermi level. By evaluating the relevant hopping amplitudes, we provide the magnitude and origin of the nearest-neighbor (NN) and next nearest-neighbor (NNN) hopping matrices in BaCoTe$_2$O$_7$. With this information, we constructed a three-orbital electronic Hubbard model for this zigzag chain system, and studied two cases: with only a NN hopping matrix, and with NN plus NNN hopping matrices. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group method, we constructed the ground-state phase diagram. A robust staggered antiferromagnetic (AFM) region was found when only the NN hopping matrix in the chain direction was employed. However, for the realistic case where the NNN hopping matrix is also included, the dominant state becomes instead a block AFM order, in agreement with experiments. The system displays Mott insulator characteristics with three half-filled orbitals, when the block AFM order is stable. Our results for BaCoTe$_2$O$_7$ provide guidance to experimentalists and theorists working on this zigzag one-dimensional chain and related materials. △ Less

Submitted 6 May, 2024; v1 submitted 7 February, 2024; originally announced February 2024.

Comments: 9 pages, 10 figures. arXiv admin note: text overlap with arXiv:2106.02753

Luther-Emery liquid and dominant singlet superconductivity in the hole-doped Haldane spin-1 chain

Authors: Pontus Laurell, Jacek Herbrych, Gonzalo Alvarez, Elbio Dagotto

Abstract: We investigate the pairing tendencies in the hole-doped Haldane spin-1 chain. To allow for doping, we extend the original spin chain Hamiltonian into a fermionic model involving a two-orbital Hubbard chain at intermediate or strong repulsive interaction strengths $U$, and for degenerate orbitals. At half-filling and large $U$, the ferromagnetic Hund's coupling, $J_\mathrm{H}$, generates effective… ▽ More We investigate the pairing tendencies in the hole-doped Haldane spin-1 chain. To allow for doping, we extend the original spin chain Hamiltonian into a fermionic model involving a two-orbital Hubbard chain at intermediate or strong repulsive interaction strengths $U$, and for degenerate orbitals. At half-filling and large $U$, the ferromagnetic Hund's coupling, $J_\mathrm{H}$, generates effective spin-$1$ moments, with antiferromagnetic correlations between sites. Using large-scale density matrix renormalization group calculations, we study accurately the system's behavior under light hole-doping. For $U=1.6$ in units of the non-interacting bandwidth and for $J_\mathrm{H}/U\gtrsim 0.275$ we find that singlet pairing dominates the long-distance physics, establishing this system as a promising platform for repulsively mediated superconductivity. We provide concrete examples of materials that could realize the physics described here. We also provide evidence that the system approaches a Luther-Emery liquid state at large system sizes, reminiscent of the behavior of doped one-orbital two-leg ladders at weak coupling, which also have superconducting tendencies. The numerically calculated central charge approaches one in the thermodynamic limit, indicating a single gapless mode as is expected for the Luther-Emery state. Exponents characterizing the power-law decays of singlet pair-pair and charge density-density correlations are determined, and found to approximately satisfy the Luther-Emery identity. △ Less

Submitted 27 August, 2024; v1 submitted 22 November, 2023; originally announced November 2023.

Comments: Main text: 10 pages, 9 figures. Supplementary: 4 pages, 2 figures

Journal ref: Phys. Rev. B 110, 064515 (2024)

Topological and magnetic properties of the interacting Bernevig-Hughes-Zhang model

Authors: Rahul Soni, Harini Radhakrishnan, Bernd Rosenow, Gonzalo Alvarez, Adrian Del Maestro

Abstract: We investigate the effects of electronic correlations on the Bernevig-Hughes-Zhang model using the real-space density matrix renormalization group (DMRG) algorithm. We introduce a method to probe topological phase transitions in systems with strong correlations using DMRG, substantiated by an unsupervised machine learning methodology that analyzes the orbital structure of the real-space edges. Inc… ▽ More We investigate the effects of electronic correlations on the Bernevig-Hughes-Zhang model using the real-space density matrix renormalization group (DMRG) algorithm. We introduce a method to probe topological phase transitions in systems with strong correlations using DMRG, substantiated by an unsupervised machine learning methodology that analyzes the orbital structure of the real-space edges. Including the full multi-orbital Hubbard interaction term, we construct a phase diagram as a function of a gap parameter ($m$) and the Hubbard interaction strength ($U$) via exact DMRG simulations on $N\times 4$ cylinders. Our analysis confirms that the topological phase persists in the presence of interactions, consistent with previous studies, but it also reveals an intriguing phase transition from a paramagnetic to an antiferromagnetic topological insulator. The combination of the magnetic structure factor, strength of magnetic moments, and the orbitally resolved density, provides real-space information on both topology and magnetism in a strongly correlated system. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: Main text: 11 pages, 12 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-interactingBHZmodel

Spinon continuum in the Heisenberg quantum chain compound Sr$_2$V$_3$O$_9$

Authors: Shang Gao, Ling-Fang Lin, Pontus Laurell, Qiang Chen, Qing Huang, Clarina dela Cruz, Krishnamurthy V. Vemuru, Mark D. Lumsden, Stephen E. Nagler, Gonzalo Alvarez, Elbio Dagotto, Haidong Zhou, Andrew D. Christianson, Matthew B. Stone

Abstract: Magnetic excitations in the spin chain candidate Sr$_2$V$_3$O$_9$ have been investigated by inelastic neutron scattering on a single crystal sample. A spinon continuum with a bandwidth of $\sim22$ meV is observed along the chain formed by alternating magnetic V$^{4+}$ and nonmagnetic V$^{5+}$ ions. Incipient magnetic Bragg peaks due to weak ferromagnetic interchain couplings emerge when approachin… ▽ More Magnetic excitations in the spin chain candidate Sr$_2$V$_3$O$_9$ have been investigated by inelastic neutron scattering on a single crystal sample. A spinon continuum with a bandwidth of $\sim22$ meV is observed along the chain formed by alternating magnetic V$^{4+}$ and nonmagnetic V$^{5+}$ ions. Incipient magnetic Bragg peaks due to weak ferromagnetic interchain couplings emerge when approaching the magnetic transition at $T_N\sim 5.3$ K while the excitations remain gapless within the instrumental resolution. Comparisons to the Bethe ansatz, density matrix renormalization group (DMRG) calculations, and effective field theories confirm Sr$_2$V$_3$O$_9$ as a host of weakly coupled $S = 1/2$ chains dominated by antiferromagnetic intrachain interactions of $\sim7.1$(1) meV. △ Less

Submitted 25 July, 2023; v1 submitted 22 July, 2023; originally announced July 2023.

Comments: 19 pages, 9 figures

Journal ref: Phys. Rev. B 109, L020402 (2024)

Majorana zero modes in Y-shape interacting Kitaev wires

Authors: Bradraj Pandey, Nitin Kaushal, Gonzalo Alvarez, Elbio Dagotto

Abstract: Motivated by the recent experimental realization of minimal Kitaev chains using quantum dots, we investigate the Majorana zero modes (MZM) in $Y$-shape Kitaev wires. We solve the associated Kitaev models analytically at the sweet spot ($t_h=Δ$) and derive the exact form of MZM wave-functions in this geometry. The novelty of our result is the observation of multi-site MZMs located near the junction… ▽ More Motivated by the recent experimental realization of minimal Kitaev chains using quantum dots, we investigate the Majorana zero modes (MZM) in $Y$-shape Kitaev wires. We solve the associated Kitaev models analytically at the sweet spot ($t_h=Δ$) and derive the exact form of MZM wave-functions in this geometry. The novelty of our result is the observation of multi-site MZMs located near the junction center on the nearby edge sites of each leg. This result is important for potential braiding of Majoranas and the performance of $Y$-junctions made from arrays of quantum dots. Furthermore, we study the stability of local (single-site) and multi-site MZMs modes in the presence of Coulomb repulsion, using density matrix renormalization group theory. Our local density-of-states calculation shows that these multi-site MZMs are as equally topologically protected as the single-site MZMs when in the presence of Coulomb repulsion or when away from the sweet-spot. △ Less

Submitted 11 August, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: 8 pages and 5 figures

Journal ref: npj Quantum Mater. 8, 51 (2023)

Structural, optical, and thermal properties of BN thin films grown on diamond via pulsed laser deposition

Authors: Abhijit Biswas, Gustavo A. Alvarez, Tao Li, Joyce Christiansen-Salameh, Eugene Jeong, Anand B. Puthirath, Sathvik Ajay Iyengar, Chenxi Li, Tia Gray, Xiang Zhang, Tymofii S. Pieshkov, Harikishan Kannan, Jacob Elkins, Robert Vajtai, A. Glen Birdwell, Mahesh R. Neupane, Elias J. Garratt, Bradford B. Pate, Tony G. Ivanov, Yuji Zhao, Zhiting Tian, Pulickel M. Ajayan

Abstract: Heterostructures based on ultrawide-bandgap (UWBG) semiconductors (bandgap >4.0 eV), boron nitride (BN) and diamond are important for next-generation high-power electronics. However, in-situ hetero-epitaxy of BN/diamond or vice-versa remains extremely challenging, due to their non-trivial growth kinetics. Here, we have grown BN thin film on (100) single crystal diamond by pulsed laser deposition a… ▽ More Heterostructures based on ultrawide-bandgap (UWBG) semiconductors (bandgap >4.0 eV), boron nitride (BN) and diamond are important for next-generation high-power electronics. However, in-situ hetero-epitaxy of BN/diamond or vice-versa remains extremely challenging, due to their non-trivial growth kinetics. Here, we have grown BN thin film on (100) single crystal diamond by pulsed laser deposition and investigated its structural and magnetic properties, optical refractive index, and thermal conductivity. Structural characterizations confirm the mixed (stable hexagonal and metastable cubic) phase growth. Film shows diamagnetic behavior at room temperature. It displays anisotropic refractive index within the visible-to-near-infrared wavelength range. The room temperature cross-plane thermal conductivity of BN is ~1.53 W/(mK), and the thermal conductance of the BN/diamond interface is ~20 MW/(m2K). Our findings are useful for various device related applications based on UWBG BN/diamond heterostructures. △ Less

Submitted 20 September, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: 16 pages, 4 figures

Journal ref: Phys. Rev. Materials 7, 094602 (2023)

Non-linear optics at twist interfaces in h-BN/SiC heterostructures

Authors: Abhijit Biswas, Rui Xu, Gustavo A. Alvarez, Jin Zhang, Joyce Christiansen-Salameh, Anand B. Puthirath, Kory Burns, Jordan A. Hachtel, Tao Li, Sathvik Ajay Iyengar, Tia Gray, Chenxi Li, Xiang Zhang, Harikishan Kannan, Jacob Elkins, Tymofii S. Pieshkov, Robert Vajtai, A. Glen Birdwell, Mahesh R. Neupane, Elias J. Garratt, Tony Ivanov, Bradford B. Pate, Yuji Zhao, Hanyu Zhu, Zhiting Tian , et al. (2 additional authors not shown)

Abstract: Understanding the emergent electronic structure in twisted atomically thin layers has led to the exciting field of twistronics. However, practical applications of such systems are challenging since the specific angular correlations between the layers must be precisely controlled and the layers have to be single crystalline with uniform atomic ordering. Here, we suggest an alternative, simple and s… ▽ More Understanding the emergent electronic structure in twisted atomically thin layers has led to the exciting field of twistronics. However, practical applications of such systems are challenging since the specific angular correlations between the layers must be precisely controlled and the layers have to be single crystalline with uniform atomic ordering. Here, we suggest an alternative, simple and scalable approach where nanocrystalline two-dimensional (2D) film on three-dimensional (3D) substrates yield twisted-interface-dependent properties. Ultrawide-bandgap hexagonal boron nitride (h-BN) thin films are directly grown on high in-plane lattice mismatched wide-bandgap silicon carbide (4H-SiC) substrates to explore the twist-dependent structure-property correlations. Concurrently, nanocrystalline h-BN thin film shows strong non-linear second-harmonic generation and ultra-low cross-plane thermal conductivity at room temperature, which are attributed to the twisted domain edges between van der Waals stacked nanocrystals with random in-plane orientations. First-principles calculations based on time-dependent density functional theory manifest strong even-order optical nonlinearity in twisted h-BN layers. Our work unveils that directly deposited 2D nanocrystalline thin film on 3D substrates could provide easily accessible twist-interfaces, therefore enabling a simple and scalable approach to utilize the 2D-twistronics integrated in 3D material devices for next-generation nanotechnology. △ Less

Submitted 4 November, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Comments: 33 pages, 4 figures

Journal ref: Adv. Mater. 2023, 2304624

Transition to the Haldane phase driven by electron-electron correlations

Authors: A. Jażdżewska, M. Mierzejewski, M. Środa, A. Nocera, G. Alvarez, E. Dagotto, J. Herbrych

Abstract: One of the most famous quantum systems with topological properties, the spin $\mathcal{S}=1$ antiferromagnetic Heisenberg chain, is well-known to display exotic $\mathcal{S}=1/2$ edge states. However, this spin model has not been analyzed from the more general perspective of strongly correlated systems varying the electron-electron interaction strength. Here, we report the investigation of the eme… ▽ More One of the most famous quantum systems with topological properties, the spin $\mathcal{S}=1$ antiferromagnetic Heisenberg chain, is well-known to display exotic $\mathcal{S}=1/2$ edge states. However, this spin model has not been analyzed from the more general perspective of strongly correlated systems varying the electron-electron interaction strength. Here, we report the investigation of the emergence of the Haldane edge in a system of interacting electrons -- the two-orbital Hubbard model -- with increasing repulsion strength $U$ and Hund interaction $J_\mathrm{H}$. We show that interactions not only form the magnetic moments but also form a topologically nontrivial fermionic many-body ground-state with zero-energy edge states. Specifically, upon increasing the strength of the Hubbard repulsion and Hund exchange, we identify a sharp transition point separating topologically trivial and nontrivial ground-states. Surprisingly, such a behaviour appears already at rather small values of the interaction, in a regime where the magnetic moments are barely developed. △ Less

Submitted 22 December, 2023; v1 submitted 21 April, 2023; originally announced April 2023.

Journal ref: Nat. Commun. 14, 8524 (2023)

Phase Stability of Hexagonal/cubic Boron Nitride Nanocomposites

Authors: Abhijit Biswas, Rui Xu, Joyce Christiansen-Salameh, Eugene Jeong, Gustavo A. Alvarez, Chenxi Li, Anand B. Puthirath, Bin Gao, Arushi Garg, Tia Gray, Harikishan Kannan, Xiang Zhang, Jacob Elkins, Tymofii S. Pieshkov, Robert Vajtai, A. Glen Birdwell, Mahesh R. Neupane, Bradford B. Pate, Tony Ivanov, Elias J. Garratt, Pengcheng Dai, Hanyu Zhu, Zhiting Tian, Pulickel M. Ajayan

Abstract: Boron nitride (BN) is an exceptional material and among its polymorphs, two-dimensional (2D) hexagonal and three-dimensional (3D) cubic BN (h-BN and c-BN) phases are most common. The phase stability regimes of these BN phases are still under debate and phase transformations of h-BN/c-BN remain a topic of interest. Here, we investigate the phase stability of 2D/3D h-BN/c-BN nanocomposites and show… ▽ More Boron nitride (BN) is an exceptional material and among its polymorphs, two-dimensional (2D) hexagonal and three-dimensional (3D) cubic BN (h-BN and c-BN) phases are most common. The phase stability regimes of these BN phases are still under debate and phase transformations of h-BN/c-BN remain a topic of interest. Here, we investigate the phase stability of 2D/3D h-BN/c-BN nanocomposites and show that the co-existence of two phases can lead to strong non-linear optical properties and low thermal conductivity at room temperature. Furthermore, spark-plasma sintering of the nanocomposite shows complete phase transformation to 2D h-BN with improved crystalline quality, where 3D c-BN grain sizes governs the nucleation and growth kinetics. Our demonstration might be insightful in phase engineering of BN polymorphs based nanocomposites with desirable properties for optoelectronics and thermal energy management applications. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 29 pages, 5 figures

Journal ref: Nano Lett. 2023, 23, 15, 6927

Stability of the novel interorbital-hopping mechanism for ferromagnetism in multi-orbital Hubbard models

Authors: Ling-Fang Lin, Yang Zhang, Gonzalo Alvarez, Michael A. McGuire, Andrew F. May, Adriana Moreo, Elbio Dagotto

Abstract: Recently, it was argued that a ferromagnetic (FM) insulating phase can be induced by a novel {\it interorbital} hopping mechanism. Here, we study the stability range of this novel FM phase under modifications in the crystal fields and electronic correlation strength, constructing a theoretical phase diagram. A plethora of states is unveiled, including the FM Mott insulator (MI), a FM orbital-selec… ▽ More Recently, it was argued that a ferromagnetic (FM) insulating phase can be induced by a novel {\it interorbital} hopping mechanism. Here, we study the stability range of this novel FM phase under modifications in the crystal fields and electronic correlation strength, constructing a theoretical phase diagram. A plethora of states is unveiled, including the FM Mott insulator (MI), a FM orbital-selective Mott phase (OSMP), several anferromagnetic (AFM) MI phases, an AFM metallic state, and a FM metal as well. Our most interesting result is that the FM regime, either in MI or OSMP forms, is shown to be stable in {\it large} portions of the phase diagram, at both intermediate and strong electronic correlations, respectively. Our results demonstrate via a detailed example that the recently proposed novel mechanism to stabilize FM insulators is not fragile but instead robust, and may enlarge substantially the relatively small family of known FM insulators. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 10 pages, 5 figures

Gene Expression Programming for Quantum Computing

Authors: Gonzalo Alvarez, Ryan Bennink, Stephan Irle, Jacek Jakowski

Abstract: We introduce QuantumGEP, a scientific computer program that uses gene expression programming (GEP) to find a quantum circuit that either (i) maps a given set of input states to a given set of output states, or (ii) transforms a fixed initial state to minimize a given physical quantity of the output state. QuantumGEP is a driver program that uses evendim, a generic computational engine for GEP, bot… ▽ More We introduce QuantumGEP, a scientific computer program that uses gene expression programming (GEP) to find a quantum circuit that either (i) maps a given set of input states to a given set of output states, or (ii) transforms a fixed initial state to minimize a given physical quantity of the output state. QuantumGEP is a driver program that uses evendim, a generic computational engine for GEP, both of which are free and open source. We apply QuantumGEP as a powerful solver for MaxCut in graphs, and for condensed matter quantum many-body Hamiltonians. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 10 pages, 5 figures

Hund bands in spectra of multiorbital systems

Authors: M. Środa, J. Mravlje, G. Alvarez, E. Dagotto, J. Herbrych

Abstract: Spectroscopy experiments are routinely used to characterize the behavior of strongly correlated systems. An in-depth understanding of the different spectral features is thus essential. Here, we show that the spectrum of the multiorbital Hubbard model exhibits unique Hund \ms{bands} that occur at energies given only by the Hund coupling $J_\mathrm{H}$, as distinct from the Hubbard satellites follow… ▽ More Spectroscopy experiments are routinely used to characterize the behavior of strongly correlated systems. An in-depth understanding of the different spectral features is thus essential. Here, we show that the spectrum of the multiorbital Hubbard model exhibits unique Hund \ms{bands} that occur at energies given only by the Hund coupling $J_\mathrm{H}$, as distinct from the Hubbard satellites following the interaction $U$. We focus on experimentally relevant single-particle and optical spectra that we calculate for a model related to iron chalcogenide ladders. The calculations are performed via the density-matrix renormalization group and Lanczos methods. The generality of the implications is verified by considering a generic multiorbital model within dynamical mean-field theory. △ Less

Submitted 3 August, 2023; v1 submitted 20 October, 2022; originally announced October 2022.

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

Spin dynamics of the generalized quantum spin compass chain

Authors: Pontus Laurell, Gonzalo Alvarez, Elbio Dagotto

Abstract: We calculate the dynamical spin structure factor of the generalized spin-$1/2$ compass spin chain using the density matrix renormalization group. The model, also known as the twisted Kitaev spin chain, was recently proposed to be relevant for the description of the spin chain compound CoNb$_2$O$_6$. It features bond-dependent interactions and interpolates between an Ising chain and a one-dimension… ▽ More We calculate the dynamical spin structure factor of the generalized spin-$1/2$ compass spin chain using the density matrix renormalization group. The model, also known as the twisted Kitaev spin chain, was recently proposed to be relevant for the description of the spin chain compound CoNb$_2$O$_6$. It features bond-dependent interactions and interpolates between an Ising chain and a one-dimensional variant of Kitaev's honeycomb spin model. The structure factor, in turn, is found to interpolate from gapped and non-dispersive in the Ising limit to gapless with non-trivial continua in the Kitaev limit. In particular, the component of the structure factor perpendicular to the Ising directions changes abruptly at the Kitaev point into a dispersionless continuum due to the emergence of an extensive groundstate degeneracy. We show this continuum is consistent with analytical Jordan-Wigner results. We also discuss implications for future inelastic scattering experiments and applications to materials, particularly CoNb$_2$O$_6$. △ Less

Submitted 27 March, 2023; v1 submitted 1 October, 2022; originally announced October 2022.

Comments: Main text: 10 pages, 6 figures. Supplemental material: 5 pages, 2 figures

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

Root-N Krylov-space correction-vectors for spectral functions with the density matrix renormalization group

Authors: Alberto Nocera, Gonzalo Alvarez

Abstract: We propose a method to compute spectral functions of generic Hamiltonians using the density matrix renormalization group (DMRG) algorithm directly in the frequency domain, based on a modified Krylov space decomposition to compute the correction-vectors. Our approach entails the calculation of the root-N (N=2 is the standard square root) of the Hamiltonian propagator using Krylov space decompositio… ▽ More We propose a method to compute spectral functions of generic Hamiltonians using the density matrix renormalization group (DMRG) algorithm directly in the frequency domain, based on a modified Krylov space decomposition to compute the correction-vectors. Our approach entails the calculation of the root-N (N=2 is the standard square root) of the Hamiltonian propagator using Krylov space decomposition, and repeating this procedure N times to obtain the actual correction-vector. We show that our method greatly alleviates the burden of keeping a large bond dimension at large target frequencies, a problem found with conventional correction-vector DMRG, while achieving better computational performance at large N. We apply our method to spin and charge spectral functions of t-J and Hubbard models in the challenging two-leg ladder geometry, and provide evidence that the root-N approach reaches a much improved resolution compared to conventional correction-vector. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 11 pages, 5 figures, Supplementary Material

Electronic structure, magnetic properties and pairing tendencies of the copper-based honeycomb lattice Na$_2$Cu$_2$TeO$_6$

Authors: Ling-Fang Lin, Rahul Soni, Yang Zhang, Shang Gao, Adriana Moreo, Gonzalo Alvarez, Andrew D. Christianson, Matthew B. Stone, Elbio Dagotto

Abstract: Spin-$1/2$ chains with alternating antiferromagnetic and ferromagnetic couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density matrix renormalization group methods, we have systematically studied the dimerized chain system Na$_2$Cu$_2$TeO$_6$. Near the Fermi level, the dominant states are mainly co… ▽ More Spin-$1/2$ chains with alternating antiferromagnetic and ferromagnetic couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density matrix renormalization group methods, we have systematically studied the dimerized chain system Na$_2$Cu$_2$TeO$_6$. Near the Fermi level, the dominant states are mainly contributed by the Cu $3d_{x^2-y^2}$ orbitals highly hybridized with the O $2p$ orbitals in the nonmagnetic phase, leading to an "effective" single-orbital low-energy model. Furthermore, the bandwidth of the Cu $3d_{x^2-y^2}$ states is small ($\sim 0.8$ eV), suggesting that electronic correlations will strongly affect this system. By introducing such electronic correlations, we found this system is a Mott insulator. Moreover, by calculating the magnetic exchange interactions ($J_1$, $J_2$ and $J_3$), we explained the size and sign of the exchange interactions in Na$_2$Cu$_2$TeO$_6$, in agreement with neutron experiments. In addition, we constructed a single-orbital Hubbard model for this dimerized chain system, where the quantum fluctuations are taken into account. Both AFM and FM coupling ($\uparrow$-$\downarrow$-$\downarrow$-$\uparrow$) along the chain were found in our DMRG and Lanczos calculations, in agreement with DFT and neutron results. We also calculated the hole pairing binding energy $ΔE$ which becomes negative at Hubbard $U \sim 11$ eV, indicating incipient pairing tendencies. Finally, we also looked at various cases of hole doping that always exhibit tight pairs. Thus, we believe our results for Na$_2$Cu$_2$TeO$_6$ could provide guidance to experimentalists and theorists working on this dimerized chain system, such as short-range magnetic coupling, doping effects, and possible pairing tendencies. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 10 pages, 12 figures

Magnetic excitations, non-classicality and quantum wake spin dynamics in the Hubbard chain

Authors: Pontus Laurell, Allen Scheie, D. Alan Tennant, Satoshi Okamoto, Gonzalo Alvarez, Elbio Dagotto

Abstract: Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations $G(r,t)$, a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor $S(k,ω)$. Here we apply theoretically these i… ▽ More Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations $G(r,t)$, a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor $S(k,ω)$. Here we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in $S(k,ω)$ going from the non-interacting limit ($U=0$) to strong coupling ($U\rightarrow \infty$), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group (DMRG) to find $S(k,ω)$, from which QFI is then calculated. We find that QFI grows with $U$. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above $U=2.5$ (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth $W(U)$ and the half-chain von Neumann entanglement entropy. We compute $G(r,t)$ by Fourier-transforming $S(k,ω)$. The results indicate a crossover in the short-time short-distance dynamics at low $U$ characterized by ferromagnetic lightcone wavefronts, to a Heisenberg-like behavior at large $U$ featuring antiferromagnetic lightcones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by $U=3$. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the $G(r,t)$ dynamics found at low $U$. △ Less

Submitted 5 August, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 15 pages, 10 figures

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

Path integral framework for characterizing and controlling decoherence induced by non-stationary environments on a quantum probe

Authors: Martin Kuffer, Analia Zwick, Gonzalo A. Alvarez

Abstract: Reliable processing of quantum information is a milestone to achieve for the deployment of quantum technologies. Uncontrolled, out-of-equilibrium sources of decoherence need to be characterized in detail for designing the control of quantum devices to mitigate the loss of quantum information. However, quantum sensing of such environments is still a challenge due to their non-stationary nature that… ▽ More Reliable processing of quantum information is a milestone to achieve for the deployment of quantum technologies. Uncontrolled, out-of-equilibrium sources of decoherence need to be characterized in detail for designing the control of quantum devices to mitigate the loss of quantum information. However, quantum sensing of such environments is still a challenge due to their non-stationary nature that in general can generate complex high-order correlations. We here introduce a path integral framework to characterize non-stationary environmental fluctuations by a quantum probe. We found the solution for the decoherence decay of non-stationary, generalized Gaussian processes that induce pure dephasing. This dephasing when expressed in a suitable basis, based on the non-stationary noise eigenmodes, is defined by the overlap of a generalized noise spectral density and a filter function that depends on the control fields. This result thus extends the validity to out-of-equilibrium environments, of the similar general expression for the dephasing of open quantum systems coupled to stationary noises. We show physical insights for a broad subclass of non-stationary noises that are local-in-time, in the sense that the noise correlation functions contain memory based on constraints of the derivatives of the fluctuating noise paths. Spectral and non-Markovian properties are discussed together with implementations of the framework to treat paradigmatic environments that are out-of-equilibrium, e.g. due to a quench and a pulsed noise. We show that our results provide tools for probing the spectral and time-correlation properties, and for mitigating decoherence effects of out-of-equilibrium -- non-stationary -- environments. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: 24 pages, 4 figures

Journal ref: PRX Quantum 3, 020321 (2022)

Strongly anisotropic electronic and magnetic structures in oxide dichlorides RuOCl$_2$ and OsOCl$_2$

Authors: Yang Zhang, Ling-Fang Lin, Adriana Moreo, Thomas A. Maier, Gonzalo Alvarez, Elbio Dagotto

Abstract: Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl$_2$ and OsOCl$_2$ with $d^4$ electronic configurations. Different from a previous study using VOI$_2$ with $d^1$ configuration, these systems with $4d^4$ or $5d^4$ do not exhibit a ferroelectric instability along the $a$-axis. Due to the fully-occup… ▽ More Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl$_2$ and OsOCl$_2$ with $d^4$ electronic configurations. Different from a previous study using VOI$_2$ with $d^1$ configuration, these systems with $4d^4$ or $5d^4$ do not exhibit a ferroelectric instability along the $a$-axis. Due to the fully-occupied $d_{xy}$ orbital in RuOCl$_2$ and OsOCl$_2$, the Peierls instability distortion disappears along the $b$-axis, leading to an undistorted I${\rm mmm}$ phase (No. 71). Furthermore, we observe strongly anisotropic electronic and magnetic structures along the $a$-axis. The large crystal-field splitting energy (between $d_{xz/yz}$ and $d_{xy}$ orbitals) and large hopping between nearest-neighbor Ru and Os atoms suppresses the spin-orbital effect in $M$OCl$_2$ ($M$ = Ru or Os) with electronic density $n = 4$, resulting in a spin-1 system instead of a $J = 0$ singlet ground state. Moreover, we find staggered antiferromagnetic order with $π$ wavevector along the $M$-O chain direction ($a$-axis) while the magnetic coupling along the $b$-axis is weak. Based on Wannier functions from first-principles calculations, we calculated the relevant hopping amplitudes and crystal-field splitting energies of the $t_{2g}$ orbitals for the Os atoms to construct a multi-orbital Hubbard model for the $M$-O chains. Staggered AFM with $\uparrow$-$\downarrow$-$\uparrow$-$\downarrow$ spin structure dominates in our DMRG calculations, in agreement with DFT calculations. △ Less

Submitted 14 April, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 14 pages, 15 figures

Journal ref: Phys. Rev. B 105, 174410 (2022)

Prediction of orbital selective Mott phases and block magnetic states in the quasi-one-dimensional iron chain Ce$_2$O$_2$FeSe$_2$ under hole and electron doping

Authors: Ling-Fang Lin, Yang Zhang, Gonzalo Alvarez, Jacek Herbrych, Adriana Moreo, Elbio Dagotto

Abstract: The recent detailed study of quasi-one-dimensional iron-based ladders, with the $3d$ iron electronic density $n = 6$, has unveiled surprises, such as orbital-selective phases. However, similar studies for $n=6$ iron chains are still rare. Here, a three-orbital electronic Hubbard model was constructed to study the magnetic and electronic properties of the quasi-one-dimensional $n=6$ iron chain Ce… ▽ More The recent detailed study of quasi-one-dimensional iron-based ladders, with the $3d$ iron electronic density $n = 6$, has unveiled surprises, such as orbital-selective phases. However, similar studies for $n=6$ iron chains are still rare. Here, a three-orbital electronic Hubbard model was constructed to study the magnetic and electronic properties of the quasi-one-dimensional $n=6$ iron chain Ce$_2$O$_2$FeSe$_2$, with focus on the effect of doping. Specifically, introducing the Hubbard $U$ and Hund $J_{H}$ couplings and studying the model via the density matrix renormalization group, we report the ground-state phase diagram varying the electronic density away from $n=6$. For the realistic Hund coupling $J_{H}/U = 1/4$, several electronic phases were obtained, including a metal, orbital-selective Mott, and Mott insulating phases. Doping away from the parent phase, the competition of many tendencies leads to a variety of magnetic states, such as ferromagnetism, as well as several antiferromagnetic and magnetic "block" phases. In the hole-doping region, two different interesting orbital-selective Mott phases were found: OSMP1 (with one localized orbital and two itinerant orbitals) and OSMP2 (with two localized orbitals and one itinerant orbital). Moreover, charge disproportionation phenomena were found in special doping regions. We argue that our predictions can be tested by simple modifications in the original chemical formula of Ce$_2$O$_2$FeSe$_2$. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 9 pages, 6 figures

Magnetization dynamics fingerprints of an excitonic condensate $t_{2g}^{4}$ magnet

Authors: Nitin Kaushal, Jacek Herbrych, Gonzalo Alvarez, Elbio Dagotto

Abstract: The competition between spin-orbit coupling $λ$ and electron-electron interaction $U$ leads to a plethora of novel states of matter, extensively studied in the context of $t_{2g}^4$ and $t_{2g}^5$ materials, such as ruthenates and iridates. Excitonic magnets -- the antiferromagnetic state of bounded electron-hole pairs -- is a prominent example of phenomena driven by those competing energy scales.… ▽ More The competition between spin-orbit coupling $λ$ and electron-electron interaction $U$ leads to a plethora of novel states of matter, extensively studied in the context of $t_{2g}^4$ and $t_{2g}^5$ materials, such as ruthenates and iridates. Excitonic magnets -- the antiferromagnetic state of bounded electron-hole pairs -- is a prominent example of phenomena driven by those competing energy scales. Interestingly, recent theoretical studies predicted that excitonic magnets can be found in the ground-state of spin-orbit-coupled $t_{2g}^4$ Hubbard models. Here, we present a detailed computational study of the magnetic excitations in that excitonic magnet, employing one-dimensional chains (via density matrix renormalization group) and small two-dimensional clusters (via Lanczos). Specifically, first we show that the low-energy spectrum is dominated by a dispersive (acoustic) magnonic mode, with extra features arising from the $λ=0$ state in the phase diagram. Second, and more importantly, we found a novel magnetic excitation forming a high-energy optical mode with the highest intensity at wavevector $q\to 0$. In the excitonic condensation regime at large $U$, we also have found a novel high-energy $π$-mode composed solely of orbital excitations. These unique fingerprints of the $t_{2g}^4$ excitonic magnet are important in the analysis of neutron and RIXS experiments. △ Less

Submitted 22 October, 2021; originally announced October 2021.

Excitonic wave-packet evolution in a two-orbital Hubbard model chain: A real-time real-space study

Authors: Bradraj Pandey, Gonzalo Alvarez, Elbio Dagotto

Abstract: Motivated by experimental developments introducing the concept of spin-orbit separation, we study the real-space time evolution of an excitonic wave-packet using a two-orbital Hubbard model. The exciton is created by exciting an electron from a lower energy half-filled orbital to a higher energy empty orbital. We carry out the real-time dynamics of the resulting excitonic wave-packet, using the ti… ▽ More Motivated by experimental developments introducing the concept of spin-orbit separation, we study the real-space time evolution of an excitonic wave-packet using a two-orbital Hubbard model. The exciton is created by exciting an electron from a lower energy half-filled orbital to a higher energy empty orbital. We carry out the real-time dynamics of the resulting excitonic wave-packet, using the time-dependent density matrix renormalization group method. We find clear evidence of charge-spin and spin-orbit separation in real-space, by tracking the time evolution of local observables. We show that the velocity of the orbiton can be tuned by varying the inter-orbital interactions. We also present a comparative study of a hole (in one orbital) and exciton (in two orbitals) dynamics in one-dimensional systems. Moreover, we analyze the dynamics of an exciton with spin-flip excitation, where we observe fractionalized spinons induced by Hund's interaction. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 7 pages, 5 figures

Journal ref: Phys. Rev. B 104, L220302 (2021)

Dynamics of quantum information scrambling under decoherence effects measured via active spins clusters

Authors: Federico D. Domínguez, Gonzalo A. Álvarez

Abstract: Developing quantum technologies requires the control and understanding of the non-equilibrium dynamics of quantum information in many-body systems. Local information propagates in the system by creating complex correlations known as information scrambling, as this process prevents extracting the information from local measurements. In this work, we develop a model adapted from solid-state NMR meth… ▽ More Developing quantum technologies requires the control and understanding of the non-equilibrium dynamics of quantum information in many-body systems. Local information propagates in the system by creating complex correlations known as information scrambling, as this process prevents extracting the information from local measurements. In this work, we develop a model adapted from solid-state NMR methods, to quantify the information scrambling. The scrambling is measured via time-reversal Loschmidt echoes (LE) and Multiple Quantum Coherences experiments that intrinsically contain imperfections. Considering these imperfections, we derive expressions for out-of-time-order correlators (OTOCs) to quantify the observable information scrambling based on measuring the number of active spins where the information was spread. Based on the OTOC expressions, decoherence effects arise naturally by the effects of the nonreverted terms in the LE experiment. Decoherence induces localization of the measurable degree of information scrambling. These effects define a localization cluster size for the observable number of active spins that determines a dynamical equilibrium. We contrast the model's predictions with quantum simulations performed with solid-state NMR experiments, that measure the information scrambling with time-reversal echoes with controlled imperfections. An excellent quantitative agreement is found with the dynamics of quantum information scrambling and its localization effects determined from the experimental data. The presented model and derived OTOCs set tools for quantifying the quantum information dynamics of large quantum systems (more than $10^{4}$ spins) consistent with experimental implementations that intrinsically contain imperfections. △ Less

Submitted 3 December, 2021; v1 submitted 8 July, 2021; originally announced July 2021.

Comments: 18 pages, 7 figures

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

Magnetic states of quasi-one-dimensional iron chalcogenide Ba$_2$FeS$_3$

Authors: Yang Zhang, Ling-Fang Lin, Gonzalo Alvarez, Adriana Moreo, Elbio Dagotto

Abstract: Quasi-one-dimensional iron-based ladders and chains, with the 3$d$ iron electronic density $n = 6$, are attracting considerable attention. Recently, a new iron chain system Ba$_2$FeS$_3$, also with $n = 6$, was prepared under high-pressure and high-temperature conditions. Here, the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on fi… ▽ More Quasi-one-dimensional iron-based ladders and chains, with the 3$d$ iron electronic density $n = 6$, are attracting considerable attention. Recently, a new iron chain system Ba$_2$FeS$_3$, also with $n = 6$, was prepared under high-pressure and high-temperature conditions. Here, the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered $\uparrow$-$\downarrow$-$\uparrow$-$\downarrow$ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard $U$ coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large $U/W$ ($W$ = bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba$_2$FeS$_3$ provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material. △ Less

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

Comments: 13 pages, 15 figures

Journal ref: Phys. Rev. B 104, 125122 (2021)

Origin of insulating ferromagnetism in iron oxychalcogenide Ce$_2$O$_2$FeSe$_2$

Authors: Ling-Fang Lin, Yang Zhang, Gonzalo Alvarez, Adriana Moreo, Elbio Dagotto

Abstract: An insulating ferromagnetic (FM) phase exists in the quasi-one-dimensional iron chalcogenide Ce$_2$O$_2$FeSe$_2$ but its origin is unknown. To understand the FM mechanism, here a systematic investigation of this material is provided, analyzing the competition between ferromagnetic and antiferromagnetic tendencies and the interplay of hoppings, Coulomb interactions, Hund's coupling, and crystal-fie… ▽ More An insulating ferromagnetic (FM) phase exists in the quasi-one-dimensional iron chalcogenide Ce$_2$O$_2$FeSe$_2$ but its origin is unknown. To understand the FM mechanism, here a systematic investigation of this material is provided, analyzing the competition between ferromagnetic and antiferromagnetic tendencies and the interplay of hoppings, Coulomb interactions, Hund's coupling, and crystal-field splittings. Our intuitive analysis based on second-order perturbation theory shows that large entanglements between doubly-occupied and half-filled orbitals play a key role in stabilizing the FM order in Ce$_2$O$_2$FeSe$_2$. In addition, via many-body computational techniques applied to a multi-orbital Hubbard model, the phase diagram confirms the proposed FM mechanism, in agreement with experiments. △ Less

Submitted 13 August, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: 6 pages, 5 figures

Journal ref: Phys. Rev. Lett. 127, 077204 (2021)

Intertwined charge, spin, and pairing orders in doped iron ladders

Authors: Bradraj Pandey, Rahul Soni, Ling-Fang Lin, Gonzalo Alvarez, Elbio Dagotto

Abstract: Motivated by recent experimental progress on iron-based ladder compounds, we study the doped two-orbital Hubbard model for the two-leg ladder BaFe$_2$S$_3$. The model is constructed by using {\it ab initio} hopping parameters and the ground state properties are investigated using the density matrix renormalization group method. We show that the $(π,0)$ magnetic ordering at half-filling, with ferro… ▽ More Motivated by recent experimental progress on iron-based ladder compounds, we study the doped two-orbital Hubbard model for the two-leg ladder BaFe$_2$S$_3$. The model is constructed by using {\it ab initio} hopping parameters and the ground state properties are investigated using the density matrix renormalization group method. We show that the $(π,0)$ magnetic ordering at half-filling, with ferromagnetic rungs and antiferromagnetic legs, becomes incommensurate upon hole doping. Moreover, depending on the strength of the Hubbard $U$ coupling, other magnetic patterns, such as $(0,π)$, are also stabilized. We found that the binding energy for two holes becomes negative for intermediate Hubbard interaction strength, indicating hole pairing. Due to the crystal-field split among orbitals, the holes primarily reside in one orbital, with the other one remaining half-filled. This resembles orbital selective Mott states. The formation of tight hole pairs continues with increasing hole density, as long as the magnetic order remains antiferromagnetic in one direction. The study of pair-pair correlations indicates the dominance of the intra-orbital spin-singlet channel, as opposed to other pairing channels. Although in a range of hole doping pairing correlations decay slowly, our results can also be interpreted as corresponding to a charge-density-wave made of pairs, a precursor of eventual superconductivity after interladder couplings are included. Such scenario of intertwined orders has been extensively discussed before in the cuprates, and our results suggest a similar physics could exist in ladder iron-based superconductors. Finally, we also show that a robust Hund's coupling is needed for pairing to occur. △ Less

Submitted 23 March, 2021; v1 submitted 10 March, 2021; originally announced March 2021.

Comments: 13 pages 10 figures

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

Witnessing entanglement in quantum magnets using neutron scattering

Authors: A. Scheie, Pontus Laurell, A. M. Samarakoon, B. Lake, S. E. Nagler, G. E. Granroth, S. Okamoto, G. Alvarez, D. A. Tennant

Abstract: We demonstrate how quantum entanglement can be directly witnessed in the quasi-1D Heisenberg antiferromagnet KCuF$_3$. We apply three entanglement witnesses --- one-tangle, two-tangle, and quantum Fisher information --- to its inelastic neutron spectrum, and compare with spectra simulated by finite-temperature density matrix renormalization group (DMRG) and classical Monte Carlo methods. We find t… ▽ More We demonstrate how quantum entanglement can be directly witnessed in the quasi-1D Heisenberg antiferromagnet KCuF$_3$. We apply three entanglement witnesses --- one-tangle, two-tangle, and quantum Fisher information --- to its inelastic neutron spectrum, and compare with spectra simulated by finite-temperature density matrix renormalization group (DMRG) and classical Monte Carlo methods. We find that each witness provides direct access to entanglement. Of these, quantum Fisher information is the most robust experimentally, and indicates the presence of at least bipartite entanglement up to at least 50 K, corresponding to around 10% of the spinon zone-boundary energy. We apply quantum Fisher information to higher spin-S Heisenberg chains, and show theoretically that the witnessable entanglement gets suppressed to lower temperatures as the quantum number increases. Finally, we outline how these results can be applied to higher dimensional quantum materials to witness and quantify entanglement. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 9 pages and 5 figures, four pages and six figures of appendices, and five pages supplemental information

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

Peierls transition, ferroelectricity, and spin-singlet formation in the monolayer VOI$_2$

Authors: Yang Zhang, Ling-Fang Lin, Adriana Moreo, Gonzalo Alvarez, Elbio Dagotto

Abstract: Using {\it ab initio} density functional theory and single-orbital Hubbard model calculations via the density matrix renormalization group method, we systematically studied the monolayer VOI$_2$ with a $3d^1$ electronic configuration. Our phonon calculations indicate that the orthorhombic $Pmm2$ FE-II phase is the most likely ground state, involving a ferroelectric distortion along the $a$-axis an… ▽ More Using {\it ab initio} density functional theory and single-orbital Hubbard model calculations via the density matrix renormalization group method, we systematically studied the monolayer VOI$_2$ with a $3d^1$ electronic configuration. Our phonon calculations indicate that the orthorhombic $Pmm2$ FE-II phase is the most likely ground state, involving a ferroelectric distortion along the $a$-axis and V-V dimerization along the $b$-axis. Specifically, the "pseudo Jahn-Teller" effect caused by the coupling between empty V ($d_{xz/yz}$ and $d_{3z^2-r^2}$) and O $2p$ states is proposed as the mechanism that stabilizes the ferroelectric distortion from the paraelectric phase. Moreover, the half-filled metallic $d_{xy}$ band displays a Peierls instability along the $b$-axis, inducing a V-V dimerization. We also found very short-range antiferromagnetic coupling along the V-V chain due to the formation of nearly-decoupled spin singlets in the ground state. △ Less

Submitted 22 February, 2021; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: 7 pages 5 figures

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

Gapless to gapless phase transitions in quantum spin chains

Authors: Shi Feng, Gonzalo Alvarez, Nandini Trivedi

Abstract: We investigate spin chains with bilinear-biquadratic (BLBQ) spin interactions as a function of an applied magnetic field $h$. At the Uimin-Lai-Sutherland (ULS) critical point we find a gapless to gapless transition revealed by the dynamical structure factor $S(q,ω)$ as a function of $h$. At $h=0$, the envelope of the lowest energy excitations goes soft at {\it two points} $q_1=2π/3$ and… ▽ More We investigate spin chains with bilinear-biquadratic (BLBQ) spin interactions as a function of an applied magnetic field $h$. At the Uimin-Lai-Sutherland (ULS) critical point we find a gapless to gapless transition revealed by the dynamical structure factor $S(q,ω)$ as a function of $h$. At $h=0$, the envelope of the lowest energy excitations goes soft at {\it two points} $q_1=2π/3$ and $q_2=4π/3$, dubbed the phase A. With increasing field, the spectral peaks at each of the gapless points bifurcate, making in total {\it four} soft modes, and combine to form a new set of excitations that soften at a {\it single} point $q=π$ at $h_{c1}\approx 0.94$. Beyond $h_{c1}$ the system enters another gapless B-phase until the transition at $h_{c2}=4$ to the fully polarized phase. We compare the ULS model results with those for the AKLT model as a representative of gapped Haldane phase. We explain the mechanism of the gapless to gapless transition in the ULS model using its conserved charges and a spinon band picture. We also discuss the universality of central charges of the BLBQ family of models subjected to a magnetic field. △ Less

Submitted 20 January, 2022; v1 submitted 1 December, 2020; originally announced December 2020.

Comments: 26 pages, 16 figures

Journal ref: Phys. Rev. B 105, 014435 (2022)

Interaction-induced topological phase transition and Majorana edge states in low-dimensional orbital-selective Mott insulators

Authors: J. Herbrych, M. Środa, G. Alvarez, M. Mierzejewski, E. Dagotto

Abstract: Topological phases of matter are among the most intriguing research directions in Condensed Matter Physics. It is known that superconductivity induced on a topological insulator's surface can lead to exotic Majorana modes, the main ingredient of many proposed quantum computation schemes. In this context, the iron-based high critical temperature superconductors are a promising platform to host such… ▽ More Topological phases of matter are among the most intriguing research directions in Condensed Matter Physics. It is known that superconductivity induced on a topological insulator's surface can lead to exotic Majorana modes, the main ingredient of many proposed quantum computation schemes. In this context, the iron-based high critical temperature superconductors are a promising platform to host such an exotic phenomenon in real condensed-matter compounds. The Coulomb interaction is commonly believed to be vital for the magnetic and superconducting properties of these systems. This work bridges these two perspectives and shows that the Coulomb interaction can also drive a canonical superconductor with orbital degrees of freedom into the topological state. Namely, we show that above a critical value of the Hubbard interaction the system simultaneously develops spiral spin order, a highly unusual triplet amplitude in superconductivity, and, remarkably, Majorana fermions at the edges of the system. △ Less

Submitted 19 May, 2021; v1 submitted 11 November, 2020; originally announced November 2020.

Journal ref: Nat. Commun. 12, 2955 (2021)

Simulating Spin Waves in Entropy Stabilized Oxides

Authors: Tom Berlijn, Gonzalo Alvarez, David S. Parker, Raphaël P. Hermann, Randy S. Fishman

Abstract: The entropy stabilized oxide Mg$_{0.2}$Co$_{0.2}$Ni$_{0.2}$Cu$_{0.2}$Zn$_{0.2}$O exhibits antiferromagnetic order and magnetic excitations, as revealed by recent neutron scattering experiments. This observation raises the question of the nature of spin wave excitations in such disordered systems. Here, we investigate theoretically the magnetic ground state and the spin-wave excitations using linea… ▽ More The entropy stabilized oxide Mg$_{0.2}$Co$_{0.2}$Ni$_{0.2}$Cu$_{0.2}$Zn$_{0.2}$O exhibits antiferromagnetic order and magnetic excitations, as revealed by recent neutron scattering experiments. This observation raises the question of the nature of spin wave excitations in such disordered systems. Here, we investigate theoretically the magnetic ground state and the spin-wave excitations using linear spin-wave theory in combination with the supercell approximation to take into account the extreme disorder in this magnetic system. We find that the experimentally observed antiferromagnetic structure can be stabilized by a rhombohedral distortion together with large second nearest neighbor interactions. Our calculations show that the spin-wave spectrum consists of a well-defined low-energy coherent spectrum in the background of an incoherent continuum that extends to higher energies. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Journal ref: Phys. Rev. Research 3, 033273 (2021)

Quantifying and controlling entanglement in the quantum magnet Cs$_2$CoCl$_4$

Authors: Pontus Laurell, Allen Scheie, Chiron J. Mukherjee, Michael M. Koza, Mechtild Enderle, Zbigniew Tylczynski, Satoshi Okamoto, Radu Coldea, D. Alan Tennant, Gonzalo Alvarez

Abstract: The lack of methods to experimentally detect and quantify entanglement in quantum matter impedes our ability to identify materials hosting highly entangled phases, such as quantum spin liquids. We thus investigate the feasibility of using inelastic neutron scattering (INS) to implement a model-independent measurement protocol for entanglement based on three entanglement witnesses: one-tangle, two-… ▽ More The lack of methods to experimentally detect and quantify entanglement in quantum matter impedes our ability to identify materials hosting highly entangled phases, such as quantum spin liquids. We thus investigate the feasibility of using inelastic neutron scattering (INS) to implement a model-independent measurement protocol for entanglement based on three entanglement witnesses: one-tangle, two-tangle, and quantum Fisher information (QFI). We perform high-resolution INS measurements on Cs$_2$CoCl$_4$, a close realization of the $S=1/2$ transverse-field XXZ spin chain, where we can control entanglement using the magnetic field, and compare with density-matrix renormalization group calculations for validation. The three witnesses allow us to infer entanglement properties and make deductions about the quantum state in the material. We find QFI to be a particularly robust experimental probe of entanglement, whereas the one- and two-tangles require more careful analysis. Our results lay the foundation for a general entanglement detection protocol for quantum spin systems. △ Less

Submitted 20 May, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: Main text: 7 pages, 4 figures. Supplementary Information: 15 pages, 15 figures

Journal ref: Phys. Rev. Lett. 127, 037201 (2021)

Origin of the Magnetic and Orbital ordering in $α$-Sr$_2$CrO$_4$

Authors: Bradraj Pandey, Yang Zhang, Nitin Kaushal, Rahul Soni, Ling-Fang Lin, Wen-Jun Hu, Gonzalo Alvarez, Elbio Dagotto

Abstract: Motivated by recent experimental progress in transition metal oxides with the K$_2$NiF$_4$ structure, we investigate the magnetic and orbital ordering in $α$-Sr$_2$CrO$_4$. Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the {\it ab initio} band structure near the Fermi level. The unique reverse splitting of $t_{2g}$ orbitals in $α$-Sr$_2$CrO… ▽ More Motivated by recent experimental progress in transition metal oxides with the K$_2$NiF$_4$ structure, we investigate the magnetic and orbital ordering in $α$-Sr$_2$CrO$_4$. Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the {\it ab initio} band structure near the Fermi level. The unique reverse splitting of $t_{2g}$ orbitals in $α$-Sr$_2$CrO$_4$, with the $3d^2$ electronic configuration for the Cr$^{4+}$ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground-state phase diagram for the two-dimensional compound $α$-Sr$_2$CrO$_4$. We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of $α$-Sr$_2$CrO$_4$, we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound $α$-Sr$_2$CrO$_4$, in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with $d_{xy}$ driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the $d_{xz}$ and $d_{yz}$ orbitals driving the orbital order. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: 13 pages, 12 figures

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

Block orbital-selective Mott insulators: a spin excitation analysis

Authors: J. Herbrych, G. Alvarez, A. Moreo, E. Dagotto

Abstract: We present a comprehensive study of the spin excitations - as measured by the dynamical spin structure factor $S(q,ω)$ - of the so-called block-magnetic state of low-dimensional orbital-selective Mott insulators. We realize this state via both a multi-orbital Hubbard model and a generalized Kondo-Heisenberg Hamiltonian. Due to various competing energy scales present in the models, the system devel… ▽ More We present a comprehensive study of the spin excitations - as measured by the dynamical spin structure factor $S(q,ω)$ - of the so-called block-magnetic state of low-dimensional orbital-selective Mott insulators. We realize this state via both a multi-orbital Hubbard model and a generalized Kondo-Heisenberg Hamiltonian. Due to various competing energy scales present in the models, the system develops periodic ferromagnetic islands of various shapes and sizes, which are antiferromagnetically coupled. The 2$\times$2 particular case was already found experimentally in the ladder material BaFe$_2$Se$_3$ that becomes superconducting under pressure. Here we discuss the electronic density as well as Hubbard and Hund coupling dependence of $S(q,ω)$ using density matrix renormalization group method. Several interesting features were identified: (1) An acoustic (dispersive spin-wave) mode develops. (2) The spin-wave bandwidth establishes a new energy scale that is strongly dependent on the size of the magnetic island and becomes abnormally small for large clusters. (3) Optical (dispersionless spin excitation) modes are present for all block states studied here. In addition, a variety of phenomenological spin Hamiltonians have been investigated but none matches entirely our results that were obtained primarily at intermediate Hubbard $U$ strengths. Our comprehensive analysis provides theoretical guidance and motivation to crystal growers to search for appropriate candidate materials to realize the block states, and to neutron scattering experimentalists to confirm the exotic dynamical magnetic properties unveiled here, with a rich mixture of acoustic and optical features. △ Less

Submitted 16 September, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. B 102, 115134 (2020)

Prediction of exotic magnetic states in the alkali metal quasi-one-dimensional iron selenide compound Na$_2$FeSe$_2$

Authors: Bradraj Pandey, Ling-Fang Lin, Rahul Soni, Nitin Kaushal, Jacek Herbrych, Gonzalo Alvarez, Elbio Dagotto

Abstract: The magnetic and electronic phase diagram of a model for the quasi-one-dimensional alkali metal iron selenide compound Na$_2$FeSe$_2$ is presented. The novelty of this material is that the valence of iron is Fe$^{2+}$ contrary to most other iron-chain compounds with valence Fe$^{3+}$. Using first-principles techniques, we developed a three-orbital tight-binding model that reproduces the {\it ab in… ▽ More The magnetic and electronic phase diagram of a model for the quasi-one-dimensional alkali metal iron selenide compound Na$_2$FeSe$_2$ is presented. The novelty of this material is that the valence of iron is Fe$^{2+}$ contrary to most other iron-chain compounds with valence Fe$^{3+}$. Using first-principles techniques, we developed a three-orbital tight-binding model that reproduces the {\it ab initio} band structure near the Fermi level. Including Hubbard and Hund couplings and studying the model via the density matrix renormalization group and Lanczos methods, we constructed the ground state phase diagram. A robust region where the block state $\uparrow \uparrow \downarrow \downarrow \uparrow \uparrow \downarrow \downarrow$ is stabilized was unveiled. The analog state in iron ladders, employing 2$\times$2 ferromagnetic blocks, is by now well-established, but in chains a block magnetic order has not been observed yet in real materials. The phase diagram also contains a large region of canonical staggered spin order $\uparrow \downarrow \uparrow \downarrow \uparrow \downarrow \uparrow$ at very large Hubbard repulsion. At the block to staggered transition region, a novel phase is stabilized with a mixture of both states: an inhomogeneous orbital-selective charge density wave with the exotic spin configuration $\uparrow \uparrow \downarrow \uparrow \downarrow \downarrow \uparrow \downarrow$. Our predictions for Na$_2$FeSe$_2$ may guide crystal growers and neutron scattering experimentalists towards the realization of block states in one-dimensional iron-selenide chain materials. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: 10 pages, 5 figures

Journal ref: Phys. Rev. B 102, 035149 (2020)

Decoherence scaling transition in the dynamics of quantum information scrambling

Authors: Federico D. Dominguez, Maria Cristina Rodriguez, Robin Kaiser, Dieter Suter, Gonzalo A. Alvarez

Abstract: Reliable processing of quantum information for developing quantum technologies requires precise control of out-of-equilibrium many-bodysystems. This is a highly challenging task as the fragility of quantum states to external perturbations increases with the system-size. Here, we report on a series of experimental quantum simulations that quantify the sensitivity of a controlled Hamiltonian evoluti… ▽ More Reliable processing of quantum information for developing quantum technologies requires precise control of out-of-equilibrium many-bodysystems. This is a highly challenging task as the fragility of quantum states to external perturbations increases with the system-size. Here, we report on a series of experimental quantum simulations that quantify the sensitivity of a controlled Hamiltonian evolution to perturbations that drive the system away from the targeted evolution. Based on out-of-time ordered correlations, we demonstrate that the decay-rate of the process fidelity increases with the effective number $K$ of correlated qubits as $K^α$. As a function of the perturbation strength, we observe a decoherence scaling transition of the exponent $α$ between two distinct dynamical regimes. In the limiting case below the critical perturbation strength, the exponent $α$ drops sharply below 1, and there is no inherent limit to the number of qubits that can be controlled. This resilient quantum feature of the controlled dynamics of quantum information is promising for reliable control of large quantum systems. △ Less

Submitted 18 June, 2021; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: 5 pages, 3 figures, material and methods

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

BCS-BEC crossover in a $(t_{2g})^4$ Excitonic Magnet

Authors: Nitin Kaushal, Rahul Soni, Alberto Nocera, Gonzalo Alvarez, Elbio Dagotto

Abstract: The condensation of spin-orbit-induced excitons in $(t_{2g})^4$ electronic systems is attracting considerable attention. In the large Hubbard U limit, antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ($J_{\textrm{eff}} = 1$). In this publication, we show that even for the weak and intermediate U regimes, the spin-orbit exciton condensation is possible… ▽ More The condensation of spin-orbit-induced excitons in $(t_{2g})^4$ electronic systems is attracting considerable attention. In the large Hubbard U limit, antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ($J_{\textrm{eff}} = 1$). In this publication, we show that even for the weak and intermediate U regimes, the spin-orbit exciton condensation is possible leading also to staggered magnetic order. The canonical electron-hole excitations (excitons) transform into local triplon excitations at large U , and this BEC strong coupling regime is smoothly connected to the intermediate U excitonic insulator region. We solved the degenerate three-orbital Hubbard model with spin-orbit coupling ($λ$) in one-dimensional geometry using the Density Matrix Renormalization Group, while in two-dimensional square clusters we use the Hartree-Fock approximation (HFA). Employing these techniques, we provide the full $λ$ vs U phase diagrams for both one- and two- dimensional lattices. Our main result is that at the intermediate Hubbard U region of our focus, increasing $λ$ at fixed U the system transitions from an incommensurate spin-density-wave metal to a Bardeen-Cooper-Schrieffer (BCS) excitonic insulator, with coherence length r coh of O(a) and O(10a) in 1d and 2d, respectively, with a the lattice spacing. Further increasing $λ$, the system eventually crosses over to the BEC limit (with r coh << a). △ Less

Submitted 17 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. B 101, 245147 (2020)

Emergence of Superconductivity in Doped Multiorbital Hubbard Chains

Authors: Niravkumar D. Patel, Nitin Kaushal, Alberto Nocera, Gonzalo Alvarez, Elbio Dagotto

Abstract: We introduce a variational state for one-dimensional two-orbital Hubbard models that intuitively explains the recent computational discovery of pairing in these systems when hole doped. Our Ansatz is an optimized linear superposition of Affleck-Kennedy-Lieb-Tasaki valence bond states, rendering the combination a valence bond liquid dubbed Orbital Resonant Valence Bond. We show that the undoped (on… ▽ More We introduce a variational state for one-dimensional two-orbital Hubbard models that intuitively explains the recent computational discovery of pairing in these systems when hole doped. Our Ansatz is an optimized linear superposition of Affleck-Kennedy-Lieb-Tasaki valence bond states, rendering the combination a valence bond liquid dubbed Orbital Resonant Valence Bond. We show that the undoped (one electron/orbital) quantum state of two sites coupled into a global spin singlet is exactly written employing only spin-1/2 singlets linking orbitals at nearest-neighbor sites. Generalizing to longer chains defines our variational state visualized geometrically expressing our chain as a two-leg ladder, with one orbital per leg. As in Anderson's resonating valence-bond state, our undoped variational state contains preformed singlet pairs that via doping become mobile leading to superconductivity. Doped real materials with one-dimensional substructures, two near-degenerate orbitals, and intermediate Hubbard U/W strengths -- W the carrier's bandwidth -- could realize spin-singlet pairing if on-site anisotropies are small. If these anisotropies are robust, spin-triplet pairing emerges. △ Less

Submitted 25 May, 2020; v1 submitted 2 January, 2020; originally announced January 2020.

Journal ref: NPJ Quantum Mater. 5, 27 (2020)

Block-spiral magnetism: An exotic type of frustrated order

Authors: J. Herbrych, J. Heverhagen, G. Alvarez, M. Daghofer, A. Moreo, E. Dagotto

Abstract: Competing interactions in Quantum Materials induce novel states of matter such as frustrated magnets, an extensive field of research both from the theoretical and experimental perspectives. Here, we show that competing energy scales present in the low-dimensional orbital-selective Mott phase (OSMP) induce an exotic magnetic order, never reported before. Earlier neutron scattering experiments on ir… ▽ More Competing interactions in Quantum Materials induce novel states of matter such as frustrated magnets, an extensive field of research both from the theoretical and experimental perspectives. Here, we show that competing energy scales present in the low-dimensional orbital-selective Mott phase (OSMP) induce an exotic magnetic order, never reported before. Earlier neutron scattering experiments on iron-based 123 ladder materials, where OSMP is relevant, already confirmed our previous theoretical prediction of block-magnetism (magnetic order of the form $\uparrow\uparrow\downarrow\downarrow$). Now we argue that another novel phase can be stabilized in multi-orbital Hubbard models, the {\it block-spiral state}. In this state, the magnetic islands form a spiral propagating through the chain but with the blocks maintaining their identity, namely rigidly rotating. This new spiral state is stabilized without any apparent frustration, the common avenue to generate spiral arrangements in multiferroics. By examining the behaviour of the electronic degrees of freedom, parity breaking quasiparticles are revealed. Finally, a simple phenomenological model that accurately captures the macroscopic spin spiral arrangement is also introduced, and fingerprints for the neutron scattering experimental detection of our new state are provided. △ Less

Submitted 15 July, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Journal ref: Proc. Natl Acad. Sci. USA 117, 16226 (2020)

Two-particle Correlation Functions in Cluster Perturbation Theory: Hubbard Spin Susceptibilities

Authors: P. T. Raum, G. Alvarez, Thomas Maier, V. W. Scarola

Abstract: Cluster Perturbation Theory (CPT) is a computationally economic method commonly used to estimate the momentum and energy resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-… ▽ More Cluster Perturbation Theory (CPT) is a computationally economic method commonly used to estimate the momentum and energy resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-particle correlation functions. CPT can be extended to compute two-particle correlation functions by approximately solving the Bethe-Salpeter equation. We implement this method and focus on the transverse spin-susceptibility, measurable via inelastic neutron scattering or with optical probes of atomic gases in optical lattices. We benchmark the method with the one-dimensional Fermi-Hubbard model at half filling by comparing with known results. △ Less

Submitted 18 February, 2020; v1 submitted 27 August, 2019; originally announced August 2019.

Journal ref: Phys. Rev. B 101, 075122 (2020)

Targeting Multiple States in the Density Matrix Renormalization Group with The Singular Value Decomposition

Authors: E. F. D'Azevedo, W. R. Elwasif, N. D. Patel, G. Alvarez

Abstract: In the Density Matrix Renormalization Group (DMRG), multiple states must be included in the density matrix when properties beyond ground state are needed, including temperature dependence, time evolution, and frequency-resolved response functions. How to include these states in the density matrix has been shown in the past. But it is advantageous to replace the density matrix by a singular value d… ▽ More In the Density Matrix Renormalization Group (DMRG), multiple states must be included in the density matrix when properties beyond ground state are needed, including temperature dependence, time evolution, and frequency-resolved response functions. How to include these states in the density matrix has been shown in the past. But it is advantageous to replace the density matrix by a singular value decomposition (SVD) instead, because of improved performance, and because it enables multiple targeting in the matrix product state description of the DMRG. This paper shows how to target multiple states using the SVD; it analyzes the implication of local symmetries, and discusses typical performance improvements using the example of the Hubbard model's photo-emission spectra on a ladder geometry. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 7 pages, 4 figures, 1 table

Novel Block Excitonic Condensate at $n=3.5$ in a Spin-Orbit Coupled $t_{2g}$ Multiorbital Hubbard Model

Authors: Nitin Kaushal, Alberto Nocera, Gonzalo Alvarez, Adriana Moreo, Elbio Dagotto

Abstract: Theoretical studies recently predicted the condensation of spin-orbit excitons at momentum $q$=$π$ in $t_{2g}^4$ spin-orbit coupled three-orbital Hubbard models at electronic density $n=4$. In parallel, experiments involving iridates with non-integer valence states for the Ir ions are starting to attract considerable attention. In this publication, using the density matrix renormalization group te… ▽ More Theoretical studies recently predicted the condensation of spin-orbit excitons at momentum $q$=$π$ in $t_{2g}^4$ spin-orbit coupled three-orbital Hubbard models at electronic density $n=4$. In parallel, experiments involving iridates with non-integer valence states for the Ir ions are starting to attract considerable attention. In this publication, using the density matrix renormalization group technique we present evidence for the existence of a novel excitonic condensate at $n=3.5$ in a one-dimensional Hubbard model with a degenerate $t_{2g}$ sector, when in the presence of spin-orbit coupling. At intermediate Hubbard $U$ and spin-orbit $λ$ couplings, we found an excitonic condensate at the unexpected momentum $q$=$π/2$ involving $j_{\textrm{eff}}=3/2,m=\pm1/2$ and $j_{\textrm{eff}}=1/2,m=\pm1/2$ bands in the triplet channel, coexisting with an also unexpected block magnetic order. We also present the entire $λ$ vs $U$ phase diagram, at a fixed and robust Hund coupling. Interestingly, this new `block excitonic phase' is present even at large values of $λ$, unlike the $n=4$ excitonic phase discussed before. Our computational study helps to understand and predict the possible magnetic phases of materials with $d^{3.5}$ valence and robust spin-orbit coupling. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Journal ref: Phys. Rev. B 99, 155115 (2019)

DCA++: A software framework to solve correlated electron problems with modern quantum cluster methods

Authors: Urs R. Hähner, Gonzalo Alvarez, Thomas A. Maier, Raffaele Solcà, Peter Staar, Michael S. Summers, Thomas C. Schulthess

Abstract: We present the first open release of the DCA++ project, a high-performance research software framework to solve quantum many-body problems with cutting edge quantum cluster algorithms. DCA++ implements the dynamical cluster approximation (DCA) and its DCA$^+$ extension with a continuous self-energy. The algorithms capture nonlocal correlations in strongly correlated electron systems, thereby givin… ▽ More We present the first open release of the DCA++ project, a high-performance research software framework to solve quantum many-body problems with cutting edge quantum cluster algorithms. DCA++ implements the dynamical cluster approximation (DCA) and its DCA$^+$ extension with a continuous self-energy. The algorithms capture nonlocal correlations in strongly correlated electron systems, thereby giving insight into high-$T_c$ superconductivity. The code's scalability allows efficient usage of systems at all scales, from workstations to leadership computers. With regard to the increasing heterogeneity of modern computing machines, DCA++ provides portable performance on conventional and emerging new architectures, such as hybrid CPU-GPU, sustaining multiple petaflops on ORNL's Titan and CSCS' Piz Daint supercomputers. Moreover, we show how sustainable and scalable development of the code base has been achieved by adopting standard techniques of the software industry. These include employing a distributed version control system, applying test-driven development and following continuous integration. △ Less

Submitted 2 September, 2019; v1 submitted 4 January, 2019; originally announced January 2019.

Comments: 32 pages, 17 figures, 1 table, 3 code listings; submitted to Comput. Phys. Commun

Journal ref: Comput. Phys. Commun. 246 (2020) 106709

Novel Magnetic Block States in Low-Dimensional Iron-Based Superconductors

Authors: J. Herbrych, J. Heverhagen, N. D. Patel, G. Alvarez, M. Daghofer, A. Moreo, E. Dagotto

Abstract: Inelastic neutron scattering recently confirmed the theoretical prediction of a $\uparrow\uparrow\downarrow\downarrow$-magnetic state along the legs of quasi-one-dimensional (quasi-1D) iron-based ladders in the orbital-selective Mott phase (OSMP). We show here that electron-doping of the OSMP induces a whole class of novel block-states with a variety of periodicities beyond the previously reported… ▽ More Inelastic neutron scattering recently confirmed the theoretical prediction of a $\uparrow\uparrow\downarrow\downarrow$-magnetic state along the legs of quasi-one-dimensional (quasi-1D) iron-based ladders in the orbital-selective Mott phase (OSMP). We show here that electron-doping of the OSMP induces a whole class of novel block-states with a variety of periodicities beyond the previously reported $π/2$ pattern. We discuss the magnetic phase diagram of the OSMP regime that could be tested by neutrons once appropriate quasi-1D quantum materials with the appropriate dopings are identified. △ Less

Submitted 10 July, 2019; v1 submitted 1 December, 2018; originally announced December 2018.

Journal ref: Phys. Rev. Lett. 123, 027203 (2019)