Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spin wave theory”

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 325 records · Page 18

Precursor of pair-density wave in doping Kitaev spin liquid on the honeycomb lattice

Abstract We study the effects of doping the Kitaev model on the honeycomb lattice where the spins interact via the bond-directional interaction J K , which is known to have a quantum spin liquid as its exact ground state. The effect of hole doping is studied within the t - J K model on a three-leg cylinder using density-matrix renormalization group. Upon light doping, we find that the ground state of the system has a dominant quasi-long-range charge-density-wave correlations but short-range single-particle correlations. In the pairing channel, the even-parity superconducting correlation is dominant with d -wave-like symmetry, which oscillates in sign as a function of separation with a period equal to that of the spin-density wave and two times the charge-density wave. Although these correlations fall rapidly (possibly exponentially) at long distances, this is never-the-less the example where a pair-density wave is the leading instability in the pairing channel on the honeycomb lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets

Natural orbitals are often used to achieve a more compact representation of correlated wave-functions. Using natural orbitals computed as eigenstates of the virtual–virtual block of the state density matrix instead of the canonical Hartree–Fock orbitals results in smaller errors when the same fraction of virtual space is frozen. This strategy, termed frozen natural orbital (FNO) approach, is effective in reducing the cost of regular coupled-cluster (CC) calculations and some multistate methods, such as EOM-IP-CC (equation-of-motion CC for ionization potentials). This contribution extends the FNO approach to the EOM-SF-CC ansatz (EOM-CC with spin-flip). In contrast to EOM-IP-CCSD, EOM-SF-CCSD relies on high-spin open-shell references. Using FNOs computed for an open-shell reference leads to an erratic behavior of the EOM-SF-CC energies and properties due to an inconsistent truncation of the α and β orbital spaces. A general solution to problems arising in the EOM-CC calculations utilizing open-shell references, termed OSFNO (open-shell FNO), is proposed. Furthermore, by means of singular value decomposition (SVD) of the overlap matrix of the α and β orbitals, the OSFNO algorithm identifies the corresponding orbitals and determines virtual orbitals corresponding to the singly occupied space. This is followed by SVD of the singlet part of the state density matrix in the remaining virtual orbital subspace. The so-computed FNOs preserve the spin purity of the open-shell orbital subspace to the extent allowed by the original reference, thus facilitating a safe truncation of the virtual space. The performance of OSFNO is benchmarked for selected diradicals and triradicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effective Hamiltonians derived from equation-of-motion coupled-cluster wave functions: Theory and application to the Hubbard and Heisenberg Hamiltonians

Effective Hamiltonians, which are commonly used for fitting experimental observables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not trivial. In this contribution, we apply Bloch’s formalism to equation-of-motion coupled-cluster (EOM-CC) wave functions to rigorously derive effective Hamiltonians in the Bloch’s and des Cloizeaux’ forms. We report the key equations and illustrate the theory by application to systems with two or three unpaired electrons, which give rise to electronic states of covalent and ionic characters. We show that the Hubbard and Heisenberg Hamiltonians can be extracted directly from the so-obtained effective Hamiltonians. By establishing a quantitative connection between many-body states and simple models, the approach facilitates the analysis of the correlated wave functions. We propose a simple diagnostic for assessing the validity of the model space choice based on the overlaps between the target- and model-space states. Finally, artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lattice collective modes from a continuum model of magic-angle twisted bilayer graphene

Here, we show that the insulating states of magic-angle twisted bilayer graphene support a series of collective modes corresponding to local particle-hole excitations on triangular lattice sites. Our theory is based on a continuum model of the magic angle flat bands. When the system is insulating at moiré band filling ν = -3, our calculations show that the ground state supports seven low-energy modes that lie well below the charge gap throughout the moiré Brillouin zone, one of which couples strongly to THz photons. The low-energy collective modes are faithfully described by a model with a local SU(8) degree of freedom in each moiré unit cell that we identify as the direct product of spin, valley, and an orbital pseudospin. Apart from spin and valley-wave modes, the collective mode spectrum includes a low-energy intraflavor exciton mode associated with transitions between flat valence and conduction band orbitals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structural Studies of Alloyed and Nanoparticle Transition Metal Dichalcogenides by Selenium-77 Solid-State Nuclear Magnetic Resonance Spectroscopy

Layered transition metal dichalcogenides (TMDCs) such as MoS 2 , MoSe 2 , and WSe 2 are under intense investigation because they are atomically thin semiconductors with photophysical properties that can be tuned by changing their composition or morphology. Mechanochemical processing has been proposed as a method to obtain alloyed TMDCs in the series Mo 1–x W x S y Se 2–y (x = 0–1; y = 0–2). However, elucidating the chemical transformations occurring at the atomic scale following mechanochemical processing can be challenging because the products are often amorphous or microcrystalline. To address this challenge, we probe TMDC mixing and alloying by using a combination of powder X-ray diffraction, Raman spectroscopy, diffuse reflectance spectroscopy, 77 Se solid-state nuclear magnetic resonance (SSNMR) spectroscopy, and planewave density functional theory (DFT) calculations. The nature of the milling material and reaction atmosphere are shown to be essential factors in limiting the formation of undesired oxide byproducts. We demonstrate acquisition of 77 Se SSNMR spectra using different combinations of Carr-Purcell Meiboom-Gill acquisition (CPMG) pulse sequences, magic angle spinning (MAS), and MAS dynamic nuclear polarization. Further, the combination of SSNMR with the other characterization methods clearly demonstrates that high energy impact ball milling induces molecular level alloying of Mo, W and chalcogen atoms in the family Mo 1–x W x S y Se 2–y . Gauge including projector augmented wave DFT calculations yield accurate 77 Se chemical shift (CS) tensor components. 77 Se SSNMR spectroscopy was also applied to study the structure of WSe 2 nanocrystals intercalated with ethylenediamine. The intercalated WSe 2 nanocrystals exhibit a more positive isotropic 77 Se CS as compared to bulk WSe 2 , however, the 77 Se CS anisotropy is the same, confirming the WSe 2 layers have a similar structure as in their bulk counterparts.

36 MATERIALS SCIENCE↗

Measures of complexity and entanglement in many-fermion systems

There is no unique and widely accepted definition of the complexity measure (CM) of a many-fermion wave function in the presence of interactions. The simplest many-fermion wave function is a Slater determinant. In shell-model or configuration interaction (CI) and other related methods, the state is represented as a superposition of a large number of Slater determinants, which in the case of CI calculations reaches about 20 billion terms [Johnson, arXiv:1809.07869]. Although in practice this number has been used as a CM for decades, it is ill defined: it is not unique, and it depends on the particular type and the number of single-particle wave functions used to construct the Slater determinants. Further, the canonical wave functions and/or natural orbitals [Löwdin, Adv. Phys. 5, 1 (1956); Löwdin and Shull, Phys. Rev. 101, 1730 (1956); Bardeen et al., Phys. Rev. 108, 1175 (1957); N. N. Bogoljubov, Il Nuovo Cimento 7, 794 (1958); Valatin, Il Nuovo Cimento 7, 843 (1958); de Gennes, Superconductivity of Metals and Alloys (CRC Press, Boca Raton, FL, 1999); Ring and Schuck, The Nuclear Many-Body Problem, 1st ed. (Springer-Verlag, Berlin, 2004)] and their corresponding occupation probabilities are intrinsic properties of any many-body wave function, irrespective of the representation, and they provide a unique solution to characterize the CM. The non-negative orbital entanglement entropy, which vanishes for a Slater determinant, provides the simplest CM, while a more complete measure of complexity is the entanglement spectrum. We illustrate these aspects in the case of a complex nonequilibrium time-dependent process, induced nuclear fission described within a real-time density functional theory framework extended to superfluid systems, which can describe simultaneously the long-range and the short-range correlations between fermions. The orbital entanglement entropy of the fissioning nucleus illustrates the localization mechanism of the many-body wave function in Fock and/or Hilbert space. The (minimal) number of Slater determinants required to represent such a complex many-body wave function with a well-defined number of particles in the case presented here is about 10 500 . The realistic case of the highly nonequilibrium nuclear fission process illustrated here is equivalent to a system of 23.328×10 9 interacting quantum spin-1/2 particles, a very large system for the study of quantum entanglement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calculated electron paramagnetic resonance g-tensor and hyperfine parameters for zinc vacancy and N related defects in ZnO

Various defects in ZnO, focused on substitutional N O and N 2 in various sites, O-site, interstitial and Zn-site are studied using first-principles calculations with the goal of understanding the electron paramagnetic resonance (EPR) center reported for N 2 in ZnO and substitutional N on the O-site. The g tensors are calculated using the gauge including projector augmented wave (GIPAW) method and compared with experiments. The g-tensor of the free $N^{+}_{2}$ and $N^{–}_{2}$ radicals and their various contributions within the GIPAW theory are analyzed first to provide a baseline reference for the accuracy of the method and for understanding the N 2 behavior in ZnO. Previous controversies on the site location of N 2 in ZnO for this EPR center and on the shallow or deep nature and donor or acceptor nature of this center are resolved. Here, we find that the N 2 on the Zn site is mostly zinc-vacancy like in its spin density and g-tensor, while for the O-site, a model with the N 2 axis lying in-the basal plane and the singly occupied π g -orbital along the c axis provides good agreement with experiment. For the interstitial location, if the N 2 is not strongly interacting with the surroundings, no levels in the gap are found and hence also no possible EPR center. The calculated g-tensors for N O and V Zn are also found to be in good agreement with experiment. The effects of different functionals affecting the localization of the spin density are shown to affect the g-tensor values.

36 MATERIALS SCIENCE↗

Ab initio calculations for void swelling bias in α - and δ -plutonium

Void swelling can develop in materials under persistent irradiation when nonequilibrium vacancy and self-interstitial populations migrate under sufficiently asymmetric interaction biases. In conventional metals, the propensity is determined to the first approximation by comparing point-defect relaxation strains. We thus present DFT-based calculations of structures and formation energies and volumes of point defects in the α and the δ phases of plutonium. Here we discuss the pros and cons of various levels of electronic structure theory: spin polarization, spin-orbit coupling, and orbital polarization. Our results show that lattice defects in δ-Pu, in contrast to most fcc metals, have surprisingly small formation volumes. Equally unexpected are the large defect formation volumes found in the low-symmetry α-Pu phase. Both these unusual properties can be satisfactorily explained from defect-induced spin/orbital moment formation and destruction in the Pu phases. Surprisingly, the point defects in α-Pu are found to induce far larger transformation of the local electronic structure than in δ-Pu. When we use the calculated defect properties to estimate the classic void swelling bias in each of the phases, we find it to be unusually small in δ-Pu but likely much larger in α-Pu. Hence, swelling rates and mechanisms can diverge dramatically between the different phases of Pu. Especially in the transient regime before the formation of large defect clusters, the swelling rate of α-Pu can reliably be expected to be much larger than δ-Pu. However, accurate forecasts over longer times will require the conventional void-swelling theory to be modified to handle the complexities presented by the different Pu phases. As a case in point, we show the possible anomalous temperature dependence of vacancy properties in δ-Pu, caused by entropic contributions from defect-induced spin-lattice fluctuations. Such complications may affect defect-defect interactions and thus alter the void swelling bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Partial densities of states from x-ray absorption and Auger electron spectroscopy: Use of core-hole memory

Ag L 3 edge x-ray absorption and L 3 −M 4,5 M 4,5 Auger electron spectra have been measured and simulated for a variety of x-ray electric-field polarization and Auger electron emission directions. The theory relies on density-functional theory, use of the Bethe–Salpeter equation, atomic multiplet theory, and a simplified model for the Auger line shape and its dependence on photon energy. We also demonstrate that partial densities of states for dT 2g , dE g , and sA g symmetry partial-wave channels at the Ag site in the solid can be deduced from the experimental measurements with only atomic theoretical input, i.e., with no solid-state theory other than assuming a Fermi edge.

36 MATERIALS SCIENCE↗

Fluctuating intertwined stripes in the strange metal regime of the Hubbard model

Strongly correlated electron systems host a variety of poorly understood correlations in their high-temperature normal state. Unlike ordered phases defined by order parameters, regions of the normal state are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground-state orders. Here, in this paper, we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short ranged; yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Helical magnetic order and Fermi surface nesting in noncentrosymmetric ScFeGe

Here, an investigation of the structural, magnetic, thermodynamic, and charge transport properties of noncentrosymmetric hexagonal ScFeGe reveals it to be an anisotropic metal with a transition to a weak itinerant incommensurate helimagnetic state below T N =36 K. Neutron diffraction measurements discovered a temperature and field independent helical wave vector k = (0 0 0.193) with magnetic moments of 0.53 μ B per Fe confined to the ab-plane. Density functional theory calculations are consistent with these measurements and find several bands that cross the Fermi level along the c-axis with a nearly degenerate set of flat bands just above the Fermi energy. The anisotropy found in the electrical transport is reflected in the calculated Fermi surface, which consists of several warped flat sheets along the c-axis with two regions of significant nesting, one of which has a wave vector that closely matches that found in the neutron diffraction. The electronic structure calculations, along with a strong anomaly in the c-axis conductivity at TN, signal a Fermi surface driven magnetic transition, similar to that found in spin density wave materials. Magnetic fields applied in the ab-plane result in a metamagnetic transition with a threshold field of ≈ 6.7 T along with a sharp, strongly temperature dependent discontinuity and a change in sign of the magnetoresistance for in-plane currents. Thus, ScFeGe is an ideal system to investigate the effect of in-plane magnetic fields on a helimagnet with a c-axis propagation vector, where the relative strength of the magnetic interactions and anisotropies determine the topology and magnetic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Calculating microscopic optical potentials for elastic scattering at intermediate energies from light nuclei in an ab initio fashion within the Watson expansion has been established within the last few years. Based on the Watson expansion of the multiple scattering series, we employ a nonlocal translationally invariant nuclear density derived within the symmetry-adapted no-core shell model (SA-NCSM) framework from a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction and the very same interaction for a consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for medium-heavy nuclei. The leading order effective (optical) folding potential is computed by integrating over a translationally invariant SA-NCSM one-body scalar density, spin-projected momentum distribution, and the Wolfenstein amplitudes 𝐴, 𝐶, and 𝑀. The resulting nonlocal potentials serve as input for a momentum space Lippmann-Schwinger equation. In the SA-NCSM, the model space is systematically up-selected using Sp⁡(3,ℝ) symmetry considerations. For the light nucleus of 6 He, we establish a systematic selection scheme in the SA-NCSM for scattering observables. Then, we apply this scheme to calculations of scattering observables, such as differential cross sections, analyzing powers, and spin rotation functions for elastic proton scattering from 20 Ne and 40 Ca in the energy regime between 65 and 200 MeV, and compare to available data. Furthermore, our calculations show that the leading order effective nucleon-nucleus potential in the Watson expansion of multiple scattering theory obtained from an up-selected SA-NCSM model space describes 40 Ca elastic scattering observables reasonably well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the nucleon-nucleon amplitudes and the one-body scalar and spin nuclear densities.

Ab initio calculations↗

Thermodynamic evidence for a two-component superconducting order parameter in Sr 2 RuO 4

Sr 2 RuO 4 has stood as the leading candidate for a spin-triplet superconductor for 26 years. However, recent NMR experiments have cast doubt on this candidacy and it is difficult to find a theory of superconductivity that is consistent with all experiments. Here, the order parameter symmetry for this material therefore remains an open question. Symmetry-based experiments are needed that can rule out broad classes of possible superconducting order parameters. Here, we use resonant ultrasound spectroscopy to measure the entire symmetry-resolved elastic tensor of Sr 2 RuO 4 through the superconducting transition. We observe a thermodynamic discontinuity in the shear elastic modulus c 66 , which implies that the superconducting order parameter has two components. A two-component p-wave order parameter, such as p x + ip y , naturally satisfies this requirement. As this order parameter appears to have been precluded by recent NMR experiments, we suggest that two other two-component order parameters, namely {d xz ,d yz } and {d x 2 –y 2 ,g xy (x 2 –y 2 ) }, are now the prime candidates for the order parameter of Sr 2 RuO 4 .

36 MATERIALS SCIENCE↗

Bardeen-Cooper-Schrieffer pairing of composite fermions

Topological pairing of composite fermions has led to remarkable ideas, such as excitations obeying non-Abelian braid statistics and topological quantum computation. Here, we construct a p-wave paired Bardeen-Cooper-Schrieffer (BCS) wave function for composite fermions in the torus geometry, which is a convenient geometry for formulating momentum space pairing as well as for revealing the underlying composite-fermion Fermi sea. Following the standard BCS approach, we minimize the Coulomb interaction energy at half filling in the lowest and the second Landau levels, which correspond to filling factors ν = 1/2 and ν = 5/2 in GaAs quantum wells, by optimizing two variational parameters that are analogous to the gap and the Debye cut-off energy of the BCS theory. Our results show no evidence for pairing at ν = 1/2 but a clear evidence for pairing at ν = 5/2. To a good approximation, the highest overlap between the exact Coulomb ground state at ν = 5/2 and the BCS state is obtained for parameters that minimize the energy of the latter, thereby providing support for the physics of composite-fermion pairing as the mechanism for the 5/2 fractional quantum Hall effect. We discuss the issue of modular covariance of the composite-fermion BCS wave function, and calculate its Hall viscosity and pair correlation function. By similar methods, we look for but do not find an instability to s-wave pairing for a spin-singlet composite-fermion Fermi sea at half-filled lowest Landau level in a system where the Zeeman splitting has been set to zero.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An Exact Separation of the Spin-Free and Spin-Dependent Terms of the Dirac-Coulomb-Breit Hamiltonian

The Dirac Hamiltonian is transformed by extracting the operator (sigma x p)/2mc from the small component of the wave function and applying it to the operators of the original Hamiltonian. The resultant operators contain products of Paull matrices that can be rearranged to give spin-free and spin-dependent operators. These operators are the ones encountered in the Breit-Pauli Hamiltonian, as well as some of higher order in alpha(sup 2). However, since the transformation of the original Dirac Hamiltonian is exact, the new Hamiltonian can be used in variational calculations, with or without the spin-dependent terms. The new small component functions have the same symmetry properties as the large component. Use of only the spin-free terms of the new Hamiltonian permits the same factorization over spin variables as in nonrelativistic theory, and therefore all the post-Self-Consistent Field (SCF) machinery of nonrelativistic calculations can be applied. However, the single-particle functions are two-component orbitals having a large and small component, and the SCF methods must be modified accordingly. Numerical examples are presented, and comparisons are made with the spin-free second-order Douglas-Kroll transformed Hamiltonian of Hess.

Dyall, Kenneth G.↗

Two-color polarization spectroscopy measurements of Zeeman state-to-state collision induced transitions of nitric oxide in binary gas mixtures

We investigated collision induced transitions in the (0, 0) band of the A2Σ+–X2Π electronic transition of nitric oxide (NO) using two-color polarization spectroscopy (TCPS). Two sets of TCPS spectra for 1% NO, diluted in different buffer gases at 295 K and 1 atm, were obtained with the pump beam tuned to the R11(11.5) and OP12(1.5) transitions. The buffer gases were He, Ar, and N2. The probe was scanned while the pump beam was tuned to the line center. Theoretical TCPS spectra, calculated by solving the density matrix formulation of the time-dependent Schrödinger wave equation, were compared with the experimental spectra. A collision model based on the modified exponential-gap law was used to model the rotational level-to-rotational level collision dynamics. A model for collisional transfer from an initial to a final Zeeman state was developed based on the difference in cosine of the rotational quantum number J projection angle with the z-axis for the two Zeeman states. Rotational energy transfer rates and Zeeman state collisional dynamics were varied to obtain good agreement between theory and experiment for the two different TCPS pump transitions and for the three different buffer gases. One key finding, in agreement with quasi-classical trajectory calculations, is that the spin-rotation changing transition rate in the A2Σ+ level of NO is almost zero for rotational quantum numbers ≥8. It was necessary to set this rate to near zero to obtain agreement with the TCPS spectra.

Chemistry↗

Unlocking New Regimes in Fractional Quantum Hall Effect with Quaternions

We demonstrate that formulating the composite-fermion theory of the fractional quantum Hall (FQH) effect in terms of quaternions greatly expands its reach and opens the door into many interesting issues that were previously not amenable to quantitative theoretical investigation. As an illustration, we explore the possibility of a nematic or a charge-density wave instability of the composite-fermion Fermi sea at half-filled Landau level and of the nearby FQH states by looking for a gap closing instability of the neutral magneto-roton excitation. As a result, our quaternion formulation of the FQH effect has been inspired by mathematical developments in the theoretical analyses of gravitational wave modes and cosmic microwave background radiation, where an important role is played by spin-weighted spherical harmonics that are nothing but monopole harmonics appearing in the spherical geometry for the FQH effect.

Composite fermions↗

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗