Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELECTROSTATIC CHARGING”

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

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Superhydrophobic Phosphonium Modified Robust 3D Covalent Organic Framework for Preferential Trapping of Charge Dispersed Oxoanionic Pollutants

Abstract Charge dispersed and less hydrophilic anionic pollutants are often difficult to be preferentially captured by common cationic framework materials. A gamma radiation modification approach is employed to anchor a superhydrophobic phosphonium‐containing coating on a robust crystalline 3D covalent organic framework (COF). By regulating the surface hydrophobicity, the prepared 3DCOF‐ g ‐VBPPh 3 Cl is endowed with a strong affinity for charge dispersed and less hydrophilic oxoanions, such as MnO 4 − , TcO 4 − and ReO 4 − , to surmount the Hofmeister bias, which much favors remediation of oxoanionic pollutants in complex aqueous systems. Batch and column experiments with 3DCOF‐g‐VBPPh 3 Cl in both tap water and simulated groundwater are performed, and rapid sequestration is achieved with removal efficacy up to 99.995%, record‐high distribution coefficient of 1.0 × 10 8 mL g −1 , and desirable priority over competing anions such as Cl − , SO 4 2− , HCO 3 − , and NO 3 − , confirming promise for remediation of charge dispersed anionic pollutants. Theoretical calculations reveal a mechanism of preferential capture based on electrostatic and dispersion interaction between charge dispersed anions and 3DCOF‐g‐VBPPh 3 Cl that effectively competes with solvation of the aqueous anions.

Wang, Yue↗
DOE OSTI2022

Electrical Gating of the Charge-Density-Wave Phases in Two-Dimensional h -BN/1T-TaS 2 Devices

Here we report on the electrical gating of the chargedensity- wave phases and current in h-BN-capped three-terminal 1T-TaS 2 heterostructure devices. It is demonstrated that the application of a gate bias can shift the source–drain current– voltage hysteresis associated with the transition between the nearly commensurate and incommensurate charge-density-wave phases. The evolution of the hysteresis and the presence of abrupt spikes in the current while sweeping the gate voltage suggest that the effect is electrical rather than self-heating. We attribute the gating to an electric-field effect on the commensurate charge-density-wave domains in the atomic planes near the gate dielectric. The transition between the nearly commensurate and incommensurate charge-density-wave phases can be induced by both the source– drain current and the electrostatic gate. Since the charge-densitywave phases are persistent in 1T-TaS 2 at room temperature, one can envision memory applications of such devices when scaled down to the dimensions of individual commensurate domains and few-atomic plane thicknesses.

2D van der Waals materials↗
DOE OSTI2024

A Gaussian field approach to the solvation of spherical ions in electrolyte solutions

In this paper, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory.

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

Understanding attenuated solvent reorganization energies near electrode interfaces

In this work, we examine the role of image charge effects on the electrostatic potential fluctuations experienced by ionic species in the vicinity of an electrode surface. We combine simulation and theory to quantify these fluctuations and how they vary with distance from the electrode surface. We observe that the potential distribution narrows significantly for species within a few electrolyte screening lengths of the electrode. We attribute this narrowing to the effects of image charge fluctuations originating from the polarization response of the electrode. We show that the physical consequences of these image charge effects can be captured in the context of a simple analytical field theory with anti-symmetric boundary conditions. We contextualize these results by discussing their implications for rates of Marcus-like outer-sphere interfacial electron transfer.

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

Gating Single Molecules with Counterions

We report atomic-scale gating and visualization of local charge distribution within individual rare-earth-based molecular complexes on a metallic surface. The complexes are formed by a positively charged lanthanum ion coordinated to a (pcam)3 molecule and a negatively charged counterion trapped underneath via electrostatic interactions on a Au(111) surface. Local gating is performed by adding an additional negatively charged counterion to one side of the complex, which results in the redistribution of charges within the complex and a positive shift of the frontier orbitals. This is caused by the internal Stark effect induced by the added counterion. This effect is directly captured using tunneling spectroscopy and spectroscopic mapping at 5 K substrate temperature. The polarizability of the complex is corroborated by density functional theory and analytical calculations based on experimental findings. Furthermore, the influence of charge polarization on nearby complexes is investigated in a cluster purposely assembled using three complexes, which reveals maintaining the charge states as in single complexes. These findings will enable the design of robust charged rare-earth complexes to be tailored for potential solid-state applications.

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 OSTI2020

Friction in a strongly magnetized neutral plasma

A charged projectile generates an electrostatic wake as it moves through a background plasma. This wake normally gives rise to a drag force that acts antiparallel to the velocity vector of the projectile, and is commonly referred to as the stopping power. A recent theory predicted that an additional component of the friction force that is perpendicular to both the projectile velocity and the Lorentz force arises when a plasma is strongly magnetized. We extend the previous analysis of this transverse friction force, which was based on the one-component plasma model, to treat a two-component charge-neutral electron-ion plasma. The direction and magnitude of the transverse force is found to depend on the projectile speed relative to the electron and ion thermal velocities, and to change sign three times. For projectile speeds below a certain threshold, the transverse force is dominated by collisions with ions, while above this threshold, only electron collisions contribute. Here, the average trajectory of the projectile is significantly altered by the transverse force as it slows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Decay of Electrostatic Force of Dust Particles on Photovoltaic Modules

Photovoltaic (PV) energy yield loss due to solar module soiling has become increasingly important as solar module deployment is now at the hundreds of gigawatts scale and continues to grow rapidly. We have reported on direct measurements-using atomic force microscopy (AFM)-of strong electric-field-induced attraction and adhesion force (Fes) of dust particles onto solar panel that are 1 to 2 orders of magnitude stronger than the van der Waals and water capillary forces, corroborated by observing the increase in system voltage-induced soiling rate. Here, we report another characteristic of Fes on soiling-long lasting or slow decay after turning off the high voltage applied to solar panels. The Fes decay time varies in a wide time range of 1 to 10 hours, depending on two factors: 1) either/both the cell or/and particle were charged with high voltage before the voltages were turned off; and 2) how the cell was connected to the ground after the voltage was turned off-either connected through the power supply electronics, directly connected to the ground, or electrically floated. The Fes decay is understood in terms of 1) net electrical charge dissipations in both particle and cell, 2) thermal disordering of dipole polarization in the module glass dielectrics, and 3) charge redistribution by the electrostatic interaction of particle and module glass. This long-lasting Fes for hours can significantly affect the solar panel soiling after sunset, especially combined with water condensation.

14 SOLAR ENERGY↗
DOE OSTI2021

Scaling of silicon nanoparticle growth in low temperature flowing plasmas

Low temperature plasmas are an emerging method to synthesize high quality nanoparticles (NPs). An established and successful technique to produce NPs is using a capacitively coupled plasma (CCP) in cylindrical geometry. Although a robust synthesis technique, optimizing or specifying NP properties using CCPs, is challenging. In this paper, results from a computational investigation for the growth of silicon NPs in flowing inductively coupled plasmas (ICPs) using Ar/SiH 4 gas mixtures of up to a few Torr are discussed. ICPs produce more locally constrained and quiescent plasma potentials. These positive plasma potentials produce an electrostatic trap for negatively charged NPs, which can significantly extend the residence time of NPs in the plasma, which in turn provides a controllable period for particle growth. The computational platforms used in this study consist of a two-dimensional plasma hydrodynamics model, a three-dimensional nanoparticle growth and trajectory tracking model, and a molecular dynamics simulation for deriving reactive sticking coefficients of silane radicals on Si NPs. Trends for the nanoparticle growth as a function of SiH 4 inlet fraction, gas residence time, energy deposition per particle, pressure, and reactor diameter are discussed. The general path for particle synthesis is the trapping of small NPs in the positive electrostatic potential, followed by entrainment in the gas flow upon reaching a critical particle size. Optimizing or controlling NP synthesis then depends on the spatial distribution of plasma potential, the density of growth species, and the relative time that particles spend in the electrostatic trap and flowing through higher densities of growth species upon leaving the trap.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Long-lasting strong electrostatic attraction and adhesion forces of dust particles on photovoltaic modules

Photovoltaic (PV) energy-yield loss due to solar module soiling has become increasingly important as solar module deployment is now at the hundreds of gigawatts scale and continues to grow rapidly. The electrostatic attraction and adhesion force (F es ) of dust particles by the high voltages from solar panels have been reported to be 1 to 2 orders of magnitude stronger than the van der Waals and water capillary forces, which has been corroborated by observing the increase in system voltage-induced soiling rate by setting up an outdoor test. Here, we report another characteristic of F es on soiling - long-lasting or slow decay after turning off the high voltage applied to solar panels. The F es decay time varies across a wide time range of 1 to 10 h, depending on two factors: 1) whether the cell, the particle, or both were charged with high voltage before the voltages were turned off; and 2) how the cell was connected to the ground after the voltage was turned off - either connected through the power-supply electronics, directly connected to the ground, or electrically floated. The Fes decay is understood in terms of 1) net electrical charge dissipations in both particle and cell, 2) thermal disordering of dipole polarization in the module glass dielectrics, and 3) charge redistribution by the electrostatic interaction of particle and module glass. This long-lasting F es for hours can affect the solar panel soiling after sunset, and it can have an even greater effect when combined with water condensation at night.

14 SOLAR ENERGY↗
DOE OSTI2023

DNA-caged nanoparticles via electrostatic self-assembly

DNA-modified nanoparticles enable DNA sensing and therapeutics in nanomedicine and are also crucial for nanoparticle self-assembly with DNA-based materials. However, methods to conjugate DNA to nanoparticle surfaces are limited, inefficient, and lack control. Inspired by DNA tile nanotechnology, we demonstrate a new approach to nanoparticle modification based on electrostatic attraction between negatively charged DNA tiles and positively charged nanoparticles. This approach does not disrupt nanoparticle surfaces and leverages the programmability of DNA nanotechnology to control DNA presentation. We demonstrated this approach using a variety of nanoparticles, including polymeric micelles, polystyrene beads, gold nanoparticles, and superparamagnetic iron oxide nanoparticles with sizes ranging from 5–20 nm in diameter. DNA cage formation was confirmed through transmission electron microscopy (TEM), neutralization of zeta potential, and a series of fluorescence experiments. DNA cages present “handle” sequences that can be used for reversible target attachment or self-assembly. Handle functionality was verified in solution, at the solid–liquid interface, and inside fixed cells, corresponding to applications in biosensing, DNA microarrays, and erasable immunocytochemistry. Finally, these experiments demonstrate the versatility of the electrostatic DNA caging approach and provide a new pathway to nanoparticle modification with DNA that will empower further applications of these materials in medicine and materials science.

60 APPLIED LIFE 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↗
DOE OSTI2022

Pulsed power to control growth of silicon nanoparticles in low temperature flowing plasmas

Low-temperature plasmas have seen increasing use for synthesizing high-quality, mono-disperse nanoparticles (NPs). Recent work has highlighted that an important process in NP growth in plasmas is particle trapping—small, negatively charged nanoparticles become trapped by the positive electrostatic potential in the plasma, even if only momentarily charged. In this article, results are discussed from a computational investigation into how pulsing the power applied to an inductively coupled plasma (ICP) reactor may be used for controlling the size of NPs synthesized in the plasma. The model system is an ICP at 1 Torr to grow silicon NPs from an Ar/SiH 4 gas mixture. This system was simulated using a two-dimensional plasma hydrodynamics model coupled to a three-dimensional kinetic NP growth and trajectory tracking model. The effects of pulse frequency and pulse duty cycle are discussed. We identified separate regimes of pulsing where particles become trapped for one pulsed cycle, a few cycles, and many cycles—each having noticeable effects on particle size distributions. For the same average power, pulsing can produce a stronger trapping potential for particles when compared to continuous wave power, potentially increasing particle mono-dispersity. Pulsing may also offer a larger degree of control over particle size for the same average power. Experimental confirmation of predicted trends is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Cellular Phosphate Sensing and Anion Binding by an Azacrown‐Calixpyrrole Hybrid

A hybrid receptor-sensor for anions originating from the merging of positively charged ammonium moieties for electrostatic attraction/stronger binding of azacrowns with directionality of calixpyrrole hydrogen bond donors for selectivity is investigated. As demonstrated this hybrid receptor-sensor shows a remarkable selectivity for orthophosphate even in the presence of other phosphates and anions found in cellular materials (K assoc H 2 PO 4 − >H 2 P 2 O 7 2− >AMP − ≫ADP 2− or ATP 3− over halides, nitrate, or hydrogen sulfate; all Na + salts in water) but also cellular polyphosphate or phospholipids. This selectivity is harnessed in a real-time monitoring of cell lysis by lysozyme, which releases orthophosphate and other phosphates and anions from the cells. This sensitive (LOD 0.4 μM) fluorescence-based microscale method compares favorably with the state-of-the-art techniques but can easily be practiced in a high-throughput screening (HTS) manner. The anion binding and selectivity in aqueous solutions were investigated by NMR and put in context with phosphate binding of the parent calix[4]pyrrole. The microscopic understanding of anion binding by the hybrid receptor was then obtained from a combination of density functional theory (DFT), classical molecular dynamics (MD) with explicit water solvation, and ab initio MD (AIMD) simulations. Correlating the NMR and fluorescence binding data with studies of solvation of the receptor, phosphate anion, and the resulting complex confirms the binding is largely driven by entropic component (TΔS) associated with receptor and anion desolvation.

Anions↗
DOE OSTI2022

Impact of hydrophobic pendant phenyl groups on transport and co-transport of methanol and acetate in PEGDA-SPMAK cation exchange membranes

CO 2 reduction cells are innovative devices that reduce CO 2 into valuable chemicals (i.e. methanol (MeOH) and acetate (OAc – )) on the cathode. A major challenge with such devices is to develop ion exchange membranes that allow ion-selective transport (i.e., protons for cation exchange membranes, CEM) and suppress the crossover of CO 2 reduction products. To design such membranes, it is important to understand multi-solute transport behavior of these solutes. Previously, our group reported acetate diffusivities in co-transport with MeOH increase in sulfonated CEMs, where we speculated charge screening of the electrostatic interactions by co-diffusing MeOH has an impact. Here, crosslinked membranes fabricated by photopolymerization of poly(ethylene glycol) diacrylate (PEGDA), 3-sulfopropyl methacrylate potassium (SPMAK, SO 3 – -containing ionomer), and a phenyl-containing comonomer either phenoxyethyl acrylate (PEA) or poly(ethylene glycol) phenyl ether acrylate (PEGPEA)) are investigated. We observe OAc – diffusivities to both (1) PEGDA-SPMAK and (2) PEGDA-PEGPEA increase in co-diffusion with MeOH and those to (3) PEGDA-SPMAK/PEGPEA decrease. To rationalize this emergent co-transport behavior, we speculate (1) electrostatic interactions are interfered with by co-diffusing MeOH, (2) chain mobility (segmental dynamics) increases in the presence of MeOH and (3) for films with both hydrophobic (PEGPEA) and hydrophilic (SPMAK) comonomers chain mobility (segmental dynamics) decrease due to interactions between the comonomer sidechains depressing overall solute diffusivities. Finally, while further investigations are needed, this work contributes to improving our fundamental understanding of the relationships between polymer film chemistry, solute chemistry, and emergent cotransport behavior observed and described in this work (and others) towards enabling the design of improved membrane materials.

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

Molecular Cage Reports on Its Contents: Spectroscopic Signatures of Cryo-Cooled K + - and Ba 2+ -Benzocryptand Complexes

UV photofragment spectroscopy and IR–UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K + or Ba 2+ imbedded inside it (labeled as K + -BzCrypt, Ba 2+ -BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH 2 -CH 2 -O-CH 2 -CH 2 -O-CH 2 -CH 2 -N chains provides a UV chromophore whose S 0 –S 1 transition we probe. K + -BzCrypt and Ba 2+ -BzCrypt have their S 0 –S 1 origin transitions at 35,925 and 36,446 cm –1 , respectively, blue-shifted by 174 and 695 cm –1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR–UV dip spectroscopy. The alkyl CH stretch region (2800–3000 cm –1 ) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K + -BzCrypt and Ba 2+ -BzCrypt, identifying two conformations (A and B) that lie below all others in energy. Here, we extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH 2 groups adjacent to these atoms. The localized CH 2 scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH 2 wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm –1 . As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH 2 symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗