Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Polar system modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

All-epitaxial long-range surface plasmon polariton structures with integrated active materials

We demonstrate all-epitaxial structures capable of supporting short- and long-range surface plasmon polariton (SRSPP and LRSPP) modes in the long-wave infrared region of the electromagnetic spectrum. The SRSPP and LRSPP modes are bound to the interfaces of a buried heavily doped (n ++ ) semiconductor layer and surrounding quantum-engineered type-II superlattice (T2SL) materials. The surrounding T2SLs are designed to allow optical transitions across the frequency dispersion of the SPP modes. We map the SPP dispersion in our structure using grating-coupled angle- and polarization-dependent reflection and photoluminescence spectroscopy. The epitaxial structures are analytically described using a simplified three-layer system (T2SL/n ++ /T2SL) and modeled using rigorous coupled wave analysis with excellent agreement to our experimental results. The presented structures offer the potential to serve as long-range interconnects or waveguides in all-epitaxial plasmonic/optoelectronic systems operating in the long-wave infrared.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gapless to gapless phase transitions in quantum spin chains

We investigate spin chains with bilinear-biquadratic (BLBQ) spin interactions as a function of an applied magnetic field h. At the Uimin-Lai-Sutherland (ULS) critical point we find a gapless to gapless transition revealed by the dynamical structure factor S(q,ω) as a function of h. At h=0, the envelope of the lowest-energy excitations goes soft at two points, q 1 =2π/3 and q 2 =4π/3, dubbed the phase A. With increasing field, the spectral peaks at each of the gapless points bifurcate, making in total four soft modes, and combine to form a new set of excitations that soften at a single point q=π at h c1 ≈0.94. Beyond h c1 the system enters another gapless B phase until the transition at h c2 =4 to the fully polarized phase. We compare the ULS model results with those for the Affleck-Kennedy-Lieb-Tasaki model as a representative of the gapped Haldane phase. We explain the mechanism of the gapless to gapless transition in the ULS model using its conserved charges and a spinon band picture. We also discuss the universality of central charges of the BLBQ family of models subjected to a magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Elucidating Electric Field-Induced Rate Promotion of Brønsted Acid-Catalyzed Alcohol Dehydration

Applied potentials have been demonstrated as a powerful tool to promote heterogeneous Brønsted acid catalysis by orders of magnitude, leveraging interfacial electric fields to stabilize protonated intermediates. However, the use of flat two-dimensional electrodes with inherently low active site densities limits the application of conventional thermochemical characterization techniques that can probe the nature of catalytic active sites. Here, we use kinetic analyses with an electrostatics-based model to elucidate the intricacies of potential-induced rate promotion, employing liquid-phase dehydration of 1-methylcyclopentanol catalyzed by carboxylic acid groups on carbon nanotubes as a probe system. By using a basket electrode to directly polarize catalyst powder, we demonstrate that thermocatalytic reaction rates can be promoted by 100,000-fold, exhibiting a log–linear dependence on applied potential with rate-potential scalings as high as 125 ± 4 mV per 10-fold rate increase. In agreement with model predictions, we show that lower ionic strengths attenuate potential sensitivity, resulting from a weakening of the interfacial electric field that interacts with the acidic proton. Furthermore, we experimentally confirm the model-predicted “isokinetic potential” (at ∼0.6 V vs Ag/AgCl)─the potential at which all rate scaling lines at various ionic strengths intersect, making the rate independent of ionic strength. Base titrations reveal that only ∼8% of the carboxylic acid sites are catalytically active, yet these same active sites are operational at the highest and lowest potentials. Collectively, our results provide a key methodology for modeling catalytic effects of electric fields, quantifying active sites under applied potential, and demonstrating fundamental principles of electric field-induced rate promotion.

Catalysts↗
DOE OSTI2023

Development of polarization-sensitive optical coherence tomography imaging platform and metrics to quantify electrostimulation-induced peripheral nerve injury in vivo in a small animal model

Neuromodulation devices are rapidly evolving for the treatment of neurological diseases and conditions. Injury from implantation or long-term use without obvious functional losses is often only detectable through terminal histology. New technologies are needed that assess the peripheral nervous system (PNS) under normal and diseased or injured conditions. We aim to demonstrate an imaging and stimulation platform that can elucidate the biological mechanisms and impacts of neurostimulation in the PNS and apply it to the sciatic nerve to extract imaging metrics indicating electrical overstimulation. A sciatic nerve injury model in a 15-rat cohort was observed using a newly developed imaging and stimulation platform that can detect electrical overstimulation effects with polarization-sensitive optical coherence tomography. The sciatic nerve was electrically stimulated using a custom-developed nerve holder with embedded electrodes for 1 h, followed by a 1-h recovery period, delivered at above-threshold Shannon model k-values in experimental groups: sham control (SC, n = 5, 0.0 mA / 0 Hz), stimulation level 1 (SL1, n = 5, 3.4 mA / 50 Hz, and k = 2.57), and stimulation level 2 (SL2, n = 5, 6.8 mA / 100 Hz, and k = 3.17). The stimulation and imaging system successfully captured study data across the cohort. When compared to a SC after a 1-week recovery, the fascicle closest to the stimulation lead showed an average change of +4 % / – 309 % (SL1/SL2) in phase retardation and –79 % / – 148 % in optical attenuation relative to SC. Analysis of immunohistochemistry (IHC) shows a +1 % / – 36 % difference in myelin pixel counts and –13 % / + 29 % difference in axon pixel counts, and an overall increase in cell nuclei pixel count of +20 % / + 35 % . These metrics were consistent with IHC and hematoxylin/eosin tissue section analysis. The poststimulation changes observed in our study are manifestations of nerve injury and repair, specifically degeneration and angiogenesis. Optical imaging metrics quantify these processes and may help evaluate the safety and efficacy of neuromodulation devices.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Cathodic Protection Modeling for Hanford Underground Double-Shell Tank Farms

Hanford stores millions of gallons of radioactive and chemically hazardous waste from the production of weapon materials in tank farms consisting of underground carbon-steel storage tanks surrounded by reinforced concrete. Six of these Hanford tank farms use double-shell storage tanks (DSTs). The DST farms were constructed from 1968 to 1986 with a planned 40–50 year design life, so some are already operating beyond their initial life expectancy. Ultrasonic testing (UT) has indicated significant thinning on the bottom of the secondary (outer) liner of these tanks, believed to arise from groundwater intrusion driving concrete side corrosion. There is no direct access to the steel/concrete interface between the tank and the concrete pad, making it difficult to apply a chemical-based mitigation strategy or to conduct repairs, but cathodic protection (CP) is a possible method to inhibit further concrete-side corrosion. Hanford already uses CP to protect below grade steel piping within the tank farms and connected to the tanks, but this system was not designed to protect the tank bottoms. CP design must account for the structures surrounding the DSTs, including the steel reinforcing bars (rebar) within the concrete pad and vault, various process lines, and the existing CP system. In this study, finite element analysis (FEA) modeling was carried out to simulate CP protection of 1) a single tank and CP anode to develop options for modeling the rebar and to compare to a simpler circuit model and 2) the entire Hanford AN tank farm as a representative example consisting of seven tanks, associated piping, and both existing and new CP anodes. Both circuit and FEA models predict that significant protective current could be delivered to the bottoms of the tanks with the addition of tank-protection anodes below the depth of the tanks. Simulations with only the existing pipe-protection anodes active confirmed that only a very small current to the tank bottoms is predicted under present conditions. Multiple simplified representations of the dome and wall rebar were tested to reduce the computational complexity of the tank-farm simulations, resulting in modeling the rebar as edge elements with a prescribed effective circumference that matches the real rebar surface area. The geometry of the rebar is also simplified into horizontal hoops around the tank walls and radial rebar over the dome with increased effective circumference to retain the target surface area. This simplification was found to greatly reduce the complexity and solution time of the models without large changes in current distributions, especially to the tank bottom. A range of values were tested for model parameters such as soil and concrete resistivities and polarization resistance to investigate their impact on the current and electric potential distributions. Depending on the parameters used, FEA simulations predict some risk of overprotection, particularly on the piping system; since overprotection can also lead to surface damage associated with hydrogen gas generation at the interface (e.g. hydrogen embrittlement or damage to coatings), this needs to be considered when refining the design of the new CP system. Comparison between the FEA models and the circuit model representation demonstrated that the circuit model could not match the predicted FEA current distribution, even when using the exact same surface areas. This discrepancy appeared to be at least partly attributable to the impact of the relative positions of the tank components and anodes to each other and to the ground surface. The FEA model accounts for the relative positions since it solves the governing equations in three dimensions, but the circuit model cannot account for the positioning. In particular, the circuit model underpredicts the current to the tank bottom and overpredicts the current to the dome compared to FEA for the baseline geometry. The FEA models omitted the electrically isolated rebar in the bottom concrete slab. However, a circuit based stray current model estimated that only 2.1% of the total current through the slab would stray into the rebar, corresponding to ~0.21 A for a target current density of 2 mA/ft2 to the tank bottom. The estimated corrosion driven by this amount of stray current is predicted to yield a lifetime of >400 years for the minimum rebar diameter, assuming an acceptable cross-section area loss of 10%.

d'Entremont, Anna [Savannah River National Laborat↗
DOE OSTI2025

Emergent topological polarization textures in relaxor ferroelectrics

Relaxor ferroelectrics underpin high-performance actuators and sensors, yet the nature of polar heterogeneities driving their broadband dielectric response remains debated. Using a unified, multimodal structural refinement framework— simultaneously fitting complementary X-ray and neutron total scattering, X-ray absorption spectra, and diffuse scattering—we reconstruct 3D mesoscale polarization maps in the classic relaxor system PbMg 1/3 Nb 2/3 O 3 –PbTiO 3 . We uncover self-organized swirling polarization textures with half-skyrmion (meron) vortices, challenging models of independent polar nanoregions. These textures, characterized by smooth changes in the polarization direction, originate from overlapping volumes in which the projections of locally correlated polarization vectors onto each volume’s long axis share the same sign. Vortex cores correlate strongly with local charge and strain gradients imposed by compositional heterogeneities. In this work, our results suggest that chemical disorder, acting via depolarizing and strain fields, stabilizes topological vortex textures of the polarization field, offering a route for engineering new dielectric and ferroelectric functionalities.

Eremenko, Maksim [National Institute of Standards ↗
DOE OSTI2021

Enhanced electro-osmosis in propylene carbonate salt solutions

Properties of solid–liquid interfaces and surface charge characteristics mediate ionic and molecular transport through porous systems, affecting many processes such as separations. Herein, we report experiments designed to probe the electrochemical properties of solid–liquid interfaces using a model system of a single polyethylene terephthalate (PET) pore in contact with aqueous and propylene carbonate solutions of LiClO 4 . First, the existence and polarity of surface charges were inferred from current–voltage curves recorded when a pore was placed in contact with a LiClO4 concentration gradient. Second, the electro-osmotic transport of uncharged polystyrene particles through the PET pore provided information on the polarity and the magnitude of the pore walls’ zeta potential. Our experiments show that the PET pores become effectively positively charged when in contact with LiClO 4 solutions in propylene carbonate, even though in aqueous LiClO 4 , the same pores are negatively charged. Additionally, the electro-osmotic velocity of the particles revealed a significantly higher magnitude of the positive zeta potential of the pores in propylene carbonate compared to the magnitude of the negative zeta potential in water. The presented methods of probing the properties of solid–liquid interfaces are expected to be applicable to a wide variety of solid and liquid systems.

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

Changes in environmental and engineered conditions alter the plasma membrane lipidome of fractured shale bacteria

ABSTRACT Microorganisms that persist in fractured shale reservoirs cause several problems including secreting foul gases and forming biofilms. Current biocontrol measures often fail due to limited knowledge of their in situ activities. The plasma membrane protects the cell, mediates many of its critical functions, and responds to intracellular cues and ecological perturbations through physicochemical modifications. As such, it provides valuable insight into the physiological adaptation of microorganisms in disturbed environmental systems. Here, we (i) demonstrate how changes in salinity and hydraulic retention time (HRT) influence the plasma membrane intact polar lipid (IPL) chemistry of model bacterium, Halanaerobium congolense WG10, and mixed microbial consortia enriched from shale-produced fluids and (ii) elucidate adjustments in membrane IPL chemistry during biofilm growth relative to planktonic cells. We incubated H. congolense WG10 in chemostats under three salinities (7%, 13%, and 20% NaCl), operated under three HRTs (19.2, 24, and 48 h), and in drip flow biofilm reactors under the same salinity gradients. Also, mixed microbial consortia in produced fluids were enriched in triplicate chemostat vessels under three HRTs (19.2, 24, and 72 h) and biofilm reactors. Lipids were analyzed by ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Our results show that phosphatidylglycerols, cardiolipins, and phosphatidylethanolamines were predominantly enriched in planktonic H. congolense WG10 cells grown at hypersalinity (20%) compared to optimum (13%). In addition, several zwitterionic phosphatidylcholines and phosphatidylethanolamines were higher in abundance during biofilm growth. These observations suggest that microbial adaptation and biofilm formation in fractured shale are enabled by strategic plasma membrane IPL chemistry adjustments. IMPORTANCE Microorganisms inadvertently introduced into the shale reservoir during fracturing face multiple stressors including brine-level salinities and starvation. However, some anaerobic halotolerant bacteria adapt and persist for long periods of time. They produce hydrogen sulfide, which sours the reservoir and corrodes engineering infrastructure. In addition, they form biofilms on rock matrices, which decrease shale permeability and clog fracture networks. These reduce well productivity and increase extraction costs. Under stress, microbes remodel their plasma membrane to optimize its roles in protection and mediating cellular processes such as signaling, transport, and energy metabolism. Hence, by observing changes in the membrane lipidome of model shale bacteria, Halanaerobium congolense WG10, and mixed consortia enriched from produced fluids under varying subsurface conditions and growth modes, we provide insight that advances our knowledge of the fractured shale biosystem. We also offer data-driven recommendations for improving biocontrol efficacy and the efficiency of energy recovery from unconventional formations.

03 NATURAL GAS↗
DOE OSTI2020

Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

Long range-corrected (LRC) or range-separated hybrid (RSH) functionals where the long-range (LR) limit of electronic interactions is set to the exact exchange have been shown to correct the tendency of traditional density functional theory (DFT) to underestimate the frontier orbital gap. Consequently, the use of such functionals in calculating electronic excited states using linear response based time-dependent DFT (TDDFT) has been successful in correcting the tendency for underestimating the energies of charge transfer states by DFT-based calculations. More recently formulations of functionals that attenuate the LR limit to address condensed-phase effects to polarize the electronic density have been reported. In particular screened RSH (SRSH) combined with polarizable continuum model (PCM) was benchmarked successfully in reproducing the fundamental gap and charge transfer state energies of molecular systems in the condensed phase. Here we use SRSH-PCM to address triplet excited states, and show its success in obtaining correspondence of the low-lying triplet states to the singlet–triplet gap in a similar way that the fundamental orbital gap corresponds to electron removal and addition energies. Importantly, the accuracy of the SRSH-PCM in calculating triplet excitations stands on the polarization consistent framework in addressing the scalar dielectric constant and without affecting the optimal tuning by triplet energies. Finally, the prospect of even further improving the SRSH-PCM accuracy in calculating triplet states can be achieved by optimal tuning on the basis of the spin multiplicity gap.

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

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

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

SQuaD: Smart Quantum Detection for Photon Recognition and Dark Count Elimination

Quantum detectors of single photons are an essential component for quantum information processing across computing, communication and networking. Today's quantum detection system, which consists of single photon detectors, timing electronics, control and data processing software, is primarily used for counting the number of single photon detection events. However, it is largely incapable of extracting other rich physical characteristics of the detected photons, such as their wavelengths, polarization states, photon numbers, or temporal waveforms. This work, for the first time, demonstrates a smart quantum detection system, SQuaD, which integrates a field programmable gate array (FPGA) with a neural network model, and is designed to recognize the features of photons and to eliminate detector dark-count. The SQuaD is a fully integrated quantum system with high timing-resolution data acquisition, onboard multi-scale data analysis, intelligent feature recognition and extraction, and feedback-driven system control. Our \name experimentally demonstrates 1) reliable photon counting on par with the state-of-the art commercial systems; 2) high-throughput data processing for each individual detection events; 3) efficient dark count recognition and elimination; 4) up to 100% accurate feature recognition of photon wavelength and polarization. Additionally, we deploy the SQuaD to an atomic (erbium ion) photon emitter source to realize noise-free control and readout of a spin qubit in the telecom band, enabling critical advances in quantum networks and distributed quantum information processing.

Linne, Karl C. [U. Chicago (main)] (ORCID:00090009↗
DOE OSTI2024

Optimal control of polar sea-ice near its tipping points

Abstract Several Earth system components are at a high risk of undergoing rapid, irreversible qualitative changes or “tipping” with increasing climate warming. It is therefore necessary to investigate the feasibility of arresting or even reversing the crossing of tipping thresholds. Here, we study feedback control of an idealized energy balance model (EBM) for Earth’s climate, which exhibits a “small icecap” instability responsible for a rapid transition to an ice-free climate under increasing greenhouse gas forcing. We develop an optimal control strategy for the EBM under different forcing scenarios to reverse sea-ice loss while minimizing costs. Control is achievable for this system, but the cost nearly quadruples once the system tips. While thermal inertia may delay tipping, leading to an overshoot of the critical forcing threshold, this leeway comes with a steep rise in requisite control once tipping occurs. Additionally, we find that the optimal control is localized in the polar region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Predictive Design of Hybrid Improper Ferroelectric Double Perovskite Oxides

The computational design of suitable multiferroic double perovskite oxides requires finding materials that exhibit sizable polarization, magnetization, and coupling between them. Oxides with the chemical formula of AA'BB'O 6 with building blocks of ABO 3 single perovskite oxides in centrosymmetric Pnma symmetry are strong candidates that have been reported to satisfy such criteria. The system lowers to noncentrosymmetric, polar P2 1 symmetry if A/A' layered and B/B' rocksalt cation orderings are imposed. A detailed compositional search over a variety of chemical spaces followed by evaluating their polarization may lead to the identification of more of these compounds with ferroelectric ordering. The standard density functional theory practices to estimate polarization within the Berry phase formalism require the systems to be perfectly insulating. The number of compounds that can be evaluated using this method is therefore limited. In this work, we introduce a predictive learning strategy based on importance sampling to build a series of machine learning models using results from first-principles simulations to predict polarization and the corresponding switching barrier. The geometry-driven features related to charge states and cationic radii play key roles in predicting the switching barrier with complementary contributions from the key structural mode-based order parameters. These modes become important to draw reasonable predictions of polarization components from machine learning models. In conclusion, our predictive models identify candidates with high polarizations and low switching barriers from a pool of double perovskite oxides, suitable for future investigation for their potential applications in spintronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Coriolis forces modify magnetostatic ponderomotive potentials

It is possible to produce a ponderomotive effect in a plasma system without time-varying fields, if the plasma flows over spatial oscillations in the field. This can be achieved by superimposing a spatially oscillatory perturbation on a guide field, then setting up an electric field perpendicular to the guide field to drive flow over the perturbation. However, subtle distinctions in the structure of the resulting electric field can entirely change the behavior of the resulting ponderomotive force. Previous work has shown that, in slab models, these distinctions can be explained in terms of the polarization of the effective wave that appears in the co-moving frame. Here, we consider what happens to this picture in a cylindrical system, where the transformation to the co-moving (rotating) frame is not inertial. It turns out that the non-inertial nature of this frame transformation can lead to counterintuitive behavior, partly due to the appearance of parallel (magnetic-field-aligned) electric fields in the rotating frame even in cases where none existed in the laboratory frame. Apart from the academic interest of this study, the practical impact lies in being better able to anticipate the antenna configuration on the plasma periphery of a cylindrical plasma that will lead to optimal ponderomotive barrier formation in the interior plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Non-volatile electric-field control of inversion symmetry

Competition between ground states at phase boundaries can lead to significant changes in properties under stimuli, particularly when these ground states have different crystal symmetries. A key challenge is to stabilize and control the coexistence of symmetry-distinct phases. Using BiFeO 3 layers confined between layers of dielectric TbScO 3 as a model system, we stabilize the mixed-phase coexistence of centrosymmetric and non-centrosymmetric BiFeO 3 phases at room temperature with antipolar, insulating and polar semiconducting behaviour, respectively. Application of orthogonal in-plane electric (polar) fields results in reversible non-volatile interconversion between the two phases, hence removing and introducing centrosymmetry. Counterintuitively, we find that an electric field 'erases' polarization, resulting from the anisotropy in octahedral tilts introduced by the interweaving TbScO 3 layers. Consequently, this interconversion between centrosymmetric and non-centrosymmetric phases generates changes in the non-linear optical response of over three orders of magnitude, resistivity of over five orders of magnitude and control of microscopic polar order. Here our work establishes a platform for cross-functional devices that take advantage of changes in optical, electrical and ferroic responses, and demonstrates octahedral tilts as an important order parameter in materials interface design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns

Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew–Burke–Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.

Chemistry↗
DOE OSTI2020

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electromagnetic Waves in Cold Plasmas

A collisionless fluid model of a plasma is presented and solved to first order in the cold, high frequency limit. From this model, a homogeneous system of linear equations is obtained for a harmonic electric field perturbation from which the Appleton-Hartree equation is derived. Moreover, the polarization states of electromagnetic waves in cold, collisionless plasmas are examined. Complete sets of equations governing the reflection and transmission of electromagnetic waves in both non magnetized and magnetized plasmas are obtained. In solving the boundary problem, the Booker quartic is derived and the nature of its roots discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗