Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FACTS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Nematic superconductivity in twisted bilayer graphene

Twisted bilayer graphene displays insulating and superconducting phases caused by exceptional flattening of its lowest energy bands. Superconductivity with highest T c appears near half filling of the valence band (n ≈ –2). The data show that in the large part of the superconducting dome near n = –2, threefold lattice rotation symmetry is broken in the superconducting phase, i.e., a superconductor is also a nematic. We perform a comprehensive analysis of superconductivity in twisted-bilayer graphene within an itinerant approach and present a mechanism for nematic superconductivity. We take as an input the fact that at dopings, where superconductivity has been observed, the Fermi energy lies in the vicinity of twist-induced Van Hove singularities in the density of states. We argue that the low-energy physics can be properly described by patch models with either six or twelve Van Hove points. We obtain pairing interactions for the patch models in terms of parameters of the microscopic model for the flat bands, which contains both local and twist-induced nonlocal interactions. Here, we show that the latter give rise to attraction in different superconducting channels. For six Van Hove points, there is just one attractive d-wave channel, and we find chiral d ± id superconducting order, which breaks time-reversal symmetry but leaves the lattice rotation symmetry intact. For twelve Van Hove points, we find two attractive channels, g and i waves, with almost equal coupling constants. We show that both order parameters are nonzero in the ground state and explicitly demonstrate that in this co-existence state the threefold lattice rotation symmetry is broken, i.e., the superconducting state is also a nematic. We find two possible nematic states, one is time-reversal symmetric, the other additionally breaks time-reversal symmetry. We apply our scenario to twisted bilayer graphene near n = –2 and argue that it is applicable also to other systems with two (or more) attractive channels with similar couplings as our reasoning for a nematic superconductivity is based on generic symmetry considerations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Trivial to nontrivial topology transition in rare-earth pnictides with epitaxial strain

The combination of magnetotransport and topological properties has brought great attention to rare-earth monopnictides semimetals. For some of them, like LaSb, it is unclear whether they show nontrivial topology or not based on density functional theory calculations and angular resolved photoemission spectroscopy measurements. Here we use hybrid density functional theory to demonstrate that LaSb is in fact a trivial topological semimetal, in agreement with experiments, but on the verge of a transition to a topological phase. In this work, we show that under compressive epitaxial strain, the La $\textit{d}$ band crosses the Sb $\textit{p}$ band near the $\textit{Z}$ point in the Brillouin zone, stabilizing a topologically nontrivial phase, opening unique opportunities to probe the inter-relation between magnetotransport properties and the effects of band topology by examining epitaxially strained and unstrained thin films of the same material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Symmetry-breaking polymorphous descriptions for correlated materials without interelectronic U

Correlated materials with open-shell d- and f-ions having degenerate band edge states show a rich variety of interesting properties ranging from metal-insulator transition to unconventional superconductivity. The textbook view for the electronic structure of these materials is that mean-field approaches are inappropriate, as the interelectronic interaction U is required to open a band gap between the occupied and unoccupied degenerate states while retaining symmetry. We show that the latter scenario often defining what Mott insulators are, is in fact not needed for the 3d binary oxides MnO, FeO, CoO, and NiO. The mean-field-like band theory can indeed lift such degeneracies in the binaries when nontrivial unit cell representations (polymorphous networks) are allowed to break symmetries, in conjunction with a recently developed non-empirical exchange and correlation density-functional without an on-site interelectronic interaction U. We explain how density-functional theory (DFT) in the polymorphous representation achieves band gap opening in correlated materials through a separate mechanism to the Mott-Hubbard approach. Here, we show the method predicts magnetic moments and gaps for the four binary monoxides in both the antiferromagnetic and paramagnetic phases, offering an effective alternative to symmetry-conserving approaches for studying a range of functionalities in open d- and f-shell complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

fcc lattice, checkerboards, fractons, and quantum field theory

In this work, we consider XY-spin degrees of freedom on an FCC lattice, such that the system respects some subsystem global symmetry. We then gauge this global symmetry and study the corresponding $U(1)$ gauge theory on the FCC lattice. Surprisingly, this $U(1)$ gauge theory is dual to the original spin system. We also analyze a similar $\mathbb{Z}_N$ gauge theory on that lattice. All these systems are fractonic. The $U(1)$ theories are gapless and the $\mathbb{Z}_N$ theories are gapped. We analyze the continuum limits of all these systems and present free continuum Lagrangians for their low-energy physics. Our $\mathbb{Z}_2$ FCC gauge theory is the continuum limit of the well known checkerboard model of fractons. Our continuum analysis leads to a straightforward proof of the known fact that this theory is dual to two copies of the $\mathbb{Z}_2$ X-cube model. We find new models and new relations between known models. The $\mathbb{Z}_N$ FCC gauge theory can be realized by coupling three copies of an anisotropic model of lineons and planons to a certain exotic $\mathbb{Z}_2$ gauge theory. Also, although for $N=2$ this model is dual to two copies of the $\mathbb{Z}_2$ X-cube model, a similar statement is not true for higher $N$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Competing energetic states in γ – Fe 2 W O 6 with strong spin-charge-lattice coupling

We report magnetic and electronic properties of γ – Fe 2 W O 6 via neutron powder-diffraction measurements and first-principles density-function theory calculations. We reveal a magnetic ground state with two phases coexisting in which the minor phase is the same as an earlier report, despite the fact that both materials studied have the same space group but with slightly different atomic positions. Interestingly, both magnetic phases are well captured by first-principles calculations. Furthermore, we show that these two magnetic phases are correlated with electronic properties, with one being insulating and the other being metallic. These features suggest that γ – Fe 2 W O 6 exhibits competing energetic states in which spin, charge, and lattice degrees of freedom are strongly coupled to each other.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

First-principles modeling of conductivity at the (001), (110), and (111) SrTiO 3 /LaAlO 3 heterointerfaces

Here, the complex polar oxide heterojunction of SrTiO 3 /LaAlO 3 (STO/LAO) is of great interest due to the emergent physical phenomena observed at the interface. STO and LAO separately are wide band-gap insulators. However, upon joining them at the 001, 110, and 111 crystallographic planes, the interface undergoes a transition to a conductive state. Although first-principles modeling of the 001 plane interface has been widely studied, there is a lack of reports regarding the 110 and 111. This paper expands the theoretical model of the STO/LAO heterointerface to the three crystallographic planes (001, 110, and 111) where the conductivity has been experimentally reported. The calculations showed that whereas at the 001 interface the conductivity appears at a critical thickness of 4 monolayers of LAO, the 110 and 111 planes have no clear critical thickness; these two interfaces were always conductive. Nevertheless, the number of conductive electrons per unit cell increases with the thickness of the LAO layer in the 110 and 111 interfaces. This is related to the energy levels downshifting due to the electrostatic potential buildup (which was in the opposite direction respect to the 001 interface), increasing the number of conductive sates below the Fermi level. Given the absence of a critical thickness and the fact that chemical intermixing and oxygen vacancies at the interface were not considered, the main mechanism responsible for the conductivity in the 110 and 111 planes was attributed to the large structure reconstruction that locally changes the energy levels at the interface causing charge transfer and accumulation at the layers close to the interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Full versus quasiparticle self-consistency in vertex-corrected GW approaches

Here, using seven semiconductors/insulators with band gaps covering the range from 1 eV to 10 eV we systematically explore the performance of two different variants of self-consistency associated with famous Hedin's system of equations: the full self-consistency and the so called quasi-particle approximation to it. The pros and cons of these two variants of self-consistency are sufficiently well documented in literature for the simplest GW approximation to the Hedin's equations. Our study, therefore, aims primarily at the level of theory beyond GW approximation, i.e. at the level of theory which includes vertex corrections. Whereas quasi-particle self-consistency has certain advantages at GW level (well known fact), the situation becomes quite different when vertex corrections are included. In the variant with full self-consistency, vertex corrections (both for polarizability and for self energy) systematically reduce the calculated band gaps making them closer to the experimental values. In the variant with quasi-particle self-consistency, however, an inclusion of the same diagrams has considerably larger effect and calculated band gaps become severely underestimated. Different effect of vertex corrections in two variants of self-consistency can be related to the Z-factor cancellation which plays positive role in quasi-particle self-consistency at GW level of theory but appears to be destructive for the quasi-particle approximation when higher order diagrams are included. Second result of our study is that we were able to reproduce the results obtained with the Questaal code using our FlapwMBPT code when the same variant of self-consistency (quasi-particle) and the same level of vertex corrections (for polarizability only, static approximation for screened interaction, and Tamm-Dancoff approximation for the Bethe-Salpeter equation) are used.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Complex structure due to As bonding and interplay with electronic structure in superconducting BaNi 2 As 2

BaNi 2 As 2 is a superconductor chemically related to the Fe-based superconductors, with a complex and poorly understood structural phase transition. We show, based on first-principles calculations, that, in fact, there are two distinct competing structures. These structures are different from the electronic, transport, and bonding points of view but are close in energy. These arise due to complex As bonding patterns and drive distortions of the Ni layers. This is supported by photoemission experiments. This is very distinct from views of the distortion as being primarily driven by electrons at the Fermi surface as in a classical charge density wave, or by correlated electron physics associated with the Ni d electrons and the Ni d orbitals. The structural distortion, although connected with As does lead to an interplay of electronic and structural behavior including induced anisotropic electronic transport. Here, the local bonding nature of the instabilities and the competition between distortions is discussed in the context of the complex behavior observed in BaNi 2 As 2 samples.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Possible quadrupole-order-driven commensurate-incommensurate phase transition in B20 CoGe

For this study, the B20-type cobalt germanide CoGe was investigated by measuring the specific heat, resistivity, and 59 Co nuclear magnetic resonance (NMR). We observed a phase transition at T Q = 13.7 K, evidenced by a very narrow peak of the specific heat and sharp changes of the nuclear spin-spin (T 2 –1 ) and spin-lattice ( T 1 – 1 ) relaxation rates. The fact that the entropy release is extremely small and the Knight shift is almost independent of temperature down to low temperatures as anticipated in a paramagnetic metal indicates that the T Q transition is of nonmagnetic origin. In addition, we detected a crossover scale T 0 ~ 30 K below which the resistivity and the NMR linewidth increase, and T 1 –1 is progressively distributed in space, that is, a static and dynamical spatial inhomogeneity develops. While the order parameter for the T Q transition remains an open question, a group-theoretical analysis suggests that the finite electric quadrupole density arising from the low local site symmetry at cobalt sites could drive the crystal symmetry lowering from the P2 1 3 symmetry that is commensurate to the R3 symmetry with an incommensurate wave vector, which fairly well accounts for the T Q transition. The quadrupole-order-driven commensurate-incommensurate phase transition may be another remarkable phenomenon arising from the structural chirality inherent in the noncentrosymmetric B20 family.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Phonon-induced rotation of the electronic nematic director in superconducting Bi 2 Se 3

The doped topological insulator A x Bi 2 Se 3 , with A = {Cu,Sr,Nb}, becomes a nematic superconductor below T c ≈ 3–4 K. The associated electronic nematic director is described by an angle α and is experimentally manifested in the elliptical shape of the in-plane critical magnetic field H c2 . Because of the threefold rotational symmetry of the lattice, α is expected to align with one of three high-symmetry directions corresponding to the in-plane nearest-neighbor bonds, consistent with a Z 3 -Potts nematic transition. Here, we show that the nematic coupling to the acoustic phonons, which makes the nematic correlation length tend to diverge along certain directions only, can fundamentally alter this phenomenology in trigonal lattices. Compared to hexagonal lattices, the former possesses a sixth independent elastic constant c14 due to the fact that the in-plane shear strain doublet (ε xx – ε yy , – 2ε xy ) and the out-of-plane shear strain doublet (2ε yz , – 2ε zx ) transform as the same irreducible representation. We find that, when c 14 overcomes a threshold value, which is expected to be the case in doped Bi 2 Se 3 , the nematic director α unlocks from the high-symmetry directions due to the competition between the quadratic phonon-mediated interaction and the cubic nematic anharmonicity. This implies the breaking of the residual in-plane twofold rotational symmetry (C 2x ), resulting in a triclinic phase. Furthermore, we discuss the implications of these findings for the structure of nematic domains, for the shape of the in-plane H c2 in A x Bi 2 Se 3 , and for the presence of nodes inside the superconducting state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Dynamical exponent of a quantum critical itinerant ferromagnet: A Monte Carlo study

Here we consider the effect of the coupling between two-dimensional (2D) quantum rotors near an XY ferromagnetic quantum critical point and spins of itinerant fermions. We analyze how this coupling affects the dynamics of rotors and the self-energy of fermions. A common belief is that near a q = 0 ferromagnetic transition, fermions induce an Ω/q Landau damping of rotors (i.e., the dynamical critical exponent is z = 3) and Landau overdamped rotors give rise to non-Fermi liquid fermionic self-energy Σ∝ω 2/3 . This behavior has been confirmed in previous quantum Monte Carlo (QMC) studies. Here we show that for the XY case the behavior is different. We report the results of large-scale quantum Monte Carlo simulations, which show that at small frequencies z = 2 and Σ∝ω 1/2 . We argue that the new behavior is associated with the fact that a fermionic spin is by itself not a conserved quantity due to spin-spin coupling to rotors, and a combination of self-energy and vertex corrections replaces 1/q in the Landau damping by a constant. We discuss the implication of these results to experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Strong increase in ultrasound attenuation below T c in Sr 2 RuO 4 : Possible evidence for domains

Recent experiments suggest that Sr 2 RuO 4 has a two-component superconducting order parameter (OP). A two-component OP has multiple degrees of freedom in the superconducting state that can result in low-energy collective modes or the formation of domain walls—a possibility that would explain a number of experimental observations including the smallness of the signature of time reversal symmetry breaking at T c and telegraph noise in critical current experiments. We use resonant ultrasound spectroscopy to perform ultrasound attenuation measurements across the superconducting T c of Sr 2 RuO 4 . We find that compressional sound attenuation increases by a factor of 7 immediately below T c , in sharp contrast with what is found in both conventional (s-wave) and high-T c (d-wave) superconductors. Our observations are most consistent with the presence of domain walls that separate different configurations of the superconducting OP. Here, the fact that we only observe an increase in sound attenuation for compressional strains, and not for shear strains, suggests an inhomogeneous superconducting state formed of two distinct, accidentally degenerate superconducting OPs that are not related to each other by symmetry. Whatever the mechanism, a factor of 7 increase in sound attenuation is a singular characteristic that must be reconciled with any potential theory of superconductivity in Sr 2 RuO 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Critical fields of superconductors with magnetic impurities

The upper critical field H c2 , the field H c3 for nucleation of the surface superconductivity, and the thermodynamic field H c are evaluated within the weak-coupling theory for the isotropic s-wave case with arbitrary transport, and pair-breaking scattering. In this study, we find that, for the standard geometry of a half-space sample in a magnetic field parallel to the surface, the ratio $\mathscr{R}$=H c3 /H c2 is within the window 1.55≲$\mathscr{R}$≲2.34 , regardless of temperature or the scattering type. While the nonmagnetic impurities tend to flatten the $\mathscr{R}$($\mathrm{T}$) variation, magnetic scattering merely shifts the maximum of $\mathscr{R}$($\mathrm{T}$) to lower temperatures. Surprisingly, while reducing the transition temperature, magnetic scattering has a milder impact on $\mathscr{R}$ than nonmagnetic scattering. The surface superconductivity is quite robust; in fact, the ratio $\mathscr{R}$≈1.7 even in the gapless state. We used Eilenberger's energy functional to evaluate the condensation energy F c and the thermodynamic critical field H c for any temperature and scattering parameters. By comparing H c2 and H c , we find that, unlike transport scattering, the pair-breaking pushes materials toward type-I behavior. Furthermore, we find a peculiar behavior of F c as a function of the pair-breaking scattering parameter at the low-$\mathrm{T}$ transition from gapped to gapless phases, which has recently been associated with the topological transition in the superconducting density of states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Suppression of ferromagnetism governed by a critical lattice parameter in CeTiGe 3 with hydrostatic pressure or V substitution

We combine structural and magnetic measurements to compare the different magnetic phase diagrams between the pressure and substitution studies in CeTiGe 3 . We report on the structural, magnetic, and electrical transport properties of single crystals of CeTi 1–x V x Ge 3 (x = 0, 0.1, 0.2, 0.3, 0.4, 0.9, and 1), and of polycrystalline samples (x = 0.5, 0.6, 0.7, 0.8), as well as structural properties of CeTiGe 3 under pressure up to 9 GPa. The ferromagnetic ordering in CeTiGe 3 is suppressed with V doping in CeTi 1–x V x Ge 3 , and suggests a possible ferromagnetic quantum critical point near x = 0.45. We perform a detailed crystalline electric field (CEF) analysis, and the magnetic susceptibility data in pure CeTiGe 3 and CeVGe 3 can be well explained by the CEF model. The proposed CEF energy levels suggest that there is a gradual change of the ground state from |±5/2 in CeTiGe 3 to |±1/2 in CeVGe 3 , and a suppression of CEF splitting energies near the quantum critical region. When hydrostatic pressure is used instead of chemical substitution, the quantum critical point is avoided by the appearance of magnetic phases above around 4.1 GPa. In the substitution study, the ferromagnetic and antiferromagnetic regions are well separated, whereas they touch in the pressure study. We observe a different trend in the temperature dependence of the resistivity maximum in both studies, suggesting that the CEF splitting energy is suppressed by V substitution but enhanced by pressure. Here, we also observe different responses in lattice constants between the two studies, highlighting the fact that substitution effects cannot be reduced to chemical pressure effects only. Nevertheless, when the magnetic phase diagrams of both hydrostatic pressure and substitution are compared, we find a common critical lattice constant c = 5.78 Å where the ferromagnetic ordering is suppressed in both studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fractionalization of subsystem symmetries in two dimensions

The fractionalization of global symmetry charges is a striking hallmark of topological quantum order. Here, we discuss the fractionalization of subsystem symmetries in two-dimensional topological phases. In line with previous no-go arguments, we show that subsystem symmetry fractionalization is not possible in many cases due to the additional rigid geometric structure of the symmetries. However, we identify a mechanism that allows fractionalization, involving global relations between macroscopically many symmetry generators. We find that anyons can fractionalize such relations, meaning that the total charge carried under all generators involved in the global relation is nontrivial, despite the fact that these generators multiply to the identity. We first discuss the general algebraic framework needed to characterize this type of fractionalization, and then explore this framework using a number of exactly solvable models with Z 2 topological order, including models having line and fractal symmetries. These models all showcase another necessary property of subsystem symmetry fractionalization: Fractionalized anyons must have restricted mobility when the symmetry is enforced, such that they are confined to a single line or point in the case of line and fractal symmetries, respectively. Looking forward, we expect that our identification of the importance of global relations in fractionalization will hold significance for the classification of phases with subsystem symmetries in all dimensions.

2-dimensional systems↗
DOE OSTI2022

Role of electromagnetic gauge-field fluctuations in the selection between chiral and nematic superconductivity

Motivated by the observation of nematic superconductivity in several systems, we revisit the problem of the leading pairing instability of two-component unconventional superconductors on the triangular lattice—such as (p x , p y )-wave and (d x 2 –y 2 , d xy )-wave superconductors. Such a system has two possible superconducting states: the chiral state (e.g., p + ip or d + id), which breaks time-reversal symmetry, and the nematic state (e.g., p + p or d + d), which breaks the threefold rotational symmetry of the lattice. Weak-coupling calculations generally favor the chiral over the nematic superconducting state, raising the question of what mechanism can stabilize the latter. Here, we show that the electromagnetic field fluctuations can play a crucial role in selecting between these two states. Specifically, we derive and analyze the effective free energy for the two-component superconducting order parameter after integrating out the gauge-field fluctuations, which is formally justified if the spatial order parameter fluctuations can be neglected. A nonanalytic cubic term arises, as in the case of a conventional s -wave superconductor. However, unlike the latter, the cubic term depends on the relative phase and on the relative amplitudes between the two order parameter components, in such a way that it generally favors the nematic state. This result is a direct consequence of the fact that the stiffness of the superconducting order parameter is not isotropic. Competition with the quartic term, which favors the chiral state, leads to a renormalized phase diagram in which the nematic state displaces the chiral state over a wide region in the parameter space. Here, we analyze the stability of the fluctuation-induced nematic phase, generalize our results to tetragonal lattices, and discuss their applicability to candidate nematic superconductors, including twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Lefschetz thimble quantum Monte Carlo for spin systems

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Furthermore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗