Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Search for efficient formulations for Hamiltonian simulation of non-Abelian lattice gauge theories

Hamiltonian formulation of lattice gauge theories (LGTs) is the most natural framework for the purpose of quantum simulation, an area of research that is growing with advances in quantum-computing algorithms and hardware. It, therefore, remains an important task to identify the most accurate, while computationally economic, Hamiltonian formulation(s) in such theories, considering the necessary truncation imposed on the Hilbert space of gauge bosons with any finite computing resources. This paper is a first step toward addressing this question in the case of non-Abelian LGTs, which further require the imposition of non-Abelian Gauss’s laws on the Hilbert space, introducing additional computational complexity. Focusing on the case of SU(2) LGT in 1+1 dimensions coupled to one flavor of fermionic matter, a number of different formulations of the original Kogut-Susskind framework are analyzed with regard to the dependence of the dimension of the physical Hilbert space on boundary conditions, system’s size, and the cutoff on the excitations of gauge bosons. The impact of such dependencies on the accuracy of the spectrum and dynamics obtained from a Hamiltonian computation is examined, and the (classical) computational-resource requirements given these considerations are studied. Besides the well-known angular-momentum formulation of the theory, the cases of purely fermionic and purely bosonic formulations (with open boundary conditions), and the Loop-String-Hadron formulation are analyzed, along with a brief discussion of a Quantum-Link-Model formulation of the same theory. Clear advantages are found in working with the Loop-String-Hadron framework which implements non-Abelian Gauss’s laws a priori using a complete set of gauge-invariant operators. Although small lattices are studied in the numerical analysis of this work, and only the simplest algorithms are considered, a range of conclusions will be applicable to larger systems and potentially to higher dimensions. Future studies will extend this investigation to the analysis of resource requirements for quantum simulations of non-Abelian LGT, with the goal of shedding light on the most efficient Hamiltonian formulation of gauge theories of relevance in nature

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Orbital Chern Insulator and Quantum Phase Diagram of a Kagome Electron System with Half-Filled Flat Bands

In this work, we study the quantum phase diagram of electrons on kagome lattice with half-filled lowest flat bands by considering the antiferromagnetic Heisenberg interaction J, and short-range Coulomb interaction V. In the weak J regime, we identify a fully spin-polarized phase. The presence of finite V drives a spontaneous chiral current, which makes the system an orbital Chern insulator by contributing an orbital magnetization. Such an out-of-plane orbital magnetization allows the presence of a Chern insulating phase independent of the spin orientation in contrast to the spin-orbit coupling induced Chern insulator that disappears with inplane ferromagnetism constrained by symmetry. Such a symmetry difference provides a criterion to distinguish the physical origin of topological responses in kagome systems. The orbital Chern insulator is robust against small coupling J. By further increasing J, we find that the ferromagnetic topological phase is suppressed, which first becomes partially polarized and then enters a nonmagnetic phase with spin and charge nematicity. The frustrated flat band allows the spin and Coulomb interaction to play an essential role in determining the quantum phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Hall effect in LaCrGe 3 at pressures up to 24 kbar and a 𝐻−𝑇 phase diagram of the pressure-induced magnetic phase

LaCrGe 3 is a metallic, itinerant ferromagnet with 𝑇 𝐶 ≈ 86K and an easy magnetic direction along its hexagonal 𝑐 axis. Here we present studies of the 𝐻 ∥ 𝑐 Hall effect in LaCrGe 3 at ambient pressure and under pressures up to 24 kbar. At ambient pressure the anomalous Hall conductivity in the ferromagnetic state reaches a large value of ∼520 Ω − 1 cm −1 . The overall data are consistent with the behavior expected for a ferromagnet with magnetic field applied along its easy axis. Measurements under pressure yield an 𝐻 = 0, 𝑃−𝑇 phase diagram consistent with earlier publications with 𝑇 𝐶 being suppressed with pressure and nonferromagnetic phases appearing for 𝑃 >15 kbar. The data suggest that the observed 𝑇 𝐶 suppression is due to the pressure-dependent exchange interactions rather than suppression of the size of the magnetic moment. Field-dependent Hall resistivity in the higher-pressure, nonferromagnetic state is consistent with the net moment being zero for 𝐻=0. In a finite field there is a sharp anomaly comparable to the one associated with a spin-flip metamagnetic transition. The 𝐻−𝑇 phase diagrams corresponding to a pressure-induced magnetic phase are constructed for several pressures. In conclusion, the Hall data suggest the possibility of the pressure-induced nonferromagnetic phase being a combination of multiple small regions of antiferromagnetic phases with different wave vectors.

Bud'ko, Sergey L. [Ames Laboratory, and Iowa State↗
DOE OSTI2025

A Misalignment Tolerant Magnetic Coupler Design and Implementation for Wireless Charging of Unmanned Aerial Vehicles (UAVs)

Wireless power transfer technology is getting more attention these days as a part of expanded battery charging infrastructure and network at a wide range of power levels. Charging Unmanned aerial vehicles (UAVs) is one wireless charging application, and charging their batteries can be challenging because they have short flight times and are frequently charged. For wireless charging of UAVs, the misalignment performance of the magnetic coupler pads is critical because landing on a charging pad at a perfectly aligned position is difficult and can be time-consuming. Different coil shapes and system structures have been proposed in the literature to overcome this challenge. A new design with a high misalignment tolerance is needed; therefore, to fill this gap, a wireless power transfer system has been designed for a 250 W UAV battery charging system with a new magnetic coupler design. A 3D model of the proposed design was created, and finite element simulations were performed. To prove the proposed design, the magnetic coupler prototypes were developed, and experimental results are presented herein. The misalignment tolerance of the proposed system was observed for various rotational positions of the secondary coil, from 0° to 360°. The results demonstrate that DC–DC efficiency varies between 80% and 85% within full 360° rotation of the receiver coil. This result enables the proposed system to address the misalignment problem for wireless charging of UAVs.

UAVs↗
DOE OSTI2020

Estimation, Minimization, and Validation of Commutation Loop Inductance for a 135-kW SiC EV Traction Inverter

With growing interests in low inductance SiC based power module packaging, it is vital to focus on system level design aspects to facilitate easy integration of the modules and reap system level benefits. To effectively utilize the low inductance modules, busbar and interconnects should also be designed with low stray inductances. A wholistic investigation of the flux path and flux cancellations in the module-busbar assembly which as differentially coupled series inductors is thus mandatory for a system level design. This work presents a busbar design which can be adopted to effectively integrate the CREE’s low inductance 1.2 kV/1.7 kV SiC power modules. The paper also proposes a novel measurement technique to measure the inductance of the modulebusbar assembly as a whole rather than deducing it from individual components. The inductance of the overall commutation loop of the inverter which encompasses the SiC power module, interconnects and PCB busbar have been estimated using finite element analysis (FEA). Furthermore, insights gained from FEA provided the guidelines to decide on the placement of the decoupling capacitors in the busbar to minimize the overall commutation loop inductance from 12.8 nH to 7.4 nH which resulted in significant reduction in the device voltage overshoot. The simulation results have been validated through measurements using an impedance analyzer with less than 5% difference between extracted loop inductance from FEA and measurements. The bus bar design study and the measurement technique discussed in this paper, can be easily extended to other power module packages. Finally, the 135 kW inverter has been compared to a similar highpower inverter utilizing a laminated busbar to highlight the performance of the former.

42 ENGINEERING↗
DOE OSTI2023

Sustainable Manufacturing and the Circular Economy

The term "zero waste" can be somewhat misleading, as it implies an idealization of a concept that may never be achieved in real-world practice. According to Waste360 , "many in the solid waste industry think zero waste is an impossible goal. (Zero-waste advocates) believe just the opposite that a zero-waste approach makes the most economical and environmental sense...Most everyone agrees, however, that actually reaching the point where absolutely no waste is created is unlikely" (Leroux 2001). "Zero waste is a goal. It's completely unachievable...It's not about perfection; it's about making better choices" (Kellogg 2018). The immediate value of moving towards a zero-waste society is not immediately clear, but since the world we live in is finite in its resources, at what point does that value become evident? How do we broaden our value analysis to better incorporate the longer-term view? What are strategies that we can easily incorporate now that will significantly increase that longer term value and mitigate the impacts for future generations? This report considers these questions through a series of case studies from different parts of the U.S. manufacturing sector.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

High-Order Multipole and Binary Love Number Universal Relations

Using a data set of approximately 2 million phenomenological equations of state consistent with observational constraints, we construct new equation-of-state-insensitive universal relations that exist between the multipolar tidal deformability parameters of neutron stars, Λ l , for several high-order multipoles (l = 5,6,7,8), and we consider finite-size effects of these high-order multipoles in waveform modeling. We also confirm the existence of a universal relation between the radius of the 1.4M ⊙ NS, R 1.4 and the reduced tidal parameter of the binary, Λ˜, and the chirp mass. We extend this relation to a large number of chirp masses and to the radii of isolated NSs of different mass M, R M . We find that there is an optimal value of M for every M such that the uncertainty in the estimate of R M is minimized when using the relation. We discuss the utility and implications of these relations for the upcoming LIGO O4 run and third-generation detectors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Slip localization behavior at triple junctions in nickel-base superalloys

Incipient slip localization in the vicinity of hundreds of grain boundary triple junctions (TJs) in a lightly deformed nickel-base superalloy IN718 is studied using a combination of three-dimensional (3D) crystal plasticity finite element (CPFE) modeling, high resolution digital image correlation (HR-DIC) and 3D electron back-scatter diffraction tomography (3D EBSD). A 3D reconstruction method enables identification of thousands of TJs and correspondence of any observed slip bands with their originating TJ lines below the specimen surface. Here, we present a large-scale CPFE model of the experimental 3D microstructure composed of high-fidelity representation of the TJ lines and the boundaries and interiors of the parent grains and use it to calculate the local micromechanical response and slip activity of all TJs at the onset of macroscopic yielding. Statistical analysis of the calculated quantities reveal TJs develop larger stress concentration and grain-average re-orientation than grain interiors and grain boundaries, however no substantial differences in cumulative slip were found among these microstructural regions. We find that TJs with observed slip bands generate lower grain-average re-orientation, fewer active slip systems, and more localized slip on a single system than those without. The distinctions in the reorientation and slip activity are stronger in TJs that experience more intense slip.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comments on the Transplanckian Censorship Conjecture

We consider some aspects of the Transplanckian Censorship Conjecture (TCC), which states that for theories of quantum gravity there is a limit on the lifetime of dS or quintessence states not too different than the current Hubble horizon. If one accepts the de Sitter Swampland conjecture, then the former are ruled out. We consider some aspects of tunneling to an isolated ground state in the presence of time-varying fields, in quantum mechanics and quantum field theory in the absence of gravitation, and note that lifetimes are typically enormous; in fact, there is often a finite probability for the system to remain eternally in its original state. With gravity in a universe with superluminal expansion, while the field evolution may be slowed, Planck scale fluctuations would seem likely to grow to superhorizon size long before the universe decays. We argue that the TCC, if it is correct, requires that superluminal expansion occur only for a brief period in the history of the universe, and will be followed by a p = ρ phase.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Assessing parallel path cooling tower performance via artificial neural networks

Real-time monitoring of a research nuclear reactor, a system in which all generated power is dissipated to the environment, can be performed via analysis of the heat rejection from the cooling system. Given an inlet water temperature and flow rate, the reactor power can be well-approximated from the outlet water temperature; however, the instrumentation to measure outlet conditions may not be robust or accurate. If we know how a cooling tower performs from historical data, but cannot measure the outlet temperature, a mathematical representation of the system can be inverted to obtain the outlet water temperature that describes the cooling capacity. Unfortunately, model inversion processes are computationally expensive. To address this, an artificial neural network (ANN) is implemented to assess the performance of a multi-cell cooling tower for a nuclear reactor. This approach leverages the Merkel model to obtain an extensive data set describing performance of the cooling tower cells throughout a wide array of potential operating conditions. The Merkel model is expressed as a function of four parameters: the inlet and outlet water temperatures, inlet air wet bulb temperature, and ratio of liquid-to-gas mass flow rates (L/G), which together provide a non-dimensional number indicative of cooling tower performance, called the Merkel integral. Computing a 4-dimensional data structure that describes finite combinations of the Merkel integral, an inverse model is then generated using an ANN to determine the cell outlet water temperature from the other three model parameters along with the computed Merkel integral. Compared to traditional model inversion methods, the ANN reduces the computational time by approximately 4 orders of magnitude, with effectively no sacrifice to solution accuracy, and could be applied for different cooling towers in the event the performance curve is known. Finally, three use cases of the ANN are then reviewed: (1) determining the cell outlet water temperatures when gas flow at rated conditions (GFRC) is known, (2) performing the prior case without knowledge of the GRFC, and (3) assessing performance differences between the individual tower cells.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Mechanistic verification of empirical UO 2 fuel fracture models

Standard UO 2 fuel pellets used in light-water reactors fracture during irradiation due to the large thermal gradient in the radial direction. Over the decades, numerous researchers have explored fuel cracking from experimental and modeling points of view. To date, there have been both empirical and mechanistic approaches to predict the number of fragments that form in UO 2 . The empirical models only consider maximum power and burnup as inputs. Existing mechanistic approaches for normal operation have not accounted for irradiation effects. Here, this work employs a mechanistic fuel cracking model using the extended finite element method to explore radial crack formation while including a sensitivity analysis that accounts for the randomization of tensile strength within the fuel, the strength randomization criteria (uniform or volume-weighted Weibull), power ramping rates, computational mesh density, maximum power level, and irradiation (burnup) effects. The results indicate that the uncertainty in this mechanistic modeling approach envelopes the predicted values from three different empirical correlations in almost all cases. This means that, for computationally intensive analyses involving UO 2 fragmentation, the empirical correlations can be used. However, since the mechanistic calculations bound those of the empirical correlations, there is confidence in the applicability of the mechanistic approach developed in this work to generate a correlation for fuel types where limited data exists (e.g., doped-UO 2 , U 3 Si 2 ).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Identification Uncertainty in Inverse Material Model Parameter Determination: A Sensitivity‐Based Decision Process for Load Path Selection

This research proposes a sensitivity-based framework for selecting the optimal prescribed loading path for a biaxial cruciform specimen. Optimality here is determined by the direction and magnitude of the prescribed displacement that minimizes the influence of random noise on the material model parameter identification. Using simulated experimental data based on finite element simulation, in this work, we identify the material model parameters of a Ludwik hardening model and plane stress implementation of the Hill-48 yield criterion using finite element model updating (FEMU). Our analysis reveals that the identification (or estimator) uncertainty of model parameters depends on the displacement boundary conditions (i.e., loading sequence) and the ground-truth value of the individual parameters. Optimal experimental design (OED) criteria based on the Fisher information matrix were investigated to mitigate indecision in the choice of optimal load path when the identification uncertainty of different material model parameters optimized at different load paths. The determinant of the Fisher information matrix was chosen here as the more useful metric due to its ability to capture uncertainty of the most influential material model parameters. The proposed framework demonstrates potential for real-time automated load step selection using scalar criteria derived prior to mechanical loading. The framework can be generalized to other geometries, boundary conditions and material models, allowing this procedure to be utilized for different experimental configurations and materials.

Fayad, Samuel S. [University of Illinois at Urbana↗
DOE OSTI2022

Microstructure-Sensitive Uncertainty Quantification for Crystal Plasticity Finite Element Constitutive Models Using Stochastic Collocation Methods

Uncertainty quantification (UQ) plays a major role in verification and validation for computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allows numerical predictions, providing the bridge from microstructure to materials properties and performances. Several constitutive models have been proposed in the last few decades to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method, which is mathematically rigorous and has been widely used in the field of UQ, to quantify the uncertainty of three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curve, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hygro-thermo-mechanical modeling of partially saturated air-entrained concrete containing dissolved salt and exposed to freeze-thaw cycles

Highlights: • A model for cyclic freezing of air-entrained concrete containing salt is proposed. • The freeze-thaw hysteresis, also including the effect of dissolved salt, is modeled. • The diffusive and convective flux of salt is coupled to the freeze-thaw processes. • The freeze-thaw response of air-entrained concrete containing salt is reproduced. • Liquid absorption from reservoirs increases with an increased concentration of salt. In cold regions, understanding the freeze-thaw behavior of air-entrained concrete is important for designing durable structures and assessing the remaining service life of existing structures. This study presents a hygro-thermo-mechanical multiphase model that describes the cyclic freeze-thaw behavior of partially saturated air-entrained concrete containing dissolved salt. An equilibrium and a non-equilibrium approach are adopted to model the ice formation, including the freeze-thaw hysteresis, inside the porous network. The model also considers the diffusive and convective transport of the dissolved salt coupled to the freeze-thaw processes. Two examples are presented to verify and highlight the capabilities of the model. The first example shows that the model is capable of reproducing the experimentally observed mechanical response of specimens containing NaC1-solutions of different concentrations. In the second example, a larger absorption of liquid from an external reservoir is obtained with an increasing salt concentration in the reservoir, which is consistent with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identification of $^{3}$He–$^{3}$H clusters in the $^{6}$Li+$^{89}$Y experiment using particle-$\gamma$ coincidence measurement

The 6 Li+ 89 Y experiment was performed to explore the reaction mechanism induced by a weakly bound nucleus 6 Li and its cluster configuration. Here, the particle-$\gamma$ coincidence method was used to identify the different reaction channels. The $\gamma$-rays coincident with 3 He/ 3 H indicate that the 3 H/ 3 He stripping reaction plays a significant role in the formation of Zr/Nb isotopes. The obtained results support the existence of a 3 He- 3 H cluster in 6 Li. Direct and sequential transfer reactions are adequately discussed, and the FRESCO code is used to perform precise finite-range cyclic redundancy check calculations. In the microscopic calculation, direct cluster transfer is more predominant than sequential transfer in 3 H transfer. However, the direct cluster transfer is of comparable magnitude to the sequential transfer in the 3 He transfer.

CRC calculations↗
DOE OSTI2020

Three-State Majority-vote Model on Scale-Free Networks and the Unitary Relation for Critical Exponents

Abstract We investigate the three-state majority-vote model for opinion dynamics on scale-free and regular networks. In this model, an individual selects an opinion equal to the opinion of the majority of its neighbors with probability 1 − q , and different to it with probability q . The parameter q is called the noise parameter of the model. We build a network of interactions where z neighbors are selected by each added site in the system, a preferential attachment network with degree distribution k − λ , where λ = 3 for a large number of nodes N . In this work, z is called the growth parameter. Using finite-size scaling analysis, we obtain that the critical exponents $$\beta /\bar{\nu }$$ β / ν ¯ and $$\gamma /\bar{\nu }$$ γ / ν ¯ associated with the magnetization and the susceptibility, respectively. Using Monte Carlo simulations, we calculate the critical noise parameter q c as a function of z for the scale-free networks and obtain the phase diagram of the model. We find that the critical exponents add up to unity when using a special volumetric scaling, regardless of the dimension of the network of interactions. We verify this result by obtaining the critical noise and the critical exponents for the two and three-state majority-vote model on cubic lattice networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A phaseless auxiliary-field quantum Monte Carlo perspective on the uniform electron gas at finite temperatures: Issues, observations, and benchmark study

We investigate the viability of the phaseless finite-temperature auxiliary-field quantum Monte Carlo (ph-FT-AFQMC) method for ab initio systems using the uniform electron gas as a model. Through comparisons with exact results and FT coupled cluster theory, we find that ph-FT-AFQMC is sufficiently accurate at high to intermediate electronic densities. We show, both analytically and numerically, that the phaseless constraint at FT is fundamentally different from its zero-temperature counterpart (i.e., ph-ZT-AFQMC), and generally, one should not expect ph-FT-AFQMC to agree with ph-ZT-AFQMC in the low-temperature limit. With an efficient implementation, we are able to compare exchange-correlation energies to the existing results in the thermodynamic limit and find that the existing parameterizations are highly accurate. In particular, we found that ph-FT-AFQMC exchange-correlation energies are in better agreement with a known parameterization than is restricted path-integral MC in the regime of Θ ≤ 0.5 and rs ≤ 2, which highlights the strength of ph-FT-AFQMC.

Lee, Joonho (ORCID:0000000296671081)↗
DOE OSTI2024

Origin of spin-driven ferroelectricity and effect of external pressure on the complex magnetism of the 6⁢ H perovskite Ba 3 ⁢Ho⁢Ru 2⁢ O 9

The compound Ba 3 HoRu 2 O 9 magnetically orders at 50 K (T N1 ), followed by another complex magnetic ordering at 10.2 K (T N2 ). The second magnetic phase transition was characterized by the coexistence of two competing magnetic ground states associated with two different magnetic wave vectors (K 1 = 0.5 0 0 and K 2 = 0.25 0.25 0). The multiferroicity and magnetoelectric coupling were predicted below T N2 in these 4d-based materials. Here, in this work, we have discussed the origin of spin-driven ferroelectricity, which is not known yet, and the nature of magnetoelectric domains. We have investigated the compound through time-of-flight neutron diffraction, synchrotron x-ray diffraction (XRD), ac susceptibility, frequency-dependent complex dielectric spectroscopy, and dc magnetization under external pressure. We have demonstrated that the noncollinear structure involving two different magnetic ions, Ru (4d) and Ho (4f), breaks the spatial inversion symmetry via inverse Dzyaloshinskii-Moriya (DM) interaction through strong 4d-4f magnetic correlation, which shifts the oxygen atoms and results in nonzero polarization. Such an observation of inverse DM interaction from two different magnetic ions which cause ferroelectricity is rarely observed. The stronger spin-orbit coupling of 4d orbital might play a major role in creating DM interaction of noncollinear spins. We have systematically studied the spin and dipolar dynamics, which exhibit intriguing behavior with shorter coherence lengths of second magnetic phase associated with the k 2 wave vector. The results manifest the development of finite-size magnetoelectric domains instead of true long-range ordering, which justifies the experimentally obtained low value of ferroelectric polarization. The lattice parameters and volume show a sharp anomaly at T N2 obtained by analyzing the temperature-dependence XRD, which is consistent with the ferroelectric transition, predicting a noncentrosymmetric space group, $P\bar{6}2c$, for this compound. Furthermore, we have investigated the effect of external pressure on this complex magnetism. The result reveals an enhancement of ordering temperature by the application of external pressure (~1.6 K/GPa). The external pressure might favor stabilizing the magnetic ground state associated with second magnetic phase. Our study shows an unconventional mechanism of spin-driven ferroelectricity involving inverse DM interaction between Ru (4d) and Ho (4f) magnetic ions due to strong 4d-4f cross coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗