Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate computing”

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 829 records · Page 46

HEPTAPOD: Orchestrating High Energy Physics Workflows Towards Autonomous Agency

Many workflows in high-energy-physics (HEP) stand to benefit from recent advances in transformer-based large language models (LLMs). While early applications of LLMs focused on text generation and code completion, modern LLMs now support orchestrated agency: the coordinated execution of complex, multi-step tasks through tool use, structured context, and iterative reasoning. We introduce the HEP Toolkit for Agentic Planning, Orchestration, and Deployment (HEPTAPOD), an orchestration framework designed to bring this emerging paradigm to HEP pipelines. The framework enables LLMs to interface with domain-specific tools, construct and manage simulation workflows, and assist in common utility and data analysis tasks through schema-validated operations and run-card-driven configuration. To demonstrate these capabilities, we consider a representative Beyond the Standard Model (BSM) Monte Carlo validation pipeline that spans model generation, event simulation, and downstream analysis within a unified, reproducible workflow. HEPTAPOD provides a structured and auditable layer between human researchers, LLMs, and computational infrastructure, establishing a foundation for transparent, human-in-the-loop systems.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗

Toward a coordinated understanding of hydro‐biogeochemical root functions in tropical forests for application in vegetation models

Summary Tropical forest root characteristics and resource acquisition strategies are underrepresented in vegetation and global models, hampering the prediction of forest–climate feedbacks for these carbon‐rich ecosystems. Lowland tropical forests often have globally unique combinations of high taxonomic and functional biodiversity, rainfall seasonality, and strongly weathered infertile soils, giving rise to distinct patterns in root traits and functions compared with higher latitude ecosystems. We provide a roadmap for integrating recent advances in our understanding of tropical forest belowground function into vegetation models, focusing on water and nutrient acquisition. We offer comparisons of recent advances in empirical and model understanding of root characteristics that represent important functional processes in tropical forests. We focus on: (1) fine‐root strategies for soil resource exploration, (2) coupling and trade‐offs in fine‐root water vs nutrient acquisition, and (3) aboveground–belowground linkages in plant resource acquisition and use. We suggest avenues for representing these extremely diverse plant communities in computationally manageable and ecologically meaningful groups in models for linked aboveground–belowground hydro‐nutrient functions. Tropical forests are undergoing warming, shifting rainfall regimes, and exacerbation of soil nutrient scarcity caused by elevated atmospheric CO 2 . The accurate model representation of tropical forest functions is crucial for understanding the interactions of this biome with the climate.

Plant Sciences↗

Effects of negative triangularity on microinstabilities in a low-recycling lithium-wall spherical tokamak

In this work, we present a linear gyrokinetic study of the impact of negative triangularity (NT) on microinstabilities in the Lithium Tokamak eXperiment-β (LTX-β), a low-recycling spherical tokamak with liquid lithium plasma-facing components that produce flat electron-temperature profiles [Elliott et al., IEEE Trans. Plasma Sci. 48, 1382 (2020)]. While NT is widely recognized as a stabilizing mechanism and often associated with improved confinement in conventional tokamaks, this study reveals that its effect is not universally stabilizing in the parameter regime of LTX-β and is shown to be highly sensitive to local equilibria. Using local linear simulations with the GS2 code [Kotschenreuther et al., Comput. Phys. Commun. 88, 128 (1995)] at ρ=0.3, 0.5, and 0.8 for two representative discharges (#103955 and #109355), and employing the Miller equilibrium model to isolate shaping effects, we find that NT can transition from stabilizing to destabilizing depending on magnetic shear, safety factor, and electron-temperature gradient. In shot #103955, NT reduces growth rates across radii, with strongest stabilization at the edge, whereas in shot #109355, it is stabilizing in the core but destabilizing at mid-radius and edge under experimental conditions. Parametric scans show that flattening the electron-temperature profile, increasing magnetic shear, and reducing the safety factor recover NT stabilization. These results demonstrate that NT stabilization is tunable rather than intrinsic and requires coordinated control of magnetic geometry and gradient drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A symmetry algebra in double-scaled SYK

The double-scaled limit of the Sachdev-Ye-Kitaev (SYK) model takes the number of fermions and their interaction number to infinity in a coordinated way. In this limit, two entangled copies of the SYK model have a bulk description of sorts known as the “chord Hilbert space”. We analyze a symmetry algebra acting on this Hilbert space, generated by the two Hamiltonians together with a two-sided operator known as the chord number. This algebra is a deformation of the JT gravitational algebra, and it contains a subalgebra that is a deformation of the \mathfrak{sl}_2 𝔰 𝔩 2 near-horizon symmetries. The subalgebra has finite-dimensional unitary representations corresponding to matter moving around in a discrete Einstein-Rosen bridge. In a semiclassical limit the discreteness disappears and the subalgebra simplifies to \mathfrak{sl}_2 𝔰 𝔩 2 , but with a non-standard action on the boundary time coordinate. One can make the action of \mathfrak{sl}_2 𝔰 𝔩 2 algebra more standard at the cost of extending the boundary circle to include some “fake” portions. Such fake portions also accommodate certain subtle states that survive the semi-classical limit, despite oscillating on the scale of discreteness. We discuss applications of this algebra, including sub-maximal chaos, the traversable wormhole protocol, and a two-sided OPE.

97 MATHEMATICS AND COMPUTING↗

Energy Technology Innovation Partnership Project Strategic Energy Plan for Washington County, Maine

This Strategic Energy Plan (SEP) for Washington County, Maine was developed through the U.S. Department of Energy's Energy Technology Innovation Partnership Project (ETIPP) to address local challenges related to electricity reliability, affordability, and safety. Grounded in a "relational infrastructure" approach, the plan recognizes that improved energy outcomes depend not only on physical grid investments but also on strengthened coordination and trust between utilities and communities. Using utility-provided reliability data, public datasets, and community input, the plan identifies key risks, including vegetation-related outages, high energy burden, and vulnerability among elderly and medically dependent populations. The SEP establishes three priority areas: hazard mitigation for energy infrastructure, outage safety and emergency readiness, and energy assistance and thermal safety. It outlines community-driven actions designed to complement utility-led system improvements, including targeted outreach, municipal coordination, and expanded access to energy programs. The plan also situates local conditions within broader state policy, market dynamics, and funding opportunities. Intended as a living document, this SEP provides a scalable framework for collaborative action to improve resilience, safety, and affordability in rural energy systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Solutal-buoyancy-driven intertwining and rotation of patterned elastic sheets

Abstract The intertwining of strands into 3D spirals is ubiquitous in biology, enabling functions from information storage to maintenance of cell structure and directed locomotion. In synthetic systems, entwined fibers can provide superior mechanical properties and act as artificial muscle or structural reinforcements. Unlike structures in nature, the entwinement of synthetic materials typically requires application of an external stimulus, such as mechanical actuation, light, or a magnetic field. Herein, we use computational modeling to design microscale sheets that mimic biology by transducing chemical energy into mechanical action, and thereby self-organize and interlink into 3D spirals, which spontaneously rotate. These flexible sheets are immersed in a fluid-filled microchamber that encompasses an immobilized patch of catalysts on the bottom wall. The sheets themselves can be passive or active (coated with catalyst). Catalytic reactions in the solution generate products that occupy different volumes than the reactants. The resulting density variations exert a force on the fluid (solutal buoyancy force) that causes motion, which in turn drives the interlinking and collective swirling of the sheets. The individual sheets do not rotate; rotation only occurs when the sheets are interlinked. This level of autonomous, coordinated 3D structural organization, intertwining, and rotation is unexpected in synthetic materials systems operating without external controls. Using physical arguments, we identify dimensionless ratios that are useful in scaling these ideas to other systems. These findings are valuable for creating materials that act as “machines”, and directing soft matter to undergo self-sustained, multistep assembly that is governed by intrinsic chemical reactions.

Manna, Raj Kumar (ORCID:0000000322729487)↗

GPS Spoofing Mitigation and Timing Risk Analysis in Networked Phasor Measurement Units via Stochastic Reachability

To address phasor measurement unit (PMU) vulnerability to spoofing, we propose the use of a set-valued state estimation technique known as stochastic reachability (SR)-based distributed Kalman filter (DKF) that computes secure global positioning system (GPS) timing across a network of receivers. Utilizing SR, we estimate not only GPS time but also its stochastic reachable set, which is parameterized by probabilistic zonotope (p-Zonotope). While requiring known measurement error bounds in only non-spoofed conditions, we designed a two-tiered approach. We first performed measurement-level spoofing mitigation via deviation of a measurement innovation from its expected p-Zonotope. We then performed state-level timing risk analysis via a determination of the intersection probability of the estimated p-Zonotope with an unsafe set that violates IEEE C37.118.1a-2014 standards. Finally, we validated our SR-DKF algorithm by subjecting it to a simulated receiver network to coordinate signal-level spoofing. We demonstrate improved timing accuracy and successful spoofing mitigation via the use of our SR-DKF algorithm. We also validated the robustness of the estimated timing risk as the number of receivers were varied.

47 OTHER INSTRUMENTATION↗

Native glycosylated HIV-1 Env in silico ensemble

The computationally modeled datasets generated in this study, namely the natively glycosylated SOSIP Env and the uniform mannose-9 glycosylated SOSIP Env ensembles will be submitted as Mendeley dataset server as required by the journal publishers, for public availability. These include two sets of Protein Data Bank (PDB)ensembles of 1000 structures each. The files are in the standardized PDB format. Each of the constituent ensemble structures are separated by the ‘END’ tag within the files. Snapshots from the beginning and end of each data set are given here. The column descriptions are as follows: Col1 –ATOM tag; Col 2 –atom number; COL 3 –atom name; COL4 –residue name; COL 5 –chain identifier; COL 6 –residue number; COL. 7,8,9 –X,Y,Z coordinates of atom in 3Dspace; COL 10 –Occupancy Factor; COL 11 -Temperature factor; COL 12 –segment identifier; COL 13 –element symbol.

59 BASIC BIOLOGICAL SCIENCES↗

Advances in Uncertainty and Sensitivity Analysis Methods and Applications in GDSA Framework.

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Fuel Cycle Technology (FCT) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). Two high priorities for SFWST disposal R&D are design concept development and disposal system modeling. These priorities are directly addressed in the SFWST ''Geologic Disposal Safety Assessment'' (GDSA) control account, which is charged with developing a geologic repository system modeling and analysis capability, and the associated software, ''GDSA Framework'', for evaluating disposal system performance for nuclear waste in geologic media. ''GDSA Framework'' is supported by SFWST Campaign and its predecessor the Used Fuel Disposition (UFD) campaign. This report fulfills the GDSA Uncertainty and Sensitivity Analysis Methods work package (SF-20SN01030403) level 3 milestone — ''Advances in Uncertainty and Sensitivity Analysis Methods and Applications in GDSA Framework'' (M3SF-20SN010304032). It presents high level objectives and strategy for development of uncertainty and sensitivity analysis tools, demonstrates uncertainty quantification (UQ) and sensitivity analysis (SA) tools in GDSA Framework in FY20, and describes additional UQ/SA tools whose future implementation would enhance the UQ/SA capability of ''GDSA Framework''. This work was closely coordinated with the other Sandia National Laboratory GDSA work packages: the GDSA Framework Development work package (SF- 2051\101030404), the GDSA Repository Systems Analysis work package (SF-2051\101030405), and the GDSA PFLOTRAN Development work package (SF-20SN01030406). This report builds on developments reported in previous ''GDSA Framework'' milestones, particularly M2SF- 19SNO1030403.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

X-ray Crystal Structure of Thorium Tetrahydroborate, Th(BH 4 ) 4 , and Computational Studies of An(BH 4 ) 4 (An = Th, U)

Here the crystal structure of Th(BH 4 ) 4 is described. Two of the four BH 4 – ions are terminal and tridentate (κ 3 ), whereas the other two bridge between neighboring Th IV centers in a κ 2 ,κ 2 (i.e., bis-bidentate) fashion. Thus, each thorium center is bound to six BH 4 – groups by 14 Th–H bonds. The six boron atoms describe a distorted octahedron in which the κ 3 -BH 4 – ions are mutually cis; the 14 ligating hydrogen atoms define a highly distorted bicapped hexagonal antiprism. The thorium centers are linked into a polymer consisting of interconnected helical chains wound about 4-fold screw axes. The structures of An(BH 4 ) 4 (An = Th, U) were also investigated by DFT. The geometries of [An(BH 4 ) 6 ] 2– , [An3(BH 4 ) 16 ] 4– , and [An 5 (BH 4 ) 26 ] 6– fragments of the polymeric structures were optimized at the B3LYP and/or PBE levels. Most calculated geometries are 14-coordinate and agree with the experimental structures, but isolated [Th(BH 4 ) 6 ] 2– units are predicted to feature 16-coordinate Th IV centers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elucidating the Role of Hydrogen at c‐In 2 O 3 /a‐In 2 O 3− x Interface

Ab initio molecular dynamics simulations and hybrid functional electronic structure calculations are employed to determine the formation, the structural and electronic properties, and the dynamics of covalent (In—OH) and ionic (In—H—In) hydrogen defects at crystalline-In 2 O 3 /amorphous-In 2 O 3−x interface. This comprehensive computational study considers (i) various interstitial and substitutional hydrogen site locations within the crystalline, amorphous, and interfacial regions; (ii) several oxygen-to-hydrogen ratios; and (iii) possible defect charged states. The results reveal hydrogen's inability to fully passivate the undercoordinated under-shared in atoms in amorphous highly substoichiometric oxide, giving rise to the formation of deep electron traps even in net-charge neutral structures. These trap defects are found to be sensitive to photoexcitation, in contrast to In—OH with H's electronic states located below the valence band and to In—H—In, where H is found to maintain its charge state upon illumination. Nevertheless, H plays a critical role in photoinduced conductivity and its relaxation by promoting In—O coordination transformations at the interfacial region, including deterioration of the crystalline layer. Finally, the results help identify mobile H species and metastable H defect complexes (such as In—H—H—In, In—OH–H—In, and In—OH–O—In) that are responsible for long relaxation times of the conductivity decay.

H defects↗

Axially-deformed solution of the Skyrme-Hartree-Fock-Bogoliubov equations using the transformed harmonic oscillator basis (IV) HFBTHO (v4.0): A new version of the program

We describe the new version 4.0 of the code HFBTHO that solves the nuclear Hartree-Fock-Bogoliubov problem by using the deformed harmonic oscillator basis in cylindrical coordinates. In the new version, we have implemented the restoration of rotational, particle number, and reflection symmetry for even-even nuclei. The restoration of rotational symmetry does not require using bases closed under rotation. Furthermore, we added the SeaLL1 functional and improved the calculation of the Coulomb potential. Finally, we refactored the code to facilitate maintenance and future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Computational Study of an Iron(II) Polypyridine Electrocatalyst for CO 2 Reduction: Key Roles for Intramolecular Interactions in CO 2 Binding and Proton Transfer

A solar-driven conversion of CO 2 into fuels by artificial photosynthesis would not only mitigate the greenhouse effect but also provide an alternative to obtain fuels in a renewable fashion. To this end, the new iron polypyridine catalyst [Fe(bpy NHEt PY2Me)L 2 ] 2+ (L = H 2 O, CH 3 CN) was recently developed for the electrochemical reduction of CO 2 to CO. Here, we performed density functional theory (DFT) electronic structure calculations to shed light on a possible pathway for CO 2 reduction and the origin of the selectivity between CO 2 reduction versus the hydrogen evolution reaction. The metal center remains Lewis acidic throughout the reduction process due to ligand loss and mainly ligand-based reduction stabilized by antiferromagnetic coupling to a high-spin Fe(II) center. This results in a high barrier for hydride formation but a facile addition and activation of CO 2 via an η 2 coordination and stabilizing hydrogen bonding by the amine group. The second unoccupied equatorial coordination site opens up the possibility for an intramolecular protonation with a coordinated water ligand. This facilitates protonation because not only CO 2 but also the proton source H 2 O is activated and properly aligned for a proton transfer due to the Fe-OH 2 bond; consequently, both protonation steps are facile. The moderate ligand field allows a rapid ligand exchange for a second intramolecular protonation step and facilitates an exergonic CO release. The lower selectivity of the related [Fe(bpy OH PY2Me)L 2 ] 2+ complex can be related to its more acidic second coordination sphere, which enables an intramolecular proton transfer that is kinetically competitive with CO 2 addition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Ligands on the Ion–Ion Selectivity of Ligand-Appended-Pillar[ n ]arene Channels

Ligand-appended pillar[5]arene (LAP) channels are emerging as a promising platform for ion-selective membranes. In this study, we investigated the ion selectivities of three LAP channels functionalized with distinct ligands: diglycolamine, carbamoylmethyl phosphine oxide, and propionamide phosphonic acid. These ligands were chosen based on their demonstrated affinities for lanthanides in solvent extraction studies. We examined the selectivity of each channel toward a series of monovalent, divalent, and trivalent ions: Li + , Na + , K + , Mg 2+ , Ca 2+ , La 3+ , Eu 3+ , and Yb 3+ . To quantify ion-channel interactions and the energetics of translocation, we computed the potential of mean force (PMF) profiles for each ion. These PMFs were subsequently recast into permeability coefficients, enabling a direct evaluation of ion permselectivity. Our results reveal distinct selectivity patterns governed by the chemical nature of the appended ligands. The diglycolamine-functionalized channel exhibits strong interactions with monovalent ions, leading to replication of the Li + bulk hydration shell within the channel. This promotes preferential transport Li + ions while effectively excluding divalent and trivalent species. In contrast, the carbamoylmethyl phosphine oxide-functionalized channel displays a reduced selectivity between monovalent ions while similarly exhibiting an effective rejection of divalent and trivalent ions. This arises from steric hindrance around channel-lining oxygens imposed by the ligand's bulky phenyl groups which limits heavy ion coordination and promotes retention of hydration shells effectively increasing effective ion size. The propionamide phosphonic acid-functionalized channel exhibits a different trend, favoring ions with lower hydration free energies, with K + ions showing enhanced permeation. These findings highlight the critical role of ligand-ion interactions in modulating ion selectivity within LAP channels. In particular, the strong coordination exhibited by diglycolamine and propionamide phosphonic acid ligands significantly influences the ion transport energetics and selectivity profiles.

hydration↗

Concerted Electron-Ion Transport by Polyacrylonitrile Elucidated with Reactive Deep Learning Potentials

Charge transport in polymers, such as polyacrylonitrile (PAN), is crucial for electronics and energy storage. For instance, PAN can transport cations e.g., Li + , by facilitating dynamic cation-nitrile coordination in batteries. However, little is known regarding the underlying role of complex reactive polymer configurations. Herein, we develop a deep-learning potential, trained on ab initio energies and forces of nonequilibrium reactive PAN configurations, to unravel the kinetics of PAN cyclization initiated by a nucleophile (OH – dissociated from LiOH) attacking the terminal nitrile carbon. We find, based on the reaction free-energetics, rates, and charge analysis, that the nucleophile attack producing the first ring is the rate-limiting step, which subsequently triggers Li + -coupled electron transfer along the PAN backbone, causing ∼10 4 times faster sequential ring-formation of the remaining nitriles. PAN’s extended configurations, where dipolar and H-bonding interactions are minimal, enable such rapid kinetics. By validating our computational findings with IR and NMR experiments, we establish a pathway for designing reactive polymers with enhanced charge transport for energy applications.

Chahal-Crockett, Rajni [Oak Ridge National Laborat↗

Small arms suppression project (LLNL final report)

US Special Operations Command (USSOCOM) was seeking a technological leap in small firearms weapon suppressor technology, because anticipated enemy capabilities are requiring the operators to have smaller detection cross sections to ensure the safe execution of missions. Suppressors have been developed almost exclusively through trial-and-error methods since the time of the original design by Hiram Maxim over one hundred years ago. Consequently USSOCOM deemed it prudent to perform a physicsbased study of weapon suppression to understand performance limits and possibly identify breakthrough technologies. Lawrence Livermore National Laboratory’s (LLNL’s) high performance production level computational tool called ALE3D (Arbitrary Lagrangian-Eulerian 3D and 2D) has unique physics models and numerical algorithms for modeling suppressor dynamics. The flexible and extendable code framework supports fully integrated hydrodynamics, heat transfer, solid and fluid dynamics, and chemistry that can be applied to simulating propellant-driven motion of a bullet down a gun barrel, the transfer of heat from the burning propellant to the barrel and suppressor, the chemistry of muzzle flash, and the shock/acoustic/optical signatures in the near-field. LLNL’s originally anticipated role was to augment ALE3D for this task, by developing the software and analysis methodologies specific to the simulation of blast and muzzle flash phenomena. It was believed that insights provided by our ALE3D simulations in tandem with a coordinated experimental component by our other team members from Oak Ridge National Laboratory (ORNL) and the U. S. Army Armament, Research, Development and Engineering Center (ARDEC), would have excellent prospects of yielding useful suppressor design improvements that could be transitioned to industry and utilized by US Special Operations Command. The three year effort has come to fruition with the development of revolutionary suppressor designs that far outperform any previous or current design by anyone outside this multi-lab team.

42 ENGINEERING↗

Origin of Stabilization of Ligand-Centered Mixed Valence Ruthenium Azopyridine Complexes: DFT Insights for Neuromorphic Applications

Redox-driven conductance changes are critical processes in molecular- and coordination-complex-based memristive thin films and devices that are envisioned for neuromorphic technologies, but fundamental mechanisms of conductance switching are not fully understood. Here, we explore charge disproportionation (CD) processes in [Ru II L 2 ](PF 6 ) 2 molecular systems that intrinsically involve interfragment charge transfer (IFCT). Using a combination of ab initio molecular dynamics simulation (AIMD), time-dependent density functional theory (TD-DFT), and density functional theory (DFT) calculations, we investigate the electron transfer mechanisms and the roles of temperature and cell volumetric expansion in facilitating the counterion movements and electronic transitions required for low-cost IFCT and charge redistribution. A detailed analysis of the density of states and TD-DFT calculations highlights that unpaired electrons play a crucial role in low-energy transitions, with the azo (N=N) groups of the ligand serving as the primary sites for electronic transport between molecular fragments, further stabilizing the asymmetric state. Localization of added electrons on azo ligands occurs with negligible change at the Ru centers, supported by atomic volume expansions up to +4.74 bohr 3 , and goes along with a progressive reduction of the HOMO−LUMO gap across redox states, suggesting enhanced conductivity. The TD-DFT analysis reveals a dominant IFCT excitation at 2082.76 nm in the doubly reduced (22) state, while a stabilization energy of 1.20 eV of the asymmetric (13) state relative to the symmetric (22) state is predicted by constrained DFT. Periodic DFT and AIMD simulations emulating a molecular film show that the stabilization of the asymmetric state, relative to a symmetric one, translates in net charge separation values (order of ∼0.33 e) that are strongly linked to increased counterion mobility (average counterion displacements exceeding 0.7 Å per atom during CD events) and the involvement of azo groups in electron redistribution. These findings, which align with previously reported experimental and computational data, provide key insights into the IFCT mechanisms and electronic transport facilitated by azo groups, with important implications for redox-driven memristive and neuromorphic technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-driven design of molecular nanomagnets

Three decades of research in molecular nanomagnets have raised their magnetic memories from liquid helium to liquid nitrogen temperature thanks to a wise choice of the magnetic ion and coordination environment. Still, serendipity and chemical intuition played a main role. In order to establish a powerful framework for statistically driven chemical design, here we collected chemical and physical data for lanthanide-based nanomagnets, catalogued over 1400 published experiments, developed an interactive dashboard (SIMDAVIS) to visualise the dataset, and applied inferential statistical analysis. Our analysis shows that the Arrhenius energy barrier correlates unexpectedly well with the magnetic memory. Furthermore, as both Orbach and Raman processes can be affected by vibronic coupling, chemical design of the coordination scheme may be used to reduce the relaxation rates. Indeed, only bis-phthalocyaninato sandwiches and metallocenes, with rigid ligands, consistently present magnetic memory up to high temperature. Analysing magnetostructural correlations, we offer promising strategies for improvement, in particular for the preparation of pentagonal bipyramids, where even softer complexes are protected against molecular vibrations.

74 ATOMIC AND MOLECULAR PHYSICS↗