Showing 1–50 of 68 results for author: Barros, K

Classical dynamics of the antiferromagnetic Heisenberg $S=1/2$ spin ladder

Authors: David A. Dahlbom, Jinu Thomas, Steven Johnston, Kipton Barros, Cristian D. Batista

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

Submitted 25 May, 2024; originally announced May 2024.

Comments: 15 pages, 8 figures

Kernel Polynomial Method for Linear Spin Wave Theory

Authors: Harry Lane, Hao Zhang, David Dahlbom, Sam Quinn, Rolando D. Somma, Martin Mourigal, Cristian D. Batista, Kipton Barros

Abstract: Calculating dynamical spin correlations is essential for matching model magnetic exchange Hamiltonians to momentum-resolved spectroscopic measurements. A major numerical bottleneck is the diagonalization of the dynamical matrix, especially in systems with large magnetic unit cells, such as those with incommensurate magnetic structures or quenched disorder. In this paper, we demonstrate an efficien… ▽ More Calculating dynamical spin correlations is essential for matching model magnetic exchange Hamiltonians to momentum-resolved spectroscopic measurements. A major numerical bottleneck is the diagonalization of the dynamical matrix, especially in systems with large magnetic unit cells, such as those with incommensurate magnetic structures or quenched disorder. In this paper, we demonstrate an efficient scheme based on the kernel polynomial method for calculating dynamical correlations of relevance to inelastic neutron scattering experiments. This method reduces the scaling of numerical cost from cubic to linear in the magnetic unit cell size. △ Less

Submitted 5 February, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: Submission to SciPost. 27 pages, 6 figures. Typos corrected, references added

Active Learning approach to simulations of Strongly Correlated Matter with the Ghost Gutzwiller Approximation

Authors: Marius S. Frank, Denis G. Artiukhin, Tsung-Han Lee, Yongxin Yao, Kipton Barros, Ove Christiansen, Nicola Lanatà

Abstract: Quantum embedding (QE) methods such as the Ghost Gutzwiller Approximation (gGA) offer a powerful approach to simulating strongly-correlated systems, but come with the computational bottleneck of computing the ground state of an auxiliary embedding Hamiltonian (EH) iteratively. In this work, we introduce an active learning (AL) framework integrated within the gGA to address this challenge. The meth… ▽ More Quantum embedding (QE) methods such as the Ghost Gutzwiller Approximation (gGA) offer a powerful approach to simulating strongly-correlated systems, but come with the computational bottleneck of computing the ground state of an auxiliary embedding Hamiltonian (EH) iteratively. In this work, we introduce an active learning (AL) framework integrated within the gGA to address this challenge. The methodology is applied to the single-band Hubbard model and results in a significant reduction in the number of instances where the EH must be solved. Through a principal component analysis (PCA), we find that the EH parameters form a low-dimensional structure that is largely independent of the geometric specifics of the systems, especially in the strongly-correlated regime. Our AL strategy enables us to discover this low-dimensionality structure on the fly, while leveraging it for reducing the computational cost of gGA, laying the groundwork for more efficient simulations of complex strongly-correlated materials. △ Less

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

Comments: 14 pages, 13 figures

Journal ref: Phys. Rev. Research 6, 013242 (2024)

SmoQyDQMC.jl: A flexible implementation of determinant quantum Monte Carlo for Hubbard and electron-phonon interactions

Authors: Benjamin Cohen-Stead, Sohan Malkaruge Costa, James Neuhaus, Andy Tanjaroon Ly, Yutan Zhang, Richard Scalettar, Kipton Barros, Steven Johnston

Abstract: We introduce the SmoQyDQMC.jl package, a Julia implementation of the determinant quantum Monte Carlo algorithm. SmoQyDQMC.jl supports generalized tight-binding Hamiltonians with on-site Hubbard and generalized electron-phonon interactions, including non-linear $e$-ph coupling and anharmonic lattice potentials. Our implementations use hybrid Monte Carlo methods with exact forces for sampling the ph… ▽ More We introduce the SmoQyDQMC.jl package, a Julia implementation of the determinant quantum Monte Carlo algorithm. SmoQyDQMC.jl supports generalized tight-binding Hamiltonians with on-site Hubbard and generalized electron-phonon interactions, including non-linear $e$-ph coupling and anharmonic lattice potentials. Our implementations use hybrid Monte Carlo methods with exact forces for sampling the phonon fields, enabling efficient simulation of low-energy phonon branches, including acoustic phonons. The SmoQyDQMC.jl package also uses a flexible scripting interface, allowing users to adapt it to different workflows and interface with other software packages in the Julia ecosystem. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDQMC.jl or installed using the Julia package manager. The online documentation, including examples, can be obtained from our document page at https://smoqysuite.github.io/SmoQyDQMC.jl/stable/. △ Less

Submitted 17 April, 2024; v1 submitted 15 November, 2023; originally announced November 2023.

Quantum to classical crossover in generalized spin systems -- the temperature-dependent spin dynamics of FeI$_2$

Authors: D. Dahlbom, D. Brooks, M. S. Wilson, S. Chi, A. I. Kolesnikov, M. B. Stone, H. Cao, Y. -W. Li, K. Barros, M. Mourigal, C. D. Batista, X. Bai

Abstract: Simulating quantum spin systems at finite temperatures is an open challenge in many-body physics. This work studies the temperature-dependent spin dynamics of a pivotal compound, FeI$_2$, to determine if universal quantum effects can be accounted for by a phenomenological renormalization of the dynamical spin structure factor $S(\mathbf{q}, ωおめが)$ measured by inelastic neutron scattering. Renormaliza… ▽ More Simulating quantum spin systems at finite temperatures is an open challenge in many-body physics. This work studies the temperature-dependent spin dynamics of a pivotal compound, FeI$_2$, to determine if universal quantum effects can be accounted for by a phenomenological renormalization of the dynamical spin structure factor $S(\mathbf{q}, ωおめが)$ measured by inelastic neutron scattering. Renormalization schemes based on the quantum-to-classical correspondence principle are commonly applied at low temperatures to the harmonic oscillators describing normal modes. However, it is not clear how to extend this renormalization to arbitrarily high temperatures. Here we introduce a temperature-dependent normalization of the classical moments, whose magnitude is determined by imposing the quantum sum rule, i.e. $\int dωおめがd\mathbf{q} S(\mathbf{q}, ωおめが) = N_S S (S+1)$ for $N_S$ dipolar magnetic moments. We show that this simple renormalization scheme significantly improves the agreement between the calculated and measured $S(\mathbf{q}, ωおめが)$ for FeI$_{2}$ at all temperatures. Due to the coupled dynamics of dipolar and quadrupolar moments in that material, this renormalization procedure is extended to classical theories based on SU(3) coherent states, and by extension, to any SU(N) coherent state representation of local multipolar moments. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: Associated source code for reproducing calculations available at: https://github.com/SunnySuite/SunnyContributed

Machine learning potentials with Iterative Boltzmann Inversion: training to experiment

Authors: Sakib Matin, Alice Allen, Justin S. Smith, Nicholas Lubbers, Ryan B. Jadrich, Richard A. Messerly, Benjamin T. Nebgen, Ying Wai Li, Sergei Tretiak, Kipton Barros

Abstract: Methodologies for training machine learning potentials (MLPs) to quantum-mechanical simulation data have recently seen tremendous progress. Experimental data has a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. We investigate a training procedure based on Iterative Boltzmann Invers… ▽ More Methodologies for training machine learning potentials (MLPs) to quantum-mechanical simulation data have recently seen tremendous progress. Experimental data has a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. We investigate a training procedure based on Iterative Boltzmann Inversion that produces a pair potential correction to an existing MLP, using equilibrium radial distribution function data. By applying these corrections to a MLP for pure aluminum based on Density Functional Theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require auto-differentiating through a molecular dynamics solver, and does not make assumptions about the MLP architecture. The results suggest a practical framework of incorporating experimental data into machine learning models to improve accuracy of molecular dynamics simulations. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Learning Together: Towards foundational models for machine learning interatomic potentials with meta-learning

Authors: Alice E. A. Allen, Nicholas Lubbers, Sakib Matin, Justin Smith, Richard Messerly, Sergei Tretiak, Kipton Barros

Abstract: The development of machine learning models has led to an abundance of datasets containing quantum mechanical (QM) calculations for molecular and material systems. However, traditional training methods for machine learning models are unable to leverage the plethora of data available as they require that each dataset be generated using the same QM method. Taking machine learning interatomic potentia… ▽ More The development of machine learning models has led to an abundance of datasets containing quantum mechanical (QM) calculations for molecular and material systems. However, traditional training methods for machine learning models are unable to leverage the plethora of data available as they require that each dataset be generated using the same QM method. Taking machine learning interatomic potentials (MLIPs) as an example, we show that meta-learning techniques, a recent advancement from the machine learning community, can be used to fit multiple levels of QM theory in the same training process. Meta-learning changes the training procedure to learn a representation that can be easily re-trained to new tasks with small amounts of data. We then demonstrate that meta-learning enables simultaneously training to multiple large organic molecule datasets. As a proof of concept, we examine the performance of a MLIP refit to a small drug-like molecule and show that pre-training potentials to multiple levels of theory with meta-learning improves performance. This difference in performance can be seen both in the reduced error and in the improved smoothness of the potential energy surface produced. We therefore show that meta-learning can utilize existing datasets with inconsistent QM levels of theory to produce models that are better at specializing to new datasets. This opens new routes for creating pre-trained, foundational models for interatomic potentials. △ Less

Submitted 8 July, 2023; originally announced July 2023.

Physically interpretable approximations of many-body spectral functions

Authors: Shubhang Goswami, Kipton Barros, Matthew R. Carbone

Abstract: The rational function approximation provides a natural and interpretable representation of response functions such as the many-body spectral functions. We apply the Vector Fitting (VFIT) algorithm to fit a variety of spectral functions calculated from the Holstein model of electron-phonon interactions. We show that the resulting rational functions are highly efficient in their fitting of sharp fea… ▽ More The rational function approximation provides a natural and interpretable representation of response functions such as the many-body spectral functions. We apply the Vector Fitting (VFIT) algorithm to fit a variety of spectral functions calculated from the Holstein model of electron-phonon interactions. We show that the resulting rational functions are highly efficient in their fitting of sharp features in the spectral functions, and could provide a means to infer physically relevant information from a spectral dataset. The position of the peaks in the approximated spectral function are determined by the location of poles in the complex plane. In addition, we developed an enhanced version of VFIT by introducing a regularization parameter that is slowly annealed to zero. With this new procedure, we demonstrate it is possible to achieve accurate spectral function fits that vary smoothly as a function of physical conditions. △ Less

Submitted 15 February, 2024; v1 submitted 19 June, 2023; originally announced June 2023.

Journal ref: Physical Review E 109.1 (2024): 015302

Machine Learning Framework for Modeling Exciton-Polaritons in Molecular Materials

Authors: Xinyang Li, Nicholas Lubbers, Sergei Tretiak, Kipton Barros, Yu Zhang

Abstract: A light-matter hybrid quasiparticle, called a polariton, is formed when molecules are strongly coupled to an optical cavity. Recent experiments have shown that polariton chemistry can manipulate chemical reactions. Polariton chemistry is a collective phenomenon and its effects increase with the number of molecules in a cavity. However, simulating an ensemble of molecules in the excited state coupl… ▽ More A light-matter hybrid quasiparticle, called a polariton, is formed when molecules are strongly coupled to an optical cavity. Recent experiments have shown that polariton chemistry can manipulate chemical reactions. Polariton chemistry is a collective phenomenon and its effects increase with the number of molecules in a cavity. However, simulating an ensemble of molecules in the excited state coupled to a cavity mode is theoretically and computationally challenging. Recent advances in machine learning techniques have shown promising capabilities in modeling ground state chemical systems. This work presents a general protocol to predict excited-state properties, such as energies, transition dipoles, and non-adiabatic coupling vectors with the hierarchically interacting particle neural network. Machine learning predictions are then applied to compute potential energy surfaces and electronic spectra of a prototype azomethane molecule in the collective coupling scenario. These computational tools provide a much-needed framework to model and understand many molecules' emerging excited-state polariton chemistry. △ Less

Submitted 2 January, 2024; v1 submitted 4 June, 2023; originally announced June 2023.

Renormalized Classical Theory of Quantum Magnets

Authors: David Dahlbom, Hao Zhang, Zoha Laraib, Daniel M. Pajerowski, Kipton Barros, Cristian Batista

Abstract: We derive a renormalized classical spin (RCS) theory for $S > 1/2$ quantum magnets by constraining a generalized classical theory that includes all multipolar fluctuations to a reduced CP$^1$ phase space of dipolar SU($2$) coherent states. When the spin Hamiltonian $\hat{\cal{H}}^{S}$ is linear in the spin operators $\hat{\boldsymbol{S}}_j$ for each lattice site $j$, the RCS Hamiltonian… ▽ More We derive a renormalized classical spin (RCS) theory for $S > 1/2$ quantum magnets by constraining a generalized classical theory that includes all multipolar fluctuations to a reduced CP$^1$ phase space of dipolar SU($2$) coherent states. When the spin Hamiltonian $\hat{\cal{H}}^{S}$ is linear in the spin operators $\hat{\boldsymbol{S}}_j$ for each lattice site $j$, the RCS Hamiltonian $\tilde{\cal{H}}_{\rm cl}$ coincides with the usual classical model $\cal{H}_{\rm cl} = \lim_{S\rightarrow\infty} \hat{\cal{H}}^S$. In the presence of non-linear terms, however, the RCS theory is more accurate than $\cal{H}_{\rm cl}$. For the many materials modeled by spin Hamiltonians with (non-linear) single-ion anisotropy terms, the use of the RCS theory is essential to accurately model phase diagrams and to extract the correct Hamiltonian parameters from neutron scattering data △ Less

Submitted 8 September, 2023; v1 submitted 7 April, 2023; originally announced April 2023.

Comments: 5 pages, 1 figures

Semi-Empirical Shadow Molecular Dynamics: A PyTorch implementation

Authors: Maksim Kulichenko, Kipton Barros, Nicholas Lubbers, Nikita Fedik, Guoqing Zhou, Sergei Tretiak, Benjamin Nebgen, Anders M. N. Niklasson

Abstract: Extended Lagrangian Born-Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov Subspace Approximation for the integration of the electronic equations of motion… ▽ More Extended Lagrangian Born-Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov Subspace Approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. Current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semiempirical quantum mechanical treatment. Applied to molecular dynamics simulation of 840 carbon atoms, one integration time step executes in 4 seconds on a single Nvidia RTX A6000 GPU. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: arXiv admin note: text overlap with arXiv:2003.09050

Reduction of Autocorrelation Times in Lattice Path Integral Quantum Monte Carlo via Direct Sampling of the Truncated Exponential Distribution

Authors: Emanuel Casiano-Diaz, Kipton Barros, Ying Wai Li, Adrian Del Maestro

Abstract: In Monte Carlo simulations, proposed configurations are accepted or rejected according to an acceptance ratio, which depends on an underlying probability distribution and an a priori sampling probability. By carefully selecting the probability distribution from which random variates are sampled, simulations can be made more efficient, by virtue of an autocorrelation time reduction. In this paper,… ▽ More In Monte Carlo simulations, proposed configurations are accepted or rejected according to an acceptance ratio, which depends on an underlying probability distribution and an a priori sampling probability. By carefully selecting the probability distribution from which random variates are sampled, simulations can be made more efficient, by virtue of an autocorrelation time reduction. In this paper, we illustrate how to directly sample random variates from a two dimensional truncated exponential distribution. We show that our direct sampling approach converges faster to the target distribution compared to rejection sampling. The direct sampling of one and two dimensional truncated exponential distributions is then applied to a recent Path Integral Monte Carlo (PIMC) algorithm for the simulation of Bose-Hubbard lattice models at zero temperature. The new sampling method results in improved acceptance ratios and reduced autocorrelation times of estimators, providing an effective speed up of the simulation. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: ~8 pages, 8 figures. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-truncExponSampling

Cubic Double Perovskites Host Noncoplanar Spin Textures

Authors: Joseph A. M. Paddison, Hao Zhang, Jiaqiang Yan, Matthew J. Cliffe, Seung-Hwan Do, Shang Gao, Matthew B. Stone, David Dahlbom, Kipton Barros, Cristian D. Batista, Andrew D. Christianson

Abstract: Magnetic materials with noncoplanar magnetic structures can show unusual physical properties driven by nontrivial topology. Topologically-active states are often multi-q structures, which are challenging to stabilize in models and to identify in materials. Here, we use inelastic neutron-scattering experiments to show that the insulating double perovskites Ba2YRuO6 and Ba2LuRuO6 host a noncoplanar… ▽ More Magnetic materials with noncoplanar magnetic structures can show unusual physical properties driven by nontrivial topology. Topologically-active states are often multi-q structures, which are challenging to stabilize in models and to identify in materials. Here, we use inelastic neutron-scattering experiments to show that the insulating double perovskites Ba2YRuO6 and Ba2LuRuO6 host a noncoplanar 3-q structure on the face-centered cubic lattice. Quantitative analysis of our neutron-scattering data reveals that these 3-q states are stabilized by biquadratic interactions. Our study identifies double perovskites as a highly promising class of materials to realize topological magnetism, elucidates the stabilization mechanism of the 3-q state in these materials, and establishes neutron spectroscopy on powder samples as a valuable technique to distinguish multi-q from single-q states, facilitating the discovery of topologically-nontrivial magnetic materials. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 5 pages, 4 figures

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

Authors: Owen Bradley, Benjamin Cohen-Stead, Steven Johnston, Kipton Barros, Richard T. Scalettar

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

Submitted 11 May, 2023; v1 submitted 12 December, 2022; originally announced December 2022.

Comments: 13 pages, 11 figures. v2: Minimal changes for consistency with published article

Journal ref: npj Quantum Materials 8, 21 (2023)

Lightweight and Effective Tensor Sensitivity for Atomistic Neural Networks

Authors: Michael Chigaev, Justin S. Smith, Steven Anaya, Benjamin Nebgen, Matthew Bettencourt, Kipton Barros, Nicholas Lubbers

Abstract: Atomistic machine learning focuses on the creation of models which obey fundamental symmetries of atomistic configurations, such as permutation, translation, and rotation invariances. In many of these schemes, translation and rotation invariance are achieved by building on scalar invariants, e.g., distances between atom pairs. There is growing interest in molecular representations that work intern… ▽ More Atomistic machine learning focuses on the creation of models which obey fundamental symmetries of atomistic configurations, such as permutation, translation, and rotation invariances. In many of these schemes, translation and rotation invariance are achieved by building on scalar invariants, e.g., distances between atom pairs. There is growing interest in molecular representations that work internally with higher rank rotational tensors, e.g., vector displacements between atoms, and tensor products thereof. Here we present a framework for extending the Hierarchically Interacting Particle Neural Network (HIP-NN) with Tensor Sensitivity information (HIP-NN-TS) from each local atomic environment. Crucially, the method employs a weight tying strategy that allows direct incorporation of many-body information while adding very few model parameters. We show that HIP-NN-TS is more accurate than HIP-NN, with negligible increase in parameter count, for several datasets and network sizes. As the dataset becomes more complex, tensor sensitivities provide greater improvements to model accuracy. In particular, HIP-NN-TS achieves a record mean absolute error of 0.927 kcal/mol for conformational energy variation on the challenging COMP6 benchmark, which includes a broad set of organic molecules. We also compare the computational performance of HIP-NN-TS to HIP-NN and other models in the literature. △ Less

Submitted 26 April, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

A flexible class of exact Hubbard-Stratonovich transformations

Authors: Seher Karakuzu, Benjamin Cohen-Stead, Cristian D. Batista, Steven Johnston, Kipton Barros

Abstract: We consider a class of Hubbard-Stratonovich transformations suitable for treating Hubbard interactions in the context of quantum Monte Carlo simulations. A tunable parameter $p$ allows us to continuously vary from a discrete Ising auxiliary field ($p=\infty$) to a compact auxiliary field that couples to electrons sinusoidally ($p=0$). In tests on the single-band Hubbard model, we find that the sev… ▽ More We consider a class of Hubbard-Stratonovich transformations suitable for treating Hubbard interactions in the context of quantum Monte Carlo simulations. A tunable parameter $p$ allows us to continuously vary from a discrete Ising auxiliary field ($p=\infty$) to a compact auxiliary field that couples to electrons sinusoidally ($p=0$). In tests on the single-band Hubbard model, we find that the severity of the sign problem decreases systematically with increasing $p$. Selecting $p$ finite, however, enables continuous sampling methods like the Langevin or Hamiltonian Monte Carlo methods. We explore the tradeoffs between various simulation methods through numerical benchmarks. △ Less

Submitted 10 May, 2023; v1 submitted 9 November, 2022; originally announced November 2022.

Journal ref: Phys. Rev. E 107, 055301 (2023)

Field-induced spin level crossings within a quasi-XY antiferromagnetic state in Ba$_{2}$FeSi$_{2}$O$_{7}$

Authors: Minseong Lee, Rico Schoenemann, Hao Zhang, David Dahlbom, Tae-Hwan Jang, Seung-Hwan Do, Andrew D. Christianson, Sang-Wook Cheong, Jae-Hoon Park, Eric Brosha, Marcelo Jaime, Kipton Barros, Cristian D. Batista, Vivien S. Zapf

Abstract: We present a high-field study of the strongly anisotropic easy-plane square lattice $S$ = 2 quantum magnet Ba$_{2}$FeSi$_{2}$O$_{7}$. This compound is a rare high-spin antiferromagnetic system with very strong easy-plane anisotropy, such that the interplay between spin level crossings and antiferromagnetic order can be studied. We observe a magnetic field-induced spin level crossing occurring with… ▽ More We present a high-field study of the strongly anisotropic easy-plane square lattice $S$ = 2 quantum magnet Ba$_{2}$FeSi$_{2}$O$_{7}$. This compound is a rare high-spin antiferromagnetic system with very strong easy-plane anisotropy, such that the interplay between spin level crossings and antiferromagnetic order can be studied. We observe a magnetic field-induced spin level crossing occurring within an ordered state. This spin level crossing appears to preserve the magnetic symmetry while producing a non-monotonic dependence the order parameter magnitude. The resulting temperature-magnetic field phase diagram exhibits two dome-shaped regions of magnetic order overlapping around 30 T. The ground state of the lower-field dome is predominantly a linear combination of $| S^{z} = 0 \rangle$ and $ |S^{z} = 1 \rangle$ states, while the ground state of the higher-field dome can be approximated by a linear combination of $| S^{z} = 1 \rangle $ and $ | S^{z} = 2\rangle$ states. At 30 T, where the spin levels cross, the magnetization exhibits a slanted plateau, {\color {black}the magnetocaloric effect shows a broad hump, and the electric polarization shows a weak slope change}. We determined the detailed magnetic phase boundaries and the spin level crossings using measurements of magnetization, electric polarization, and the magnetocaloric effect in pulsed magnetic fields to 60 T. We calculate these properties using a mean field theory based on direct products of SU(5) coherent states and find good agreement. Finally, we measure and calculate the magnetically-induced electric polarization that reflects magnetic ordering and spin level crossings. This multiferroic behavior provides another avenue for detecting phase boundaries and symmetry changes. △ Less

Submitted 14 March, 2023; v1 submitted 25 October, 2022; originally announced October 2022.

Comments: 9 pages, 5 figures

Report number: LA-UR-22-31305

Langevin dynamics of generalized spins as SU($N$) coherent states

Authors: David Dahlbom, Cole Miles, Hao Zhang, Cristian D. Batista, Kipton Barros

Abstract: Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude $S$ have $N=2S+1$ levels. Alternatively, the spin state is fully characterized by a set of $N^{2}-1$ local physical observables, which we interpret as generalized spin components. For example, a spin with $S=1$ has three dipole compon… ▽ More Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude $S$ have $N=2S+1$ levels. Alternatively, the spin state is fully characterized by a set of $N^{2}-1$ local physical observables, which we interpret as generalized spin components. For example, a spin with $S=1$ has three dipole components and five quadrupole components. These components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU($N$) coherent states in the Schrödinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a non-equilibrium relaxation process that creates CP$^{2}$ Skyrmions, which are topological defects with both dipole and quadrupole character. △ Less

Submitted 29 December, 2022; v1 submitted 2 September, 2022; originally announced September 2022.

A hybrid Monte Carlo study of bond-stretching electron-phonon interactions and charge order in BaBiO$_3$

Authors: Benjamin Cohen-Stead, Kipton Barros, Richard Scalettar, Steven Johnston

Abstract: The relationship between electron-phonon ($e$-ph) interactions and charge-density-wave (CDW) order in the bismuthate family of high-temperature superconductors remains unresolved. We address this question using nonperturbative hybrid Monte Carlo calculations for the parent compound BaBiO$_3$. Our model includes the Bi $6s$ and O $2p_σしぐま$ orbitals and coupling to the Bi-O bond-stretching branch of op… ▽ More The relationship between electron-phonon ($e$-ph) interactions and charge-density-wave (CDW) order in the bismuthate family of high-temperature superconductors remains unresolved. We address this question using nonperturbative hybrid Monte Carlo calculations for the parent compound BaBiO$_3$. Our model includes the Bi $6s$ and O $2p_σしぐま$ orbitals and coupling to the Bi-O bond-stretching branch of optical phonons via modulations of the Bi-O hopping integral. We simulate three-dimensional clusters of up to 4000 orbitals, with input model parameters taken from {\it ab initio} electronic structure calculations and a phonon energy $\hbarΩおめが_0 = 60$~meV. Our results demonstrate that the coupling to the bond-stretching modes is sufficient to reproduce the CDW transition in this system, despite a relatively small dimensionless coupling. We also find that the transition deviates from the weak-coupling Peierls' picture. This work demonstrates that off-diagonal $e$-ph interactions in orbital space are vital in establishing the bismuthate phase diagram. △ Less

Submitted 10 May, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Journal ref: npj Computational Materials 9, 40 (2023)

Understanding temperature-dependent SU($3$) spin dynamics in the $S=1$ antiferromagnet Ba$_2$FeSi$_2$O$_7$

Authors: Seung-Hwan Do, Hao Zhang, David A. Dahlbom, Travis J. Williams, V. Ovidiu Garlea, Tao Hong, Tae-Hwan Jang, Sang-Wook Cheong, Jae-Hoon Park, Kipton Barros, Cristian D. Batista, Andrew D. Christianson

Abstract: Quantum magnets admit more than one classical limit and $N$-level systems with strong single-ion anisotropy are expected to be described by a classical approximation based on SU($N$) coherent states. Here we test this hypothesis by modeling finite temperature inelastic neutron scattering (INS) data of the effective spin-one antiferromagnet \bfso{}. The measured dynamic structure factor is calculat… ▽ More Quantum magnets admit more than one classical limit and $N$-level systems with strong single-ion anisotropy are expected to be described by a classical approximation based on SU($N$) coherent states. Here we test this hypothesis by modeling finite temperature inelastic neutron scattering (INS) data of the effective spin-one antiferromagnet \bfso{}. The measured dynamic structure factor is calculated with a generalized Landau-Lifshitz dynamics for SU($3$) spins. Unlike the traditional classical limit based on SU($2$) coherent states, the results obtained with classical SU($3$) spins are in good agreement with the measured temperature-dependent spectrum. The SU($3$) approach developed here provides a general framework to understand the broad class of materials comprising weakly coupled antiferromagnetic dimers, trimers, or tetramers, and magnets with strong single-ion anisotropy. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 8 pages 4 figures

Geometric integration of classical spin dynamics via a mean-field Schrödinger equation

Authors: David Dahlbom, Hao Zhang, Cole Miles, Xiaojian Bai, Cristian D. Batista, Kipton Barros

Abstract: The Landau-Lifshitz equation describes the time-evolution of magnetic dipoles, and can be derived by taking the classical limit of a quantum mechanical spin Hamiltonian. To take this limit, one constrains the many-body quantum state to a tensor product of coherent states, thereby neglecting entanglement between sites. Expectation values of the quantum spin operators produce the usual classical spi… ▽ More The Landau-Lifshitz equation describes the time-evolution of magnetic dipoles, and can be derived by taking the classical limit of a quantum mechanical spin Hamiltonian. To take this limit, one constrains the many-body quantum state to a tensor product of coherent states, thereby neglecting entanglement between sites. Expectation values of the quantum spin operators produce the usual classical spin dipoles. One may also consider expectation values of polynomials of the spin operators, leading to quadrupole and higher-order spin moments, which satisfy a dynamical equation of motion that generalizes the Landau-Lifshitz dynamics [Zhang and Batista, Phys. Rev. B 104, 104409 (2021)]. Here, we reformulate the dynamics of these $N^2-1$ generalized spin components as a mean-field Schrödinger equation on the $N$-dimensional coherent state. This viewpoint suggests efficient integration methods that respect the local symplectic structure of the classical spin dynamics. △ Less

Submitted 22 July, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

CP$^2$ Skyrmions and Skyrmion Crystals in Realistic Quantum Magnets

Authors: Hao Zhang, Zhentao Wang, David Dahlbom, Kipton Barros, Cristian D. Batista

Abstract: Magnetic skyrmions are nanoscale topological textures that have been recently observed in different families of quantum magnets. These textures are known as CP$^1$ skyrmions because the target manifold of the magnetization field is the 2D sphere $S^2 \cong$ CP$^{1}$. Here we report the emergence of magnetic CP$^2$ skyrmions in realistic spin-$1$ models. Unlike CP$^1$ skyrmions, CP$^2$ skyrmions ca… ▽ More Magnetic skyrmions are nanoscale topological textures that have been recently observed in different families of quantum magnets. These textures are known as CP$^1$ skyrmions because the target manifold of the magnetization field is the 2D sphere $S^2 \cong$ CP$^{1}$. Here we report the emergence of magnetic CP$^2$ skyrmions in realistic spin-$1$ models. Unlike CP$^1$ skyrmions, CP$^2$ skyrmions can also arise as metastable textures of quantum paramagnets, opening a new road to discover emergent topological solitons in non-magnetic materials. The quantum phase diagram of the spin-$1$ models also includes magnetic field-induced CP$^2$ skyrmion crystals that can be detected with regular momentum- (diffraction) and real-space (Lorentz transmission electron microscopy) experimental techniques. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: main text: 8 pages, 4 figures. Supplemental Info: 6 pages, 3 figures

Journal ref: Nat Commun 14, 3626 (2023)

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

Authors: Benjamin Cohen-Stead, Owen Bradley, Cole Miles, George Batrouni, Richard Scalettar, Kipton Barros

Abstract: We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Altern… ▽ More We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details. △ Less

Submitted 15 July, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

Authors: Cole Miles, Benjamin Cohen-Stead, Owen Bradley, Steven Johnston, Richard Scalettar, Kipton Barros

Abstract: We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle num… ▽ More We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. We perform a variety of tests, and in all cases find rapid convergence of the chemical potential -- inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations. △ Less

Submitted 14 April, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 10 pages, 6 figures

Machine Learning of consistent thermodynamic models using automatic differentiation

Authors: David Rosenberger, Kipton Barros, Timothy C. Germann, Nicholas Lubbers

Abstract: We propose a data-driven method to describe consistent equations of state (EOS) for arbitrary systems. Complex EOS are traditionally obtained by fitting suitable analytical expressions to thermophysical data. A key aspect of EOS are that the relationships between state variables are given by derivatives of the system free energy. In this work, we model the free energy with an artificial neural net… ▽ More We propose a data-driven method to describe consistent equations of state (EOS) for arbitrary systems. Complex EOS are traditionally obtained by fitting suitable analytical expressions to thermophysical data. A key aspect of EOS are that the relationships between state variables are given by derivatives of the system free energy. In this work, we model the free energy with an artificial neural network, and utilize automatic differentiation to directly learn the derivatives of the free energy. We demonstrate this approach on two different systems, the analytic van der Waals EOS, and published data for the Lennard-Jones fluid, and show that it is advantageous over direct learning of thermodynamic properties (i.e. not as derivatives of the free energy, but as independent properties), in terms of both accuracy and the exact preservation of the Maxwell relations. Furthermore, the method implicitly provides the free energy of a system without explicit integration. △ Less

Submitted 9 March, 2022; v1 submitted 10 August, 2021; originally announced August 2021.

Journal ref: Phys. Rev. E 105, 045301, 2022

Machine learning of Kondo physics using variational autoencoders and symbolic regression

Authors: Cole Miles, Matthew R. Carbone, Erica J. Sturm, Deyu Lu, Andreas Weichselbaum, Kipton Barros, Robert M. Konik

Abstract: We employ variational autoencoders to extract physical insight from a dataset of one-particle Anderson impurity model spectral functions. Autoencoders are trained to find a low-dimensional, latent space representation that faithfully characterizes each element of the training set, as measured by a reconstruction error. Variational autoencoders, a probabilistic generalization of standard autoencode… ▽ More We employ variational autoencoders to extract physical insight from a dataset of one-particle Anderson impurity model spectral functions. Autoencoders are trained to find a low-dimensional, latent space representation that faithfully characterizes each element of the training set, as measured by a reconstruction error. Variational autoencoders, a probabilistic generalization of standard autoencoders, further condition the learned latent space to promote highly interpretable features. In our study, we find that the learned latent variables strongly correlate with well known, but nontrivial, parameters that characterize emergent behaviors in the Anderson impurity model. In particular, one latent variable correlates with particle-hole asymmetry, while another is in near one-to-one correspondence with the Kondo temperature, a dynamically generated low-energy scale in the impurity model. Using symbolic regression, we model this variable as a function of the known bare physical input parameters and "rediscover" the non-perturbative formula for the Kondo temperature. The machine learning pipeline we develop suggests a general purpose approach which opens opportunities to discover new domain knowledge in other physical systems. △ Less

Submitted 16 December, 2021; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: Update to match PRB publication + typo fixes + minor edits

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

Quantum-based Molecular Dynamics Simulations Using Tensor Cores

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

Abstract: Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency, but with the caveat that they were designed for tensor contractions (matrix-matrix multiplications) using only low-precision f… ▽ More Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency, but with the caveat that they were designed for tensor contractions (matrix-matrix multiplications) using only low-precision floating point operations. In spite of this, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born-Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on the tensor cores of a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born-Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics. △ Less

Submitted 10 September, 2021; v1 submitted 6 July, 2021; originally announced July 2021.

Comments: 20 pages, 9 figures

Report number: LA-UR-21-25427

Quantum Monte Carlo study of an anharmonic Holstein model

Authors: G. Paleari, F. Hébert, B. Cohen-Stead, K. Barros, R. T. Scalettar, G. G. Batrouni

Abstract: We study the effects of anharmonicity on the physics of the Holstein model, which describes the coupling of itinerant fermions and localized quantum phonons, by introducing a quartic term in the phonon potential energy. We find that the presence of this anharmonic term reduces the extent of the charge density wave phase (CDW) at half-filling as well as the transition temperature to this phase. Dop… ▽ More We study the effects of anharmonicity on the physics of the Holstein model, which describes the coupling of itinerant fermions and localized quantum phonons, by introducing a quartic term in the phonon potential energy. We find that the presence of this anharmonic term reduces the extent of the charge density wave phase (CDW) at half-filling as well as the transition temperature to this phase. Doping away from half-filling, we observe a first order phase transition between the CDW and a homogeneous phase which is also present in the harmonic model. In addition, we study the evolution of the superconducting susceptibility in the doped region and show that anharmonicity can enhance the superconducting response. △ Less

Submitted 31 May, 2021; v1 submitted 20 January, 2021; originally announced January 2021.

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

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

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

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

Submitted 16 January, 2021; originally announced January 2021.

Comments: 10 pages, 6 figures

Report number: LA-UR-21-20350

Phase Diagram of a Spin-ice Kondo Lattice Model in a Breathing Pyrochlore Lattice

Authors: Munir Shahzad, Kipton Barros, Stephanie H. Curnoe

Abstract: We study a spin-ice Kondo lattice model on a breathing pyrochlore lattice with classical localized spins. The highly efficient kernel polynomial expansion method, together with a classical Monte Carlo method, is employed in order to study the magnetic phase diagram at four representative values of the number density of itinerant electrons. We tune the breathing mode by varying the hopping ratio --… ▽ More We study a spin-ice Kondo lattice model on a breathing pyrochlore lattice with classical localized spins. The highly efficient kernel polynomial expansion method, together with a classical Monte Carlo method, is employed in order to study the magnetic phase diagram at four representative values of the number density of itinerant electrons. We tune the breathing mode by varying the hopping ratio -- the ratio of hopping parameters for itinerant electrons along inequivalent paths. Several interesting magnetic phases are stabilized in the phase diagram parameterized by the hopping ratio, Kondo coupling, and electronic filling fraction, including an "all-in/all-out" ordered spin configuration phase, spin-ice, ordered phases containing $16$ and $32$ spin sites in the magnetic unit cell, as well as a disordered phase at small values of the hopping ratio. △ Less

Submitted 30 July, 2020; originally announced July 2020.

Comments: 9 pages and 9 figures. Any comments or suggestions will be appreciated

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

Simple and efficient algorithms for training machine learning potentials to force data

Authors: Justin S. Smith, Nicholas Lubbers, Aidan P. Thompson, Kipton Barros

Abstract: Abstract Machine learning models, trained on data from ab initio quantum simulations, are yielding molecular dynamics potentials with unprecedented accuracy. One limiting factor is the quantity of available training data, which can be expensive to obtain. A quantum simulation often provides all atomic forces, in addition to the total energy of the system. These forces provide much more information… ▽ More Abstract Machine learning models, trained on data from ab initio quantum simulations, are yielding molecular dynamics potentials with unprecedented accuracy. One limiting factor is the quantity of available training data, which can be expensive to obtain. A quantum simulation often provides all atomic forces, in addition to the total energy of the system. These forces provide much more information than the energy alone. It may appear that training a model to this large quantity of force data would introduce significant computational costs. Actually, training to all available force data should only be a few times more expensive than training to energies alone. Here, we present a new algorithm for efficient force training, and benchmark its accuracy by training to forces from real-world datasets for organic chemistry and bulk aluminum. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Langevin Simulations of the Half-Filled Cubic Holstein Model

Authors: B. Cohen-Stead, Kipton Barros, Zi Yang Meng, Chuang Chen, R. T Scalettar, G. G. Batrouni

Abstract: Over the past several years, reliable Quantum Monte Carlo results for the charge density wave transition temperature $T_{cdw}$ of the half-filled two dimensional Holstein model in square and honeycomb lattices have become available for the first time. Exploiting the further development of numerical methodology, here we present results in three dimensions, which are made possible through the use of… ▽ More Over the past several years, reliable Quantum Monte Carlo results for the charge density wave transition temperature $T_{cdw}$ of the half-filled two dimensional Holstein model in square and honeycomb lattices have become available for the first time. Exploiting the further development of numerical methodology, here we present results in three dimensions, which are made possible through the use of Langevin evolution of the quantum phonon degrees of freedom. In addition to determining $T_{cdw}$ from the scaling of the charge correlations, we also examine the nature of charge order at general wave vectors for different temperatures, couplings, and phonon frequencies, and the behavior of the spectral function and specific heat. △ Less

Submitted 2 May, 2020; originally announced May 2020.

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

Automated discovery of a robust interatomic potential for aluminum

Authors: Justin S. Smith, Benjamin Nebgen, Nithin Mathew, Jie Chen, Nicholas Lubbers, Leonid Burakovsky, Sergei Tretiak, Hai Ah Nam, Timothy Germann, Saryu Fensin, Kipton Barros

Abstract: Accuracy of molecular dynamics simulations depends crucially on the interatomic potential used to generate forces. The gold standard would be first-principles quantum mechanics (QM) calculations, but these become prohibitively expensive at large simulation scales. Machine learning (ML) based potentials aim for faithful emulation of QM at drastically reduced computational cost. The accuracy and rob… ▽ More Accuracy of molecular dynamics simulations depends crucially on the interatomic potential used to generate forces. The gold standard would be first-principles quantum mechanics (QM) calculations, but these become prohibitively expensive at large simulation scales. Machine learning (ML) based potentials aim for faithful emulation of QM at drastically reduced computational cost. The accuracy and robustness of an ML potential is primarily limited by the quality and diversity of the training dataset. Using the principles of active learning (AL), we present a highly automated approach to dataset construction. The strategy is to use the ML potential under development to sample new atomic configurations and, whenever a configuration is reached for which the ML uncertainty is sufficiently large, collect new QM data. Here, we seek to push the limits of automation, removing as much expert knowledge from the AL process as possible. All sampling is performed using MD simulations starting from an initially disordered configuration, and undergoing non-equilibrium dynamics as driven by time-varying applied temperatures. We demonstrate this approach by building an ML potential for aluminum (ANI-Al). After many AL iterations, ANI-Al teaches itself to predict properties like the radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. To demonstrate transferability, we perform a 1.3M atom shock simulation, and show that ANI-Al predictions agree very well with DFT calculations on local atomic environments sampled from the nonequilibrium dynamics. Interestingly, the configurations appearing in shock appear to have been well sampled in the AL training dataset, in a way that we illustrate visually. △ Less

Submitted 24 August, 2020; v1 submitted 10 March, 2020; originally announced March 2020.

Fast and stable deep-learning predictions of material properties for solid solution alloys

Authors: Massimiliano Lupo Pasini, Ying Wai Li, Junqi Yin, Jiaxin Zhang, Kipton Barros, Markus Eisenbach

Abstract: We present a novel deep learning (DL) approach to produce highly accurate predictions of macroscopic physical properties of solid solution binary alloys and magnetic systems. The major idea is to make use of the correlations between different physical properties in alloy systems to improve the prediction accuracy of neural network (NN) models. We use multitasking NN models to simultaneously predic… ▽ More We present a novel deep learning (DL) approach to produce highly accurate predictions of macroscopic physical properties of solid solution binary alloys and magnetic systems. The major idea is to make use of the correlations between different physical properties in alloy systems to improve the prediction accuracy of neural network (NN) models. We use multitasking NN models to simultaneously predict the total energy, charge density and magnetic moment. These physical properties mutually serve as constraints during the training of the multitasking NN, resulting in more reliable DL models because multiple physics properties are correctly learned by a single model. Two binary alloys, copper-gold (CuAu) and iron-platinum (FePt), were studied. Our results show that once the multitasking NN's are trained, they can estimate the material properties for a specific configuration hundreds of times faster than first-principles density functional theory calculations while retaining comparable accuracy. We used a simple measure based on the root-mean-squared errors (RMSE) to quantify the quality of the NN models, and found that the inclusion of charge density and magnetic moment as physical constraints leads to more stable models that exhibit improved accuracy and reduced uncertainty for the energy predictions. △ Less

Submitted 6 November, 2020; v1 submitted 18 December, 2019; originally announced December 2019.

Comments: 17 pages, 8 figures, 3 tables

Journal ref: J Phys Condens Matter . 2020 Nov 17

Anomalous Magnetoresistance due to Longitudinal Spin Fluctuations in a Jeff = 1/2 Mott Semiconductor

Authors: Lin Hao, Zhentao Wang, Junyi Yang, D. Meyers, Joshua Sanchez, Gilberto Fabbris, Yongseong Choi, Jong-Woo Kim, Daniel Haskel, Philip J. Ryan, Kipton Barros, Jiun-Haw Chu, M. P. M. Dean, Cristian D. Batista, Jian Liu

Abstract: As a hallmark of electronic correlation, spin-charge interplay underlies many emergent phenomena in doped Mott insulators, such as high-temperature superconductivity, whereas the half-filled parent state is usually electronically frozen with an antiferromagnetic order that resists external control. We report on the observation of a new positive magnetoresistance that probes the staggered susceptib… ▽ More As a hallmark of electronic correlation, spin-charge interplay underlies many emergent phenomena in doped Mott insulators, such as high-temperature superconductivity, whereas the half-filled parent state is usually electronically frozen with an antiferromagnetic order that resists external control. We report on the observation of a new positive magnetoresistance that probes the staggered susceptibility of a pseudospin-half square-lattice Mott insulator built as an artificial SrIrO3/SrTiO3 superlattice. Its size is particularly large in the high-temperature insulating paramagnetic phase near the Néel transition. This novel magnetoresistance originates from a collective charge response to the large longitudinal spin fluctuations under a linear coupling between the external magnetic field and the staggered magnetization enabled by strong spin-orbit interaction. Our results demonstrate a magnetic control of the binding energy of the fluctuating particle-hole pairs in the Slater-Mott crossover regime analogous to the BCS-to-Bose-Einstein condensation crossover of ultracold-superfluids. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: Accepted by Nature communications

Journal ref: Nature Communications (2019) 10:5301

Machine Learned Hückel Theory: Interfacing Physics and Deep Neural Networks

Authors: Tetiana Zubatyuk, Ben Nebgen, Nicholas Lubbers, Justin S. Smith, Roman Zubatyuk, Guoqing Zhou, Christopher Koh, Kipton Barros, Olexandr Isayev, Sergei Tretiak

Abstract: The Hückel Hamiltonian is an incredibly simple tight-binding model famed for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions… ▽ More The Hückel Hamiltonian is an incredibly simple tight-binding model famed for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these traditionally static parameters with dynamically predicted values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability while the deep neural network parameterization is smooth, accurate, and reproduces insightful features of the original static parameterization. Finally, we demonstrate that the Hückel model, and not the deep neural network, is responsible for capturing intricate orbital interactions in two molecular case studies. Overall, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Machine Learning Assisted Insight to Spin Ice Dy$_2$Ti$_2$O$_7$

Authors: Anjana M Samarakoon, Kipton Barros, Ying Wai Li, Markus Eisenbach, Qiang Zhang, Feng Ye, Z. L. Dun, Haidong Zhou, Santiago A. Grigera, Cristian D. Batista, D. Alan Tennant

Abstract: Complex behavior poses challenges in extracting models from experiment. An example is spin liquid formation in frustrated magnets like Dy$_2$Ti$_2$O$_7$. Understanding has been hindered by issues including disorder, glass formation, and interpretation of scattering data. Here, we use a novel automated capability to extract model Hamiltonians from data, and to identify different magnetic regimes. T… ▽ More Complex behavior poses challenges in extracting models from experiment. An example is spin liquid formation in frustrated magnets like Dy$_2$Ti$_2$O$_7$. Understanding has been hindered by issues including disorder, glass formation, and interpretation of scattering data. Here, we use a novel automated capability to extract model Hamiltonians from data, and to identify different magnetic regimes. This involves training an autoencoder to learn a compressed representation of three-dimensional diffuse scattering, over a wide range of spin Hamiltonians. The autoencoder finds optimal matches according to scattering and heat capacity data and provides confidence intervals. Validation tests indicate that our optimal Hamiltonian accurately predicts temperature and field dependence of both magnetic structure and magnetization, as well as glass formation and irreversibility in Dy$_2$Ti$_2$O$_7$. The autoencoder can also categorize different magnetic behaviors and eliminate background noise and artifacts in raw data. Our methodology is readily applicable to other materials and types of scattering problems. △ Less

Submitted 11 November, 2020; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: 18 pages, 6 figures

Machine learning for molecular dynamics with strongly correlated electrons

Authors: Hidemaro Suwa, Justin S. Smith, Nicholas Lubbers, Cristian D. Batista, Gia-Wei Chern, Kipton Barros

Abstract: We use machine learning to enable large-scale molecular dynamics (MD) of a correlated electron model under the Gutzwiller approximation scheme. This model exhibits a Mott transition as a function of on-site Coulomb repulsion $U$. The repeated solution of the Gutzwiller self-consistency equations would be prohibitively expensive for large-scale MD simulations. We show that machine learning models o… ▽ More We use machine learning to enable large-scale molecular dynamics (MD) of a correlated electron model under the Gutzwiller approximation scheme. This model exhibits a Mott transition as a function of on-site Coulomb repulsion $U$. The repeated solution of the Gutzwiller self-consistency equations would be prohibitively expensive for large-scale MD simulations. We show that machine learning models of the Gutzwiller potential energy can be remarkably accurate. The models, which are trained with $N=33$ atoms, enable highly accurate MD simulations at much larger scales ($N\gtrsim10^{3}$). We investigate the physics of the smooth Mott crossover in the fluid phase. △ Less

Submitted 10 April, 2019; v1 submitted 5 November, 2018; originally announced November 2018.

Comments: 6 pages, 4 figures

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

Earthquake catalog-based machine learning identification of laboratory fault states and the effects of magnitude of completeness

Authors: Nicholas Lubbers, David C. Bolton, Jamaludin Mohd-Yusof, Chris Marone, Kipton Barros, Paul A. Johnson

Abstract: Machine learning regression can predict macroscopic fault properties such as shear stress, friction, and time to failure using continuous records of fault zone acoustic emissions. Here we show that a similar approach is successful using event catalogs derived from the continuous data. Our methods are applicable to catalogs of arbitrary scale and magnitude of completeness. We investigate how machin… ▽ More Machine learning regression can predict macroscopic fault properties such as shear stress, friction, and time to failure using continuous records of fault zone acoustic emissions. Here we show that a similar approach is successful using event catalogs derived from the continuous data. Our methods are applicable to catalogs of arbitrary scale and magnitude of completeness. We investigate how machine learning regression from an event catalog of laboratory earthquakes performs as a function of the catalog magnitude of completeness. We find that strong model performance requires a sufficiently low magnitude of completeness, and below this magnitude of completeness model performance saturates. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Transferable Molecular Charge Assignment Using Deep Neural Networks

Authors: Ben Nebgen, Nick Lubbers, Justin S. Smith, Andrew Sifain, Andrey Lokhov, Olexandr Isayev, Adrian Roitberg, Kipton Barros, Sergei Tretiak

Abstract: We use HIP-NN, a neural network architecture that excels at predicting molecular energies, to predict atomic charges. The charge predictions are accurate over a wide range of molecules (both small and large) and for a diverse set of charge assignment schemes. To demonstrate the power of charge prediction on non-equilibrium geometries, we use HIP-NN to generate IR spectra from dynamical trajectorie… ▽ More We use HIP-NN, a neural network architecture that excels at predicting molecular energies, to predict atomic charges. The charge predictions are accurate over a wide range of molecules (both small and large) and for a diverse set of charge assignment schemes. To demonstrate the power of charge prediction on non-equilibrium geometries, we use HIP-NN to generate IR spectra from dynamical trajectories on a variety of molecules. The results are in good agreement with reference IR spectra produced by traditional theoretical methods. Critically, for this application, HIP-NN charge predictions are about 104 times faster than direct DFT charge calculations. Thus, ML provides a pathway to greatly increase the range of feasible simulations while retaining quantum-level accuracy. In summary, our results provide further evidence that machine learning can replicate high-level quantum calculations at a tiny fraction of the computational cost. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Nonequilibrium dynamics of superconductivity in the attractive Hubbard model

Authors: Gia-Wei Chern, Kipton Barros

Abstract: We present a framework of semiclassical superconductivity (SC) dynamics that properly includes effects of spatial fluctuations for the attractive Hubbard model. We consider both coherent and adiabatic limits. To model the coherent SC dynamics, we develop a real-space von~Neumann equation based on the time-dependent Hartree-Fock-Bogoliubov theory. Applying our method to interaction quenches in the… ▽ More We present a framework of semiclassical superconductivity (SC) dynamics that properly includes effects of spatial fluctuations for the attractive Hubbard model. We consider both coherent and adiabatic limits. To model the coherent SC dynamics, we develop a real-space von~Neumann equation based on the time-dependent Hartree-Fock-Bogoliubov theory. Applying our method to interaction quenches in the negative-$U$ Hubbard model, we show that the relaxation of SC order at weak coupling is dominated by Landau-damping. At strong coupling, we find a two-stage relaxation of the pairing field: a collapse of the synchronized oscillation of Cooper pairs due to spatial inhomogeneity, followed by a slow relaxation to a quasi-stationary state. SC dynamics in adiabatic limit is described by a quantum Landau-Lifshitz equation with Ginzburg-Landau relaxation. Numerical simulations of the pump-probe process show that long time recovery of the pairing field is dominated by defects dynamics. Our results demonstrate the important role of spatial fluctuations in both limits. △ Less

Submitted 29 October, 2018; v1 submitted 12 March, 2018; originally announced March 2018.

Comments: 6 pages, 4 figures

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

Gradient-based stochastic estimation of the density matrix

Authors: Zhentao Wang, Gia-Wei Chern, Cristian D. Batista, Kipton Barros

Abstract: Fast estimation of the single-particle density matrix is key to many applications in quantum chemistry and condensed matter physics. The best numerical methods leverage the fact that the density matrix elements $f(H)_{ij}$ decay rapidly with distance $r_{ij}$ between orbitals. This decay is usually exponential. However, for the special case of metals at zero temperature, algebraic decay of the den… ▽ More Fast estimation of the single-particle density matrix is key to many applications in quantum chemistry and condensed matter physics. The best numerical methods leverage the fact that the density matrix elements $f(H)_{ij}$ decay rapidly with distance $r_{ij}$ between orbitals. This decay is usually exponential. However, for the special case of metals at zero temperature, algebraic decay of the density matrix appears and poses a significant numerical challenge. We introduce a gradient-based probing method to estimate all local density matrix elements at a computational cost that scales linearly with system size. For zero-temperature metals the stochastic error scales like $S^{-(d+2)/2d}$, where $d$ is the dimension and $S$ is a prefactor to the computational cost. The convergence becomes exponential if the system is at finite temperature or is insulating. △ Less

Submitted 13 April, 2018; v1 submitted 28 November, 2017; originally announced November 2017.

Hierarchical modeling of molecular energies using a deep neural network

Authors: Nicholas Lubbers, Justin S. Smith, Kipton Barros

Abstract: We introduce the Hierarchically Interacting Particle Neural Network (HIP-NN) to model molecular properties from datasets of quantum calculations. Inspired by a many-body expansion, HIP-NN decomposes properties, such as energy, as a sum over hierarchical terms. These terms are generated from a neural network--a composition of many nonlinear transformations--acting on a representation of the molecul… ▽ More We introduce the Hierarchically Interacting Particle Neural Network (HIP-NN) to model molecular properties from datasets of quantum calculations. Inspired by a many-body expansion, HIP-NN decomposes properties, such as energy, as a sum over hierarchical terms. These terms are generated from a neural network--a composition of many nonlinear transformations--acting on a representation of the molecule. HIP-NN achieves state-of-the-art performance on a dataset of 131k ground state organic molecules, and predicts energies with 0.26 kcal/mol mean absolute error. With minimal tuning, our model is also competitive on a dataset of molecular dynamics trajectories. In addition to enabling accurate energy predictions, the hierarchical structure of HIP-NN helps to identify regions of model uncertainty. △ Less

Submitted 29 September, 2017; originally announced October 2017.

Semiclassical dynamics of spin density waves

Authors: Gia-Wei Chern, Kipton Barros, Zhentao Wang, Hidemaro Suwa, Cristian D. Batista

Abstract: We present a theoretical framework for equilibrium and nonequilibrium dynamical simulation of quantum states with spin-density-wave (SDW) order. Within a semiclassical adiabatic approximation that retains electron degrees of freedom, we demonstrate that the SDW order parameter obeys a generalized Landau-Lifshitz equation. With the aid of an enhanced kernel polynomial method, our linear-scaling qua… ▽ More We present a theoretical framework for equilibrium and nonequilibrium dynamical simulation of quantum states with spin-density-wave (SDW) order. Within a semiclassical adiabatic approximation that retains electron degrees of freedom, we demonstrate that the SDW order parameter obeys a generalized Landau-Lifshitz equation. With the aid of an enhanced kernel polynomial method, our linear-scaling quantum Landau-Lifshitz dynamics (QLLD) method enables dynamical SDW simulations with $N \simeq 10^5$ lattice sites. Our real-space formulation can be used to compute dynamical responses, such as dynamical structure factor, of complex and even inhomogeneous SDW configurations at zero or finite temperatures. Applying the QLLD to study the relaxation of a noncoplanar topological SDW under the excitation of a short pulse, we further demonstrate the crucial role of spatial correlations and fluctuations in the SDW dynamics. △ Less

Submitted 13 December, 2017; v1 submitted 27 August, 2017; originally announced August 2017.

Comments: 17 pages, 5 figures

Journal ref: Phys. Rev. B 97, 035120 (2018)

Domain engineering of the metastable domains in the 4f-uniaxial-ferromagnet CeRu$_2$Ga$_2$B

Authors: D. Wulferding, H. Kim, I. Yang, J. Jeong, K. Barros, Y. Kato, I. Martin, O. E. Ayala-Valenzuela, M. Lee, H. C. Choi, F. Ronning, L. Civale, R. E. Baumbach, E. D. Bauer, J. D. Thompson, R. Movshovich, J. Kim

Abstract: In search of novel, improved materials for magnetic data storage and spintronic devices, compounds that allow a tailoring of magnetic domain shapes and sizes are essential. Good candidates are materials with intrinsic anisotropies or competing interactions, as they are prone to host various domain phases that can be easily and precisely selected by external tuning parameters such as temperature an… ▽ More In search of novel, improved materials for magnetic data storage and spintronic devices, compounds that allow a tailoring of magnetic domain shapes and sizes are essential. Good candidates are materials with intrinsic anisotropies or competing interactions, as they are prone to host various domain phases that can be easily and precisely selected by external tuning parameters such as temperature and magnetic field. Here, we utilize vector magnetic fields to visualize directly the magnetic anisotropy in the uniaxial ferromagnet CeRu$_2$Ga$_2$B. We demonstrate a feasible control both globally and locally of domain shapes and sizes by the external field as well as a smooth transition from single stripe to bubble domains, which opens the door to future applications based on magnetic domain tailoring. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Journal ref: Sci. Rep. 7, 46296 (2017)

Machine Learning Predicts Laboratory Earthquakes

Authors: Bertrand Rouet-Leduc, Claudia Hulbert, Nicholas Lubbers, Kipton Barros, Colin Humphreys, Paul A. Johnson

Abstract: Forecasting fault failure is a fundamental but elusive goal in earthquake science. Here we show that by listening to the acoustic signal emitted by a laboratory fault, machine learning can predict the time remaining before it fails with great accuracy. These predictions are based solely on the instantaneous physical characteristics of the acoustical signal, and do not make use of its history. Surp… ▽ More Forecasting fault failure is a fundamental but elusive goal in earthquake science. Here we show that by listening to the acoustic signal emitted by a laboratory fault, machine learning can predict the time remaining before it fails with great accuracy. These predictions are based solely on the instantaneous physical characteristics of the acoustical signal, and do not make use of its history. Surprisingly, machine learning identifies a signal emitted from the fault zone previously thought to be low-amplitude noise that enables failure forecasting throughout the laboratory quake cycle. We hypothesize that applying this approach to continuous seismic data may lead to significant advances in identifying currently unknown signals, in providing new insights into fault physics, and in placing bounds on fault failure times. △ Less

Submitted 19 February, 2017; originally announced February 2017.

Comments: 17 pages, 4 figures

Report number: LAUR 1808407

Learning molecular energies using localized graph kernels

Authors: G. Ferré, T. Haut, K. Barros

Abstract: Recent machine learning methods make it possible to model potential energy of atomic configurations with chemical-level accuracy (as calculated from ab-initio calculations) and at speeds suitable for molecular dynam- ics simulation. Best performance is achieved when the known physical constraints are encoded in the machine learning models. For example, the atomic energy is invariant under global t… ▽ More Recent machine learning methods make it possible to model potential energy of atomic configurations with chemical-level accuracy (as calculated from ab-initio calculations) and at speeds suitable for molecular dynam- ics simulation. Best performance is achieved when the known physical constraints are encoded in the machine learning models. For example, the atomic energy is invariant under global translations and rotations, it is also invariant to permutations of same-species atoms. Although simple to state, these symmetries are complicated to encode into machine learning algorithms. In this paper, we present a machine learning approach based on graph theory that naturally incorporates translation, rotation, and permutation symmetries. Specifically, we use a random walk graph kernel to measure the similarity of two adjacency matrices, each of which represents a local atomic environment. This Graph Approximated Energy (GRAPE) approach is flexible and admits many possible extensions. We benchmark a simple version of GRAPE by predicting atomization energies on a standard dataset of organic molecules. △ Less

Submitted 25 April, 2017; v1 submitted 1 December, 2016; originally announced December 2016.

Journal ref: The Journal of Chemical Physics, 146(11), 114107 (2017)

Inferring low-dimensional microstructure representations using convolutional neural networks

Authors: Nicholas Lubbers, Turab Lookman, Kipton Barros

Abstract: We apply recent advances in machine learning and computer vision to a central problem in materials informatics: The statistical representation of microstructural images. We use activations in a pre-trained convolutional neural network to provide a high-dimensional characterization of a set of synthetic microstructural images. Next, we use manifold learning to obtain a low-dimensional embedding of… ▽ More We apply recent advances in machine learning and computer vision to a central problem in materials informatics: The statistical representation of microstructural images. We use activations in a pre-trained convolutional neural network to provide a high-dimensional characterization of a set of synthetic microstructural images. Next, we use manifold learning to obtain a low-dimensional embedding of this statistical characterization. We show that the low-dimensional embedding extracts the parameters used to generate the images. According to a variety of metrics, the convolutional neural network method yields dramatically better embeddings than the analogous method derived from two-point correlations alone. △ Less

Submitted 30 November, 2018; v1 submitted 8 November, 2016; originally announced November 2016.

Comments: 14 Pages, 12 Figures

Journal ref: Phys. Rev. E 96, 052111 (2017)

Shape of magnetic domain walls formed by coupling to mobile charges

Authors: Ryo Ozawa, Satoru Hayami, Kipton Barros, Yukitoshi Motome

Abstract: Magnetic domain walls, which are crucially important in both fundamental physics and technical applications, often have a preference in their form due to many different origins, such as the crystalline shape, lattice symmetry, and magnetic anisotropy. We theoretically investigate yet another origin stemming from the coupling to mobile charges in itinerant magnets. Performing a large-scale numerica… ▽ More Magnetic domain walls, which are crucially important in both fundamental physics and technical applications, often have a preference in their form due to many different origins, such as the crystalline shape, lattice symmetry, and magnetic anisotropy. We theoretically investigate yet another origin stemming from the coupling to mobile charges in itinerant magnets. Performing a large-scale numerical simulation in a minimal model for itinerant magnets, i.e., the Kondo lattice model with classical localized spins, we show that the shape of magnetic domain walls depends on the electronic band structure and electron filling. While Neel and 120 antiferromagnetic states do not show a strong preference in the shape of domain walls, noncoplanar spin states with scalar chiral ordering have distinct directional preferences of the domain walls depending on the electron filling. We find that the directional preference is rationalized by the wave-number dependence of the effective magnetic interactions induced by the mobile charges, which are set by the band structure and electron filling. We also observe that, in the noncoplanar chiral states, an electric current is induced along the domain walls owing to the spin Berry phase mechanism, with very different spatial distributions depending on whether the bulk state is metallic or insulating. △ Less

Submitted 22 September, 2016; originally announced September 2016.

Comments: 12 pages, 4 figures

Journal ref: Phys. Rev. B 96, 094417 (2017)

Searching for cosmological signatures of the Einstein equivalence principle breaking

Authors: R. F. L. Holanda, K. N. N. O. Barros

Abstract: Modifications of gravity generated by a multiplicative coupling of a scalar field to the electromagnetic Lagrangian lead to a breaking of Einstein equivalence principle (EEPB) as well as to variations of fundamental constants. In these theoretical frameworks, deviations of standard values of the fine structure constant, $Δでるたαあるふぁ/αあるふぁ=φふぁい$, and of the cosmic distance duality relation,… ▽ More Modifications of gravity generated by a multiplicative coupling of a scalar field to the electromagnetic Lagrangian lead to a breaking of Einstein equivalence principle (EEPB) as well as to variations of fundamental constants. In these theoretical frameworks, deviations of standard values of the fine structure constant, $Δでるたαあるふぁ/αあるふぁ=φふぁい$, and of the cosmic distance duality relation, $D_L(1+z)^{-2}/D_A=ηいーた=1$, where $D_L$ and $D_A$ are the luminosity and angular diameter distances, respectively, are unequivocally linked. In this paper, we search for cosmological signatures of the EEPB by using angular diameter distance from galaxy clusters, obtained via their Sunyaev-Zeldovich effect (SZE) and X-ray observations, and distance modulus of type Ia supernovae (SNe Ia). The crucial point here is that we take into account the dependence of the SZE/X-ray technique with $φふぁい$ and $ηいーた$. Our new results show no indication of the EEPB. △ Less

Submitted 5 July, 2016; v1 submitted 25 June, 2016; originally announced June 2016.

Comments: 2 figures, 5pages, To appear in Physical Review D

Journal ref: Phys. Rev. D 94, 023524 (2016)