Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explainable”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Unification of nonequilibrium molecular dynamics and the mode-resolved phonon Boltzmann equation for thermal transport simulations

Nano-size confinement induces many intriguing non-Fourier heat conduction phenomena, such as nonlinear temperature gradients, temperature jumps near the contacts, and size-dependent thermal conductivity. Over the past decades, these effects have been studied and interpreted by nonequilibrium molecular dynamics (NEMD) and phonon Boltzmann transport equation (BTE) simulations separately, but no theory that unifies these two methods has ever been established. In this work, we unify these methods using a quantitative mode-level comparison and demonstrate that they are equivalent for various thermostats. We show that different thermostats result in different non-Fourier thermal transport characteristics due to the different mode-level phonon excitations inside the thermostats, which explains the different size-dependent thermal conductivities calculated using different reservoirs, even though they give the same bulk thermal conductivity. Specifically, the Langevin thermostat behaves like a thermalizing boundary in phonon BTE and provides mode-level thermal-equilibrium phonon outlets, while the Nose-Hoover chain thermostat and velocity rescaling method behave like biased reservoirs, which provide a spatially uniform heat generation and mode-level nonequilibrium phonon outlets. These findings explain why different experimental measurement methods can yield different size-dependent thermal conductivity. They also indicate that the thermal conductivity of materials can be tuned for various applications by specifically designing thermostats. Finally, the unification of NEMD and phonon BTE will largely facilitate the study of thermal transport in complex systems in the future by, e.g., replacing computationally unaffordable first-principles NEMD simulations with computationally less expensive spectral BTE simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Low-energy physics of isotropic spin-1 chains in the critical and Haldane phases

Using a matrix product state algorithm with infinite boundary conditions, we compute high-resolution dynamic spin and quadrupolar structure factors in the thermodynamic limit to explore the low-energy excitations of isotropic bilinear-biquadratic spin-1 chains. Haldane mapped the spin-1 Heisenberg antiferromagnet to a continuum field theory, the nonlinear sigma model (NLσM). In this work, we find that the NLσM fails to capture the influence of the biquadratic term and provides only an unsatisfactory description of the Haldane phase physics. But several features in the Haldane phase can be explained by noninteracting multimagnon states. The physics at the Uimin-Lai-Sutherland point is characterized by multisoliton continua. Moving into the extended critical phase, we find that these excitation continua contract, which we explain using a field-theoretic description. New excitations emerge at higher energies and, in the vicinity of the purely biquadratic point, they show simple cosine dispersions. Using block fidelities, we identify them as elementary one-particle excitations and relate them to the integrable Temperley-Lieb chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electrically controlled two-dimensional electron-hole fluids

We study the electronic properties of dual-gated electron-hole bilayers in which the two layers are separated by a perfectly opaque tunnel barrier. Combining an electrostatic and thermodynamic analysis with mean-field theory estimates of interacting system chemical potentials, we explain the dependence of the electron and hole densities on the two gate voltages. Because chemical potential jumps occur for both electrons and holes at neutrality, there is a finite area in gate voltage parameter space over which electron and hole densities are equal. In that regime the electron-hole pair density depends only on the sum of the two gate voltages. As a result, we are able to explain a recent experimental study of electrically controlled bilayers by allowing for interlayer tunneling and using a nonequilibrium steady-state transport picture.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spin-orbit exciton in a honeycomb lattice magnet CoTiO 3 : Revealing a link between magnetism in d - and f -electron systems

We carried out inelastic neutron scattering to study the spin-orbital (SO) exciton in a single crystal sample of CoTiO 3 as a function of temperature. CoTiO 3 is a honeycomb magnet with dominant XY-type magnetic interaction and an A-type antiferromagnetic order below T N ≈ 38 K. We found that the SO exciton becomes softer, but acquires a larger bandwidth in the paramagnetic phase, compared to that in the magnetically ordered phase. Moreover, an additional mode is only observed in the intermediate temperature range, as the sample is warmed up above the lowest accessible temperature below T N . Such an unusual temperature dependence observed in this material suggests that its ground states (an S eff = 1 2 doublet) and excited states multiplets are strongly coupled, and therefore cannot be treated independently, as often done in a pseudo-spin model. Furthermore, our observations can be explained by a multi-level theory within random phase approximation that explicitly takes into account both the ground and excited multiplets. The success of our theory, which is originally developed to explain temperature dependence of magnetic excitations in the rare-earth magnets, highlight the similarity between the magnetic excitations in rare-earth systems and those in transition metal systems with strong spin orbit coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Ba ( Ti 1 - x , Zr x ) O 3 relaxors: Dynamic ferroelectrics in the gigahertz frequency range

Relaxor ferroelectrics are of great scientific and technological significance as they exhibit large and unusual responses to external stimuli. Their hallmark features are broadness and frequency dispersion of the peak in the temperature dependence of the dielectric constant. Both are believed to originate from the dynamics of polar nanoregions. Here, we apply first-principles-based molecular dynamics to resolve the gigahertz electric response of dynamically poled Ba (Ti 1 - x , Zr x ) O 3 ferroelectric relaxors that remained overlooked so far. We find that (i) the hallmark relaxor features continue to persist even in the dynamically poled structures, but do not necessarily originate from the polar nanoregions dynamics; (ii) dynamically poled samples exhibit polarization aging which leads to the frequency dependence of both remnant polarization and the Curie temperature; (iii) “dynamical” nature of the latter naturally explains the frequency dependence of the dielectric susceptibility maximum in dynamically poled Ba (Ti 1 - x , Zr x) O 3 ; and (iv) incorporation of the polarization aging into dielectric susceptibility expression explains its enhancing contribution in an intuitive way.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Oxygen vacancy induced site-selective Mott transition in LaNiO 3

While defects such as oxygen vacancies in correlated materials can modify their electronic properties dramatically, understanding the microscopic origin of electronic correlations in materials with defects has been elusive. Lanthanum nickelate with oxygen vacancies, LaNiO 3 - x , exhibits the metal-to-insulator transition as the oxygen vacancy level x increases from the stoichiometric LaNiO 3 . In particular, LaNiO 2.5 exhibits a paramagnetic insulating phase, also stabilizing an antiferromagnetic state below T N ≃ 152 K . Here, we study the electronic structure and energetics of LaNiO 3 - x using first principles. We find that LaNiO 2.5 exhibits a “site-selective” paramagnetic Mott insulating state at T ≃ 290 K as obtained using density functional theory plus dynamical mean field theory ( DFT + DMFT ). The Ni octahedron site develops a Mott insulating state with strong correlations as the Ni e g orbital is half-filled while the Ni square-planar site with apical oxygen vacancies becomes a band insulator. Overall, our oxygen vacancy results cannot be explained by the pure change of the Ni oxidation state alone within the rigid band-shift approximation. Our DFT + DMFT density of states explains that the peak splitting of unoccupied states in LaNiO 3 - x measured by the experimental x-ray absorption spectra originates from two nonequivalent Ni ions in the vacancy-ordered structure.

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

Collective excitations in the tetravalent lanthanide honeycomb antiferromagnet

Thermomagnetic and inelastic neutron-scattering measurements on Na 2 PrO 3 are reported. This material is an antiferromagnetic honeycomb magnet based on the tetravalent lanthanide Pr 4+ and has been proposed to host dominant antiferromagnetic Kitaev interactions. These measurements reveal magnetic fluctuations in Na 2 PrO 3 below an energy of 2 meV as well as crystal-field excitations around 230 meV. Furthermore, the latter energy is comparable to the scale of the spin-orbit interaction and explains both the very small effective moment of around 1.0μB per Pr 4+ and the difficulty to uncover any static magnetic scattering below the ordering transition at T N = 4.6 K. By comparing the low-energy magnetic excitations in Na 2 PrO 3 to those of the isostructural spin-only compound, Na 2 TbO 3 , a microscopic model of exchange interactions is developed that implicates dominant and surprisingly large Heisenberg exchange interactions J≈1.1(1) meV. Although antiferromagnetic Kitaev interactions with K≤0.2J cannot be excluded, the inelastic neutron-scattering data of Na 2 PrO 3 is best explained with a Δ=1.22 easy-axis XXZ exchange anisotropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Role of carbon and hydrogen in limiting n -type doping of monoclinic ( Al x Ga 1 – x ) 2 O 3

In this study, we use hybrid density functional calculations to assess n -type doping in monoclinic ( Al x Ga 1 – x ) 2 O 3 alloys. We focus on silicon, the most promising donor dopant, and study the structural properties, formation energies, and charge-state transition levels of its various configurations. We also explore the impact of carbon and hydrogen, which are common impurities in metal-organic chemical vapor deposition (MOCVD). In Ga 2 O 3 , Si Ga is an effective shallow donor, but in Al 2 O 3 Si Al acts as a DX center with a ( + / – ) transition level in the band gap. Interstitial hydrogen acts as a shallow donor in Ga 2 O 3 but behaves as a compensating acceptor in n -type Al 2 O 3 . Interpolation indicates that Si is an effective donor in ( Al x Ga 1 – x ) 2 O 3 up to 70% Al, but it can be compensated by hydrogen already at 1% Al. We also assess the diffusivity of hydrogen and study complex formation. Si cation – H complexes have relatively low binding energies. Substitutional carbon on a cation site acts as a shallow donor in Ga 2 O 3 , but can be stable in a negative charge state in ( Al x Ga 1 – x ) 2 O 3 when x > 5 % . Substitutional carbon on an oxygen site ( C O ) always acts as an acceptor in n -type ( Al x Ga 1 – x ) 2 O 3 , but will incorporate only under relatively oxygen-poor conditions. C O – H complexes can actually incorporate more easily, explaining observations of carbon-related compensation in Ga 2 O 3 grown by MOCVD. We also investigate C cation – H complexes, finding they have high binding energies and act as compensating acceptors when x > 56 % ; otherwise the hydrogen just passivates the unintentional carbon donors. C-H complex formation explains why MOCVD-grown Ga 2 O 3 can exhibit record-low free-carrier concentrations, in spite of the unavoidable incorporation of carbon. Our study highlights that, while Si is in principle a suitable shallow donor in ( Al x Ga 1 – x ) 2 O 3 alloys up to high Al compositions, control of unintentional impurities is essential to avoid compensation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Critical role of magnetic moments in heavy-fermion materials: Revisiting mysteries of SmB 6

Heavy-fermion family exhibits fascinating and often puzzling properties due to the presence of open-shell f ions and the complexity of the associated charge, orbital, and spin degrees of freedom. SmB 6 is a prototypical heavy-fermion compound that is electrically insulating but yet it displays quantum oscillations, which are a telltale signature of the metallic state. Adding to the enigma is the possibility that SmB 6 is a topological Kondo insulator. Here, by treating the spin degree of freedom on an equal footing with other degrees of freedom using the parameter-free strongly-constrained and appropriately-normed (SCAN) density functional, we explore the ground-state electronic structure of SmB6. A number of competing magnetic phases lying very closely in energy are found, indicating the key role of spin fluctuations in the material. The computed band structure, crystal-field splittings in the f-electron complex, the heavy effective electron mass at the Fermi energy, and the large specific heat are all in good agreement with the corresponding experimental results. In particular, our predicted FS explains the experimentally observed bulk quantum oscillations as well as the low electrical conductivity of SmB 6 . The topological Kondo state of SmB6 is shown to be robust regardless of its magnetic configuration. The excellent performance of SCAN in heavy-fermion systems is explained in terms of its ability to treat self-interaction errors and symmetry breaking within the framework of the density functional theory. Our study provides a new approach for modeling heavy-fermion materials

36 MATERIALS SCIENCE↗
DOE OSTI2022

Trigonal symmetry breaking and its electronic effects in the two-dimensional dihalides MX 2 and trihalides MX 3

We study the consequences of the approximately trigonal (D 3d ) point symmetry of the transition metal (M) site in two-dimensional van der Waals MX 2 dihalides and MX 3 trihalides. The trigonal symmetry leads to a 2-2-1 orbital splitting of the transition metal d shell, which is best represented by d orbitals that describe the trigonal rather than Oh symmetry. The ligand-ligand bond length differences (rather than metal-ligand) take the role of a Jahn-Teller-like mode and in combination with interlayer distances and dimensionality effects tune the crystal field splittings and bandwidths. These effects, in turn, are amplified by electronic correlation effects—leading to crystal field splittings of the order of 0.1–1 eV between the singlet and the lowest orbital doublet—as opposed to the often assumed degenerate 3 t 2g states. Our calculations explain why most of the materials in this family are insulating, and why these considerations have to be taken into account in realistic models of them. Further, orbital order coupled to various lower symmetry lattice modes may lift the remaining orbital degeneracies, and we explain how these may support unique electronic states using ZrI 2 and CuCl 2 as examples, and offer a brief overview of possible electronic configurations in this class of materials. By building and analyzing Wannier models adapted to the appropriate symmetry we examine how the interplay among trigonal symmetry, electronic correlation effects, and – d orbital charge transfer leads to insulating, orbitally polarized magnetic and/or orbital-selective Mott states, and we provide a simple framework and numerical tool for others. Our paper establishes a rigorous framework to understand, control, and tune the electronic states in low-dimensional correlated halides. Furthermore, our analysis shows that trigonal symmetry and its breaking is a key feature of the two-dimensional halides that needs to be accounted for in search of novel electronic states in materials ranging from CrI 3 to α–RuCl 3 .

2-dimensional systems↗
DOE OSTI2022

Effect of realistic out-of-plane dopant potentials on the superfluid density of overdoped cuprates

Recent experimental papers on hole-doped overdoped cuprates have argued that a series of observations showing unexpected behavior in the superconducting state imply the breakdown of the quasiparticle-based Landau–BCS paradigm in that doping range. In contrast, some of the present authors have argued that a phenomenological “dirty d-wave” theoretical analysis explains essentially all aspects of thermodynamic and transport properties in the superconducting state, provided the unusual effects of weak, out-of-plane dopant impurities are properly accounted for. Here we attempt to place this theory on a more quantitative basis by performing ab initio calculations of dopant impurity potentials for LSCO and Tl-2201. These potentials are more complex than the pointlike impurity models considered previously, and require calculation of forward scattering corrections to transport properties. Including realistic, ARPES-derived band structures, Fermi liquid renormalizations, and vertex corrections, we show that the theory can explain semiquantitatively the unusual superfluid density measurements of the two most studied overdoped materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Composite fermion mass: Experimental measurements in ultrahigh quality two-dimensional electron systems

Composite fermions (CFs), exotic quasiparticles formed by pairing an electron and an even number of magnetic flux quanta, emerge at high magnetic fields in an interacting electron system, and can explain phenomena such as the fractional quantum Hall state (FQHS) and other many-body phases. CFs possess an effective mass (m CF ) whose magnitude is inversely related to the most fundamental property of a FQHS, namely its energy gap. Here we present here experimental measurements of m CF in ultrahigh quality two-dimensional electron systems confined to GaAs quantum wells of varying thickness. An advantage of measuring m CF over gap measurements is that mass values are insensitive to disorder and are therefore ideal for comparison with theoretical calculations, especially for high-order FQHS. Our data reveal that m CF increases with increasing well width, reflecting a decrease in the energy gap as the electron layer becomes thicker and the in-plane Coulomb energy softens. Comparing our measured masses with available theoretical results, we find significant quantitative discrepancies, highlighting that more rigorous and accurate calculations are needed to explain the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interplay between magnetism and band topology in the kagome magnets R Mn 6 Sn 6

Kagome-lattice magnets R Mn 6 Sn 6 recently emerged as a new platform to exploit the interplay between magnetism and topological electronic states. Some of the most exciting features of this family are the dramatic dependence of the easy magnetization direction on the rare-earth specie, despite other magnetic and electronic properties being essentially unchanged, and the kagome geometry of the Mn planes that in principle can generate flat bands and Dirac points; gapping of the Dirac points by spin-orbit coupling has been suggested recently to be responsible for the observed anomalous Hall response in the member TbMn 6 Sn 6 . In this paper, we address both issues with ab initio calculations. We have discovered the significant role played by higher-order crystal-field parameters and rare-earth magnetic anisotropy constants in these systems. Further, we demonstrate that the microscopic origin of rare-earth magnetic anisotropy can also be quantified and understood at various levels: ab initio, phenomenological, and analytical. In particular, using a simple and physically transparent analytical model based on perturbation theory, we are able to explain, with full quantitative agreement, the evolution of rare-earth magnetic anisotropy across the series. We analyze in detail the topological properties of Mn-dominated bands and demonstrate how they emerge from the multiorbital planar kagome model. We further show that, despite this fact, most of the topological features at the Brillouin zone corner K are strongly 3D and therefore cannot explain the observed quasi-2D anomalous Hall effect, while the most pronounced quasi-2D dispersion are too far removed from the Fermi level. By employing self-consistent calculations with ab initio many-body approaches, we demonstrate that the exchange-correlation effects beyond the density functional theory for itinerant Mn- d electrons do not significantly alter the obtained electronic and magnetic structure. Therefore, we conclude that, contrary to previous claims, the most pronounced 2D kagome-derived topological band features bear little relevance to transport in RMn 6 Sn 6 , albeit they may possibly be brought to focus by electron or hole doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Symmetry breaking forms split-off flat bands in quantum oxides controlling metal versus insulator phases

The crystal structure used as input to electronic structure calculation of conventionally bonded solids consists primarily of the standard crystallographic degrees of freedom. Quantum solids sometimes have additional microscopic degrees of freedom (m-DOF) nested within the crystallographic structure. These might include local motifs including positional (Peierls dimers, deformed octahedra, Jahn-Teller distortions), magnetic (local moment configurations) and dipolar (local ferroelectric configurations). Such motifs can be observed experimentally via local probes that avoid averaging, and theoretically as distinct total energy lowering features relative to more simplified “average crystallographic structures”. Here we examine the ability of electronic structure methods independent of strong correlation physics to explain the broad phenomenology of metal-insulator selectivity in quantum oxides by energy-lowering symmetry breaking. We do this by avoiding the restriction of considering only (i) electron-electron interactions in a fixed unresponsive lattice—as done in Mott strong correlation explanations—allowing, however, (ii) local positional and magnetic motifs that coexist within a crystallographic landscape. We find in a broad range of quantum oxides with different symmetry breaking modes the formation of split-off flat bands that can also control metallic versus insulating characteristics. The split-off band effect is common to magnetic and non-magnetic materials, including both binary and ternary oxides. One finds that when such local symmetry breaking motifs are considered in mean-field-like (e.g. density functional theory) electronic structure calculations of quantum oxides, they explain many of the observed trends in metal vs insulator phases discussed otherwise in prior literature via strong correlation in symmetry unbroken view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Role of magnetic and structural symmetry breaking in forming the Mott insulating gap in Nb 3 ⁢Cl 8

The α-phase of the gapped insulator Nb 3 Cl 8 has recently emerged as the long-sought critical testing bed for examining the importance of strong interelectronic correlation vs symmetry breaking in understanding insulation of such Mott compounds. Structural symmetry breaking detected by density functional theory (DFT) energy lowering (such as dimer formation, disproportionation, or Jahn-Teller distortions) explains insulation in both d-electron Mott-like systems and in non-d-electron cases without recourse to strong correlation. Yet, in Nb 3 Cl 8 , structural symmetry breaking alone (viz. formation of Nb trimers) fails to explain insulation, leading instead to a partially occupied metallic flat band, in contrast with experimental observations. We examine the role of magnetic symmetry breaking, noting that Nb 3 Cl 8 is an observed paramagnet (not an antiferromagnet), thus potentially carrying also short-range ordered magnetic moments. Describing the latter as a polymorphous distribution of nonzero local moments with total zero net magnetization is demonstrated to lower the DFT total energy, while gapping the system without recourse to strong correlation or long-range magnetic order. This suggests that degeneracy removal by symmetry breaking in mean-field-like approaches—either structural, or magnetic, or both—can reduce or eliminate the need for strong correlation, allowing the use of DFT for such Mott systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Tensor force role in β decays analyzed within the Gogny-interaction shell model

The half-life of the famous C 14 β decay is anomalously long, with different mechanisms: the tensor force, cross-shell mixing, and three-body forces, proposed to explain the cancellations that lead to a small transition matrix element. In this study, we revisit and analyze the role of the tensor force for the β decay of C 14 as well as of neighboring isotopes. We add a tensor force to the Gogny interaction, and derive an effective Hamiltonian for shell-model calculations. The calculations were carried out in a p – s d model space to investigate cross-shell effects. Furthermore, we decompose the wave functions according to the total orbital angular momentum L in order to analyze the effects of the tensor force and cross-shell mixing. The inclusion of the tensor force significantly improves the shell-model calculations of the β -decay properties of carbon isotopes. In particular, the anomalously slow β decay of C 14 can be explained by the isospin T = 0 part of the tensor force, which changes the components of N 14 with the orbital angular momentum L = 0 , 1 , and results in a dramatic suppression of the Gamow-Teller transition strength. At the same time, the description of other nearby β decays are improved. Decomposition of wave function into L components illuminates how the tensor force modifies nuclear wave functions, in particular suppression of β -decay matrix elements. Cross-shell mixing also has a visible impact on the β -decay strength. Inclusion of the tensor force does not seem to significantly change, however, binding energies of the nuclei within the phenomenological interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Probing quarkyonic matter in neutron stars with the Bayesian nuclear-physics multimessenger astrophysics framework

The interiors of neutron stars contain matter at the highest densities realized in our Universe. Interestingly, theoretical studies of dense matter, in combination with the existence of two-solar-mass neutron stars, indicate that the speed of sound $c_s$ has to increase to values well above the conformal limit ($c_s^2$ = 1/3) before decreasing again at higher densities. Further, the decrease could be explained by either a strong first-order phase transition or a crossover transition from hadronic to quark matter. The latter scenario leads to a pronounced peak in the speed of sound, reaching values above the conformal limit, naturally explaining the inferred behavior. In this work, we use the nuclear-physics multimessenger astrophysics (NMMA) framework to compare predictions of the quarkyonic matter model with astrophysical observations of neutron stars, with the goal of constraining model parameters. Assuming quarkyonic matter to be realized within neutron stars, we find that there can be a significant amount of quarks inside the cores of neutron stars with masses in the two-solar-mass range, amounting to up to ≈0.13$M$ ⊙ , contributing ≈ 5.9% of the total mass. Furthermore, for the quarkyonic matter model investigated here, the radius of a 1.4$M$ ⊙ neutron star would be $13.44_{–1.54}^{+1.69}(13. 54_{–1.04}^{+1.02})$ km, at 95% credibility, without (with) the inclusion of AT2017gfo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗