Materials Science (cond-mat.mtrl-sci)

Abstract: Second-harmonic generation (SHG) has emerged as a promising tool for detecting electronic and magnetic structures in noncentrosymmetric materials, but 100$\%$ spin-polarized SHG has not been reported. In this work, we demonstrate nearly 100$\%$ spin-polarized SHG from half-metallic ferrimagnet Mn$_{2}$RuGa. A band gap in the spin-down channel suppresses SHG, so the spin-up channel contributes nearly all the signal, as large as 3614 pm/V about 10 times larger than that of GaAs. In the spin-up channel, $\chi_{xyz}^{(2)}$ is dominated by the large intraband current in three highly dispersed bands near the Fermi level. With the spin-orbit coupling (SOC), the reduced magnetic point group allows additional SHG components, where the interband contribution is enhanced. Our finding is important as it predicts a large and complete spin-polarized SHG in a all-optical spin switching ferrimagnet. This opens the door for future applications.

Abstract: Needle-like twins are observed experimentally within the transition layer at the martensite-twinned martensite interface. We utilize a phase-field approach to investigate this microstructure. Our goal is to simulate the morphology of the transition layer and to perform a detailed analysis to characterize its interfacial and elastic micro-strain energy. To illustrate the micromechanical framework developed for that purpose, sample computations are carried out for a CuAlNi shape memory alloy undergoing the cubic-to-orthorhombic martensitic transformation. A particular focus of the study is on size-dependent morphology through examining the impact of twin spacing. Additionally, our results reveal that certain twin volume fractions lead to the emergence of twin branching, as a way to minimize the total free energy stored in the microstructure.

Abstract: Ferroelectric domain walls are a rich source of emergent electronic properties and unusual polar order. Recent studies showed that the configuration of ferroelectric walls can go well beyond the conventional Ising-type structure. N\'eel-, Bloch-, and vortex-like polar patterns have been observed, displaying strong similarities with the spin textures at magnetic domain walls. Here, we report the discovery of antiferroelectric domain walls in the uniaxial ferroelectric Pb$_{5}$Ge$_{3}$O$_{11}$. We resolve highly mobile domain walls with an alternating displacement of Pb atoms, resulting in a cyclic 180$^{\circ}$ flip of dipole direction within the wall. Density functional theory calculations reveal that Pb$_{5}$Ge$_{3}$O$_{11}$ is hyperferroelectric, allowing the system to overcome the depolarization fields that usually suppress antiparallel ordering of dipoles along the longitudinal direction. Interestingly, the antiferroelectric walls observed under the electron beam are energetically more costly than basic head-to-head or tail-to-tail walls. The results suggest a new type of excited domain-wall state, expanding previous studies on ferroelectric domain walls into the realm of antiferroic phenomena.

Abstract: The PAW method has been used to compute the iron L$_{2,3}$ edges of x-ray absorption spectra (XAS) and x-ray magnetic circular dichroism (XMCD) of the spin-crossover Fe(phen)$_{2}$(NCS)$_{2}$ molecule when adsorbed on Cu(001) surface and in the gas phase, for both the high spin (HS) and low spin (LS) states. It is found that the calculated XAS and XMCD with the static core hole or the Slater transition state half hole are in less good agreement with experiment than those using the so called initial state. This disagreement is due to the reduction of the iron spin magnetic moment caused by the static screening of the core hole by the photo-electron. The L$_{2,3}$ XAS formula is found to be directly related to the unoccupied $3d$ density of states (DOS), and hence the symmetry broken $e_g$ and the $t_{2g}$ iron DOS are used to explain the XAS and XMCD results. It is demonstrated that the dependence of the HS XMCD on the direction of incident x-ray circularly polarized light with respect to the magnetization direction can be used to determine the iron octahedron deformation, while the XMCD for various magnetization directions is directly related to the anisotropy of the orbital magnetic moment and the magneto-crystalline energy. It is also shown that the magnetic dipole moment $T_z$ is very large due to the strong distortion of the iron octahedron and is necessary for an accurate determination of the sum rule computed spin magnetic moment.

Abstract: We establish the sign of the linear magnetoelectric (ME) coefficient, $\alpha$, in chromia, Cr$_2$O$_3$. Cr$_2$O$_3$ is the prototypical linear ME material, in which an electric (magnetic) field induces a linearly proportional magnetization (polarization), and a single magnetic domain can be selected by annealing in combined magnetic (H) and electric (E) fields. Opposite antiferromagnetic domains have opposite ME responses, and which antiferromagnetic domain corresponds to which sign of response has previously been unclear. We use density functional theory (DFT) to calculate the magnetic response of a single antiferromagnetic domain of Cr$_2$O$_3$ to an applied in-plane electric field at 0 K. We find that the domain with nearest neighbor magnetic moments oriented away from (towards) each other has a negative (positive) in-plane ME coefficient, $\alpha_{\perp}$, at 0 K. We show that this sign is consistent with all other DFT calculations in the literature that specified the domain orientation, independent of the choice of DFT code or functional, the method used to apply the field, and whether the direct (magnetic field) or inverse (electric field) ME response was calculated. Next, we reanalyze our previously published spherical neutron polarimetry data to determine the antiferromagnetic domain produced by annealing in combined E and H fields oriented along the crystallographic symmetry axis at room temperature. We find that the antiferromagnetic domain with nearest-neighbor magnetic moments oriented away from (towards) each other is produced by annealing in (anti-)parallel E and H fields, corresponding to a positive (negative) axial ME coefficient, $\alpha_{\parallel}$, at room temperature. Since $\alpha_{\perp}$ at 0 K and $\alpha_{\parallel}$ at room temperature are known to be of opposite sign, our computational and experimental results are consistent.

Abstract: Friction and wear of metals are critically influenced by the microstructures of the bodies constituting the tribological contact. Understanding the microstructural evolution taking place over the lifetime of a tribological system therefore is crucial for strategically designing tribological systems with tailored friction and wear properties. Here, we focus on single-crystalline CoCrFeMnNi that is prone to form twins at room temperature. Deformation twins feature a pronounced orientation dependence with a tension-compression anisotropy, a distinct strain release in an extended volume and robust onset stresses. This makes deformation twinning an ideal probe to experimentally investigate the complex stress fields occurring in a tribological contact. Our results clearly show a grain orientation dependence of twinning under tribological load. Unexpectedly, neither the crystal direction parallel to the sliding nor the normal direction are solely decisive for twinning. This experimental approach is ideal to experimentally validate tribological stress field models, as is demonstrates here.

Abstract: Using both time-dependent density functional theory (TDDFT) and the ``single-shot" $GW$ plus Bethe-Salpeter equation ($GW$-BSE) approach, we compute optical band gaps and optical absorption spectra from first principles for eight common binary and ternary closed-shell metal oxides (MgO, Al$_2$O$_3$, CaO, TiO$_2$, Cu$_2$O, ZnO, BaSnO$_3$, and BiVO$_4$), based on the non-empirical Wannier-localized optimally-tuned screened range-separated hybrid functional. Overall, we find excellent agreement between our TDDFT and $GW$-BSE results and experiment, with a mean absolute error less than 0.4 eV, including for Cu$_2$O and ZnO, traditionally considered to be challenging for both methods.

Abstract: Circular dichroism in angle-resolved photoemission (CD-ARPES) is one of the promising techniques for obtaining experimental insight into topological properties of novel materials, in particular to the orbital angular momentum (OAM) in dispersive bands, which might be related, albeit certainly in a non-trivial way, to the momentum resolved Berry curvature of the bands. Therefore, it is important to understand how non-vanishing CD-ARPES signal arises in graphene, a material where Dirac bands are made from C $|2p_z\rangle$ orbitals that carry zero OAM, spin-orbit-coupling (SOC) can be neglected, and Berry curvature effectively vanishes. Dubs et al., Phys. Rev. B 32, 8389 (1985) have demonstrated non-vanishing cricular dichroism in angular distribution (CDAD) from an oriented $p_z$ orbital, and this process can be responsible for the experimentally observed CD-ARPES in graphene. In this paper, we derive the CD-ARPES from $p_z$ orbitals by elementary means, using only simple algebraic formulas and tabulated numerical values, and show that it leads to significant CD-ARPES signal over the entire vacuum ultraviolet and soft x-ray energy range, with an exception of the photon energy region near $h\nu \approx 40$ eV. We also demonstrate that another process, emerging from the finite electron inelastic mean free path, also leads to CD-ARPES of the potentially similar order of magnitude, as previously discussed by Moser, J. Electron Spectrosc. Relat. Phenom. 214, 29 (2017). We present calculated CDAD maps for selected orbitals and briefly discuss the consequences of the findings for CD-ARPES, focusing on graphene, graphite and WSe$_2$.

Abstract: Triphosphides, with a chemical formula of XP$_3$ (X is a group IIIA, IVA, or VA element), have recently attracted much attention due to their great potential in several applications. Here, using density functional theory calculations, we describe for the first time the structural and electronic properties of the bulk bismuth triphosphide (BiP$_3$). Phonon spectra and molecular dynamics simulations confirm that the 3D crystal of BiP$_3$ is a metal thermodynamically stable with no bandgap. Unlike the bulk, the mono-, bi-, tri-, and tetra-layers of BiP$_3$ are semiconductors with a bandgap ranging from 1.4 to 0.06 eV. However, stackings with more than five layers exhibit metallic behavior equal to the bulk. The results show that quantum confinement is a powerful tool for tuning the electronic properties of BiP$_3$ triphosphide, making it suitable for technological applications. Building on this, the electronic properties of van der Waals heterostructure constructed by graphene (G) and the \bip~monolayer (m-\bip) were investigated. Our results show that the Dirac cone in graphene remains intact in this heterostructure. At the equilibrium interlayer distance, the G/m-BiP$_3$ forms an n-type contact with a Schottky barrier height of 0.5 eV. It is worth noting that the SHB in the G/m-BiP$_3$ heterostructure can be adjusted by changing the interlayer distance or applying a transverse electric field. Thus, we show that few-layers \bip~is an interesting material for realizing nanoelectronic and optoelectronic devices and is an excellent option for designing Schottky nanoelectronic devices.

Abstract: Atomic simulations of material microstructure require significant resources to generate, store and analyze. Here, atomic descriptor functions are proposed as a general latent space to compress atomic microstructure, ideal for use in large-scale simulations. Descriptors can regress a broad range of properties, including character-dependent dislocation densities, stress states or radial distribution functions. A vector autoregressive model can generate trajectories over yield points, resample from new initial conditions and forecast trajectory futures. A forecast confidence, essential for practical application, is derived by propagating forecasts through the Mahalanobis outlier distance, providing a powerful tool to assess coarse-grained models. Application to nanoparticles and yielding of dislocation networks confirms low uncertainty forecasts are accurate and resampling allows for the propagation of smooth microstructure distributions. Yielding is associated with a collapse in the intrinsic dimension of the descriptor manifold, which is discussed in relation to the yield surface.

Abstract: Upon film growth by physical vapor deposition, the preferential orientation of polar organic molecules can result in a non-zero permanent dipole moment (PDM) alignment, causing a macroscopic film polarization. This effect, known as spontaneous orientation polarization (SOP), was studied in the case of different phosphine oxides. We investigate the control of SOP by molecular design and film-growth conditions. Our results show that using less polar phosphine oxides with just one phosphor-oxygen bond yields an exceptionally high degree of SOP with the so-called giant surface potential (slope) reaching more than 150mV/nm in a neat BCPO film grown at room temperature. Additionally, by altering the evaporation rate and the substrate temperature, we are able to control the SOP magnitude over a broad range from 0 to almost 300mV/nm. Diluting BCPO in a non-polar host enhances the PDM alignment only marginally, but combining temperature control together with dipolar doping can result in almost perfectly aligned molecules with more than 80% of their PDMs standing upright on the substrate on average.

Abstract: In this paper, we report the growth of pure {\alpha}-MoO3 micro-flakes by CVD technique and their structural, electronic, optical, and magnetic properties. Samples are annealed at various temperatures in an H2 atmosphere to induce ferromagnetism. All the samples exhibit ferromagnetism at room temperature, and 250oC annealed sample shows the highest magnetic moment of 0.087 emu/g. It is evident from PL data that pristine as well as annealed samples contain different types of defects like oxygen vacancies, surface defects, interstitial oxygen, etc. It is deduced from the analysis of Mo3d and O1s core-level XPS spectra that oxygen vacancies increase up to an annealing temperature of 250oC that correlates with the magnetic moment. Significant changes in the total density of states and also in the magnetic moment for two and three oxygen vacancies are noticed through first-principle-based calculations. It is concluded that the magnetic moment is produced by oxygen vacancies or vacancy clusters, which is consistent with our experimental findings.

Abstract: Controlled doping and understanding its underlying microscopic mechanisms is crucial for advancement of nanoscale electronic technologies, especially in semiconducting single-wall carbon nanotubes (s-SWNTs), where adsorbed counterions are known to govern redox-doping levels. However, modeling the associated 'Coulomb defects' is challenging due to the need for large-scale simulations at low doping levels. Using modified H\"uckel calculations on 120 nm long s-SWNTs with adsorbed $\rm Cl^-$ ions, we study the scaling properties of shallow Coulomb defect states at the valence band edge and quantum well (QW) states in the conduction band. Interestingly, the QW states may underlie observed exciton band shifts of inhomogeneously doped semiconductors. Binding energies of Coulomb defects are found to scale with counterion distance, effective band mass, relative permittivity and counterion charge according to $d^{\alpha-2}m^{\alpha-1}\epsilon_r^{-\alpha}|z_j|^{\alpha}$, with $\alpha$ as an empirical parameter, deepening our understanding of s-SWNT doping.

Abstract: We present a simple Landau phenomenology for plastic-to-crystal phase transitions and use the resulting model to calculate barocaloric effects in plastic crystals that are driven by hydrostatic pressure. The essential ingredients of the model are (i) a multipole-moment order parameter that describes the orientational ordering of the constituent molecules, (ii) coupling between such order parameter and elastic strains, and (iii) the thermal expansion of the solid. The model captures main features of plastic-to-crystal phase transitions, namely large volume and entropy changes at the transition, and strong dependence of the transition temperature with pressure. Using solid C$_{60}$ under $0.60\,$GPa as case example, we show that calculated peak isothermal entropy changes of $\sim 58 \,{\rm J K^{-1} kg^{-1}}$ and peak adiabatic entropy changes of $\sim 23 \,{\rm K}$ agree well with experimental values.

Abstract: An essential aspect of extending safe operation of the active nuclear reactors is understanding and predicting the embrittlement that occurs in the steels that make up the Reactor pressure vessel (RPV). In this work we integrate state of the art machine learning methods using ensembles of neural networks with unprecedented data collection and integration to develop a new model for RPV steel embrittlement. The new model has multiple improvements over previous machine learning and hand-tuned efforts, including greater accuracy (e.g., at high-fluence relevant for extending the life of present reactors), wider domain of applicability (e.g., including a wide-range of compositions), uncertainty quantification, and online accessibility for easy use by the community. These improvements provide a model with significant new capabilities, including the ability to easily and accurately explore compositions, flux, and fluence effects on RPV steel embrittlement for the first time. Furthermore, our detailed comparisons show our approach improves on the leading American Society for Testing and Materials (ASTM) E900-15 standard model for RPV embrittlement on every metric we assessed, demonstrating the efficacy of machine learning approaches for this type of highly demanding materials property prediction.

Abstract: RhPb was initially recognized as one of a CoSn-like compounds with $P6/mmm$ symmetry, containing an ideal kagome lattice of $d$-block atoms. However, theoretical calculations predict the realization of the phonon soft mode which leads to the kagome lattice distortion and stabilization of the structure with $P\bar{6}2m$ symmetry [A. Ptok et al., Phys. Rev. B 104, 054305 (2021)]. Here, we present the single crystal x-ray diffraction results supporting this prediction. Furthermore, we discuss the main dynamical properties of RhPb with $P\bar{6}2m$ symmetry. The bulk phononic dispersion curves contain several flattened bands, Dirac nodal lines, and triple degenerate Dirac points. As a consequence, the phononic drumhead surface state is realized for the (100) surface, terminated by the zigzag-like edge of Pb honeycomb sublattice.

Abstract: Ruthenium oxide with the rutile structure is one of example of altermagnets. These systems are characterized by compensated magnetic moments (typical for antiferromagnets) and strong time reversal symmetry breaking (typical for ferromagnets). However, in such cases, the electronic band structure exhibit strong spin splitting along some directions in the momentum space. Occurrence of the compensated magnetic textures allows for realization of surfaces with specific magnetization, which dependent on the surface orientation and/or its termination. Here, we study interplay between the electronic surface states and the surface magnetization. We show that the spin-resolved spectra strongly depends on a direction in reciprocal space. Such properties can be used for the experimental confirmation of the altermagnetism in RuO$_{2}$ within the spectroscopic techniques. Additionally, we show that the most modified orbitals in the system are $d_{z^{2}}$ and $d_{xy}$ orbitals of Ru. Similarly, the Ru $e_{g}$ states are most sensitive on epitaxial strain, what can suggest some link between altermagnetism and strain.