Search NASA⌕ Search

SEARCH · Search NASA

Results for “force field development”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

The False Dilemma: Rethinking AF Science and Technology Officer Talent Management

In 2019, Secretary of the Air Force (SECAF) Heather Wilson launched the 2030 Science and Technology (S&T) Strategy by stating, “The advantage will go to those who create the best technologies and who integrate and field them in creative operational ways that provide military advantages.” In 2021, 15% of Air Force general officers responsible for creating and integrating this technological edge had a science, technology, engineering, and math (STEM) graduate degree, and less than 1% of generals had a STEM doctorate. As one point of comparison, at least 32% and 13% of founding CEOs of Fortune 500 technology companies had a master’s and doctoral STEM degree, respectively. The gulf is larger when you consider that these CEOs used their technical degrees throughout their careers whereas most general officers do not. In large organizations, if the leader does not possess technical knowledge, it is difficult to drive innovation systems and connect ideas to reality. The evidence for this disconnect is abundant in acquisition challenges for high-tech systems, significant pushback to the “revolution in military affairs”, and failure of the “Third Offset” to take hold and deliver capabilities to offset Chinese and Russian capabilities. Indeed, when the former Air Force and Space Force Chief Software Officer, Nicolas M. Chaillan, offered his resignation, one of the main reasons for leaving was “the failure of OSD and the Joint Staff to deliver on their own alleged top ‘priority’, JADC2 – they couldn’t ‘walk the walk.’” Senior leaders have recognized this disconnect, and several inquiries and studies regarding STEM competency have been conducted. Many have posed a false dilemma: a large, STEM-cognizant or small, STEM-expert force. However, current and future STEM human capital needs are managed at the unit and functional area level. This leads to “silos” and the tactical “needs of the present” dominating the strategic “needs of the future”, highlighted by only two functional areas systematically tracking future STEM needs. Therefore, the decision on force structure appears to have been historically made by default through these and other perceived structural realities, and a vast majority of recommendations have focused on the small STEM community – an exercise in diminishing returns. In essence, the solution often boils down to finding unicorns – officers that have chartered the unforgiving pathway through the traditional “gates” to general officer while obtaining sufficient STEM proficiency along the way. Fortunately for the Air Force, these officers exist, albeit at a rate below what the evidence would suggest is necessary. However, this approach fundamentally limits the ability to develop a deep pool of officers with the ability to create and integrate technologies and ensures the Dunning-Kruger effect is prevalent. The Dunning-Kruger effect is a cognitive bias where those with limited knowledge in an area lack the expertise necessary to recognize their lack of expertise and consequently are prone to overestimate their knowledge and performance. In short, we need less leaders overestimating their STEM knowledge and more with hard-earned STEM competency required to “walk the walk” to create and integrate technologies for military advantage. This will not be accomplished by restricting the emphasis to the traditional scientists and engineering (S&E) career fields – the pool, only 10% of line officers, and general officer progression is just too small. Fortunately, the recent SECAF’s Management Initiatives and Chief of Staff’s (CSAF) Action Orders, both containing emphasis on organic expertise to “accrue advantage in military-technological competition,” present an opportunity to truly develop a framework for the force of the future.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Solving the problem of overdetermination of quasisymmetric equilibrium solutions by near-axis expansions. I. Generalized force balance

It is well known that the process of construction of quasisymmetric magnetic fields in magnetostatic equilibrium with isotropic pressure suffers from the problem of overdetermination. This has led to the widespread belief that global quasisymmetric solutions are likely not to exist. Here we develop a general near-axis expansion procedure that does not rely on the assumption of magnetostatic equilibria with isotropic pressure. We then demonstrate that in equilibria with anisotropic pressure, it is possible to circumvent the problem of overdetermination and carry out the power-series solutions to higher order. This suggests, contrary to current belief, that the existence of globally quasisymmetric fields is likely if one relaxes the assumption of magnetostatic equilibria with isotropic pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nuclear-matter saturation and symmetry energy within Δ -full chiral effective field theory

Nuclear saturation and the symmetry energy are key properties of low-energy nuclear physics that depend on fine details of the nuclear interaction. The equation of state around saturation is also an important anchor for extrapolations to higher densities and studies of neutron stars. Here we develop a unified statistical framework that uses realistic nuclear forces to link the theoretical modeling of finite nuclei and infinite nuclear matter. We construct fast and accurate emulators for nuclear-matter observables and employ an iterative history-matching approach to explore and reduce the enormous parameter domain of Δ -full chiral interactions. We perform rigorous uncertainty quantification and find that model calibration including O 16 observables gives saturation predictions that are more precise than those that only use few-body data. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

1. Physical Security Engineering by Design for Nuclear Facilities; 2. Nuclear Power Plant Site Security Management – A Security Strategy; 3. The UAE Women in Nuclear Energy Security

1. Security by design, or SeBD, is a comprehensive approach that integrates the physical protection system of a nuclear plant into every stage of its existence. This includes planning, designing, constructing, commissioning, and operating the facility, using a combination of analytical, physical, technological, and procedural measures. Essentially, SeBD involves intentionally applying and incorporating security into all aspects of design and operation throughout the entire lifecycle of a facility. By implementing this methodology throughout various phases such as program development, process implementation, staff training and procedures management in conjunction with plant equipment, facilities can be optimized to minimize security risks without compromising functional design requirements. This ultimately improves the overall security posture of the site and reduces the need for costly modifications or additional security resources post-design. 2. A site security strategy is a living document that is revised on a periodic or event-driven basis, ensuring that site security operations and corresponding procedures provide long-term, effective protection for the entire nuclear power plant (NPP) site. A site security strategy aims to mitigate threats across the entire NPP site via in-depth defense approaches and mutual support; therefore, if a layer is omitted or altered, then the effect across all layers must be re-evaluated. Therefore, the aim of the site security strategy is to provide an appropriate, scalable security regime that deters, denies, delays, and detects incidents and, equally importantly, reassures legitimate users and the regulator that due diligence and regulatory compliance have been achieved, ensuring that the site is safe and secure. Robust access control for vehicles and pedestrians is at the heart of the strategy. Vehicle and pedestrian searching and screening are seen as the strongest mitigation methods against vehicle- and pedestrian-borne attacks. The security strategy must also be supported through comprehensive staff training and the development of robust processes, procedures, and planning. If all these measures are to be effective, then training must be implemented during each phase of construction, partial operation/commissioning, and full operation. No single element of site security is completely isolated from the influence of other elements. Ideally, consideration of all key elements will result in a security strategy that is integrated and proportional to the threat and that does not over specify individual security solutions through the application of isolated measures but rather applies a holistic, all-encompassing approach. 3. When women enter the labor force, numerous positive outcomes emerge, including increased GDP, educational gains, and decreased maternal mortality. Despite these benefits, women's employment rates and equality vary significantly worldwide as does support for women in the workforce. This paper will explore the multifaceted benefits of women's employment, the factors influencing labor force participation rates, and the urgency to achieve gender equality as outlined in the 2015 United Nations Sustainable Development Goals (SDGs). It will then examine the emerging presence of women in the traditionally male-dominated nuclear field, specifically within the United Arab Emirates (UAE) as a testament to their resilience and determination to break social norms and advance gender equality.

Zineddin, Dr. Z.↗
DOE OSTI2025

Electromagnetic Analysis of ITER Electron Cyclotron Emission Components and Support Structures for the Preliminary Design Review

The ITER electron cyclotron emission (ECE) diagnostic system is located at Design Shield Module (DSM) 2, Equatorial Port (EP) 9 to measure electron temperature profile and electron temperature fluctuations and also to assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two couples of mirrors for two different optical views. This ECE diagnostic system shall be exposed to significant power due to unabsorbed electron cyclotron heating power in the plasma. It shall also receive large electromagnetic (EM) loads up to 100 MN/m 3 force density due to the eddy currents generated by the short 16-ms transient plasma disruption. The global EM models of EP12 and EP11 with the worst case of plasma disruption Major Disruption Downward Exponential 16 ms Category III (MD_DW_EXP16MS_CATIII) have been provided by ITER. However, the support structure for these critical components, including four mirrors, is still under the design and development phase, especially during the period of preliminary design reviews (PDRs). For those small in-vessel components that are not modeled explicitly for the PDRs, we can extract the magnetic field (B) data and flux variations (dB/dt) at those locations from the global model and then use the B and dB/dt method to calculate the eddy currents, EM forces, and torques on those components in the local submodels. In conclusion, this paper will present the detailed PDR EM analysis results of the ECE components and support structures.

Electromagnetic↗
DOE OSTI2022

Ponderomotive force driven density modifications parallel to B on the LAPD

The ponderomotive force has previously been identified as a possible driver of observed density modifications close to radio frequency actuators during operation. This nonlinear force redistributes density in regions of gradients in the magnitude of an oscillating electric field and describes the influence of the fast time scale RF wave dynamics on slow time scale plasma transport. Depletion of the saturation current (a proxy for the density) measured at the Large Plasma Device (LAPD) was 30–35% during ion cyclotron range of frequencies operation. A coupled 1D plasma transport and cold plasma frequency domain wave solver was developed to self-consistently describe ponderomotive effects and was used to compare with results obtained from the LAPD experiment. The scaled current density driver for the wave model yielded an RF B field in close agreement with two components of the experimental data. However, the 1D parallel model did not accurately reproduce the amplitude or spatial distribution observed in experimental measurements of By. Within the limitations of the 1D model, initial simulation results showed that the ponderomotive force depleted up to 8% for high power (1 MW) and around 1% for the experimental power of 120 kW. This could suggest that the ponderomotive force is not the main driver of density modification for the LAPD experiments presented in this paper. Higher fidelity tools of at least 2D will be required to give a more realistic description of the RF E fields and the effect of the ponderomotive force on the LAPD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development of a military-specific mesh-type computational phantom library and its application to internal dosimetry and in-field radiological triage screening

Estimates of organ-absorbed and committed doses to individuals exposed to radioactive materials via acute inhalation often rely on internal dose coefficients and detector responses from reference human computational models. To achieve more accurate dose assessments to United States Armed Forces service members exposed in-field, computational models with varying morphometric parameters representative of this population are necessary. The International Commission on Radiological Protection (ICRP) Publication 145 provides detailed mesh reference computational phantoms (MRCPs) for adult males and females, with morphometric parameters matched to the 50th percentile. Previously, these phantoms were 2D and 3D scaled to match desired height, mass, and secondary anthropomorphic parameters in the creation of the University of Florida / Memorial Sloan Kettering (UF/MSK) computational phantom library. To achieve body fat percentage targets required for accession into the US Armed Forces, muscle and fat volumes were adjusted accordingly, thus, creating the UF/Department of Defence computational phantom library presented in this study. A comprehensive library of mesh-type computational human phantoms was created, including 57 adult males and 49 adult females with morphometric parameters aligned with United States Armed Forces service members. Phantoms were restricted to a body mass index between 19 and 27.5, with body fat percentages below 26% for males and 36% for females. Specific absorbed fractions were computed for selected source and target combinations, demonstrating how variations in height and body mass influence energy absorption in target regions relative to the ICRP MRCPs. Radiation detector responses were also computed, revealing that higher body masses resulted in decreased registered counts in the detection volume. These findings highlight the importance of incorporating morphometric variability in computational phantoms to achieve more accurate dose assessments and radiation detection responses for United States Armed Forces service members who inhale radioactive materials in-field.

computational phantoms↗
DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

Database and deep-learning scalability of anharmonic phonon properties by automated brute-force first-principles calculations

Understanding the anharmonic phonon properties of crystal compounds—such as phonon lifetimes and thermal conductivities—is essential for investigating and optimizing their thermal transport behaviors. These properties also impact optical, electronic, and magnetic characteristics through interactions between phonons and other quasiparticles and fields. In this study, we develop an automated first-principles workflow to calculate anharmonic phonon properties and build a comprehensive database encompassing more than 6500 inorganic compounds. Utilizing this dataset, we train a graph neural network model to predict thermal conductivity values and spectra from structural parameters, demonstrating a scaling law in which prediction accuracy improves with increasing training data size. High-throughput screening with the model enables the identification of materials exhibiting extreme thermal conductivities—both high and low. The resulting database offers valuable insights into the anharmonic behavior of phonons, thereby accelerating the design and development of advanced functional materials.

Ohnishi, Masato [University of Tokyo (Japan); Inst↗
DOE OSTI2023

Toward prediction of turbulent atmospheric flows over propagating oceanic waves via machine-learning augmented large-eddy simulation

Wind-wave interactions have important effects on the energy harvesting of offshore wind farms. High-fidelity large-eddy simulation (LES) is a powerful approach for investigating wind-wave interactions in turbulent oceanic environments. Due to the large scale of the flow domain and the high grid resolution required to resolve multi-scale flow motions, however, brute-force LES of wind-wave interactions is computationally very expensive. We propose augmenting brute-force LES via machine-learning data-driven modeling (ML-LES) to dramatically reduce the computational time required to obtain converged turbulence statistics when the brute-force approach is employed. Namely, we employ a convolutional neural network (CNN) autoencoder trained and validated with LES data sets to develop a highly efficient ML-LES approach for computing turbulence statistics from just a few snapshots of instantaneous LES flow fields. Further, our results demonstrate the accuracy and efficiency of ML-LES in predicting the mean velocity, velocity fluctuations, and turbulence kinetic energy in highly stretched computational grid systems required to carry out simulations in real-life oceanic environments.

42 ENGINEERING↗
DOE OSTI2020

Seeing is believing: visualization of antiferromagnetic domains

Abstract Understanding and utilizing novel antiferromagnetic (AFM) materials has been recently one of the central issues in condensed matter physics, as well as in materials science and engineering. The relevant contemporary topics include multiferroicity, topological magnetism and AFM spintronics. The ability to image magnetic domains in AFM materials is of key importance for the success of these exciting fields. While imaging techniques of magnetic domains on the surfaces of ferro-(ferri)magnetic materials with, for example, magneto-optical Kerr microscopy and magnetic force microscopy have been available for a number of decades, AFM domain imaging is a relatively new development. We review various experimental techniques utilizing scanning, optical, and synchrotron X-ray probes to visualize AFM domains and domain walls, and to unveil their physical properties. We also discuss the existing challenges and opportunities in these techniques, especially with further increase of spatial and temporal resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Relativistic magnetized jets, such as those from AGN, GRBs, and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and nonthermal radiation. Here, we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron–positron pair plasma. Generalizing recent treatments by Alves et al. and Davelaar et al., we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases, a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Interactions of particles with “continuous spin” fields

Powerful general arguments allow only a few families of long-range interactions, exemplified by gauge field theories of electromagnetism and gravity. However, all of these arguments presuppose that massless fields have zero spin scale (Casimir invariant) and hence exactly boost invariant helicity. This misses the most general behavior compatible with Lorentz symmetry. We present a Lagrangian formalism describing interactions of matter particles with bosonic “continuous spin” fields with arbitrary spin scale ρ. Remarkably, physical observables are well approximated by familiar theories at frequencies larger than ρ, with calculable deviations at low frequencies and long distances. For example, we predict specific ρ-dependent modifications to the Lorentz force law and the Larmor formula, which lay the foundation for experimental tests of the photon’s spin scale. We also reproduce known soft radiation emission amplitudes for nonzero ρ. The particles’ effective matter currents are not fully localized to their worldlines when ρ ≠ 0, which motivates investigation of manifestly local completions of our theory. Our results also motivate the development of continuous spin analogues of gravity and non-Abelian gauge theories. Given the correspondence with familiar gauge theory in the small ρ limit, we conjecture that continuous spin particles may in fact mediate known long-range forces, with testable consequences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Comparison of three DWM-based wake models at above-rated wind speeds

In this study we investigate three mid-fidelity wind turbine wake models based on the dynamic wake meandering (DWM) model principle, and compare their performance with a reference dataset, produced with large-eddy simulations using the actuator line model. The models are compared with respect to flow field, power, and loads on a row of four 5MW reference turbines experiencing above-rated wind conditions. In general, the DWM models show fairly good agreement with large-eddy simulation for the time-averaged flow fields, blade forces and power, with increasing differences along the turbine row. Also when comparing fatigue loads of blade root moments, the differences between the models increase further into the row, with deviations up to 25 % of the reference case. However, while the development in blade root moment fatigue along the turbine row is predominantly driven by the energy content at the frequency corresponding to the turbine's rotational period (1 P ) for the DWM models, the large-eddy simulation results suggest that the key drivers for the blade root and tower loads are the increase in meandering and energy at higher frequencies (> 1 P ) deeper into the turbine row. For the tower loads, the DWM models highly underestimate the fatigue for the waked turbines. From these results, we suggest priorities for future model developments so that robust model implementations can be used in wind farm design and operation.

17 WIND ENERGY↗
DOE OSTI2021

Control of polarization in bulk ferroelectrics by mechanical dislocation imprint

Defects are essential to engineering the properties of functional materials ranging from semiconductors and superconductors to ferroics. Whereas point defects have been widely exploited, dislocations are commonly viewed as problematic for functional materials and not as a microstructural tool. In this work, we developed a method for mechanically imprinting dislocation networks that favorably skew the domain structure in bulk ferroelectrics and thereby tame the large switching polarization and make it available for functional harvesting. The resulting microstructure yields a strong mechanical restoring force to revert electric field–induced domain wall displacement on the macroscopic level and high pinning force on the local level. This induces a giant increase of the dielectric and electromechanical response at intermediate electric fields in barium titanate [electric field–dependent permittivity (ε 33 ) ≈ 5800 and large-signal piezoelectric coefficient ( d 33 *) ≈ 1890 picometers/volt]. Dislocation-based anisotropy delivers a different suite of tools with which to tailor functional materials.

36 MATERIALS SCIENCE↗