SEARCH · Engineering Papers
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.
The quantum anomalous Hall effect in graphene, predicted by first-principles calculations, is induced by the proximity effect of a magnetic topological insulator. Our model further reveals a rich phase diagram in graphene/MnBi 2 Te 4 heterostructures.
Multi-twinned structures have been observed in technologically important crystal systems, for example diamond cubic and face-centered cubic (FCC) lattices, that include materials such as diamond, silicon, a wide range of noble metals, and their nanoscale counterparts. Beyond atomic building blocks, the special arrangements also occur in the self-assembly of nanoparticles (NP) and μm-sized colloidal particles and occupy parts of their phase diagrams. Spanning a wide range of length scales, the universality of the structures arises when the systems attempt to achieve multitwinned structures by overcoming geometric misfits during minimizing surface energies with entirely {111} or close-packing facets. While it is fundamental to understand how strain is sustained upon twinned structures and symmetry breaking, the knowledge will be paramount in practical aspects such as guiding and controlling the thin film growth, anisotropic NP growth, and self-assembly of NPs. Au decahedral (Dh) NP, as the most prevalent multi-twinned model system, fits five tetrahedral motifs into a circle by sharing an axis resulting in a geometric misfit angle of 7.35°, or a disclination with power of -7.35°. Postulating how Au FCC lattice adopts the misfit, theoretical models have been developed to address the underlying lattice symmetry and inhomogeneous strain distribution separately. Yet, experimental reports regarding the former have been limited due to the relatively large X-ray beam sizes that do not fit the sizes of NPs. On the other hand, though the latter has been widely adapted in the thermodynamics of small (<10 nm) multi-twinned nanoparticles, previous literature has shown that, at edge length of 17 nm, the theory’s is invalidated by shear strain that is observed in a defect-free Au Dh NP by high-resolution transmission electron microscopy (HR-TEM) imaging. Though the advancement of aberrationcorrected scanning transmission electron microscopy (AC-STEM) imaging and ab initio calculation techniques brings new opportunities, along with challenges in complicated image analysis and limitation in particle size (usually below 10 nm), the gap between nanoscale and mesoscale has never been extended to gain insight from atomic system with straightforward interaction potentials.
Recent theory progress in (3+1)D dynamical descriptions of relativistic nuclear collisions at finite baryon density are reviewed. Heavy-ion collisions at different collision energies produce strongly coupled nuclear matter to probe the phase structure of Quantum Chromodynamics (QCD). Dynamical frameworks serve as a quantitative tool to study properties of hot QCD matter and map collisions to the QCD phase diagram. Outstanding challenges are highlighted when confronting theoretical models with the current and forthcoming experimental measurements from the RHIC beam energy scan program.
Subensemble Acceptance Method (SAM) [1, 2] is an essential link between measured event-by-event fluctuations and their grand canonical theoretical predictions such as lattice QCD. The method allows quantifying the global conservation law effects in fluctuations. In its basic formulation, SAM requires a sufficiently large system such as created in central nucleus-nucleus collisions and sufficient space-momentum correlations. Directly in the spinodal region of the First Order Phase Transition (FOPT) different approximations should be used that account for finite size effects. Thus, we present the generalization of SAM applicable in both the pure phases, metastable and unstable regions of the phase diagram [3]. Obtained analytic formulas indicate the enhancement of fluctuations due to crossing the spinodal region of FOPT and are tested using molecular dynamics simulations. A rather good agreement is observed. Using transport model calculations with interaction potential we show that the spinodal enhancement of fluctuations survives till the later stages of collision via the memory effect [4]. However, at low collision energies the space-momentum correlation is not strong enough for this signal to be transferred to second and third order cumulants measured in momentum subspace. This result agrees well with recent HADES data on proton number fluctuations at $\sqrt{S_{NN}}$ = 2.4 GeV which are found to be consistent with the binomial momentum space acceptance [5].
Resistive heating of a sample in a diamond anvil cell (DAC) can generate a homogeneous temperature field across the sample chamber with reliable temperatures measured by a thermocouple. It is of importance in experiments aiming at exploring phase diagrams and quantifying thermoelastic properties of materials. Here, we present a ring-heater design developed for BX90 diamond anvil cells (DACs). Furthermore, it is made of a ring-shaped aluminum oxide holder hosting a tungsten wire coil inside and coupled with Ar + 2% H 2 gas to prevent oxidation during experiment. This modular plug-and-play design enables in situ studies of samples via x-ray diffraction up to a temperature of 1700 K. Temperature in the BX90 sample volume as measured through a thermocouple was calibrated using the melting point of gold. As an application of this design, we report the thermal expansion coefficient of MgO at 9.5(1) GPa.
In bilayers of ferromagnets and heavy metals, which form the so-called spintronic emitters, the phenomena of ultrafast demagnetization and the inverse spin Hall effect (ISHE) conspire to yield remarkably efficient emission of electric pulses in the THz band. Light-induced demagnetization of the ferromagnet launches a pulse of spin current into the heavy metal, wherein it bifurcates into a radiative charge transient due to the ISHE. The influence of temperature on this combined effect should depend on both the magnetic phase diagram and the microscopic origin of spin Hall conductivity, but its exact dependence remains to be clarified. Here, we experimentally study the temperature dependence of an archetypal spintronic emitter, the Co/Pt bilayer, using electro-optic sampling of the emitted THz pulses in the time domain. The emission amplitude is attenuated with decreasing temperature, consistent with an inverse spin Hall effect in platinum of predominantly intrinsic origin.
We present a systematic analysis of state-of-the-art polarizable and flexible water models from a many-body perspective, with a specific focus on their ability to represent the Born–Oppenheimer potential energy surface of water from the gas to the liquid phase. Using coupled cluster data in the completed basis set limit as a reference, we examine the accuracy of the polarizable models in reproducing individual many-body contributions to interaction energies and harmonic frequencies of water clusters and compare their performance with that of MB-pol, an explicit many-body model that has been shown to correctly predict the properties of water across the entire phase diagram. Based on these comparisons, we use MB-pol as a reference to analyze the ability of the polarizable models to reproduce the energy landscape of liquid water under ambient conditions. We find that, while correctly reproducing the energetics of minimum-energy structures, the polarizable models examined in this study suffer from inadequate representations of many-body effects for distorted configurations. To investigate the role played by geometry-dependent representations of 1-body charge distributions in reproducing coupled cluster data for both interaction and many-body energies, we introduce a simplified version of MB-pol that adopts fixed atomic charges and demonstrate that the new model retains the same accuracy as the original MB-pol model. Based on the analyses presented in this study, we believe that future developments of both polarizable and explicit many-body models should continue in parallel and would benefit from synergistic efforts aimed at integrating the best aspects of the two theoretical/computational frameworks.
Statistical physics models ranging from simple lattice to complex quantum Hamiltonians are one of the mainstays of modern physics that have allowed both decades of scientific discovery and provided a universal framework to understand a broad range of phenomena from alloying to frustrated and phase separated materials to quantum systems. Traditionally, exploration of the phase diagrams corresponding to multidimensional parameter spaces of Hamiltonians was performed using a combination of basic physical principles, analytical approximations, and extensive numerical modeling. However, exploration of complex multidimensional parameter spaces is subject to the classic dimensionality problem, and the behaviors of interest concentrated on low dimensional manifolds remain undiscovered. Here, we demonstrate that a combination of exploration and exploration–exploitation with Gaussian process modeling and Bayesian optimization allows effective exploration of the parameter space for lattice Hamiltonians and effectively maps the regions at which specific macroscopic functionalities or local structures are maximized. We argue that this approach is general and can be further extended well beyond the lattice Hamiltonians to effectively explore the parameter space of more complex off-lattice and dynamic models.
Based on density functional theory (DFT) calculations, we perform an extensive investigation of intrinsic point defects and Xe impurities in uranium monocarbide (UC). The DFT calculations involve both the conventional generalized gradient approximation (GGA) and the GGA+U approach with the Hubbard parametric term (U), using up to 5×5×5 supercells. GGA calculations for the formation energy of intrinsic defects demonstrate the significant effect of using larger supercells than in previous studies. Results confirm that the ⟨111⟩ and ⟨100⟩ dumbbell interstitials are the most stable interstitial configurations for U and C, respectively. The interstitial mechanisms are favored for self-diffusion of both uranium and carbon and diffusion of Xe under equilibrium conditions. Calculations also reveal that the Xe substitutional defect at the C lattice site tends to adopt an off-site configuration, which can be interpreted as a Xe interstitial–C vacancy complex. We also utilize GGA+U to assess the impact of effective U parameter (Ueff) on the results. Moreover, we introduce a method to estimate the carbon chemical potential by fitting the phase diagram composition data and propose a selection Ueff=1.25 eV based on the experimental Xe diffusion activation energy. With this approach, GGA+U calculations reproduce the available experimental data for the formation energy of the carbon Frenkel pair and can explain the stepwise recovery of intrinsic properties and burst Xe release behavior in UC observed in annealing experiments.
Finite temperature auxiliary field-based quantum Monte Carlo methods, including determinant quantum Monte Carlo and Auxiliary Field Quantum Monte Carlo (AFQMC), have historically assumed pivotal roles in the investigation of the finite temperature phase diagrams of a wide variety of multidimensional lattice models and materials. Despite their utility, however, these techniques are typically formulated in the grand canonical ensemble, which makes them difficult to apply to condensates such as superfluids and difficult to benchmark against alternative methods that are formulated in the canonical ensemble. Working in the grand canonical ensemble is furthermore accompanied by the increased overhead associated with having to determine the chemical potentials that produce desired fillings. Given this backdrop, in this work, we present a new recursive approach for performing AFQMC simulations in the canonical ensemble that does not require knowledge of chemical potentials. To derive this approach, we exploit the convenient fact that AFQMC solves the many-body problem by decoupling many-body propagators into integrals over one-body problems to which non-interacting theories can be applied. We benchmark the accuracy of our technique on illustrative Bose and Fermi–Hubbard models and demonstrate that it can converge more quickly to the ground state than grand canonical AFQMC simulations. We believe that our novel use of HS-transformed operators to implement algorithms originally derived for non-interacting systems will motivate the development of a variety of other methods and anticipate that our technique will enable direct performance comparisons against other many-body approaches formulated in the canonical ensemble.
Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.
Temperature-induced metal–insulator transition (MIT) in vanadium dioxide (VO2) has been under intense research interest for decades both theoretically and experimentally. Due to the complex nature of electron correlations, the underlying physics behind the MIT in VO2 has yet to be fully grasped. In this work, we utilize the fine resolution of the scattering-type scanning near-field optical microscope to investigate the MIT in an epitaxial VO2 thin film on the (100)R TiO2 substrate with mid-infrared light. Bidirectional tweed-like metal–insulator phase coexistence patterns are observed and understood under the Landau free energy paradigm. More interestingly, delayed metallic nucleation is observed near the surface cracks due to local strain relief. This research proposes ideas in investigating the temperature–pressure phase diagram and tuning the interplay between local strain and MIT in oxide thin films.
The Jordan-Wigner transformation establishes a duality between $su(2)$ and fermionic algebras. We present qualitative arguments and numerical evidence that when mapping spins to fermions, the transformation makes strong correlation weaker, as demonstrated by the Hartree-Fock approximation to the transformed Hamiltonian. This result can be rationalized in terms of rank reduction of spin shift terms when transformed to fermions. Conversely, the mapping of fermions to qubits makes strong correlation stronger, complicating its solution when one uses qubit-based correlators. The presence of string operators poses challenges to the implementation of quantum chemistry methods on classical computers, but these can be dealt with using established techniques of low computational cost. Here, our proof of principle results for XXZ and J1-J2 Heisenberg (in 1D and 2D) indicate that the JW transformed fermionic Hamiltonian has reduced complexity in key regions of their phase diagrams, and provides a better starting point for addressing challenging spin problems.
Nickel- and cobalt-based superalloys are commonly used as turbine materials for high-temperature applications. However, their maximum operating temperature is limited to about 1100 °C. Therefore, to improve turbine efficiency, current research is focused on designing materials that can withstand higher temperatures. Niobium-based alloys can be considered as promising candidates because of their exceptional properties at elevated temperatures. The conventional approach to alloy design relies on phase diagrams and structure–property data of limited alloys and extrapolates this information into unexplored compositional space. In this work, we harness machine learning and provide an efficient design strategy for finding promising niobium-based alloy compositions with high yield and ultimate tensile strength. Unlike standard composition-based features, we use domain knowledge-based custom features and achieve higher prediction accuracy. We apply Bayesian optimization to screen out novel Nb-based quaternary and quinary alloy compositions and find these compositions have superior predicted strength over a range of temperatures. We develop a detailed design flow and include Python programming code, which could be helpful for accelerating alloy design in a limited alloy data regime.
Sb thin films have attracted wide interest due to their tunable band structure, topological phases, high electron mobility, and thermoelectric properties. We successfully grow epitaxial Sb thin films on a closely lattice-matched GaSb(001) surface by molecular beam epitaxy. We find a novel anisotropic directional dependence on their structural, morphological, and electronic properties. The origin of the anisotropic features is elucidated using first-principles density functional theory (DFT) calculations. The growth regime of crystalline and amorphous Sb thin films was determined by mapping the surface reconstruction phase diagram of the GaSb(001) surface under Sb 2 flux, with confirmation of structural characterizations. Crystalline Sb thin films show a rhombohedral crystal structure along the rhombohedral (211) surface orientation parallel to the cubic (001) surface orientation of the GaSb substrate. At this coherent interface, Sb atoms are aligned with the GaSb lattice along the [1̄10] crystallographic direction but are not aligned well along the [110] crystallographic direction, which results in anisotropic features in reflection of high-energy electron diffraction patterns, misfit dislocation formation, surface morphology, and transport properties. Our DFT calculations show that the preferential orientation of the rhombohedral Sb (211) plane may originate from the GaSb surface, where Sb atoms align with the Ga and Sb atoms on the reconstructed surface. The formation energy calculations confirm the stability of the experimentally observed structures. Our results provide optimal film growth conditions for further studies of novel properties of Bi 1-x Sb x thin films with similar lattice parameters and an identical crystal structure, as well as functional heterostructures of them with III–V semiconductor layers along the (001) surface orientation, supported by a theoretical understanding of the anisotropic film orientation.
The (PbTiO3)n/(SrTiO3)n (PTO/STO) superlattice system has been shown to exhibit interesting topological phases (e.g., vortices and skyrmions) in addition to normal ferroelectric domain states. Existing studies are mostly focused on the dependence of topological polar distributions and properties of PTO/STO superlattice on its periodicity. Here, we study the strain effect on the topological phase transitions and ferroelectric domain structures employing phase-field simulations. We summarized in an isotropic strain (in-plane misfit strain along the x direction is equal to that along the y direction) periodicity phase diagram displaying the stability regions of different polar topological states, including normal ferroelectric twins, vortices, skyrmions, and mixtures of vortices and twins. We also analyzed the polarization configurations under anisotropic in-plane strains (in-plane misfit strain along the x direction is not equal to that along the y direction) and demonstrated that the strain anisotropy can be used to tune the directions of vortex arrays along either the [100] pc or [010] pc directions or labyrinth vortex arrays. This work offers guidance to manipulating polar structures in the PTO/STO superlattices via strain engineering.
This contribution provides a thorough examination of the structural characteristics of pyrochlore-type lanthanide titanates and zirconates Ln2Ti2O7 and Ln2Zr2O7, across various length scales. This paper also examines their processing, interesting physical properties (electrical, magnetic, and thermal characteristics), and responses to high pressure and ion irradiation. Brief sections on the elemental oxides' crystal chemistry, pertinent phase diagrams, and energetics of defect formation are also provided. Pyrochlore-type Ln2Ti2O7 and Ln2Zr2O7 stand out as truly multifunctional materials. Moreover, they have emerged as fascinating materials due to magnetic geometrical frustration, arising from the ordering of magnetic Ln3+ and non-magnetic Ti4+ (or Zr4+) cations into separate, interpenetrating lattices of corner-sharing tetrahedra. This results in a diverse array of exotic magnetic ground states, such as spin-ice (e.g., Dy2Ti2O7 or Ho2Ti2O7) or quantum spin ice (e.g., Tb2Ti2O7), observed at both low and room temperatures. They also exhibit varied electrical and electrochemical characteristics. Some members such as Gd2Zr2O7, function as fast ion conductors with a conductivity (σ) of ≈10−2 S·cm−1 at 800 °C and activation energy (Ea) ranging from 0.85 to 1.52 eV, depending on the degree of structural disorder. Others, such as Gd2TiMoO7, are mixed ionic-electronic conductors with σ ≈ 25 S·cm−1 at 1000 °C, making them promising candidate materials for applications in energy conversion and storage devices and oxygen separation membranes. Their exceptionally low thermal conductivity (e.g., κ ∼ 1.1–1.7 W·m−1·K−1 between 700 and 1200 °C for Ln2Zr2O7), close to the glass-like lower limit of highly disordered solids, positions them as valuable materials for thermal barrier coatings. They can also effectively accommodate actinides (e.g., Pu, Np, Cm, Am) in solid solutions and sustain prolonged exposure to radiation due to alpha-decay events, while preserving the integrity of the periodic atomic structure. Proposed as major components in actinide-bearing ceramics, they contribute to the long-term immobilization and disposal of long-lived waste radionuclides from nuclear programs. Some of these properties are displayed simultaneously, opening avenues for new applications. Despite the wealth of data available in the literature, this review highlights the need for a better understanding of order/disorder processes in pyrochlore-type materials and the influence of the structural length scale on their physical and chemical properties. Recent experimental evidence has revealed that pyrochlore short-range structure is far more complex than originally thought. Moreover, pyrochlore local structure is now believed to include short-range, lower symmetry, ordered domains, such as the orthorhombic weberite-type of structure. Notably, short- and long-range structures appear decoupled across different length scales and temperature regimes, and these differences persist even in well-ordered samples. We believe that the pyrochlore structure offers a unique opportunity for examining the interplay between chemical composition, defect chemistry, and properties. In Memoriam: Rodney C. Ewing, Fondly Remembered.
Lead-free piezoelectric ceramics of (1 − x)K0.5Na0.5NbO3-xCaTiO3 were fabricated, and their crystal structure, microstructure, and electrical properties were systematically studied. Rietveld refinement of the x-ray diffraction data and Raman spectroscopic analyses revealed a composition-dependent structural phase transition: three phase transitions, namely, from a pure orthorhombic phase for x ≤ 0.02 to a mixed phase of orthorhombic and tetragonal phases (0.03 ≤ x ≤ 0.08) and finally another mixed phase of tetragonal + cubic for x = 0.10 and 0.15 at room temperature (RT). The morphological study reveals a decrease in grain size along with a more uniform distribution of grains as the concentration of CaTiO3 (CT) increases; notably, a homogeneous distribution of grains is observed for x = 0.05. The temperature-dependent dielectric properties show two phase transitions, from orthorhombic to tetragonal (TO-T) and tetragonal to cubic (TC), for unmodified K0.5Na0.5NbO3 (KNN). However, both the phase transition temperatures (TO-T and TC) decrease, and the transition peaks broaden with an increase in CT substitution, and for x > 0.06, the TO-T shifted below RT. The broadening of the transition peak at TO-T may be due to the relaxation behavior. Among the prepared samples, the 5 mol. % CT-modified KNN shows the optimum electrical properties (d33 = 114 pC/N, ɛr = 412, and 2Pr = 15.25 μC/cm2) at RT. The enhanced electrical properties for x = 0.05 are due to the coexistence of orthorhombic and tetragonal phases, facilitating easy polarization rotation and flattening of the free energy profile. A phase diagram has been constructed based on the information gathered from the temperature-dependent dielectric measurements, RT x-ray diffraction, and Raman spectroscopy data and is discussed in detail.