Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex quantum materials”

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

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Collective Phonon–Polaritonic Modes in Silicon Carbide Subarrays

Localized surface phonon polaritons (LSPhPs) can be implemented to engineer light–matter interactions through nanoscale patterning for a range of midinfrared application spaces. However, the polar material systems studied to date have mainly focused on simple designs featuring a single element in the periodic unit cell. Increasing the complexity of the unit cell can serve to modify the resonant near-fields and intra- and inter-unit-cell coupling as well as to dictate spectral tuning in the far-field. In this work, we exploit more complicated unit-cell structures to realize LSPhP modes with additional degrees of design freedom, which are largely unexplored. Collectively excited LSPhP modes with distinctly symmetric and antisymmetric near-fields are supported in these subarray designs, which are based on nanopillars that are scaled by the number of subarray elements to ensure a constant unit-cell size. Moreover, we observe an anomalous mode-matching of the collective symmetric mode in our fabricated subarrays that is robust to changing numbers of pillars within the subarrays as well as to defects intentionally introduced in the form of missing pillars. This work therefore illustrates the hierarchical design of tailored LSPhP resonances and modal near-field profiles simultaneously for a variety of IR applications such as surface-enhanced spectroscopies and biochemical sensing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Spin-orbital liquids and insulator-metal transitions on the pyrochlore lattice

The two orbital Hubbard model, with the electrons additionally coupled to a complex magnetic background, arises in the pyrochlore molybdates. The background involves local moments Hund's coupled to the electrons, driving double exchange ferromagnetism, and antiferromagnetic (AF) tendency arising from competing superexchange. The key scales include the Hubbard repulsion and the superexchange, both of which can be tuned in these materials. They control the phase transition from a ferromagnetic metal to a spin glass metal and then to a spin glass (Mott) insulator. We provide a comprehensive description of the ground state of this model using an unrestricted Hartree-Fock scheme implemented via a simulated annealing procedure and establish the metal-insulator transition line for varying Hubbard interaction and superexchange. The electrons see an effective disorder, due to orbital frustration, already in the ferromagnetic phase. The disorder is further enhanced by antiferromagnetic coupling and the resulting magnetic disorder. As a result, increasing AF coupling shifts the metal-insulator transition to lower Hubbard interaction and gives it an additional "Anderson" character. In conclusion, we provide detailed results on the magnetic and orbital correlations, the density of states, and the optical conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Magnetic anisotropy and spin scattering in (La 2/3 Sr 1/3 )MnO 3 /CaRuO 3 bilayers

Complex ferromagnetic oxides have been identified as possible candidate materials for sources of spin currents. Here we study bilayers of ferromagnetic (La 2/3 Sr 1/3 )MnO 3 (LSMO) and metallic CaRuO 3 (CRO) on LSAT substrates as a model system for spin pumping. Ferromagnetic resonance (FMR) measurements of these bilayers show evidence of spin pumping across the interface in the form of an increase in Gilbert damping with the addition of CRO. FMR indicates that the presence of CRO modifies the magnetic anisotropy of the LSMO. By increasing CRO thickness, we find a reduction of the out-of-plane anisotropy and simultaneous rotation of the easy axis within the plane, from the $\langle$110$\rangle$ to $\langle$100$\rangle$ axis. The evolution of magnetic anisotropy determined by FMR disagrees with that measured by bulk SQUID magnetometry and is accompanied by structural distortions in the LSMO layer as measured by x-ray diffraction, thus suggesting a change in magnetic anisotropy attributed to structural distortions imposed on LSMO by CRO. These results suggest that while LSMO and CRO remain promising candidates for efficient pure spin current generation and detection, respectively, epitaxial integration of perovskites will cause additional changes which must be accounted for in spintronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Correlating Mechanical Sensitivity with Spin Transition in the Explosive Spin Crossover Complex [Fe(Htrz) 3 ] n [ClO 4 ] 2 n

Spin crossover complexes are known to undergo bond length, volume, and enthalpy changes during spin transition. In an explosive spin crossover complex, these changes could affect the mechanical and initiation sensitivity of the explosive and lead to the development of a new class of sensitivity switchable materials. To explore this relationship, the well-known spin crossover compound [Fe(Htrz) 3 ] n [ClO 4 ] 2n (1) was re-evaluated for its explosive properties, and its mechanical impact sensitivity was correlated to spin transition. A variable temperature impact test was developed and used to evaluate the impact sensitivity of 1 in the low spin (LS, S = 0), thermally accessed high spin (HS, S = 2), and mixed LS and HS states. For comparison, the structurally similar Ni compound, [Ni(Htrz) 3 ] n [ClO 4 ] 2n (2), which does not undergo a spin transition at accessible temperatures, was synthesized and characterized, and its explosive properties and variable temperature impact sensitivity measured. These results reveal a correlation between impact sensitivity and spin transition, where 1 exhibits lower impact sensitivity in the LS state and increases in sensitivity upon transition to the HS state. Here, density functional theory was used to predict structural changes that occur upon spin transition that correlate to the change in sensitivity. This demonstrates, for the first time, an explosive spin crossover compound (ExSCO) that exhibits switchable impact sensitivity with a fully reversible internal switching mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Dual-Wavelength Simultaneous Patterning of Degradable Thermoset Supports for One-Pot Embedded 3D Printing

Vat photopolymerization (VP) techniques have enabled the fabrication of complex geometries while balancing high precision and fast processing times. 3D printed objects are traditionally built layer-by-layer with newly cured layers being structurally supported by previous ones. Fabricating unsupported features such as overhangs and arches risks misalignment and sagging, limiting the range of accessible designs. To overcome this issue, support structures are fabricated along with the primary object as temporary scaffolds that provide stability and conserve print fidelity. For VP specifically, patterning dissolvable sacrificial supports is attractive to avoid manual removal after printing. In this study, we demonstrate a base-degradable thermoset to pattern print supports in a one-pot formulation along with the primary structural material. Efficient printing is enabled using a dual-wavelength negative imaging (DWNI) DLP printer that patterns the degradable thermoset with visible light and the permanent network with UV light, which are simultaneously projected using a single digital micromirror device (DMD). Printed objects undergo thermal postprocessing to enhance the final conversion of the primary material, after which thermoset supports are degraded in a basic, aqueous solution. This approach provides a robust method for the dual-wavelength patterning of sacrificial thermoset supports, broadening the range of accessible 3D printable materials and geometries.

3D printing↗
DOE OSTI2020

Feasibility of Realizing a Universal Identifier Based on Fiber Bragg Grating Technology

Over the past year we have collected data on different fiber Bragg gratings (FBG) with the end goal of understanding their suitability for use in a new type of unique identifier (UID) that could be used in both tags and seals for high-value, high-risk assets. The ideal UID will have a complex signature that allows items to be uniquely identified by that signature. Tamper indication must be an integral part of the tag/seal process with measurable changes to the device signature indicative of any attempt to tamper with the monitored item. Further, the signature must be so difficult to reproduce that even state-actors with significant resources at their disposal will be discouraged from attempting it. Finally, readout verification of the product must be simple and secure, while providing red-light/green-light results in real time. To evaluate how well an FBG-based UID could meet these goals, we collected data on a large number of FBG samples from different vendors. Measurements included the response of the devices to uniform tension, tension when wrapped on a mandrel, and the response to changes in temperature and humidity, as well as the effect of simple potting in expanding and nonexpanding potting materials. We tested FBG samples written into single-mode fibers that were coated with the two commonly used materials, polyimide and acrylate and gratings that were both uniform and those with an apodized response (i.e., the intensity of the grating was varied to minimize sidelobes). The result of the work to date did not find any significant issues with the design concept for a UID based on this technology. Specific findings include a baseline design comprising five FBGs with distinct spectral returns analyzed using a difference-correlation approach. To realize this device, careful attention must be given to the impact of environmental factors on the UID response, limiting the properties of the potting compounds and FBG coating materials that can be used.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Tunable correlated Chern insulator and ferromagnetism in a moiré superlattice

Studies of two-dimensional electron systems in a strong magnetic field revealed the quantum Hall effect, a topological state of matter featuring a finite Chern number C and chiral edge states. Haldane later theorized that Chern insulators with integer quantum Hall effects could appear in lattice models with complex hopping parameters even at zero magnetic field. The ABC-trilayer graphene/hexagonal boron nitride (ABC-TLG/hBN) moiré superlattice provides an attractive platform with which to explore Chern insulators because it features nearly flat moiré minibands with a valley-dependent, electrically tunable Chern number. In this work, we report the experimental observation of a correlated Chern insulator in an ABC-TLG/hBN moiré superlattice. We show that reversing the direction of the applied vertical electric field switches the moiré minibands of ABC-TLG/hBN between zero and finite Chern numbers, as revealed by large changes in magneto-transport behaviour. For topological hole minibands tuned to have a finite Chern number, we focus on quarter filling, corresponding to one hole per moiré unit cell. The Hall resistance is well quantized at h/2e 2 (where h is Planck’s constant and e is the charge on the electron), which implies C = 2, for a magnetic field exceeding 0.4 tesla. The correlated Chern insulator is ferromagnetic, exhibiting substantial magnetic hysteresis and a large anomalous Hall signal at zero magnetic field. Finally, our discovery of a C = 2 Chern insulator at zero magnetic field should open up opportunities for discovering correlated topological states, possibly with topological excitations, in nearly flat and topologically nontrivial moiré minibands.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental investigation of laser ablation of stone polycrystalline targets

We report the results of an experimental investigation of ablation of stone polycrystalline targets of complex multicomponent composition, which imitate the substance of asteroids. The targets were irradiated by nanosecond pulses of a neodymium laser at an energy density Φ{sub L} of up to 5 × 10{sup 4} J cm{sup −2}. The experiments demonstrated the existence of two ablation regimes, with the boundary between them lying at Φ{sub L} ≈ 4000 J cm{sup −2}. The regime change is characterised by a change in the form of the dependence of the surface mass density of removed target material on the laser energy density and by the appearance of a minimum in the dependence of specific energy of destruction on Φ{sub L}. This is supposedly related to the passage from a one-dimensional plasma plume expansion to the three-dimensional one and the corresponding decrease in the efficiency of energy transfer from the laser beam to the target due to a lowering of laser-produced plasma density. Our experiments also showed the existence of a maximum in impulse coupling coefficient C{sub m} as a function of laser energy density (C{sub m} ≈ 6.3 × 10{sup −5} N W{sup −1} for Φ{sub L} = Φ{sub opt} ≈ 25 J cm{sup −2}). Maxima were also recorded in the dependences of the ablation efficiency and average ablation flow velocity on Φ{sub L}. For Φ{sub L} > Φ{sub opt}, the decrease in the function C{sub m}(Φ{sub L}) turns out to be much steeper than for metals and polymer materials. The difference is supposedly due to the lower strength and lower plasticity of the polycrystalline stone targets. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Colossal Strain Tuning of Ferroelectric Transitions in KNbO 3 Thin Films

Strong coupling between polarization ( P ) and strain (ɛ) in ferroelectric complex oxides offers unique opportunities to dramatically tune their properties. Here colossal strain tuning of ferroelectricity in epitaxial KNbO 3 thin films grown by sub‐oxide molecular beam epitaxy is demonstrated. While bulk KNbO 3 exhibits three ferroelectric transitions and a Curie temperature ( T c ) of ≈676 K, phase‐field modeling predicts that a biaxial strain of as little as −0.6% pushes its T c > 975 K, its decomposition temperature in air, and for −1.4% strain, to T c > 1325 K, its melting point. Furthermore, a strain of −1.5% can stabilize a single phase throughout the entire temperature range of its stability. A combination of temperature‐dependent second harmonic generation measurements, synchrotron‐based X‐ray reciprocal space mapping, ferroelectric measurements, and transmission electron microscopy reveal a single tetragonal phase from 10 K to 975 K, an enhancement of ≈46% in the tetragonal phase remanent polarization ( P r ), and a ≈200% enhancement in its optical second harmonic generation coefficients over bulk values. These properties in a lead‐free system, but with properties comparable or superior to lead‐based systems, make it an attractive candidate for applications ranging from high‐temperature ferroelectric memory to cryogenic temperature quantum computing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Relative permeability as a stationary process: Energy fluctuations in immiscible displacement

Relative permeability is commonly used to model immiscible fluid flow through porous materials. In this work, we derive the relative permeability relationship from conservation of energy, assuming that the system to be non-ergodic at large length scales and relying on averaging in both space and time to homogenize the behavior. Explicit criteria are obtained to define stationary conditions: (1) there can be no net change for extensive measures of the system state over the time averaging interval; (2) the net energy inputs into the system are zero, meaning that the net rate of work done on the system must balance with the heat removed; and (3) there is no net work performed due to the contribution of internal energy fluctuations. Results are then evaluated based on direct numerical simulation. Dynamic connectivity is observed during steady-state flow, which is quantitatively assessed based the Euler characteristic. We show that even during steady-state flow at low capillary number ([Formula: see text]), typical flow processes will explore multiple connectivity states. The residence time for each connectivity state is captured based on the time-and-space average. The distribution for energy fluctuations is shown to be multi-modal and non-Gaussian when terms are considered independently. However, we demonstrate that their sum is zero. Given an appropriate choice of the thermodynamic driving force, we show that the conventional relative permeability relationship is sufficient to model the energy dissipation in systems with complex pore-scale dynamics that routinely alter the structure of fluid connected pathways.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Real-space observation of ergodicity transitions in artificial spin ice

Abstract Ever since its introduction by Ludwig Boltzmann, the ergodic hypothesis became a cornerstone analytical concept of equilibrium thermodynamics and complex dynamic processes. Examples of its relevance range from modeling decision-making processes in brain science to economic predictions. In condensed matter physics, ergodicity remains a concept largely investigated via theoretical and computational models. Here, we demonstrate the direct real-space observation of ergodicity transitions in a vertex-frustrated artificial spin ice. Using synchrotron-based photoemission electron microscopy we record thermally-driven moment fluctuations as a function of temperature, allowing us to directly observe transitions between ergodicity-breaking dynamics to system freezing, standing in contrast to simple trends observed for the temperature-dependent vertex populations, all while the entropy features arise as a function of temperature. These results highlight how a geometrically frustrated system, with thermodynamics strictly adhering to local ice-rule constraints, runs back-and-forth through periods of ergodicity-breaking dynamics. Ergodicity breaking and the emergence of memory is important for emergent computation, particularly in physical reservoir computing. Our work serves as further evidence of how fundamental laws of thermodynamics can be experimentally explored via real-space imaging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Improved loss functions for machine-learned atomic potentials

Machine learning (ML) has become an invaluable tool across a wide array of domains in science as researchers find new ways to leverage its predictive power. This is especially true in chemistry, where ML is used to fit chemical properties or desirable attributes to the local structure of molecules and materials. In the pursuit of greater accuracy, it is relatively simple to increase the size or complexity of such models, although this often requires simultaneously seeking larger datasets in order to both fit and interpret the larger number of parameters. However, it is equally important to assess the quality and relative importance of the data and how these factors impact the training process. We, therefore, investigate the impact of using different loss functions for training neural network potentials (NNPs), as the loss function defines the error and parameter gradients used to train the NNP. In particular, we test the mean-squared error and Huber loss functions and, using insight from these functions, derive a new loss function based on the Asinh function, which yields significant improvement in the accuracy and generality of NNPs. We show that by discounting/minimizing errors and anomalies in the optimization process, both the Huber and Asinh loss functions improve the training of NNPs, leading to a final potential with a greater effective dimensionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A Novel Spatio-Temporal Regime Tracking Method for Impact Simulations

In this proposal, we present a novel method of tracking the rheological regimes activated during impact cratering events that will allow researchers to gain new insights into cratering mechanics. Rheology describes the stress-strain response of rocks to different conditions. Planetary impact cratering events often occur on too large of a scale to be feasibly captured in controlled experiments. Instead, these dynamic events are primarily studied using multi-physics codes equipped with complex material models that enable calculations of the impact event at scale. However, determining which physical processes are required for the problem of interest is challenging. Because the dominant rheological regimes change with space and time during crater formation, it is difficult to link numerical simulations with observable features of craters at the end of the event. The basis of this work is the implementation of numerical flags that track the activation of each rheological regime throughout impact simulations. We demonstrate this with the ‘Rock Model’ implemented in the CTH shock-physics code. We will use this rheology tracking method to ’zoom in’ on a specific region within an event and track the conditions the rock experiences over time. This work will develop community benchmarks to validate and distribute our implemented rheological models. We will focus on improving the melt models used by the planetary impact modeling community by developing an EOS-aware rheological transition from solid to melt. Through analysis of cell and tracer-particle based tracking data, we will study the effects of different rheologies on modeled outcomes, particularly on the volume and distributions of impacts melts. We will also use this method to link observable features with the rheological mechanisms responsible for them. The deliverables (peer-reviewed papers) from the proposed work are (i) tests of the implemented rheologic processes and demonstrations of the tracking flags; (ii) the first calculations of the spatio-temporal evolution of the dominant rheologies during impact cratering events; and (iii) application to delivery of impactor iron during basin-scale impacts.

58 GEOSCIENCES↗
DOE OSTI2022

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrically controlled interlayer trion fluid in electron-hole bilayers

Here, the combination of repulsive and attractive Coulomb interactions in a quantum electron-hole (e-h) fluid can produce correlated phases of multiparticle charge complexes, such as excitons, trions, and biexcitons. We report an experimental realization of an electrically controlled interlayer trion fluid in van der Waals heterostructures. In strongly coupled e-h bilayers, electrons and holes spontaneously form three-particle trion bound states. The interlayer trions can assume 1e-2h and 2e-1h configurations. We show that the two holes in 1e-2h trions form a spin-singlet with a spin gap of approximately one milli–electron volt. By electrostatic gating, the equilibrium state can be continuously tuned into an exciton fluid, a trion fluid, an exciton-trion mixture, or a trion-charge mixture. Our work demonstrates a platform to study correlated phases of tunable Bose-Fermi mixtures.

Qi, Ruishi [University of California, Berkeley, CA↗
DOE OSTI2022

Strongly Correlated Quantum Spin Liquids versus Heavy Fermion Metals: A Review

This review considers the topological fermion condensation quantum phase transition (FCQPT) that explains the complex behavior of strongly correlated Fermi systems, such as frustrated insulators with quantum spin liquid and heavy fermion metals. The review contrasts theoretical consideration with recent experimental data collected on both heavy fermion metals (HF) and frustrated insulators. Such a method allows to understand experimental data. We also consider experimental data collected on quantum spin liquid in Lu 3 Cu 2 Sb 3 O 14 and quasi-one dimensional (1D) quantum spin liquid in both YbAlO 3 and Cu(C 4 H 4 N 2 )(NO 3 ) 2 with the aim to establish a sound theoretical explanation for the observed scaling laws, Landau Fermi liquid (LFL) and non-Fermi-liquid (NFL) behavior exhibited by these frustrated insulators. The recent experimental data on the heavy-fermion metal α-YbAl 1-x Fe x B 4 , with x=0.014, and on its sister compounds β-YbAlB 4 and YbCo 2 Ge 4 , carried out under the application of magnetic field as a control parameter are analyzed. We show that the thermodynamic and transport properties as well as the empirical scaling laws follow from the fermion condensation theory. We explain how both the similarity and the difference in the thermodynamic and transport properties of α-YbAl 1-x Fe x B 4 and in its sister compounds β-YbAlB 4 and YbCo 2 Ge 4 emerge, as well as establish connection of these (HF) metals with insulators Lu 3 Cu 2 Sb 3 O 14 , Cu(C 4 H 4 N 2 )(NO 3 ) 2 and YbAlO 3 . We demonstrate that the universal LFL and NFL behavior emerge because the HF compounds and the frustrated insulators are located near the topological FCQPT or are driven by the application of magnetic fields.

36 MATERIALS SCIENCE↗