Search NASA⌕ Search

SEARCH · Search NASA

Results for “canonization”

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 361 records · Page 20

Fluctuations in heavy ion collisions and global conservation effects

Subensemble is a type of statistical ensemble which is the generalization of grand canonical and canonical ensembles. The subensemble acceptance method (SAM) provides general formulas to correct the cumulants of distributions in heavy-ion collisions for the global conservation of all QCD charges. The method is applicable for an arbitrary equation of state and sufficiently large systems, such as those created in central collisions of heavy ions. The new fluctuation measures insensitive to global conservation effects are presented. The main results are illustrated in the hadron resonance gas and van der Waals fluid frameworks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born–Oppenheimer approximation

In this work, we develop a formalism to accurately account for the renormalization of the electronic structure due to quantum and thermal nuclear motions within the Born–Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from quantum Monte Carlo calculations in either the canonical or grand-canonical ensembles. The formalism applies as well to effective single electron theories such as those based on density functional theory. We show that the electronic (Bloch) crystal momentum can be restored by marginalizing the total electron–ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born–Oppenheimer approximation. Based on the Kubo–Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and, in general, differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200 K and 250 GPa and discuss the optical absorption profile of hydrogen crystals at 200 K and carbon diamond at 297 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA–ellipticine and buckycatcher–fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA–ellipticine intercalation complex, and the buckycatcher–C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller−Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA–ellipticine, and buckycatcher–C60, respectively. Such sizable discrepancy above “chemical accuracy” for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree–Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Ballesteros, Francisco (ORCID:0000000325732307)↗

Constrained non-negative matrix factorization enabling real-time insights of in situ and high-throughput experiments

Non-negative matrix factorization (NMF) is an appealing class of methods for performing unsupervised learning on streaming spectral data, particularly in time-sensitive applications such as in situ characterization of materials. These methods seek to decompose a dataset into a small number of components and weights that can compactly represent the underlying signal while effectively reconstructing the observations with minimal error. However, canonical NMF methods have no underlying requirement that the reconstruction uses components or weights that are representative of the true physical processes. In this work, we demonstrate how constraining a subset of the NMF weights or components as rigid priors, provided as known or assumed values, can provide significant improvement in revealing true underlying phenomena. We present a PyTorch-based method for efficiently applying constrained NMF and demonstrate its application to several synthetic examples. Our implementation allows an expert researcher-in-the-loop to provide and dynamically adjust the constraints during a live experiment involving streaming spectral data. Such interactive priors allow researchers to specify known or identified independent components, as well as functional expectations about the mixing or transitions between the components. We further demonstrate the application of this method to measured synchrotron x-ray total scattering data from in situ beamline experiments. In such a context, constrained NMF can result in a more interpretive and scientifically relevant decomposition than canonical NMF or other decomposition techniques. As a result, the details of the method are provided, along with general guidance for employing constrained NMF in the extraction of critical information and insights during time-sensitive experimental applications.

36 MATERIALS SCIENCE↗

Magnetic field properties in non-axisymmetric divertors

Stellarator power plants require a plan for the removal of the particles and the heat that are exhausted across the plasma edge. Unless a flowing liquid metal can be used to carry the helium exhaust to places where it can be removed from the plasma chamber, the particle exhaust must be magnetically diverted into pumping chambers. Studies are required to determine how magnetic features relate to the required divertor properties, how these magnetic features can be produced, and how they can be controlled. General studies are clarified and simplified by the use of the magnetic field line Hamiltonian $ψ_p(ψ,θ,φ)$ and a vector $\vec{x}(ψ,θ,φ)$ that gives the point in space associated with each point in the $(ψ,θ,φ)$ canonical coordinates, a flux and two angles. The non-resonant Fourier terms in ψp can be removed by a canonical transformation, so only resonant Fourier terms can determine the field line properties in the plasma edge and divertor. This paper discusses the important divertor properties and explains how $ψ_p(ψ,θ,φ)$ and $\vec{x}(ψ,θ,φ)$ can be obtained numerically in a special form for any stellarator magnetic field, $\vec{B}(\vec{x})$. This form holds between an arbitrary magnetic surface and the chamber walls with the non-resonant terms eliminated. Studies based on variations in the terms in such derived field-line Hamiltonians can determine what magnetic features are mathematically possible and how they could be produced and controlled by the external magnetic field coils.

Physics↗

Benchmarking the exponential ansatz for the Holstein model

Polarons are quasiparticles formed as a result of lattice distortions induced by charge carriers. The single-electron Holstein model captures the fundamentals of single polaron physics. We examine the power of the exponential ansatz for the polaron ground-state wavefunction in its coupled cluster, canonical transformation, and (canonically transformed) perturbative variants across the parameter space of the Holstein model. Our benchmark serves to guide future developments of polaron wavefunctions beyond the single-electron Holstein model.

Chemistry↗

Theoretical assessment of the transition between electron emission mechanisms in a collisional crossed-field gap

Understanding the transition to the limiting current under vacuum and collisional conditions for thermal-field emission in a crossed-field gap is critical for characterizing the operating regime for numerous applications, such as high-power microwaves and Hall thrusters. By deriving the exact solution for the current density as a function of cathode temperature (and, by extension, initial velocity), cathode electric field, applied voltage, collision frequency, and magnetic field, we determine the transition from the canonical current source (thermal-field emission) to the vacuum and collisional limiting currents under different operating conditions. Here, we further show that the exact solution recovers the asymptotic solutions for the canonical and limiting currents for vacuum and collisional crossed-field gaps under appropriate limits. The exact solution demonstrates that the critical current decreases with both increasing magnetic field strength and increasing collisionality. Nexus phase space plots constructed by matching the individual asymptotic solutions provide the conditions for the transitions between mechanisms, elucidating when a single mechanistic equation would suffice and when theories combining multiple mechanisms are necessary.

Breen, Lorin I. [Purdue Univ., West Lafayette, IN ↗

Mapping the Temperature-dependent and network site-specific onset of spectral diffusion at the surface of a water cluster cage

We explore the kinetic processes that sustain equilibrium in a microscopic, finite system. This is accomplished by monitoring the spontaneous, time-dependent frequency evolution (the frequency autocorrelation) of a single OH oscillator, embedded in a water cluster held in a temperature-controlled ion trap. The measurements are carried out by applying two-color, IR-IR photodissociation mass spectrometry to the D3O+?(HDO)(D2O)19 isotopologue of the “magic number” protonated water cluster, H+?(H2O)21. The OH group can occupy any one of the five spectroscopically distinct sites in the distorted pentagonal dodecahedron cage structure. The OH frequency is observed to evolve over tens of milliseconds in the temperature range (90-120 K). Starting at 100 K, large “jumps” are observed between two OH frequencies separated by ~300 cm-1 indicating migration of the OH group from the bound OH site at 3350 cm-1 to the free position at 3686 cm-1. Increasing the temperature to 110 K leads to partial interconversion among many sites. All sites are observed to interconvert at 120 K such that the distribution of the unique OH group among them adopts the form one would expect for a canonical ensemble. The spectral dynamics displayed by the clusters thus offer an unprecedented view into the molecular-level processes that drive spectral diffusion in an extended network of water molecules.

Yang, Nan↗

Structure of saguaro cactus virus 3′ translational enhancer mimics 5′ cap for eIF4E binding

The genomes of several plant viruses contain RNA structures at their 3' ends called cap-independent translation enhancers (CITEs) that bind the host protein factors such as mRNA 5' cap-binding protein eIF4E for promoting cap-independent genome translation. However, the structural basis of such 5' cap-binding protein recognition by the uncapped RNA remains largely unknown. Here, we have determined the crystal structure of a 3' CITE, panicum mosaic virus-like translation enhancer (PTE) from the saguaro cactus virus (SCV), using a Fab crystallization chaperone. The PTE RNA folds into a three-way junction architecture with a pseudoknot between the purine-rich R domain and pyrimidine-rich Y domain, which organizes the overall structure to protrude out a specific guanine nucleotide, G18, from the R domain that comprises a major interaction site for the eIF4E binding. The superimposable crystal structures of the wild-type, G18A, G18C, and G18U mutants suggest that the PTE scaffold is preorganized with the flipped-out G18 ready to dock into the eIF4E 5' cap-binding pocket. The binding studies with wheat and human eIF4Es using gel electrophoresis and isothermal titration calorimetry, and molecular docking computation for the PTE–eIF4E complex demonstrated that the PTE structure essentially mimics the mRNA 5' cap for eIF4E binding. Such 5' cap mimicry by the uncapped and structured viral RNA highlights how viruses can exploit RNA structures to mimic the host protein-binding partners and bypass the canonical mechanisms for their genome translation, providing opportunities for a better understanding of virus-host interactions and non-canonical translation mechanisms found in many pathogenic RNA viruses.

3' cap-independent translation enhancers↗

Insights into an unusual Auxiliary Activity 9 family member lacking the histidine brace motif of lytic polysaccharide monooxygenases

Lytic polysaccharide monooxygenases (LPMOs) are redox-enzymes involved in biomass degradation. All characterized LPMOs possess an active site of two highly conserved histidine residues coordinating a copper ion (the histidine brace), which are essential for LPMO activity. However, some protein sequences that belong to the AA9 LPMO family display a natural N-terminal His to Arg substitution (Arg-AA9). These are found almost entirely in the phylogenetic fungal class Agaricomycetes, associated with wood decay, but no function has been demonstrated for any Arg-AA9. Through bioinformatics, transcriptomic, and proteomic analyses we present data, which suggest that Arg-AA9 proteins could have a hitherto unidentified role in fungal degradation of lignocellulosic biomass in conjunction with other secreted fungal enzymes. We present the first structure of an Arg-AA9, LsAA9B, a naturally occurring protein from Lentinus similis The LsAA9B structure reveals gross changes in the region equivalent to the canonical LPMO copper-binding site, whereas features implicated in carbohydrate binding in AA9 LPMOs have been maintained. We obtained a structure of LsAA9B with xylotetraose bound on the surface of the protein although with a considerably different binding mode compared with other AA9 complex structures. In addition, we have found indications of protein phosphorylation near the N-terminal Arg and the carbohydrate-binding site, for which the potential function is currently unknown. Finally, our results are strong evidence that Arg-AA9s function markedly different from canonical AA9 LPMO, but nonetheless, may play a role in fungal conversion of lignocellulosic biomass.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural insight into metallocofactor maturation in carbon monoxide dehydrogenase

The nickel-dependent carbon monoxide dehydrogenase (CODH) employs a unique heterometallic nickel–iron–sulfur cluster, termed the C-cluster, to catalyze the interconversion of CO and CO 2 . Like other complex metalloenzymes, CODH requires dedicated assembly machinery to form the fully intact and functional C-cluster. In particular, nickel incorporation into the C-cluster depends on the maturation factor CooC; however, the mechanism of nickel insertion remains poorly understood. Here, we compare X-ray structures (1.50–2.48 Å resolution) of CODH from Desulfovibrio vulgaris ( Dv CODH) heterologously expressed in either the absence ( Dv CODH –CooC ) or presence ( Dv CODH +CooC ) of co-expressed CooC. Here, we find that the C-cluster of Dv CODH –CooC is fully loaded with iron but does not contain any nickel. Interestingly, the so-called unique iron ion (Fe u ) occupies both its canonical site (80% occupancy) and the nickel site (20% occupancy), with addition of reductant causing further mismetallation of the nickel site (60% iron occupancy). We also demonstrate that a Dv CODH variant that lacks a surface-accessible iron–sulfur cluster (the D-cluster) has a C-cluster that is also replete in iron but lacks nickel, despite co-expression with CooC. In this variant, all Fe u is in its canonical location, and the nickel site is empty. This D-cluster–deficient CODH is inactive despite attempts to reconstitute it with nickel. Taken together, these results suggest that an empty nickel site is not sufficient for nickel incorporation. Based on our findings, we propose a model for C-cluster assembly that requires both CooC and a functioning D-cluster, involves precise redox-state control, and includes a two-step nickel-binding process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Field theory and a structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence

A field theory and the associated structure-preserving geometric particle-in-cell (PIC) algorithm are developed to study low frequency electrostatic perturbations with fully kinetic ions and adiabatic electrons in magnetized plasmas. The algorithm is constructed by geometrically discretizing the field theory using discrete exterior calculus, high-order Whitney interpolation forms, and non-canonical Hamiltonian splitting method. The discretization preserves the non-canonical symplectic structure of the particle-field system, as well as the electromagnetic gauge symmetry. As a result, the algorithm is charge-conserving and possesses long-term conservation properties. Because drift wave turbulence and anomalous transport intrinsically involve multi time-scales, simulation studies using fully kinetic particles demand algorithms with long-term accuracy and fidelity. The structure-preserving geometric PIC algorithm developed adequately serves this purpose. The algorithm has been implemented in the code, tested and benchmarked using the examples of ion Bernstein waves and drift waves. We apply the algorithm to study the ion temperature gradient (ITG) instability and turbulence in a 2D slab geometry. Simulation results show that at the early stage of the turbulence, the energy diffusion is between the Bohm scaling and gyro-Bohm scaling. Finally, at later time, the observed diffusion is closer to the gyro-Bohm scaling, and density blobs generated by the rupture of unstable modes are the prominent structures of the fully developed ITG turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ACTINN: automated identification of cell types in single cell RNA sequencing

Abstract Motivation Cell type identification is one of the major goals in single cell RNA sequencing (scRNA-seq). Current methods for assigning cell types typically involve the use of unsupervised clustering, the identification of signature genes in each cluster, followed by a manual lookup of these genes in the literature and databases to assign cell types. However, there are several limitations associated with these approaches, such as unwanted sources of variation that influence clustering and a lack of canonical markers for certain cell types. Here, we present ACTINN (Automated Cell Type Identification using Neural Networks), which employs a neural network with three hidden layers, trains on datasets with predefined cell types and predicts cell types for other datasets based on the trained parameters. Results We trained the neural network on a mouse cell type atlas (Tabula Muris Atlas) and a human immune cell dataset, and used it to predict cell types for mouse leukocytes, human PBMCs and human T cell sub types. The results showed that our neural network is fast and accurate, and should therefore be a useful tool to complement existing scRNA-seq pipelines. Availability and implementation The codes and datasets are available at https://figshare.com/articles/ACTINN/8967116. Tutorial is available at https://github.com/mafeiyang/ACTINN. All codes are implemented in python. Supplementary information Supplementary data are available at Bioinformatics online.

Ma, Feiyang (ORCID:0000000262600787)↗

Non-universal stellar initial mass functions: large uncertainties in star formation rates at z ≈ 2–4 and other astrophysical probes

ABSTRACT We explore the assumption, widely used in many astrophysical calculations, that the stellar initial mass function (IMF) is universal across all galaxies. By considering both a canonical broken-power-law IMF and a non-universal IMF, we are able to compare the effect of different IMFs on multiple observables and derived quantities in astrophysics. Specifically, we consider a non-universal IMF that varies as a function of the local star formation rate, and explore the effects on the star formation rate density (SFRD), the extragalactic background light, the supernova (both core-collapse and thermonuclear) rates, and the diffuse supernova neutrino background. Our most interesting result is that our adopted varying IMF leads to much greater uncertainty on the SFRD at $z \approx 2-4$ than is usually assumed. Indeed, we find an SFRD (inferred using observed galaxy luminosity distributions) that is a factor of $\gtrsim 3$ lower than canonical results obtained using a universal IMF. Secondly, the non-universal IMF we explore implies a reduction in the supernova core-collapse rate of a factor of $\sim 2$, compared against a universal IMF. The other potential tracers are only slightly affected by changes to the properties of the IMF. We find that currently available data do not provide a clear preference for universal or non-universal IMF. However, improvements to measurements of the star formation rate and core-collapse supernova rate at redshifts $z \gtrsim 2$ may offer the best prospects for discernment.

79 ASTRONOMY AND ASTROPHYSICS↗

How to optimally constrain galaxy assembly bias: supplement projected correlation functions with count-in-cells statistics

Most models for the statistical connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than halo mass, a phenomenon known as galaxy assembly bias. And yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data that are analysed using traditional halo occupation models. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy-halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function w p ( r p ), the galaxy-galaxy lensing signal ΔΣ( r p ), the void probability function VPF( r ), the distributions of counts-in-cylinders P(NCIC), and counts-in-annuli P(NCIA), and the distribution of the ratio of counts in cylinders of different sizes P ( N 2 /N 5 ). We find that despite the frequent use of the combination w p ( r p ) + ΔΣ( r p ) in interpreting galaxy data, the count statistics, P ( N CIC ) and P ( N CIA ), are generally more efficient in constraining galaxy assembly bias when combined with w p ( r p ). Constraints based upon w p ( r p ) and ΔΣ(rp) share common degeneracy directions in the parameter space, while combinations of r p w p ( r p ) with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy-halo connection.

galaxies: evolution↗

DXO/Rai1 enzymes remove 5′-end FAD and dephospho-CoA caps on RNAs

Abstract In eukaryotes, the DXO/Rai1 enzymes can eliminate most of the incomplete and non-canonical NAD caps through their decapping, deNADding and pyrophosphohydrolase activities. Here, we report that these enzymes can also remove FAD and dephospho-CoA (dpCoA) non-canonical caps from RNA, and we have named these activities deFADding and deCoAping. The crystal structures of mammalian DXO with 3′-FADP or CoA and fission yeast Rai1 with 3′-FADP provide elegant insight to these activities. FAD and CoA are accommodated in the DXO/Rai1 active site by adopting folded conformations. The flavin of FAD and the pantetheine group of CoA contact the same region at the bottom of the active site tunnel, which undergoes conformational changes to accommodate the different cap moieties. We have developed FAD-capQ to detect and quantify FAD-capped RNAs and determined that FAD caps are present on short RNAs (with less than ∼200 nucleotides) in human cells and that these RNAs are stabilized in the absence of DXO.

59 BASIC BIOLOGICAL SCIENCES↗

Uncovering translation roadblocks during the development of a synthetic tRNA

Ribosomes are remarkable in their malleability to accept diverse aminoacyl-tRNA substrates from both the same organism and other organisms or domains of life. This is a critical feature of the ribosome that allows the use of orthogonal translation systems for genetic code expansion. Optimization of these orthogonal translation systems generally involves focusing on the compatibility of the tRNA, aminoacyl-tRNA synthetase, and a non-canonical amino acid with each other. As we expand the diversity of tRNAs used to include non-canonical structures, the question arises as to the tRNA suitability on the ribosome. Specifically, we investigated the ribosomal translation of allo-tRNA UTu1 , a uniquely shaped (9/3) tRNA exploited for site-specific selenocysteine insertion, using single-molecule fluorescence. With this technique we identified ribosomal disassembly occurring from translocation of allo-tRNA UTu1 from the A to the P site. Using cryo-EM to capture the tRNA on the ribosome, we pinpointed a distinct tertiary interaction preventing fluid translocation. Through a single nucleotide mutation, we disrupted this tertiary interaction and relieved the translation roadblock. With the continued diversification of genetic code expansion, our work highlights a targeted approach to optimize translation by distinct tRNAs as they move through the ribosome.

59 BASIC BIOLOGICAL SCIENCES↗