Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite size”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Universal relation between doping content and normal-state resistance in gate voltage tuned ultrathin Bi 2 Sr 2 CaCu 2 O 8 + x flakes

We report gate voltage tunable ultra-thin high-T c cuprates supply a unique platform to investigate the electronic phase diagram and superconductor-insulator transition. One of the challenges in this field is the precise determination of the doping content in the underdoped non-superconducting region. Here we report the discovery of a universal relation between the doping content p and the normal-state resistance at a fixed temperature R(T f ), p = α + β ln[1/R(T f )], in the ultra-thin Bi 2 Sr 2 CaCu 2 O 8+x flakes. The in-depth analysis shows that the evolution of carrier scattering probability with doping content and the change of effective mass caused by superconductor-insulator transition are two key factors leading to this logarithmic relation. Based on our finding, the more precise electronic phase diagram can be established. In addition, the superconductor-insulator transition is verified to be a quantum phase transition using a finite size scaling analysis. The scaling exponent z v is found to have a close correlation with the disorder levels. The present result provides an important foundation to investigate the fascinating electronic states in the ultra-thin cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Bilayer Hubbard model: Analysis based on the fermionic sign problem

The bilayer Hubbard model describes the antiferromagnet to spin singlet transition and, potentially, aspects of the physics of unconventional superconductors. Despite these important applications, significant aspects of its phase diagram in the interplane hopping t ⊥ versus on-site interaction U parameter space, at half filling, are largely in disagreement. Here we provide an analysis making use of the average sign of weights over the course of the importance sampling in quantum Monte Carlo simulations to resolve several central open questions. Specifically, this metric of the weights clarifies the finite-sized metallic regimes at small U. Furthermore, at strong interactions, it points to the existence of a crossover from a correlated to uncorrelated band insulator not yet explored in a variety of existing, unbiased numerical methods. Our paper demonstrates the versatility of using properties of the weights in quantum Monte Carlo simulations to reveal important physical characteristics of the models under study.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fully self-consistent finite-temperature $GW$ in Gaussian Bloch orbitals for solids

In this work, we present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated, finite-temperature sc ⁢$GW$ implementations when finite-size corrections and basis-set errors are taken into account. By migrating computationally intensive calculations to graphics processing units, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc⁢ $GW$ for ab initio correlated material simulations and provides a sound starting point for embedding methods built on top of $GW$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Plasmons in Z 2 topological insulators

Here, we study plasmonic excitations in the Kane-Mele model, a two-dimensional Z 2 topological insulator on the honeycomb lattice, using the random phase approximation (RPA). In the topologically nontrivial phase, the model has conducting edge states that traverse the bulk energy gap and display spin-momentum locking. Such a state of matter is called the quantum spin hall (QSH) phase, which is robust against time-reversal (TR) invariant perturbations. We find that in the QSH phase, gapless spin-polarized plasmons can be excited on the edges of the system. The propagation of these plasmons is chiral for each individual spin component and shows spin-momentum locking for both spin components on the same edge. Moreover, we study the effect of external magnetic fields on the gapless edge plasmons. Specifically, out-of-plane magnetic fields delocalize edge plasmons propagating in one direction without affecting the other one, while an in-plane magnetic field can be applied to selectively excite a specific spin-plasmon branch with proper doping or gating to the system. Our findings may have potential applications in novel plasmonic and spintronic devices. We also investigate plasmons in the Kane-Mele model on a finite-sized diamond-shaped nanoflake and observe low-energy plasmons circulating the boundary of the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Dipole condensates in tilted Bose-Hubbard chains

Here we study the quantum phase diagram of a Bose-Hubbard chain whose dynamics conserves both a boson number and boson dipole moment, a situation which can arise in strongly tilted optical lattices. The conservation of the dipole moment has a dramatic effect on the phase diagram, which we analyze by combining a field theory analysis with DMRG simulations. In the thermodynamic limit, the phase diagram is dominated by various types of incompressible dipolar condensates. In finite-sized systems, however, it may be possible to stabilize a Bose-Einstein insulator: an exotic compressible phase which is insulating, despite the absence of a charge gap. We suggest several ways by which these exotic phases can be identified in near-term cold-atom experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Analytic calculation of the vison gap in the Kitaev spin liquid

Although the ground-state energy of the Kitaev spin liquid can be calculated exactly, the associated vison gap energy has to date only been calculated numerically from finite size diagonalization. Here we show that the phase shift for scattering Majorana fermions off a single bond-flip can be calculated analytically, leading to a closed-form expression for the vison gap energy Δ = 0.2633 J . Generalizations of our approach can be applied to Kitaev spin liquids on more complex lattices such as the three dimensional hyper-octagonal lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Itinerant ferromagnetism in transition metal dichalcogenide moiré superlattices

Moiré materials are artificial crystals formed at van der Waals heterojunctions that have emerged as a highly tunable platform that is able to realize much of the rich quantum physics of electrons in atomic scale solids, and in several cases even new quantum phases of matter. Here we use finite-size exact diagonalization methods to explore the physics of single-band itinerant electron ferromagnetism in semiconductor moiré materials. As a result, we predict where ferromagnetism is likely to occur in triangular-lattice moiré systems, and where it is likely to yield the highest Curie temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Origin of spin-driven ferroelectricity and effect of external pressure on the complex magnetism of the 6⁢ H perovskite Ba 3 ⁢Ho⁢Ru 2⁢ O 9

The compound Ba 3 HoRu 2 O 9 magnetically orders at 50 K (T N1 ), followed by another complex magnetic ordering at 10.2 K (T N2 ). The second magnetic phase transition was characterized by the coexistence of two competing magnetic ground states associated with two different magnetic wave vectors (K 1 = 0.5 0 0 and K 2 = 0.25 0.25 0). The multiferroicity and magnetoelectric coupling were predicted below T N2 in these 4d-based materials. Here, in this work, we have discussed the origin of spin-driven ferroelectricity, which is not known yet, and the nature of magnetoelectric domains. We have investigated the compound through time-of-flight neutron diffraction, synchrotron x-ray diffraction (XRD), ac susceptibility, frequency-dependent complex dielectric spectroscopy, and dc magnetization under external pressure. We have demonstrated that the noncollinear structure involving two different magnetic ions, Ru (4d) and Ho (4f), breaks the spatial inversion symmetry via inverse Dzyaloshinskii-Moriya (DM) interaction through strong 4d-4f magnetic correlation, which shifts the oxygen atoms and results in nonzero polarization. Such an observation of inverse DM interaction from two different magnetic ions which cause ferroelectricity is rarely observed. The stronger spin-orbit coupling of 4d orbital might play a major role in creating DM interaction of noncollinear spins. We have systematically studied the spin and dipolar dynamics, which exhibit intriguing behavior with shorter coherence lengths of second magnetic phase associated with the k 2 wave vector. The results manifest the development of finite-size magnetoelectric domains instead of true long-range ordering, which justifies the experimentally obtained low value of ferroelectric polarization. The lattice parameters and volume show a sharp anomaly at T N2 obtained by analyzing the temperature-dependence XRD, which is consistent with the ferroelectric transition, predicting a noncentrosymmetric space group, $P\bar{6}2c$, for this compound. Furthermore, we have investigated the effect of external pressure on this complex magnetism. The result reveals an enhancement of ordering temperature by the application of external pressure (~1.6 K/GPa). The external pressure might favor stabilizing the magnetic ground state associated with second magnetic phase. Our study shows an unconventional mechanism of spin-driven ferroelectricity involving inverse DM interaction between Ru (4d) and Ho (4f) magnetic ions due to strong 4d-4f cross coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Classical dynamics of the antiferromagnetic Heisenberg spin ladder

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. Here, 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.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Properties of the neutron star crust: Quantifying and correlating uncertainties with improved nuclear physics

A compressible liquid-drop model (CLDM) is used to correlate uncertainties associated with the properties of the neutron star (NS) crust with theoretical estimates of the uncertainties associated with the equation of state (EOS) of homogeneous neutron and nuclear matter. For the latter, we employ recent calculations based on Hamiltonians constructed using chiral effective-field theory (χEFT). Fits to experimental nuclear masses are employed to constrain the CLDM further, and we find that they disfavor some of the χEFT Hamiltonians. The CLDM allows us to study the complex interplay between bulk, surface, curvature, and Coulomb contributions, and their impact on the NS crust. It also reveals how the curvature energy alters the correlation between the surface energy and the bulk symmetry energy. Here our analysis quantifies how the uncertainties associated with the EOS of homogeneous matter implies significant uncertainties for the composition of the crust, its proton fraction, and the volume fraction occupied by nuclei. We find that the finite-size effects impact the crust composition but have a negligible effect on the net isospin asymmetry of matter. The isospin asymmetry is largely determined by the bulk properties and the isospin dependence of the surface energy. The most significant uncertainties associated with matter properties in the densest regions of the crust, the precise location of the crust-core transition, are found to be strongly correlated with uncertainties associated with the Hamiltonians. By adopting a unified model to describe the crust and the core of NSs, we tighten the correlation between their global properties such as their mass-radius relationship, moment of inertia, crust thickness, and tidal deformability with uncertainties associated with the nuclear Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Relativistic mean-field theories for neutron-star physics based on chiral effective field theory

We describe and implement a procedure for determining the couplings of a relativistic mean-field theory (RMFT) that is optimized for application to neutron star phenomenology. In the standard RMFT approach, the couplings are constrained by comparing the theory's predictions for symmetric matter at saturation density with measured nuclear properties. The theory is then applied to neutron stars which consist of neutron-rich matter at densities ranging up to several times saturation density, which allows for additional astrophysical constraints. In our approach, rather than using the RMFT to extrapolate from symmetric to neutron-rich matter and from finite-sized nuclei to uniform matter, we fit the RMFT to properties of uniform pure neutron matter obtained from chiral effective field theory. Chiral effective field theory incorporates the experimental data for nuclei in the framework of a controlled expansion for nuclear forces valid at nuclear densities and enables us to account for theoretical uncertainties when fitting the RMFT. We construct four simple RMFTs that span the uncertainties provided by chiral effective field theory for neutron matter, and are consistent with current astrophysical constraints on the equation of state. Lastly, our RMFTs can be used to model the properties of neutron-rich matter across the vast range of densities and temperatures encountered in neutron stars and their mergers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Lattice sine-Gordon model

We obtain nonperturbative results on the sine-Gordon model using the lattice field technique. In particular, we employ the Fourier accelerated hybrid Monte Carlo algorithm for our studies. We find the critical temperature of the theory based on autocorrelation time, as well as the finite size scaling of the "thickness" observable used in an earlier lattice study by Hasenbusch et al. We study the entropy, which is smooth across all temperatures, supportive of an infinite order transition. This system has a well-known duality with the massive Thirring model, which can play the role of a toy models for Montonen-Olive duality in $\mathcal {N}$ = 4 super-Yang-Mills theory, since it relates solitons to elementary field excitations. Our research lays a groundwork for such study on the lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamics of non-Gaussian fluctuations in model A

Motivated by the experimental search for the QCD critical point, we perform simulations of a stochastic field theory with purely relaxational dynamics (model A). We verify the expected dynamic scaling of correlation functions. Using a finite size scaling analysis, we obtain the dynamic critical exponent z = 2.026(56). We investigate time dependent correlation functions of higher moments M n (t) of the order parameter M(t) for n = 1, 2, 3, 4. We obtain dynamic scaling with the same critical exponent z for all n, but the relaxation constant depends on n. We also study the relaxation of M n (t) after a quench, where the simulation is initialized in the high temperature phase, and the dynamics is studied at the critical temperature T c . Finally, we find that the evolution does not follow simple scaling with the dynamic exponent z, and that it involves an early time rise followed by late stage relaxation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Detectability of subsolar mass neutron stars through a template bank search

Here, we study the detectability of gravitational-wave signals from subsolar-mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as 78.4% for an equal mass binary system of chirp mass 0.17M ⊙ , in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for subsolar-mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of subsolar-mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for subsolar-mass binary neutron star signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Nonperturbative running of the tensor operator for N f = 3 QCD from the chirally rotated Schrödinger functional

We study the renormalization group (RG) running of the nonsinglet tensor operator, for N f = 3 QCD with Wilson fermions in a mixed action setup, with standard Schrödinger functional (SF) boundary conditions for sea quarks and chirally rotated Schrödinger functional ( χ SF ) boundary conditions for valence quarks. Based on a recursive finite-size scaling technique we compute nonperturbatively the tensor step-scaling function for an energy range between a hadronic scale and an electroweak scale, above which perturbation theory may be safely applied. Our result is expressed as the RG-running factor T RGI / [ T ( μ had ) ] R , where the numerator is the scale independent (renormalization group invariant—RGI) tensor operator and the denominator is its renormalized counterpart at a hadronic scale μ had = 233 ( 8 ) MeV in a given scheme. We determine the step-scaling function in four distinct renormalization schemes. We also compute the renormalization parameters of these schemes at μ had which, combined with the RG-running factor, gives the scheme-independent quantity Z T RGI ( g 0 2 ) in four schemes and for a range of bare gauge couplings in which large volume hadronic matrix element simulations are performed by the CLS consortium in N f = 2 + 1 QCD. All four results are compatible and also agree with a recent determination based on a unitary setup for Wilson quarks with Schrödinger functional boundary conditions [arXiv:2309.04314]. This provides a strong universality test. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗