Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo method”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hierarchy among the crystal lattice, charge density wave, and superconducting orders in transition metal dichalcogenides

Using high-energy x-ray scattering and large-scale three-dimensional (3D) structure modeling, in this study we investigate the relationship between the crystal lattice, charge density wave (CDW), and superconducting (SC) orders in transition metal dichalcogenides (TMDs). In particular, we systematically substitute Te for Se in Ta-Se-Te solid solutions, determine changes in their crystal lattice, and relate them to changes in the CDW transition temperature, $T_{CDW}$, and SC critical temperature, $T_c$ . We find that strong lattice distortions such as buckling of Ta layers are detrimental to the CDW and SC orders. The presence of a perfect lattice order in two dimensions is a prerequisite to the emergence of CDWs but insufficient to achieve a SC ordered state. For the SC order to emerge, the Ta sublattice should also appear periodic in 3D. Local chemical disorder may promote the SC order, and the perfectness of Ta coordination polyhedra is a factor contributing to its strength. A hierarchical relationship among the crystal lattice, CDW, and SC orders thus appears to exist in TMDs in a sense that different degrees of crystal lattice order-disorder promote and maintain the CDW and SC orders to a different extent, offering an opportunity to control the latter through modifying the former by rational design. Our findings are a step towards a better understanding of the correlation between lattice and electronic degrees of freedom in TMDs. We also demonstrate an efficient experimental approach to study them in fine detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Starting-point-independent quantum Monte Carlo calculations of iron oxide

Quantum Monte Carlo (QMC) methods are useful for studies of strongly correlated materials because they are many body in nature and use the physical Hamiltonian. Typical calculations assume as a starting point a wave function constructed from single-particle orbitals obtained from one-body methods, e.g., density functional theory. However, mean-field-derived wave functions can sometimes lead to systematic QMC biases if the mean-field result poorly describes the true ground state. In this study, we examine the accuracy and flexibility of QMC trial wave functions using variational and fixed-node diffusion QMC estimates of the total spin density and lattice distortion of antiferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. We found that for relatively simple wave functions the predicted lattice distortion was controlled by the choice of single-particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. By optimizing the orbitals with QMC, we then demonstrate starting-point independence of the trial wave function with respect to the method by which the orbitals were constructed by demonstrating convergence of the energy, spin density, and predicted lattice distortion for two qualitatively different sets of orbitals. The results suggest that orbital optimization is a promising method for accurate many-body calculations of strongly correlated condensed phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

From hard spheres to hard-core spins

A system of hard spheres exhibits physics that is controlled only by their density. This comes about because the interaction energy is either infinite or zero, so all allowed configurations have exactly the same energy. The low-density phase is liquid, while the high-density phase is crystalline, an example of “order by disorder” as it is driven purely by entropic considerations. Here we study a family of hard spin models, which we call hard-core spin models, where we replace the translational degrees of freedom of hard spheres with the orientational degrees of freedom of lattice spins. Their hard-core interaction serves analogously to divide configurations of the many spin system into allowed and disallowed sectors. We present detailed results on the square lattice in d = 2 for a set of models with $\mathbb{Z}_n$ symmetry, which generalize Potts models, and their U(1) limits, for ferromagnetic and antiferromagnetic senses of the interaction, which we refer to as exclusion and inclusion models. As the exclusion and inclusion angles are varied, we find a Kosterlitz-Thouless phase transition between a disordered phase and an ordered phase with quasi-long-ranged order, which is the form order by disorder takes in these systems. These results follow from a set of height representations, an ergodic cluster algorithm, and transfer matrix calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Oxygen magnetic polarization, nodes in spin density, and zigzag spin order in oxides

Recent studies on Ba 2 CoO 4 (BCO) and SrRuO 3 (SRO) have unveiled a variety of intriguing phenomena, such as magnetic polarization on oxygens, unexpected nodes in the spin density profile along bonds, and unusual zigzag spin patterns in triangular lattices. Here, using simple model calculations supplemented by density functional theory we explain the presence of nodes based on the antibonding character of the dominant singly occupied molecular orbitals along the transition metal (TM) to oxygen bonds. Our simple model also allows us to explain the net polarization on oxygen as originated from the hybridization between atoms and mobility of the electrons with spins opposite to those of the closest TM atoms. Overall, our results are not limited to BCO and SRO, but they are generic and qualitatively predict the net polarization expected on any ligands, according to the spin order of the closest TM atoms and the number of intermediate ligand atoms. Finally, we propose that a robust easy-axis anisotropy would suppress the competing 120° antiferromagnetic order to stabilize the zigzag pattern order as ground state in a triangular lattice. Our generic predictions should be applicable to any other compound with characteristics similar to those of BCO and SRO.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Defect-driven ferrimagnetism and hidden magnetization in MnBi 2 Te 4

MnBi 2 Te 4 (MBT) materials are promising antiferromagnetic topological insulators in which field-driven ferromagnetism is predicted to cause a transition between axion insulator and Weyl semimetallic states. However, the presence of antiferromagnetic coupling between Mn/Bi antisite defects and the main Mn layer can reduce the low-field magnetization, and it has been shown that such defects are more prevalent in the structurally identical magnetic insulator MnSb 2 Te 4 (MST). We use high-field magnetization measurements to show that the magnetization of MBT and MST occur in stages and full saturation requires fields of ~60 T. As a consequence, the low-field magnetization plateau state in MBT, where many determinations of the quantum anomalous Hall state are studied, actually consists of ferrimagnetic septuple blocks containing both uniform and staggered magnetization components.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Temperature evolution of the phonon dynamics in the Kitaev spin liquid

In this work, we present a study of the phonon dynamics in the honeycomb Kitaev spin model at finite temperatures. We show that the fractionalized spin excitations of the Kitaev spin liquid, the itinerant Majorana fermions, and static Z 2 fluxes have distinct effects on the phonon dynamics, which makes the phonon dynamics a promising tool for exploring the Kitaev spin liquid candidate materials. In particular, we focus on the signature of the fractionalized excitations in the thermodynamic behavior of the sound attenuation and the phonon Hall viscosity: the former describes the phonon decay into the fractionalized excitations, and the latter is the leading order time-reversal symmetry breaking effect on the acoustic phonon. We find that the angular dependence of the attenuation coefficient and its magnitude is modified by the thermal excitation of the Z 2 fluxes. The strength of this effect strongly depends on the relative magnitude of the sound velocity and the Fermi velocity characterizing the low-energy Majorana fermions. We also show that the Hall viscosity is strongly suppressed by the increase of the density of the Z 2 fluxes at finite temperatures. All these observations reflect the effects of the emergent disorder on the Majorana fermions introduced by the Z 2 fluxes. Our analysis is based on the complementary analytical calculations in the low-temperature zero-flux sector, and numerical calculations in the inhomogeneous flux sectors at intermediate and high temperatures with stratified Monte Carlo (strMC) method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Spin dynamics of the antiferromagnetic Heisenberg model on a kagome bilayer

We research the spin dynamics of a classical Heisenberg antiferromagnet with nearest-neighbor interactions on a quasi-two-dimensional kagome bilayer. This geometrically frustrated lattice consists of two kagome layers connected by a triangular-lattice linking layer. By combining Monte Carlo with precessional spin dynamics simulations, we compute the dynamical structure factor of the classical spin liquid in the kagome bilayer and investigate the thermal and dilution effects. While the low-frequency and long-wavelength dynamics of the cooperative paramagnetic phase is dominated by spin diffusion, weak magnon excitations persist at higher energies, giving rise to the half-moon pattern in the dynamical structure factor. In the presence of spin vacancies, the dynamical properties of the diluted system can be understood within the two-population picture. The spin diffusion of the “correlated” spin clusters is mainly driven by the zero-energy weather-vane modes, giving rise to an autocorrelation function that decays exponentially with time. On the other hand, the diffusive dynamics of the quasifree “orphan” spins leads to a distinctive longer-time power-law tail in the autocorrelation function. We discuss the implications of our work for the glassy behaviors observed in the archetypal frustrated magnet SrCr 9p Ga 12-9 pO 19 (SCGO).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Quantifying multipoint ordering in alloys

A central problem in multicomponent lattice systems is to systematically quantify multipoint ordering. Ordering in such systems is often described in terms of pairs, even though this is not sufficient when three-point and higher-order interactions are included in the Hamiltonian. Current models and parameters for multipoint ordering are often only applicable for very specific cases or require approximating a subset of correlated occupational variables on a lattice as being uncorrelated. In this paper, cluster order parameters are introduced to systematically quantify arbitrary multipoint ordering motifs in substitutional systems through direct calculations of normalized cluster probabilities. These parameters can describe multipoint chemical ordering in crystal systems with multiple sublattices, multiple components, and systems with reduced symmetry. These are defined in this paper and applied to quantify four-point chemical ordering motifs in platinum/palladium alloy nanoparticles that are of practical interest to the synthesis of catalytic nanocages. Impacts of chemical ordering on nanocage stability are discussed. It is demonstrated that approximating four-point probabilities from superpositions of lower-order pair probabilities is not sufficient in cases where three- and four-body terms are included in the energy expression. Furthermore, conclusions about the formation mechanisms of nanocages may change significantly when using common pair approximations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Importance of bulk excitations and coherent electron-photon-phonon scattering in photoemission from PbTe(111): Ab initio theory with experimental comparisons

This paper presents a fully ab initio many-body photoemission framework that includes coherent three-body electron-photon-phonon scattering to predict the transverse momentum distributions and the mean transverse energies (MTEs) of bulk photoelectrons from single-crystal photocathodes. The need to develop such a theory stems from the lack of studies that provide complete understanding of the underlying fundamental processes governing the transverse momentum distribution of photoelectrons emitted from single crystals. For example, initial predictions based on density-functional theory calculations of effective electron masses suggested that the (111) surface of PbTe would produce very small MTEs (≤ 15 meV), whereas our experiments yielded MTEs ten to twenty times larger than these predictions and also exhibited a lower photoemission threshold than predicted. The ab initio framework presented in this paper correctly reproduces both the magnitude of the MTEs from our measurements in PbTe(111) and the observed photoemission below the predicted threshold. Our results show that both photoexcitations into states that propagate in the bulk of the material and coherent many-body electron-photon-phonon scattering processes, which initial predictions ignored, play surprisingly important roles in photoemission from PbTe(111). Lastly, from the lessons learned, we recommend a procedure for rapid computational screening of potential single-crystal photocathodes for applications in next-generation ultrafast electron diffraction and X-ray free-electron lasers, which will enable significant advances in condensed matter research.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fixed-phase diffusion Monte Carlo study of activation gap and skyrmion excitations of a ν = 1 system in the presence of charged impurities

The discrepancy between the theoretically calculated and experimentally measured activation gaps in quantum Hall effect has long been a puzzle. We revisit this issue in the context of the v = 1 quantum Hall state, while also incorporating the skyrmion physics. We find that the finite width and the Landau level mixing effects are not sufficient to explain the observed activation gap. We further show that the presence of charged impurities located adjacent to the quantum well can cause a significant reduction in the activation gap, while also causing a suppression of the skyrmion size.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Trapping interlayer excitons in van der Waals heterostructures by potential arrays

Transition metal dichalcogenide heterostructures can host interlayer excitons (IXs), which consist of electrons and holes spatially separated in different layers. IXs possess permanent dipoles and have proven to offer a wealth of novel physics. Here we develop a discrete, random-walk model which includes annihilation and repulsion interactions among IXs. Using this model, we simulate the trapping of IXs in traps of different depths, densities, and shapes. Our results show that dipole-dipole interactions play an important role in regulating IX trapping. The effects of dipole interactions can be mitigated with small, deep traps which are realizable with atomic defects and moire potentials.

2-dimensional systems↗
DOE OSTI2021

Physical properties of the quasi-two-dimensional square lattice antiferromagnet Ba 2 FeSi 2 O 7

In this work, we report magnetization (χ, M ), magnetic specific heat ( C M ), and neutron powder diffraction results on a quasi-two-dimensional (2D) S = 2 square lattice antiferromagnet Ba 2 FeSi 2 O 7 consisting of FeO 4 tetrahedrons with highly compressive tetragonal distortion (27%). Despite of the quasi-2D lattice structure, both χ and C M present three-dimensional magnetic long-range ordering below the Néel temperature T N = 5.2 K . Neutron diffraction data show a collinear Q m = (1,0,1/2) antiferromagnetic (AFM) structure below T N but the ordered moment aligned in the a b plane is suppressed by 26% from the ionic spin S = 2 value ( 4 μ B ). Both the AFM structure and the suppressed moments are well explained by using Monte Carlo simulations with a large single-ion in-plane anisotropy D = 1 .4 meV and a rather small Heisenberg exchange J intra = 0 .15 meV in the plane. The characteristic 2D spin fluctuations are recognized in the magnetic entropy release and diffuse scattering above T N . This new quasi-2D magnetic system also displays unusual nonmonotonic dependence of T N as a function of magnetic field H .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reaching the equilibrium state of the frustrated triangular Ising magnet Ca 3 Co 2 O 6

Ca 3 Co 2 O 6 is a frustrated magnet consisting of a triangular arrangement of chains of Ising spins. It shows regular magnetization steps versus magnetic field every 1.2 T that are metastable with very slow dynamics. This has puzzled the community for many years and given rise to numerous potential theories. Here we approach the problem by seeking the elusive magnetic equilibrium state at T=2 K. To this end, we explore two approaches: (1) bypassing the slow dynamics produced by changing fields by instead field-cooling directly to the target temperature and (2) quantum annealing in transverse magnetic fields. While we observe no measurable effect of the quantum annealing in fields up to 7 T, which is likely due to the large Ising anisotropy of Co spins in this material, we find that for the field cooling in longitudinal fields we achieve the predicted equilibrium 1/3 magnetization. We perform Monte Carlo simulations of the ground state phase diagram and we also simulate the quantum annealing process and find good agreement between experiment and theory. Thus we present an investigation of the elusive ground state properties of the canonical frustrated triangular system Ca 3 Co 2 O 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

High-pressure phase transitions and melt structure of Pb O 2 : An analog for silica

In situ x-ray diffraction measurements and inverse Monte Carlo simulations of pair distribution functions were used to study the structural response of PbO 2 under pressure and to characterize the local structure of liquid PbO 2 . Two phase transitions are observed upon room-temperature compression of crystalline PbO 2 up to ~65 GPa. Here, the starting mixture of rutile structured β– PbO 2 and orthorhombic α– PbO 2 undergoes a transition to a ZrO 2 -type orthorhombic phase with space group Pbca at ~21 GPa. Above 42 GPa, the Pbca phase transforms to a cotunnite-type phase with space group Pnam. In this paper, we also report short-range order in liquid PbO 2 , showing ~8–9-fold and 4–5-fold coordination around Pb and O atoms, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnetism and magnetotransport in the kagome antiferromagnet Mn 3 Ge

We perform classical Monte Carlo and stochastic Landau-Lifshitz-Gilbert simulations to study the temperature-dependent magnetism of the kagome antiferromagnet Weyl metal Mn3Ge, and we find that a long-range chiral order sets in at a transition temperature well below the Neel temperature (T N ). Based on the crystalline symmetries imposed by the chiral magnetic order, we argue for the presence of multiple isoenergetic Weyl nodes (nodes that are at the same energy and with a congruent Fermi surface around them) near the chemical potential. Using the semiclassical Boltzmann equations, we show that the combined contribution to the net longitudinal magnetoconductance (LMC) and the planar Hall conductance (PHC) from tilted Weyl nodes can lead to signatures that are qualitatively distinct from that of a single pair of Weyl nodes. In particular, we show that magnetic orders with different chiralities can give rise to different periods in LMC and PHC as a function of the in-plane magnetic field direction. As a result, this is ultimately related to differences in the symmetry-imposed constraints on the Weyl nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnetic structure and exchange interactions in the Heisenberg pyrochlore antiferromagnet Gd 2 Pt 2 O 7

In this report the Heisenberg pyrochlore antiferromagnet Gd 2 Pt 2 O 7 is one of a series of gadolinium pyrochlore compounds with a variety of B-site cations. Despite the expected simplicity of a spin-only Gd 3+ Heisenberg interaction model, the gadolinium pyrochlore series exhibits various complex magnetic ground states at low temperature. Gd 2 Pt 2 O 7 displays the highest temperature magnetic order of the series with $T_N$=1.6 K, which has been attributed to enhanced superexchange pathways facilitated by empty $5d$ $e_g$ Pt orbitals. In this study, we use various neutron scattering techniques on an isotopically enriched polycrystalline 160 Gd 2 Pt 2 O 7 sample to examine the magnetic structure and spin-wave excitation spectrum below $T_N$ in order to extract the dominant exchange interactions. We find that the ground-state magnetic structure is the Palmer-Chalker state previously seen in Gd 2 Sn 2 O 7 with an associated gapped excitation spectrum consistent with enhanced exchange interactions between further near-neighbor Gd 3+ ions. We confirm this exchange model with analysis of the magnetic diffuse scattering in the paramagnetic regime using polarized neutrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effect of emitters on quantum state transfer in coupled cavity arrays

Over the last decade, conditions for perfect state transfer in quantum spin chains have been discovered and their experimental realizations addressed, as have their extensions to more complex geometries of coupled cavity-emitter arrays. In this paper we further consider such studies and situations in which quantum state transfer can occur with high fidelity, even when the cavity-cavity coupling rates and cavity-emitter interaction rates are comparable. This is accomplished through the development and use of a Monte Carlo approach to the inverse eigenvalue problem, which allows the determination of coupling rates which optimize quantum state transfer fidelity and subsequent time evolution of the polariton wave function through exact diagonalization of the resulting Jaynes-Cummings-Hubbard Hamiltonian. The effect of inhomogeneous emitter locations is also evaluated. Furthermore, our key results include the demonstration that our methodology can be used successfully to establish Hamiltonian parameters for high-fidelity state transfer in more general lattice geometries and excitation number sectors, and also a determination of the effects of fluctuations in those parameters about their optimal values.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Revisiting excitation gaps in the fractional quantum Hall effect

Recent systematic measurements of the quantum well width dependence of the excitation gaps of fractional quantum Hall states in high mobility samples open the possibility of a better quantitative understanding of this important issue. Here we present what we believe to be accurate theoretical gaps including the effects of finite width and Landau level (LL) mixing. While theory captures the width dependence, there still remains a deviation between the calculated and the measured gaps, presumably caused by disorder. It is customary to model the experimental gaps of the n/(2n ± 1) states as Δ n/(2n+1) = Ce 2 /[(2n ± 1)ϵl]-Γ, where ϵ is the dielectric constant of the background semiconductor and l is the magnetic length; the first term is interpreted as the cyclotron energy of composite fermions and Γ as a disorder-induced broadening of composite-fermion LLs. Fitting the gaps for various fractional quantum Hall states, we find that Γ can be nonzero even in the absence of disorder.

36 MATERIALS SCIENCE↗