Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrostatic”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

From ionic clusters dynamics to network constraints in ionic polymer solutions

Physical networks formed by ionizable polymers with ionic clusters as crosslinks are controlled by coupled dynamics that transcend from ionic clusters through chain motion to macroscopic response. Here, the coupled dynamics, across length scales, from the ionic clusters to the networks in toluene swollen polystyrene sulfonate networks, were directly correlated, as the electrostatic environment of the physical crosslinks was altered. The multiscale insight is attained by coupling neutron spin echo measurements with molecular dynamics simulations, carried out to times typical of relaxation of polymers in solutions. The experimental dynamic structure factor is in outstanding agreement with the one calculated from computer simulations, as the networks are perturbed by elevating the temperature and changing the electrostatic environment. In toluene, the long-lived clusters remain stable over hundreds of ns across a broad temperature range, while the polymer network remains dynamic. In conclusion, though the size of the clusters changes as the dielectric constant of the solvent is modified through the addition of ethanol, they remain stable but morph, enhancing the polymer chain dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

First Measurement of Drift-Alfvén Wave Polarization in Magnetically Confined Fusion Plasmas

Polarization of drift-Alfvén waves, defined as the ratio of electrostatic to electromagnetic fluctuations, has remained unmeasurable in fusion plasmas for decades, despite its pivotal role in understanding wave dynamics and their impact on plasmas. We report the first measurements of drift-Alfvén wave polarization in a hot, magnetically-confined plasma. Here, the breakthrough is enabled by a novel methodology developed from gyrokinetic theory, utilizing fluctuations of electron temperature and density. Analysis of data from the DIII-D tokamak reveals that the waves above the geodesic acoustic mode frequency exhibit dominant electromagnetic polarization, whereas lower frequency waves show a mix of electromagnetic and electrostatic polarization, indicating a strong coupling between shear Alfvén waves and drift-acoustic waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Improved chemistry restraints for crystallographic refinement by integrating the Amber force field into Phenix

The refinement of biomolecular crystallographic models relies on geometric restraints to help to address the paucity of experimental data typical in these experiments. Limitations in these restraints can degrade the quality of the resulting atomic models. Here, an integration of the full all-atom Amber molecular-dynamics force field into Phenix crystallographic refinement is presented, which enables more complete modeling of biomolecular chemistry. The advantages of the force field include a carefully derived set of torsion-angle potentials, an extensive and flexible set of atom types, Lennard–Jones treatment of nonbonded interactions and a full treatment of crystalline electrostatics. The new combined method was tested against conventional geometry restraints for over 22 000 protein structures. Structures refined with the new method show substantially improved model quality. On average, Ramachandran and rotamer scores are somewhat better, clashscores and MolProbity scores are significantly improved, and the modeling of electrostatics leads to structures that exhibit more, and more correct, hydrogen bonds than those refined using traditional geometry restraints. In general it is found that model improvements are greatest at lower resolutions, prompting plans to add the Amber target function to real-space refinement for use in electron cryo-microscopy. This work opens the door to the future development of more advanced applications such as Amber -based ensemble refinement, quantum-mechanical representation of active sites and improved geometric restraints for simulated annealing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Influence of oxygen vacancies on the domain wall stability in BiFeO 3

Oxygen vacancy is the most common type of point defects in functional oxides, and it is known to have profound influence on their properties. This is particularly true for ferroelectric oxides since their interaction with ferroelectric polarization often dictates the ferroelectric responses. Here, we study the influence of the concentration of oxygen vacancies on the stability of ferroelectric domain walls (DWs) in BiFeO 3 , a material with a relatively narrow bandgap among all perovskite oxides, which enables strong interactions among electronic charge carriers, oxygen vacancies, and ferroelectric domains. It is found that the electronic charge carriers in the absence of oxygen vacancies have essentially no influence on the spatial polarization distribution of the DWs due to their low concentrations. Upon increasing the concentration of oxygen vacancies, charge-neutral DWs with an originally symmetric polarization distribution symmetric around the center of the wall can develop a strong asymmetry of the polarization field, which is mediated by the electrostatic interaction between polarization and electrons from the ionization of oxygen vacancies. Strongly charged head-to-head DWs that are unstable without oxygen vacancies can be energetically stabilized in the off-stoichiometric BiFeO 3–δ with δ ~ 0.02 where ionization of oxygen vacancies provides sufficient free electrons to compensate the bound charge at the wall. Our results delineate the electrostatic coupling of the ionic defects and the associated free electronic charge carriers with the bound charge in the vicinity of neutral and charged DWs in perovskite ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Demonstration of p-side down green light emitting diodes with high external quantum efficiencies

GaN/(In,Ga)N heterostructure based visible Light Emitting Diodes (LED) have enabled a wide range of solid-state lighting applications through excellent efficiency and power output in the shorter wavelengths (≤ 475nm) of violet/blue emission. However, the efficiency of emitters in the longer wavelength range (≥ 500nm) drops drastically due to the need to include higher Indium-content in the InGaN quantum wells. Large average polarization fields for high Indium-content quantum wells for conventional P-up structure, opposes the depletion field leading to large electrostatic barriers for both electrons and holes injection. LEDs fabricated along the N-polar direction with a p-up orientation or Ga-polar direction with p-down orientation lower such electrostatic barriers to carrier injection due to alignment of the polarization dipole field and depletion region field. This can therefore theoretically improve the electrical injection efficiency and reduce the forward voltage of operation. Such a Ga-polar p-down LED requires a bottom buried tunnel junction to avoid current spreading issues for a buried p-GaN layer. In this report, we demonstrate for the first time Ga-polar p-down green emitting LEDs using bottom tunnel junctions and having external quantum efficiencies comparable to those of equivalent p-up LEDs grown by Metal Organic Chemical Vapor Deposition (MOCVD).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantized critical supercurrent in SrTiO 3 -based quantum point contacts

Superconductivity in SrTiO 3 occurs at remarkably low carrier densities and therefore, unlike conventional superconductors, can be controlled by electrostatic gates. Here, we demonstrate nanoscale weak links connecting superconducting leads, all within a single material, SrTiO 3 . Ionic liquid gating accumulates carriers in the leads, and local electrostatic gates are tuned to open the weak link. These devices behave as superconducting quantum point contacts with a quantized critical supercurrent. This is a milestone toward establishing SrTiO 3 as a single-material platform for mesoscopic superconducting transport experiments that also intrinsically contains the necessary ingredients to engineer topological superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Altered morphology and diffusivity of water confined in MXenes: Machine learning–accelerated computations combined with experiments

Nanoconfined water exhibits unique properties compared to bulk water due to limited quantities, frustrated hydrogen bonding, and surface interactions, which are fundamental for energy storage and transport applications. We integrate machine learning–accelerated ab initio molecular dynamics with x-ray diffraction (XRD) and inelastic neutron scattering (INS) to systematically analyze the thermodynamic and dynamic behavior of water confined between functionalized (-F, -O, and -OH) two-dimensional (2D) Ti 3 C 2 T x MXene layers. As water intercalates between layers, the interlayer spacing exhibits layer-dependent staging characteristics. The water polarization can be flipped by the count and morphology of intercalated molecules interacting with MXene surface groups, resulting in varying electrostatic potential profiles. On the basis of interfacial electrostatic potential, hydrogen bond lifetime, and molecular orientation, we establish a linear combination of exponential model describing water diffusivity. These computational insights align well with experimental x-ray and neutron measurements, suggesting strategies for tuning water morphology and transport by tailoring MXene surface chemistry and water content for electrochemical energy storage and nanofluidic applications.

Tang, Jiawei [Southeast Univ., Nanjing (China)] (O↗
DOE OSTI2024

Electrically driven proton transfer promotes Brønsted acid catalysis by orders of magnitude

Electric fields play a key role in enzymatic catalysis and can enhance reaction rates by 100,000-fold, but the same rate enhancements have yet to be achieved in thermochemical heterogeneous catalysis. In this work, we probe the influence of catalyst potential and interfacial electric fields on heterogeneous Brønsted acid catalysis. We observed that variations in applied potential of ~380 mV led to a 100,000-fold rate enhancement for 1-methylcyclopentanol dehydration, which was catalyzed by carbon-supported phosphotungstic acid. Mechanistic studies support a model in which the interfacial electrostatic potential drop drives quasi-equilibrated proton transfer to the adsorbed substrate prior to rate-limiting C–O bond cleavage. Large increases in rate with potential were also observed for the same reaction catalyzed by Ti/TiO y H x and for the Friedel Crafts acylation of anisole with acetic anhydride by carbon-supported phosphotungstic acid.

Science & Technology - Other Topics↗
DOE OSTI2020

Suppression of Penning discharges between the KATRIN spectrometers

Abstract The KArlsruhe TRItium Neutrino experiment (KATRIN) aims to determine the effective electron (anti)-neutrino mass with a sensitivity of 0.2eV/c $$^2$$ 2 by precisely measuring the endpoint region of the tritium $$\beta $$ β -decay spectrum. It uses a tandem of electrostatic spectrometers working as magnetic adiabatic collimation combined with an electrostatic (MAC-E) filters. In the space between the pre-spectrometer and the main spectrometer, creating a Penning trap is unavoidable when the superconducting magnet between the two spectrometers, biased at their respective nominal potentials, is energized. The electrons accumulated in this trap can lead to discharges, which create additional background electrons and endanger the spectrometer and detector section downstream. To counteract this problem, “electron catchers” were installed in the beamline inside the magnet bore between the two spectrometers. These catchers can be moved across the magnetic-flux tube and intercept on a sub-ms time scale the stored electrons along their magnetron motion paths. In this paper, we report on the design and the successful commissioning of the electron catchers and present results on their efficiency in reducing the experimental background.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Charging and ion ejection dynamics of large helium nanodroplets exposed to intense femtosecond soft xray pulses

Ion ejection from charged helium nanodroplets exposed to intense femtosecond soft X-ray pulses is studied by single-pulse ion time-of-flight (TOF) spectroscopy in coincidence with small-angle X-ray scattering. Scattering images encode the droplet size and absolute photon flux incident on each droplet, while ion TOF spectra are used to determine the maximum ion kinetic energy, Ekin, of He+ j fragments (j = 1-4). Measurements span HeN droplet sizes between N similar to 10(7) and similar to 10(10) (radii R-0 = 78-578 nm), and droplet charges between similar to 9x10(-5) and similar to 3x10(-3) e/atom. Conditions encompass a wide range of ionization and expansion regimes, from departure of all photoelectrons from the droplet, leading to pure Coulomb explosion, to substantial electron trapping by the electrostatic potential of the charged droplet, indicating the onset of hydrodynamic expansion. The unique combination of absolute Xray intensities, droplet sizes, and ion Ekin on an event-by-event basis reveals a detailed picture of the correlations between the ionization conditions and the ejection dynamics of the ionic fragments. The maximum E-kin of He+ is found to be governed by Coulomb repulsion from unscreened cations across all expansion regimes. The impact of ion-atom interactions resulting from the relatively low charge densities is increasingly relevant with less electron trapping. The findings are consistent with the emergence of a charged spherical shell around a quasineutral plasma core as the degree of ionization increases. The results demonstrate a complex relationship between measured ion E-kin and droplet ionization conditions that can only be disentangled through the use of coincident single-pulse TOF and scattering data.

Helium↗
DOE OSTI2021

Electromechanics of lipid-modulated gating of potassium channels

Experimental studies reveal that the anionic lipid phosphatidic acid (POPA), non-phospholipid cholesterol, and cationic lipid DOTAP inhibit the gating of voltage-sensitive potassium (Kv) channels. In this work, we develop a continuum electromechanical model to investigate the interaction of these lipids with the ion channel. Our model suggests that: (i) POPA lipids may restrict the vertical motion of the voltage-sensor domain through direct electrostatic interactions; (ii) cholesterol may oppose the radial motion of the pore domain of the channel by increasing the mechanical rigidity of the membrane; and (iii) DOTAP can reduce the effect of electrostatic forces by regulating the dielectric constant at the channel–lipid interface. The electromechanical model predictions for the three lipid types match well with the experimental observations and provide mechanistic insights into lipid-dependent gating of Kv channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 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 OSTI2024

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Redesign of the Jefferson Lab -300 kV DC Photo-Gun for High Bunch Charge Operations

Production of high bunch charge beams for the Electron-Ion Collider (EIC) is a challenging task. High bunch charge (a few nC) electron beam studies at Jefferson Lab using an inverted insulator DC high voltage photo-gun showed evidence of space charge limitations starting at 0.3 nC, limiting the maximum delivered bunch charge to 0.7 nC for beam at -225 kV, 75 ps (FWHM) pulse width, and 1.64 mm (rms) laser spot size. The low extracted charge is due to the modest longitudinal electric field (Ez) at the photocathode leading to beam loss at the anode and downstream beam pipe. To reach the few nC high bunch charge goal, and to correct the beam deflection exerted by the non-symmetric nature of the inverted insulator photo-gun the existing photo-gun was modified. This contribution discusses the electrostatic design of the modified photo-gun obtained using CST Studio Suite’s electromagnetic field solver. Beam dynamics simulations performed using General Particle Tracer (GPT) with the resulting electrostatic field map obtained from the modified electrodes confirmed the validity of the new design.

Wijethunga, S.A.K.↗
DOE OSTI2021

Nanoporous weakly coordinating anionic frameworks (Final Technical Report)

Under this award, we synthesized and characterized new weakly coordinating, anionic covalent organic frameworks. The motivation was to expand the class of covalent organic frameworks though the introduction of the first weakly coordinating anionic ones. The significance of the weakly coordinating anionic character of these frameworks is associated with the character of the counter (metal) cations in their pores. Those counter-cations, due to the inability to effectively coordinate to the pore walls, have a “naked”, and thus highly electrophilic, Lewis acidic character. Such properties are expected to translate to a number of interesting physical properties that can be potentially used for a range of energy-relevant applications. One of these applications is gas adsorption, separation, and storage. The electrophilic, Lewis acidic character of the counter-cations are expected to strongly bind to gas molecules entering the pores leading to higher gas sorption capacities. Other applications include heterogeneous catalysis, ion exchange, and solid state ion conduction. Within the project we explored a range of chemical blocks, and tested synthesis conditions to polymerize these building blocks to form the target materials. In particular, we focused on materials with group 15 elements like phosphorus in the oxidation state +5 and the coordination number 6, that have a formal -1 charge. The synthesis of such materials is challenging because its success depends on the absence of residual monomeric and oligomeric side-products in the pores, and the absence of significant framework interpenetration. Crystalline forms of these materials require sufficient reversibility of the polymerization reaction (which is difficult to achieve because group 15 elements like P tend to make strong covalent bonds with neighbored elements), as well as suitable polymerization rates that prevent the formation of kinetically controlled products. We had success with the preparation of anionic porous organic frameworks with anionic building blocks containing P(V)O6 octahedra synthesized through Yamamoto coupling reactions. These materials are non-crystalline because of the irreversible nature of the Yamamoto coupling reaction, but show porosity. The porosity and surface area can be significantly increased by co-reaction with electrostatically neutral co-reactants. Very high porosity and surface areas can be achieved when the materials are washed with hydrochloric acid. This method is commonly used to free products from Yamamoto coupling reactions from by-products. Indeed, this method removed by products from the pores of our materials, however it also led to the hydrolysis of the PO6 octahedra to form pending phosphonic acid groups (-PO3H)-groups in the frameworks (PA-POFs) which was actually not desired. Nonetheless, these materials showed high CO2 sorption capacities. In addition, they show extremely high (actually record) sorption capacities for the removal of a range of common micropollutants such as bisphenol A (BPA) and 4-nitrophenol in water. BPA is a wide-spread organic pollutant and a known endocrine disruptor. The maximum adsorption capacity of BPA at equilibrium is found to be as high as 3,366 mg g-1 by Langmuir adsorption model, which is more than 10 times greater than peer materials. The polymer also rapidly removes various other organic micropollutants with more than 90% removal efficiencies. In addition, the PA-POF material can be regenerated at least five times by mild washing using methanol without significant loss in removal efficiency. This is significant advantage over commonly used activated carbon materials (that are in Brita water faucet filters, for example). Comparison studies showed that only materials with the phosphonate groups show extreme sorption although the number of PA group is actually small according to solid state NMR and EDX data. More research would be needed to investigate why this is (a possibility is that the phosphonate groups have a multiplying effect in which a PA-induced initial adsorption leads to the adsorption of much more adsorbate). The extreme sorption capabilities can also be visually seen when colored adsorbates such as 4-nitrophenol are used, as the adsorption of the nitrobenzene leads to the decoloration of the solution. We have produced three short videos that illustrate the sorption capabilities compared to two peer materials (activated carbon from dismantled Brita water filters, and a electrostatically neutral, chemically related, porous organic polymer, respectively).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of 2D-3D heterointerface on remote epitaxial interaction through graphene

Remote epitaxy has drawn attention as it offers epitaxy of functional materials that can be released from the substrates with atomic precision. The realization of remote epitaxy leads to profound impact on production and heterointegration of flexible, transferrable, and stackable freestanding single-crystalline membranes, which can be extremely useful in manufacturing next generation electronic and photonic devices. In addition, the remote interaction of atoms and adatoms through two-dimensional (2D) materials in remote epitaxy allows investigation and utilization of electrical/chemical/physical coupling of bulk (3D) materials via 2D materials (3D-2D-3D coupling). Here, we unveil the respective roles and impacts of the substrate material, graphene, substrate-graphene interface, and epitaxial material for electrostatic coupling of these materials. The coupling governs cohesive ordering and can lead to single-crystal epitaxy in the overlying film. We show that simply coating a graphene layer on wafers does not guarantee successful implementation of remote epitaxy, since atomically precise control of the graphene-coated interface is required. Also, it provides key considerations for maximizing the remote electrostatic interaction between the substrate and adatoms. This was enabled by exploring various material systems and processing conditions, and we demonstrate that the rules of remote epitaxy vary significantly depending on the ionicity of material systems as well as the graphene-substrate interface and the epitaxy environment. The general rule of thumb discovered here enables expanding 3D material libraries that can be stacked in freestanding form.

Lu, Kuangye↗