Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular diagnostics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗
DOE OSTI2021

A CO-to-H{sub 2} Ratio of ≈10{sup -5} toward the Herbig Ae Star HK Ori

Measurements of gas mass in protoplanetary gas disks form the basis for estimating the conditions of planet formation. Among the most important constraints derived from disk diagnostics are the abundances of gas-phase species critical for understanding disk chemistry. Toward this end, we present direct line-of-sight measurements of H{sub 2} and CO, employing UV absorption spectroscopy from Hubble Space Telescope Cosmic Origins Spectrograph to characterize disk composition, molecular excitation temperatures, and spatial distribution in the circumstellar material around the Herbig Ae stars HK Ori and T Ori. We observe strong CO (N(CO) = 10{sup 15.5} cm{sup −2}; T {sub rot}(CO) = 19 K) and H{sub 2} (N(H{sub 2}) = 10{sup 20.34} cm{sup −2}; T {sub rot}(H{sub 2}) = 141 K) absorption toward HK Ori with a CO/H{sub 2} ratio of (≡N(CO)/N(H{sub 2})) = 1.3{sub −0.7}{sup +1.6} × 10{sup −5}. These measurements place direct empirical constraints on the CO-to-H{sub 2} conversion factor in the disk around a Herbig Ae star for the first time, although there is uncertainty concerning the exact viewing geometry of the disk. The spectra of T Ori show CO (N(CO) = 10{sup 14.9} cm{sup −2}; T {sub rot}(CO) = 124 K) absorption. Interestingly, we do not detect any H{sub 2} absorption toward this star (N(H{sub 2}) < 10{sup 15.9} cm{sup −2}). We discuss a potential scenario for the detection of CO without H{sub 2}, which deserves further investigation. The low abundance ratio measured around HK Ori suggests significant depletion of CO in the circumstellar gas, which conforms with the handful of other recent CO abundance measurements in protoplanetary disks.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Polyethylene Glycol Crowder’s Effect on Enzyme Aggregation, Thermal Stability, and Residual Catalytic Activity

Protein stability and performance in various natural and artificial systems incorporating many other macromolecules for therapeutic, diagnostic, sensor, and biotechnological applications attract increasing interest with the expansion of these technologies. Here we address the catalytic activity of lysozyme protein (LYZ) in the presence of a polyethylene glycol (PEG) crowder in a broad range of concentrations and temperatures in aqueous solutions of two different molecular mass PEG samples (Mw = 3350 and 10000 g/mol). The phase behavior of PEG–protein solutions is examined by using dynamic light scattering (DLS) and small-angle X-ray scattering (SAXS), while the enzyme denaturing is monitored by using an activity assay (AS) and circular dichroism (CD) spectroscopy. Molecular dynamic (MD) simulations are used to illustrate the effect of PEG concentration on protein stability at high temperatures. The results demonstrate that LYZ residual activity after 1 h incubation at 80 °C is improved from 15% up to 55% with the addition of PEG. The improvement is attributed to two underlying mechanisms. (i) Primarily, the stabilizing effect is due to the suppression of the enzyme aggregation because of the stronger PEG–protein interactions caused by the increased hydrophobicity of PEG and lysozyme at elevated temperatures. (ii) The MD simulations showed that the addition of PEG to some degree stabilizes the secondary structures of the enzyme by delaying unfolding at elevated temperatures. The more pronounced effect is observed with an increase in PEG concentration. This trend is consistent with CD and AS experimental results, where the thermal stability is strengthened with increasing of PEG concentration and molecular mass. Overall, the results show that the highest stabilizing effect is approached at the critical overlap concentration of PEG.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

A Molecular Gas Ring Hidden in the Sombrero Galaxy

We present Herschel, ALMA, and MUSE observations of the molecular ring of Messier 104, also known as the Sombrero galaxy. These previously unpublished archival data shed new light on the content of the interstellar medium of M104. In particular, molecular hydrogen measured by CO emission and dust measured by far-infrared light are uniformly distributed along the ring. The ionized gas revealed by Hα and [C II] emission is distributed in knots along the ring. Despite being classified as an SAa galaxy, M104 displays features typical of early-type galaxies. We therefore compared its [C II] and dust emission to a sample of early-type galaxies observed with Herschel and SOFIA. The [C II]/FIR ratio of M104 is much lower than that of typical star-forming galaxies and is instead much more similar to that of early-type galaxies. By classifying regions using optical emission line diagnostics, we also find that regions classified as H II lie closer to star-forming galaxies in the [C II]/FIR diagram than those classified as low-ionization emission regions. The good match between [C II] and Hα emission, in conjunction with the lack of correlation between CO emission and star formation, suggests that there is very limited active star formation along the ring and that most of the [C II] emission is from ionized and neutral atomic gas rather than molecular gas. From the total intensity of the CO line, we estimate a molecular hydrogen mass of 0.9 × 10 9 M ⊙ , a value intermediate between those of early-type galaxies and the content of the molecular ring of our galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

ACBD3 is up-regulated in gastric cancer and promotes cell cycle G1-to-S transition in an AKT-dependent manner

Highlights: • ACBD3 is up-regulated in gastric cancer. • ACBD3 regulates cell-cycle of gastric cancer in an AKT-dependent manner. It has been reported that ACBD3 is closely related to the malignant process of cells, but its role in gastric cancer has not been elucidated. This study aims to investigate the expression and function of ACBD3 in human gastric cancer. The Cancer Genome Atlas (TCGA) database were selected to analyze mRNA levels of ACBD3 in gastric cancer tissues and normal gastric epithelial tissues. qPCR and Western blot were conducted to detect the expression of ACBD3 in two normal gastric epithelial cell lines and five gastric cancer cell lines which were cultured in our laboratory. To exclude differences in individual background between different patients, we further detected the expression of ACBD3 in 8 pairs of malignant/non-malignant clinical gastric tissues. Through the establishment of stable cells, in vitro cell experiments and in vivo xenotransplantation models in mice, the role of ACBD3 in the proliferation of gastric cancer cells has been further explored. AKT inhibitors were used to deeply explore the molecular regulation mechanism of ACBD3. The results showed that the elevated ACBD3 in gastric cancer tissue were positively correlated with the clinical grade and prognosis of gastric cancer. In terms of molecular function, we found that ACBD3 can enhance the production and growth of gastric cancer cells. At the same time, the activation of AKT kinase played an important role in ACBD3's promotion of G1-to-S transition. The experiments generally indicate that ACBD3 is expected to become a potential diagnostic molecule or therapeutic target for gastric cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Understanding the structure and composition of recalcitrant oligosaccharides in hydrolysate using high-throughput biotin-based glycome profiling and mass spectrometry

Novel Immunological and Mass Spectrometry Methods for Comprehensive Analysis of Recalcitrant Oligosaccharides in AFEX Pretreated Corn Stover. Lignocellulosic biomass is a sustainable alternative to fossil fuel and is extensively used for developing bio-based technologies to produce products such as food, feed, fuel, and chemicals. The key to these technologies is to develop cost competitive processes to convert complex carbohydrates present in plant cell wall to simple sugars such as glucose, xylose, and arabinose. Since lignocellulosic biomass is highly recalcitrant, it must undergo a combination of thermochemical treatment such as Ammonia Fiber Expansion (AFEX), dilute acid (DA), Ionic Liquid (IL) and biological treatment such as enzyme hydrolysis and microbial fermentation to produce desired products. However, when using commercial fungal enzymes during hydrolysis, only 75–85% of the soluble sugars generated are monomeric sugars, while the remaining 15–25% are soluble recalcitrant oligosaccharides that cannot be easily utilized by microorganisms. Previously, we successfully separated and purified the soluble recalcitrant oligosaccharides using a combination of charcoal and celite-based separation followed by size exclusion chromatography and studies their inhibitory properties on enzymes. We discovered that the oligosaccharides with higher degree of polymerization (DP) containing methylated uronic acid substitutions were more recalcitrant towards commercial enzyme mixtures than lower DP and neutral oligosaccharides. Here, we report the use of several complementary techniques that include glycome profiling using plant biomass glycan specific monoclonal antibodies (mAbs) to characterize sugar linkages in plant cell walls and enzymatic hydrolysate, matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS) using structurally-informative diagnostic peaks offered by negative ion post-secondary decay spectra, gas chromatography followed by mass spectrometry (GC–MS) to characterize oligosaccharide sugar linkages with and without derivatization. Since oligosaccharides (DP 4–20) are small, it is challenging to mobilize these molecules for mAbs binding and characterization. To overcome this problem, we have applied a new biotin-coupling based oligosaccharide immobilization method that successfully tagged most of the low DP soluble oligosaccharides on to a micro-plate surface followed by specific linkage analysis using mAbs in a high-throughput system. This new approach will help develop more advanced versions of future high throughput glycome profiling methods that can be used to separate and characterize oligosaccharides present in biomarkers for diagnostic applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Experimental measurements of soot formation in fuel-rich homogeneous mixtures using an optical rapid compression machine

Advanced combustion strategies are necessary for the use of more environmentally sustainable fuels than traditional diesel. Alcohol fuels and alcohol/gasoline blends are of particular interest as they are readily available in the marketplace. Heavy-duty engines typically use compression ignited, conventional diesel mixing controlled combustion. Mixing controlled combustion features a non-premixed diffusion flame with a wide range of local equivalence ratios, leading to potentially high rates of soot formation. This work studies the sooting behavior of iso-octane and ethanol as a function of equivalence ratio. Measurements are carried out in a rapid compression machine (RCM) and are reported for pre-ignition conditions of 10–30 bar and temperatures of 650–800 K. Theoretical equilibrium and bulk gas temperatures are calculated for both fuels. These data are used to identify the critical equivalence ratio, the lowest equivalence ratio where soot is detected with a single-pass laser extinction diagnostic. The critical equivalence ratio for iso-octane varies between 1.82 and 1.77 for compressed pressures of 10 and 20 bar, respectively. Ethanol, sometimes considered sootless, had a critical equivalence ratio between 2.37 and 2.12 for compressed pressures of 20 and 30 bar, respectively. When characterizing soot formation by oxygenated equivalence ratio, the critical equivalence ratios for ethanol approach those of iso-octane. This suggests the oxygenated nature of alcohol fuels reduces sooting tendency, but other factors such as fuel molecular structure and morphology may play a role. It was observed that ethanol will form soot at equivalence ratios only slightly higher than iso-octane, which could have implications in mixing controlled combustion. It was seen for both fuels that soot formation is pressure sensitive, with the critical equivalence ratio being inversely proportional to compressed pressure and the rate of soot formation. Future work will investigate the sooting behavior of gasoline/ethanol blends.

02 PETROLEUM↗
DOE OSTI2025

Spectroscopic Identification of Carbamate Formation and Synergistic Binding in Amide–CO$^{–}_{2}$ Complexes

Carbamate formation is an elementary step that governs nitrogencentered nucleophilic CO 2 capture across diverse environments, yet a direct, structure-specific understanding has been lacking. Here, we report the first gas-phase characterization of closed-shell carbamate formation by deprotonated amides, using cryogenic ion trap vibrational spectroscopy combined with quantum chemical calculations. Reactions of deprotonated benzamide and isophthalamide anions with CO 2 form carbamate species, with the amide nitrogen serving as the nucleophilic site. Diagnostic, strongly red-shifted antisymmetric CO 2 stretching vibrational bands, supported by a bonding analysis, establish chemisorption with substantial charge transfer. In the multiamide system, an intramolecular N−H···O hydrogen bond provides synergistic stabilization, correlating with larger red shifts and more exothermic binding. These structure-assigned benchmarks provide molecular-level insights into the binding motif, charge redistribution, and hydrogen bond-mediated stabilization of amide carbamates, which aid the characterization of carbamate formation in condensed phase.

Amides↗
DOE OSTI2020

Low Energy Implantation into Transition-Metal Dichalcogenide Monolayers to Form Janus Structures

Atomically-thin two-dimensional (2D) materials face significant energy barriers for synthesis and processing into functional metastable phases such as Janus structures. Here, in this study, the controllable implantation of hyperthermal species from pulsed laser deposition (PLD) plasmas is introduced as a top-down method to compositionally engineer 2D monolayers. The kinetic energies of Se clusters impinging on suspended monolayer WS 2 crystals were controlled in the <10 eV/atom range with in situ plasma diagnostics to determine the thresholds for selective top layer replacement of sulfur by selenium for the formation of high quality WSSe Janus monolayers at low (300 °C) temperatures, and bottom layer replacement for complete conversion to WSe 2 . Atomic-resolution electron microscopy and spectroscopy in tilted geometry confirm the WSSe Janus monolayer. Molecular dynamics simulations reveal that Se clusters implant to form disordered metastable alloy regions, which then recrystallize to form highly ordered structures, opening the door for low-energy implantation by PLD as a novel method to explore the synthesis of 2D Janus layers and alloys of variable composition.

2D materials↗
DOE OSTI2020

Global Issues of Radiopharmaceutical Access and Availability: a Nuclear Medicine Global Initiative Project

The Nuclear Medicine Global Initiative (NMGI) was formed in 2012 by 13 international organizations to promote human health by advancing the field of nuclear medicine and molecular imaging by supporting the practice and application of nuclear medicine. The first project focused on standardization of administered activities in pediatric nuclear medicine and resulted in two manuscripts. For its second project the NMGI chose to explore issues impacting on access and availability of radiopharmaceuticals around the world. Methods: Information was obtained by survey responses from 35 countries on available radioisotopes, radiopharmaceuticals and kits for diagnostic and therapeutic use. Issues impacting on access and availability of radiopharmaceuticals in individual countries were also identified. Results: Detailed information on radiopharmaceuticals utilized in each country, and sources of supply, was evaluated. Responses highlighted problems in access particularly due to the reliance on a sole provider, regulatory issues and reimbursement, as well as issues of facilities and workforce particularly in low- and middle-income countries. Conclusion: Finally, strategies to address access and availability of radiopharmaceuticals are outlined, to enable timely and equitable patient access to nuclear medicine procedures worldwide. In the face of disruptions to global supply chains by the COVID-19 outbreak, renewed focus on ensuring reliable supply of radiopharmaceuticals is a major priority for nuclear medicine practice globally.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

A compact stem-loop DNA aptamer targets a uracil-binding pocket in the SARS-CoV-2 nucleocapsid RNA-binding domain

Abstract SARS-CoV-2 nucleocapsid (N) protein is a structural component of the virus with essential roles in the replication and packaging of the viral RNA genome. The N protein is also an important target of COVID-19 antigen tests and a promising vaccine candidate along with the spike protein. Here, we report a compact stem-loop DNA aptamer that binds tightly to the N-terminal RNA-binding domain of SARS-CoV-2 N protein. Crystallographic analysis shows that a hexanucleotide DNA motif (5′-TCGGAT-3′) of the aptamer fits into a positively charged concave surface of N-NTD and engages essential RNA-binding residues including Tyr109, which mediates a sequence-specific interaction in a uracil-binding pocket. Avid binding of the DNA aptamer allows isolation and sensitive detection of full-length N protein from crude cell lysates, demonstrating its selectivity and utility in biochemical applications. We further designed a chemically modified DNA aptamer and used it as a probe to examine the interaction of N-NTD with various RNA motifs, which revealed a strong preference for uridine-rich sequences. Our studies provide a high-affinity chemical probe for the SARS-CoV-2 N protein RNA-binding domain, which may be useful for diagnostic applications and investigating novel antiviral agents.

Biochemistry & Molecular Biology↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-zirconium (Nb-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

Amoeba species colonizing the gills of rainbow trout ( Oncorhynchus mykiss ) in Swiss aquaculture

Nodular gill disease (NGD) is an infectious condition characterized by proliferative gill lesions leading to respiratory problems, oxygen deficiency and mortality in fish. Globally, NGD primarily impacts freshwater salmonids in intensive aquaculture systems. In recent years, numerous outbreaks of severe gill disease have affected more than half of the larger rainbow trout ( Oncorhynchus mykiss ) farms in Switzerland, mainly during spring and early summer. Mortality has reached up to 50% in cases where no treatment was administered. Freshwater amoeba are the presumed aetiologic agent of NGD. The gross gill score (GS) categorising severity of gill pathology is a valuable first-line diagnostic tool aiding fish farmers in identifying and quantifying amoebic gill disease (AGD) in farmed marine salmonids. In this study, the GS was adapted to the NGD outbreak in farmed trout in Switzerland. In addition to scoring disease severity, gill swabs from NGD-affected rainbow trout were sampled and amoeba were cultured from these swabs. Morphologic and molecular methods identified six amoeba strains: Cochliopodium sp., Naegleria sp., Vannella sp., Ripella sp., Saccamoeba sp. and Mycamoeba sp. However, the importance of the different amoeba species for the onset and progression of NGD still has to be evaluated. This paper presents the first description of NGD with associated amoeba infection in farmed rainbow trout in Switzerland.

59 BASIC BIOLOGICAL SCIENCES↗