Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “condensed phase systems”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experimental realization of classical Z 2 spin liquids in a programmable quantum device

We build and probe a Z 2 spin liquid in a programmable quantum device, the D-Wave DW-2000Q. Specifically, we observe the classical eight-vertex and six-vertex (spin ice) states and transitions between them. To realize this state of matter, we design a Hamiltonian with combinatorial gauge symmetry using only pairwise-qubit interactions and a transverse field, i.e., interactions which are accessible in this quantum device. The combinatorial gauge symmetry remains exact along the full quantum annealing path, landing the system onto the classical eight-vertex model at the endpoint of the path. The output configurations from the device allow us to directly observe the loop structure of the classical model. Moreover, we deform the Hamiltonian so as to vary the weights of the eight vertices and show that we can selectively attain the classical six-vertex (ice) model, or drive the system into a ferromagnetic state. Additionally, we present studies of the classical phase diagram of the system as a function of the eight-vertex deformations and effective temperature, which we control by varying the relative strengths of the programmable couplings, and we show that the experimental results are consistent with theoretical analysis. Finally, we identify additional capabilities that, if added to these devices, would allow us to realize Z 2 quantum spin liquids on which to build topological qubits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A machine learning inversion scheme for determining interaction from scattering

Small angle scattering techniques have now been routinely used to quantitatively determine the potential of mean force in colloidal suspensions. However the numerical accuracy of data interpretation is often compounded by the approximations adopted by liquid state analytical theories. To circumvent this long standing issue, here we outline a machine learning strategy for determining the effective interaction in the condensed phases of matter using scattering. Via a case study of colloidal suspensions, we show that the effective potential can be probabilistically inferred from the scattering spectra without any restriction imposed by model assumptions. Comparisons to existing parametric approaches demonstrate the superior performance of this method in accuracy, efficiency, and applicability. This method can effectively enable quantification of interaction in highly correlated systems using scattering and diffraction experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles

The molecular dipole moment ( μ ) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly—via the ground state electron density—from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). Here, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression framework which assigns a (vector) dipole moment to each atom, while the movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting “MuML” models are fitted together to reproduce molecular μ computed using high-level coupled-cluster theory and density functional theory (DFT) on the QM7b dataset, achieving more accurate results due to the physics-based combination of these complementary terms. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models—and the optimal weighting of their combination—depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to μ ; furthermore, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.

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

Sachdev-Ye-Kitaev models and beyond: Window into non-Fermi liquids

In this work is a review of the Sachdev-Ye-Kitaev (SYK) model of compressible quantum many-body systems without quasiparticle excitations, and its connections to various theoretical studies of non-Fermi liquids in condensed matter physics. The review is placed in the context of numerous experimental observations on correlated electron materials. Strong correlations in metals are often associated with their proximity to a Mott transition to an insulator created by the local Coulomb repulsion between the electrons. The phase diagrams of a number of models of such a local electronic correlation are explored, employing a dynamical mean-field theory in the presence of random spin exchange interactions. Numerical analyses and analytical solutions, using renormalization group methods and expansions in large spin degeneracy, lead to critical regions that display SYK physics. The models studied include the single-band Hubbard model, the t–J model, and the two-band Kondo-Heisenberg model in the presence of random spin exchange interactions. Also examined are non-Fermi liquids obtained by considering each SYK model with random four-fermion interactions to be a multiorbital atom, with the SYK atoms arranged in an infinite lattice. Connections are made to theories of sharp Fermi surfaces without any low-energy quasiparticles in the absence of spatial disorder, obtained by coupling a Fermi liquid to a gapless boson; a systematic large-N theory of such a critical Fermi surface, with SYK characteristics, is obtained by averaging over an ensemble of theories with random boson-fermion couplings. Finally, an overview of the links between the SYK model and quantum gravity is presented, and the review ends with an outlook on open questions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The phase diagram of ultra quantum liquids

In this paper, we discuss the dependence of the phase diagram of a hypothetical isotope of helium with nuclear mass less than 4 atomic mass units. We argue that with decreasing nucleus mass, the temperature of the superfluid phase transition (about 2.2 K in real 4 He) increases, while that of the liquid-gas critical point (about 5.2 K in real 4 He) decreases. We discuss various scenarios that may occur when the two temperatures approach each other and the order parameters of the superfluid and the liquid-gas phase transitions interact with each other. The simplest scenario, in which both order parameters become critical at particular values of the nuclear mass, temperature, and pressure, can be ruled out through on an analysis of the Landau theory. We argue that in the most likely scenario, as the nuclear mass decreases, first, a tricritical point appears on the line separating the superfluid and the normal fluid phase, then the critical point disappears under the first-order part of superfluid phase transition line, and in the end the tricritical point disappears. The last change in the phase diagram occurs when the two-body scattering length crosses zero, which corresponds to the nuclear mass of about 1.55 u. We develop a quantitative theory that allows one to determine the phase diagram in the vicinity of this point. Lastly, we discuss several ways to physically realize such liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Stacking faults and alternate crystal structures for the shape-memory alloy NiTi

Here, we use ab initio calculations to study the role of stacking faults in connecting the high-temperature B2 and the theoretically predicted low-temperature B33 NiTi phases. In contrast with prior work, we describe the B 2 → B 33 phase transformation in terms of alternate bilayer shifts by 1 2 [100] on the ( 011 ) B 2 plane, obtaining a viable pathway; the same mechanism could also work with the B19 parent phase. We then examine B33-like structures built from alternate stacking sequences of B19 bilayers, constructed to have monoclinic tilt angles close to the experimentally reported NiTi B19' martensite, and find four low-energy stacking-fault variants with energies 5.8–8.5 meV/atom above the calculated B19' martensite structure, suggesting that such structures might appear as a part of the NiTi martensite phase at low temperatures. Investigating further the occurrence of specific coordinated planar shifts in NiTi systems, we report a dynamically stable NiTi B27 phase and find that it is only 1.2 meV/atom above the calculated B33 ground-state structure, thus having a potential to also play a role in NiTi martensitic phase transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Correlation strength and orbital differentiation across the phase diagram of plutonium metal

We compare the trends on the strength of electronic correlations across the different phases of elemental Pu focusing on its site and orbital dependence, using a combination of density functional theory (DFT) and dynamical mean-field theory (DMFT) calculations within the vertex corrected one crossing approximation. We find that Pu-5f states are more correlated in δ -Pu, followed by some crystallographic sites in α and β phases. In addition, we observe that Pu-5f 5/2 and Pu-5f 7/2 orbital differentiation is a general feature of this material, as is site differentiation in the low-symmetry phases. The Pu-5f 5/2 states show Fermi liquid like behavior, whereas the Pu-5f 7/2 states remaining incoherent down to very low temperatures. We correlate the correlation strength in the different phases to their structure and the Pu-5 f occupancy of their crystallographic sites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Revealing the competition between charge density wave and superconductivity in CsV 3 Sb 5 through uniaxial strain

In this paper we report the effect of uniaxial strain ε applied along the crystalline a axis on the newly discovered kagome superconductor CsV 3 Sb 5 . At ambient conditions, CsV 3 Sb 5 shows a charge-density wave (CDW) transition at T CDW = 94.5K and superconducts below T c = 3.34K. In our paper, when the uniaxial strain ε is varied from –0.90% to 0.90%,T c monotonically increases by ~33% from 3.0 to 4.0 K, giving rise to the empirical relation T c (ε) = 3.4 + 0.56ε + 0.12ε 2 . On the other hand, for ε changing from –0.76% to 1.26%,TCDW decreases monotonically by ~10% from 97.5 to 87.5 K with T CDW (ε) = 94.5 – 4.72ε – 0.60ε 2 . The opposite response of T c and T CDW to the uniaxial strain suggests strong competition between these two orders. Comparison with hydrostatic pressure measurements indicate that it is the change in the c axis that is responsible for these behaviors of the CDW and superconducting transitions, and that the explicit breaking of the sixfold rotational symmetry by strain has a negligible effect. Combined with our first-principles calculations and phenomenological analysis, we conclude that the enhancement in Tc with decreasing c is caused primarily by the suppression of T CDW , rather than strain-induced modifications in the bare superconducting parameters. We propose that the sensitivity of T CDW with respect to the changes in the c axis arises from the impact of the latter on the trilinear coupling between the $M$$^{+}_{1}$ and the $L$$^{–}_{2}$ phonon modes associated with the CDW. Altogether, our paper reveals that the c-axis lattice parameter, which can be controlled by both pressure and uniaxial strain, is a powerful tuning knob for the phase diagram of CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantification of Energy Savings and Demand Reduction for a Heat Pump Integrated with Thermal Energy Storage (Final Report)

To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries or installing new peaking power plants. Stakeholders and policy makers across the United States have expressed interests in promoting TES, as demonstrated by the US Department of Energy’s Grid-Interactive Efficient Buildings program and the efforts of various state legislatures. However, the cost value provided by TES are unclear. If reliable cost benefits were determined, stakeholders would have a clearer picture of the financial returns that can be gained from their investment in TES. The study in this report is conducted by ORNL with collaboration with Emerson the Helix Innovation Center. In the first part of this report, EnergyPlus was used to perform whole-building simulations for two residential buildings in Indianapolis and Atlanta. The HVAC system in both buildings were equipped with phase change material TES. The TES tank was charged in off-peak hours and discharged in peak hours to perform load shifting. First, the economic value implied by existing time-of-use (TOU) rates offered by utility companies was analyzed via whole-building simulation. Second, existing demand reduction (DR) incentives sourced from 3 different electrical grid administrators (i.e., California, Texas, and New England region) were surveyed to determine their implied value. The study suggests that the traditional value analysis that focuses on ROI for the building owner significantly undervalues TES technology making economic viability difficult. A more comprehensive value analysis that includes peak demand management and deferred capital for peaking power plants shows that TES should be economically viable but here the value is greater for the utility and requires large market penetration and aggregation to fully realize the benefits. Therefore, to facilitate commercialization, new business models are needed that include a broader range of stakeholders and distribute the value of TES proportionally. In the 2 nd part of this project, the benefits of a novel phase change material (PCM) integrated heat pump configuration were evaluated via detailed component based simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation. To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.

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

Thermodynamic signatures of an antiferromagnetic quantum critical point inside a superconducting dome

Recent experiments in unconventional superconductors, and in particular iron-based materials, have reported evidence of an antiferromagnetic quantum critical point (AFM-QCP) emerging inside the superconducting dome of the phase diagram. Fluctuations associated with such an AFM-QCP are expected to promote unusual temperature dependencies of thermodynamic quantities. Here, we compute the T dependence of the specific heat C(T) deep inside a fully gapped s +– superconducting state as the AFM-QCP is approached. We find that at the AFM-QCP, the specific heat C(T) vanishes quadratically with temperature, as opposed to the typical exponential suppression seen in fully gapped BCS superconductors. This robust result is due to a nonanalytic contribution to the free energy arising from the general form of the bosonic (AFM) propagator in the SC state. Away from the AFM-QCP, as temperature is lowered, C(T) shows a crossover from a T2 behavior to an exponential behavior, with the crossover temperature scale set by the value of the superconducting gap and the distance to the QCP. As a result, we argue that these features in the specific heat can be used to unambiguously determine the existence of AFM-QCPs inside the superconducting domes of iron-based and other fully gapped unconventional superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multi-messenger heavy-ion collision physics

We report that this work studies the production of direct photons in relativistic nuclear collisions, along with the production of hadrons. Radiation from the very first instants to the final moments of the evolution is included. The hybrid model used here describes all stages of relativistic heavy-ion collisions. Chronologically, those are an initial state reflecting the collision of nuclei described within the Color Glass Condensate effective theory; a pre-equilibrium phase based on non-equilibrium linear response; relativistic viscous hydrodynamics, and a hadronic afterburner. The effect of the preequilibrium phase on both photonic and hadronic observables is highlighted for the first time. The potential of photon observables - spectrum, differential elliptic and triangular flow - to reveal the chemical equilibration time is studied. Finally, we consider "small collision systems", including proton+nucleus collisions and collisions of light nuclei, as probed by hadronic and electromagnetic observables. We demonstrate how photon production can signal the formation of quark-gluon plasma in such small systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

SPLIT: Stable Protein Coacervation Using a Light Induced Transition

Protein coacervates serve as hubs to concentrate and sequester proteins and nucleotides and thus function as membraneless organelles to manipulate cell physiology. We have engineered a coacervating protein to create tunable, synthetic membraneless organelles that assemble in response to a single pulse of light. Coacervation is driven by the intrinsically disordered RGG domain from the protein LAF-1, and opto-responsiveness is coded by the protein PhoCl, which cleaves in response to 405 nm light. We developed a fusion protein containing a solubilizing maltose-binding protein domain, PhoCl, and two copies of the RGG domain. Several seconds of illumination at 405 nm is sufficient to cleave PhoCl, removing the solubilization domain and enabling RGG-driven coacervation within minutes in cellularsized water-in-oil emulsions. An optimized version of this system displayed light-induced coacervation in Saccharomyces cerevisiae. Finally, the methods described here provide novel strategies for inducing protein phase separation using light.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Topological obstructed atomic limit insulators by annihilating Dirac fermions

We show that annihilating a pair of Dirac fermions implies a topological transition from the critical semimetallic phase to an obstructed atomic limit insulator phase instead of a trivial insulator. This is shown to happen because of branch cuts in the phase of the wave functions, leading to nontrivial Zak phase along certain directions. To this end, we study their ${\mathbb{Z}}_{2}$ invariant and also study the phase transition using entanglement entropy. We use low-energy Hamiltonians and numerical result from model systems to show this effect. These transitions are observed in realistic materials, including strained graphene and buckled honeycomb group V (Sb/As).

36 MATERIALS SCIENCE↗
DOE OSTI2021

van der Waals corrected density functionals for cylindrical surfaces: Ammonia and nitrogen dioxide adsorbed on a single-walled carbon nanotube

In this work, we extend the damped Zaremba-Kohn model (dZK) for long-range dispersion interaction between a molecule and a planar surface to molecules adsorbed on a curved cylindrical surface, and employ this extended model as an additive correction to the semilocal density functionals PBE (Perdew-Burke-Ernzerhof) and SCAN (strongly constrained and appropriately normed). The resulting PBE+vdW (van der Waals)-dZK and SCAN+vdW-dZK are applied to two systems, NH3 and NO2 molecules adsorbed on a single-wall carbon nanotube (CNT), for calculations of binding energies and equilibrium distances. For comparison, the results from vdW nonlocal functionals, such as SCAN+rVV10 and PBE+rVV10, are also presented. The binding energies from PBE+rVV10 (Vydrov and Van Voorhis), SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK are about 70–115 meV for the system of CNT + NH 3 and 300–500 meV for the system of CNT + NO 2 . The results from PBE+vdW-dZK and SCAN+vdW-dZK are closer to each other than those from PBE+rVV10 and SCAN+rVV10 are. The relatively closer results from PBE+vdW-dZK and SCAN+vdW-dZK indicate the consistency of our developed vdW–dZK model for cylindrical surfaces. All methods, including PBE, SCAN, PBE+rVV10, SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK, give approximately the same binding energy differences between two adsorption configurations (types I and II) for the two systems. This implies that the two adsorption sites have approximately the same adsorption stability. The exponent of the vdW interaction power law from our vdW-dZK model for the two systems is about 0 at short distance, largely due to the damping factor, and tends slowly to –4 to –4.5 at distances D about 20–50 Å. At even larger distances, the vdW power-law exponent approaches –5. This feature is very similar to the one calculated with random-phase approximation and renormalization group approaches, supporting the applicability of our methods. Our developed vdW-dZK method provides a highly efficient and reliable method for large systems with cylindrical surfaces, such as vdW interactions with nanotubes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Competing Zero-Field Chern Insulators in Superconducting Twisted Bilayer Graphene

The discovery of magic angle twisted bilayer graphene has unveiled a rich variety of superconducting, magnetic, and topologically nontrivial phases. Here, we show that the zero-field states at odd integer filling factors in h-BN nonaligned devices are consistent with symmetry broken Chern insulators, as is evidenced by the observation of the anomalous Hall effect near moiré cell filling factor ν = +1. The corresponding Chern insulator has a Chern number C = ±1 and a relatively high Curie temperature of T c ≈ 4.5 K. In a perpendicular magnetic field above B > 0.5 T we observe a transition of the ν = +1 Chern insulator from Chern number C = ±1 to C = 3, characterized by a quantized Hall plateau with R yx = h/3e 2 . These observations demonstrate that interaction-induced symmetry breaking leads to zero-field ground states that include almost degenerate and closely competing Chern insulators, and that states with larger Chern numbers couple most strongly to the B field. In addition, the device reveals strong superconducting phases with critical temperatures of up to T c ≈ 3.5 K. By providing the first demonstration of a system that allows gate-induced transitions between magnetic and superconducting phases, our observations mark a major milestone in the creation of a new generation of quantum electronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High-efficiency frost and ice control via sensing-assisted nanovibrational slippery surfaces

Frost and ice accretions on surfaces pose persistent challenges across numerous industrial, residential and transportation systems. While various removal strategies exist, they often suffer from limited effectiveness or high energy consumption, such as frosting delay, ice crack generation, and Joule heating. Here, in this work, we report a novel integrated approach combining vibrational quasi-liquid surface (QLS) and capacitive sensing for efficient condensate, frost, and ice management. Compared to Joule heating, our approach does not rely on complete melting and evaporation for removal, resulting in 68% and 95% energy savings for frost and ice removal, respectively. Our QLS coating significantly reduces surface retention forces, achieving 91% and 87% less residual mass compared to hydrophilic surfaces for frost and ice removal through surface nanovibration, respectively. The integrated capacitive sensor provides real-time detection of different phase states, enabling on-demand removal in precise timeframes. This sensor-assisted approach showed 3.8 times lower energy consumption compared to conventional Joule heating for defrosting. This synergistic integration of surface engineering, nanovibration, and intelligent sensing represents a significant advancement in phase change processes, offering an energy-efficient solution for frost and ice mitigation in energy-intensive systems.

Shen, Yuchen [Univ. of Texas at Dallas, Richardson↗
DOE OSTI2022

Effects of doping, hydrostatic pressure, and thermal quenching on the phase transitions and magnetocaloric properties in Mn 1– x Co x NiGe

The effects of doping, hydrostatic pressure, and thermal quenching on the phase transitions and magnetocaloric properties of the Mn 1–x Co x NiGe system have been investigated. Cobalt doping on the Mn site shifted the martensitic structural transition toward lower temperature until it was ultimately absent, leaving only a magnetic transition from a ferromagnetic (FM) to a paramagnetic (PM) state in the high-temperature hexagonal phase. Co-occurrence of the magnetic and structural transitions to form a first-order magnetostructural transition (MST) from the FM orthorhombic to the PM hexagonal phase was observed in samples with 0.05 < x < 0.20. An additional antiferromagnetic–ferromagnetic-like transition was observed in the martensite phase for 0.05 < x < 0.10, which gradually vanished with increasing Co concentration (x > 0.10) or magnetic field (H > 0.5 T). The application of external hydrostatic pressure shifted the structural transition to lower temperature until an MST was formed in samples with x = 0.03 and 0.05, inducing large magnetic entropy changes up to –80.3 J kg –1 K –1 (x = 0.03) for a 7-T field change under 10.6-kbar pressure. Similar to the effects of the application of hydrostatic pressure, an MST was formed near room temperature in the sample with x = 0.03 by annealing at high temperature (1200 °C) followed by quenching, resulting in a large magnetic entropy change of –56.2 J kg –1 K –1 . Furthermore, these experimental results show that the application of pressure and thermal quenching, in addition to compositional variations, are effective methods to create magnetostructural transitions in the MnNiGe system, resulting in large magnetocaloric effects.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗