Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “crystal structure prediction”

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 811 records · Page 45

Lifing of Engine Components

The successful development of advanced aerospace engines depends greatly on the capabilities of high performance materials and structures. Advanced materials, such as nickel based single crystal alloys, metal foam, advanced copper alloys, and ceramics matrix composites, have been engineered to provide higher engine temperature and stress capabilities. Thermal barrier coatings have been developed to improve component durability and fuel efficiency, by reducing the substrate hot wall metal temperature and protecting against oxidation and blanching. However, these coatings are prone to oxidation and delamination failures. In order to implement the use of these materials in advanced engines, it is necessary to understand and model the evolution of damage of the metal substrate as well as the coating under actual engine conditions. The models and the understanding of material behavior are utilized in the development of a life prediction methodology for hot section components. The research activities were focused on determining the stress and strain fields in an engine environment under combined thermo-mechanical loads to develop life prediction methodologies consistent with the observed damage formation of the coating and the substrates.

Source record↗

Identifying Order and Disorder in Double Four-Membered Rings via Raman Spectroscopy during Crystallization of LTA Zeolite

Fluoride (F – ) has been essential for the synthesis of low-defect siliceous zeolites. It has been hypothesized that F – balances the positive charges from organic structure-directing agents (OSDAs) and stabilizes key building units during zeolite crystallization such as the double four-membered ring (D4R). However, due to the lack of characterization techniques for investigating medium-range structures, including rings and cages formed during zeolite crystallization, the roles of F – in stabilizing building units and maintaining local charge balance during zeolite assembly are not yet fully understood. Here, the crystallization of siliceous Linde type A (LTA) zeolite in the presence of F – was investigated using Raman spectroscopy and periodic density functional theory (DFT) calculations. We have discovered that the F – -filled double four-membered ring (F – /D4R) and the empty D4R exhibit rather distinct Raman features. Both the F – /D4R and empty D4R are formed in the LTA zeolite synthesized in the presence of F – using 1,2-dimethyl-3-(4-methylbenzyl) imidazolium as the OSDA. Observed Raman bands of the F – /D4R and empty D4R, along with predictive DFT calculations on LTA supercells, reveal an ordered distribution of these two D4R units in the final as-made LTA zeolite. The discovery of these distinct Raman signatures of F – /D4R and empty D4R units opens an interesting new window for studying defects in the D4R during zeolite formation. In particular, we have observed variation in Raman intensities of F – /D4R and empty D4R bands during LTA crystallization; periodic DFT calculations indicate that the observed Raman behavior is consistent with empty D4R units containing one or two Si vacancies surrounded by Q 3 Si—consistent also with solid-state nuclear magnetic resonance measurement. These defects appear to heal during further crystallization, leading to the formation of defect-free LTA zeolite crystals. Altogether, our results provide deeper understanding on the roles of F – in charge balancing and stabilizing intact D4R units during zeolite formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring interlayer coupling in the twisted bilayer PtTe 2

We have investigated interlayer interactions in the bilayer PtTe 2 system, which influence the electronic energy bands near the Fermi level. Our diffusion Monte Carlo (DMC) calculations for the high-symmetry bilayer stackings (AA, AB, AC) manifest distinct interlayer binding characteristics among the stacking modes by revealing significantly different interlayer separations depending on the stacking, which is critical to understanding the interlayer coupling of the twisted bilayers consisting of various local stacking arrangements. Furthermore, a comparison between the interlayer separations obtained from DMC and density functional theory (DFT) shows that meta-generalized gradient approximation (GGA)-based van der Waals–DFT results agree with DMC for different layer stackings, including twisted bilayers, but only the ground-state AA stacking matches well with GGA-based DFT predictions. This underscores the importance of accurate exchange-correlation potentials even for capturing the stacking-dependent interlayer binding properties. We further show that the variability in DFT-predicted interlayer separations is responsible for the large discrepancy of band structures in the 21 . 79 ∘ twisted bilayer PtTe 2 , affecting its classification as metallic or insulating. These results demonstrate the importance of obtaining a correct description of stacking-dependent interlayer coupling in modeling delicate bilayer systems at finite twists. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Mechanism and dynamics of fatty acid photodecarboxylase

Photoenzymes are rare biocatalysts driven by absorption of a photon at each catalytic cycle; they inspire development of artificial photoenzymes with valuable activities. Fatty acid photodecarboxylase (FAP) is a natural photoenzyme that has potential applications in the bio-based production of hydrocarbons, yet its mechanism is far from fully understood. RATIONALE To elucidate the mechanism of FAP, we studied the wild-type (WT) enzyme from Chlorella variabilis (CvFAP) and variants with altered active-site residues using a wealth of techniques, including static and time-resolved crystallography and spectroscopy, as well as biochemical and computational approaches. RESULTS A 1.8-Å-resolution CvFAP x-ray crystal structure revealed a dense hydrogen-bonding network positioning the fatty acid carboxyl group in the vicinity of the flavin adenine dinucleotide (FAD) cofactor. Structures solved from free electron laser and low-dose synchrotron x-ray crystal data further highlighted an unusual bent shape of the oxidized flavin chromophore, and showed that the bending angle (14°) did not change upon photon absorption (step 1) or throughout the photocycle. Calculations showed that bending substantially affected the energy levels of the flavin. Structural and spectroscopic analysis of WT and mutant proteins targeting two conserved active-site residues, R451 and C432, demonstrated that both residues were crucial for proper positioning of the substrate and water molecules and for oxidation of the fatty acid carboxylate by 1 FAD* (~300 ps in WT FAP) to form FAD ∙– (step 2). Time-resolved infrared spectroscopy demonstrated that decarboxylation occured quasi-instantaneously upon this forward electron transfer, consistent with barrierless bond cleavage predicted by quantum chemistry calculations and with snapshots obtained by time-resolved crystallography. Transient absorption spectroscopy in H 2 O and D 2 O buffers indicated that back electron transfer from FAD ∙– was coupled to and limited by transfer of an exchangeable proton or hydrogen atom (step 3). Unexpectedly, concomitant with FAD ∙– reoxidation (to a red-shifted form FAD RS ) in 100 ns, most of the CO 2 product was converted, most likely into bicarbonate (as inferred from FTIR spectra of the cryotrapped FAD RS intermediate). Calculations indicated that this catalytic transformation involved an active-site water molecule. Cryo-Fourier transform infrared spectroscopy studies suggested that bicarbonate formation (step 4) was preceded by deprotonation of an arginine residue (step 3). At room temperature, the remaining CO 2 left the protein in 1.5 μs (step 4'). The observation of residual electron density close to C432 in electron density maps derived from time-resolved and cryocrystallography data suggests that this residue may play a role in stabilizing CO 2 and/or bicarbonate. Three routes for alkane formation were identified by quantum chemistry calculations; the one shown in the figure is favored by the ensemble of experimental data. CONCLUSION Finally, we provide a detailed and comprehensive characterization of light-driven hydrocarbon formation by FAP, which uses a remarkably complex mechanism including unique catalytic steps. We anticipate that our results will help to expand the green chemistry toolkit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Liquid-Liquid Transition in an Undercooled Ti-Zr-Ni Liquid

If crystallization can be avoided, liquids enter a metastable (undercooled) state below their equilibrium liquidus temperatures, TI, finally freezing into a glass below a characteristic temperature called the glass transition temperature, T,. In rare cases, the undercooled liquid may undergo a liquid-liquid phase transition (liquid polymorphism) before entering the glassy state. This has been suggested from experimental studies of HzO and Si4. Such phase transitions have been predicted in some stable liquids, i.e. above TI at atmospheric pressure, for Si02 and BeF;, but these have not been verified experimentally. They have been observed in liquids of P7, Sis and C9, but only under high pressure. All of these transitions are driven by an anomalous density change, i.e. change in local structure, with temperature or pressure. In this letter we present the first experimental evidence for a phase transition in a low viscosity liquid that is not driven by an anomalous density change, but by an approach to a constant configuration state. A maximum in the specific heat at constant pressure, similar to what is normally observed near T,, is reported here for undercooled low viscosity liquids of quasicrystal- forming Ti-Zr-Ni alloys. that includes cooperativity, by incorporating a temperature dependent excitation energy fits the data well, signaling a phase transition.

Lee, G. W.↗

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pseudospin-lattice coupling and electric control of the square-lattice iridate Sr2IrO4

Sr 2 IrO 4 is an archetypal spin-orbit-coupled Mott insulator and has been extensively studied in part because of a wide range of predicted states. Limited experimental characterization of these states thus far brings to light the extraordinary susceptibility of the physical properties to the lattice, particularly, the Ir-O-Ir bond angle. Here, we report a microscopic rotation of the IrO 6 octahedra below 50 K measured by single crystal neutron diffraction. Overall, this sharp lattice anomaly provides keys to understanding the anomalous low-temperature physics and a direct confirmation of a crucial role that the Ir-O-Ir bond angle plays in determining the ground state. Indeed, as also demonstrated in this study, applied electric current readily weakens the antiferromagnetic order via the straightening of the Ir-O-Ir bond angle, highlighting that even slight change in the local structure can disproportionately affect the physical properties in the spin-orbit-coupled system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Lead-Free Cs 4 CuSb 2 Cl 12 Layered Double Perovskite Nanocrystals

Concerns about the toxicity of lead-based perovskites have aroused great interest for the development of alternative lead-free perovskite-type materials. Recently, theoretical calculations predict that Pb 2+ cations can be substituted by a combination of Cu 2+ and Sb 3+ cations to form a vacancy-ordered layered double perovskite structure with superior optoelectronic properties. However, accessibilities to this class of perovskite-type materials remain inadequate, hindering their practical implementations in various applications. Here, we report the first colloidal synthesis of Cs 4 CuSb 2 Cl 12 perovskite-type nanocrystals (NCs). The resulting NCs exhibit a layered double perovskite structure with ordered vacancies and a direct bandgap of 1.79 eV. A composition-structure-property relationship has been established by investigating a series of Cs 4 Cu x Ag 2-2x Sb 2 Cl 12 perovskite-type NCs (0 ≤ x ≤ 1). The composition-induced crystal structure transformation, thus the electronic bandgap evolution has been explored by experimental observations and further confirmed by theoretical calculations. Taking advantages of both the unique electronic structure and solution processability, we demonstrate that the Cs 4 CuSb 2 Cl 12 NCs can be solution-processed as high-speed photodetectors with ultrafast photo-response and narrow bandwidth. We anticipate that our study will prompt future research to design and fabricate novel and high-performance lead-free perovskite-type NCs for a range of applications.

25 ENERGY STORAGE↗

Thermophysical and Optical Properties of Semiconducting Ga2Te3 Melt

The majority of bulk semiconductor single crystals are presently grown from their melts. The thermophysical and optical properties of the melts provide a fundamental understanding of the melt structure and can be used to optimize the growth conditions to obtain higher quality crystals. In this paper, we report several thermophysical and optical properties for Ga2Te3 melts, such as electrical conductivity, viscosity, and optical transmission for temperatures ranging from the melting point up to approximately 990 C. The conductivity and viscosity of the melts are determined using the transient torque technique. The optical transmission of the melts is measured between the wavelengths of 300 and 2000 nm by an dual beam reversed-optics spectrophotometer. The measured properties are in good agreement with the published data. The conductivities indicate that the Ga2Te3 melt is semiconductor-like. The anomalous behavior in the measured properties are used as an indication of a structural transformation in the Ga2Te3 melt and discussed in terms of Eyring's and Bachinskii's predicted behaviors for homogeneous melts.

Li, Chao↗

Low-frequency Raman spectrum of 2D layered perovskites: Local atomistic motion or superlattice modes?

We report the low-frequency Raman spectrum (ω = 10 cm -1 –150 cm -1 ) of a wide variety of alkylammonium iodide based 2D lead halide perovskites (2D LHPs) as a function of A-site cation (MA = methylammonium and FA = formamidinium), octahedral layer thickness (n = 2–4), organic spacer chain length (butyl-, pentyl-, hexyl-), and sample temperature (T = 77 K–293 K). Using density functional theory calculations under the harmonic approximation for n = 2 BA:MAPbI, we assign several longitudinal/transverse optical phonon modes between 30 cm -1 and 100 cm -1 , the eigen displacements of which are analogous to that observed previously for octahedral twists/distortions in bulk MAPbI. Additionally, we propose an alternative assignment for low-frequency modes below this band (<30 cm -1 ) as zone-folded longitudinal acoustic phonons corresponding to the periodicity of the entire layered structure. Furthermore, we compare measured spectra to predictions of the Rytov elastic continuum model for zone-folded dispersion in layered structures. Our results are consistent across the various 2D LHPs studied herein, with energetic shifts of optical phonons corresponding to microscopic structural differences between materials and energetic shifts of acoustic phonons according to changes in the periodicity and elastic properties of the perovskite/organic subphases. This study highlights the importance of both the local atomic order and the superlattice structure on the vibrational properties of layered 2D materials.

2D materials↗

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗

Surface Composition of NiPd Alloys

Surface segregation in Ni-Pd alloys has been studied using the BFS method for alloys. Not only does the method predict an oscillatory segregation profile but it also indicates that the number of Pd-enriched surface planes can vary as a function of orientation. The segregation profiles were computed as a function of temperature, crystal face, and composition. Pd enrichment of the first layer is observed in (111) and (100) surfaces, and enrichment of the top two layers occurs for (110) surfaces. In all cases, the segregation profile shows oscillations that are actually related to weak ordering tendencies in the bulk. An atom-by-atom analysis was performed to identify the competing mechanisms leading to the observed surface behaviors. Large-scale atomistic simulations were also performed to investigate the temperature dependence of the segregation profiles as well as for analysis of the bulk structures. Finally, the observed surface behaviors are discussed in relation to the bulk phase structure of Ni-Pd alloys, which exhibit a tendency to weakly order.

Noebe, Ronald D.↗

Anomaly Detection in Materials Digital Twins with Multiscale ICME for Additive Manufacturing

Detecting anomaly in fatigue and fracture experimental materials science is an interesting yet challenging topic. The reasons are threefold. First, the anomalous microstructure feature that gives rise to structural failure is small, sometimes in the order of 10 -7 of the interrogated volume. This, in turn, results in a highly imbalanced classification problem in machine learning (ML). Second, the consequence is high, in the sense that the test specimen is destructed in such case. Third, the convolution between microstructure stochasticity and the small probability of void nucleation, growth, and coalescence makes failure and fracture a hard-to-predict and challenging problem in materials science due to its irreproducibility, even experimentally. In this paper, we developed a materials digital twin and applied anomaly detection methods to detect voids and anomaly in additive manufacturing (AM). The materials digital twin is driven by two integrated computational materials engineering (ICME) models, which are kinetic Monte Carlo (kMC) and crystal plasticity finite element method (CPFEM). In conclusion, we demonstrated that by using anomaly detection, it is possible to detect voids and other defects in materials digital twin, which paves way for future research in integrating materials digital twin with its physical counterpart.

ICME↗

NASA Tech Briefs, May 2007

Topics include: Noise-Canceling Helmet Audio System; Program Analyzes Spacecraft/Ground Radio Links; Two-Way Communication Using RFID Equipment and Techniques; Six-Message Electromechanical Display System; Scanning Terahertz Heterodyne Imaging Systems; Master Clock and Time-Signal-Distribution System; Synchronous Phase-Resolving Flash Range Imaging; Integrated Radial Probe Transition From MMIC to Waveguide; Bar-Code System for a Microbiological Laboratory; MMIC Amplifier Produces Gain of 10 dB at 235 GHz; Mapping Nearby Terrain in 3D by Use of a Grid of Laser Spots; Digital Beam Deflectors Based Partly on Liquid Crystals; Narrow-Band WGM Optical Filters With Tunable FSRs; Better Finite-Element Analysis of Composite Shell Structures; Computing Spacecraft-Pointing Vectors for Limb Tracking; Enhanced Master Controller Unit Tester; Rover Graphical Simulator; Increasing Durability of Flame-Sprayed Strain Gauges; Multifunctional, High-Temperature Nanocomposites; Multilayer Impregnated Fibrous Thermal Insulation Tiles; Radiation-Shielding Polymer/Soil Composites; Film/Adhesive Processing Module for Fiber-Placement Processing of Composites; Fabrication of Submillimeter Axisymmetric Optical Components; Electrochemical Disposal of Hydrazines in Water; Statistical Model of Evaporating Multicomponent Fuel Drops; Resistively Heated SiC Nozzle for Generating Molecular Beams; Compact Packaging of Photonic Millimeter-Wave Receiver; Diffractive Combiner of Single-Mode Pump Laser-Diode Beams; Wide-Band, High-Quantum-Efficiency Photodetector; A Robustly Stabilizing Model Predictive Control Algorithm; Modeling Evaporation of Drops of Different Kerosenes; Development of Vapor-Phase Catalytic Ammonia Removal System; Several Developments in Space Tethers; Design Concept for a Nuclear Reactor-Powered Mars Rover; Formation-Initialization Algorithm for N Spacecraft; and DNSs of Multicomponent Gaseous and Drop-Laden Mixing Layers Achieving Transition to Turbulence.

Source record↗

Computational prediction of ferromagnetic 𝐴⁢𝑇 6 ⁢𝑋 6 kagome compounds

We present a systematic high-throughput density-functional-theory investigation of the structural and magnetic stability of 312 substitutional compounds in the magnetic kagome 𝐴⁢𝑇 6⁢ 𝑋 6 family. Our screening confirms the stability of many previously reported structures and predicts several additional stable candidates. Within collinear spin configurations, we find that Fe-based systems predominantly adopt antiferromagnetic ground states, whereas Mn-based analogs exhibit a more balanced distribution between ferromagnetic and antiferromagnetic order. For compounds exhibiting several nearly degenerate collinear configurations, we analyze the nature of their magnetic ground states, assess the possible emergence of noncollinear order, and discuss the limitations and uncertainties inherent to standard density-functional approaches. Our electronic structure analysis further reveals that predicted ferromagnetic kagome systems display characteristic features of topological metals, with rich magnetic configurations that can be tuned by chemical substitution. Altogether, these ferromagnetic kagome compounds constitute a broad and still largely unexplored materials platform for the emergence of exciting magnetotransport phenomena.

Density functional calculations↗

Atomic-Level Understanding of Surface Reconstruction Based on Li[Ni x Mn y Co 1–x–y ]O 2 Single-Crystal Studies

The stability of cathode particle surfaces that are directly exposed to the electrolyte is one of the most crucial and determining factors for cathode performance at high operating voltages. Theory has predicted a strong dependence of surface stability on chemical compositions as well as surface facets of layered oxides, yet conflicting results on the correlations exist as most experimental studies focus on cycled secondary particles recovered from composite electrodes. Here, we synthesize well-formed Li[Ni x Mn y Co 1–x–y ]O 2 (NMC) single-crystal samples, carefully define pristine surface properties, and then monitor their evolution with cycling. Atomic-resolution scanning transmission electron microscopy (STEM) imaging and electron energy loss spectroscopy (EELS) analysis show the formation of a surface reconstruction layer (SRL) as well as an extended surface reduction layer on pristine, Li-permeable non-(001) surfaces, even before cycling. We reveal a transition region with chemical gradient, in which the layered structure gradually densifies and eventually transforms into the SRL on the top surface. Contrary to these observations, no SRL is observed on pristine, Li-impermeable (001) surfaces, revealing the facet-dependent nature of surface reconstructions during particle synthesis. Upon electrochemical cycling, significant composition-and facet-dependent SRL growth is observed. The driving force and mechanism for surface reconstruction are further discussed. The present study provides insights into the origin as well as the nature of SRLs, highlighting the significance of surface engineering in cathode material optimization.

25 ENERGY STORAGE↗

Substitution-Mediated Calcination of Nickel-Based Cathodes: Decoupling Lithiation and Crystallization

Nickel-based layered cathodes such as LiNiO 2 offer high energy density for lithium-ion batteries, yet improvements in cycling performance and safety are required for practical use–often achieved through manganese and cobalt substitution as in LiNi 0.80 Mn 0.10 Co 0.10 O 2 (NMC811). However, how such substitution impacts calcination, the key process that governs lithiation, structural ordering, crystallization, and ultimately the resulting material properties, remains unclear. Here, in this study, we investigate substitution-mediated calcination dynamics in NMC811 compared to LiNiO 2 using multiscale-correlated in situ spectroscopy and atomistic-to-mesoscale modeling. While both systems progress through the same sequence of intermediates toward the thermodynamically favored layered phase, NMC811 exhibits an earlier onset of layering, concurrent with hydroxide decomposition followed by sluggish crystallization. Modeling reveals that Mn and Co lower the energy barrier for lithium incorporation and ordering but increase the penalty for interlayer gliding, thereby slowing crystal growth at elevated temperatures. This substitution-mediated decoupling of lithiation and crystallization explains the fine-grained microstructure observed in NMC811 versus coarsened particles in LiNiO 2 and establishes a mechanistic framework for predictive microstructure engineering of Ni-based cathodes.

Calcination↗

Manipulation of spin orientation via ferroelectric switching in Fe-doped Bi 2 WO 6 from first principles

Atomic-scale control of spins by electric fields is highly desirable for future technological applications. Magnetically doped Aurivillius-phase oxides present one route to achieve this, with magnetic ions substituted into the ferroelectric structure at dilute concentrations, resulting in spin-charge coupling. However, there has been minimal exploration of the ferroelectric switching pathways in this materials class, limiting predictions of the influence of an electric field on magnetic spins in the structure. Here, we determine the ferroelectric switching pathways of the end member of the Aurivillius phase family, Bi 2 WO 6 , using a combination of group theoretic analysis and density functional theory calculations. We find that in the ground state P2 1 ab phase, a two-step switching pathway via C2 and Cm intermediate phases provides the lowest energy barrier. Considering iron substitutions on the W site in Bi 2 WO 6 , we determine the spin easy axis. Finally, by tracking the change in spin directionality during ferroelectric switching, we find that a 90° switch in the polarization direction leads to a 112° reorientation of the spin easy axis. The low-symmetry crystal-field environment of Bi 2 WO 6 and magnetoelastic coupling on the magnetic dopant provide a route to spin control via an applied electric field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗