SEARCH · Engineering Papers

Results for “second order phase transitions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1990

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 140 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce thus that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism/ - Pauli matrices/. The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗

NASA NTRS1991

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 240 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce, thus, that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism (- Pauli matrices). The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗

DOE OSTI2024

Noncollinear polar magnet Fe 2 (SeO 3 ) 3 (H 2 O) 3 with inequivalent Fe 3+ sites

The emergence of novel magnetic states becomes more likely when the inversion symmetry of the crystal field, relative to the center between two spins, is broken. We propose that placing magnetic spins in inequivalent sites in a polar lattice can promote a realization of nontrivial magnetic states and associated magnetic properties. To test our hypothesis, we study Fe 2 (SeO 3 )(H 2 O) 3 as a model system that displays two distinct Fe(1) and Fe(2) magnetic sites in a polar structure (R3c space group). At low fields μ 0 H≤ 0.06 T, the material undergoes an antiferromagnetic ordering with T N1 = 77 K and a second transition at T N2 ≈ 4 K. At μ 0 H≥ 0.06 T and 74 K ≤T≤ 76 K, a positive entropy change of ∼0.12 mJ mol −1 K −1 can be associated with a metamagnetic transition to possibly nontrivial spin states. At zero field, Fe(1) is nearly fully ordered at T≈ 25 K, while Fe(2) features magnetic frustration down to T = 4 K. The magnetic ground state, a result corroborated by single-crystal neutron diffraction and 57 Fe Mössbauer spectroscopy, is a noncollinear antiparallel arrangement of ferrimagnetic Fe(1)–Fe(2) dimers along the c-axis. The results demonstrate that placing distinct magnetic sites in a polar crystal lattice can enable a new pathway to modifying spin, orbital, and lattice degrees of freedom for unconventional magnetism.

Crystal lattices↗

NASA NTRS2001

NASA Alternate Access to Station Service Concept

The evolving nature of the NASA space enterprise compels the agency to develop new and innovative space systems concepts. NASA, working with increasingly strained budgets and a declining manpower base, is attempting to transform from operational activities to procurement of commercial services. NASA's current generation reusable launch vehicle, the Shuttle, is in transition from a government owned and operated entity to a commercial venture to reduce the civil servant necessities for that program. NASA foresees its second generation launch vehicles being designed and operated by industry for commercial and government services. The "service" concept is a pioneering effort by NASA. The purpose the "service" is not only to reduce the civil servant overhead but will free up government resources for further research - and enable industry to develop a space business case so that industry can sustain itself beyond government programs. In addition, NASA desires a decreased responsibility thereby decreasing liability. The Second Generation Reusable Launch Vehicle (RLV) program is implementing NASA's Space Launch Initiative (SLI) to enable industry to develop the launch vehicles of the future. The Alternate Access to Station (AAS) project office within this program is chartered with enabling industry to demonstrate an alternate access capability for the International Space Station (ISS). The project will not accomplish this by traditional government procurement methods, not by integrating the space system within the project office, or by providing the only source of business for the new capability. The project funds will ultimately be used to purchase a service to take re-supply cargo to the ISS, much the same as any business might purchase a service from FedEx to deliver a package to its customer. In the near term, the project will fund risk mitigation efforts for enabling technologies. AAS is in some ways a precursor to the 2nd Generation RLV. By accomplishing ISS resupply with existing technologies, not only will a new category of autonomous vehicles deliver cargo, but a commercial business base will be incubated that will improve the likelihood of commercial convergence with the next generation of RLVs. Traditional paradigms in government management and acquisition philosophy are being challenged in order to bring about the objective of the AAS project. The phased procurement approach is proving to be the most questionable aspect to date. This work addresses the fresh approach AAS is adopting in management and procurement through a study of the AAS history, current solutions, key technologies, procurement complications, and an incremental forward plan leading to the purchase of a service to deliver goods to ISS. Included in this work is a discussion of the Commercial Space Act of 1998 and how it affects government purchase of space launch and space vehicle services. Industry should find these topics pertinent to their current state of business.

Bailey, Michelle D.↗

NASA NTRS2026

Spacecraft Radiator Protection from Ionizing Radiation, Dust, and Excessive Heat Loss

Under a Phase II SBIR project funded by NASA Johnson Space Center (Contract No. 80NSSC25C0088), Analytical Scientific Products LLC (ASP) has been developing an actively controlled louver to protect spacecraft radiators from degradation due to exposure to various types of environmental conditions. Of particular interest are ionizing radiation during spacecraft transit through the Van Allen belts, dust during spacecraft landing and surface operations on the moon, and excessive heat loss during the long lunar night especially near the poles where the local ambient temperatures can drop below -200°C. Exposure to ionizing radiation and dust can degrade the optical properties of the radiator coating that in turn reduces its ability to reject excess heat from the spacecraft into the environment. Exposure to the extremely low temperature conditions during the lunar night near the poles can freeze the radiator fluids that can compromise the integrity of the radiator. Passive louvers constructed using thick aluminum vanes are used currently to protect spacecraft radiators from some of the above effects, but they have several drawbacks: (i) their high aerial density makes it impractical to scale them to protect the much larger spacecraft radiators needed for future manned missions to the moon and beyond, (ii) the bimetallic actuators used to open and close the vanes in passive louvers rely on external temperature alone and so the louver cannot offer protection from dust and ionizing radiation when the ambient temperatures are high, and (iii) the bimetallic actuators need time scales of the order of hours to open and close. We have addressed the above limitations of passive louvers through the development of a low aerial density and rapidly actuating actively controlled louver. It is constructed from a custom alloy that offers much higher levels of protection against the ionizing radiation, dust and excessive heat loss at a fraction of weight compared to passive louvers while providing opening and closing time scales of the order of a second. Our modular design makes it easy to scale the system up or down depending on the application. We have recently constructed a 31 inch × 31 inch module of this louver and tested its functionality and effectiveness. These tests have shown that the louver can be opened and closed in less than 2 s. It reduces the transmission of ionizing radiation by 78% to 100% (depending on the radiation source), dust transmission by more than 93%, and heat loss by more than 97%. We are currently preparing to evaluate this device at the Johnson Space Center’s cryogenic vacuum chamber where it can be subjected to simulated lunar surface conditions. This paper provides a detailed discussion of the test designs as well as the data from tests that demonstrate the effectiveness in reducing the transmission of ionizing radiation, dust, and heat under laboratory conditions.

Radiator↗

NASA NTRSDate not supplied

Spacecraft Radiator Protection from Ionizing Radiation, Dust, and Excessive Heat Loss

Radiator↗

NASA NTRS2024

Aerocapture: An Enabling Technology for Flagship-Class Uranus Orbiter and Probe Mission

Introduction: The current decadal survey published by the National Academies of Sciences has informed National Aeronautics and Space Administration (NASA) to prioritize the study of the Ice Giants, especially Uranus. To gather the required data that addresses the science questions raised in this survey, a mission to Uranus with an orbiter and atmospheric probe must be designed. The Uranus Orbiter and Probe (UOP) study, which the survey identified as the flagship mission of this decade, proposes a 2031 launch to take advantage of a Jupiter fly-by and utilizes a fully propulsive orbit insertion design with an Earth-to-Uranus transit times ranging from 13 to 15 years. This fully propulsive orbit insertion at Uranus will be very fuel expensive (wet mass percentages of around 60-70\%) thereby leaving less mass for the scientific payload and additional planetary probes. In addition, scientists are more interested in visiting Uranus before 2049, when the Spring Equinox will occur, as it allows studying Uranus seasons not seen during Voyager 2's flyby in 1987. A NASA Flagship-Class mission would require at least 10 years of lead time prior to launch thereby making the 2031 launch to take advantage of the Jupiter fly-by extremely challenging. The consequence of missing the Jupiter fly-by and launching in the late 2030s is the challenge of a fully-propulsive mission like UOP to have a feasible alternative interplanetary trajectory that reaches Uranus before 2050. As an alternative, to address the shortcomings of the fully propulsive mission, one can design a mission to Uranus using aerocapture. What is aerocapture: Aerocapture is an atmospheric maneuver that uses aerodynamic forces generated by flight within a planetary atmosphere to decelerate and achieve orbit insertion. Using aerocapture allows one to provide the change in velocity (Delta V) needed to slow down from the approach hyperbolic trajectory to achieve the desired captured orbit around the target planet using aerodynamic forces of the vehicle (lift and drag). Using aerodynamic forces instead of fully-propulsive maneuvers results in significant savings in the fuel. Furthermore, aerocapture can also allow one to consider interplanetary trajectories with higher approach hyperbolic velocities, thus reducing the mission transit times. Aerocapture as an enabling technology: To use aerocapture as an enabling technology for the Uranus exploration mission, one would require an integrated system-level design, including a Thermal Protection System (TPS), hardware needed for aerodynamic modulation, and autonomous Guidance, Navigation, and Control (GNC) systems. Aerocapture has yet to be demonstrated, despite considering it for several past missions. Recent advancements in TPS and GNC capabilities show the potential for using rigid, heritage entry vehicle configurations already flown at other planetary bodies for Uranus aerocapture. Aerocapture can be a robust technology that can deliver spacecraft to Uranus science orbits while substantially increasing on-orbit payload mass (more than 40\%) that can enable robust atmospheric entry probes. Additionally, the aerocapture maneuver would reduce the interplanetary transit time by 2–5 years (15-30\%) relative to fully-propulsive orbit insertion. Recent work has shown that one can conduct a flagship-class mission in a shorter period than fully-propulsive missions if using aerocapture. What does aerocapture bring in for a Uranus Mission: Using aerocapture for a Uranus orbiter and atmospheric probe mission can provide one with considerable propellant savings. Spacecraft in interplanetary trajectories to Uranus typically need an Delta V in orders of kilometers per second to insert into science orbit. One would require thousands of kilograms of fuel to achieve such a Delta V using a traditional fully propulsive maneuver, leaving less mass for payload during the mission launch. Aerocapture can reduce the propulsion needs by dissipating energy in the sizable atmosphere of Uranus without a significant mass increase due to the need for an aeroshell. One can use the mass savings achieved using the aerocapture to reduce the launch vehicle requirements. In addition, one can also have additional science instruments on the orbiter or create a robust instrumentation suite on an atmospheric probe that can significantly increase the science outcome of the Uranus exploration mission. Furthermore, since aerocapture performance is relatively insensitive to increases in hyperbolic excess velocity, one can design the interplanetary trajectory to arrive at Uranus faster, reducing the interplanetary transit time and operations cost. All these savings, achieved using aerocapture, could help fit a larger class mission, such as the Uranus mission within a smaller capital, e.g., a Flagship-class orbiter mission in a New Frontiers class capital. Summary: This talk will provide an overview of how aerocapture can enable the Uranus exploration mission. Specifically, this talk will discuss the latest advancements made in the Uranus aerocapture study, such as investigating interplanetary trajectories with higher hyperbolic approach velocities and their implications on the aero heating and the TPS design, incorporating FNPAG (an advanced numerical-predictor guidance) and comparison of multiple navigation approaches. In addition, this talk will focus on mechanical design that can house more than one atmospheric probe and the six degrees of freedom simulation of aerocapture at Uranus. Findings from a recent NASA Space Technology Mission Directorate (STMD)-funded activity studying the aerocapture as an enabling technology for a Uranus orbiter will be presented. Using the science payload recommended by the Decadal Survey for Uranus exploration, this work shows many improvements over the baseline fully-propulsive mission. These improvements include a shorter cruise phase, flexibility in launch opportunities late into the 2030s while reaching Uranus before the 2050 equinox for the desired science opportunities, and lower propellant mass needs. This talk will highlight how aerocapture can be utilized for Uranus science orbit insertion using a lower-risk, heritage entry vehicle configuration used extensively as a Mars entry, descent, and landing vehicle. Furthermore, this talk will explore how the demonstration of aerocapture at Earth can benefit the aerocapture-enabled Uranus mission.

Aerocapture↗

DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2026

Altermagnetic polar metallic phase in ultrathin epitaxially strained RuO 2 films

Altermagnetism refers to a wide class of magnetic orders featuring magnetic sublattices with opposite spins related by rotational symmetries, resulting in nontrivial spin splitting and magnetic multipoles. However, the direct observation of the altermagnetic transition remains elusive. Here, by combining theoretical analysis, electrical transport, X-ray, and optical spectroscopies, we establish a phase diagram in hybrid molecular beam epitaxy-grown RuO 2 /TiO 2 (110) films, mapping symmetries along with altermagnetic/electronic/structural phase transitions as functions of film thickness and temperature. This features an altermagnetic metallic polar phase in epitaxially strained 2 nm films, suggesting a potential link between polar metals and altermagnetic materials. Such a clear signature of a magnetic phase transition at ~500 K is observed exclusively in ultrathin strained films, unlike in bulk RuO 2 single crystals. These results highlight the power of epitaxial heterostructure engineering to induce altermagnetism in systems initially nonmagnetic, opening avenues for realizing emergent quantum phases with multifunctional properties.

RuO2 thin film↗

DOE OSTI2021

Optimal surface-tension isotropy in the Rothman-Keller color-gradient lattice Boltzmann method for multiphase flow

The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Here, two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents. Here, we investigate the accuracy of the CG calculation for interfaces between fluids with several radii of curvature and find that the standard CG calculations lead to significant inaccuracy. Specifically, we observe significant anisotropy of the color gradient of order 7% for high curvature of an interface such as when a pinchout occurs. We derive a second order accurate color gradient and find that the diagonal nearest neighbors can be weighted differently than in the usual color-gradient calculation such that anisotropy is minimized to a fraction of a percent. The optimal weights that minimize anisotropy for the smallest radius of curvature interface are found to be w = (0.298, 0.284, 0.275) for diagonal nearest neighbors for the cases of the interface smoothing parameter β = (0.5, 0.7, 0.99), somewhat higher than the w = 0.25 value derived by Leclaire et al. [Leclaire, Reggio, and Trepanier, Computers and Fluids 48, 98 (2011)] based on obtaining isotropic errors to second order. We find that use of these optimal w values yields over a factor of 10 decrease in anisotropy and over a factor of 30 decrease in mean anisotropy relative to using the standard w = 1 value. And we find a factor of about 2 decrease in the anisotropic error and up to factor 15 decrease in mean anisotropic error relative to the choice of w = 0.25 for small radius of curvature interfaces. The improved CG calculations will allow the method to be more reliably applied to studies of phenomenology and pore scale processes such as viscous and capillary fingering, and droplet formation where surface-tension isotropy of narrow fingers and small droplets plays a crucial role in correctly capturing phenomenology. We present an example illustrating how different phenomena can be captured using the improved color-gradient method. Namely, we present simulations of a wetting fluid invading a fluid filled pipe where the viscosity ratio of fluids is unity in which droplets form at the transition to fingering using the improved CG calculations that are not captured using the standard CG calculations. We present an explanation of why this is so which relates to anisotropy of the surface tension, which inhibits the pinchouts needed to form droplets.

58 GEOSCIENCES↗

DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide[1]Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo2C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo2C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo2C to H2S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS2. The formation of MoS2 led to a deficiency of Mo atoms in the underlying Mo2C, resulting in an inhomogeneous phase change from α-Mo2C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo2C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo2C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional materials↗

DOE OSTI2020

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Broken symmetries associated with a Kagome chiral charge order

Chirality, or handedness, is ubiquitous in science, from cell biology to physics, and in condensed matter can underlie exotic phases such as chiral charge density waves and chiral superconductivity. However, detecting subtle broken symmetries that define such states is challenging, leading to debate and controversy. Here, using second-order optical response, we reveal the broken symmetries of a chiral charge density wave in the Kagome lattice KV 3 Sb 5 . Polarization-dependent mid-infrared photocurrent microscopy uncovers a longitudinal, helicity-dependent photocurrent associated with the charge order, indicating broken inversion and mirror symmetries. These findings, supported by theoretical analysis, directly establish the intrinsic chiral nature of the ordered state. Moreover, the absence of a circular photogalvanic effect perpendicular to the incident light imposes stringent constraints on the point-group symmetries. Our study not only visualizes the chiral nature of the Kagome charge order, but also highlights the nonlinear photogalvanic effect as a sensitive probe for detecting subtle symmetry breakings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo 2 C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo 2 C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo 2 C to H 2 S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS 2 . The formation of MoS 2 led to a deficiency of Mo atoms in the underlying Mo 2 C, resulting in an inhomogeneous phase change from α-Mo 2 C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo 2 C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo 2 C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional Materials↗

NASA NTRS2018

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗

DOE OSTI2026

Phase diagram and spectroscopic signatures of a supersolid in the quantum ising magnet K 2 Co(SeO 3 ) 2

Supersolid phases are quantum-entangled states of matter exhibiting the dual characteristics of superfluidity and solidity. Theory predicts that hard-core bosons on a triangular lattice can form such phases at half filling and near complete filling. Leveraging an exact mapping between bosons and spin-$\frac{1}{2}$ degrees of freedom, here we show that these phases are realized in the triangular-lattice antiferromagnet K 2 Co(SeO 3 ) 2 . At zero field, neutron diffraction reveals the development of quasi-two-dimensional $\sqrt3$ x $\sqrt3$ magnetic order with Z 3 translational symmetry breaking (solidity), though with reduced amplitude indicating strong quantum fluctuations. These fluctuations manifest as equidistant bands of continuum neutron scattering, where the lowest-energy mode is gapless at K ($\frac{1}{3}$ $\frac{1}{3}$), consistent with broken U(1) spin rotational symmetry (superfluidity). For c-axis-oriented magnetic fields near saturation, we find a second phase consistent with a high-field supersolid. These two supersolids are separated by a pronounced 1/3 magnetization plateau phase that supports coherent spin waves, from which we determine the underlying spin Hamiltonian.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Strangeness-correlations on the pseudocritical line in ( 2 + 1 )-flavor QCD

We present some lattice QCD results on first ( χ 1 i ) and second ( χ 2 i ) cumulants of and correlations ( χ 11 i j ) among net baryon-number ( B ), strangeness ( S ) and electric charge ( Q ) along the pseudocritical line [ T p c ( μ B ) ] in the temperature ( T )–baryon chemical potential ( μ B ) phase diagram of ( 2 + 1 )-flavor QCD. We point out that violations of sum rules among second order cumulants, which hold in the isospin symmetric limit of vanishing electric charge chemical potential, are small along the T p c ( μ B ) for the entire range of μ B covered in the RHIC beam energy scan. For the strangeness neutral matter produced in heavy-ion collisions this leads to a close relation between χ 11 B S and χ 11 Q S . We compare lattice QCD results for χ 11 B S / χ 2 S along the T p c ( μ B ) line with preliminary experimental measurements of χ 11 B S / χ 2 S for collision energies 7.7 GeV ≤ s N N ≤ 62.4 GeV . While we find good agreements for s N N ≥ 39 GeV , differences are sizeable at smaller values of s N N . Moreover, we compare lattice QCD results for the ratio of the strangeness ( μ S ) to baryon ( μ B ) chemical potentials, which define a strangeness neutral system with fixed electric charge to baryon number density, with experimental results obtained by the STAR collaboration for μ S / μ B using strange baryon yields on the freeze-out line. Finally, we determine the baryon chemical potential at the freeze-out ( μ B f ) by comparing χ 1 B / χ 2 B along the T p c ( μ B ) with the experimentally measured net-proton cumulants χ 1 p / χ 2 p . We find that { μ B f , T p c ( μ B f ) } are consistent with the freeze-out parameters of the statistical-model fits to experimentally measured hadron yields for s N N ≥ 11.5 GeV . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2002

Dynamics of Sheared Granular Materials

This work focuses on the properties of sheared granular materials near the jamming transition. The project currently involves two aspects. The first of these is an experiment that is a prototype for a planned ISS (International Space Station) flight. The second is discrete element simulations (DES) that can give insight into the behavior one might expect in a reduced-g environment. The experimental arrangement consists of an annular channel that contains the granular material. One surface, say the upper surface, rotates so as to shear the material contained in the annulus. The lower surface controls the mean density/mean stress on the sample through an actuator or other control system. A novel feature under development is the ability to 'thermalize' the layer, i.e. create a larger amount of random motion in the material, by using the actuating system to provide vibrations as well control the mean volume of the annulus. The stress states of the system are determined by transducers on the non-rotating wall. These measure both shear and normal components of the stress on different size scales. Here, the idea is to characterize the system as the density varies through values spanning dense almost solid to relatively mobile granular states. This transition regime encompasses the regime usually thought of as the glass transition, and/or the jamming transition. Motivation for this experiment springs from ideas of a granular glass transition, a related jamming transition, and from recent experiments. In particular, we note recent experiments carried out by our group to characterize this type of transition and also to demonstrate/ characterize fluctuations in slowly sheared systems. These experiments give key insights into what one might expect in near-zero g. In particular, they show that the compressibility of granular systems diverges at a transition or critical point. It is this divergence, coupled to gravity, that makes it extremely difficult if not impossible to characterize the transition region in an earth-bound experiment. In the DE modeling, we analyze dynamics of a sheared granular system in Couette geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better understand what we might encounter in a reduced-g environment, and at a deeper level to deduce the physics of sheared systems in a density regime that has not been addressed by past experiments or simulations. One aspect of the simulations addresses sheared 2D system in zero-g environment. For low volume fractions, the expected dynamics of this type of system is relatively well understood. However, as the volume fraction is increased, the system undergoes a phase transition, as explained above. The DES concentrate on the evolution of the system as the solid volume fraction is slowly increased, and in particular on the behavior of very dense systems. For these configurations, the simulations show that polydispersity of the sheared particles is a crucial factor that determines the system response. Figures 1 and 2 below, that present the total force on each grain, show that even relatively small (10 %) nonuniformity of the size of the grains (expected in typical experiments) may lead to significant modifications of the system properties, such as velocity profiles, temperature, force propagation, and formation shear bands. The simulations are extended in a few other directions, in order to provide additional insight to the experimental system analyzed above. In one direction, both gravity, and driving due to vibrations are included. These simulations allow for predictions on the driving regime that is required in the experiments in order to analyze the jamming transition. Furthermore, direct comparison of experiments and DES will allow for verification of the modeling assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of these simulations of an annular system in zero-g environment will conclude the presentation.

Kondic, Lou↗