Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coordinate transformations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

NXREFINE

NXRefine implements a complete workflow for both data acquisition and reduction of single crystal diffuse x-ray scattering collected as a series of area detector frames during continuous sample rotation. The data are stored in HDF5 files, which conform to the NeXus data format standard. These files contain a comprehensive set of metadata, including the incident wavelength, powder calibration scans used to define detector distances, beam centers, and yaw, pitch, and roll corrections. Automated peak searches of the three-dimensional data arrays define a set of Bragg peaks that are used to define an orientation matrix, defined according to Busing and Levy, which is then used to transform the data from angular coordinates to reciprocal space coordinates using the CCTW program. Multiple rotations at different detector translations and offset rotation axes are combined to ensure that there are no gaps in the reconstructed data, and allow the construction of masks to eliminate detector artifacts due to Compton scattering within the sensor layers. Finally, the data can be transformed into 3D-ΔPDF maps using the punch-and-fill method after symmetrization.

OSBORN, RAYMOND↗
DOE OSTI2023

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Compartmentalization and Photo-Regulating Pathways for Incompatible Tandem Catalysis

This contribution describes an advanced compartmentalized micellar nanoreactor that possesses a reversible photoresponsive feature and its application towards photo-regulating reaction pathways for incompatible tandem catalysis under aqueous conditions. The smart nanoreactor is based on multi-functional amphiphilic poly(2-oxazoline)s and covalently crosslinked with spiropyran upon micelle formation in water. It responds to light irradiation in a wavelength-selective manner switching its morphology as confirmed by dynamic light scattering and cryo-transition electron microscopy. The compartmental structure renders distinct nanoconfinements for two incompatible enantioselective transformations: a rhodium-diene complex-catalyzed asymmetric 1,4-addition occurs in the hydrophilic corona while a Rh-TsDPEN-catalyzed asymmetric transfer hydrogenation proceeds in the hydrophobic core. Control experiments and kinetic studies showed that the gated behavior induced by the photo-triggered reversible spiropyran to merocyanine transition in the cross-linking layer is key to discriminate among substrates/reagents during the catalysis. Here, the smart nanoreactor realized photo-regulation to direct the reaction pathway to give a multi-chiral product with high conversions and perfect enantioselectivities in aqueous media. Our SCM catalytic system, on a basic level, mimics the concepts of compartmentalization and responsiveness Nature uses to coordinate thousands of incompatible chemical transformations into streamlined metabolic processes.

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

Sulfated Zirconium Metal–Organic Frameworks as Well-Defined Supports for Enhancing Organometallic Catalysis

Understanding heterogeneous catalysts is a challenging pursuit due to surface site nonuniformity and aperiodicity in traditionally-used materials. One example is sulfated metal oxides, which function as highly-active catalysts and as supports for organometallic complexes. These applications are due to traits such as acidity, ability to act as a weakly coordinating ligand, and aptitude for promoting transformations via radical cation intermediates. Research is ongoing about the structural features of sulfated metal oxides that imbue the aforementioned properties, such as sulfate geometry and coordination. To better understand these materials, metal–organic frameworks (MOFs) have been targeted as structurally-defined analogs. Composed of inorganic nodes and organic linkers, MOFs possess features such as high porosity and crystallinity, which make them attractive for mechanistic studies of heterogeneous catalysts. In this work, a Zr 6 -based MOF NU-1000 is sulfated and characterized using atomically-precise techniques such as single crystal X-ray diffraction (SCXRD) in addition to diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). Here, the dynamic nature of the sulfate binding motif is found to transition from monodentate, to bidentate, to tridentate depending on the degree of hydration, as supported by density functional theory (DFT) calculations. Heightened Brønsted acidity compared to the parent MOF was observed upon sulfation, and probed through trimethylphosphine oxide (TMPO) physisorption, ammonia sorption, in-situ ammonia DRIFTS, and DFT studies. With the support structure benchmarked, an organoiridium complex was chemisorbed onto the sulfated MOF node and the efficacy of this supported catalyst was demonstrated for stoichiometric and catalytic activation of benzene-d6 and toluene with structure-activity relationships derived.

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

Bundling measures for food systems transformation: a global, multimodel assessment

Background Current food systems leave one in ten individuals at risk of hunger while driving unsustainable environmental impacts. Inaction risks further exacerbating negative impacts on both human and planetary health. These challenges emerge from complex system interactions, requiring approaches that engage with this complexity and consider how transformation measures interact across food systems. We aimed to quantify the magnitude and uncertainty of the impacts of key food systems transformation measures both individually and in a bundle using an ensemble of global economic models. Methods In this global multimodel assessment, we applied an ensemble of ten state-of-the-art global economic models to evaluate the potential of four key measures in transforming food systems: increasing agricultural productivity, halving food loss and waste, shifting towards healthier diets, and economy-wide climate mitigation policies aligned with limiting warming to 1·5°C. The scenarios used a middle-of-the-road shared socioeconomic pathway for population and gross domestic product growth, climate impact data from Jägermeyr and colleagues, Thornton and colleagues, and Nelson and colleagues, and dietary targets based on the EAT–Lancet healthy reference diet, with model simulations conducted from 2020 to 2050. We then assessed the effect of these measures in isolation and in combination in a bundled scenario. To further understand the interactions between these measures, we conducted a decomposition analysis that distinguishes between the individual effects of a measure (effect when implemented alone), total effects (its contribution within the bundle), and interaction effects (the difference between total and individual effects). This approach aimed to show complementarities and trade-offs that emerge when multiple measures are implemented simultaneously. Findings Our analysis showed that individual measures in isolation are insufficient to achieve high-level environmental objectives and might generate unintended consequences. In contrast, bundling measures produces co-benefits: avoiding 50% of projected agricultural greenhouse gas emissions by 2050 and almost 20% of anticipated land conversion, while moderating food price increases associated with ambitious climate change mitigation policies. Our decomposition analysis further shows that measures can have varying effects across different dimensions. Although dietary shifts and climate mitigation policies are the largest drivers of environmental benefits (each contributing to a median decline of >10 percentage points in non-CO 2 emissions and 5 percentage points in agricultural land use globally), productivity improvements and reducing food loss and waste play essential roles in moderating price increases (each contributing to a median decline of >5 percentage points in average prices). Interpretation This study highlights the importance of implementing coordinated approaches to food system transformation and climate change mitigation rather than relying on isolated interventions. Comprehensive transformation requires understanding how supply-side and demand-side changes can interact with climate mitigation policies, enabling policy makers to design intervention packages that maximise benefits while minimising trade-offs across environmental, economic, and social dimensions.

Sundiang, Marina [Cornell Univ., Ithaca, NY (Unite↗
DOE OSTI2026

Hollow-structured Ni-N-C catalysts for highly selective CO 2 electroreduction

Atomically dispersed single-atom catalysts have emerged as promising non-precious catalyst alternatives to expensive Ag and Au catalysts for electrochemical CO 2 reduction reaction (CO 2 RR). In particular, nickel-nitrogen-carbon (Ni-N-C) catalysts have demonstrated a high faradaic efficiency (FE) toward CO formation at low overpotentials. Nonetheless, the exact nature of Ni active sites under CO 2 RR remains elusive and conventional Ni-N-C catalysts are limited by microporosity and low density of Ni single atoms, hindering performance in CO 2 electrolyzers. Here, we report the synthesis of hollow-structured Ni-N-C ( hs -Ni-N-C) catalysts via a post-synthesis modification (PSM) strategy using partial ligand exchange of 2-methylimidazole with 3-amino-1,2,4-triazole. This approach enables the formation of a hollow structure, resulting in more than a twofold increase in Ni atom density compared to regular Ni-N-C (r-Ni-N-C). In a zero-gap CO 2 electrolyzer, the optimized hs -Ni-N-C allows for achieving an FE CO of 97% at a current density of > 100 mA cm⁻ 2 , while maintaining high CO selectivity with stable performance over 100 h at 2.5 V. hs-Ni-N-C shows a more than sevenfold increase in the CO partial current density relative to r-Ni-N-C resulting from the combined effects of a higher density of Ni single-atom sites, improved kinetics, and lower transport resistance under the operating conditions, as indicated by electrochemical impedance spectra and distribution of relaxation times analysis. Operando high energy-resolution X-ray absorption spectroscopy (XAS) reveals that atop-bonded CO on Ni single sites induces dynamic transformations of the Ni–N coordination environment, leading to a symmetric coordination structure of hs -Ni-N-C. Under CO 2 RR, the catalysts undergo a more pronounced structural change and form a minor fraction of Ni nanoparticles. Density functional theory calculations are consistent with the XAS results and provide molecular insights showing that the interplay between protonation and CO adsorption leads to adsorbate-induced restructuring of the Ni single atom. This work demonstrates the synergistic role of hollow structure and high-density Ni atoms in governing CO 2 RR selectivity and provides mechanistic insights into the structural dynamics of single-atom catalysts under operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Solution of the linear wave-particle kinetic equation for global modes of arbitrary frequency in a tokamak

The linear response of a plasma to perturbations of arbitrary frequency and wavelength is derived for any axisymmetric magnetized toroidal plasma. An explicit transformation to action-angle coordinates is achieved using orthogonal magnetic coordinates and the Littlejohn Lagrangian, establishing the validity of this result to arbitrary order in normalized Larmor radius. The global resonance condition for compressional modes is clarified in more detail than in previous works, confirming that the poloidal orbit-average of the cyclotron frequency gives the desired result at lowest order in Larmor radius. The global plasma response to the perturbation at each resonance is captured by a poloidal and gyroaverage of the perturbing potential. A “global gyroaveraging” of the potential is a natural by-product of this analysis which takes into account the changing of the magnetic field over an orbit. The resonance condition depends on two arbitrary integers which completely separately capture the effects poloidal non-uniformity and finite Larmor radius in generating sidebands. We learn that poloidal sidebands generated for compressional modes are dominated by the change in gyrofrequency over the orbit, which is very different to shear modes where the gyrofrequency only contributes via a finite Larmor radius effect. This increases the number of bounce harmonics required to compute the linear drive, giving a more complicated resonance map. An example calculation is given comparing resonance of shear and compressional modes in a published DIII-D case.

Compressional↗
DOE OSTI2020

Comparing light-front quantization with instant-time quantization

Here in this paper we compare light-front quantization and instant-time quantization both at the level of operators and at the level of their Feynman diagram matrix elements. At the level of operators light-front quantization and instant-time quantization lead to equal light-front time commutation (or anticommutation) relations that appear to be quite different from equal instant-time commutation (or anticommutation) relations. Despite this we show that at unequal times instant-time and light-front commutation (or anticommutation) relations actually can be transformed into each other, with it only being the restriction to equal times that makes the commutation (or anticommutation) relations appear to be so different. While our results are valid for both bosons and fermions, for fermions there are subtleties associated with tip of the light cone contributions that need to be taken care of. At the level of Feynman diagrams we show for non-vacuum Feynman diagrams that the pole terms in four-dimensional light-front Feynman diagrams reproduce the widely used three-dimensional light-front on-shell Hamiltonian Fock space formulation in which the light-front energy and light-front momentum are on shell. Moreover, we show that the contributions of pole terms in non-vacuum instant-time and non-vacuum light-front Feynman diagrams are equal. However, because of circle at infinity contributions we show that this equivalence of pole terms fails for four-dimensional light-front vacuum tadpole diagrams. Then, and precisely because of these circle at infinity contributions, we show that light-front vacuum tadpole diagrams are not only nonzero, they quite remarkably are actually equal to the pure pole term instant-time vacuum tadpole diagrams. Light-front vacuum diagrams are not correctly describable by the on-shell Hamiltonian formalism, and thus not by the closely related infinite momentum frame prescription either. Thus for the light-front vacuum sector we must use the off-shell Feynman formalism as it contains information that is not accessible in the on-shell Hamiltonian Fock space approach. We show that light-front quantization is intrinsically nonlocal, and that for fermions this nonlocality is present in Ward identities. One can project fermion spinors into so-called good and bad components, and both of these components contribute in Ward identities. Central to our analysis is that the transformation from instant-time coordinates and fields to light-front coordinates and fields is a unitary, spacetime-dependent translation. Consequently, not only are instant-time quantization and light-front quantization equivalent, because of general coordinate invariance they are unitarily equivalent.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The magnetic structure of long-wavelength magnetohydrodynamic modes in current-carrying stellarator plasmas

Magnetic field fluctuations are observed in current-carrying stellarator plasmas when the rotational transform is close to a rational value at the edge of the plasma. At low plasma pressure, these fluctuations are associated with perturbed currents parallel to the equilibrium magnetic field lines. A model for these magnetohydrodynamic modes in a low-β, three-dimensional stellarator equilibria has been developed. A set of helical current filaments are constrained to mimic the structure of magnetic field lines on rational surfaces derived from three-dimensional (3D) equilibrium reconstructions. Transformation to straight field line coordinates then allows fitting of the poloidal magnetic sensor data to a single harmonic function, which fixes the modeled toroidal mode structure via the field line flow geometry. The developed procedure accurately captures phase and amplitude variation for m/n = 3/2, 3/1, and 4/1 modes in the 3D equilibria of the compact toroidal hybrid experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Novel angular velocity estimation technique for plasma filaments

Magnetic field aligned filaments such as blobs and edge localized mode filaments carry significant amounts of heat and particles to the plasma facing components and they decrease their lifetime. The dynamics of these filaments determine at least a part of the heat and particle loads. These dynamics can be characterized by their translation and rotation. In this paper, we present an analysis method novel for fusion plasmas, which can estimate the angular velocity of the filaments on frame-by-frame time resolution. After pre-processing, the frames are two-dimensional (2D) Fourier-transformed, then the resulting 2D Fourier magnitude spectra are transformed to log-polar coordinates, and finally the 2D cross-correlation coefficient function (CCCF) is calculated between the consecutive frames. The displacement of the CCCF’s peak along the angular coordinate estimates the angle of rotation of the most intense structure in the frame. Further, the proposed angular velocity estimation method is tested and validated for its accuracy and robustness by applying it to rotating Gaussian-structures. The method is also applied to gas-puff imaging measurements of filaments in National Spherical Torus Experiment plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Disordered interfaces of alkaline aluminate salt hydrates provide glimpses of Al 3+ coordination changes

Hypothesis: The precipitation and dissolution of aluminum-bearing mineral phases in aqueous systems often proceed via changes in both aluminum coordination number and connectivity, complicating molecular-scale interpretation of the transformation mechanism. Here, the thermally induced transformation of crystalline sodium aluminum salt hydrate, a phase comprised of monomeric octahedrally coordinated aluminate which is of relevance to industrial aluminum processing, has been studied. Because intermediate aluminum coordination states during melting have not previously been detected, it is hypothesized that the transition to lower coordinated aluminum ions occurs within a highly disordered quasi-two-dimensional phase at the solid-solution interface. Experiments and simulations: In this work, in situ X-ray diffraction (XRD), Raman and 27 Al nuclear magnetic resonance (NMR) spectroscopy were used to monitor the melting transition of nonasodium aluminate hydrate (NSA, Na 9 [Al(OH) 6 ] 2 ·3(OH)·6H 2 O). A mechanistic interpretation was developed based on complementary classical molecular dynamics (CMD) simulations including enhanced sampling. A reactive forcefield was developed to bridge speciation in the solution and in the solid phase. Findings: In contrast to classical dissolution, aluminum coordination change proceeds through a dynamically stabilized ensemble of intermediate states in a disordered layer at the solid-solution interface. In both melting and dissolution of NSA, octahedral, monomeric aluminum transition through an intermediate of pentahedral coordination. The intermediate dehydroxylates to form tetrahedral aluminate (Al(OH)$_4^–$) in the liquid phase. This coordination change is concomitant with a breaking of the ionic aluminate-sodium ion linkages. The solution phase Al(OH)$_4^–$ ions subsequently polymerize into polynuclear aluminate ions. However, there are some differences between bulk melting and interfacial dissolution, with the onset of the surface-controlled process occurring at a lower temperature (~30 °C) and the coordination change taking place more gradually as a function of temperature. This work to determine the local structure and dynamics of aluminum in the disordered layer provides a new basis to understand mechanisms controlling aluminum phase transformations in highly alkaline solutions.

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

An extensible vector toolkit and parts library for advanced engineering of plant genomes

Abstract Plant biotechnology is rife with new advances in transformation and genome engineering techniques. A common requirement for delivery and coordinated expression in plant cells, however, places the design and assembly of transformation constructs at a crucial juncture as desired reagent suites grow more complex. Modular cloning principles have simplified some aspects of vector design, yet many important components remain unavailable or poorly adapted for rapid implementation in biotechnology research. Here, we describe a universal Golden Gate cloning toolkit for vector construction. The toolkit chassis is compatible with the widely accepted Phytobrick standard for genetic parts, and supports assembly of arbitrarily complex T‐DNAs through improved capacity, positional flexibility, and extensibility in comparison to extant kits. We also provision a substantial library of newly adapted Phytobricks, including regulatory elements for monocot and dicot gene expression, and coding sequences for genes of interest such as reporters, developmental regulators, and site‐specific recombinases. Finally, we use a series of dual‐luciferase assays to measure contributions to expression from promoters, terminators, and from cross‐cassette interactions attributable to enhancer elements in certain promoters. Taken together, these publicly available cloning resources can greatly accelerate the testing and deployment of new tools for plant engineering.

42 ENGINEERING↗
DOE OSTI2020

A discrete integral transform for rapid spectral synthesis

Accurate synthetic spectra that rely on large Line-By-Line (LBL)-databases are used in a wide range of applications such as high temperature combustion, atmospheric re-entry, planetary surveillance and laboratory plasmas. Conventionally synthetic spectra are calculated by computing a lineshape for every spectral line in the database and adding those together, which may take multiple hours for large databases. In this paper we propose a new approach for spectral synthesis based on an integral transform: the synthetic spectrum is calculated as the integral over the product of a Voigt profile and a newly proposed three-dimensional “lineshape distribution function”, which is a function of spectral position and Gaussian- & Lorentzian width coordinates. A fast discrete version of this transform based on the Fast Fourier Transform (FFT) is proposed, which improves performance compared to the conventional approach by several orders of magnitude while maintaining accuracy. Strategies that minimize the discretization error are discussed. A Python implementation of the method is compared against state-of-the-art spectral code RADIS, and is since adopted as RADIS's default synthesis method. The synthesis of a benchmark CO2 spectrum consisting of 1.8 M spectral lines and 200k spectral points took only 3.1 s using the proposed method (1011 lines × spectral points/s), a factor ~300 improvement over the state-of-the-art, with the relative improvement generally increasing for higher number of lines and/or number of spectral points. Finally, an experimental GPU-implementation of the method was also benchmarked, which demonstrated another 2~3 orders performance increase, achieving up to 5 ∙ 10 14 lines × spectral points/s.

42 ENGINEERING↗
DOE OSTI2023

Faster-than-real-time Simulation with Demonstration for Resilient DER Integration

The US electric grid is facing operational, stability, and security challenges. Transmission system operators need some measure of visibility into distribution system renewable generation. Distribution system generation needs to support transmission system voltage. The grid is experiencing an expansion in measurement systems. How to take full advantage of this expansion and defend against attacks, both cyber and physical, poses additional challenges. The Faster-than-real-time Simulation with demonstration for Resilient DER Integration project set out to do the following: a. Flatten the voltage profile through the feeders and system for cost saving and voltage stabilization needs. b. Increase the amount of intermittent distributed energy resources (IDERs) that could be deployed on a utility feeder and provide 100% or more energy needed for the demands on that feeder, and c. based on an accurate model (Digital Twin) of the utilities system, be able to detect any abnormalities on the utilities distribution system. To manage the voltage and increase IDER penetration (a,b), Graph Trace Analysis is employed in a time-series, optimal power flow to coordinate the time-varying feedback control setpoints of a distribution feeder’s utility control devices. Under the coordinated control are a Load Tap Changing Transformer, a voltage regulator, and five switched capacitor banks. The feeder serves over 2000 customers, the feeder secondaries are modeled, and the feeder has 2.3 MW of PV generation, corresponding to a 17.4% penetration of PV generation. The feeder model has over 12,000 components, where every customer load bus and PV generator are modeled. The accuracy of the power flow solution is compared against historical meter voltage measurements, the improvement in conservation voltage reduction energy savings as a function of the coordinated control desired voltage profile is investigated, and the increase in PV penetration of the coordinated control over the existing control is presented. To achieve improved control performance while observing system operation constraints, bellwether Advanced Metering Infrastructure (AMI) voltage measurements are used to adjust the desired voltage profile used by the optimal power flow analysis. To detect and alleviate or negate attacks or failures on the distribution and transmission utility grids (c) the grid needs to be resilient and self-healing. In this project software was designed to do just that. At the center of the software is an Integrated System Model (ISM) that spans from transmission to secondary distribution. The ISM is employed in real-time abnormality detection, voltage stability forecasting, and multi-mode control. Testing results are presented for: 1—attacks on utility infrastructure; 2—energy savings from optimal control; 3—distribution system control response during a low voltage transmission system event; 4—cyber-attacks on PV inverters, where physical inverters are used in hard-ware-in-the-simulation-loop studies. Contributions of this work include real-time analysis that spans from three-phase transmission through secondary distribution; an approach for detecting abnormalities that employs measurements from three independent measurement systems; and a multi-mode distribution system control that responds to cyber-attacks, physical attacks, equipment failures, and transmission system needs.

Integrated System Model, Graph Trace Analysis, Adv↗
DOE OSTI2025

Extending the Mott–Gurney law to one-dimensional nonplanar diodes using point transformations

Recent studies have applied variational calculus, conformal mapping, and point transformations to generalize the one-dimensional (1D) space-charge limited current density (SCLCD) and electron emission mechanisms to nonplanar geometries; however, these assessments have focused on extending the Child–Langmuir law (CLL) for SCLCD in vacuum. Since the charge in the diode is independent of the coordinate system (i.e., covariant), we apply bijective point transformations to extend the Mott–Gurney law (MGL) for the SCLCD in a collisional or semiconductor gap to nonplanar 1D geometries. This yields a modified MGL that replaces the Cartesian gap distance with a canonical gap distance that may be written generally in terms of geometric scale factors that are known for multiple geometries. We tabulate results for common geometries. Such an approach may be applied to any current density, including non-space-charge limited gaps and SCLCD that may fall between the CLL and MGL.

Garner, Allen L. [Purdue Univ., West Lafayette, IN↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Coordination environment tuning of nickel sites by oxyanions to optimize methanol electro-oxidation activity

To achieve zero-carbon economy, advanced anode catalysts are desirable for hydrogen production and biomass upgrading powered by renewable energy. Ni-based non-precious electrocatalysts are considered as potential candidates because of intrinsic redox attributes, but in-depth understanding and rational design of Ni site coordination still remain challenging. Here, we perform anodic electrochemical oxidation of Ni-metalloids (NiP x , NiS x , and NiSe x ) to in-situ construct different oxyanion-coordinated amorphous nickel oxyhydroxides (NiOOH-TO x ), among which NiOOH-PO x shows optimal local coordination environment and boosts electrocatalytic activity of Ni sites towards selective oxidation of methanol to formate. Experiments and theoretical results demonstrate that NiOOH-PO x possesses improved adsorption of OH* and methanol, and favors the formation of CH 3 O* intermediates. The coordinated phosphate oxyanions effectively tailor the d band center of Ni sites and increases Ni-O covalency, promoting the catalytic activity. This study provides additional insights into modulation of active-center coordination environment via oxyanions for organic molecules transformation.

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

Theory, observations, and simulations of kinetic entropy in a magnetotail electron diffusion region

In this work, we examine velocity-space kinetic entropy, a spatially local measure of entropy for systems out of thermal equilibrium, during an encounter of an electron diffusion region at a magnetic reconnection site in Earth’s magnetotail by the Magnetospheric Multiscale (MMS) mission. We start by generalizing the theory of kinetic entropy to the case of non-uniform velocity space grids and transforming the equations into spheri- cal energy coordinates useful to experimental plasma detectors. The theory is then applied to MMS data and compared to particle-in-cell simulations of reconnection. We demonstrate that the entropy-based non-Maxwellianity measure from the MMS data is of sufficiently high precision to reliably identify non-Maxwellian distributions and therefore the measurements when kinetic effects are most significant. By comparing two different non-Maxwellian measures, we show that total entropy density suffers from “information loss” because it lacks a dependence on the velocity space grid, and so has lost information about how well a distribution function is resolved. Local velocity-space kinetic entropy density recovers this information. We quantify information loss and argue that the considerations needed to minimize it are crucial for instruments designed to measure distribution functions in situ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗