Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electrostatic Potential”

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 523 records · Page 29

Quantitative Determination of Metal Ion Adsorption on Cellulose Nanocrystals Surfaces

Nanocellulose is a bio-based material that holds significant potential in the field of water purification. Of particular interest is their potential use as a key sorbent material for the removal of metal ions from solution. However, the structure of metal ions adsorbed onto cellulose surfaces is not well understood. The focus of this work is to determine quantitatively the three-dimensional distribution of metal ions of different valencies surrounding negatively charged carboxylate functionalized cellulose nanocrystals (CNCs) using anomalous small-angle X-ray scattering (ASAXS). These distributions can affect the water and ionic permeability in these materials. The data show that increasing the carboxylate density on the surface of the CNCs from 740 to 1100 mmol/kg changed the nature of the structure of the adsorbed ions from a monolayer into a multilayer structure. The monolayer was modeled as a Stern layer around the CNC nanoparticles, whereas the multilayer structure was modeled as a diffuse layer on top of the Stern layer around the nanoparticles. Within the Stern layer, the maximum ion density increases from 1680 to 4350 mmol of Rb + /(kg of CNC) with the increase in the carboxylate density on the surface of the nanoparticles. Additionally, the data show that CNCs can leverage multiple mechanisms, such as electrostatic attraction and the chaotropic effect, to adsorb ions of different valencies. By understanding the spatial organization of the adsorbed metal ions, the design of cellulose-based sorbents can be further optimized to improve the uptake capacity and selectivity in separation applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon nanospikes have improved sensitivity and antifouling properties for adenosine, hydrogen peroxide, and histamine

Carbon nanospikes (CNSs) are a new nanomaterial that has enhanced surface roughness and surface oxide concentration, increasing the sensitivity for dopamine detection. However, CNS-modified electrodes (CNSMEs) have not been characterized for other neurochemicals, particularly those with higher oxidation potentials. The purpose of this study was to evaluate CNSMEs for the detection of adenosine, hydrogen peroxide (H 2 O 2 ), and histamine. The sensitivity increased with CNSs, and signals at CNSMEs were about 3.3 times higher than CFMEs. Normalizing for surface area differences using background currents, CNSMEs show an increased signal of 4.8 times for adenosine, 1.5 times for H 2 O 2 , and 2 times for histamine. CNSMEs promoted the formation of secondary products for adenosine and histamine, which enables differentiation from other analytes with similar oxidation potentials. Furthermore, CNSs also selectively enhance the sensitivity for adenosine and histamine compared to H 2 O 2 . A scan rate test reveals that adenosine is more adsorption-controlled at CNS electrodes than CFMEs. CNSMEs are antifouling for histamine, with less fouling because the polymers formed after histamine electrooxidation do not adsorb due to an elevated number of edge planes. CNSMEs were useful for detecting each analyte applied in brain slices. Because of the hydrophilic surface compared to CFMEs, CNSMEs also have reduced biofouling when used in tissue. Therefore, CNSMEs are useful for tissue measurements of adenosine, hydrogen peroxide, and histamine with high selectivity and low fouling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Subtask 2.6 – Optimization of Aerosol Mitigation Technology for Postcombustion CO 2 Capture

Growing concerns over the impact of CO 2 emissions from combustion sources on global climate change have prompted numerous research and development projects aimed at developing cost-effective technologies for CO 2 capture. One family of technologies being demonstrated at pilot and full scale globally is postcombustion CO 2 capture (PCCC) systems that employ amine-based solvents. The captured CO 2 can be compressed and permanently stored underground or used for enhanced oil recovery (EOR). The proximity of North Dakota’s lignite-fired fleet of power plants to potential CO 2 storage options creates a unique atmosphere for PCCC within the state. However, the aerosols present in lignite flue gas present a challenge for large-scale PCCC at North Dakota power plants. Aerosols can negatively impact the long-term performance of amine-based solvents for CO 2 capture. Amine-based solvents are volatile, and aerosols provide nucleation sites where amine vapors can condense. Because aerosols cannot be easily captured at the column outlet using conventional technologies, the amine-laden aerosols escape the system and lead to accelerated solvent losses. Moreover, aerosol components can chemically react with amines to form degradation products that can permanently deactivate the amine, cause fouling, and lead to hazardous emissions. Many of the elements that have been shown to catalyze solvent degradation are present in lignite coals and can exacerbate solvent replacement economics. Understanding this issue is critical to the implementation of solvent-based CO 2 capture systems as applied to lignite-fired generation systems. The Energy & Environmental Research Center (EERC) designed and carried out this project to optimize aerosol mitigation technology for PCCC at a lignite-fired power plant. To meet the goal of this project, the following objectives were identified: Determine the effectiveness of a wet electrostatic precipitator (WESP) on collection of aerosols at a low-rank coal-fired power station. Determine the impact of aerosols on the efficiency and degradation products of amine-based carbon capture systems fired with low-rank fuels. Work was conducted at Minnkota Power Cooperative’s (MPC’s) Milton R. Young (MRY) Station Unit 2 using a slipstream of flue gas from the outlet of the plant’s flue gas desulfurization (FGD) unit. To gather initial data for sizing and specifying a WESP for this system, a temporary pilot-scale WESP was rented and installed on-site. Several different conditions were tested to examine the impact of flow rate, voltage, and current on WESP performance. The WESP was effective at removing large particulate (>200 nm) but caused an increase in fine particulate (<75 nm). Fine particulate material at the inlet and outlet of the WESP was collected, analyzed, and showed that crystalline sulfates carried over from the plant’s FGD unit were being converted to fine aerosols and SO 2 was being converted to SO 3 through the WESP. Additionally, the high moisture content of the flue gas stream at this sample location also contributed to an overall increase in aerosol mass under some of the test conditions. Using the results from the rented WESP, the project team installed a smaller-scale WESP upstream of the EERC’s slipstream CO 2 capture system. Flue gas was routed through a pilot-scale FGD unit to remove SO 2 to very low levels (~1 ppm) and then through a direct contact cooler (DCC) to further cool the gas and to remove moisture. The gas exiting the DCC was then routed through the new WESP before passing to the CO 2 absorber columns. Fluor’s amine-based solvent was used to scrub CO 2 from the slipstream through a set of two absorber columns. The rich solvent was regenerated in a stripper column by heating to drive off captured CO 2 . The system operated using a catch-and-release method where the CO 2 was separated to provide data on the process, but the captured CO 2 was released back into the host site stack. Aerosols and sulfur species were measured at multiple locations throughout the pilot-scale system. The inlet FGD and DCC removed much of the particulate matter and gaseous sulfur upstream of the WESP. With this configuration, the WESP achieved >95% particulate capture. The new WESP did not show any of the increases in SO 3 or other aerosol species that had been consistently observed with the larger-scale WESP installed immediately downstream of the plant’s full-scale FGD unit. Particulate samples captured and analyzed from the WESP inlet did not show any presence of crystalline sulfate materials, indicating that the pilot-scale FGD and DCC were efficient at reducing carryover from the plant’s full-scale FGD unit. A set of parametric tests were conducted on the new WESP to assess the impacts of flow rate, voltage, number of online WESP fields, and gas-phase sulfur content on aerosol and sulfur transformations. The results showed that the WESP performed similarly well at all sets of conditions. Sulfur and particulate matter exiting the WESP were further reduced through the absorber column as the amine-based solvent captured some of the residual contaminants. Solvent analysis showed that these species were slowly concentrating in the solvent over the duration of the test. When the WESP was taken offline and the sulfur slip through the FGD allowed to rise, the sulfate content in the solvent rose sharply, showing that the extra FGD and WESP were effective at reducing sulfate and cation uptake. This would be expected to extend amine-based solvent life by slowing the formation of heat-stable salts and other degradation products. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission and MPC.

01 COAL, LIGNITE, AND PEAT↗

Twisted bilayer graphene as topological heavy fermion: II. Analytical approximations of the model parameters

The recently-introduced topological heavy fermion (THF) model [1] of twisted bilayer graphene (TBG) aims to reconcile the quantum-dot-like electronic structure of the latter observed by scanning tunneling microscopy, with its electron delocalization seen in transport measurements. The THF model achieves this by coupling localized (heavy) fermions with anomalous conduction electrons. Originally, the parameters of the THF model were obtained numerically from the Bistritzer–Macdonald (BM) model of TBG [1]. In this work, we derive analytical expressions for the THF model parameters as a function of the twist angle, the ratio between the tunneling amplitudes at the AA and AB regions (w 0 /w 1 ), and the screening length of the interaction potential. By numerically computing the THF model parameters across an extensive experimentally-relevant parameter space, we show that the resulting approximations are remarkably good, i.e., within the 30% relative error for almost the entire parameter space. At the single-particle level, the THF model accurately captures the energy spectrum of the BM model over a large phase space of angles and tunneling amplitude ratios. When interactions are included, we also show that the THF description of TBG is good around the magic angle for realistic values of the tunneling amplitude ratios (0.6 ≤ w 0 /w 1 ≤ 1.0), for which the hybridization between the localized and conduction fermions γ is smaller than the onsite repulsion of the heavy fermions U 1 (i.e., |γ| < U 1 ).

Approximation methods↗

General kinetic ion-induced electron emission model for metallic walls applied to biased Z-pinch electrodes

A kinetic ion-induced electron emission (IIEE) model for general applications is developed to obtain the emitted electron energy spectrum for a distribution of ion impacts on a metallic surface. We assume an ionization cascade mechanism and use empirical models for the ion and electron stopping powers. The emission spectrum and the secondary electron yield (SEY) are validated for a variety of materials. The IIEE model is used to study the effect of IIEE on the plasma-material interactions of Z-pinch electrodes. Un-magnetized Boltzmann-Poisson simulations are performed for a Z-pinch plasma doubly bounded by two biased copper electrodes with and without IIEE at bias potentials from 0 to 9 kV. At the anode, the SEY decreases from 0 to 1 kV, but then increases at higher bias potentials. At the cathode, the SEY is much larger due to higher energy ion bombardment and grows with bias potential. As the bias potential increases, the emitted cathode electrons are accelerated to higher energies into the domain, collisionally heating the plasma. Above 1 kV, the heating is strong enough to increase the plasma potential. Despite SEY greater than 1, only a classical sheath forms as opposed to a space-charge limited or inverse sheath due to the emitted electron flux not reaching the space charge current saturation limits. Furthermore, the current in the emissionless cases saturates to a value lower than experiment. With IIEE, the current does not saturate and continues to increase with the 4 kV case, matching most closely with the experiment.

Carbon based materials↗

Correlating Protein Dynamics and Catalytic Activity of a Model Hydrogenase Using Paramagnetic and Biological Nuclear Magnetic Resonance Spectroscopy

Rational catalyst design remains a significant challenge, with electronic structure, steric, and electrostatic effects known to contribute to activity. Recently, dynamics has been recognized as another factor that impacts catalysis, though identifying and predicting these effects has remained out of reach. Nickel-substituted rubredoxin (NiRd), a protein-based mimic of a hydrogenase enzyme, serves as a model catalytic system in which dynamics can be systematically investigated with respect to activity. While over 30 secondary-sphere mutants of NiRd have been shown to be catalytically active, no significant correlation was observed between the rates and catalytic overpotential or electronic structure, prompting questions about the protein-derived factors that modulate activity. Here, in this work, NMR spectroscopy was used to investigate the roles of substrate accessibility, protein dynamics, and protein stability in controlling catalysis. Significant paramagnetic effects from the nickel center (S = 1) isolate the methylene proton resonances of the metal-coordinating cysteine residues. The sensitivity of resonance positions and linewidths to local environment offers an opportunity to study dynamical molecular changes around the metal center with high resolution. Machine learning algorithms were employed to identify correlations between the catalytic activity and the paramagnetic NMR spectra. These analyses revealed spectroscopic features of specific cysteine protons that report on catalytic overpotential and increased turnover rates, which are further supported by the results obtained using high-field NMR techniques. Collectively, these studies indicate the potential for multifrequency NMR techniques to resolve key contributors to catalytic activity and highlight the importance of local and outer-sphere dynamics.

Protein Engineering↗

Topical Report – Findings on Subtask 2.7 – Wet ESP and Aerosol Testing at Coal Creek Station

Growing concerns over the impact of CO 2 emissions from combustion sources on global climate change have prompted numerous research and development projects aimed at developing cost-effective technologies for CO 2 capture. One family of technologies being demonstrated at pilot and full scale globally is postcombustion carbon capture (PCCC) systems that employ amine-based solvents. The captured CO 2 can be compressed and permanently stored underground or used for enhanced oil recovery. The proximity of North Dakota’s lignite-fired fleet of power plants to potential CO 2 storage options creates a unique atmosphere for PCCC within the state. However, the unique components present in lignite flue gas present a challenge for large-scale PCCC at North Dakota power plants by contributing to aerosol formation. Aerosols can negatively impact the long-term performance of amine-based solvents for CO 2 capture. Amine-based solvents are volatile, and flue gas particulate provides nucleation sites where amine vapors can condense as aerosols. Because aerosols cannot be easily captured at the column outlet using conventional technologies, the amine-laden aerosols escape the system and lead to accelerated solvent losses. Moreover, particulate components can chemically react with amines to form degradation products that can permanently deactivate the amine, cause fouling, and lead to hazardous emissions. Many of the elements that have been shown to catalyze solvent degradation are present in lignite coals and can exacerbate solvent replacement economics. Understanding this issue is critical to the implementation of solvent-based CO 2 capture systems as applied to lignite-fired generation systems. The Energy & Environmental Research Center (EERC) designed and carried out this project to fully characterize aerosol behavior with various control technologies installed to better optimize aerosol mitigation technology for CO 2 capture. To meet the goal of this project, the following objectives were identified: Determine the effectiveness of a wet electrostatic precipitator (WESP) on mitigating formation of problematic aerosols at Great River Energy’s Coal Creek Station, upstream of the PCCC system. Determine the effectiveness of the Mitsubishi Heavy Industries (MHI) proprietary amine emission reduction unit (AERU) as a postcapture solvent recovery system for reducing aerosol emissions and extending solvent life downstream of the PCCC system. Determine the impact of aerosols on the efficiency and degradation products of both commercial and advanced solvents within the PCCC system. Work was conducted at Coal Creek Station Unit 1 using a slipstream of flue gas from the outlet of the plant’s flue gas desulfurization (FGD) unit. Flue gas was routed through a pilot-scale FGD unit to remove SO 2 to very low levels (~1 ppm) and then through a direct contact cooler (DCC) to further cool the gas and to remove moisture. The gas exiting the DCC was then optionally routed through a WESP before passing to the CO 2 absorber columns. The MHI solvent was used to scrub CO 2 from the slipstream through a set of two absorber columns. The rich solvent was regenerated in a stripper column by heating to drive off captured CO 2 . Flue gas exiting the absorber column was routed to MHI’s proprietary AERU to recover entrained solvent. The system operated using a catch-and-release method where the CO 2 was separated to provide data on the process, but the captured CO 2 was released back into the host site stack. Particulate was measured, collected, and analyzed from multiple locations throughout the pilot-scale system. Unlike the performance observed in prior work, the inlet FGD and DCC did not remove significant particulate matter from the flue gas. This appears to be due to a difference in the nature of the particulate. The DCC seemed to increase particulate size and count, most likely owing to water condensing onto the surfaces of fly ash particles. When the WESP was operated, it achieved >95% particulate capture. Very little particulate matter or indications of solvent were detected at the AERU outlet. When operating with advanced KS 21 solvent, the particulate material at the AERU outlet was even further decreased. Solvent analysis showed that some species derived from flue gas and ash were slowly concentrating in the solvent over the duration of the test. The levels observed were reported to be within expected ranges and were not of concern to MHI. A high-level techno-economic assessment of installing CO 2 capture at Coal Creek Station suggested that, when using a standard monoethanolamine (MEA)-based solution with simple heat integration, the energy penalty to net generation would be 34%. The bulk of this was due to steam losses for regenerating solvent, followed by parasitic electrical demand for CO 2 compression and then by increased parasitic load for pushing flue gas through the absorber column. These demands could be decreased with a more advanced solvent that exhibits lower heat of regeneration and lower pressure drop than does a simple MEA solution. Further energy could be saved with more thorough heat integration to recover useful energy from the steam used for solvent regeneration. Installing a WESP was predicted to increase the cost of electricity by nearly $5/MWh. This would become cost-effective if solvent losses were roughly 10 times the baseline estimate when not using a WESP but could be returned to baseline by installing the WESP. Piping CO 2 for storage in more favorable geology could help with carbon capture and storage economics. Although storing off-site would necessitate construction of a CO 2 transport pipeline, the cost of this pipeline might be more than offset by reducing the number of wells required, the depths of the wells required, and the electrical demand for the CO 2 compressor. Additional factors that favor off-site storage costs include smaller expected CO 2 plume sizes, which translates to less monitoring and fewer landowner agreements. More detailed assessment of the specific geology in the region under and around Coal Creek Station would be needed to accurately assess the costs and benefits of different storage site options. This subtask was cofunded through the Energy & Environmental Research Center–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission.

Strege, Joshua↗

Roles for Countercharge in the Voltage Sensor Domain of Ion Channels

Voltage-gated ion channels share a common structure typified by peripheral, voltage sensor domains. Their S4 segments respond to alteration in membrane potential with translocation coupled to ion permeation through a central pore domain. The mechanisms of gating in these channels have been intensely studied using pioneering methods such as measurement of charge displacement across a membrane, sequencing of genes coding for voltage-gated ion channels, and the development of all-atom molecular dynamics simulations using structural information from prokaryotic and eukaryotic channel proteins. One aspect of this work has been the description of the role of conserved negative countercharges in S1, S2, and S3 transmembrane segments to promote sequential salt-bridge formation with positively charged residues in S4 segments. These interactions facilitate S4 translocation through the lipid bilayer. In this review, we describe functional and computational work investigating the role of these countercharges in S4 translocation, voltage sensor domain hydration, and in diseases resulting from countercharge mutations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗

Structure and Dynamics of Domains in Ferroelectric Nanostructures – Phase-Field Modeling (Final Technical Report)

Domain pattern formation is one of the most common phenomena in nature and is a topic of immense interest in many fields ranging from materials science and physics to chemistry and biology. This DOE sponsored research project explored the basic science concerning the thermodynamic stability of mesoscale polarization domain patterns and their temporal evolution mechanisms during formation and subsequent switching in ferroelectric nanostructures and heterostructures. The project employed the computational phase-field method in combination of microelasticity and electrostatic theories. The research was carried out in close collaboration with a number of experimental groups who used High Resolution Transmission Electron Microscopy (HRTEM), in situ TEM with Scanning Probe Microscopy (SPM), or Piezoresponse Force Microscopy (PFM) to characterize the domain structures and dynamics of ferroelectric thin films and heterostructures and who grow high-quality ferroelectric and multiferroic thin films using advanced growth techniques such as Molecular Beam Epitaxy (MBE), Pulsed Laser Deposition (PLLD), and sputtering. The research efforts of the project help establish the phase-field method as the most powerful method for understanding and predicting domain structures in ferroelectric thin films and nanostructures. The findings of the project led to the basic understanding of stability and switching mechanisms of ferroelectric domains under different mechanical boundary conditions and under either homogeneous capacitor configurations or local fields using metallic probes as electrodes and the emergence of charged domain walls during domain switching. The project predicted the spatial length scales, temperature ranges, and electromechanical conditions for different polar states in heterostructures and guided the discovery of both transient and stable novel polarization states containing vortex lattices in oxide superlattices. The project resulted in 134 journal publications and 6 PhD theses with all the PhD graduates currently working in either academia or industry within the United States. The basic understanding on the stability of mesoscale polar states and pattern evolution achieved by the project improved our ability to control and engineer properties of ferroelectric thin films and heterostructures for potential applications in nanoscale electronic devices.

36 MATERIALS SCIENCE↗

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Foundations of Molecular 'Isotomics'

The naturally occurring rare isotopes are versions of common elements, such as hydrogen, carbon and oxygen, that contain a larger than usual number of neutrons in their atomic nuclei and therefore are higher in mass than the common atoms of that element. Isotopes exist for most elements and are found in most natural and synthetic materials, but are uneven in their distribution because chemical and physical processes are isotope-selective (e.g., a chemical reaction may proceed more rapidly for one isotope than for another). For this reason, abundances of isotopes in a material of interest can provide a record, or ‘signature’ of various features of that material’s origin and history. These signatures have been used in the geo, life, chemical and physical sciences in a wide variety of ways over close to 8 decades. However, many such applications struggle to reach unique interpretations of isotopic data because multiple factors combine to control a given sample’s overall isotopic content. That is, the factors controlling isotopic content are too numerous and complex to fully constrain from a simple measurement of a material’s isotope abundances. However, the distribution of isotopes within materials, at molecular scales potentially provides a vastly larger number and diversity of constraints on the chemical and physical processes that comprise a material’s history. The rare isotopes may be concentrated into one atomic position in a molecule relative to another, some proportion of molecules in a sample may contain two or more rare isotopes, and those multiply-isotope-substituted forms of molecules may also have uneven distributions of those isotopes across individual atomic sites. For these reasons, even small, seemingly simple molecules, such as sugars, amino acids or drug compounds, actually exist in a vast number of isotopically unique forms (often millions or more), and each one of those forms is in some sense an independent ‘vote’ on that sample’s history. This project has focused on opening this rich archive of information by enabling the creation of routinely and widely applicable ways of measuring and interpreting isotopic structures of molecules. This work has included the development of core technologies and analytical methods, advancing fundamental understanding of the physical and chemical properties of isotopic versions of molecules, and conducting proof of concept studies of illustrative geochemical, cosmochemical and forensic problems in order to show how these technologies, methods and principles come together to solve problems in new ways. A key to the success of this project was the adaptation of ‘Fourier transform mass spectrometry’ (FTMS) to the task of precisely measuring proportions of the rare, naturally occurring isotopic forms of molecules. FTMS is a highly specialized form of mass spectrometry that traps ions within magnetic or electrostatic cavities and, effectively, ‘listens’ (through registering of subtle electrical signals) to the harmonic signals they make while rapidly orbiting within those cavities. These signals have periods that are a function of their mass and strength (or ‘loudness’) that is proportional to their abundances. Thus, these signals constrain relative amounts of molecules that differ in their mass due to various isotopic substitutions. This technology has been essential to the identification of organic molecules in the life, chemical and environmental sciences for over 4 decades, but generally has lacked the control, stability and precision to meaningfully measure rare isotope forms of molecules. This project’s most fundamental contribution has been to modify FTMS, both in terms of hardware and methods, to enable such measurements. The raw data of molecular isotopic structure is tremendously voluminous and complex, so another important activity of this project has been developing the theoretical and data-science tools needed to interpret the data generated by this new form of isotopic measurement. A particularly challenging part of this task has been predicting molecular isotopic structure, as only through the comparison of measurements with predictions can we make progress on hypothesis driven research questions. We have attacked this this prediction task through a combination of first-principles chemical-physics models of the effects of isotope substitution on molecule properties and data-science models that permit us to generalize that chemical physics to cases that have not yet been studied by detailed chemical physics theory. The proof of concept applications we have pursued over the course of this study include biological reactions of amino acids and other biomolecules, non-biological synthesis of organic molecules in extra-terrestrial settings such as meteorites, petroleum geoscience questions concerning the origin and evolution of natural gas, oil and kerogen compounds, and forensic questions such as the sourcing of chemical weapons. The successes of these applications have laid the groundwork for the next phase of this field’s development, which will include larger scale and more ambitious studies of molecular isotopic structure as a means of diagnosing human diseases, such as cancer, and reconstructing detailed interpretations of the origin and evolution of organic molecules in modern and geological environments.

Cesar, Jaime↗

Few-mode and anisotropic quantum transport in InSb nanoribbons using an all-van der Waals material-based gate

High-quality electrostatic gating is a fundamental ingredient for successful semiconducting device physics, and a key element of realizing clean quantum transport. Inspired by the widespread improvement of transport quality when two-dimensional van der Waals (vdW) materials are gated exclusively by other vdW materials, we have developed a method for gating non-vdW materials with an all-vdW gate stack, consisting of a hexagonal boron nitride dielectric layer and a few-layer graphite gate electrode. We demonstrate this gating approach on MOVPE-grown InSb nanoribbons (NRs), a novel variant of the InSb nanowire, with a flattened cross-section. In our all-vdW gated NR devices we observe conductance features that are reproducible and have low- to near-zero gate hysteresis. We also report quantized conductance, which persists to lower magnetic fields and longer channel lengths than typical InSb nanowire devices reported to date. Additionally, we observe level splitting that is highly anisotropic in an applied magnetic field, which we attribute to the ribbon cross-section. The performance of our devices is consistent with the reduced disorder expected from the all-vdW gating scheme, and marks the first report of ballistic, few-modes quantum transport in a non-vdW material with an all-vdW gate. Our results establish all-vdW gating as a promising approach for high-quality gating of non-vdW materials for quantum transport, which is in principle applicable generically, beyond InSb systems. In addition, the work showcases the specific potential of all-vdW gate/InSb NR devices for enabling clean quantum devices that may be relevant for spintronics and topological superconductivity studies.

InSb↗

Laser aberration signatures in expelled electrons from a tenuous gas

We describe a mechanism by which the aberration content of a focused, multiterawatt laser is imprinted upon the forward angular distribution of electrons ionized within and ponderomotively expelled from the focal volume. In our experiments, the laser aberration type and magnitude are controllably varied, and the measured electron distributions are correspondingly modified in a way consistent with predictions from numerical simulations. This imprint mechanism shows potential for enabling the development of accurate focal-spot characterization of intense lasers when fired at full power and is being developed for extension to petawatt laser systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of an Ultrahigh-bandwidth Phase Contrast Imaging System for detection of electron scale turbulence and Gigahertz Radio-Frequency Waves

The study of waves and turbulence is vital to the development of future reactor-grade plasma devices developed in the quest for fusion energy. These fluctuations are responsible for moving heat and particles across the magnetic field, and a predictive understanding of them is needed to achieve the density and temperature required to sustain a plasma fusion reaction. While many techniques have been developed for measuring waves and fluctuations, every measurement method has limitations. There are relatively few techniques for measuring very high frequency fluctuations, such as radio frequency waves injected to heat the plasma, unstable waves driven by suprathermal particles, or short wavelength electrostatic waves driven by electron temperature or density gradients. The present project builds upon the proven phase contrast imaging (PCI) technique to extend the response of the diagnostic by orders of magnitude in frequency and almost a factor of ten in spatial resolution. PCI provides a measurement of electron density based on small angle scattering of a CO-2 laser beam by using optical techniques to render a phase shift as an intensity change on a detector. Due to the telecommunications revolution, technological development by manufacturers has focused on components in the near -infrared, so that high- quality lasers and detectors at 1.55 µm are readily available. Shifting PCI design to a new, shorter wavelength has numerous advantages and challenges. Similar detector performance is available with room-temperature arrays with GHz bandwidth, while the detector arrays for 10.6 µm required liquid nitrogen cooling and were therefore limited to a bandwidth of about 1 MHz. Shifting to a shorter wavelength reduces the angle at which the laser beam scatters off of plasma waves, allowing more such scattered components to pass through the aperture of the vacuum vessel port and be collected by the PCI, which increases the spatial resolution. Concomitant with these benefits, various questions of performance arise. At shorter laser wave- length, the sensitivity to mirror and lens quality is increased, the contribution of the laser to the overall system noise is increased, and the sensitivity to vibrations is increased. The custom optical components at the heart of the PCI technique were required to be properly scaled for the shorter wavelength, so fabrication technologies needed to be explored. This project was designed to show that a low noise, high response PCI system at 1.55 µm could be constructed and operated, and then to quantify potential issues to allow extrapolation to a full-size production system providing physics measurements on a large plasma device. The first stage, producing the custom optical component called a Phase Plate, was successfully achieved using two methods. First, an easily reproduced masking and coating technique was able to produce good phase plates with the required parameters. Second, a nanofabrication technique was found to produce extremely high quality phase plates at a competitive price. The PCI constructed with the new phase plates and 1.55 µm laser was found to provide excellent wavelength measurements with the theoretically expected response. The sensitivity to optical surface quality was found to be in line with previous measurements at 10.6 µm. The observed signal-to-noise ratio was similar to the theoretically expected value. The effect of vibrations on PCI was studied with the first measurement of the effect of beam motion on PCI response and comparison to theory, allowing for a quantitative prediction of the effect of vibrations on a production PCI system and the requirements for improved beam stabilization. PCI is an extremely cost-effective method to provide a low noise, absolutely calibrated measurement of plasma fluctuations across a wide spatial scale. This project has shown that a 1.55 µm PCI using modern techniques and components is less expensive than the 10.6 µm PCI of a few years ago, with the largest savings in phase plate fabrication and the infrared detector array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗