Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive model”

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 757 records · Page 42

Paths Forward for Nuclear Energy: Using a Nationwide Post-Stratified Hierarchical Model to Facilitate Matching of New Nuclear Technologies to Receptive Host Communities

This project was a collaboration between the University of Oklahoma (OU) and the University of Michigan (UMich). The overall objective of the project was to address a critical problem facing the siting of nuclear facilities, where there is no mechanism for incorporating public attitudes at a state and local level into decision-making. Utilizing recent advances in data modeling over local spatial scales, the research team created a state and county level map of public attitudes towards nuclear energy. Current estimates for public acceptance of nuclear technologies primarily exists on a national level, and sub-national estimates of public opinion about reactor siting, research facilities, transport routes, or storage and disposal facilities are very scarce. Model results from this project provide a systematic basis for technology developers to evaluate public acceptance across alternative siting options as a critical component for weighing potential benefit-cost tradeoffs for alternative paths to deploy new nuclear facilities. The project utilized data from the largest database of nuclear attitudes in the US, coupled with Census and elections-related data as well as data on social vulnerability and proximity to current nuclear facilities. In addition to modeling public support, the team also developed an interactive map to showcase results, using an online tool called Public and Local Attitudes about Nuclear Energy Technologies (PLANET). The tool allows technology developers and stakeholders to interact with the model results and better understand variation in support for nuclear facility siting on a state and county level. The interactive nature of the tool also allows users to analyze results in a targeted way, by comparing levels of support across different states and/or counties as desired.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Accurate Simulation of Both Sensitivity and Variability for Amazonian Photosynthesis: Is It Too Much to Ask?

Abstract Estimates of Amazon rainforest gross primary productivity (GPP) differ by a factor of 2 across a suite of three statistical and 18 process models. This wide spread contributes uncertainty to predictions of future climate. We compare the mean and variance of GPP from these models to that of GPP at six eddy covariance (EC) towers. Only one model's mean GPP across all sites falls within a 99% confidence interval for EC GPP, and only one model matches EC variance. The strength of model response to climate drivers is related to model ability to match the seasonal pattern of the EC GPP. Models with stronger seasonal swings in GPP have stronger responses to rain, light, and temperature than does EC GPP. The model to data comparison illustrates a trade‐off inherent to deterministic models between accurate simulation of a mean (average) and accurate responsiveness to drivers. The trade‐off exists because all deterministic models simplify processes and lack at least some consequential driver or interaction. If a model's sensitivities to included drivers and their interactions are accurate, then deterministically predicted outcomes have less variability than is realistic. If a GPP model has stronger responses to climate drivers than found in data, model predictions may match the observed variance and seasonal pattern but are likely to overpredict GPP response to climate change. High or realistic variability of model estimates relative to reference data indicate that the model is hypersensitive to one or more drivers.

54 ENVIRONMENTAL SCIENCES↗

Understanding aerosol–cloud interactions using a single-column model for a cold-air outbreak case during the ACTIVATE campaign

Marine boundary layer clouds play a critical role in Earth's energy balance. Their microphysical and radiative properties are highly impacted by ambient aerosols and dynamic forcings. In this study, we evaluate the representation of these clouds and related aerosol–cloud interaction processes in the single-column version of the E3SM climate model (SCM) against field measurements collected during the NASA ACTIVATE campaign over the western North Atlantic, as well as intercompare results with high-resolution process level models. We show that E3SM SCM reproduces the macrophysical properties of post-frontal boundary layer clouds in a cold-air outbreak (CAO) case well. However, it generates fewer but larger cloud droplets compared to aircraft measurements. Further sensitivity tests show that the underestimation of both aerosol number concentration and vertical velocity variance contributes to this bias. Aerosol–cloud interactions are examined by perturbing prescribed aerosol properties in E3SM SCM with fixed dynamics. Higher aerosol number concentration or hygroscopicity leads to more numerous but smaller cloud droplets, resulting in a stronger cooling via shortwave cloud forcing. This apparent Twomey effect is consistent with prior climate model studies. The cloud liquid water path shows a weakly positive relation with cloud droplet number concentration due to precipitation suppression. This weak aerosol effect on cloud macrophysics may be attributed to the dominant impact of strong dynamical forcing associated with the CAO. Our findings indicate that the SCM framework is a key tool to bridge the gap between climate models, process level models, and field observations to facilitate process level understanding.

54 ENVIRONMENTAL SCIENCES↗

Fingerprinting Interactions between Proteins and Ligands for Facilitating Machine Learning in Drug Discovery

Molecular recognition is fundamental in biology, underpinning intricate processes through specific protein–ligand interactions. This understanding is pivotal in drug discovery, yet traditional experimental methods face limitations in exploring the vast chemical space. Computational approaches, notably quantitative structure–activity/property relationship analysis, have gained prominence. Molecular fingerprints encode molecular structures and serve as property profiles, which are essential in drug discovery. While two-dimensional (2D) fingerprints are commonly used, three-dimensional (3D) structural interaction fingerprints offer enhanced structural features specific to target proteins. Machine learning models trained on interaction fingerprints enable precise binding prediction. Recent focus has shifted to structure-based predictive modeling, with machine-learning scoring functions excelling due to feature engineering guided by key interactions. Notably, 3D interaction fingerprints are gaining ground due to their robustness. Various structural interaction fingerprints have been developed and used in drug discovery, each with unique capabilities. This review recapitulates the developed structural interaction fingerprints and provides two case studies to illustrate the power of interaction fingerprint-driven machine learning. The first elucidates structure–activity relationships in β2 adrenoceptor ligands, demonstrating the ability to differentiate agonists and antagonists. The second employs a retrosynthesis-based pre-trained molecular representation to predict protein–ligand dissociation rates, offering insights into binding kinetics. Despite remarkable progress, challenges persist in interpreting complex machine learning models built on 3D fingerprints, emphasizing the need for strategies to make predictions interpretable. Binding site plasticity and induced fit effects pose additional complexities. Interaction fingerprints are promising but require continued research to harness their full potential.

3D structural interaction fingerprints↗

Quantifying structural errors in cloud condensation nuclei activity from reduced representation of aerosol size distributions

Aerosol effects on clouds and radiation are the dominant contribution to uncertainty in radiative forcing relative to the pre-industrial atmosphere. While previous studies have assessed the impact of parametric uncertainty on modeled forcing, structural errors from the numerical representation of particle distributions have not been well quantified. Here we present a framework for quantifying error in aerosol size distributions and cloud condensation nuclei activity, which we apply to the widely used 4-mode version of the Modal Aerosol Module (MAM4). Box model predictions from the MAM4 are evaluated against the Particle Monte Carlo Model for Simulating Aerosol Interactions and Chemistry (PartMC-MOSAIC), a benchmark model that tracks the evolution of individual particles. We show that size distributions simulated by MAM4 diverge from those simulated by PartMC-MOSAIC after only a few hours of aging by condensation and coagulation in polluted conditions, which leads to large errors in modeled cloud condensation nuclei concentrations. We find that differences between MAM4 and PartMC-MOSAIC are largest under polluted conditions, where the size distribution evolves rapidly though aging by condensation of semi-volatile substances and coagulation among particles. These findings suggest that structural error in modeled aerosol properties contributes to the large inter-model variability in aerosol radiative forcing.

Fierce, Laura M.↗

Integrated Experimental and Modeling Study of Geochemical Reactions of Simple Fracturing Fluids with Caney Shale

Interactions between rock minerals and hydraulic fluids directly impact the geochemical and geomechanical properties of shale formations. However, the mechanisms of geochemical reactions in shale unconventional reservoirs remain poorly understood. Here, to investigate the geochemical reactions between shale and hydraulic fracturing fluids, a series of batch reactor experiments were undertaken. Three rock samples with different mineralogical compositions and three fluid samples of different compositions (deionized water, deionized water + 2% potassium chloride (KCl), and deionized water + 0.5% choline chloride (C 5 H 14 ClNO) were used. Experiments were undertaken at reservoir temperature and atmospheric pressure. Elemental compositions of effluents after 1, 3, 7, 14, 28 days were analyzed using Inductively Couple Plasma Mass Spectrometry (ICP-MS). Medical Computed Tomography (CT) scan and X-ray Fluorescence (XRF) were conducted on the entire core run to help upscale results obtained from rock-fluid interaction experiments. Geochemical modeling using a reactive simulator, TOUGHREACT, was undertaken to corroborate experimental results. Results show that lower pH triggered high dissolution rates in the rock samples, especially the carbonate components. As pH increased, the rate of dissolution declined significantly, though for most cases dissolution still continued. Observed dissolved silica concentrations were much higher than the quartz solubility, suggesting that much of the silica originates from more soluble silica polymorphs and possibly desorption from clay mineral exchange sites. Concentration of most elemental species in solution increased but aluminium and magnesium concentrations declined rapidly following initial entry into solution. Geochemical modeling corroborated the conclusions regarding mineral dissolution and precipitation observed from experiments, notably; the dissolution of calcite and oxidation of pyrite in reacted shale samples, the likely presence of silica polymorphs such as opal, chalcedony or amorphous silica in these samples, and the reduction of Al and Mg concentrations in solution by precipitation of secondary aluminosilicate phases. The de-flocculation of clay minerals during reaction implies fines migration after hydraulic fracturing. This is detrimental to reservoir productivity as clay fines are displaces and lodged within the micro and nano-fractures created during fracturing. The immediate consumption of aluminium and magnesium also has implications on blockage of hydrocarbon pathways due to precipitation of new minerals in these locations.

04 OIL SHALES AND TAR SANDS↗

Assembly of JAZ–JAZ and JAZ–NINJA complexes in jasmonate signaling

Jasmonates (JAs) are plant hormones with crucial roles in development and stress resilience. They activate MYC transcription factors by mediating the proteolysis of MYC inhibitors called JAZ proteins. In the absence of JA, JAZ proteins bind and inhibit MYC through the assembly of MYC–JAZ–Novel Interactor of JAZ (NINJA)–TPL repressor complexes. However, JAZ and NINJA are predicted to be largely intrinsically unstructured, which has precluded their experimental structure determination. Through a combination of biochemical, mutational, and biophysical analyses and AlphaFold-derived ColabFold modeling, we characterized JAZ–JAZ and JAZ–NINJA interactions and generated models with detailed, high-confidence domain interfaces. We demonstrate that JAZ, NINJA, and MYC interface domains are dynamic in isolation and become stabilized in a stepwise order upon complex assembly. By contrast, most JAZ and NINJA regions outside of the interfaces remain highly dynamic and cannot be modeled in a single conformation. Our data indicate that the small JAZ Zinc finger expressed in Inflorescence Meristem (ZIM) motif mediates JAZ–JAZ and JAZ–NINJA interactions through separate surfaces, and our data further suggest that NINJA modulates JAZ dimerization. This study advances our understanding of JA signaling by providing insights into the dynamics, interactions, and structure of the JAZ–NINJA core of the JA repressor complex.

59 BASIC BIOLOGICAL SCIENCES↗

CFD modeling of near-wall combustion and unburned methane prediction in natural gas spark ignition engines

Natural gas-powered engines play a critical role in gas drilling, compression, and transmission sectors, but methane (CH 4 ) from engine combustion slip can be significant over their lifespan, contributing to atmospheric pollution and signaling reduced engine efficiency. Here, to address this challenge, computational fluid dynamics (CFD) simulations offer valuable insights into the in-cylinder combustion process, enabling the optimization of combustion strategies and engine designs to minimize unburned CH 4 slip. This study aims to evaluate and improve combustion models for simulating the combustion process and predicting unburned CH 4 concentrations in natural gas spark-ignition (SI) engines, including engines that are part of combined reformer-engine systems. Specifically, the performance of two flamelet-based combustion models—the Extended Coherent Flame Model (ECFM) and the G-equation model—was assessed using experimental engine data collected under varying excess-air ratio (λ) conditions and fuel compositions, including natural gas and syngas blends. In addition, to enhance the predictive capabilities of the G-equation model, a flame-wall interaction (FWI) sub-model was integrated into its framework. The effects of its model parameters, such as quenching and influence distance, on combustion behavior and unburned methane predictions were analyzed in detail. The ECFM tended to predict delayed combustion phasing under diluted mixture conditions, resulting in overprediction of unburned CH 4 concentrations. In contrast, the G-equation model provided reasonable predictions of combustion pressure, while representing higher the CH 4 reduction rate across the operating condition compared to experimental data. Incorporating the FWI sub-model—with the quenching distance calculated based on a pressure-dependent relation (P -0.48 ) and a fixed influence distance of 1.5 mm—further improved the G-equation model’s accuracy in predicting CH 4 reduction rates without compromising its ability to simulate the combustion process.

Combustion model↗

Gauge invariance of light-matter interactions in first-principle tight-binding models

In this work, we study the different ways of introducing light-matter interaction in first-principle tight-binding (TB) models. The standard way of describing optical properties is the velocity gauge, defined by linear coupling to the vector potential. In finite systems a transformation to represent the electromagnetic radiation by the electric field instead is possible, albeit subtleties arise in periodic systems. The resulting dipole gauge is a multi-orbital generalization of the Peierls substitution. In this work we investigate the accuracy for both pathways, with particular emphasis on gauge invariance, for TB models constructed from maximally localized Wannier functions. This approach accurately captures the light-matter interaction close to the Fermi level. Focusing on paradigmatic two-dimensional materials, we construct first-principle models and calculate the response to electromagnetic fields in linear response and for strong excitations. Benchmarks against fully converged first-principle calculations allow for ascertaining the accuracy of the TB models. We find that the dipole gauge provides a more accurate description than the velocity gauge in all cases. The main deficiency of the velocity gauge is an imperfect cancellation of paramagnetic and diamagnetic current. Formulating a corresponding sum rule however provides a way to explicitly enforce this cancellation. This procedure corrects the TB models in the velocity gauge, yielding excellent agreement with dipole gauge and thus restoring gauge invariance.

2-dimensional systems↗

Triaxiality and the nature of low-energy excitations in Ge 76

The deformation properties of the low-lying states in 76 Ge have been investigated following a safe energy Coulomb excitation measurement with the GRETINA tracking array and CHICO2 heavy-ion counter at the ATLAS accelerator facility at Argonne National Laboratory. A comprehensive set of transition and static E2 matrix elements were extracted from the measured differential Coulomb cross-sections, and compared with results of configuration interaction shell-model calculations and computations carried out within the framework of the generalized triaxial rotor model. The remarkable agreement between the calculated and experimental data supports a near-maximum triaxial deformation for the ground state of 76 Ge. Additionally, the degree of softness of the asymmetry in 76 Ge and 76 Se was investigated using rotational invariants generated from configuration interaction shell-model wave functions computed with the jj44b and JUN45 effective interactions. The resulting invariants are shown to be consistent with a stiff triaxial deformation in 76 Ge and a predominantly soft triaxial potential for 76 Se, in agreement with the conclusions of recent works by this collaboration.

59 ≤ A ≤ 89↗

The impact of coupled air–sea interaction on extreme East Asian summer monsoon simulation in CMIP5 models

Abstract In this study, the relationship between the ability to simulate air–sea interactions over the western North Pacific (WNP), and to reproduce the extreme East Asian summer monsoon (EASM), were investigated by comparing the performances of several global climate models (GCMs). High ranked in air–sea interaction simulation (HRA) and low ranked in air–sea interaction simulation (LRA) models were selected, according to their performance in simulating relations between sea surface temperature (SST) and precipitation over the WNP, from the ensemble of models that participated in the third and fifth phases of the Coupled Model Intercomparison Project (CMIP3, CMIP5). Compared with CMIP3 models, CMIP5 models exhibited improved simulations of the distinctive air–sea interaction over the WNP, namely, the strong atmospheric forcing on the ocean. Among CMIP5 models, HRA models, which reproduced intrinsic negative correlations between precipitation and SST over the WNP, could simulate the extreme EASM better than LRA models. In particular, HRA models generated a more realistic spatial distribution of the extreme EASM compared with LRA models. The defects of the LRA models resulted from distorted synoptic fields, including underestimated geopotential height and overestimated low‐level wind over the WNP, inducing unrealistic moisture supply and convection due to the exaggerated SST forcing. In contrast, reasonable air–sea interactions represented in HRA models lead to realistic synoptic fields over the WNP, and proper simulation of the extreme EASM.

Kim, Taehyung↗

Plant-microbe interactions: from genes to ecosystems using Populus as a model system

Plant-microbe symbioses span a continuum from pathogenic to mutualistic with functional consequences for both organisms in the symbiosis. In order to increase sustainable food and fuel production in the future, it is imperative that we harness these symbioses. The tree genus Populus is an excellent model system for studies examining plant-microbe interactions due to the wealth of genomic information available and the molecular tools that have been developed to manipulate Populus-microbe symbioses. In this review, we highlight how Populus can serve as a model system to explore plant-microbe interactions. Specifically, highlighting research linking Populus-microbe interactions from the gene to the ecosystem level. We explore why Populus is an excellent model for perennial plant systems, the molecular underpinnings of Populus-microbe interactions, how host genetics influence microbial community composition, and how microbial communities vary at fine spatial scales and between Populus species. Further, we explore how the patterns of the microbiome may affect ecosystem level functions in managed and natural ecosystems. Understanding and manipulating these interactions in Populus has the potential to improve plant health and impact ecosystem sustainability and processes as Populus trees function as foundational species in many natural ecosystems and are also deployed in managed ecosystems for various agroforestry applications.

59 BASIC BIOLOGICAL SCIENCES↗

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development and Evaluation of Chemistry‐Aerosol‐Climate Model CAM5‐Chem‐MAM7‐MOSAIC: Global Atmospheric Distribution and Radiative Effects of Nitrate Aerosol

Abstract An advanced aerosol treatment, with a focus on semivolatile nitrate formation, is introduced into the Community Atmosphere Model version 5 with interactive chemistry (CAM5‐chem) by coupling the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the 7‐mode Modal Aerosol Module (MAM7). An important feature of MOSAIC is dynamic partitioning of all condensable gases to the different fine and coarse mode aerosols, as governed by mode‐resolved thermodynamics and heterogeneous chemical reactions. Applied in the free‐running mode from 1995 to 2005 with prescribed historical climatological conditions, the model simulates global distributions of sulfate, nitrate, and ammonium in good agreement with observations and previous studies. Inclusion of nitrate resulted in ∼10% higher global average accumulation mode number concentrations, indicating enhanced growth of Aitken mode aerosols from nitrate formation. While the simulated accumulation mode nitrate burdens are high over the anthropogenic source regions, the sea‐salt and dust modes respectively constitute about 74% and 17% of the annual global average nitrate burden. Regional clear‐sky shortwave radiative cooling of up to −5 W m −2 due to nitrate is seen, with a much smaller global average cooling of −0.05 W m −2 . Significant enhancements in regional cloud condensation nuclei (at 0.1% supersaturation) and cloud droplet number concentrations are also attributed to nitrate, causing an additional global average shortwave cooling of −0.8 W m −2 . Taking into consideration of changes in both longwave and shortwave radiation under all‐sky conditions, the net change in the top of the atmosphere radiative fluxes induced by including nitrate aerosol is −0.7 W m −2 .

54 ENVIRONMENTAL SCIENCES↗

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗

Identification of a parasitic symbiosis between respiratory metabolisms in the biogeochemical chlorine cycle

A key step in the chlorine cycle is the reduction of perchlorate (ClO 4 - ) and chlorate (ClO 3 - ) to chloride by microbial respiratory pathways. Perchlorate-reducing bacteria and chlorate-reducing bacteria differ in that the latter cannot use perchlorate, the most oxidized chlorine compound. However, a recent study identified a bacterium with the chlorate reduction pathway dominating a community provided only perchlorate. Here we confirm a metabolic interaction between perchlorate- and chlorate-reducing bacteria and define its mechanism. Perchlorate-reducing bacteria supported the growth of chlorate-reducing bacteria to up to 90% of total cells in communities and co-cultures. Chlorate-reducing bacteria required the gene for chlorate reductase to grow in co-culture with perchlorate-reducing bacteria, demonstrating that chlorate is responsible for the interaction, not the subsequent intermediates chlorite and oxygen. Modeling of the interaction suggested that cells specialized for chlorate reduction have a competitive advantage for consuming chlorate produced from perchlorate, especially at high concentrations of perchlorate, because perchlorate and chlorate compete for a single enzyme in perchlorate-reducing cells. We conclude that perchlorate-reducing bacteria inadvertently support large populations of chlorate-reducing bacteria in a parasitic relationship through the release of the intermediate chlorate. An implication of these findings is that undetected chlorate-reducing bacteria have likely negatively impacted efforts to bioremediate perchlorate pollution for decades.

59 BASIC BIOLOGICAL SCIENCES↗

Free-energy analysis of the nonhysteretic first-order phase transition of Eu 2 In

Binary intermetallic Eu 2 In was recently reported to exhibit a giant anhysteretic magnetocaloric effect due to a first-order magnetic phase transition between paramagnetic and ferromagnetic states. Experimentally, the transition occurs with a small phase volume change, ΔV/V, of approximately 0.1% around T C of ca. 55 K. We representmagnetic and computemagnetocaloric properties of a Eu 2 In compound using a microscopic description based on a model Hamiltonian that takes into account magnetic exchange and magnetoelastic interactions. In the model the thermodynamic nature of the transition is conveniently represented by a single magnetoelastic interaction parameter. We report a good agreement between the theoretical results and earlier published experimental data confirms the effectiveness of our approach.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗