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.
Paramagnetic single-molecule magnets (SMMs) interacting with the ferromagnetic electrodes of a magnetic tunnel junction (MTJ) produce new molecular spintronics testbed and highly ordered magnetic metamaterial promising for room temperature.
We used angle-resolved photoemission spectroscopy (ARPES) and density functional theory calculations to study the electronic properties of MnBi 2 Te 4 , a material that was predicted to be an intrinsic antiferromagnetic (AFM) topological insulator. In striking contrast to earlier literature showing a full gap opening between two surface band manifolds on the (0001) surface, we observed a gapless Dirac surface state with a Dirac point sitting at E B = -280 meV . Furthermore, our ARPES data revealed the existence of a second Dirac cone sitting closer to the Fermi level. Surprisingly, these surface states remain intact across the AFM transition. The presence of gapless Dirac states in this material may be caused by different ordering at the surface from the bulk or weaker magnetic coupling between the bulk and surface. Whereas the surface Dirac cones seem to be remarkably insensitive to the AFM ordering most likely due to weak coupling to magnetism, we did observe a splitting of the bulk band accompanying the AFM transition. Finally, with a moderately high ordering temperature and interesting gapless Dirac surface states, MnBi 2 Te 4 provides a unique platform for studying the interplay between magnetic ordering and topology.
The Fe-based superconductor Fe (Se,Te) combines non-trivial topology with unconventional superconductivity and may be an ideal platform to realize exotic states such as high-order topological corner modes and Majorana modes. Thin films of Fe (Se,Te) are particularly important for device fabrication and phase sensitive transport measurements. While bulk Fe (Se,Te) has been extensively studied, the nature of the superconducting order parameter in the monolayer limit has not yet been explored. In this work, we study monolayer films of Fe (Se,Te) on Bi 2 Te 3 with scanning tunneling spectroscopy. Monolayer Fe (Se,Te)/Bi 2 Te 3 heterostructures host a multigap superconducting state that strongly resembles the bulk. Analysis of the phase-referenced quasiparticle interference signal reveals a sign-changing s-wave order parameter similar to the bulk as well as a unique pattern of sign changes which have not been observed in the bulk. Our work establishes monolayer Fe (Se,Te)/Bi 2 Te 3 as a robust multi-band unconventional superconductor and sets the stage for explorations of non-trivial topology in this highly-tunable system.
Utilizing the self-assembly of block copolymers with large Flory–Huggins interaction parameters (χ) for nanofabrication is a formidable challenge due to the attendant large surface energy differences between the blocks. This work reports a robust protocol for the fabrication of thin films with highly ordered cylindrical nanopore arrays via the self-assembly of an asymmetric poly(styrene-block-4-vinylpyridine) (PS-b-P4VP) diblock copolymer blended with a P4VP homopolymer. The desired vertical domain orientation is achieved at the air–polymer interface by controlled solvent vapor annealing (SVA) using acetone, a solvent with weak selectivity for PS over P4VP, and at the substrate interface by functionalization using a hydroxy-terminated poly(2-vinylpyridine) (P2VP-OH) homopolymer brush. In contrast, the vertical cylinder orientation is unstable during acetone SVA on substrates functionalized using hydroxy-terminated poly(methyl methacrylate) (PMMA-OH). Although PMMA exhibits more balanced interfacial energies between PS and P4VP than P2VP in the dry state, it is also swollen more selectively by acetone. We hypothesize that the nearly balanced solvent swelling of the three polymers (P2VP, P4VP, and PS) stabilizes the vertical cylinder orientation, while unbalanced swelling (PMMA > P4VP and PS) does not. Here we further characterize pore formation by addition of a P4VP homopolymer and its postassembly extraction using ethanol, revealing a narrow window of pore size tunability. Notably, minimal differences in nanopore morphologies are observed for P4VP volume fractions as high as 0.1, regardless of the P4VP molar mass. However, further increasing the P4VP volume fraction results in domain reorientation or macrophase separation when its molar mass is less than or greater than the P4VP block molar mass, respectively. Using a P4VP homopolymer that is nearly equal in length to the P4VP block enables the fabrication of well-ordered arrays of vertical, through-film nanopores with high aspect ratios (>10), small periods (<23 nm), and diameters less than 10 nm.
The practical application of Li-rich Mn-based layered oxides (LLO) cathode is hindered by severe capacity and voltage degradation resulting from severe oxygen release and irreversible phase transition. In this work, the reaction heterogeneity, which describes the spatially resolved electrochemical divergence within individual cathode particles, is engineered through compositional gradient design to couple Li + transport kinetics and the anion redox activity between particle interiors and surfaces. It is revealed that Co/Mn concentration gradient within particles creates heterogeneous phase content distribution and structural ordering, inducing surface-bulk reaction heterogeneity that significantly impacts the overall electrochemical performance. Specifically, Li 2 MnO 3 -poor and Co-enriched surface effectively mitigates the oxygen loss and enhances electrochemical reaction kinetics, benefited from the reduced surface redox reactivity and induced highly ordered intra-layered cationic arrangement. Meanwhile, the Li 2 MnO 3 -enriched core with slight Li/Ni intermixing provides high reversible capacity and strong mechanical stability. Consequently, the greatly enhanced anion redox reversibility, Li + diffusion dynamics, and structure stability endow LLO with exceptional electrochemical properties, showing a capacity retention of 86.0% and a reduced voltage decay of 0.518 mV per cycle after 500 cycles at 1 C. This work provides a valuable strategy to tailor the redox chemistry and achieve robust LLO.
Covalent organic frameworks (COFs) are an emerging type of porous crystalline material with highly ordered aperture size and tunable structures with designer properties. COFs have been proposed as promising materials for water treatment because of their notable intrinsic properties like excellent chemostability, high surface area, abundant functional sites, and uniform adjustable aperture size. This review focuses on fundamental COF design principles for water treatment (stability, aperture size, and surface functionalization) and the state-of-the-art application of COFs in desalination, organic contaminant sorption, and ion capture. Additional potential promising applications of COFs for water treatment, including solar steam generation, photocatalysis for degradation of organic contaminants, and capacitive deionization are also presented along with an outlook toward future opportunities in the field.
The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab aims to measure the muon anomalous magnetic moment $a_{\mu}$ with unprecedented precision, potentially uncovering physics beyond the Standard Model of particle physics. The result based on Runs 1-3, released in 2023, achieved a precision of 0.20 ppm. The experiment circulates muons in a storage ring, measuring $a_{\mu}$ from decay positron time and energy measurements collected with calorimeters. To achieve the required accuracy, it is crucial to measure and control the magnetic field in the ring with high precision. Beam dynamics corrections are necessary for muons not orbiting exactly in the midplane, for their oscillations, and for electric field effects. Highly accurate beam dynamics simulations are instrumental for quantifying and validating the beam dynamics corrections, ultimately improving the precision of the $a_{\mu}$ measurement and facilitating the achievement of the ambitious $70\:\mathrm{ppb}$ systematic uncertainty goal. The measured field data was incorporated into models for simulations using three codes: \texttt{gm2ringsim} (an internal Geant4-based code), \textit{COSY INFINITY}, and \textit{BMAD}. The advantages of \texttt{gm2ringsim} include using CAD-based geometry and modelling the detector effects. \textit{COSY INFINITY} is a highly accurate and efficient code that uses high-order differential-algebraic transfer maps, precise fringe field calculations, and advanced symplectification methods. Symplectification is important for maintaining the physical correctness of the muon beam behaviour with high precision over the storage time, ensuring conservation of phase space volume and preventing artificial damping or excitation of particle motion. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous BNL experiment. Analyses of data from Runs 4-6 are ongoing, with results planned for release in 2025, potentially resolving the current tension between experiment and theory.
This project investigated multifluid models for simulating interpenetrating plasmas. These models solve the Euler equations for each plasma stream along with an appropriate model for electron dynamics. The project developed a novel code, Euclid, that simulates these models using high-order conservative spatial discretization methods and explicit and semi-implicit time integration methods. Euclid was used to solve several multifluid test cases representative of plasma dynamics in high-energy density physics experiments. Kinetic simulations of interpenetrating flows using the LOKI code revealed novel physical effects of the ion-ion streaming instability. A new multi-flow collision operator was developed for the LOKI code. A simple model of the ion-ion streaming instability was also developed. Several presentations were given on the project results; three journal articles have been published, one has been submitted, and two are planned.
The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.
Due to the strong coupling between magnetism and ferroelectricity, (ND 4 ) 2 FeCl 5 ·D 2 O exhibits several intriguing magnetic and electric phases. In this Letter, we include higher-order onsite anisotropic spin interactions to explain the ferroelectric phase transitions of (ND 4 ) 2 FeCl 5 ·D 2 O in a magnetic field and to produce the large weights of high-order harmonic components in the cycloid structure that are observed from neutron diffraction experiments. Moreover, we predict a new ferroelectric phase sandwiched between the FE II and FE III phases in a magnetic field. Furthermore, by emphasizing the importance of the higher-order spin anisotropic interactions, our work provides a framework to understand multiferroic materials with rich phase diagrams.
We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].
Many applications in magnetic confinement fusion require the efficient calculation of surface integrals with singular integrands. The singularity subtraction approaches typically used to handle such singularities are complicated to implement and low-order accurate. In contrast, we demonstrate that the Kapur–Rokhlin quadrature scheme is well-suited for the logarithmically singular integrals encountered for a toroidally axisymmetric confinement system, is easy to implement and is high-order accurate. As an illustration, we show how to apply this quadrature scheme for the efficient and accurate calculation of the normal component of the magnetic field due to the plasma current on the plasma boundary, via the virtual-casing principle.
Higher symmetries in interacting many-body systems often give rise to new phases and unexpected dynamical behavior. Here, we theoretically investigate a variant of the Dicke model with higher-order discrete symmetry, resulting from complex-valued coupling coefficients between quantum emitters and a bosonic mode. We propose a driven-dissipative realization of this model focusing on optomechanical response of a driven atom tweezer array comprised of 𝑛 subensembles and placed within an optical cavity, with the phase of the driving field advancing stepwise between subensembles. Examining stationary points and their dynamical stability, we identify a phase diagram for 𝑛≥3 with three distinctive features: a ℤ 𝑛 (ℤ 2𝑛 ) symmetry-breaking superradiant phase for even (odd) 𝑛, a normal unbroken-symmetry phase that is dynamically unstable due to nonreciprocal forces between emitters, and a first-order phase transition separating these phases. This 𝑛-phase Dicke model may be equivalently realized in a variety of optomechanical or optomagnonic settings, where it can serve as a test bed for studying high-order symmetry breaking and nonreciprocal interactions in open systems.
The increasing penetration of distributed energy resources (DERs) highlights the growing importance of microgrids (MGs) in enhancing power system reliability. Employing electromagnetic transient (EMT) analysis in MGs becomes crucial for controlling the rapid transients. However, this requires an accurate but high-order model of power electronics and their underlying control loops, complexifying the stability analysis from the viewpoint of a higher control level. To overcome these challenges, this paper proposes a large-signal order reduction (LSOR) method for MGs with considerations of external control inputs and the detailed dynamics of underlying control levels based on singular perturbation theory (SPT). Specially, we innovatively proposed and strictly proved a general stability and accuracy assessment theorem that allows us to analyze the dynamic stability of a full-order nonlinear system by only leveraging our derived reduced-order model (ROM) and boundary layer model (BLM). Furthermore, this theorem furnishes a set of conditions that determine the accuracy of the developed ROM. Lastly, by embedding such a theorem into the SPT, we propose a novel LSOR approach with guaranteed accuracy and stability analysis equivalence. Case studies are conducted on MG systems to show the effectiveness of the proposed approach.
The self-assembly of colloidal nanoparticles into ordered superlattices typically uses dynamic interactions to govern particle crystallization, as these non-permanent bonds prevent the formation of kinetically trapped, disordered aggregates. However, while the use of reversible bonding is critical in the formation of highly ordered particle arrangements, dynamic interactions also inherently make the structures more prone to disassembly or disruption when subjected to different environmental stimuli. Thus, there is typically a trade-off between the ability to initially form an ordered colloidal material and the ability of that material to retain its order under different conditions. Here, we present a method for embedding colloidal nanoparticle superlattices into a polymer gel matrix. This encapsulation strategy physically prevents the nanoparticles from dissociating upon heating, drying, or the introduction of chemicals that would normally disrupt the lattice. However, the use of a gel as the embedding medium still permits further modification of the colloidal nanoparticle lattice by introducing stimuli that deform the gel network (as this deformation in turn alters the nanoparticle lattice structure in a predictable manner). Moreover, encapsulation of the lattice within a gel permits further stabilization into fully solid materials by removing the solvent from the gel or by replacing the solvent with a liquid monomer that can be photopolymerized. This embedding method therefore makes it possible to incorporate ordered colloidal arrays into a polymer matrix as either dynamic or static structures, expanding their potential for use in responsive materials.
Abstract Nonfused ring acceptors (NFRAs) have attracted significant attention for nonfullerene organic solar cells (OSCs) owing to their chemical tunability and facile synthesis. In this study, a benzotriazole‐based NFRA with chlorinated end groups (Triazole‐4Cl) is developed to realize highly efficient and thermally stable NFRA‐based OSCs; an analogous NFRA with nonchlorinated end groups (Triazole‐H) is synthesized for comparison. Triazole‐4Cl film exhibits the high‐order packing structure and the near‐infrared absorption capability, which are advantageous in charge transport and light harvesting of the resulting OSCs. In particular, the strong crystalline behavior of Triazole‐4Cl results in enhanced self‐aggregation, leading to high charge carrier mobility. Owing to these properties, a PBDB‐T (polymer donor):Triazole‐4Cl OSC demonstrates a high short‐circuit current, fill factor, and power conversion efficiency (PCE = 10.46%), outperforming a PBDB‐T:Triazole‐H OSC (PCE = 7.65%). In addition, the thermal stability of a PBDB‐T:Triazole‐4Cl OSC at an elevated temperature of 120 °C exceeds that of a PBDB‐T:Triazole‐H OSC. This is mainly attributed to the significantly higher cold crystallization temperature of Triazole‐4Cl (205.9 °C). This work provides useful guidelines for the design of NFRAs to achieve efficient and thermally stable NFRA‐based OSCs.
We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.
Here, this paper proposes a semi-analytical approach for efficient and accurate electromagnetic transient (EMT) simulation of a power grid. The approach first derives a high-order semi-analytical solution (SAS) of the grid’s state-space EMT model using the differential transformation (DT), and then evaluates the solution over enlarged, variable time steps to significantly accelerate the simulations while maintaining its high accuracy on detailed fast EMT dynamics. The approach also addresses switches during large time steps by using a limit violation detection algorithm with a binary search-enhanced quadratic interpolation. Case studies are conducted on EMT models of the IEEE 39-bus system and large-scale systems to demonstrate the merits of the new simulation approach against traditional numerical methods.