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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2022

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

Lu, KUANGYE↗

DOE OSTI2023

Transformer Analysis with ANSYS Maxwell and COMSOL

This is an extension of work described by Rodriguez et al . (2021). It continues analyses of a generic transformer design by Wes Greenwood. In this report, we summarize that work and add comparable results using the ANSYS Maxwell software (henceforward, “ANSYS”), and with COMSOL . We found the ANSYS and COMSOL calculations of inductance agreed well with previous results for simplified coils in air, and with a ferromagnetic core. We then describe the ANSYS and COMSOL approach and show results for a full transformer model based on magnetic field analyses. Finally, we present electrostatic analyses of E field enhancement, once again resolving individual wires. The purpose is to assess the electrostatic fields in order to locate where electric breakdown is likely to originate. We found the maximum enhancement between the secondary and either the primary or the tertiary at the end of the windings with a large potential difference.

42 ENGINEERING↗

DOE OSTI2025

dynamics of organic-mineral interactions at the metal oxide-solution interface as studied via binding energetics (Final report)

This project focused on addressing longstanding fundamental and experimental uncertainties on how dissolved organic substances (DOS) interact with metal oxide surface under environmentally relevant conditions. By leveraging a custom-built real-time, in-tandem flow adsorption microcalorimetry-UV-Vis/fluorescence spectroscopy platform, we characterized the binding energetics, kinetics and mechanistic pathways driving DOS-metal oxide interactions at temporal resolution on the order of 1-5 seconds. We studied a diverse suite of model organic compounds/substances – including monocarboxylates (e.g. acetate and benzoate), di-carboxylates (oxalate and succinate), amino acids, amino-based nanparticles and natural organic matter – interacting at the mineral-water interface of structurally- and/or chemically distinct metal oxides (including SiO2, boehmite, ferrihydrite, and γ-Al2O3). Our results indicated that DOS-metal oxide interactions are governed by multi-step reaction pathways, often switching between distinct, resolvable enthalpy- and entropy-driven non-electrostatic or electrostatic configurations. To quantify these interactions, we developed and implemented an analytical workflow that integrates peak deconvolution and Monte-Carlo based error propagation to determine site-specific thermodynamic and kinetic parameters for individual binding/debinding events. In addition to resolving apparent first-order rate constants of each event, we were able to quantify associated apparent equilibrium constants as well as free energy, enthalpy and entropy contribution to the activation and subsequent progression of the binding/debinding process across compounds, compound class and metal oxide surfaces. The kinetic-thermodynamic data produced in this study captured how the interplay between oxide surface reactivity and DOS molecular structure jointly drives binding-debinding dynamics. Notably, that at pH below PZC of the oxide surface, neutral species were heavily involved in monocarboxylate binding, while anionic species drove dicarboxylate binding. Also, that among amino acids 1) positional isomers show distinctive binding characteristics to each other while enantiomers show no significant differences in binding characteristics, 2) molecules that bind via outer-sphere complexation show a larger entropic shift between binding and debinding with no impact on oxide surface while 3) inner-sphere interactions increased anion exchange capacity of the oxide surface. The new insights and data from this work has great potential for improving predictive modeling of carbon dynamics and specifically organic-mineral interactions in environmental and industrial systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Tau and Membranes: Interactions That Promote Folding and Condensation

Tau misfolding and assembly is linked to a number of neurodegenerative diseases collectively described as tauopathies, including Alzheimer’s disease (AD) and Parkinson’s disease. Anionic cellular membranes, such as the cytosolic leaflet of the plasma membrane, are sites that concentrate and neutralize tau, primarily due to electrostatic interactions with tau’s microtubule binding repeat domain (RD). In addition to electrostatic interactions with lipids, tau also has interactions with membrane proteins, which are important for tau’s cellular functions. Tau also interacts with lipid tails to facilitate direct translocation across the membrane and can form stable protein-lipid complexes involved in cell-to-cell transport. Concentrated tau monomers at the membrane surface can form reversible condensates, change secondary structures, and induce oligomers, which may eventually undergo irreversible crosslinking and fibril formation. These β-sheet rich tau structures are capable of disrupting membrane organization and are toxic in cell-based assays. Given the evidence for relevant membrane-based tau assembly, we review the emerging hypothesis that polyanionic membranes may serve as a site for phase-separated tau condensation. Membrane-mediated phase separation may have important implications for regulating tau folding/misfolding, and may be a powerful mechanism to spatially direct tau for native membrane-mediated functions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Acid-Base Equilibrium and Dielectric Environment Regulate Charge in Supramolecular Nanofibers

Peptide amphiphiles are a class of molecules that can self-assemble into a variety of supramolecular structures, including high-aspect-ratio nanofibers. It is challenging to model and predict the charges in these supramolecular nanofibers because the ionization state of the peptides are not fixed but liable to change due to the acid-base equilibrium that is coupled to the structural organization of the peptide amphiphile molecules. Here, we have developed a theoretical model to describe and predict the amount of charge found on self-assembled peptide amphiphiles as a function of pH and ion concentration. In particular, we computed the amount of charge of peptide amphiphiles nanofibers with the sequence C 16 − V 2 A 2 E 2 . In our theoretical formulation, we consider charge regulation of the carboxylic acid groups, which involves the acid-base chemical equilibrium of the glutamic acid residues and the possibility of ion condensation. The charge regulation is coupled with the local dielectric environment by allowing for a varying dielectric constant that also includes a position-dependent electrostatic solvation energy for the charged species. We find that the charges on the glutamic acid residues of the peptide amphiphile nanofiber are much lower than the same functional group in aqueous solution. There is a strong coupling between the charging via the acid-base equilibrium and the local dielectric environment. Our model predicts a much lower degree of deprotonation for a position-dependent relative dielectric constant compared to a constant dielectric background. Furthermore, the shape and size of the electrostatic potential as well as the counterion distribution are quantitatively and qualitatively different. These results indicate that an accurate model of peptide amphiphile self-assembly must take into account charge regulation of acidic groups through acid–base equilibria and ion condensation, as well as coupling to the local dielectric environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗