Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orientational order”

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 433 records · Page 24

Numerical Study of a Passive Wedge Absorber System for Momentum Selection of Muon Beams

The Muon g-2 experiment aims to measure the anomalous magnetic moment of the muon with higher precision than previous experiments. In order to do this muons are created from protons hitting the target and are then sent to the M4 and M5 lines to be injected into the storage ring. The storage ring only accepts muons with momentum within 0.2% of 3094 MeV/c. In order to decrease the number of lost muons a wedge was placed in the M5 beam line to utilize the properties of ionization cooling to create a more uniform momentum spread. Previous simulations predicted a much larger gain of storable muons than the experimental results showed. In order to determine the cause of this discrepancy, simulations were run using a GEANT4 based code to track muons through the M4 and M5 lines, the inflector, and the storage ring while adjusting the orientation, material, geometrical parameters and studying the performance of this wedge for a wide range of incoming muon energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

4D-STEM of Beam-Sensitive Materials

We report that a scanning electron nanobeam diffraction, or 4D-STEM (four-dimensional scanning transmission electron microscopy), is a flexible and powerful approach to elucidate structure from "soft" materials that are challenging to image in the transmission electron microscope because their structure is easily damaged by the electron beam. In a 4D-STEM experiment, a converged electron beam is scanned across the sample, and a pixelated camera records a diffraction pattern at each scan position. This four-dimensional data set can be mined for various analyses, producing maps of local crystal orientation, structural distortions, crystallinity, or different structural classes. Holding the sample at cryogenic temperatures minimizes diffusion of radicals and the resulting damage and disorder caused by the electron beam. The total fluence of incident electrons can easily be controlled during 4D-STEM experiments by careful use of the beam blanker, steering of the localized electron dose, and by minimizing the fluence in the convergent beam thus minimizing beam damage. This technique can be applied to both organic and inorganic materials that are known to be beam-sensitive; they can be highly crystalline, semicrystalline, mixed phase, or amorphous. One common example is the case for many organic materials that have a π-π stacking of polymer chains or rings on the order of 3.4-4.2 Å separation. If these chains or rings are aligned in some regions, they will produce distinct diffraction spots (as would other crystalline spacings in this range), though they may be weak or diffuse for disordered or weakly scattering materials. We can reconstruct the orientation of the π-π stacking, the degree of π-π stacking in the sample, and the domain size of the aligned regions. This Account summarizes illumination conditions and experimental parameters for 4D-STEM experiments with the goal of producing images of structural features for materials that are beam-sensitive. We will discuss experimental parameters including sample cooling, probe size and shape, fluence, and cameras. 4D-STEM has been applied to a variety of materials, not only as an advanced technique for model systems, but as a technique for the beginning microscopist to answer materials science questions. It is noteworthy that the experimental data acquisition does not require an aberration-corrected TEM but can be produced on a variety of instruments with the right attention to experimental parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multimodal Approach Reveals the Symmetry-Breaking Pathway to the Broken Helix in EuIn 2 ⁢As 2

Understanding and manipulating emergent phases, which are themes at the forefront of quantum-materials research, rely on identifying their underlying symmetries. This general principle has been particularly prominent in materials with coupled electronic and magnetic degrees of freedom, in which magnetic order influences the electronic band structure and can lead to exotic topological effects. However, identifying symmetry of a magnetically ordered phase can pose a challenge, particularly in the presence of small domains. Here we introduce a multimodal approach for determining magnetic structures, which combines symmetry-sensitive optical probes, scattering, and group-theoretical analysis. We apply it to EuIn 2 ⁢As 2 , a material that has received attention as a candidate axion insulator. While first-principles calculations predict this state on the assumption of a simple collinear antiferromagnetic structure, subsequent neutron-scattering measurements reveal a much more intricate magnetic ground state characterized by two coexisting magnetic wave vectors reached by successive thermal phase transitions. The proposed high- and low-temperature phases are a spin helix and a state with interpenetrating helical and Néel antiferromagnetic order termed a “broken helix,” respectively. Employing a multimodal approach, we identify the magnetic structure associated with these two phases of EuIn 2 ⁢As 2 . We find that the higher-temperature phase is characterized by a variation of the magnetic moment amplitude from layer to layer, with the moment vanishing entirely in every third Eu layer. The lower-temperature structure is similar to the broken helix, with one important difference: Because of local strain, the relative orientation of the magnetic structure and the lattice is not fixed. Consequently, the symmetry required to protect the axion phase is not generically protected in EuIn 2 ⁢As 2 , but we show that it can be restored if the magnetic structure is tuned with uniaxial strain. Finally, we present a spin Hamiltonian that identifies the spin interactions that account for the complex magnetic order in EuIn 2 ⁢As 2 . Our work highlights the importance of a multimodal approach in determining the symmetry of complex order parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Seeing structural evolution of organic molecular nano-crystallites using 4D scanning confocal electron diffraction (4D-SCED)

Abstract Direct observation of organic molecular nanocrystals and their evolution using electron microscopy is extremely challenging, due to their radiation sensitivity and complex structure. Here, we introduce 4D-scanning confocal electron diffraction (4D-SCED), which enables direct in situ observation of bulk heterojunction (BHJ) thin films. 4D-SCED combines confocal electron optic setup with a pixelated detector to record focused spot-like diffraction patterns with high angular resolution, using an order of magnitude lower dose than previous methods. We apply it to study an active layer in organic solar cells, namely DRCN5T:PC 71 BM BHJ thin films. Structural details of DRCN5T nano-crystallites oriented both in- and out-of-plane are imaged at ~5 nm resolution and dose budget of ~5 e − /Å 2 . We use in situ annealing to observe the growth of the donor crystals, evolution of the crystal orientation, and progressive enrichment of PC 71 BM at interfaces. This highly dose-efficient method opens more possibilities for studying beam sensitive soft materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetic properties of the Ising-like rare earth pyrosilicate: D-Er 2 Si 2 O 7

Ising-like spin-1/2 magnetic materials are of interest for their ready connection to theory, particularly in the context of quantum critical behavior. In this work we report detailed studies of the magnetic properties of a member of the rare earth pyrosilicate family, D-Er 2 Si 2 O 7 , which is known to display a highly anisotropic Ising-like g-tensor and effective spin-1/2 magnetic moments. We used powder neutron diffraction, powder inelastic neutron spectroscopy (INS), and single crystal AC susceptibility to characterize its magnetic properties. Neutron diffraction enabled us to determine the magnetic structure below the known transition temperature (T N = 1.9 K) in zero field, confirming that the magnetic state is a four-sublattice antiferromagnetic structure with two non-collinear Ising axes, as was previously hypothesized. Our powder INS data revealed a gapped excitation at zero field, consistent with anisotropic (possibly Ising) exchange. An applied field of 1 T produces a mode softening, which is consistent with a field-induced second order phase transition. To assess the relevance of D-Er 2 Si 2 O 7 to the transverse field Ising model, we performed AC susceptibility measurements on a single crystal with the magnetic field oriented in the direction transverse to the Ising axes. This revealed a transition at 2.65 T at 0.1 K, a field significantly higher than the mode-softening field observed by powder INS, showing that the field-induced phase transitions are highly field-direction dependent as expected. Furthermore, these measurements suggest that D-Er 2 Si 2 O 7 may be a candidate for further exploration related to the transverse field Ising model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NSUF RTE completion report for 23-4780: Microstructural Defect Induced Thermal Conductivity Reduction in Uranium Nitride and Thorium Nitride

Uranium nitride (UN) is known to have a higher thermal conductivity than traditional oxide fuels, which could lead to a more efficient energy transport and lower local temperature during its lifetime. But the thermal transport performance of UN in extreme environments has not been systematic studied. This study investigated the irradiation induced microstructural defects in UN and the impacts on thermal conductivity. The samples were produced by spark plasma sintering at University of Texas-San Antonio (UTSA), Los Alamos National Laboratory (LANL), and Idaho National Laboratory (INL). Thermal conductivity of UN before and after 2MeV proton irradiation were measured by using laser metrology at INL in a temperature range of 77-295K. Thermal conductivity measurements in this wide, cryogenic temperature range is critical to understand the phonon scattering mechanisms between the thermal energy carrier, phonons, and different types of irradiation-induced defects. In order to ensure the measurements were conducted in the same grains and minimize the impact of the local heterogeneities, the measurement locations were highlighted by fiducial marks using the focused ion beam (FIB) with grain orientation identified using the electron backscatter diffraction (EBSD) at the Center of Advanced Energy Studies (CAES). Irradiation experiment was conducted at the Ion Beam Laboratory at Texas A&M University (TAMU). A total of 6 samples were irradiated with different irradiation doses and temperatures. After irradiation, the microstructure was characterized by using Transmission Electron Microscope (TEM) at CAES (also in the FIB marks).

36 - MATERIALS SCIENCE↗

CondiDiag1.0: a flexible online diagnostic tool for conditional sampling and budget analysis in the E3SM atmosphere model (EAM)

Abstract. Numerical models used in weather and climate prediction take into account a comprehensive set of atmospheric processes (i.e., phenomena) such as the resolved and unresolved fluid dynamics, radiative transfer, cloud and aerosol life cycles, and mass or energy exchanges with the Earth's surface. In order to identify model deficiencies and improve predictive skills, it is important to obtain process-level understanding of the interactions between different processes. Conditional sampling and budget analysis are powerful tools for process-oriented model evaluation, but they often require tedious ad hoc coding and large amounts of instantaneous model output, resulting in inefficient use of human and computing resources. This paper presents an online diagnostic tool that addresses this challenge by monitoring model variables in a generic manner as they evolve within the time integration cycle. The tool is convenient to use. It allows users to select sampling conditions and specify monitored variables at run time. Both the evolving values of the model variables and their increments caused by different atmospheric processes can be monitored and archived. Online calculation of vertical integrals is also supported. Multiple sampling conditions can be monitored in a single simulation in combination with unconditional sampling. The paper explains in detail the design and implementation of the tool in the Energy Exascale Earth System Model (E3SM) version 1. The usage is demonstrated through three examples: a global budget analysis of dust aerosol mass concentration, a composite analysis of sea salt emission and its dependency on surface wind speed, and a conditionally sampled relative humidity budget. The tool is expected to be easily portable to closely related atmospheric models that use the same or similar data structures and time integration methods.

58 GEOSCIENCES↗

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey. II. Anisotropic Large-scale Coherence in Hot Gas, Galaxies, and Dark Matter

Statistics that capture the directional dependence of the baryon distribution in the cosmic web enable unique tests of cosmology and astrophysical feedback. We use constrained oriented stacking of thermal Sunyaev–Zel’dovich (tSZ) maps to measure the anisotropic distribution of hot gas 2.5–40 Mpc away from galaxy clusters embedded in massive filaments and superclusters. The cluster selection and orientation (at a scale of ∼15 Mpc) use Dark Energy Survey (DES) Year 3 data, while expanded tSZ maps from the Atacama Cosmology Telescope Data Release 6 enable a ∼3× more significant measurement of the extended gas compared to the technique’s proof-of-concept. Decomposing stacks into cosine multipoles of order m , we detect a dipole ( m = 1) and quadrupole ( m = 2) at 8σ–10σ, as well as evidence for m = 4 signal at up to 6σ, indicating sensitivity to late-time non-Gaussianity. We compare to Cardinal simulations with spherical gas models pasted onto dark matter halos. The fiducial tSZ data can discriminate between two models that deplete pressure differently in low-mass halos (mimicking astrophysical feedback), preferring higher average pressure in extended structures. However, uncertainty in the amount of cosmic infrared background contamination reduces the constraining power. Additionally, we apply the technique to DES galaxy density and weak lensing to study for the first time their oriented relationships with tSZ. In the tSZ-to-lensing relation, averaged on 7.5 Mpc (transverse) scales, we observe dependence on redshift but not shape or radial distance. Thus, on large scales, the superclustering of gas pressure, galaxies, and total matter is coherent in shape and extent.

79 ASTRONOMY AND ASTROPHYSICS↗

Evidence for Interfacial Octahedral Coupling as a Route to Enhance Magnetoresistance in Perovskite Oxide Superlattices

Engineering octahedral rotations in oxide heterostructures is a promising route for controlling magnetic properties in perovskites, with recent work focusing on magnetic-ordering temperatures and magnetic anisotropies. Here the effects of interfacial octahedral coupling on magnetoresistance are demonstrated in a series of (La 0.7 Sr 0.3 MnO 3 ) n /(LaFeO 3 ) 10 superlattices grown on (001)- and (111)-oriented SrTiO 3 substrates. The different crystallographic orientations allow for the interfacial octahedral connectivity to be tuned, with weaker interfacial coupling present at the (001)-oriented than the (111)-oriented structures as revealed by density functional theory calculations. In n = 14 superlattices, the effect of orientation on the physical properties is minimal with both (001)- and (111)-oriented samples exhibiting similar magnetoresistance. As the fraction of interfacial volume within the LSMO layers is increased by decreasing n, the magnetoresistive behavior of the samples diverges with significantly larger magnetoresistance magnitudes present in the (111)-oriented superlattices. Furthermore, the results are consistent with octahedral coupling playing a greater role in the functional properties at (111)-heterointerfaces and demonstrate a structure-driven approach to tuning interfacial magnetoresistance in complex-oxide heterostructures.

36 MATERIALS SCIENCE↗

Single crystal yield criterion for chromium based on atomistic studies of isolated 1/2[111] screw dislocations

Using the recently developed Bond Order Potential for Cr, we investigate the mobility of isolated 1/2[111] screw dislocations under stress at 0 K. These studies provide the dependence of the critical resolved shear stress (i.e., the Peierls stress) on the orientation of the maximum resolved shear stress plane and on the shear stress perpendicular to the slip direction. Similar studies are made for loading in tension and compression to show that the onset of yielding in Cr single crystals can be described using a combination of the shear stresses parallel and perpendicular to the slip direction. Based on these results, we develop a single crystal yield criterion for Cr that includes all non-Schmid contributions to the Peierls stress. This criterion is used to predict the shape of the yield surface and compare it with the regular Tresca's hexagon arising from the Schmid law. For all directions of uniaxial load, the criterion is used to quantify the character and magnitude of the tension-compression asymmetry.

36 MATERIALS SCIENCE↗

The microstructure and He + ion irradiation behavior of novel low-activation W-Ta-Cr-V refractory high entropy alloy for nuclear applications

Microstructure and nanohardness of a nearly equimolar W-Ta-Cr-V high entropy alloy (HEA), as well as its irradiation response under He + irradiation, were investigated. The single-phase body-centered cubic nanostructured alloy with a 1 µm thick layer was fabricated on a silicon substrate using a magnetron sputtering method. The HEA film has a complex microstructure consisting of micrometric domains that exhibit internal nanostructure controlled by their crystal orientation. The measured nanohardness of the W-Ta-Cr-V alloy is 13 ± 2 GPa, which significantly exceeds the hardness of nanocrystalline tungsten as a result of the high solid-solution strengthening effect. In order to evaluate the irradiation resistance of the HEA film, the material was irradiated with 200 keV He + ions at room temperature, with two different ion fluences: 1 × 10 16 and 5 × 10 16 ions/cm 2 . Using transmission electron microscopy, a high density of extremely fine He bubbles is observed that were uniformly distributed in the matrix. The increase of He + ion fluence increased the density of bubbles, whereas their size remained at a similar level, which indicates that the damage proceeds by the nucleation of additional He bubbles, not by their growth.

36 MATERIALS SCIENCE↗

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗

Electronically Coupled 2D Polymer/MoS 2 Heterostructures

Emergent quantum phenomena in electronically coupled two-dimensional heterostructures are central to next-generation optical, electronic, and quantum information applications. Tailoring electronic band gaps in coupled heterostructures would permit control of such phenomena and is the subject of significant research interest. Two-dimensional polymers (2DPs) offer a compelling route to tailored band structures through the selection of molecular constituents. However, despite the promise of synthetic flexibility and electronic design, fabrication of 2DPs that form electronically coupled 2D heterostructures remains an outstanding challenge. Here, we report the rational design and optimized synthesis of electronically coupled semiconducting 2DP/2D transition metal dichalcogenide van der Waals heterostructures, demonstrate direct exfoliation of the highly crystalline and oriented 2DP films down to a few nanometers, and present the first thickness-dependent study of 2DP/MoS 2 heterostructures. In this work, control over the 2DP layers reveals enhancement of the 2DP photoluminescence by two orders of magnitude in ultrathin sheets and an unexpected thickness-dependent modulation of the ultrafast excited state dynamics in the 2DP/MoS 2 heterostructure. These results provide fundamental insight into the electronic structure of 2DPs and present a route to tune emergent quantum phenomena in 2DP hybrid van der Waals heterostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning-enhanced model-based scenario optimization for DIII-D

Abstract Scenario development in tokamaks is an open area of investigation that can be approached in a variety of different ways. Experimental trial and error has been the traditional method, but this required a massive amount of experimental time and resources. As high fidelity predictive models have become available, offline development and testing of proposed scenarios has become an option to reduce the required experimental resources. The use of predictive models also offers the possibility of using a numerical optimization process to find the controllable inputs that most closely achieve the desired plasma state. However, this type of optimization can require as many as hundreds or thousands of predictive simulation cases to converge to a solution; many of the commonly used high fidelity models have high computational burdens, so it is only reasonable to run a handful of predictive simulations. In order to make use of numerical optimization approaches, a compromise needs to be found between model fidelity and computational burden. This compromise can be achieved using neural networks surrogates of high fidelity models that retain nearly the same level of accuracy as the models they are trained to replicate while reducing the computation time by orders of magnitude. In this work, a model-based numerical optimization tool for scenario development is described. The predictive model used by the optimizer includes neural network surrogate models integrated into the fast Control-Oriented Transport simulation framework. This optimization scheme is able to converge to the optimal values of the controllable inputs that produce the target plasma scenario by running thousands of predictive simulations in under an hour without sacrificing too much prediction accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Dark Energy Spectrographic Instrument (DESI) Guide-Focus Assembly

The specific research of the DESI project is to study in detail the expansion history of the Universe over the past 10 billion years. In order to do this, the project designed, fabricated, tested and commissioned the DESI instrument which has been deployed at the Mayall Telescope at the Kitt Peak National Observatory near Tucson, Arizona. The DESI project was conducted by approximately 30 US and foreign national laboratories. Our CRADA with the The Institut de Física d’Altes Energies (IFAE) in Barcelona (Spain) involved both research planning for eventual observing on the DESI telescope; and in providing hardware jointly designed by LBNL and IFAE; ten guide and focus sensor systems (GFAs) that either provide feedback to the Mayall telescope tracking system or provide focus position information to the DESI corrector hexapod. The key major goals of this CRADA were achieved.For survey planning, Drs. Fernandez, Fosalba, Garcia-Bellido, and Gaztañaga were all important contributors to the working groups to which they were assigned. Their application of redshift techniques, development of simulations to support target selection, exploration of the use of alternative observables to extract cosmological information like AP and void statistics, and the use of PAU data to assist in the optimization of targeting selections have been vital in the early planning for targets during commissioning. The major role in the project involved Drs. Castander, Miquel, Padilla, and Sànchez. They were all directly involved with the development of the Guide-Focus Assembly (GFA) hardware, software and filter validation. In order to accomplish the goals of the DESI Project, the Mayall telescope must track DESI targeted galaxies for up to 15 minutes and the DESI corrector must be repositioned on its hexapod system for different telescope orientations. The guide and focus assemblies (GFA) provide in this information. There are ten GFA modules on the focal plane, six providing tracking information and four providing focus information. A common GFA design supports both functions by a simple switching of its operating mode. IFAE designed, fabricated, tested and delivered twelve GFAs. Once received at LBNL, the GFAs were integrated with the focal plane. This was an iterative process as the GFA’s were delivered as they were made. Our working relationship with the IFAE scientists and engineers was remarkable given the distance and time zone issues. The GFAs have performed extremely well during the commissioning phase of the instrument.

79 ASTRONOMY AND ASTROPHYSICS↗

The Origin of Anomalous Density Behavior of Silica Glass

The anomalous density–temperature relationship of vitreous silica with low hydroxyl content is explained by the formation of medium-range ordering structure in the glass transition process. The ordered medium-range structure has the shape of a “nanoflake” and consists of two layers of SiO 4 tetrahedra, bonded by O atoms located in the middle of the structure. The nanoflakes interact with their surrounding structures through both covalent chemical bonds and van der Waals bonds. In the formation of the van der Waals bonds, the orientation of SiO 4 tetrahedra can change, which results in an increase in distance between the nanoflakes and their surrounding structures. Thus, there is a slight volume enlargement associated with the formation of nanoflakes. Since the nanoflakes’ formation starts at a temperature near 1480 °C, and the population of the nanoflakes grows continuously as temperature decreases until about 950 °C, the bulk volume of silica glass increases in the temperature range from about 1480 °C to 950 °C. Therefore, the density anomaly of silica glass can be explained as a byproduct of forming of medium-range ordering structure in the silica glass transition.

36 MATERIALS SCIENCE↗

Goal-oriented a-posteriori estimation of model error as an aid to parameter estimation

In this work, a Bayesian model calibration framework is presented that utilizes goal-oriented a-posterior error estimates in quantities of interest (QoIs) for classes of high-fidelity models characterized by PDEs. It is shown that for a large class of computational models, it is possible to develop a computationally inexpensive procedure for calibrating parameters of high-fidelity models of physical events when the parameters of low-fidelity (surrogate) models are known with acceptable accuracy. The main ingredients in the proposed model calibration scheme are goal-oriented a-posteriori estimates of error in QoIs computed using a so-called lower fidelity model compared to those of an uncalibrated higher fidelity model. The estimates of error in QoIs are used to define likelihood functions in Bayesian inversion analysis. A standard Bayesian approach is employed to compute the posterior distribution of model parameters of high-fidelity models. As applications, parameters in a quasi-linear second-order elliptic boundary-value problem (BVP) are calibrated using a second-order linear elliptic BVP. In a second application, parameters of a tumor growth model involving nonlinear time-dependent PDEs are calibrated using a lower fidelity linear tumor growth model with known parameter values.

A-posterior estimates↗

Noncollinear magnetic configurations and substrate-mediated interactions in Mn trimers on the GaN (000$\bar{1}$) surface

Collinear and noncollinear calculations based on density functional theory are carried out to elucidate the magnetic ordering of Mn trimers on a GaN ($000\bar{1}$) substrate. These trimers had previously been observed in $3a \times 3a$ surface reconstructions through Mn deposition onto the N-polar face of wurtzite GaN($000\bar{1}$). In this work, we start off by studying the effect of spin orbit coupling for the case of monomers and dimers of Mn atoms on top of a GaN surface. Based on an effective spin Hamiltonian, we estimate the magnetic anisotropy energy (MAE) for those cases and found that it is four orders of magnitude weaker than the exchange magnetic coupling between Mn adatoms. In the Mn trimer case, the magnetic ground state has the Mn spins in-plane with the GaN surface in which the relative spin orientation within each trimer is noncollinear due to the competition between the two antiferromagnetic interactions that affect each Mn spin in the trimer, which leads to the found energy minimum with 120 degree angles between the spins. By exploring the nature and fundamental mechanisms for the magnetic interaction among the Mn trimers, we find that the surface states of the substrate play a key role, involving a Ruderman-Kittel-Kasuya-Yosida (RKKY)-type interaction. We report on an electron-mediated long-distance exchange coupling between localized magnetic moments on a GaN($000\bar{1}$) surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗